Hybrid correlation models based on active-space partitioning: correcting second-order Moller-Plesset perturbation theory for bond-breaking reactions.

Moller-Plesset second-order (MP2) perturbation theory breaks down at molecular geometries which are far away from equilibrium. We decompose the MP2 energy into contributions from different orbital subspaces and show that the divergent behavior of the MP2 energy comes from the excitations located within a small (or sometimes even the minimal) active space. The divergent behavior of the MP2 energy at large interfragment distances may be corrected by replacing a small number of terms by their more robust counterparts from coupled-cluster (CCSD) theory. We investigated several schemes of such a substitution, and we find that a coupling between the active-space CCSD and the remaining MP2 amplitudes is necessary to obtain the best results. This naturally leads us to an approach which has previously been examined in the context of cost-saving approximations to CCSD for equilibrium properties by Nooijen [J. Chem. Phys. 111, 10815 (1999)]. The hybrid MP2-CCSD approach, which has the same formal scaling as conventional MP2 theory, provides potential curves with a correct shape for bond-breaking reactions of BH, CH(4), and HF. The error of the MP2-CCSD method (measured against full configuration-interaction data) is smaller than that of MP2 at all interfragment separations and is qualitatively similar to that of full CCSD.

[1]  Jeppe Olsen,et al.  Full configuration interaction benchmarking of coupled-cluster models for the lowest singlet energy surfaces of N2 , 2000 .

[2]  Marcel Nooijen,et al.  Combining coupled cluster and perturbation theory , 1999 .

[3]  Antara Dutta,et al.  Full configuration interaction potential energy curves for breaking bonds to hydrogen: An assessment of single-reference correlation methods , 2003 .

[4]  J. Paldus,et al.  Reduced multireference coupled cluster method IV: open-shell systems , 2000 .

[5]  Jae Shin Lee,et al.  Critical comparison of single-reference and multireference coupled-cluster methods: Geometry, harmonic frequencies, and excitation energies of N2O2 , 2003 .

[6]  T. Crawford,et al.  On the energy invariance of open‐shell perturbation theory with respect to unitary transformations of molecular orbitals , 1996 .

[7]  Martin Head-Gordon,et al.  A second-order perturbative correction to the coupled-cluster singles and doubles method: CCSD(2) , 2001 .

[8]  Michael J. McGuire,et al.  Recent advances in electronic structure theory: Method of moments of coupled-cluster equations and renormalized coupled-cluster approaches , 2002 .

[9]  Henry F. Schaefer,et al.  Compact Variational Wave Functions Incorporating Limited Triple and Quadruple Substitutions , 1996 .

[10]  Hans-Joachim Werner,et al.  Local treatment of electron correlation in coupled cluster theory , 1996 .

[11]  Isaiah Shavitt,et al.  The Method of Configuration Interaction , 1977 .

[12]  Ludwik Adamowicz,et al.  A state-selective multireference coupled-cluster theory employing the single-reference formalism , 1993 .

[13]  Micah L. Abrams,et al.  General-order single- and multi-reference configuration interaction and coupled-cluster theory: Symmetric dissociation of water , 2005 .

[14]  P. Löwdin Quantum Theory of Many-Particle Systems. III. Extension of the Hartree-Fock Scheme to Include Degenerate Systems and Correlation Effects , 1955 .

[15]  Anna I. Krylov,et al.  A spin-complete version of the spin-flip approach to bond breaking: What is the impact of obtaining spin eigenfunctions? , 2003 .

[16]  H. Schaefer Methods of Electronic Structure Theory , 1977 .

[17]  J. Olsen,et al.  Erratum: “A full configuration interaction and coupled-cluster study of the potential-energy surfaces for the lowest singlet excited state of N2” [J. Chem. Phys. 113, 6677 (2000)] , 2001 .

[18]  Rodney J. Bartlett,et al.  Full configuration-interaction and state of the art correlation calculations on water in a valence double-zeta basis with polarization functions , 1996 .

[19]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[20]  T. Daniel Crawford,et al.  Locally correlated equation-of-motion coupled cluster theory for the excited states of large molecules , 2002 .

[21]  S. Kucharski,et al.  Method of Moments of Coupled-Cluster Equations: A New Theoretical Framework for Designing "Black-Box" Approaches for Molecular Potential Energy Surfaces , 2002 .

[22]  M. Hoffmann,et al.  Low-lying potential energy surfaces , 2002 .

[23]  A. Szabó,et al.  Modern quantum chemistry : introduction to advanced electronic structure theory , 1982 .

[24]  Trygve Helgaker,et al.  Multiple basis sets in calculations of triples corrections in coupled-cluster theory , 1997 .

[25]  Peter Pulay,et al.  The local correlation treatment. II. Implementation and tests , 1988 .

[26]  Anna I. Krylov,et al.  Second order perturbation corrections to singles and doubles coupled-cluster methods: General theory and application to the valence optimized doubles model , 2000 .

[27]  Anna I. Krylov,et al.  Size-consistent wave functions for bond-breaking: the equation-of-motion spin-flip model , 2001 .

[28]  C. David Sherrill,et al.  Chapter 4 Bond Breaking in Quantum Chemistry , 2005 .

[29]  Jeppe Olsen,et al.  Determinant based configuration interaction algorithms for complete and restricted configuration interaction spaces , 1988 .

[30]  Anna I. Krylov,et al.  The spin–flip approach within time-dependent density functional theory: Theory and applications to diradicals , 2003 .

[31]  R. Bartlett,et al.  A direct product decomposition approach for symmetry exploitation in many-body methods. I. Energy calculations , 1991 .

[32]  J. Paldus,et al.  Reduced multireference couple cluster method. II. Application to potential energy surfaces of HF, F2, and H2O , 1998 .

[33]  Anna I. Krylov,et al.  Perturbative corrections to the equation-of-motion spin–flip self-consistent field model: Application to bond-breaking and equilibrium properties of diradicals , 2002 .

[34]  C. David Sherrill,et al.  Full configuration interaction potential energy curves for the X 1Σg+, B 1Δg, and B′ 1Σg+ states of C2: A challenge for approximate methods , 2004 .

[35]  Rodney J. Bartlett,et al.  Coupled-cluster methods with internal and semi-internal triply and quadruply excited clusters: CCSDt and CCSDtq approaches , 1999 .