Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2210.05702

Strongly Contracted N-Electron Valence State Perturbation Theory Using Reduced Density Matrices from a Quantum Computer

Michał Krompiec, David Muñoz Ramo

Published Oct 11, 2022• 10 citations
arXivHTMLPDF

Abstract

We introduce QRDM-NEVPT2: a hybrid quantum-classical implementation of strongly-contracted N-electron Valence State $2^{nd}$-order Perturbation Theory (SC-NEVPT2), in which the Complete Active Space Configuration Interaction (CASCI) step, capturing static correlation effects, is replaced by a simulation performed on a quantum computer. Subsequently, n-particle Reduced Density Matrices (n-RDMs) measured on a quantum device are used directly in a classical SC-NEVPT2 calculation, which recovers remaining dynamic electron correlation effects approximately. We also discuss the use of the cumulant expansion to approximate the whole 4-RDM matrix or only its zeros. In addition to noiseless state-vector quantum simulations, we demonstrate, for the first time, a hybrid quantum-classical multi-reference perturbation theory calculation, with the quantum part performed on a quantum computer.

DOI
10.48550/arxiv.2210.05702

Related papers

Quantum-centric strong and dynamical electron correlation: A resource-efficient second-order $N$-electron valence perturbation theory formulation for near-term quantum devices

Aaron Fitzpatrick, N. Walter Talarico, Roberto Di Remigio +1 more

We present a measurement-cost efficient implementation of Strongly-Contracted $N$-Electron Valence Perturbation Theory (SC-NEVPT2) for use on near-term quantum devices. At the heart of our algorithm we exploit the properties of adaptive Informationally Complete positive operator valued measures (IC-POVMs) to recycle the measurement outcomes from a ground state energy estimation on a quantum device to reconstruct the matrix elements of the three- and four-body reduced density matrices for use in a subsequent CPU-driven NEVPT2 calculation. The proposed scheme is capable of producing results in good agreement with corresponding conventional NEVPT2 simulations, while significantly reducing the cost of quantum measurements and allowing for embarrassingly parallel estimations of higher-order RDMs in classical post-processing. Our scheme shows favourable scaling of the total number of shots with respect to system size. This paves the way for routine inclusion of dynamic electron correlation effects in hybrid quantum-classical computing pipelines.

arXiv:2405.15422•May 24, 2024•0 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.08137•Journal of Chemical Theory and Computation•Feb 25, 2016•187 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 Physics•Apr 11, 2025•7 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 Physics•May 26, 2016•58 citations

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, 2017•0 citations

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 +1 more

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.

arXiv:1703.10830•The Journal of Chemical Physics•Jun 22, 2017•66 citations
We introduce QRDM-NEVPT2: a hybrid quantum-classical implementation of strongly-contracted N-electron Valence State 2nd2^{nd}2nd-order Perturbation Theory (SC-NEVPT2), in which the Complete Active Space Configuration Interaction (CASCI) step, capturing static correlation effects, is replaced by a simulation performed on a quantum computer. Subsequently, n-particle Reduced Density Matrices (n-RDMs) measured on a quantum device are used directly in a classical SC-NEVPT2 calculation, which recovers remaining dynamic electron correlation effects approximately. We also discuss the use of the cumulant expansion to approximate the whole 4-RDM matrix or only its zeros. In addition to noiseless state-vector quantum simulations, we demonstrate, for the first time, a hybrid quantum-classical multi-reference perturbation theory calculation, with the quantum part performed on a quantum computer.
Quantum-centric strong and dynamical electron correlation: A resource-efficient second-order NNN-electron valence perturbation theory formulation for near-term quantum devices
We present a measurement-cost efficient implementation of Strongly-Contracted NNN-Electron Valence Perturbation Theory (SC-NEVPT2) for use on near-term quantum devices. At the heart of our algorithm we exploit the properties of adaptive Informationally Complete positive operator valued measures (IC-POVMs) to recycle the measurement outcomes from a ground state energy estimation on a quantum device to reconstruct the matrix elements of the three- and four-body reduced density matrices for use in a subsequent CPU-driven NEVPT2 calculation. The proposed scheme is capable of producing results in good agreement with corresponding conventional NEVPT2 simulations, while significantly reducing the cost of quantum measurements and allowing for embarrassingly parallel estimations of higher-order RDMs in classical post-processing. Our scheme shows favourable scaling of the total number of shots with respect to system size. This paves the way for routine inclusion of dynamic electron correlation effects in hybrid quantum-classical computing pipelines.
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)).
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.
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.
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.
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.