We present a quantum-classical algorithm for electronic structure calculations that dramatically reduces the quantum measurement cost of variational quantum eigensolver (VQE) approaches. While conventional VQE methods require measurements scaling as O(M^4) with system size M, the proposed occupation-number VQE (ON-VQE) reduces this cost to O(M/2) by avoiding reduced density matrix (RDM) measurements and relying exclusively on ONs. The method exploits only the diagonal elements of the one-particle RDM in the natural orbital representation, where occupations are obtained directly from computational-basis measurement outcomes. By restricting the variational ansatz to double excitations within orbital subspaces associated with electron pairs, the required measurements can be grouped into a small number of qubit-wise commuting observables, yielding an efficient and scalable measurement strategy. The approach is validated through simulations and executions on quantum hardware for the cubic H$_8$ cluster, demonstrating the feasibility of extracting accurate ONs from quantum measurements and evaluating electronic energies within the natural orbital functional (NOF) framework. Across representative molecular systems, the extracted ONs enable accurate energy evaluation with state-of-the-art NOFs while maintaining a dramatically reduced measurement cost. These results establish a scalable route toward quantum simulation of strongly correlated electronic systems, demonstrating that accurate electronic energies can be obtained from quantum measurements of ONs alone.
We present a quantum-classical algorithm for electronic structure calculations that dramatically reduces the quantum measurement cost of variational quantum eigensolver (VQE) approaches. While conventional VQE methods require measurements scaling as O(M^4) with system size M, the proposed occupation-number VQE (ON-VQE) reduces this cost to O(M/2) by avoiding reduced density matrix (RDM) measurements and relying exclusively on ONs. The method exploits only the diagonal elements of the one-particle RDM in the natural orbital representation, where occupations are obtained directly from computational-basis measurement outcomes. By restricting the variational ansatz to double excitations within orbital subspaces associated with electron pairs, the required measurements can be grouped into a small number of qubit-wise commuting observables, yielding an efficient and scalable measurement strategy. The approach is validated through simulations and executions on quantum hardware for the cubic H8 cluster, demonstrating the feasibility of extracting accurate ONs from quantum measurements and evaluating electronic energies within the natural orbital functional (NOF) framework. Across representative molecular systems, the extracted ONs enable accurate energy evaluation with state-of-the-art NOFs while maintaining a dramatically reduced measurement cost. These results establish a scalable route toward quantum simulation of strongly correlated electronic systems, demonstrating that accurate electronic energies can be obtained from quantum measurements of ONs alone.
In this work, we explore the use of the one-particle reduced density matrix (1RDM) to streamline energy measurements of chemical systems on quantum computers, particularly within the variational quantum eigensolver (VQE) framework. This approach leverages the existence of an exact energy functional of the 1RDM, enabling a reduction in both the number of expectation values to measure and the number of circuits to execute, thereby optimizing quantum resource usage. Specifically, sampling the 1RDM involves measuring only elements, which is significantly fewer than the required for the Hamiltonian’s expectation value ⟨ Ĥ ⟩. We demonstrate our approach by harnessing the well-established natural orbital functional (NOF) theory, using the natural orbitals and occupation numbers derived from the diagonalization of the 1RDM measured from the quantum computer. Starting with the H 2 system, we validate the accuracy of our method by comparing the energy derived from NOF approximations applied to the exact wave function with the value obtained from ⟨ Ĥ ⟩. This is followed by an optimization of the gate parameters by minimizing the energy using the NOF approximations as the objective function. The analysis is extended to larger systems, such as LiH, Li 2, OH –, FH, NeH +, and F 2 using a wave function ansatz with single and double excitation gates. This NOF-based method reduces the scaling cost of circuit executions compared to standard VQE implementations, achieving around 90% savings in the systems used in this work. Overall, by using a well-performing NOF as the objective function, the proposed NOF-VQE demonstrates the viability of NOF approximations for obtaining accurate energies in the noisy intermediate-scale quantum era and underscores the potential for developing new functionals tailored to quantum computing applications.
Abstract Variational quantum eigensolvers (VQE) are among the most promising approaches for solving electronic structure problems on near-term quantum computers. A critical challenge for VQE in practice is that one needs to strike a balance between the expressivity of the VQE ansatz versus the number of quantum gates required to implement the ansatz, given the reality of noisy quantum operations on near-term quantum computers. In this work, we consider an orbital-optimized pair-correlated approximation to the unitary coupled cluster with singles and doubles (uCCSD) ansatz and report a highly efficient quantum circuit implementation for trapped-ion architectures. We show that orbital optimization can recover significant additional electron correlation energy without sacrificing efficiency through measurements of low-order reduced density matrices (RDMs). In the dissociation of small molecules, the method gives qualitatively accurate predictions in the strongly-correlated regime when running on noise-free quantum simulators. On IonQ’s Harmony and Aria trapped-ion quantum computers, we run end-to-end VQE algorithms with up to 12 qubits and 72 variational parameters—the largest full VQE simulation with a correlated wave function on quantum hardware. We find that even without error mitigation techniques, the predicted relative energies across different molecular geometries are in excellent agreement with noise-free simulators.
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.
ABSTRACT 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 a measurement‐reduced orbital‐optimized (OO) 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.
Improving the Accuracy of the Variational Quantum Eigensolver for Molecular Systems by the Explicitly-Correlated Perturbative [2]R12-Correction
We provide an integration of the universal, perturbative explicitly correlated [2]R12-correction in the context of the Variational Quantum Eigensolver (VQE). This approach is able to increase the accuracy of the underlying reference method significantly while requiring no additional quantum resources. Our proposed approach only requires knowledge of the one- and two-particle reduced density matrices (RDMs) of the reference wavefunction; these can be measured after having reached convergence in VQE. The RDMs are then combined with a set of molecular integrals. This computation comes at a cost that scales as the sixth power of the number of electrons. We explore the performance of the VQE+[2]R12 approach using both conventional Gaussian basis sets and our recently proposed directly determined pair-natural orbitals obtained by multiresolution analysis (MRA-PNOs). Both Gaussian orbital and PNOs are investigated as a potential set of complementary basis functions in the computation of [2]R12. In particular the combination of MRA-PNOs with [2]R12 has turned out to be very promising -- persistently throughout our data, this allowed very accurate simulations at a quantum cost of a minimal basis set. Additionally, we found that the deployment of PNOs as complementary basis can greatly reduce the number of complementary basis functions that enter the computation of the correction at a cubic complexity.
Accurately predicting the properties of molecules and materials from first principles requires solving the electronic structure problem. However, the computational cost of solving the electronic structure problem scales prohibitively in the number of atoms, and the usual approximate methods designed to reduce the overhead are often inaccurate for strongly correlated systems. For this reason, quantum algorithms have been developed as a promising alternative. In this work, we develop algebraic tools to improve the efficiency of three quantum algorithms for the electronic structure problem, addressing key bottlenecks for both near-term and fault-tolerant quantum computers. The first is the variational quantum eigensolver (VQE), which suffers from a large measurement overhead needed to estimate the Hamiltonian expectation value. To this end, we introduce a general class of measurable fragments of the electronic Hamiltonian, which enable partitioning schemes that significantly reduce the sampling cost for energy estimation. The second are quantum subspace methods (QSMs). We introduce a new QSM, the quantum seniority-based subspace expansion algorithm, which is a scalable alternative to VQE that alleviates VQE's state preparation, measurement, and optimization challenges while reducing qubit requirements. The third is quantum phase estimation (QPE), which requires an encoding of the electronic Hamiltonian directly on the quantum computer, the cost of which scales indirectly with the Hamiltonian spectral range. To address this, we develop linear programming approaches to implement the Block Invariant Symmetry Shift protocol, a method which reduces the spectral range of the Hamiltonian, improving its block-encoding cost for QPE. Collectively, these methods demonstrate that one can exploit the Hamiltonian's algebraic structure to mitigate computational costs, accelerating the path toward practical quantum advantage for electronic structure.