An improved method of moving the parameters in variational quantum Monte Carlo (VMC) calculations is described, and then applied to the calculation of atomic forces for first-row hydrides. We propose use of Newton's method with singular-value decomposition (SVD) to update the variational parameters, combined with steepest descent (SD) updates along directions rejected by the SVD. Using this approach to minimize the energy expectation value, we compute atomic forces [J. Chem. Phys. 118, 7193 (2003)] for first-row monohydrides in VMC and fixed-node diffusion Monte Carlo (DMC). The forces agree quantitatively with potential-energy surfaces, demonstrating the accuracy of this force procedure. Properties obtained from these calculations, such as harmonic vibrational frequencies and anharmonicity constants, agree well with experimental values.
An improved method of moving the parameters in variational quantum Monte Carlo (VMC) calculations is described, and then applied to the calculation of atomic forces for first-row hydrides. We propose use of Newton's method with singular-value decomposition (SVD) to update the variational parameters, combined with steepest descent (SD) updates along directions rejected by the SVD. Using this approach to minimize the energy expectation value, we compute atomic forces [J. Chem. Phys. 118, 7193 (2003)] for first-row monohydrides in VMC and fixed-node diffusion Monte Carlo (DMC). The forces agree quantitatively with potential-energy surfaces, demonstrating the accuracy of this force procedure. Properties obtained from these calculations, such as harmonic vibrational frequencies and anharmonicity constants, agree well with experimental values.
Ab initio quantum Monte Carlo (QMC) methods are state-of-the-art electronic structure calculations based on highly parallelizable stochastic frameworks for accurate solutions of the many-body Schrödinger equation, suitable for modern many-core supercomputer architectures. Despite its potential, one of the major drawbacks that still hinders QMC applications, especially when targeting dynamical properties of large systems or extensive datasets, is the lack of an affordable method to compute atomic forces that are consistent with the corresponding potential energy surfaces (PESs), also known as unbiased atomic forces. Recently, one of the authors in the present paper proposed a way to obtain unbiased forces with the Jastrow-correlated Slater determinant Ansatz, where the determinant part is frozen to the values obtained by a mean-field method, such as density functional theory [K. Nakano, M. Casula, and G. Tenti, Phys. Rev. B 109, 205151 (2024)]. However, the proposed method has a significant drawback for its applications: for a system with N nuclei, one requires 6N additional density functional theory (DFT) calculations to get unbiased forces, which is not negligible as the system size increases. This paper presents a way to replace the 6N DFT calculations with a single coupled-perturbed Kohn-Sham calculation, following the so-called Lagrangian technique established in quantum chemistry. This improves the computational cost and scalability of the method. We also demonstrate that the developed unbiased variational Monte Carlo (VMC) force calculation improves not only the consistency with PESs but also its accuracy, by investigating three molecules from the rMD17 benchmark set, and comparing the unbiased VMC forces with those obtained by the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] calculations. We found that the bare VMC forces are biased from the CCSD(T) ones, while the unbiased ones give values closer to those of the CCSD(T) ones. Our benchmark test also reveals that the unbiased VMC forces yield very consistent values with hybrid and meta generalized gradient approximations (e.g., ωB97X-D3BJ and ωB97M-D3BJ), but do not necessarily yield values that are very close to those of CCSD(T). Our finding paves the way to generate machine learning interatomic potentials based on VMC forces more efficiently and accurately.
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.
In this work, we introduce three algorithmic improvements to reduce the cost and improve the scaling of orbital space variational Monte Carlo (VMC). First, we show that, by appropriately screening the one- and two-electron integrals of the Hamiltonian, one can improve the efficiency of the algorithm by several orders of magnitude. This improved efficiency comes with the added benefit that the cost of obtaining a constant error per electron scales as the second power of the system size O ( N 2 ), down from the fourth power O ( N 4 ). Using numerical results, we demonstrate that the practical scaling obtained is, in fact, O ( N 1.5 ) for a chain of hydrogen atoms. Second, we show that, by using the adaptive stochastic gradient descent algorithm called AMSGrad, one can optimize the wave function energies robustly and efficiently. Remarkably, AMSGrad is almost as inexpensive as the simple stochastic gradient descent but delivers a convergence rate that is comparable to that of the Stochastic Reconfiguration algorithm, which is significantly more expensive and has a worse scaling with the system size. Third, we introduce the use of the rejection-free continuous time Monte Carlo (CTMC) to sample the determinants. Unlike the first two improvements, CTMC does come at an overhead that the local energy must be calculated at every Monte Carlo step. However, this overhead is mitigated to a large extent because of the reduced scaling algorithm, which ensures that the asymptotic cost of calculating the local energy is equal to that of updating the walker. The resulting algorithm allows us to calculate the ground state energy of a chain of 160 hydrogen atoms using a wave function containing ∼2 × 10 5 variational parameters with an accuracy of 1 mE h /particle at a cost of just 25 CPU h, which when split over 2 nodes of 24 processors each amounts to only about half hour of wall time. This low cost coupled with embarrassing parallelizability of the VMC algorithm and great freedom in the forms of usable wave functions, represents a highly effective method for calculating the electronic structure of model and ab initio systems.
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.
Quantum Monte Carlo (QMC) methods are powerful, stochastic techniques for computing the properties of interacting electrons and nuclei with an accuracy comparable to the standard post-Hartree--Fock methods of quantum chemistry. Whilst the favourable scaling of QMC methods enables a quantum, many-body treatment of much larger systems, the lack of accurate and efficient total energy derivatives, required to compute atomic forces, has hindered their widespread adoption. The work contained within this thesis provides an efficient procedure for calculating exact derivatives of QMC results. This procedure uses the programming technique of algorithmic differentiation (AD), which allows access to the derivatives of a computed function by applying chain rule differentiation to the underlying source code. However, this thesis shows that a straightforward differentiation of a stochastic function fails to capture the important contribution to the derivative from probabilistic decisions. A general approach for calculating the derivatives of a stochastic function is presented, where a similar adaptation of AD applied to the diffusion Monte Carlo (DMC) algorithm yields exact DMC atomic forces. The approach is validated by performing the largest ever DMC force calculations, which demonstrate the feasibility of treating systems containing thousands of electrons. The efficiency of AD also enables molecular dynamics simulations driven entirely by DMC, adding new functionality to the QMC toolkit. Another focus of this thesis is using the phenomenon of stochastic coherence to correlate DMC simulations, allowing finite difference derivatives to be obtained with a small error. Whilst this method is far easier to implement than AD, preliminary results show an instability when treating larger systems. A different approach is obtained from extrapolating this method to a finite difference step size of zero, producing algebraic expressions for a direct differentiation of the DMC algorithm.