The restricted active space followed by second-order perturbation theory method: theory and application to the study of CuO2 and Cu2O2 systems

A multireference second-order perturbation theory using a restricted active space self-consistent field wave function as reference (RASPT2/RASSCF) is described. This model is particularly effective for cases where a chemical system requires a balanced orbital active space that is too large to be addressed by the complete active space self-consistent field model with or without second-order perturbation theory (CASPT2 or CASSCF, respectively). Rather than permitting all possible electronic configurations of the electrons in the active space to appear in the reference wave function, certain orbitals are sequestered into two subspaces that permit a maximum number of occupations or holes, respectively, in any given configuration, thereby reducing the total number of possible configurations. Subsequent second-order perturbation theory captures additional dynamical correlation effects. Applications of the theory to the electronic structure of complexes involved in the activation of molecular oxygen by mono- and binuclear copper complexes are presented. In the mononuclear case, RASPT2 and CASPT2 provide very similar results. In the binuclear cases, however, only RASPT2 proves quantitatively useful, owing to the very large size of the necessary active space.     

Time-dependent quasirelativistic density-functional theory based on the zeroth-order regular approximation

A time-dependent quasirelativistic density-functional theory for excitation energies of systems containing heavy elements is developed, which is based on the zeroth-order regular approximation (ZORA) for the relativistic Hamiltonian and a noncollinear form for the adiabatic exchange-correlation kernel. To avoid the gauge dependence of the ZORA Hamiltonian a model atomic potential, instead of the full molecular potential, is used to construct the ZORA kinetic operator in ground-state calculations. As such, the ZORA kinetic operator no longer responds to changes in the density in response calculations. In addition, it is shown that, for closed-shell ground states, time-reversal symmetry can be employed to simplify the eigenvalue equation into an approximate form that is similar to that of time-dependent nonrelativistic density-functional theory. This is achieved by invoking an independent-particle approximation for the induced density matrix. The resulting theory is applied to investigate the global potential-energy curves of low-lying LambdaS- and omega omega-coupled electronic states of the AuH molecule. The derived spectroscopic parameters, including the adiabatic and vertical excitation energies, equilibrium bond lengths, harmonic and anharmonic vibrational constants, fundamental frequencies, and dissociation energies, are in good agreement with those of time-dependent four-component relativistic density-functional theory and ab initio multireference second-order perturbation theory. Nonetheless, this two-component relativistic version of time-dependent density-functional theory is only moderately advantageous over the four-component one as far as computational efforts are concerned.     

A CSF-based multireference coupled pair approximation

Summary MRCPA (multireference coupled pair approximation) is formulated by the use of the wave operator formalism and Rayleigh-Schrödiner perturbation theory with special selection of the unperturbed part of the electronic Hamiltonian. By considering super molecule, it is shown that the theory is size consistent with the help of the new formalism. The method has been tested for three simple systems, H₂O, FH, and O₂.     

Simplified reference wave functions for multireference perturbation theory

A number of simplifications in defining the reference wave functions used in multireference second-order M?ller-Plesset perturbation theory (MRMP2) calculations are studied. The usual multiconfigurational orbital optimization is avoided by using Hartree-Fock or Kohn-Sham orbitals; the complete configuration expansion in the active-space orbitals is replaced by a severely truncated expansion, and the spin-component-scaling idea is applied to the multireference perturbation expansion. We assess these approximations to the full procedure by calculating the barrier heights for 15 processes taken from the Zhao-Gonzalez-Garcia-Truhlar database. Our results suggest that reliable and relatively cheap reference wave functions for MRMP2 calculations can be obtained from the simplifications introduced here. We hope that this will enable the application of the MRMP2 method to a larger range of chemical systems.     

Communication: extended multi-state complete active space second-order perturbation theory: energy and nuclear gradients

The extended multireference quasi-degenerate perturbation theory, proposed by Granovsky [J. Chem. Phys. 134, 214113 (2011)], is combined with internally contracted multi-state complete active space second-order perturbation theory (XMS-CASPT2). The first-order wavefunction is expanded in terms of the union of internally contracted basis functions generated from all the reference functions, which guarantees invariance of the theory with respect to unitary rotations of the reference functions. The method yields improved potentials in the vicinity of avoided crossings and conical intersections. The theory for computing nuclear energy gradients for MS-CASPT2 and XMS-CASPT2 is also presented and the first implementation of these gradient methods is reported. A number of illustrative applications of the new methods are presented.     

On multireference superdirect configuration interaction in third order

The superdirect configuration interaction (Sup-CI ) method has the usual versatility and stability of the CI methods with computational efficiency typical to that of the many-body methods, such as the many-body perturbation theory (MBPT ). The Hamilton operator is projected into a space of a few trial vectors, such as Krylov, Nesbet, or Møller–Plesset correction vectors. In this space, Hamiltonian matrix elements may be directly computed in the many-body fashion, as weighted sums of integral products over orbital indices. The variation-perturbation method based on the first-order wave function is equivalent to the Sup-CI method with a single correction vector of the Møller–Plesset type. Different points of view on the superdirect CI method are discussed and a version in which third-order contributions are computed for a relatively small (10–100) space of reference and correction vectors is tested. Selection of the best “effective first-order spaces” and size-extensivity corrections in Sup-CI are briefly discussed. Møoller–Plesset, Epstein–Nesbet, and other correction vectors are included in the model calculations on the symmetric stretch of bonds in water, acetylene, and the NH₂ molecule. Errors are almost independent of molecular geometry and the method appears to be superior than the multireference second-order perturbation methods. © 1994 John Wiley & Sons, Inc.     

General orbital invariant MP2-F12 theory

A general form of orbital invariant explicitly correlated second-order closed-shell Moller-Plesset perturbation theory (MP2-F12) is derived, and compact working equations are presented. Many-electron integrals are avoided by resolution of the identity (RI) approximations using the complementary auxiliary basis set approach. A hierarchy of well defined levels of approximation is introduced, differing from the exact theory by the neglect of terms involving matrix elements over the Fock operator. The most accurate method is denoted as MP2-F12/3B. This assumes only that Fock matrix elements between occupied orbitals and orbitals outside the auxiliary basis set are negligible. For the chosen ansatz for the first-order wave function this is exact if the auxiliary basis is complete. In the next lower approximation it is assumed that the occupied orbital space is closed under action of the Fock operator [generalized Brillouin condition (GBC)]; this is equivalent to approximation 2B of Klopper and Samson [J. Chem. Phys. 116, 6397 (2002)]. Further approximations can be introduced by assuming the extended Brillouin condition (EBC) or by neglecting certain terms involving the exchange operator. A new approximation MP2-F12/3C, which is closely related to the MP2-R12/C method recently proposed by Kedzuch et al. [Int. J. Quantum Chem. 105, 929 (2005)] is described. In the limit of a complete RI basis this method is equivalent to MP2-F12/3B. The effect of the various approximations (GBC, EBC, and exchange) is tested by studying the convergence of the correlation energies with respect to the atomic orbital and auxiliary basis sets for 21 molecules. The accuracy of relative energies is demonstrated for 16 chemical reactions. Approximation 3C is found to perform equally well as the computationally more demanding approximation 3B. The reaction energies obtained with smaller basis sets are found to be most accurate if the orbital-variant diagonal Ansatz combined with localized orbitals is used for the first-order wave function. This unexpected result is attributed to geminal basis set superposition errors present in the formally more rigorous orbital invariant methods.     

One-electron properties of molecules calculated using second-order multireference perturbation theory

A new form of second-order multireference perturbation theory coupled with finite-field perturbation theory is applied to evaluate some one-electron molecular properties. Several possible definitions of the zeroth-order Hamiltonian are considered and results tested against bench-mark full CI calculations. © 1995 John Wiley & Sons, Inc.     

Elimination of translational and rotational motions in nuclear orbital plus molecular orbital theory: application of Moller-Plesset perturbation theory

The translation- and rotation-free nuclear orbital plus molecular orbital (TRF-NOMO) theory was developed to determine the nonadiabatic nuclear and electronic wave functions. This study presents a formulation of TRF-NOMO second-order Moller-Plesset (MP2) perturbation and Epstein-Nesbet (EN) theory with the use of the TRF Hamiltonian. Numerical assessment of the TRF-NOMO/MP2 and EN is performed for several molecules. We confirm the importance of the elimination of translational and rotational motions in the many-body calculations.     

Complete-active-space self-consistent-field/Amber parameterization of the Lys296–retinal–Glu113 rhodopsin chromophore-counterion system

A special hybrid quantum mechanics/molecular mechanics forcefield is defined, parameterized and validated for studying the photoisomerization path of the retinal chromophore in the rhodopsin protein. It couples a multireference ab initio Hamiltonian (CASSCF and second-order multireference many-body perturbation theory using a CASSCF reference) to describe the chromophore while the rest of the protein is approximated with the Amber forcefield. The frontier has been carefully parameterized in order to reproduce full quantum mechanics torsional energy profiles, for both the ground state and the first excited state. It is also shown that replacing the chromophore counterion with point charges is a valid approximation. This result is interpreted in terms of a cancellation effect for which a possible explanation is given.     

Second order explicitly correlated R12 theory revisited: a second quantization framework for treatment of the operators' partitionings

Second order R12 theory is presented and derived alternatively using the second quantized hole-particle formalism. We have shown that in order to ensure the strong orthogonality between the R12 and the conventional part of the wave function, the explicit use of projection operators can be easily avoided by an appropriate partitioning of the involved operators to parts which are fully describable within the computational orbital basis and complementary parts that involve imaginary orbitals from the complete orbital basis. Various Hamiltonian splittings are discussed and computationally investigated for a set of nine molecules and their atomization energies. If no generalized Brillouin condition is assumed, with all relevant partitionings the one-particle contribution arising in the explicitly correlated part of the first order wave function has to be considered and has a significant role when smaller atomic orbital basis sets are used. The most appropriate Hamiltonian splitting results if one follows the conventional perturbation theory for a general non-Hartree-Fock reference. Then, no couplings between the R12 part and the conventional part arise within the first order wave function. The computationally most favorable splitting when the whole complementary part of the Hamiltonian is treated as a perturbation fails badly. These conclusions also apply to MP2-F12 approaches with different correlation factors.     

Multireference second-order perturbation theory: how size consistent is "almost size consistent"

A systematic study of the deviation from size consistency of the multireference second-order Moller-Plesset perturbation theory (MRMP2) method is presented. The size-consistency error is shown to depend on the number of monomers in a supermolecule calculation, size of basis set, number of correlated valence electrons, and size of active space. HF, F(2), and N(2) are used as test cases, with stretched bonds, to include simple, well-defined multireference character. This is essential in ensuring that MRMP2 is being tested as a multireference method. It is concluded that the MRMP2 and other multireference perturbation theory methods can exhibit significant size-consistency errors, and that the size of the error depends on the manner in which the perturbation theory is implemented.     

New perspectives in multireference perturbation theory: the n-electron valence state approach

The n-electron valence state perturbation theory (NEVPT) is a form of multireference perturbation theory which is based on a zero order reference wavefunction of CAS-CI type (complete active space configuration interaction) and which is characterized by the utilization of correction functions (zero order wavefunctions external to the CAS) of multireference nature, obtained through the diagonalization of a suitable two-electron model Hamiltonian (Dyall’s Hamiltonian) in some well defined determinant spaces. A review of the NEVPT approach is presented, starting from the original second order state-specific formulation, going through the quasidegenerate multi-state extension and arriving at the recent implementations of the third order in the energy and of the internally contracted configuration interaction. The chief properties of NEVPT—size consistence and absence of intruder states—are analyzed. Finally, an application concerning the calculation of the vertical spectrum of the biologically important free base porphin molecule, is presented.     

Influence of multi-atom bridging ligands on the electronic structure and magnetic properties of homodinuclear titanium molecules

The electronic structure and magnetic properties of homodinuclear titanium(III) molecules with bridging ligands from groups 14, 15, and 16 are examined. Single- and multireference methods with triple-zeta plus polarization basis sets are employed. Dynamic electron correlation effects are included via second-order multireference perturbation theory. Isotropic interaction parameters are calculated, and two of the complexes studied are predicted to be ferromagnetic based on multireference second-order perturbation (MRMP2) theory, using the TZVP(fg) basis set. Zero-field splitting parameters are determined using spin-orbit coupling obtained from complete active space (CAS) self-consistent field (SCF) and multiconfigurational quasi-degenerate perturbation theory (MCQDPT) wave functions. Three Breit-Pauli-based spin coupling methods were employed: full Breit-Pauli (HSO2), the partial two-electron method (P2E), and the semiempirical one-electron method (HSO1).     

The centrifugally decoupled exponential distorted wave (CDEDW) approximation for the calculation of rotationally inelastic molecular collision cross se

A new approximate method is presented for the rapid calculation of rotationally inelastic molecular collision cross sections. The method is called the centrifugally decoupled exponential distorted wave (CDEDW) approximation and involves the combination of two well known approximations. The first approximation is the neglect of the off-diagonal coupling terms which arise from the orbital angular momentum operator in the coupled differential equations in the body-fixed axis system. The second approximation is to treat the remaining coupling terms, which arise from the interaction potential, using a unitary perturbation approximation. The CDEDW method is applied to the calculation of total and partial rotationally inelastic cross sections in the ArN₂ system, and detailed comparisons are made with exact and several other types of approximate calculations. Agreement with exact calculations is good and often comparable with the coupled states and p-helicity decoupled approximations. The CDEDW method requires a similar amount of computational effort to the infinite order sudden (IOS) approximation, and we show that for the present system the CDEDW method gives more reliable results.     

Multireference composite approaches for the accurate study of ground and excited electronic states: C2, N2, and O2

A multireference analog of the correlation consistent composite approach (MR-ccCA) based on complete active space with second-order perturbation theory (CASPT2) has been utilized in an investigation of the ground and valence excited states of C(2), N(2), and O(2). The performance of different second-order multireference perturbation theory methods including second-order n-electron valence state perturbation theory, second-order multireference M?ller-Plesset, and second-order generalized van Vleck perturbation theory has been analyzed as potential alternatives to CASPT2 within MR-ccCA. The MR-ccCA-P predicts spectroscopic constants with overall mean absolute deviations from experimental values of 0.0006 A?, 7.0 cm(-1), and 143 cm(-1) for equilibrium bond length (r(e)), harmonic frequency (ω(e)), and term values (T(e)), respectively, which are comparable to the predictions by more computationally costly multireference configuration interaction-based methods.