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

Abstract

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.

References in corpus (6)

• The density-matrix renormalization group in the age of matrix product states
• Low-energy spectrum of iron-sulfur clusters directly from many-particle quantum mechanics
• The density matrix renormalization group for ab initio quantum chemistry
• A spin-adapted Density Matrix Renormalization Group algorithm for quantum chemistry
• Spectroscopic accuracy directly from quantum chemistry: application to ground and excited states of beryllium dimer
• Variational determination of the second-order density matrix for the isoelectronic series of beryllium, neon and silicon

Cited by in corpus (2)

• Multi-reference perturbation theory with Cholesky decomposition for the density matrix renormalization group
• A nonorthogonal state-interaction approach for matrix product state wave functions