Communication: A flexible multi-reference perturbation theory by minimizing the Hylleraas functional with matrix product states.

We describe a formulation of multi-reference perturbation theory that obtains a rigorous upper bound to the second order energy by minimizing the Hylleraas functional in the space of matrix product states (MPS). The first order wavefunctions so obtained can also be used to compute the third order energy with little overhead. Our formulation has several advantages including (i) flexibility with respect to the choice of zeroth order Hamiltonian, (ii) recovery of the exact uncontracted multi-reference perturbation theory energies in the limit of large MPS bond dimension, (iii) no requirement to compute high body density matrices, (iv) an embarrassingly parallel algorithm (scaling up to the number of virtual orbitals, squared, processors). Preliminary numerical examples show that the MPS bond dimension required for accurate first order wavefunctions scales sub-linearly with the size of the basis.

[1]  T. Yanai,et al.  High-performance ab initio density matrix renormalization group method: applicability to large-scale multireference problems for metal compounds. , 2009, The Journal of chemical physics.

[2]  George H. Booth,et al.  Spectroscopic accuracy directly from quantum chemistry: application to ground and excited states of beryllium dimer. , 2014, The Journal of chemical physics.

[3]  Garnet Kin-Lic Chan,et al.  An algorithm for large scale density matrix renormalization group calculations. , 2004, The Journal of chemical physics.

[4]  G. Chan,et al.  Entangled quantum electronic wavefunctions of the Mn₄CaO₅ cluster in photosystem II. , 2013, Nature chemistry.

[5]  Bjoern O. Roos.  Theoretical Studies of Electronically Excited States of Molecular Systems Using Multiconfigurational Perturbation Theory , 1999 .

[6]  E. Hylleraas.  Über den Grundterm der Zweielektronenprobleme von H−, He, Li+, Be++ usw. , 1930 .

[7]  T. Yanai,et al.  Canonical transcorrelated theory with projected Slater-type geminals. , 2012, The Journal of chemical physics.

[8]  Takeshi Yanai,et al.  Multireference configuration interaction theory using cumulant reconstruction with internal contraction of density matrix renormalization group wave function. , 2013, The Journal of chemical physics.

[9]  Ericka Stricklin-Parker,et al.  Ann , 2005 .

[10]  Roman Orus,et al.  A Practical Introduction to Tensor Networks: Matrix Product States and Projected Entangled Pair States , 2013, 1306.2164.

[11]  Sandeep Sharma,et al.  The density matrix renormalization group in quantum chemistry. , 2011, Annual review of physical chemistry.

[12]  Edward F. Valeev,et al.  Explicitly correlated R12/F12 methods for electronic structure. , 2012, Chemical reviews.

[13]  Ali Alavi,et al.  An explicitly correlated approach to basis set incompleteness in full configuration interaction quantum Monte Carlo. , 2012, The Journal of chemical physics.

[14]  Alistair P. Rendell,et al.  The restricted active space self-consistent-field method, implemented with a split graph unitary group approach , 1990 .

[15]  Markus Reiher,et al.  New electron correlation theories for transition metal chemistry. , 2011, Physical chemistry chemical physics : PCCP.

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

[17]  Seiichiro Ten-no,et al.  Explicitly correlated wave functions: summary and perspective , 2012, Theoretical Chemistry Accounts.

[18]  Xiang.  Density-matrix renormalization-group method in momentum space. , 1996, Physical review. B, Condensed matter.

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

[20]  Ali Alavi,et al.  Linear-scaling and parallelisable algorithms for stochastic quantum chemistry , 2013, 1305.6981.

[21]  Frank Neese,et al.  Low-energy spectrum of iron-sulfur clusters directly from many-particle quantum mechanics. , 2014, Nature chemistry.

[22]  Celestino Angeli,et al.  Introduction of n-electron valence states for multireference perturbation theory , 2001 .

[23]  J. Sólyom,et al.  Applications of Quantum Information in the Density-Matrix Renormalization Group , 2008 .

[24]  Ali Alavi,et al.  Fermion Monte Carlo without fixed nodes: a game of life, death, and annihilation in Slater determinant space. , 2009, The Journal of chemical physics.

[25]  Kenneth G. Dyall,et al.  The choice of a zeroth-order Hamiltonian for second-order perturbation theory with a complete active space self-consistent-field reference function , 1995 .

[26]  O. Sǐnanoğlu.  VARIATION-PERTURBATION METHOD FOR EXCITED STATES , 1961 .

[27]  Garnet Kin-Lic Chan,et al.  Low entanglement wavefunctions , 2012 .

[28]  David W. Small,et al.  Post-modern valence bond theory for strongly correlated electron spins. , 2011, Physical chemistry chemical physics : PCCP.

[29]  Kerstin Andersson,et al.  Second-order perturbation theory with a CASSCF reference function , 1990 .

[30]  Frederick R Manby,et al.  Local explicitly correlated second-order perturbation theory for the accurate treatment of large molecules. , 2009, The Journal of chemical physics.

[31]  Helen van Aggelen,et al.  Variational determination of the second-order density matrix for the isoelectronic series of beryllium, neon, and silicon , 2009, 0909.0225.

[32]  F. Neese,et al.  Assessment of n-Electron Valence State Perturbation Theory for Vertical Excitation Energies. , 2013, Journal of chemical theory and computation.

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

[34]  Edward F. Valeev,et al.  Variational formulation of perturbative explicitly-correlated coupled-cluster methods. , 2008, Physical chemistry chemical physics : PCCP.

[35]  Legeza,et al.  Accuracy of the density-matrix renormalization-group method. , 1996, Physical review. B, Condensed matter.

[36]  K. Morokuma,et al.  Multireference Ab Initio Density Matrix Renormalization Group (DMRG)-CASSCF and DMRG-CASPT2 Study on the Photochromic Ring Opening of Spiropyran. , 2013, Journal of chemical theory and computation.

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

[38]  Takeshi Yanai,et al.  Second-order perturbation theory with a density matrix renormalization group self-consistent field reference function: theory and application to the study of chromium dimer. , 2011, The Journal of chemical physics.

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

[40]  U. Schollwoeck.  The density-matrix renormalization group in the age of matrix product states , 2010, 1008.3477.

[41]  Martin Schütz,et al.  Molpro: a general‐purpose quantum chemistry program package , 2012 .

[42]  Garnet Kin-Lic Chan,et al.  Spin-adapted density matrix renormalization group algorithms for quantum chemistry. , 2012, The Journal of chemical physics.

[43]  M. Head‐Gordon,et al.  Highly correlated calculations with a polynomial cost algorithm: A study of the density matrix renormalization group , 2002 .

[44]  Frederick R. Manby,et al.  R12 methods in explicitly correlated molecular electronic structure theory , 2006 .

[45]  Richard L. Martin,et al.  Ab initio quantum chemistry using the density matrix renormalization group , 1998 .

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

[47]  O. Sǐnanoğlu.  Relation of Perturbation Theory to Variation Method , 1961 .

[48]  Jeppe Olsen,et al.  The initial implementation and applications of a general active space coupled cluster method , 2000 .

[49]  Mihály Kállay,et al.  Higher excitations in coupled-cluster theory , 2001 .

[50]  Giovanni Li Manni,et al.  The generalized active space concept in multiconfigurational self-consistent field methods. , 2011, The Journal of chemical physics.