Abstract Variational quantum eigensolvers are considered promising approaches for quantum chemistry on near-term quantum hardware, yet their performance is strongly affected by sampling noise. Despite the widespread use of classical optimizers, their robustness under realistic finite-shot sampling conditions remains insufficiently understood.The optimization of the truncated Variational Hamiltonian Ansatz for H2, H4, and LiH is investigated using eight classical optimization algorithms. A total of 1,280 independent optimization runs are performed under both exact statevector and finite-shot sampling conditions.Finite-shot sampling is found to reshape optimization dynamics and reverse optimizer rankings observed in noiseless simulations. While gradient-based methods achieve near-machine precision in exact simulations, their performance deteriorates substantially under sampling noise. In contrast, population-based approaches, particularly the Covariance Matrix Adaptation Evolution Strategy, remain robust, achieving corrected errors below 5 x 10-4 Ha for H2 and below 10-2 Ha for LiH.Finite-shot sampling is further shown to frequently produce apparent violations of the variational principle. By reevaluating candidate solutions using 105 measurements, sampling-induced fluctuations are exploited to obtain energy estimates beyond the intrinsic sampling limit.These results provide practical guidelines for optimizer selection and highlight the importance of high-shot reevaluation in noisy variational quantum chemistry.
Abstract Variational quantum eigensolvers are considered promising approaches for quantum chemistry on near-term quantum hardware, yet their performance is strongly affected by sampling noise. Despite the widespread use of classical optimizers, their robustness under realistic finite-shot sampling conditions remains insufficiently understood.The optimization of the truncated Variational Hamiltonian Ansatz for H2, H4, and LiH is investigated using eight classical optimization algorithms. A total of 1,280 independent optimization runs are performed under both exact statevector and finite-shot sampling conditions.Finite-shot sampling is found to reshape optimization dynamics and reverse optimizer rankings observed in noiseless simulations. While gradient-based methods achieve near-machine precision in exact simulations, their performance deteriorates substantially under sampling noise. In contrast, population-based approaches, particularly the Covariance Matrix Adaptation Evolution Strategy, remain robust, achieving corrected errors below 5 x 10-4 Ha for H2 and below 10-2 Ha for LiH.Finite-shot sampling is further shown to frequently produce apparent violations of the variational principle. By reevaluating candidate solutions using 105 measurements, sampling-induced fluctuations are exploited to obtain energy estimates beyond the intrinsic sampling limit.These results provide practical guidelines for optimizer selection and highlight the importance of high-shot reevaluation in noisy variational quantum chemistry.
The optimization of Variational Quantum Eigensolver is severely challenged by finite-shot sampling noise, which distorts the cost landscape, creates false variational minima, and induces statistical bias called winner's curse. We investigate this phenomenon by benchmarking eight classical optimizers spanning gradient-based, gradient-free, and metaheuristic methods on quantum chemistry Hamiltonians H2, H4 chain, LiH (in both full and active spaces) using the truncated Variational Hamiltonian Ansatz. We analyze difficulties of gradient-based methods (e.g., SLSQP, BFGS) in noisy regimes, where they diverge or stagnate. We show that the bias of estimator can be corrected by tracking the \textit{population mean}, rather than the biased best individual when using population based optimizer. Our findings, which are shown to generalize to hardware-efficient circuits and condensed matter models, identify adaptive metaheuristics (specifically CMA-ES and iL-SHADE) as the most effective and resilient strategies. We conclude by presenting a set of practical guidelines for reliable VQE optimization under noise, centering on the co-design of physically motivated ansatz and the use of adaptive 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.
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.
Quantum chemistry is a near-term application for quantum computers. This application may be facilitated by variational quantum-classical algorithms (VQCAs), although a concern for VQCAs is the large number of measurements needed for convergence, especially for chemical accuracy. Here we introduce a strategy for reducing the number of measurements (i.e., shots) by randomly sampling operators hi from the overall Hamiltonian H=∑icihi. In particular, we employ weighted sampling, which is important when the ci's are highly non-uniform, as is typical in chemistry. We integrate this strategy with an adaptive optimizer developed recently by our group to construct an improved optimizer called Rosalin (Random Operator Sampling for Adaptive Learning with Individual Number of shots). Rosalin implements stochastic gradient descent while adapting the shot noise for each partial derivative and randomly assigning the shots amongst the hi according to a weighted distribution. We implement this and other optimizers to find the ground states of molecules H2, LiH, and BeH2, without and with quantum hardware noise, and Rosalin outperforms other optimizers in most cases.
Variational quantum algorithms (VQAs) offer a promising path toward using near-term quantum hardware for applications in academic and industrial research. These algorithms aim to find approximate solutions to quantum problems by optimizing a parametrized quantum circuit using a classical optimization algorithm. A successful VQA requires fast and reliable classical optimization algorithms. Understanding and optimizing how off-the-shelf optimization methods perform in this context is important for the future of the field. In this work, we study the performance of four commonly used gradient-free optimization methods [sequential least-squares quadratic programming, constrained optimization by linear approximations, the covariance matrix adaptation evolutionary strategy (CMA-ES), and the simultaneous perturbation stochastic approximation (SPSA)] to find ground-state energies of a range of small chemistry and material science problems. We test a telescoping sampling scheme (where the accuracy of the cost-function estimate provided to the optimizer is increased as the optimization converges) for all methods, demonstrating mixed results across our range of optimizers and problems chosen. We further hyperparameter tune two of the four optimizers (CMA-ES and SPSA) across a large range of models and demonstrate that with appropriate hyperparameter tuning, CMA-ES is competitive with and sometimes outperforms SPSA (which is not observed in the absence of hyperparameter tuning). Finally, we investigate the ability of an optimizer to beat the ``sampling-noise floor'' given by the sampling noise of each cost-function estimate provided to the optimizer. Our results demonstrate the necessity for tailoring and hyperparameter tuning known optimization techniques for inherently noisy variational quantum algorithms and that the variational landscape that one finds in a VQA is highly problem and system dependent. This provides guidance for future implementations of these algorithms in experiments.
As a hybrid quantum-classical algorithm, the variational quantum eigensolver is widely applied in quantum chemistry simulations, especially in computing the electronic structure of complex molecular systems. However, on existing noisy intermediate-scale quantum devices, some factors such as quantum decoherence, measurement errors, and gate operation imprecisions are unavoidable. To overcome these challenges, this study proposes an efficient noise-mitigating variational quantum eigensolver for accurate computation of molecular ground state energies in noisy environments. We design the quantum circuit with reference to the structure of matrix product states and utilize it to pre-train the circuit parameters, which ensures circuit stability and mitigates fluctuations caused by initialization. We also employ zero-noise extrapolation to mitigate quantum noise and combine it with neural networks to improve the accuracy of the noise-fitting function, which significantly eliminates noise interference. Furthermore, we implement an intelligent grouping strategy for measuring Hamiltonian Pauli strings, which not only reduces measurement errors but also improves sampling efficiency. We perform numerical simulations to solve the ground state energy of the H4molecule by using MindSpore Quantum framework, and the results demonstrate that our algorithm can constrain noise errors within the range of O(10−2)∼O(10−1), outperforming mainstream variational quantum eigensolvers. This work provides a new strategy for high-precision quantum chemistry calculations on near-term noisy quantum hardware.