arXiv:1703.10830 · The Journal of Chemical Physics

Time-dependent N-electron valence perturbation theory with matrix product state reference wavefunctions for large active spaces and basis sets: Applications to the chromium dimer and all-trans polyenes

Alexander Yu. Sokolov, Sheng Guo, Enrico Ronca, Garnet Kin-Lic Chan

Published Jun 22, 2017• 66 citations

Abstract

In earlier work [A. Y. Sokolov and G. K.-L. Chan, J. Chem. Phys. 144, 064102 (2016)], we introduced a time-dependent formulation of the second-order N-electron valence perturbation theory (t-NEVPT2) which (i) had a lower computational scaling than the usual internally contracted perturbation formulation and (ii) yielded the fully uncontracted NEVPT2 energy. Here, we present a combination of t-NEVPT2 with a matrix product state (MPS) reference wavefunction (t-MPS-NEVPT2) that allows us to compute uncontracted dynamic correlation energies for large active spaces and basis sets, using the time-dependent density matrix renormalization group algorithm. In addition, we report a low-scaling MPS-based implementation of strongly contracted NEVPT2 (sc-MPS-NEVPT2) that avoids computation of the four-particle reduced density matrix. We use these new methods to compute the dissociation energy of the chromium dimer and to study the low-lying excited states in all-trans polyenes (C4H6 to C24H26), incorporating dynamic correlation for reference wavefunctions with up to 24 active electrons and orbitals.

DOI: 10.1063/1.4986975

Related papers

Time-dependent and strongly-contracted N-electron valence perturbation theories with matrix product state reference wave functions for large active spaces and basis sets: Applications to the chromium dimer and all-trans polyenes

Alexander Yu. Sokolov, Sheng Guo, Enrico Ronca +1 more

In earlier work [J. Chem. Phys. 144, 064102 (2016)], we introduced a time-dependent formulation of the second-order N-electron valence perturbation theory (t-NEVPT2) which (i) had a lower computational scaling than the usual internally-contracted perturbation formulation, and (ii) yielded the fully uncontracted NEVPT2 energy. Here, we present a combination of t-NEVPT2 with a matrix product state (MPS) reference wavefunction (t-MPS-NEVPT2) that allows to compute uncontracted dynamic correlation energies for large active spaces and basis sets, using the time-dependent density matrix renormalization group (td-DMRG) algorithm. In addition, we report a low-scaling MPS-based implementation of strongly-contracted NEVPT2 (sc-MPS-NEVPT2) that avoids computation of the four-particle reduced density matrix. We use these new methods to compute the dissociation energy of the chromium dimer and to study the low-lying excited states in all-trans polyenes (C_4H_6 to C_(24)H_(26)), incorporating dynamic correlation for reference wavefunctions with up to 24 active electrons and orbitals.

Mar 31, 20170 citations

N-Electron Valence State Perturbation Theory Based on a Density Matrix Renormalization Group Reference Function, with Applications to the Chromium Dimer and a Trimer Model of Poly(p-Phenylenevinylene)

Sheng Guo, Mark A. Watson, Weifeng Hu +2 more

The strongly contracted variant of second-order N-electron valence state perturbation theory (NEVPT2) is an efficient perturbative method to treat dynamic correlation without the problems of intruder states or level shifts, while the density matrix renormalization group (DMRG) provides the capability to address static correlation in large active spaces. We present a combination of the DMRG and strongly contracted NEVPT2 (DMRG-SC-NEVPT2) that uses an efficient algorithm to compute high-order reduced-density matrices from DMRG wave functions. The capabilities of DMRG-SC-NEVPT2 are demonstrated on calculations of the chromium dimer potential energy curve at the basis set limit, and the excitation energies of a trimer model of poly(p-phenylenevinylene) (PPV(n = 3)).

arXiv:1512.08137Journal of Chemical Theory and ComputationFeb 25, 2016187 citations

Multi-reference electron correlation with large active spaces from time-dependent perturbation theory with density matrix renormalization group

Alexander Yu. Sokolov, Garnet Kin‐Lic Chan

In the recent years there has been a remarkable progress in the development of methods for the description of strong electron correlation in mols. There are now many methods that describe strong correlation in active spaces with up to 60 orbitals. However, capturing the remaining weak (dynamic) correlation outside of the active space remains a challenge. Conventional multi-ref. methods for dynamic correlation (such as internally-contracted perturbation theory, e.g. CASPT2 or NEVPT2) require computation of up to fourth-order reduced d. matrixes, which prohibits their application to systems with active spaces with more than 24 \norbitals. In this talk I will describe a new method which provides an efficient description of static and dynamic \nelectron correlation by combining the formalism of d. matrix renormalization group (DMRG) and timedependent \nmulti-ref. perturbation theory. This approach has a no. of important advantages:. 1) It is equiv. to fully uncontracted perturbation theory, but has a lower computational scaling than the contracted approaches;. 2) It avoids computation of the three- and four-particle reduced d. matrixes;. 3) It has a polynomial scaling with the no. of active orbitals and can be applied to active spaces with much more than 24 orbitals. In my talk, I will present a brief overview of multi-ref. time-dependent perturbation theory, discuss its implementation in combination with DMRG, and demonstrate its performance for strongly correlated systems with large active spaces.

Apr 1, 20170 citations

Approximations of density matrices in N-electron valence state second-order perturbation theory (NEVPT2). III. Large active space calculations with selected configuration interaction reference

Yang Guo, Kantharuban Sivalingam, Vijay Gopal Chilkuri +1 more

In multi-reference (MR) methods, addressing systems with large active spaces remains a challenge in the field. In Papers I and II of this series, we demonstrated that full rank N-electron valence state second-order perturbation theory (FR-NEVPT2) is a robust MR perturbation theory capable of computing strongly correlated systems with approximate density matrices. However, the previous FR-NEVPT2 implementation requires the computation and storage of fifth-order reduced density matrices (RDMs), limiting the usage of FR-NEVPT2 for systems with large active spaces. In the present work, as Paper III of the series, we report a new FR-NEVPT2 algorithm to handle systems with large active spaces. In the new algorithm, an approximate complete active space (CAS) self-consistent field (SCF) method, iterative configuration expansion (ICE) SCF, is employed to compute the reference wave functions for FR-NEVPT2. Then, the necessary Koopmans matrices of FR-NEVPT2 involving various RDMs are constructed using the intermediates designed by Kollmar et al. [J. Chem. Phys. 155, 234104 (2021)] to avoid storage bottlenecks. The performance of the new FR-NEVPT2 algorithm for systems with large active spaces is evaluated. Our results show that even with aggressive truncation parameters to truncate the ICE-SCF reference wave function, FR-NEVPT2 effectively recovers the missing static correlations of ICE-SCF. Several interesting systems with active spaces up to CAS(34,34) are studied using FR-NEVPT2 with ICE-SCF reference.

The Journal of Chemical PhysicsApr 11, 20257 citations

DMRG-CASPT2 study of the longitudinal static second hyperpolarizability of all-trans polyenes

Sebastian Wouters, Véronique Van Speybroeck, Dimitri Van Neck

We have implemented internally contracted complete active space second order perturbation theory (CASPT2) with the density matrix renormalization group (DMRG) as active space solver [Y. Kurashige and T. Yanai, J. Chem. Phys. 135, 094104 (2011)]. Internally contracted CASPT2 requires to contract the generalized Fock matrix with the 4-particle reduced density matrix (4-RDM) of the reference wavefunction. The required 4-RDM elements can be obtained from 3-particle reduced density matrices (3-RDM) of different wavefunctions, formed by symmetry-conserving single-particle excitations op top of the reference wavefunction. In our spin-adapted DMRG code chemps2 https://github.com/sebwouters/chemps2, we decompose these excited wavefunctions as spin-adapted matrix product states and calculate their 3-RDM in order to obtain the required contraction of the generalized Fock matrix with the 4-RDM of the reference wavefunction. In this work, we study the longitudinal static second hyperpolarizability of all-trans polyenes C2nH2n+2 [n = 4-12] in the cc-pVDZ basis set. DMRG-SCF and DMRG-CASPT2 yield substantially lower values and scaling with system size compared to RHF and MP2, respectively.

arXiv:1605.05526The Journal of Chemical PhysicsAug 5, 201675 citations

A projected approximation to strongly contracted N-electron valence perturbation theory for DMRG wavefunctions

Michael Roemelt, Sheng Guo, Garnet Kin‐Lic Chan

A novel approach to strongly contracted N-electron valence perturbation theory (SC-NEVPT2) as a means of describing dynamic electron correlation for quantum chemical density matrix renormalization group (DMRG) calculations is presented. In this approach the strongly contracted perturber functions are projected onto a renormalized Hilbert space. Compared to a straightforward implementation of SC-NEVPT2 with DMRG wavefunctions, the computational scaling and storage requirements are reduced. This favorable scaling opens up the possibility of calculations with larger active spaces. A specially designed renormalization scheme ensures that both the electronic ground state and the perturber functions are well represented in the renormalized Hilbert space. Test calculations on the N2 and [Cu2O2(en)2](2+) demonstrate some key properties of the method and indicate its capabilities.

The Journal of Chemical PhysicsMay 26, 201658 citations