Equilibrium Geometries of Noncovalently Bound Intermolecular Complexes Derived from Subsystem Formulation of Density Functional Theory

Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces

The theory of analytic gradients is presented for the projector-based density functional theory (DFT) embedding approach utilizing the Huzinaga-equation. The advantages of the Huzinaga-equation-based formulation are demonstrated. In particular, it is shown that the projector employed does not appear in the Lagrangian, and the potential risk of numerical problems is avoided at the evaluation of the gradients. The efficient implementation of the analytic gradient theory is presented for approaches where hybrid DFT, second-order Møller-Plesset perturbation theory, or double hybrid DFT are embedded in lower-level DFT environments. To demonstrate the applicability of the method and to gain insight into its accuracy, it is applied to equilibrium geometry optimizations, transition state searches, and potential energy surface scans. Our results show that bond lengths and angles converge rapidly with the size of the embedded system. While providing structural parameters close to high-level quality for the embedded atoms, the embedding approach has the potential to relax the coordinates of the environment as well. Our demonstrations on a 171-atom zeolite and a 570-atom protein system show that the Huzinaga-equation-based embedding can accelerate (double) hybrid gradient computations by an order of magnitude with sufficient active regions and enables affordable force evaluations or geometry optimizations for molecules of hundreds of atoms.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Analytical Nuclear Excited-State Gradients for the Second-Order Approximate Coupled-Cluster Singles and Doubles (CC2) Method Employing Uncoupled Frozen-Density Embedding

We report the derivation and implementation of analytical orbital-relaxed properties and nuclear gradients for excited states using the second-order approximate coupled-cluster singles and doubles (CC2) model combined with uncoupled frozen-density embedding (FDEu). An implementation of the algebraic diagrammatic construction through second-order ADC(2), which arises from simplification of RICC2 FDEu, is also presented. In order to ensure a RICC2 FDEu Lagrange functional that is linear in the Lagrange multipliers, the Hartree-Fock density is employed for the target subsystem in the embedding contributions. The accuracy of the new scheme is assessed using the carbon monoxide molecule, 4-aminophthalimide, and a benzonitrile dimer, revealing that the obtained errors are below the method error of RICC2. Using density functional theory for the environment, the efficiency of the new method is illustrated by computing the perturbed excited-state dipole moment of a chromophore in a biological environment. For this system, comprising 32 molecules consisting of 366 atoms in total, the computation requires only a couple of days on a standard compute node. RICC2 FDEu thus enables large-scale calculations of ab initio wave functions for molecules in complex environments as routine applications.

Explicit vs. implicit electronic polarisation of environment of an embedded chromophore in frozen-density embedding theory

A comparison of strategies to account for environment polarisation in Frozen Density Embedding Theory (FDET).

Analytical nuclear excited-state gradients for the Tamm-Dancoff approximation using uncoupled frozen-density embedding

We report the derivation and implementation of analytical nuclear gradients for excited states using time-dependent density functional theory using the Tamm-Dancoff approximation combined with uncoupled frozen-density embedding using density fitting. Explicit equations are presented and discussed. The implementation is able to treat singlet as well as triplet states and functionals using the local density approximation, the generalized gradient approximation, combinations with Hartree-Fock exchange (hybrids), and range-separated functionals such as CAM-B3LYP. The new method is benchmarked against supermolecule calculations in two case studies: The solvatochromic shift of the (vertical) fluorescence energy of 4-aminophthalimide on solvation, and the first local excitation of the benzonitrile dimer. Whereas for the 4-aminophthalimide-water complex deviations of about 0.2 eV are obtained to supermolecular calculations, for the benzonitrile dimer the maximum error for adiabatic excitation energies is below 0.01 eV due to a weak coupling of the subsystems. © 2017 Wiley Periodicals, Inc.

Analytical gradients for excitation energies from frozen-density embedding

Analytical excitation-energy gradients from frozen-density embedding–time-dependent density functional theory are derived and implemented, which are important for photochemistry in complex systems.

Subsystem-DFT potential-energy curves for weakly interacting systems

Subsystem DFT for weakly interacting systems: good at reference geometries, but not necessarily good for potential-energy curves.

Frozen density embedding with non-integer subsystems' particle numbers

We extend the frozen density embedding theory to non-integer subsystems' particles numbers. Different features of this formulation are discussed, with special concern for approximate embedding calculations. In particular, we highlight the relation between the non-integer particle-number partition scheme and the resulting embedding errors. Finally, we provide a discussion of the implications of the present theory for the derivative discontinuity issue and the calculation of chemical reactivity descriptors.

On the accuracy of frozen density embedding calculations with hybrid and orbital-dependent functionals for non-bonded interaction energies

We analyze the accuracy of the frozen density embedding (FDE) method, with hybrid and orbital-dependent exchange-correlation functionals, for the calculation of the total interaction energies of weakly interacting systems. Our investigation is motivated by the fact that these approaches require, in addition to the non-additive kinetic energy approximation, also approximate non-additive exact-exchange energies. Despite this further approximation, we find that the hybrid/orbital-dependent FDE approaches can reproduce the total energies with the same accuracy (about 1 mHa) as the one of conventional semi-local functionals. In many cases, thanks to error cancellation effects, hybrid/orbital-dependent approaches yield even the smallest error. A detailed energy-decomposition investigation is presented. Finally, the Becke-exchange functional is found to reproduce accurately the non-additive exact-exchange energies also for non-equilibrium geometries. These performances are rationalized in terms of a reduced-gradient decomposition of the non-additive exchange energy.

Multilevel X-Pol: a fragment-based method with mixed quantum mechanical representations of different fragments

The explicit polarization (X-Pol) method is a fragment-based quantum mechanical model, in which a macromolecular system or other large or complex system in solution is partitioned into monomeric fragments. The present study extends the original X-Pol method, where all fragments are treated using the same electronic structure theory, to multilevel representations, called multilevel X-Pol, in which different electronic structure methods are used to describe different fragments. The multilevel X-Pol method has been implemented into a locally modified version of Gaussian 09. A key ingredient that is used to couple interfragment electrostatic interactions at different levels of theory is the use of the response density for the post-self-consistent-field energy. (The response density is also called the generalized density.) The method is useful for treating fragments in a small region of the system such as a solute molecule or the substrate and amino acids in the active site of an enzyme with a high-level theory, and the fragments in the rest of the system by a lower-level and computationally more efficient method. Multilevel X-Pol is illustrated here by applications to hydrogen bonding complexes in which one fragment is treated with the hybrid M06 density functional, Møller-Plesset perturbation theory, or coupled cluster theory, and the other fragments are treated by Hartree-Fock theory or the B3LYP or M06 hybrid density functionals.

Energy decomposition analysis based on a block-localized wavefunction and multistate density functional theory

An interaction energy decomposition analysis method based on the block-localized wavefunction (BLW-ED) approach is described. The first main feature of the BLW-ED method is that it combines concepts of valence bond and molecular orbital theories such that the intermediate and physically intuitive electron-localized states are variationally optimized by self-consistent field calculations. Furthermore, the block-localization scheme can be used both in wave function theory and in density functional theory, providing a useful tool to gain insights on intermolecular interactions that would otherwise be difficult to obtain using the delocalized Kohn-Sham DFT. These features allow broad applications of the BLW method to energy decomposition (BLW-ED) analysis for intermolecular interactions. In this perspective, we outline theoretical aspects of the BLW-ED method, and illustrate its applications in hydrogen-bonding and π-cation intermolecular interactions as well as metal-carbonyl complexes. Future prospects on the development of a multistate density functional theory (MSDFT) are presented, making use of block-localized electronic states as the basis configurations.

Block-Localized Density Functional Theory (BLDFT), Diabatic Coupling, and Their Use in Valence Bond Theory for Representing Reactive Potential Energy Surfaces

A multistate density functional theory in the framework of the valence bond model is described. The method is based on a block-localized density functional theory (BLDFT) for the construction of valence-bond-like diabatic electronic states and is suitable for the study of electron transfer reactions and for the representation of reactive potential energy surfaces. The method is equivalent to a valence bond theory with the treatment of the localized configurations by using density functional theory (VBDFT). In VBDFT, the electron densities and energies of the valence bond states are determined by BLDFT. A functional estimate of the off-diagonal matrix elements of the VB Hamiltonian is proposed, making use of the overlap integral between Kohn-Sham determinants and the exchange-correlation functional for the ground state substituted with the transition (exchange) density. In addition, we describe an approximate approach, in which the off-diagonal matrix element is computed by wave function theory using block-localized Kohn-Sham orbitals. The key feature is that the electron density of the adiabatic ground state is not directly computed nor used to obtain the ground-state energy; the energy is determined by diagonalization of the multistate valence bond Hamiltonian. This represents a departure from the standard single-determinant Kohn-Sham density functional theory. The multistate VBDFT method is illustrated by the bond dissociation of H2+ and a set of three nucleophilic substitution reactions in the DBH24 database. In the dissociation of H2+, the VBDFT method yields the correct asymptotic behavior as the two protons stretch to infinity, whereas approximate functionals fail badly. For the S(N)2 nucleophilic substitution reactions, the hybrid functional B3LYP severely underestimates the barrier heights, while the approximate two-state VBDFT method overcomes the self-interaction error, and overestimates the barrier heights. Inclusion of the ionic state in a three-state model, VBDFT(3), significantly improves the computed barrier heights, which are found to be in accord with accurate results. The BLDFT method is a versatile theory that can be used to analyze conventional DFT results to gain insight into chemical bonding properties, and it is illustrated by examining the intricate energy contributions to the ion-dipole complex stabilization.