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:2405.15422

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, Stefan Knecht

Published May 24, 2024• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.2405.15422

Related papers

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

Michał Krompiec, David Muñoz Ramo

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.

arXiv:2210.05702•Oct 11, 2022•10 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

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

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

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

Jan-Niklas Boyn, Лыхин Александр Олегович, Scott E. Smart +2 more

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. We obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

arXiv:2106.11972•The Journal of Chemical Physics•Dec 28, 2021•32 citations
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.
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.
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.
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 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.
While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. We obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.