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

Physical Chemistry Chemical Physics

Statistical errors in reduced density matrices sampled from quantum circuit simulation and the impact on multireference perturbation theory

Soichiro Nishio, Yuki Oba, Yuki Kurashige

Published Jan 1, 2023• 5 citations
Publisher

Abstract

In this work, we present a detailed analysis of statistical errors in reduced density matrices (RDMs) of active space wavefunctions sampled from quantum circuit simulation and the impact on results obtained by the multireference theories.
DOI
10.1039/d3cp03520d

Related papers

Multireference error mitigation for quantum computation of chemistry

Hang Zou, Erika Magnusson, Hampus Brunander +2 more

MREM extends reference-state error mitigation by incorporating multireference wavefunctions, enabling accurate and noise-resilient quantum simulations of strongly correlated molecules.

arXiv:2505.08291•Digital Discovery•Jan 1, 2025•5 citations

Multireference configuration interaction and perturbation theory without reduced density matrices

Ankit Mahajan, Nick S. Blunt, Iliya Sabzevari +1 more

The computationally expensive evaluation and storage of high-rank reduced density matrices (RDMs) has been the bottleneck in the calculation of dynamic correlation for multireference wave functions in large active spaces. We present a stochastic formulation of multireference configuration interaction and perturbation theory that avoids the need for these expensive RDMs. The algorithm presented here is flexible enough to incorporate a wide variety of active space reference wave functions, including selected configuration interaction, matrix product states, and symmetry-projected Jastrow mean field wave functions. It enjoys the usual attractive features of Monte Carlo methods, such as embarrassing parallelizability and low memory costs. We find that the stochastic algorithm is already competitive with the deterministic algorithm for small active spaces, containing as few as 14 orbitals. We illustrate the utility of our stochastic formulation using benchmark applications.

arXiv:1909.06935•The Journal of Chemical Physics•Dec 2, 2019•30 citations

Density-matrix spectra for two-dimensional quantum systems

Ming-Chiang Chung, Ingo Peschel

For a two-dimensional system of coupled oscillators, the spectra of reduced density matrices can be obtained analytically. This provides an example where the features of these quantities, which are of central importance in numerical studies using the density-matrix renormalization-group method, can be seen

arXiv:cond-mat/0004222•Physical review. B, Condensed matter•Aug 15, 2000•93 citations

Perturbative Variational Quantum Eigensolver via Reduced Density Matrices

Yuhan Zheng, Yibin Guo, Huili Zhang +5 more

Current noisy intermediate-scale quantum (NISQ) devices remain limited in their ability to perform accurate quantum chemistry simulations because of restricted numbers of high-fidelity qubits and short coherence times. To overcome these challenges, we introduce a reduced-density-matrix (RDM)-based perturbative variational quantum eigensolver (VQE) framework that augments active-space VQE with perturbation theory to recover electron correlation beyond the active space without increasing the qubit count or variational circuit depth. We formulated a fully coupled approach (VQE-PTs) and a diagonal approximation (VQE-PT). The former retains couplings among orthonormalized perturbers, whereas the latter neglects these couplings to simplify the classical postprocessing. Numerical simulations of HF, N2, and F2 show that VQE-PTs provides a robust formulation across different molecular systems, while VQE-PT offers an efficient approximation. We further experimentally implement VQE-PT on the Quafu superconducting quantum processor for F2, achieving a mean absolute error of 1.2 mHa along the potential energy surface after error mitigation. These results demonstrate perturbative VQE as a practical framework for incorporating dynamic correlations in quantum chemistry simulations.

arXiv:2504.02340•The Journal of Physical Chemistry A•Oct 7, 2026•0 citations

The Rigorous Stochastic Matrix Product Method for the Calculations of Reduced Equilibrium Density Matrices of Multilevel Open Quantum Systems

Xin Chen

•Oct 28, 2013•0 citations

Combining Internally Contracted States and Matrix Product States To Perform Multireference Perturbation Theory

Sandeep Sharma, Gerald Knizia, Sheng Guo +1 more

We present two efficient and intruder-free methods for treating dynamic correlation on top of general multiconfiguration reference wave functions — including such as obtained by the density matrix renormalization group (DMRG) with large active spaces. The new methods are the second order variant of the recently proposed multireference linearized coupled cluster method (MRLCC) [Sharma, S.; Alavi, A. J. Chem. Phys. 2015, 143, 102815] and of N-electron valence perturbation theory (NEVPT2), with expected accuracies similar to MRCI+Q and (at least) CASPT2, respectively. Great efficiency gains are realized by representing the first order wave function with a combination of internal contraction (IC) and matrix product state perturbation theory (MPSPT). With this combination, only third order reduced density matrices (RDMs) are required. Thus, we obviate the need for calculating (or estimating) RDMs of fourth or higher order; these had so far posed a severe bottleneck for dynamic correlation treatments involving the large active spaces accessible to DMRG. Using several benchmark systems, including first and second row containing small molecules, Cr 2, pentacene, and oxo-Mn(Salen), we show that active spaces containing at least 30 orbitals can be treated using this method. On a single node, MRLCC2 and NEVPT2 calculations can be performed with over 550 and 1100 virtual orbitals, respectively. We also critically examine the errors incurred due to the three sources of errors introduced in the present implementation – calculating second order instead of third order energy corrections, use of internal contraction, and approximations made in the reference wave function due to DMRG.

arXiv:1609.03496•Journal of Chemical Theory and Computation•Jan 7, 2017•67 citations
The computationally expensive evaluation and storage of high-rank reduced density matrices (RDMs) has been the bottleneck in the calculation of dynamic correlation for multireference wave functions in large active spaces. We present a stochastic formulation of multireference configuration interaction and perturbation theory that avoids the need for these expensive RDMs. The algorithm presented here is flexible enough to incorporate a wide variety of active space reference wave functions, including selected configuration interaction, matrix product states, and symmetry-projected Jastrow mean field wave functions. It enjoys the usual attractive features of Monte Carlo methods, such as embarrassing parallelizability and low memory costs. We find that the stochastic algorithm is already competitive with the deterministic algorithm for small active spaces, containing as few as 14 orbitals. We illustrate the utility of our stochastic formulation using benchmark applications.
Current noisy intermediate-scale quantum (NISQ) devices remain limited in their ability to perform accurate quantum chemistry simulations because of restricted numbers of high-fidelity qubits and short coherence times. To overcome these challenges, we introduce a reduced-density-matrix (RDM)-based perturbative variational quantum eigensolver (VQE) framework that augments active-space VQE with perturbation theory to recover electron correlation beyond the active space without increasing the qubit count or variational circuit depth. We formulated a fully coupled approach (VQE-PTs) and a diagonal approximation (VQE-PT). The former retains couplings among orthonormalized perturbers, whereas the latter neglects these couplings to simplify the classical postprocessing. Numerical simulations of HF, N2, and F2 show that VQE-PTs provides a robust formulation across different molecular systems, while VQE-PT offers an efficient approximation. We further experimentally implement VQE-PT on the Quafu superconducting quantum processor for F2, achieving a mean absolute error of 1.2 mHa along the potential energy surface after error mitigation. These results demonstrate perturbative VQE as a practical framework for incorporating dynamic correlations in quantum chemistry simulations.
We present two efficient and intruder-free methods for treating dynamic correlation on top of general multiconfiguration reference wave functions — including such as obtained by the density matrix renormalization group (DMRG) with large active spaces. The new methods are the second order variant of the recently proposed multireference linearized coupled cluster method (MRLCC) [Sharma, S.; Alavi, A. J. Chem. Phys. 2015, 143, 102815] and of N-electron valence perturbation theory (NEVPT2), with expected accuracies similar to MRCI+Q and (at least) CASPT2, respectively. Great efficiency gains are realized by representing the first order wave function with a combination of internal contraction (IC) and matrix product state perturbation theory (MPSPT). With this combination, only third order reduced density matrices (RDMs) are required. Thus, we obviate the need for calculating (or estimating) RDMs of fourth or higher order; these had so far posed a severe bottleneck for dynamic correlation treatments involving the large active spaces accessible to DMRG. Using several benchmark systems, including first and second row containing small molecules, Cr 2, pentacene, and oxo-Mn(Salen), we show that active spaces containing at least 30 orbitals can be treated using this method. On a single node, MRLCC2 and NEVPT2 calculations can be performed with over 550 and 1100 virtual orbitals, respectively. We also critically examine the errors incurred due to the three sources of errors introduced in the present implementation – calculating second order instead of third order energy corrections, use of internal contraction, and approximations made in the reference wave function due to DMRG.