Simulation of Potential Energy Curves for H₂ and LiH Molecules Using the Variational Quantum Eigensolver: A Comparative Study of Ideal and Noisy Quantum Hardware Conditions on Qiskit
This study investigates the Variational Quantum Eigensolver (VQE), a hybrid quantum-classical algorithm that partitions the computational workload between a quantum processing unit and a classical optimiser to determine the ground-state energies of molecular systems. The Qiskit SDK v1.x framework is employed, utilising the PySCFDriver module within the qiskit-nature package to construct second-quantised molecular electronic Hamiltonians for two representative systems: Hydrogen ( ) and Lithium Hydride ( ). For the complete active space is retained without spatial reduction, mapping the full Hamiltonian onto a 4-qubit register. However for the full STO-3G representation demands 12 qubits, rendering direct classical parameter optimisation computationally intractable. This bottleneck is resolved by applying an ActiveSpaceTransformer (num_electrons=2, num_spatial_orbitals=3) to freeze the chemically inactive lithium core electrons, reducing the active register to 6 qubits. Fermionic operators are subsequently converted to Pauli strings via the Jordan-Wigner transformation ( ) enabling the quantum processor to evaluate the molecular Hamiltonian as a circuit while preserving the fermionic anti-commutation relations. Under ideal, noise-free simulation using the StatevectorEstimator primitive paired with the COBYLA optimiser (maxiter=200), all VQE energy points converge precisely onto the Full Configuration Interaction (FCI) classical baseline across the full range of interatomic separations. This confirms that the hardware-efficient TwoLocal ansatz ( rotation layers and entanglement layers with reps=3) possesses sufficient expressive power to capture the molecular ground-state wavefunction. Introducing gate-level noise through AerSimulator backends calibrated from the FakeManilaV2 and FakeNairobiV2 hardware profiles—with optimisation driven by the SPSA algorithm (maxiter=150)—reveals two critical failure modes. A systematic upward energy bias is observed, attributable to depolarising gate errors and measurement readout faults that prevent the optimiser from reaching the true minimum energy eigenvalue. Additionally, severe stochastic oscillations between 1.0 Å and 1.5 Å on the curve induce premature convergence of the SPSA optimiser to spurious local minima. These results demonstrate that unmitigated hardware noise renders VQE incapable of reliably predicting equilibrium bond lengths on current NISQ devices, motivating the future integration of Zero-Noise Extrapolation (ZNE), readout error mitigation, and Virtual Distillation into the computational pipeline. Keywords—Variational Quantum Eigensolver (VQE); Qiskit 1.x; Potential Energy Curve; Molecular Ground-State Energy; Hardware-Efficient Ansatz; Gate-Level Noise Simulation; ActiveSpaceTransformer; NISQ Devices.
This study investigates the Variational Quantum Eigensolver (VQE), a hybrid quantum-classical algorithm that partitions the computational workload between a quantum processing unit and a classical optimiser to determine the ground-state energies of molecular systems. The Qiskit SDK v1.x framework is employed, utilising the PySCFDriver module within the qiskit-nature package to construct second-quantised molecular electronic Hamiltonians for two representative systems: Hydrogen ( ) and Lithium Hydride ( ). For the complete active space is retained without spatial reduction, mapping the full Hamiltonian onto a 4-qubit register. However for the full STO-3G representation demands 12 qubits, rendering direct classical parameter optimisation computationally intractable. This bottleneck is resolved by applying an ActiveSpaceTransformer (num_electrons=2, num_spatial_orbitals=3) to freeze the chemically inactive lithium core electrons, reducing the active register to 6 qubits. Fermionic operators are subsequently converted to Pauli strings via the Jordan-Wigner transformation ( ) enabling the quantum processor to evaluate the molecular Hamiltonian as a circuit while preserving the fermionic anti-commutation relations. Under ideal, noise-free simulation using the StatevectorEstimator primitive paired with the COBYLA optimiser (maxiter=200), all VQE energy points converge precisely onto the Full Configuration Interaction (FCI) classical baseline across the full range of interatomic separations. This confirms that the hardware-efficient TwoLocal ansatz ( rotation layers and entanglement layers with reps=3) possesses sufficient expressive power to capture the molecular ground-state wavefunction. Introducing gate-level noise through AerSimulator backends calibrated from the FakeManilaV2 and FakeNairobiV2 hardware profiles—with optimisation driven by the SPSA algorithm (maxiter=150)—reveals two critical failure modes. A systematic upward energy bias is observed, attributable to depolarising gate errors and measurement readout faults that prevent the optimiser from reaching the true minimum energy eigenvalue. Additionally, severe stochastic oscillations between 1.0 Å and 1.5 Å on the curve induce premature convergence of the SPSA optimiser to spurious local minima. These results demonstrate that unmitigated hardware noise renders VQE incapable of reliably predicting equilibrium bond lengths on current NISQ devices, motivating the future integration of Zero-Noise Extrapolation (ZNE), readout error mitigation, and Virtual Distillation into the computational pipeline. Keywords—Variational Quantum Eigensolver (VQE); Qiskit 1.x; Potential Energy Curve; Molecular Ground-State Energy; Hardware-Efficient Ansatz; Gate-Level Noise Simulation; ActiveSpaceTransformer; NISQ Devices.
molecule, emphasizing practical implementation and performance on current quantum hardware. Our research presents a comparative study of HEA and UCCSD ansätze on noiseless and noisy simulations and implements VQE on recent IBM quantum computer noise models and a real quantum computer, IBM Fez, providing a fully functional code employing Qiskit 1.2. Our experiments confirm UCCSD's reliability in ideal conditions, while the HEA demonstrates greater robustness to hardware noise, achieving chemical accuracy on state-vector simulation (SVS). The results reveal that achieving ground-state energy within chemical accuracy is feasible without error mitigation during VQE convergence. We demonstrate that current quantum devices effectively optimize circuit parameters despite misestimating simulated energies. The SVS-evaluated energies provide a more accurate representation of the solution quality compared to QPU-estimated energy values, indicating that VQE converges to the correct ground state despite quantum noise. Our study also applies noise mitigation as a postprocessing technique, using zero-noise extrapolation (ZNE) on a real quantum computer. The detailed methodologies presented in this study, including Hamiltonian construction and Fermionic-to-qubit transformations, facilitate flexible adaptation of the VQE approach for various algorithm variants and across different levels of algorithmic implementation.
Classical simulation of molecular electronic structure is often limited by the rapid increase in computational cost associated with describing electron correlation. The Variational Quantum Eigensolver (VQE) is one hybrid quantum-classical approach ideally adapted to the Noisy Intermediate-Scale Quantum (NISQ) regime. This work applies the VQE algorithm to compute ground state energies of a set of molecules:H2,H3+,LiH,HF,OH−,H2O,BH3,CH4, andNH3. The approach utilizes the Unitary Coupled Cluster for Single and Double excitations (UCCSD) ansatz along with the Bravyi-Kitaev transformation for fermionic-to-qubit encoding. Quantum simulation outputs are systematically benchmarked against Full Configuration Interaction (FCI), a numerically exact traditional approach, and Unrestricted Hartree-Fock (UHF). The results establish a robust consistency of the VQE-obtained molecular ground state energies against standard FCI and UHF benchmarks. This research provides a critical benchmark, confirming that VQE has the capability for high-precision molecular simulations and playing the role of a pivotal stepping stone toward revolutionizing applications in drug discovery and materials science.
While numerical simulations are presented in most papers introducing new methods to enhance the VQE performance, comprehensive, comparative, and applied studies remain relatively rare. We present a comprehensive, yet concise guide for the implementation of the VQE for molecular problems on NISQ devices, specifically applied to estimate the ground state energy of the BeH2 molecule using hardware-efficient and chemically informed ansätze. This work clarifies several under-documented aspects in the literature, such as the construction of the electronic Hamiltonian, the transformation of fermionic operators into qubit operators via second quantization, and the mathematical framework's details for the unitary coupled cluster single and double (UCCSD) ansatz. Our methodology, implemented using Qiskit 1.2, the latest release as of the date of this writing, is demonstrated on a noiseless simulator and further tested with noisy quantum circuits. The resilience of the VQE to quantum noise remains an open question. This study compares the computational accuracy of ground state energy estimations for molecules using the VQE across three different current quantum hardware noise models. Furthermore, our experiment on IBM's 156-qubit actual quantum computer revealed valuable insights on the real performance of the VQE on current quantum hardware.
Abstract We present a comprehensive quantum chemical study of nine molecules H₂, HeH⁺, LiH, BeH₂, H₂O, NH₃, CH₄, CO₂, and N₂ using the Variational Quantum Eigensolver (VQE) algorithm implemented in Qiskit 2.4.1. Molecular Hamiltonians are constructed in the STO-3G minimal basis set, transformed to qubit operators via the Jordan-Wigner mapping, and optimized using the n-local hardware-efficient ansatz with the SLSQP optimizer. Ground state energies are benchmarked against Full Configuration Interaction (FCI)/STO-3G reference values from authoritative literature sources. VQE achieves accuracies from 96.46% (HeH⁺) to 99.99% (H₂), average 98.68%. Potential energy surfaces are generated for all nine molecules and compared with high level ab initio literature PES references including the Kolos-Wolniewicz H₂ curve, the Polyansky water surface, and the Bytautas-Ruedenberg N₂ FCI surface. Equilibrium bond lengths agree with experiment to within 0.000–0.006 Å. Quantum ansatz circuits are analyzed in detail for all molecules. These results confirm the viability of near-term VQE algorithms for molecular quantum chemistry across diverse bond types and geometries.
Abstract Finding the ground-state energy of molecules is an important and challenging computational problem for which quantum computing can potentially find efficient solutions. The variational quantum eigensolver (VQE) is a quantum algorithm that tackles the molecular groundstate problem and is regarded as one of the flagships of quantum computing. Yet, to date, only very small molecules were computed via VQE, due to high noise levels in current quantum devices. Here we present an alternative variational quantum scheme that requires significantly less qubits than VQE. The reduction in the qubit number allows for shallower circuits to be sufficient, rendering the method more resistant to noise. The proposed algorithm, termed variational quantum selected-configuration-interaction (VQ-SCI), is based on: (a) representing the target groundstate as a superposition of Slater determinant configurations, encoded directly upon the quantum computational basis states; and (b) selecting a-priory only the most dominant configurations. This is demonstrated through a set of groundstate calculations of the H2, LiH, BeH2, H2O, NH3 and C2H4 molecules in the sto-3g basis set, performed on IBM quantum devices. We show that the VQ-SCI reaches the full configuration interaction energy within chemical accuracy using the lowest number of qubits reported to date. Moreover, when the SCI matrix is generated ‘on the fly’, the VQ-SCI requires exponentially less memory than classical SCI methods. This offers a potential remedy to a severe memory bottleneck problem in classical SCI calculations. Finally, the proposed scheme is general and can be straightforwardly applied for finding the groundstate of any Hermitian matrix, outside the chemical context.
Accurate determination of ground-state energies for molecules remains a challenge in quantum chemistry and a cornerstone for progress in fields such as drug discovery and materials design. The Variational Quantum Eigensolver (VQE) represents a leading hybrid quantum-classical paradigm for addressing this challenge; however, its widespread realization is limited by noise and the restricted scalability of current quantum hardware. Achieving efficient simulations on Noisy Intermediate-Scale Quantum (NISQ) devices and forthcoming Fault-Tolerant Application-Scalable Quantum (FASQ) systems demands a detailed understanding of how computational resources scale with molecular complexity and fermion-to-qubit encodings. In this study, resource requirements for VQE implementations employing the Unitary Coupled Cluster Singles and Doubles (UCCSD) ansatz are systematically analyzed. The molecular Hamiltonian is formulated in second quantization and mapped to qubit operators through the Jordan-Wigner (JW), Bravyi-Kitaev (BK), and Parity (Pa) transformations. Hamiltonian reduction strategies, including Z2 tapering and frozen-core approximations, are examined to assess their effect on quantum resource scaling. The analysis reveals that appropriate transformations, when combined with symmetry-based reductions, can substantially reduce qubit counts by up to ≈50% and quantum gate counts by up to ≈27.5× and Hamiltonian Pauli string counts by up to ≈2.75×, relative to the corresponding unreduced Hamiltonian representations within the same active-space configuration for the representative set of molecular systems under study. These findings provide practical circuit-level insights for executing chemically relevant simulations on NISQ hardware, while establishing physical-resource baselines that may inform future logical-level analyses targeting FASQ systems.