Abstract:The high cost of common multireference second-order perturbation theory (MRPT2) methods compared with the single-reference variant (MP2) arises from the expensive complete active space self-consistent field (CASSCF) orbital optimization step. Furthermore, the use of conventional four-index electron repulsion integrals (ERIs) prevents their application to larger molecular systems due to expensive I/O procedures. To address these bottlenecks of the multiconfigurational second-order quasidegenerate perturbation theory (MC-QDPT2), an efficient implementation of QDPT2 with the density-fitting (DF) and Cholesky decomposition (CD) approximations, denoted by DF-QDPT2 and CD-QDPT2, is reported. For the DF/CD-QDPT2 methods, the Hose--Kaldor approach is used. The DF-QDPT2 method, with the cc-pwCVTZ basis set, dramatically reduces the computational cost compared to conventional multiconfigurational QDPT2 (MC-QDPT2, from the {\sc Gamess 2017.R2} package), with a more than 122-fold reduction for the largest member of the diradical test set considered. The DF approximation enables substantially accelerated energies to be obtained for the QDPT2 approach due to the significantly reduced I/O time. The performance of the DF-QDPT2 and CD-QDPT2 methods is compared with that of CASSCF, the multireference second-order perturbation theory (MRMP2), MC-QDPT2, and CASSCF-based second-order perturbation theory (CASPT2) methods for singlet-triplet energy splitting ($E_{ST}$) in \ce{O2} and \ce{C2} molecules and for the dissociation energy of \ce{F2}. For the \ce{O2} and \ce{C2} molecules, the performance of the DF-QDPT2 and CD-QDPT2 methods is significantly better than that of CASSCF, MRMP2, MC-QDPT2, and CASPT2; while for the \ce{F2} case, the results of DF-QDPT2, CD-QDPT2, MRMP2, MC-QDPT2, and CASPT2 are similar and remarkably better than that of CASSCF, which fails dramatically. Moreover, the DF-QDPT2, CASSCF, CASPT2, and MRCI+Q methods are applied to potential energy curves (PECs) for \ce{N2}, \ce{CH4}, and \ce{F2} molecules. Our results demonstrate that the performance of DF-QDPT2 is substantially better than that of CASSCF, and comparable with that of CASPT2 for the molecules considered. Overall, the present application results demonstrate that the DF-QDPT2 and CD-QDPT2 methods are very promising for electronically challenging molecular systems suffering from (quasi)degeneracy problems, where the single-reference methods cannot provide an accurate electronic description, but the DF-QDPT2 and CD-QDPT2 methods can do so at significantly reduced computational costs.

Extended Koopmans' Theorem at the Second-Order Perturbation Theory.

Integration Approach at the Second-Order Perturbation Theory: Applications to Ionization Potential and Electron Affinity Calculations

Second-Order Perturbation Theory for Fractional Occupation Systems: Applications to Ionization Potential and Electron Affinity Calculations.

A Correlation-Relaxation-balanced Direct Method at the Second Order Perturbation Theory for Accurate Ionization Potential Predictions.

Perturbation Theory Made Efficient and Effective for Predictions of Ionization Potential and Electron Affinity.

Extended Koopmans' Theorem in the Adiabatic Connection Formalism: Applied to Doubly Hybrid Density Functionals.

Attaining high accuracy for charge-transfer excitations in non-covalent complexes at second-order perturbation cost: the importance of state-specific self-consistency

Single-reference Moller-Plesset perturbation extrapolation technique in the accurate model chemistry is unreliable for non-equilibrium molecules

Koopmans' theorem for large molecular systems within density functional theory.

Chemical Potentials and the One-Electron Hamiltonian of the Second-Order Perturbation Theory from the Functional Derivative Approach

Reaching high accuracy for energetic properties at second-order perturbation cost by merging self-consistency and spin-opposite scaling

Molecular Orbital Analysis In Evaluation Of Electron-Transfer Matrix Element By Koopmans' Theory

Second-order Many-Body Perturbation Theory: an Eternal Frontier.

Numerical Analysis of the Complete Active-Space Extended Koopmans's Theorem

Numerical analysis of the complete active-space extended Koopmans’s theorem

Insights into Direct Methods for Predictions of Ionization Potential and Electron Affinity in Density Functional Theory.

Linear and Nonlinear Optical Properties from TDOMP2 Theory

An Efficient Implementation of the Second-Order Quasidegenerate Perturbation Theory with Density-Fitting and Cholesky Decomposition Approximations: Is It Possible to Use Hartree-Fock Orbitals for a Multiconfigurational Perturbation Theory?

Theoretical Study on Low-Lying Electronic States of Kr2(+), Xe2(+), and Rn2(+).

A simple scheme for the direct calculation of ionization potentials with coupled-cluster theory that exploits established excitation energy methods

Koopmans-compliant functionals and potentials and their application to the GW100 test-set