Orbital-optimized third-order Møller-Plesset perturbation theory and its spin-component and spin-opposite scaled variants: application to symmetry breaking problems.

In this research, orbital-optimized third-order Møller-Plesset perturbation theory (OMP3) and its spin-component and spin-opposite scaled variants (SCS-OMP3 and SOS-OMP3) are introduced. Using a Lagrangian-based approach, an efficient, quadratically convergent algorithm for variational optimization of the molecular orbitals (MOs) for third-order Møller-Plesset perturbation theory (MP3) is presented. Explicit equations for response density matrices, the MO gradient, and Hessian are reported in spin-orbital form. The OMP3, SCS-OMP3, and SOS-OMP3 approaches are compared with the second-order Møller-Plesset perturbation theory (MP2), MP3, coupled-cluster doubles (CCD), optimized-doubles (OD), and coupled-cluster singles and doubles (CCSD) methods. All these methods are applied to the O(4)(+), O(3), and seven diatomic molecules. Results demonstrate that the OMP3 and its variants provide significantly better vibrational frequencies than MP3, CCSD, and OD for the molecules where the symmetry-breaking problems are observed. For O(4)(+), the OMP3 prediction, 1343 cm(-1), for ω(6) (b(3u)) mode, where symmetry-breaking appears, is even better than presumably more reliable methods such as Brueckner doubles (BD), 1194 cm(-1), and OD, 1193 cm(-1), methods (the experimental value is 1320 cm(-1)). For O(3), the predictions of SCS-OMP3 (1143 cm(-1)) and SOS-OMP3 (1165 cm(-1)) are remarkably better than the more robust OD method (1282 cm(-1)); the experimental value is 1089 cm(-1). For the seven diatomics, again the SCS-OMP3 and SOS-OMP3 methods provide the lowest average errors, ∣Δω(e)∣ = 44 and ∣Δω(e)∣ = 35 cm(-1), respectively, while for OD, ∣Δω(e)∣ = 161 cm(-1)and CCSD ∣Δω(e)∣ = 106 cm(-1). Hence, the OMP3 and especially its spin-scaled variants perform much better than the MP3, CCSD, and more robust OD approaches for considered test cases. Therefore, considering both the computational cost and the reliability, SCS-OMP3 and SOS-OMP3 appear to be the best methods for the symmetry-breaking cases, based on present application results. The OMP3 method offers certain advantages: it provides reliable vibrational frequencies in case of symmetry-breaking problems, especially with spin-scaling tricks, its analytic gradients are easier to compute since there is no need to solve the coupled-perturbed equations for the orbital response, and the computation of one-electron properties are easier because there is no response contribution to the particle density matrices. The OMP3 has further advantages over standard MP3, making it promising for excited state properties via linear response theory.

[1]  Matthew L. Leininger,et al.  PSI3: An open‐source Ab Initio electronic structure package , 2007, J. Comput. Chem..

[2]  M. Jacox,et al.  Evidence for the stabilization of rectangular O+4 in solid neon , 1994 .

[3]  H. Monkhorst,et al.  Calculation of properties with the coupled-cluster method , 2009 .

[4]  Henry F. Schaefer,et al.  Large multiconfiguration self-consistent-field wave functions for the ozone molecule , 1981 .

[5]  Rodney J. Bartlett,et al.  Similarity transformed equation-of-motion coupled-cluster theory: Details, examples, and comparisons , 1997 .

[6]  Poul Jo,et al.  Transition moments and dynamic polarizabilities in a second order polarization propagator approach , 1980 .

[7]  K. Andersson,et al.  Vibrational frequencies of ozone : a multiconfigurational approach , 1992 .

[8]  J. Olsen,et al.  Orbital-optimized coupled-cluster theory does not reproduce the full configuration-interaction limit. , 2005, The Journal of chemical physics.

[9]  M. Jacox,et al.  The vibrational spectra of molecular ions isolated in solid neon. II. O+4 and O−4 , 1989 .

[10]  R. Lindh,et al.  Symmetry breaking in O+4: an application of the Brueckner coupled-cluster method , 1994 .

[11]  S. Grimme,et al.  Is spin-component scaled second-order Møller-Plesset perturbation theory an appropriate method for the study of noncovalent interactions in molecules? , 2007, The journal of physical chemistry. A.

[12]  Thomas Bondo Pedersen,et al.  Gauge invariant coupled cluster response theory using optimized nonorthogonal orbitals , 2001 .

[13]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

[14]  A. Barbe,et al.  Infrared spectra of 16O3 and 18O3: Darling and Dennison resonance and anharmonic potential function of ozone , 1974 .

[15]  M. Head‐Gordon,et al.  Approaching closed-shell accuracy for radicals using coupled cluster theory with perturbative triple substitutions , 2003 .

[16]  Frank Neese,et al.  Assessment of Orbital-Optimized, Spin-Component Scaled Second-Order Many-Body Perturbation Theory for Thermochemistry and Kinetics. , 2009, Journal of chemical theory and computation.

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

[18]  G. Scuseria,et al.  The optimization of molecular orbitals for coupled cluster wavefunctions , 1987 .

[19]  Trygve Helgaker,et al.  Mo/ller–Plesset energy derivatives , 1988 .

[20]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[21]  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 .

[22]  Krishnan Raghavachari,et al.  Highly correlated systems: Structure, binding energy and harmonic vibrational frequencies of ozone , 1989 .

[23]  Poul Jo,et al.  Optimization of orbitals for multiconfigurational reference states , 1978 .

[24]  E. Davidson,et al.  Symmetry breaking in polyatomic molecules: real and artifactual , 1983 .

[25]  H. Schaefer,et al.  C3+ is bent , 1990 .

[26]  Takehiko Tanaka,et al.  Coriolis interaction and anharmonic potential function of ozone from the microwave spectra in the excited vibrational states , 1970 .

[27]  William A. Goddard,et al.  Configuration interaction studies of O3 and O+3. Ground and excited states , 1975 .

[28]  S. Grimme.  Improved second-order Møller–Plesset perturbation theory by separate scaling of parallel- and antiparallel-spin pair correlation energies , 2003 .

[29]  Poul Jørgensen,et al.  The second-order approximate coupled cluster singles and doubles model CC2 , 1995 .

[30]  Ivan Hubač,et al.  Size-extensivity correction for the state-specific multireference Brillouin–Wigner coupled-cluster theory , 2000 .

[31]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[32]  Josef Paldus,et al.  Simultaneous handling of dynamical and nondynamical correlation via reduced multireference coupled cluster method: Geometry and harmonic force field of ozone , 1999 .

[33]  Henry F. Schaefer,et al.  On the evaluation of analytic energy derivatives for correlated wave functions , 1984 .

[34]  John F. Stanton,et al.  Many-body perturbation theory with a restricted open-shell Hartree—Fock reference , 1991 .

[35]  Stefan Grimme,et al.  Improved third‐order Møller–Plesset perturbation theory , 2003, J. Comput. Chem..

[36]  Poul Jørgensen,et al.  Response functions from Fourier component variational perturbation theory applied to a time-averaged quasienergy , 1998 .

[37]  W. D. Allen,et al.  Hartree-Fock orbital instability envelopes in highly correlated single-reference wave functions , 1997 .

[38]  M. Head‐Gordon,et al.  Violations of N-representability from spin-unrestricted orbitals in Møller–Plesset perturbation theory and related double-hybrid density functional theory , 2009 .

[39]  W. D. Allen,et al.  The analytic evaluation of energy first derivatives for two‐configuration self‐consistent‐field configuration interaction (TCSCF‐CI) wave functions. Application to ozone and ethylene , 1987 .

[40]  A. Varandas.  Spin-component-scaling second-order Møller-Plesset theory and its variants for economical correlation energies: unified theoretical interpretation and use for quartet N3. , 2010, The Journal of chemical physics.

[41]  Henry F. Schaefer,et al.  Ordering of the O-O stretching vibrational frequencies in ozone , 1989 .

[42]  K. Kowalski.  Implementation of the locally renormalized CCSD(T) approaches for arbitrary reference function. , 2005, The Journal of chemical physics.

[43]  Trygve Helgaker,et al.  Analytical Calculation of Geometrical Derivatives in Molecular Electronic Structure Theory , 1988 .

[44]  Frank Neese,et al.  Correlated ab initio spin densities for larger molecules: orbital-optimized spin-component-scaled MP2 method. , 2010, The journal of physical chemistry. A.

[45]  Henrik Koch,et al.  Coupled cluster response functions , 1990 .

[46]  Ron Shepard,et al.  The Analytic Gradient Method for Configuration Interaction Wave Functions , 1995 .

[47]  W. D. Allen,et al.  Is the oxywater radical cation more stable than neutral oxywater , 1996 .

[48]  Anna I. Krylov,et al.  Excited states theory for optimized orbitals and valence optimized orbitals coupled-cluster doubles models , 2000 .

[49]  Ove Christiansen,et al.  Response functions in the CC3 iterative triple excitation model , 1995 .

[50]  Anna I. Krylov,et al.  Size-consistent wave functions for nondynamical correlation energy: The valence active space optimized orbital coupled-cluster doubles model , 1998 .

[51]  James A Platts,et al.  Spin-Component Scaling Methods for Weak and Stacking Interactions. , 2007, Journal of chemical theory and computation.

[52]  Martin Head-Gordon,et al.  Optimization of wave function and geometry in the finite basis Hartree-Fock method , 1988 .

[53]  R. Bartlett,et al.  The equilibrium structure and harmonic vibrational frequencies of ozone: coupled cluster results including triple excitations , 1989 .

[54]  Garnet Kin-Lic Chan,et al.  Tailored coupled cluster singles and doubles method applied to calculations on molecular structure and harmonic vibrational frequencies of ozone. , 2006, The Journal of chemical physics.

[55]  Jürgen Gauss,et al.  Calculation of excited-state properties using general coupled-cluster and configuration-interaction models. , 2004, The Journal of chemical physics.

[56]  Martin Head-Gordon,et al.  The perfect quadruples model for electron correlation in a valence active space. , 2009, The Journal of chemical physics.

[57]  Martin Head-Gordon,et al.  Scaled opposite-spin second order Møller-Plesset correlation energy: an economical electronic structure method. , 2004, The Journal of chemical physics.

[58]  David H. Magers,et al.  Highly correlated single‐reference studies of the O3 potential surface. I. Effects of high order excitations on the equilibrium structure and harmonic force field of ozone , 1989 .

[59]  Francesco A Evangelista,et al.  Coupling term derivation and general implementation of state-specific multireference coupled cluster theories. , 2007, The Journal of chemical physics.

[60]  Rudolph C. Mayrhofer,et al.  Complete active space self‐consistent field potential energy surfaces, dipole moment functions, and spectroscopic properties of O3, CF2, NO−2, and NF+2 , 1991 .

[61]  Björn O. Roos,et al.  The Multiconfigurational (MC) SCF Method , 1983 .

[62]  R. Bartlett,et al.  Connected quadruples for the frequencies of O3 , 1999 .

[63]  C. Sherrill.  Computations of Noncovalent π Interactions , 2009 .

[64]  Rodney J. Bartlett,et al.  Analytic energy derivatives in many‐body methods. I. First derivatives , 1989 .

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

[66]  Alistair P. Rendell,et al.  The analytic configuration interaction gradient method: The calculation of one electron properties , 1987 .

[67]  Martin Head-Gordon,et al.  The Theoretical Prediction of Molecular Radical Species: a Systematic Study of Equilibrium Geometries and Harmonic Vibrational Frequencies , 2001 .

[68]  Trygve Helgaker,et al.  The CC3 model: An iterative coupled cluster approach including connected triples , 1997 .

[69]  Thom H. Dunning,et al.  Gaussian basis sets for use in correlated molecular calculations. V. Core-valence basis sets for boron through neon , 1995 .

[70]  R. Fink.  Spin-component-scaled Møller-Plesset (SCS-MP) perturbation theory: a generalization of the MP approach with improved properties. , 2010, The Journal of chemical physics.

[71]  H. Schaefer,et al.  Molecular geometry and vibrational frequencies of ozone from compact variational wave functions explicitly including triple and quadruple substitutions , 1997 .

[72]  J. Linderberg,et al.  State vectors and propagators in many‐electron theory. A unified approach , 1977 .

[73]  Anna I. Krylov,et al.  Energies and analytic gradients for a coupled-cluster doubles model using variational Brueckner orbitals: Application to symmetry breaking in O4+ , 1998 .

[74]  Gustavo E. Scuseria,et al.  The vibrational frequencies of ozone , 1990 .

[75]  Martin Head-Gordon,et al.  A tractable and accurate electronic structure method for static correlations: the perfect hextuples model. , 2010, The Journal of chemical physics.

[76]  Rohini C. Lochan,et al.  Orbital-optimized opposite-spin scaled second-order correlation: an economical method to improve the description of open-shell molecules. , 2007, The Journal of chemical physics.

[77]  Trygve Helgaker,et al.  Excitation energies from the coupled cluster singles and doubles linear response function (CCSDLR). Applications to Be, CH+, CO, and H2O , 1990 .

[78]  S. Grimme,et al.  Spin-component scaled second-order Møller–Plesset perturbation theory for the calculation of molecular geometries and harmonic vibrational frequencies , 2004 .

[79]  John F. Stanton,et al.  A benchmark coupled-cluster single, double, and triple excitation (CCSDT) study of the structure and harmonic vibrational frequencies of the ozone molecule☆ , 1991 .

[80]  John F. Stanton,et al.  Restricted open-shell Hartree-Fock-based many-body perturbation theory: Theory and application of energy and gradient calculations , 1992 .

[81]  D. A. Horner,et al.  The lithium superoxide radical: Symmetry breaking phenomena and potential energy surfaces , 1989 .

[82]  R. Lindh,et al.  The fraternal twins of quartet O+4 , 1994 .

[83]  N. Handy,et al.  Large basis set calculations using Brueckner theory , 1994 .

[84]  Michael T. Heath,et al.  Scientific Computing: An Introductory Survey , 1996 .

[85]  Hua Guo,et al.  Accurate ab initio near-equilibrium potential energy and dipole moment functions of the ground electronic state of ozone , 2000 .

[86]  P. Jeffrey Hay,et al.  Geometries and energies of the excited states of O3 from ab initio potential energy surfaces , 1977 .

[87]  Charles W. Bauschlicher,et al.  Theoretical calculation of ozone vibrational infrared intensities , 1985 .

[88]  H. Koch,et al.  Gauge invariant coupled cluster response theory , 1999 .

[89]  A. D. McLean,et al.  Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z=11–18 , 1980 .

[90]  R. Bartlett,et al.  On the choice of orbitals for symmetry breaking problems with application to NO3 , 1992 .

[91]  Hans-Joachim Werner,et al.  Calculation of intermolecular interactions in the benzene dimer using coupled-cluster and local electron correlation methods. , 2006, Physical chemistry chemical physics : PCCP.

[92]  R. Bartlett,et al.  COUPLED-CLUSTER CALCULATIONS OF STRUCTURE AND VIBRATIONAL FREQUENCIES OF OZONE : ARE TRIPLE EXCITATIONS ENOUGH? , 1998 .

[93]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[94]  Uğur Bozkaya,et al.  Quadratically convergent algorithm for orbital optimization in the orbital-optimized coupled-cluster doubles method and in orbital-optimized second-order Møller-Plesset perturbation theory. , 2011, The Journal of chemical physics.

[95]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .