Quantum computers are reaching one crucial milestone after another. Motivated by their progress in quantum chemistry, we performed an extensive series of simulations of quantum-computer runs that were aimed at inspecting the best-practice aspects of these calculations. In order to compare the performance of different setups, the ground-state energy of the hydrogen molecule was chosen as a benchmark for which the exact solution exists in the literature. Applying the variational quantum eigensolver (VQE) to a qubit Hamiltonian obtained by the Bravyi-Kitaev transformation, we analyzed the impact of various computational technicalities. These included (i) the choice of the optimization methods, (ii) the architecture of the quantum circuits, as well as (iii) the different types of noise when simulating real quantum processors. On these, we eventually performed a series of experimental runs as a complement to our simulations. The simultaneous perturbation stochastic approximation (SPSA) and constrained optimization by linear approximation (COBYLA) optimization methods clearly outperformed the Nelder-Mead and Powell methods. The results obtained when using the Ry variational form were better than those obtained when the RyRz form was used. The choice of an optimum entangling layer was sensitively interlinked with the choice of the optimization method. The circular entangling layer was found to worsen the performance of the COBYLA method, while the full-entangling layer improved it. All four optimization methods sometimes led to an energy that corresponded to an excited state rather than the ground state. We also show that a similarity analysis of measured probabilities can provide a useful insight.
Quantum computers are reaching one crucial milestone after another. Motivated by their progress in quantum chemistry, we performed an extensive series of simulations of quantum-computer runs that were aimed at inspecting the best-practice aspects of these calculations. In order to compare the performance of different setups, the ground-state energy of the hydrogen molecule was chosen as a benchmark for which the exact solution exists in the literature. Applying the variational quantum eigensolver (VQE) to a qubit Hamiltonian obtained by the Bravyi-Kitaev transformation, we analyzed the impact of various computational technicalities. These included (i) the choice of the optimization methods, (ii) the architecture of the quantum circuits, as well as (iii) the different types of noise when simulating real quantum processors. On these, we eventually performed a series of experimental runs as a complement to our simulations. The simultaneous perturbation stochastic approximation (SPSA) and constrained optimization by linear approximation (COBYLA) optimization methods clearly outperformed the Nelder-Mead and Powell methods. The results obtained when using the Ry variational form were better than those obtained when the RyRz form was used. The choice of an optimum entangling layer was sensitively interlinked with the choice of the optimization method. The circular entangling layer was found to worsen the performance of the COBYLA method, while the full-entangling layer improved it. All four optimization methods sometimes led to an energy that corresponded to an excited state rather than the ground state. We also show that a similarity analysis of measured probabilities can provide a useful insight.
Classical optimizers play a crucial role in determining the accuracy and convergence of variational quantum algorithms; leading algorithms use a near-term quantum computer to solve the ground state properties of molecules, simulate dynamics of different quantum systems, and so on. In the literature, many optimizers, each having its own architecture, have been employed expediently for different applications. In this work, we consider a few popular and efficacious optimizers and assess their performance in variational quantum algorithms for applications in quantum chemistry in a realistic noisy setting. We benchmark the optimizers with critical analysis based on quantum simulations of simple molecules, such as hydrogen, lithium hydride, beryllium hydride, water, and hydrogen fluoride. The errors in the ground state energy, dissociation energy, and dipole moment are the parameters used as yardsticks. All the simulations were carried out with an ideal quantum circuit simulator, a noisy quantum circuit simulator, and finally a noisy simulator with noise embedded from the IBM Cairo quantum device to understand the performance of the classical optimizers in ideal and realistic quantum environments. We used the standard unitary coupled cluster ansatz for simulations, and the number of qubits varied from two starting from the hydrogen molecule to ten qubits in hydrogen fluoride. Based on the performance of these optimizers in the ideal quantum circuits, the conjugate gradient, limited-memory Broyden-Fletcher-Goldfarb-Shanno bound, and sequential least squares programming optimizers are found to be the best-performing gradient-based optimizers. While constrained optimization by linear approximation (COBYLA) and Powell's conjugate direction algorithm for unconstrained optimization (POWELL) perform most efficiently among the gradient-free methods, in noisy quantum circuit conditions, simultaneous perturbation stochastic approximation, POWELL, and COBYLA are among the best-performing optimizers.
New approaches into computational quantum chemistry can be developed through the use of quantum computing. While universal, fault-tolerant quantum computers are still not available, and we want to utilize today's noisy quantum processors. One of their flagship applications is the variational quantum eigensolver (VQE)-an algorithm for calculating the minimum energy of a physical Hamiltonian. In this study, we investigate how various types of errors affect the VQE and how to efficiently use the available resources to produce precise computational results. We utilize a simulator of a noisy quantum device, an exact statevector simulator, and physical quantum hardware to study the VQE algorithm for molecular hydrogen. We find that the optimal method of running the hybrid classical-quantum optimization is to: (i) allow some noise in intermediate energy evaluations, using fewer shots per step and fewer optimization iterations, but ensure a high final readout precision; (ii) emphasize efficient problem encoding and ansatz parametrization; and (iii) run all experiments within a short time-frame, avoiding parameter drift with time. Nevertheless, current publicly available quantum resources are still very noisy and scarce/expensive, and even when using them efficiently, it is quite difficult to perform trustworthy calculations of molecular energies.
The Variational Quantum Eigensolver (VQE) is a promising quantum algorithm for applications in chemistry within the Noisy Intermediate-Scale Quantum (NISQ) era. The ability for a quantum computer to simulate electronic structures with high accuracy would have a profound impact on material and biochemical science with potential applications e.g., to the development of new drugs. However, considering the variety of quantum hardware architectures, it is still uncertain which hardware concept is most suited to execute the VQE for e.g., the simulation of molecules. Aspects to consider here are the required connectivity of the quantum circuit used, the size and the depth and thus the susceptibility to noise effects. Besides theo-retical considerations, empirical studies using available quantum hardware may help to clarify the question of which hardware technology might be better suited for a certain given application and algorithm. Going one step into this direction, within this work, we present results using the VQE for the simulation of the hydrogen molecule, comparing superconducting and ion trap quantum computers. The experiments are carried out with a standardized setup of ansatz and optimizer, selected to reduce the number of required iterations. The findings are analyzed considering different quantum processor types, calibration data as well as the depth and gate counts of the circuits required for the different hardware concepts after transpilation.
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 work investigates the performance of numerical optimization algorithms applied to the State-Averaged Orbital-Optimized Variational Quantum Eigensolver for the H2 molecule under various quantum noise conditions. The goal is to assess the stability, accuracy, and computational efficiency of commonly used gradient-based, gradient-free, and global optimization strategies within the Noisy Intermediate-Scale Quantum regime. We systematically compare six representative optimizers, BFGS, SLSQP, Nelder-Mead, Powell, COBYLA, and iSOMA, under ideal, stochastic, and decoherence noise models, including phase damping, depolarizing, and thermal relaxation channels. Each optimizer was tested over multiple noise intensities and measurement settings to characterize convergence behavior and sensitivity to noise-induced landscape distortions. The results show that BFGS consistently achieves the most accurate energies with minimal evaluations, maintaining robustness even under moderate decoherence. COBYLA performs well for low-cost approximations, while SLSQP exhibits instability in noisy regimes. Global approaches such as iSOMA show potential but are computationally expensive. These findings provide practical guidance for selecting suitable optimizers in variational quantum simulations, highlighting the importance of noise-aware optimization strategies for reliable and efficient quantum chemistry computations on current hardware.
The variational quantum eigensolver (or VQE), first developed by Peruzzo et al. (2014), has received significant attention from the research community in recent years. It uses the variational principle to compute the ground state energy of a Hamiltonian, a problem that is central to quantum chemistry and condensed matter physics. Conventional computing methods are constrained in their accuracy due to the computational limits facing exact modeling of the exponentially growing electronic wavefunction for these many-electron systems. The VQE may be used to model these complex wavefunctions in polynomial time, making it one of the most promising near-term applications for quantum computing. One important advantage is that variational algorithms have been shown to present some degree of resilience to the noise in the quantum hardware. Finding a path to navigate the relevant literature has rapidly become an overwhelming task, with many methods promising to improve different parts of the algorithm, but without clear descriptions of how the diverse parts fit together. The potential practical advantages of the algorithm are also widely discussed in the literature, but with varying conclusions. Despite strong theoretical underpinnings suggesting excellent scaling of individual VQE components, studies have pointed out that their various pre-factors could be too large to reach a quantum computing advantage over conventional methods. This review aims at disentangling the relevant literature to provide a comprehensive overview of the progress that has been made on the different parts of the algorithm, and to discuss future areas of research that are fundamental for the VQE to deliver on its promises. All the different components of the algorithm are reviewed in detail. These include the representation of Hamiltonians and wavefunctions on a quantum computer, the optimization process to find ground state energies, the post processing mitigation of quantum errors, and suggested best practices. We identify four main areas of future research: (1) optimal measurement schemes for reduction of circuit repetitions required; (2) large scale parallelization across many quantum computers; (3) ways to overcome the potential appearance of vanishing gradients in the optimization process for large systems, and how the number of iterations required for the optimization scales with system size; (4) the extent to which VQE suffers for quantum noise, and whether this noise can be mitigated in a tractable manner. The answers to these open research questions will determine the routes for the VQE to achieve quantum advantage as the quantum computing hardware scales up and as the noise levels are reduced.