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:$\mathrm{H}_{2}, \mathrm{H}_{3}^{+}, \text{LiH}, \text{HF}, \text{OH}^{-}, \mathrm{H}_{2} \mathrm{O}, \text{BH}_{3}, \text{CH}_{4}$, and$\text{NH}_{3}$. 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.
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.
Accurate prediction of binding energies between drug candidates and target molecules is essential in computational drug discovery, as stronger (more negative) binding energies suggest higher efficacy. This energy corresponds to the ground state of the molecular system and can be computed using methods such as Full Configuration Interaction (FCI), Hartree-Fock (HF), Variational Quantum Eigensolver (VQE), and Quantum Phase Estimation (QPE). While FCI offers benchmark accuracy, its computational cost renders it impractical for larger systems. HF scales efficiently but at the expense of precision. QPE, though theoretically robust, is currently constrained by hardware limitations. VQE emerges as a practical alternative for Noisy Intermediate-Scale Quantum (NISQ) devices, balancing accuracy and feasibility. Nevertheless, quantum noise significantly affects VQE performance. To address this, we incorporate Zero Noise Extrapolation (ZNE) to mitigate errors and improve energy estimates. This study compares these approaches and demonstrates the efficacy of the VQEZNE framework in generating reliable binding energy calculations for drug design applications.
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.
The Variational Quantum Eigensolver (VQE) is a hybrid algorithm that combines quantum and classical computing to determine the ground-state energy of molecular systems. In this context, this study applies VQE to investigate the ground state of protocatechuic acid, analyzing its performance with various Ansatzes and active spaces. Subsequently, all VQE results were compared to those obtained with the CISD and FCI methods. The results demonstrate that Ansatzes, like Unitary Coupled Cluster Singles and Doubles (UCCSD) and variations of Hardware-Efficient Ansatzes, generally achieve accuracy close to that of FCI. In conclusion, this study highlights the effectiveness of VQE as a robust method for investigating molecular ground-state energies. Additionally, the findings emphasize the pivotal role of Ansatz design and active space selection in optimizing VQE performance, offering meaningful insights into its capabilities and constraints.
We present computational chemistry data for small molecules (CO, HCl, F2, NH4+, CH4, NH3, H3O+, H2O, BeH2, LiH, OH−, HF, HeH+, H2), obtained by implementing the Unitary Coupled Cluster method with Single and Double excitations (UCCSD) on a quantum computer simulator. We have used the Variational Quantum Eigensolver (VQE) algorithm to extract the ground state energies of these molecules. This energy data represents the expected ground state energy that a quantum computer will produce for the given molecules, on the STO-3G basis. Since there is a lot of interest in the implementation of UCCSD on quantum computers, we hope that our work will serve as a benchmark for future experimental implementations.
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.
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.