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.
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.
Variational quantum eigensolver (VQE) is a hybrid quantum-classical algorithm designed for noisy intermediate-scale quantum (NISQ) computers. It is promising for quantum chemical calculations (QCC) because it can calculate the ground-state energy of molecules. Although VQE has a theoretical potential to achieve a higher accuracy than classical approximation methods, it is challenging to achieve it on current NISQ computers due to the significant impact of noise. Density matrix embedding theory (DMET) is a well-known technique that can be used to decompose a molecule into multiple fragments, which is available to mitigate the noise impact on VQE. However, our preliminary evaluation shows that the naive combination of DMET and VQE hardly mitigates it. In this work, we present two proposals to mitigate the noise impact on the DMET+VQE combination. (1) The size of quantum circuits used by VQE is decreased by reducing the number of bath orbitals which represent interactions between multiple fragments in DMET. (2) While VQE is executed on a NISQ computer, reduced density matrices (RDMs), which are used to calculate the energy of fragments in DMET, are calculated accurately by executing quantum circuits using a noise-less quantum computer simulator. The evaluation using a noisy quantum computer simulator shows that our proposals improve the accuracy of the DMET+VQE combination by an order of magnitude.
Variational quantum eigensolver (VQE) is a hybrid quantum-classical algorithm designed for noisy intermediate-scale quantum (NISQ) computers. It is promising for quantum chemical calculations (QCC) because it can calculate the ground-state energy of a target molecule. Although VQE has a potential to achieve a higher accuracy than classical approximation methods in QCC, it is challenging to achieve it on current NISQ computers due to the significant impact of noises. Density matrix embedding theory (DMET) is a well-known technique to divide a molecule into multiple fragments, which is available to mitigate the noise impact on VQE. However, our preliminary evaluation shows that the naive combination of DMET and VQE does not outperform a gold standard classical method. In this work, we present three approaches to mitigate the noise impact for the DMET+VQE combination. (1) The size of quantum circuits used by VQE is decreased by reducing the number of bath orbitals which represent interactions between multiple fragments in DMET. (2) Reduced density matrices (RDMs), which are used to calculate a molecular energy in DMET, are calculated accurately based on expectation values obtained by executing quantum circuits using a noise-less quantum computer simulator. (3) The parameters of a quantum circuit optimized by VQE are refined with mathematical post-processing. The evaluation using a noisy quantum computer simulator shows that our approaches significantly improve the accuracy of the DMET+VQE combination. Moreover, we demonstrate that on a real NISQ device, the DMET+VQE combination applying our three approaches achieves a higher accuracy than the gold standard classical method.
Quantum hardware noise remains a major obstacle to accurate molecular simulations using the variational quantum eigensolver (VQE). We introduce a novel error-mitigation strategy based on the extended random phase approximation (ERPA), which reconstructs the two-particle reduced density matrix (2-RDM) from approximate transition density matrices, effectively filtering errors in quantum measurements. Unlike quantum subspace expansion, the approach requires only the one- and two-particle RDMs. Since ERPA preserves the input 1-RDM, we further propose a hybrid quantum-classical scheme employing an inexpensive density matrix renormalization group (DMRG) calculation to provide a noise-free approximation to this quantity. Benchmarks on water double-bond stretching and cyclobutadiene automerization, performed on IBM quantum hardware, demonstrate reductions in energy errors by orders of magnitude. Combined with multireference adiabatic connection (AC0) corrections, the approach additionally recovers dynamical correlation beyond the active space, yielding accurate molecular potential energy curves and reaction barriers. Our results establish ERPA as a promising, measurement-efficient route toward reliable quantum-classical electronic structure calculations.
The variational quantum eigensolver (VQE) has emerged as one of the leading quantum algorithms for solving electronic structure problems on near-term noisy intermediate-scale quantum devices. However, its practical application to quantum chemistry remains challenging due to the limited coherence time, imperfect quantum gate fidelity, and the large number of measurements required, which together confine current electronic structure simulations to relatively small active spaces. In this work, we present an orbital-optimized VQE framework based on orbital compression, designed to improve the accuracy of electronic structure calculations while maintaining relatively small active spaces. Frozen natural orbitals (FNO) and split virtual orbitals (SVO) are first employed to construct compact active spaces for VQE simulations, leading to the FNO/SVO-VQE approach. Orbital optimization is then incorporated to further recover electron correlation effects, resulting in the FNO/SVO-OO-VQE methods. We apply the proposed method to simulate potential energy surfaces for molecular dissociation and the activation energy of formaldehyde decomposition. Numerical results demonstrate that both FNO-OO-VQE and SVO-OO-VQE improve the variational accuracy while substantially reducing measurement cost.
The variational quantum eigensolver (VQE) is a hybrid quantum-classical algorithm designed for noisy intermediate-scale quantum (NISQ) hardware to estimate eigenvalues of many-body Hamiltonians. Unlike fully quantum approaches such as quantum phase estimation (QPE), VQE trades deep coherent circuits for repeated state preparation, measurement, and classical optimization, making it more compatible with limited qubit counts and finite coherence times. Recent developments have focused on reducing quantum resource requirements while retaining chemically meaningful wavefunction structure. In this paper, we examine recent progress in VQE methods for molecular electronic structure with an emphasis on three themes: (i) strategies for circuit and ansatz complexity reduction, including adaptive and selectively screened approaches, (ii) chemically motivated workflows that combine VQE with orbital optimization, fragmentation, and localized active-space ideas to better address strong correlation, and (iii) extensions of VQE to excited-state calculations. Throughout, we emphasize the tradeoffs among parameter count, gate depth, symmetry preservation, measurement overhead, and classical preprocessing, and discuss where these approaches may become most useful for chemically challenging active spaces. We also highlight benchmarking considerations for assessing both accuracy and resource requirements, and conclude with a perspective on regimes in which VQE may offer the greatest long-term value, particularly multireference active spaces and low-lying excited-state manifolds.
Recent practical approaches for the use of current generation noisy quantum devices in the simulation of quantum many-body problems have been dominated by the use of a variational quantum eigensolver (VQE). These coupled quantum-classical algorithms leverage the ability to perform many repeated measurements to avoid the currently prohibitive gate depths often required for exact quantum algorithms, with the restriction of a parametrized circuit to describe the states of interest. In this work, we show how the calculation of zero-temperature dynamic correlation functions defining the linear response characteristics of quantum systems can also be recast into a modified VQE algorithm, which can be incorporated into the current variational quantum infrastructure. This allows for these important physical expectation values describing the dynamics of the system to be directly converged on the frequency axis, and they approach exactness over all frequencies as the flexibility of the parametrization increases. The frequency resolution hence does not explicitly scale with gate depth, which is approximately twice as deep as a ground-state VQE. We apply the method to compute the single-particle Green's function of ab initio dihydrogen and lithium hydride molecules, and demonstrate the use of a practical active space embedding approach to extend to larger systems. While currently limited by the fidelity of two-qubit gates, whose number is increased compared to the ground-state algorithm on current devices, we believe the approach shows potential for the extraction of frequency dynamics of correlated systems on near-term quantum processors.