Optimization of configuration interaction coefficients in multideterminant Jastrow–Slater wave functions

Optimization of quantum Monte Carlo wave functions by energy minimization

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.

Quantum Monte Carlo study of the Ne atom and the Ne+ ion

We report all-electron and pseudopotential calculations of the ground-state energies of the neutral Ne atom and the Ne(+) ion using the variational and diffusion quantum Monte Carlo (DMC) methods. We investigate different levels of Slater-Jastrow trial wave function: (i) using Hartree-Fock orbitals, (ii) using orbitals optimized within a Monte Carlo procedure in the presence of a Jastrow factor, and (iii) including backflow correlations in the wave function. Small reductions in the total energy are obtained by optimizing the orbitals, while more significant reductions are obtained by incorporating backflow correlations. We study the finite-time-step and fixed-node biases in the DMC energy and show that there is a strong tendency for these errors to cancel when the first ionization potential (IP) is calculated. DMC gives highly accurate values for the IP of Ne at all the levels of trial wave function that we have considered.

Energy and variance optimization of many-body wave functions

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.

Delayed rejection variational Monte Carlo

An acceleration algorithm to address the problem of multiple time scales in variational Monte Carlo simulations is presented. After a first attempted move has been rejected, the delayed rejection algorithm attempts a second move with a smaller time step, so that even moves of the core electrons can be accepted. Results on Be and Ne atoms as test cases are presented. Correlation time and both average accepted displacement and acceptance ratio as a function of the distance from the nucleus evidence the efficiency of the proposed algorithm in dealing with the multiple time scales problem.

Optimized Jastrow–Slater wave functions for ground and excited states: Application to the lowest states of ethene

A quantum Monte Carlo method is presented for determining multideterminantal Jastrow-Slater wave functions for which the energy is stationary with respect to the simultaneous optimization of orbitals and configuration interaction coefficients. The approach is within the framework of the so-called energy fluctuation potential method which minimizes the energy in an iterative fashion based on Monte Carlo sampling and a fitting of the local energy fluctuations. The optimization of the orbitals is combined with the optimization of the configuration interaction coefficients through the use of additional single excitations to a set of external orbitals. A new set of orbitals is then obtained from the natural orbitals of this enlarged configuration interaction expansion. For excited states, the approach is extended to treat the average of several states within the same irreducible representation of the pointgroup of the molecule. The relationship of our optimization method with the stochastic reconfiguration technique by Sorella et al. is examined. Finally, the performance of our approach is illustrated with the lowest states of ethene, in particular with the difficult case of the 1(1)B(1u) state.