The electronic structure of the triiodide ion from relativistic correlated calculations: A comparison of different methodologies J. Chem. Phys. 133, 064305 (2010); 10.1063/1.3474571 Relativistic quasidegenerate perturbation theory with four-component general multiconfiguration reference functionsThe 1 2 A 1 , 1 2 B 2 , and 1 2 A 2 states of the SO 2 + ion studied using multiconfiguration second-order perturbation theory Generalized molecular mechanics including quantum electronic structure variation of polar solvents. I. Theoretical formulation via a truncated adiabatic basis set description

First results on electron densities and energies for a number of molecular complexes with different interaction strengths (ranging from ca. 0.3 to 40 kcal/mol), obtained using our recently introduced DFT-in-DFT embedding equations (i.e., Kohn-Sham equations with constrained electron density (KSCED) and external orbital orthogonality (ext orth), KSCED(x, ext orth), where x denotes the single particle support: monomer (m); supermolecular (s); or extended monomer (e)) are compared with densities from supermolecular Kohn-Sham (KS)-DFT calculations and traditional DFT-in-DFT results. Because our methodology does not rely on error-prone potentials that are not present in supermolecular KS-DFT calculations, it allows DFT-in-DFT calculations to achieve much higher accuracy than previous protocols of DFT-in-DFT that employed such potentials. It is shown that whereas conventional DFT-in-DFT embedding theory leads to errors in the electron density at the boundary between subsystems, the situation is remedied when orbital orthogonality between subsystems (i.e., external orthogonality) is enforced. Our approach reproduces KS-DFT total energies at least to the seventh decimal place (and exactly at most geometries) for the tested systems. Potential energy curves (PECs) of the separation of some of the tested systems into fragments are calculated. PECs, obtained with the new equations, using the usual Kohn-Sham equations with constrained electron density and supermolecular basis expansion [KSCED(s, ext orth, v(T) = 0), where v(T) is the nonadditive kinetic potential] were found to be virtually identical to those from conventional KS-DFT; equilibrium distances and interaction energies were reproduced to all reported digits for both local density approximation (LDA) and generalized gradient approximation (GGA) functionals. As an additional approximation, an alternative one-particle space (to the common monomer or supermolecular spaces) in which KS orbitals of a subsystem are expanded is introduced. This expansion, which we refer to as the extended monomer expansion [e.g., KSCED(e)], includes basis functions centered on atom(s) of the complementary subsystem in the interfacial region. Density differences and PECs obtained with the new equations and new one-particle space [i.e., KSCED(e, ext orth, v(T) = 0)] were closely related to those obtained from KSCED(s, ext orth, v(T) = 0). The new approach does not require any supermolecular calculations.

No abstract

Using the technique of Perdew and Levy [Phys. Rev. B 31, 6264 (1985)], it is shown that both the density function theory (DFT)-in-DFT and wave function theory (WFT)-in-DFT embedding approaches are formally correct in studying not only the ground state but also a subset of the excited states of the total system. Without further approximations, the DFT-in-DFT embedding approach results in a pair of coupled Euler-Lagrange equations. In contrast to DFT-in-DFT, the WFT-in-DFT approach is shown to ensure a systematic description of excited states if such states are mainly related to excitations within the embedded subsystem. Possible ways for the practical realization of the WFT-in-DFT approach for studying excited states are briefly discussed.

A new, efficient, configuration-driven algorithm utilizing the unitary group approach (UGA) was developed and implemented for the generalized van Vleck perturbation theory (GVVPT) variant of multireference perturbation theory. The computational speed has been improved by 1 or 2 orders of magnitude compared to the previous implementation based on the Table-CI technique. It is shown that the reformulation is applicable to both the second-order (GVVPT2) and third-order (GVVPT3) approximations. Calculations on model problems and on a chemically realistic description of cyclobutadiene are used to illustrate the performance of the method. The calculations on cyclobutadiene, using over 2.3 billion CSFs, provide results on geometric parameters and the barrier height of the automerization reaction in good agreement with established high accuracy results.