We present a simple and efficient method to optimize within energy minimization the determinantal component of the many-body wave functions commonly used in quantum Monte Carlo calculations. The approach obtains the optimal wave function as an approximate perturbative solution of an effective Hamiltonian iteratively constructed via Monte Carlo sampling. The effectiveness of the method, as well as its ability to substantially improve the accuracy of quantum Monte Carlo calculations, are demonstrated by optimizing a large number of parameters for the ground state of acetone and the difficult case of the $1^{1}B_{u}$ state of hexatriene.
We present a simple and efficient method to optimize within energy minimization the determinantal component of the many-body wave functions commonly used in quantum Monte Carlo calculations. The approach obtains the optimal wave function as an approximate perturbative solution of an effective Hamiltonian iteratively constructed via Monte Carlo sampling. The effectiveness of the method, as well as its ability to substantially improve the accuracy of quantum Monte Carlo calculations, are demonstrated by optimizing a large number of parameters for the ground state of acetone and the difficult case of the 11Bu state of hexatriene.
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.
We present a simple, robust, and efficient method for varying the parameters in a many-body wave function to optimize the expectation value of the energy. The effectiveness of the method is demonstrated by optimizing the parameters in flexible Jastrow factors that include 3-body electron-electron-nucleus correlation terms for the NO2 and decapentaene (C10H12) molecules. The basic idea is to add terms to the straightforward expression for the Hessian of the energy that have zero expectation value, but that cancel much of the statistical fluctuations for a finite Monte Carlo sample. The method is compared to what is currently the most popular method for optimizing many-body wave functions, namely, minimization of the variance of the local energy. The most efficient wave function is obtained by optimizing a linear combination of the energy and the variance.
We present a simple, robust, and highly efficient method for optimizing all parameters of many-body wave functions in quantum Monte Carlo calculations, applicable to continuum systems and lattice models. Based on a strong zero-variance principle, diagonalization of the Hamiltonian matrix in the space spanned by the wave function and its derivatives determines the optimal parameters. It systematically reduces the fixed-node error, as demonstrated by the calculation of the binding energy of the small but challenging C(2) molecule to the experimental accuracy of 0.02 eV.
An appropriate iterative scheme for the minimization of the energy, based on the variational Monte Carlo (VMC) technique, is introduced and compared with existing stochastic schemes. We test the various methods for the one-dimensional Heisenberg ring and the two-dimensional $t\text{\ensuremath{-}}J$ model and show that, with the present scheme, very accurate and efficient calculations are possible, even for several variational parameters. Indeed, by using a very efficient statistical evaluation of the first and the second energy derivatives, it is possible to define a very rapidly converging iterative scheme that, within VMC, is much more convenient than the standard Newton method. It is also shown how to optimize simultaneously both the Jastrow and the determinantal part of the wave function.
We investigate Monte Carlo energy and variance-minimization techniques for optimizing many-body wave functions. Several variants of the basic techniques are studied, including limiting the variations in the weighting factors that arise in correlated sampling estimations of the energy and its variance. We investigate the numerical stability of the techniques and identify two reasons why variance minimization exhibits superior numerical stability to energy minimization. The characteristics of each method are studied using a noninteracting 64-electron model of crystalline silicon. While our main interest is in solid-state systems, the issues investigated are relevant to Monte Carlo studies of atoms, molecules, and solids. We identify a robust and efficient variance-minimization scheme for optimizing wave functions for large systems.
We present a new method for the optimization of large configuration interaction (CI) expansions in the quantum Monte Carlo (QMC) framework. The central idea here is to replace the nonorthogonal variational optimization of CI coefficients performed in usual QMC calculations by an orthogonal non-Hermitian optimization thanks to the so-called transcorrelated (TC) framework, the two methods yielding the same results in the limit of a complete basis set. By rewriting the TC equations as an effective self-consistent Hermitian problem, our approach requires the sampling of a single quantity per Slater determinant, leading to minimal memory requirements in the QMC code. Using analytical quantities obtained from both the TC framework and the usual CI-type calculations, we also propose improved estimators which reduce the statistical fluctuations of the sampled quantities by more than an order of magnitude. We demonstrate the efficiency of this method on wave functions containing 10 5 –10 6 Slater determinants, using effective core potentials or all-electron calculations. In all the cases, a sub-milli-Hartree convergence is reached within only two or three iterations of optimization.