Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Progress in Natural Science Materials International

Variance minimization for variational quantum Monte Carlo method

HX Huang, ZX Cao, Shubin Liu, 曹泽星

Published Oct 1, 1997• 0 citations
OpenAlex

Abstract

A new algorithm for the variational quantum Monte Carlo (VMC) is proposed. This algorithm takesthe internal structure of a local energy as starting point, and directly reduces its fluctuation, making the variance decrease to the minimum. An analytical expression of the local energy is presented. The principle of variance minimization for the VMC is deduced, and the steps of variance minimization are established. This new algorithm was used to compute the energies of the ground states for molecules H2, LiH, Li2 and H2O. The Calculated results show that this novel algorithm is superior to the usual VMC in both accuracy and statistical error, with but very little increase in computational cost.

Related papers

Variational Principles in Quantum Monte Carlo: The Troubled Story of Variance Minimization

Alice Cuzzocrea, Anthony Scemama, Wim J. Briels +2 more

We investigate the use of different variational principles in quantum Monte Carlo, namely, energy and variance minimization, prompted by the interest in the robust and accurate estimation of electronic excited states. For two prototypical, challenging molecules, we readily reach the accuracy of the best available reference excitation energies using energy minimization in a state-specific or state-average fashion for states of different or equal symmetry, respectively. On the other hand, in variance minimization, where the use of suitable functionals is expected to target specific states regardless of the symmetry, we encounter severe problems for a variety of wave functions: as the variance converges, the energy drifts away from that of the selected state. This unexpected behavior is sometimes observed even when the target is the ground state and generally prevents the robust estimation of total and excitation energies. We analyze this problem using a very simple wave function and infer that the optimization finds little or no barrier to escape from a local minimum or local plateau, eventually converging to a lower-variance state instead of the target state. For the increasingly complex systems becoming in reach of quantum Monte Carlo simulations, variance minimization with current functionals appears to be an impractical route.

arXiv:2002.05951•Journal of Chemical Theory and Computation•May 17, 2020•38 citations

A novel method for optimizing quantum Monte Carlo wave functions

Hongxin Huang, Zexing Cao

We propose an algorithm for optimizing quantum Monte Carlo wave functions. An improved steepest-descent technique is used, with step size automatically adjustable to obtain a procedure that converges superlinearly. We also propose a novel trial function, which has both correct electron–electron cusp conditions and correct electron–nucleus cusp conditions. To test the optimization procedure and the optimized trial function, the ground states for CH4 and H2O molecules were investigated using variational Monte Carlo (VMC) and fixed-node quantum Monte Carlo (FNQMC) calculations. For CH4 and H2O, the VMC recovered 73.3% and 57.9% of the correlation energy, respectively, and the FNQMC recovered 99.3% and 92.8%, respectively. The optimization procedure is three to five times faster than conventional usual steepest-descent procedures. The trial functions optimized are more accurate than prior trial functions of similar complexity.

The Journal of Chemical Physics•Jan 1, 1996•26 citations

Variance minimization variational Monte Carlo method

Imran Khan, Bo Gao

We present a variational Monte Carlo (VMC) method that works equally well for the ground and the excited states of a quantum system. The method is based on the minimization of the variance of energy, as opposed to the energy itself in standard methods. As a test, it is applied to the investigation of the universal spectrum at the van der Waals length scale for two identical Bose atoms in a symmetric harmonic trap, with results compared to the basically exact results obtained from a multiscale quantum-defect theory.

arXiv:physics/0701223•APS•Jan 19, 2007•0 citations

Accelerating variational quantum Monte Carlo using the variational quantum eigensolver

Ashley Montanaro, Stasja Stanisic

Variational Monte Carlo (VMC) methods are used to sample classically from distributions corresponding to quantum states which have an efficient classical description. VMC methods are based on performing a number of steps of a Markov chain starting with samples from a simple initial distribution. Here we propose replacing this initial distribution with samples produced using a quantum computer, for example using the variational quantum eigensolver (VQE). We show that, based on the use of initial distributions generated by numerical simulations and by experiments on quantum hardware, convergence to the target distribution can be accelerated compared with classical samples; the energy can be reduced compared with the energy of the state produced by VQE; and VQE states produced by small quantum computers can be used to accelerate large instances of VMC. Quantum-enhanced VMC makes minimal requirements of the quantum computer and offers the prospect of accelerating classical methods using noisy samples from near-term quantum computers which are not yet able to accurately represent ground states of complex quantum systems.

arXiv:2307.07719•Jul 15, 2023•3 citations

Optimization of quantum Monte Carlo wave functions by energy minimization

Julien Toulouse, Cyrus J. Umrigar

We study three wave function optimization methods based on energy minimization in a variational Monte Carlo framework: the Newton, linear, and perturbative methods. In the Newton method, the parameter variations are calculated from the energy gradient and Hessian, using a reduced variance statistical estimator for the latter. In the linear method, the parameter variations are found by diagonalizing a nonsymmetric estimator of the Hamiltonian matrix in the space spanned by the wave function and its derivatives with respect to the parameters, making use of a strong zero-variance principle. In the less computationally expensive perturbative method, the parameter variations are calculated by approximately solving the generalized eigenvalue equation of the linear method by a nonorthogonal perturbation theory. These general methods are illustrated here by the optimization of wave functions consisting of a Jastrow factor multiplied by an expansion in configuration state functions (CSFs) for the C2 molecule, including both valence and core electrons in the calculation. The Newton and linear methods are very efficient for the optimization of the Jastrow, CSF, and orbital parameters. The perturbative method is a good alternative for the optimization of just the CSF and orbital parameters. Although the optimization is performed at the variational Monte Carlo level, we observe for the C2 molecule studied here, and for other systems we have studied, that as more parameters in the trial wave functions are optimized, the diffusion Monte Carlo total energy improves monotonically, implying that the nodal hypersurface also improves monotonically.

arXiv:physics/0701039•The Journal of Chemical Physics•Feb 22, 2007•290 citations

Improved energies and wave function accuracy with Weighted Variational Monte Carlo

Zhang, Huan, Robert J. Webber, Michael Lindsey +2 more

Neural network parametrizations have increasingly been used to represent the ground and excited states in variational Monte Carlo (VMC) with promising results. However, traditional VMC methods only optimize the wave function in regions of peak probability. The wave function is uncontrolled in the tails of the probability distribution, which can limit the accuracy of the trained wavefunction approximation. To improve the approximation accuracy in the probability tails, this paper interprets VMC as a gradient flow in the space of wave functions, followed by a projection step. From this perspective, arbitrary probability distributions can be used in the projection step, allowing the user to prioritize accuracy in different regions of state space. Motivated by this theoretical perspective, the paper tests a new weighted VMC method on the antiferromagnetic Heisenberg model for a periodic spin chain. Compared to traditional VMC, weighted VMC reduces the error in the ground state energy by a factor of 2 and it reduces the errors in the local energies away from the mode by large factors of $10^2$--$10^4$.

arXiv:2507.01905•Jun 17, 2025•0 citations
A new algorithm for the variational quantum Monte Carlo (VMC) is proposed. This algorithm takesthe internal structure of a local energy as starting point, and directly reduces its fluctuation, making the variance decrease to the minimum. An analytical expression of the local energy is presented. The principle of variance minimization for the VMC is deduced, and the steps of variance minimization are established. This new algorithm was used to compute the energies of the ground states for molecules H2, LiH, Li2 and H2O. The Calculated results show that this novel algorithm is superior to the usual VMC in both accuracy and statistical error, with but very little increase in computational cost.
We investigate the use of different variational principles in quantum Monte Carlo, namely, energy and variance minimization, prompted by the interest in the robust and accurate estimation of electronic excited states. For two prototypical, challenging molecules, we readily reach the accuracy of the best available reference excitation energies using energy minimization in a state-specific or state-average fashion for states of different or equal symmetry, respectively. On the other hand, in variance minimization, where the use of suitable functionals is expected to target specific states regardless of the symmetry, we encounter severe problems for a variety of wave functions: as the variance converges, the energy drifts away from that of the selected state. This unexpected behavior is sometimes observed even when the target is the ground state and generally prevents the robust estimation of total and excitation energies. We analyze this problem using a very simple wave function and infer that the optimization finds little or no barrier to escape from a local minimum or local plateau, eventually converging to a lower-variance state instead of the target state. For the increasingly complex systems becoming in reach of quantum Monte Carlo simulations, variance minimization with current functionals appears to be an impractical route.
We propose an algorithm for optimizing quantum Monte Carlo wave functions. An improved steepest-descent technique is used, with step size automatically adjustable to obtain a procedure that converges superlinearly. We also propose a novel trial function, which has both correct electron–electron cusp conditions and correct electron–nucleus cusp conditions. To test the optimization procedure and the optimized trial function, the ground states for CH4 and H2O molecules were investigated using variational Monte Carlo (VMC) and fixed-node quantum Monte Carlo (FNQMC) calculations. For CH4 and H2O, the VMC recovered 73.3% and 57.9% of the correlation energy, respectively, and the FNQMC recovered 99.3% and 92.8%, respectively. The optimization procedure is three to five times faster than conventional usual steepest-descent procedures. The trial functions optimized are more accurate than prior trial functions of similar complexity.
We present a variational Monte Carlo (VMC) method that works equally well for the ground and the excited states of a quantum system. The method is based on the minimization of the variance of energy, as opposed to the energy itself in standard methods. As a test, it is applied to the investigation of the universal spectrum at the van der Waals length scale for two identical Bose atoms in a symmetric harmonic trap, with results compared to the basically exact results obtained from a multiscale quantum-defect theory.
Variational Monte Carlo (VMC) methods are used to sample classically from distributions corresponding to quantum states which have an efficient classical description. VMC methods are based on performing a number of steps of a Markov chain starting with samples from a simple initial distribution. Here we propose replacing this initial distribution with samples produced using a quantum computer, for example using the variational quantum eigensolver (VQE). We show that, based on the use of initial distributions generated by numerical simulations and by experiments on quantum hardware, convergence to the target distribution can be accelerated compared with classical samples; the energy can be reduced compared with the energy of the state produced by VQE; and VQE states produced by small quantum computers can be used to accelerate large instances of VMC. Quantum-enhanced VMC makes minimal requirements of the quantum computer and offers the prospect of accelerating classical methods using noisy samples from near-term quantum computers which are not yet able to accurately represent ground states of complex quantum systems.
We study three wave function optimization methods based on energy minimization in a variational Monte Carlo framework: the Newton, linear, and perturbative methods. In the Newton method, the parameter variations are calculated from the energy gradient and Hessian, using a reduced variance statistical estimator for the latter. In the linear method, the parameter variations are found by diagonalizing a nonsymmetric estimator of the Hamiltonian matrix in the space spanned by the wave function and its derivatives with respect to the parameters, making use of a strong zero-variance principle. In the less computationally expensive perturbative method, the parameter variations are calculated by approximately solving the generalized eigenvalue equation of the linear method by a nonorthogonal perturbation theory. These general methods are illustrated here by the optimization of wave functions consisting of a Jastrow factor multiplied by an expansion in configuration state functions (CSFs) for the C2 molecule, including both valence and core electrons in the calculation. The Newton and linear methods are very efficient for the optimization of the Jastrow, CSF, and orbital parameters. The perturbative method is a good alternative for the optimization of just the CSF and orbital parameters. Although the optimization is performed at the variational Monte Carlo level, we observe for the C2 molecule studied here, and for other systems we have studied, that as more parameters in the trial wave functions are optimized, the diffusion Monte Carlo total energy improves monotonically, implying that the nodal hypersurface also improves monotonically.
Neural network parametrizations have increasingly been used to represent the ground and excited states in variational Monte Carlo (VMC) with promising results. However, traditional VMC methods only optimize the wave function in regions of peak probability. The wave function is uncontrolled in the tails of the probability distribution, which can limit the accuracy of the trained wavefunction approximation. To improve the approximation accuracy in the probability tails, this paper interprets VMC as a gradient flow in the space of wave functions, followed by a projection step. From this perspective, arbitrary probability distributions can be used in the projection step, allowing the user to prioritize accuracy in different regions of state space. Motivated by this theoretical perspective, the paper tests a new weighted VMC method on the antiferromagnetic Heisenberg model for a periodic spin chain. Compared to traditional VMC, weighted VMC reduces the error in the ground state energy by a factor of 2 and it reduces the errors in the local energies away from the mode by large factors of 10210^2102--10410^4104.