Orbital-optimized density cumulant functional theory.

In density cumulant functional theory (DCFT) the electronic energy is evaluated from the one-particle density matrix and two-particle density cumulant, circumventing the computation of the wavefunction. To achieve this, the one-particle density matrix is decomposed exactly into the mean-field (idempotent) and correlation components. While the latter can be entirely derived from the density cumulant, the former must be obtained by choosing a specific set of orbitals. In the original DCFT formulation [W. Kutzelnigg, J. Chem. Phys. 125, 171101 (2006)] the orbitals were determined by diagonalizing the effective Fock operator, which introduces partial orbital relaxation. Here we present a new orbital-optimized formulation of DCFT where the energy is variationally minimized with respect to orbital rotations. This introduces important energy contributions and significantly improves the description of the dynamic correlation. In addition, it greatly simplifies the computation of analytic gradients, for which expressions are also presented. We offer a perturbative analysis of the new orbital stationarity conditions and benchmark their performance for a variety of chemical systems.

[1]  M. Nooijen,et al.  State specific equation of motion coupled cluster method in general active space. , 2009, The Journal of chemical physics.

[2]  D. Mazziotti.  Approximate solution for electron correlation through the use of Schwinger probes , 1998 .

[3]  Wilfried Meyer,et al.  PNO–CI Studies of electron correlation effects. I. Configuration expansion by means of nonorthogonal orbitals, and application to the ground state and ionized states of methane , 1973 .

[4]  Werner Kutzelnigg,et al.  Error analysis and improvements of coupled-cluster theory , 1991 .

[5]  M. Head‐Gordon,et al.  Size-consistent Brueckner theory limited to double and triple substitutions , 1990 .

[6]  M. Head‐Gordon,et al.  A fifth-order perturbation comparison of electron correlation theories , 1989 .

[7]  Edward F. Valeev,et al.  Universal perturbative explicitly correlated basis set incompleteness correction. , 2009, The Journal of chemical physics.

[8]  C. Valdemoro,et al.  The contracted Schrödinger equation methodology: study of the third-order correlation effects , 2007 .

[9]  Edward F. Valeev,et al.  SF-[2]R12: a spin-adapted explicitly correlated method applicable to arbitrary electronic states. , 2011, The Journal of chemical physics.

[10]  D. Mazziotti.  Two-electron reduced density matrices from the anti-Hermitian contracted Schrodinger equation: enhanced energies and properties with larger basis sets. , 2007, The Journal of chemical physics.

[11]  John F. Stanton,et al.  Analytic energy gradients for open-shell coupled-cluster singles and doubles (CCSD) calculations using restricted open-shell Hartree—Fock (ROHF) reference functions , 1991 .

[12]  Garnet Kin-Lic Chan,et al.  Multireference quantum chemistry through a joint density matrix renormalization group and canonical transformation theory. , 2010, The Journal of chemical physics.

[13]  E. Davidson.  N-representability of the electron pair density , 1995 .

[14]  W. Kutzelnigg.  Unconventional Aspects of Coupled-Cluster Theory , 2010 .

[15]  J. Cizek.  On the Correlation Problem in Atomic and Molecular Systems. Calculation of Wavefunction Components in Ursell-Type Expansion Using Quantum-Field Theoretical Methods , 1966 .

[16]  H. Schaefer,et al.  Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method: improving upon CCSD(T) and CCSD(T)(Λ): preliminary application. , 2012, The Journal of chemical physics.

[17]  Liguo Kong.  Extension of the cumulant definition for low particle number systems , 2010 .

[18]  F. Neese,et al.  Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method. , 2009, The Journal of chemical physics.

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

[20]  Julia E. Rice,et al.  Analytic evaluation of energy gradients for the single and double excitation coupled cluster (CCSD) wave function: Theory and application , 1987 .

[21]  Nicholas C. Handy,et al.  Size-consistent Brueckner theory limited to double substitutions , 1989 .

[22]  D. Mazziotti.  Variational method for solving the contracted Schrödinger equation through a projection of the N-particle power method onto the two-particle space , 2002 .

[23]  Werner Kutzelnigg,et al.  Quantum chemistry in Fock space. I. The universal wave and energy operators , 1982 .

[24]  Uğur Bozkaya,et al.  Orbital-optimized third-order Møller-Plesset perturbation theory and its spin-component and spin-opposite scaled variants: application to symmetry breaking problems. , 2011, The Journal of chemical physics.

[25]  John M. Herbert,et al.  Magnitude and significance of the higher-order reduced density matrix cumulants , 2007 .

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

[27]  Andrew C Simmonett,et al.  Analytic gradients for density cumulant functional theory: the DCFT-06 model. , 2012, The Journal of chemical physics.

[28]  Uğur Bozkaya,et al.  Orbital-optimized coupled-electron pair theory and its analytic gradients: accurate equilibrium geometries, harmonic vibrational frequencies, and hydrogen transfer reactions. , 2013, The Journal of chemical physics.

[29]  W. Kutzelnigg.  Almost variational coupled cluster theory , 1998 .

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

[31]  Francesco A Evangelista,et al.  Alternative single-reference coupled cluster approaches for multireference problems: the simpler, the better. , 2011, The Journal of chemical physics.

[32]  Garnet Kin-Lic Chan,et al.  Canonical transformation theory for multireference problems. , 2006, The Journal of chemical physics.

[33]  C. Valdemoro,et al.  Approximating q-order reduced density matrices in terms of the lower-order ones. II. Applications. , 1993, Physical review. A, Atomic, molecular, and optical physics.

[34]  Garnet Kin-Lic Chan,et al.  Quadratic canonical transformation theory and higher order density matrices. , 2009, The Journal of chemical physics.

[35]  M. Hanauer,et al.  Meaning and magnitude of the reduced density matrix cumulants , 2012 .

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

[37]  John F. Stanton,et al.  The accurate determination of molecular equilibrium structures , 2001 .

[38]  A. DePrince,et al.  Parametric approach to variational two-electron reduced-density-matrix theory , 2007 .

[39]  Edward F. Valeev,et al.  A novel interpretation of reduced density matrix and cumulant for electronic structure theories. , 2011, The Journal of chemical physics.

[40]  T. Crawford,et al.  An Introduction to Coupled Cluster Theory for Computational Chemists , 2007 .

[41]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[42]  Rodney J. Bartlett,et al.  Alternative ansätze in single reference coupled-cluster theory. III. A critical analysis of different methods , 1995 .

[43]  Andrew G. Taube,et al.  Rethinking linearized coupled-cluster theory. , 2009, The Journal of chemical physics.

[44]  Edward F. Valeev,et al.  Perturbative correction for the basis set incompleteness error of complete-active-space self-consistent field. , 2010, The Journal of chemical physics.

[45]  F. Coester,et al.  Bound states of a many-particle system , 1958 .

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

[47]  Hans-Joachim Werner,et al.  A comparison of the efficiency and accuracy of the quadratic configuration interaction (QCISD), coupled cluster (CCSD), and Brueckner coupled cluster (BCCD) methods , 1992 .

[48]  M. Nooijen,et al.  A state-specific partially internally contracted multireference coupled cluster approach. , 2011, The Journal of chemical physics.

[49]  R. Bartlett.  Many-Body Perturbation Theory and Coupled Cluster Theory for Electron Correlation in Molecules , 1981 .

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

[51]  Debashis Mukherjee,et al.  Normal order and extended Wick theorem for a multiconfiguration reference wave function , 1997 .

[52]  G. Scuseria.  An alternative expression for the perturbative triples correction to Brueckner coupled-cluster doubles theory , 1994 .

[53]  D. Mazziotti.  Anti-hermitian part of the contracted Schrödinger equation for the direct calculation of two-electron reduced density matrices , 2007 .

[54]  Andreas Köhn,et al.  Communication: Restoring full size extensivity in internally contracted multireference coupled cluster theory. , 2012, The Journal of chemical physics.

[55]  U. Bozkaya.  Analytic energy gradients for the orbital-optimized third-order Møller-Plesset perturbation theory. , 2013, The Journal of chemical physics.

[56]  Toru Shiozaki,et al.  Explicitly correlated multireference configuration interaction: MRCI-F12. , 2011, The Journal of chemical physics.

[57]  Yasuda,et al.  Direct determination of the quantum-mechanical density matrix using the density equation. , 1996, Physical review letters.

[58]  P. Szalay,et al.  Analytic energy derivatives for coupled‐cluster methods describing excited states: General formulas and comparison of computational costs , 1995 .

[59]  Rodney J. Bartlett,et al.  Many‐body perturbation theory, coupled‐pair many‐electron theory, and the importance of quadruple excitations for the correlation problem , 1978 .

[60]  J. Percus,et al.  Reduction of the N‐Particle Variational Problem , 1964 .

[61]  Werner Kutzelnigg,et al.  Generalized k-particle brillouin conditions and their use for the construction of correlated electronic wavefunctions , 1979 .

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

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

[64]  D. Klein.  Degenerate perturbation theory , 1974 .

[65]  W. Kutzelnigg,et al.  Irreducible Brillouin conditions and contracted Schrödinger equations for n-electron systems. IV. Perturbative analysis. , 2004, The Journal of chemical physics.

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

[67]  V. H. Smith,et al.  On Different Criteria for the Best Independent‐Particle Model Approximation , 1964 .

[68]  Werner Kutzelnigg,et al.  Density-cumulant functional theory. , 2006, The Journal of chemical physics.

[69]  Matthew L. Leininger,et al.  Psi4: an open‐source ab initio electronic structure program , 2012 .

[70]  R. Nesbet.  Brueckner's Theory and the Method of Superposition of Configurations , 1958 .

[71]  D. Cramb,et al.  The (2,5) vibrational band of the A—X system of N+2 observed using an infrared diode laser velocity modulation spectrometer , 1991 .

[72]  Bridgette Cooper,et al.  Benchmark studies of variational, unitary and extended coupled cluster methods. , 2010, The Journal of chemical physics.

[73]  R. Bartlett,et al.  Coupled‐cluster open‐shell analytic gradients: Implementation of the direct product decomposition approach in energy gradient calculations , 1991 .

[74]  F. Weinhold,et al.  Reduced Density Matrices of Atoms and Molecules. II. On the N‐Representability Problem , 1967 .

[75]  David A Mazziotti,et al.  Anti-hermitian contracted schrödinger equation: direct determination of the two-electron reduced density matrices of many-electron molecules. , 2006, Physical review letters.

[76]  F. Coester,et al.  Short-range correlations in nuclear wave functions , 1960 .

[77]  Debashis Mukherjee,et al.  Cumulant expansion of the reduced density matrices , 1999 .

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

[79]  C. Kollmar.  A size extensive energy functional derived from a double configuration interaction approach: the role of N representability conditions. , 2006, The Journal of chemical physics.

[80]  Andrew C Simmonett,et al.  Density cumulant functional theory: the DC-12 method, an improved description of the one-particle density matrix. , 2013, The Journal of chemical physics.

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

[82]  J. Stanton.  Why CCSD(T) works: a different perspective , 1997 .

[83]  K. Brueckner,et al.  NUCLEAR SATURATION AND TWO-BODY FORCES. II. TENSOR FORCES , 1954 .

[84]  R. Bartlett,et al.  Coupled-cluster theory in quantum chemistry , 2007 .

[85]  F. Neese,et al.  A comparative study of single reference correlation methods of the coupled-pair type , 2008 .

[86]  W. Kutzelnigg.  New derivation and a k‐particle generalization of SCF‐type theories , 1980 .

[87]  David A. Mazziotti,et al.  Parametrization of the two-electron reduced density matrix for its direct calculation without the many-electron wave function: Generalizations and applications , 2010 .

[88]  G. Winnewisser,et al.  Torsional Dependence of the Effective Rotational Constants of H2O2 and H2S2 , 1993 .

[89]  Debashis Mukherjee,et al.  Irreducible Brillouin conditions and contracted Schrödinger equations for n-electron systems. I. The equations satisfied by the density cumulants , 2001 .

[90]  D. Mazziotti.  Contracted Schrödinger equation: Determining quantum energies and two-particle density matrices without wave functions , 1998 .

[91]  Garnet Kin-Lic Chan,et al.  Canonical transformation theory from extended normal ordering. , 2007, The Journal of chemical physics.

[92]  Garnet Kin-Lic Chan,et al.  A review of canonical transformation theory , 2010 .

[93]  D. Mazziotti.  Two-electron reduced density matrix as the basic variable in many-electron quantum chemistry and physics. , 2012, Chemical reviews.

[94]  H. Primas.  GENERALIZED PERTURBATION THEORY IN OPERATOR FORM , 1963 .

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

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

[97]  Garnet Kin-Lic Chan,et al.  Strongly contracted canonical transformation theory. , 2010, The Journal of chemical physics.

[98]  Bernard Kirtman,et al.  Simultaneous calculation of several interacting electronic states by generalized Van Vleck perturbation theory , 1981 .

[99]  A. J. Coleman.  THE STRUCTURE OF FERMION DENSITY MATRICES , 1963 .

[100]  David A Mazziotti,et al.  Parametrization of the two-electron reduced density matrix for its direct calculation without the many-electron wave function. , 2008, Physical review letters.

[101]  Frank Neese,et al.  The coupled electron pair approximation: variational formulation and spin adaptation , 2010 .

[102]  R. Shepard.  The Multiconfiguration Self‐Consistent Field Method , 2007 .

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

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

[105]  Jeremiah J. Wilke,et al.  Density cumulant functional theory: first implementation and benchmark results for the DCFT-06 model. , 2010, The Journal of chemical physics.