We take variational and diffusion quantum Monte Carlo (VMC and DMC) methods to investigate the chemical binding of the ground state chromium dimer, employing various single determinant (SD) or multi-determinant (MD) wavefunctions multiplied by a Jastrow factor as a trial/guiding wavefunction. The molecular orbitals (MOs) entering the SD wavefunction were calculated using restricted or unrestricted Hartree-Fock or density functional theory (DFT) calculations where five commonly-used local (SVWN5), semi-local (PW91 and BLYP), and hybrid (B1LYP and B3LYP) functionals were examined. The MD expansions were obtained from the complete-active space self-consistent field, generalized valence bond, and unrestricted configuration interaction methods. We also adopted the UB3LYP MOs to construct the MD expansion and optimized their coefficients at the VMC level. The DMC binding curves have a minimum indicating a bound state, but the comparison of atomic and molecular energies gives rise to a negative binding energy for the DMC simulations.
We take variational and diffusion quantum Monte Carlo (VMC and DMC) methods to investigate the chemical binding of the ground state chromium dimer, employing various single determinant (SD) or multi-determinant (MD) wavefunctions multiplied by a Jastrow factor as a trial/guiding wavefunction. The molecular orbitals (MOs) entering the SD wavefunction were calculated using restricted or unrestricted Hartree-Fock or density functional theory (DFT) calculations where five commonly-used local (SVWN5), semi-local (PW91 and BLYP), and hybrid (B1LYP and B3LYP) functionals were examined. The MD expansions were obtained from the complete-active space self-consistent field, generalized valence bond, and unrestricted configuration interaction methods. We also adopted the UB3LYP MOs to construct the MD expansion and optimized their coefficients at the VMC level. The DMC binding curves have a minimum indicating a bound state, but the comparison of atomic and molecular energies gives rise to a negative binding energy for the DMC simulations.
This work investigates the accuracy of real-space quantum Monte Carlo (QMC) methods for calculating molecular geometries. We present the equilibrium bond lengths of a test set of 30 diatomic molecules calculated using variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC) methods. The effect of different trial wavefunctions is investigated using single determinants constructed from Hartree-Fock (HF) and Density Functional Theory (DFT) orbitals with LDA, PBE, and B3LYP functionals, as well as small multi-configurational self-consistent field (MCSCF) multi-determinant expansions. When compared to experimental geometries, all DMC methods exhibit smaller mean-absolute deviations (MADs) than those given by HF, DFT, and MCSCF. The most accurate MAD of 3 ± 2 × 10(-3) Å is achieved using DMC with a small multi-determinant expansion. However, the more computationally efficient multi-determinant VMC method has a similar MAD of only 4.0 ± 0.9 × 10(-3) Å, suggesting that QMC forces calculated from the relatively simple VMC algorithm may often be sufficient for accurate molecular geometries.
We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.
Variational quantum Monte Carlo (VMC) method applied to interacting electron systems has played a pivotal role for the studies of electronic states of atoms, molecules and condensed matters. First, the many-body wave functions optimized by the VMC method provide good guiding functions in more accurate, diffusion quantum Monte Carlo (DMC) calculations for identical systems. Second, the VMC calculations themselves can give a good upper bound for the exact ground-state energy of many-electron systems when an appropriate form for the variational wave function is employed. In addition to a conventional form for the variational wave function that is expressed by the product of the Jastrow factor and the single Slater determinant, recent investigations have found that the expression in terms of many Slater determinants can provide a very accurate description of electronic states, especially for intrinsically multi-reference systems, which corresponds to the configuration interaction (CI) expansion method in conventional quantum chemistry. Although the CI method needs an extremely large number of expansion coefficients as variational parameters, the number of parameters can be reduced efficiently by using an idea of (molecular) orbital pair correlations. The purpose of the present paper is to examine the feasibility and usefulness of this idea within the framework of VMC method with the Jastrow–Slater wave function, where the ground state of a water molecule is investigated as a test case. In this study we rely on a variational wave function composed of a product of the Jastrow and Slater parts as 1⁄4 J S; ð1Þ
Quantum chemistry is a useful tool that provides insight into the properties and behavior of chemical systems. Modern software packages have made quantum chemistry methods more easily accessible, and the continued increase in available computational resources has allowed them to be applied to larger systems at higher levels of theory. Two significant problems that the field faces are the high computational complexity of high-level methods and the shift toward parallelism in high performance computing architectures. This work examines the treatment of weakly interacting molecular systems with the fixed-node diffusion Monte Carlo (DMC) method. DMC and other quantum Monte Carlo (QMC) methods offer a possible solution to both of the aforementioned problems: they can produce near-exact results with a lower scaling (with respect to problem size) than other similarly-accurate methods, and they are inherently parallel, so there is little additional cost associated with distributing the work of a single QMC calculation across a large number of processing units. The only error in DMC that is not systematically improvable is the constraint of a fixed nodal surface of the many-particle wave function of the system being studied. There are many cases in which a single Slater determinant trial function is sufficient to obtain accurate results, but there are others in which more sophisticated multi-determinant trial functions are necessary. Furthermore, it is non-trivial to generate nodal surfaces of similar quality for isolated and interacting molecules, so cancellation of errors is not guaranteed. We examine the use of different single- and multi-determinant trial functions in DMC calculations on small chemical systems with the goal of further understanding how to construct appropriate trial functions for molecules and clusters.
Abstract In this paper Grossman reported pseudopotential fixed-node diffusion QMC calculations of atomization energies for the Gl set of 55 small molecules,0 a standard benchmark set for testing methods in quantum chemistry. To maintain the uniformity in the calculations, they were carried out with the same basis sets and orbital types for single-determinant trial functions and effective core potentials (with a few exceptions). Variational QMC was used to optimize the determinan Jastrow trial functions. The mean absolute deviation EMAD in atomization energy from the experimental values was found to be 2.9 kcal/mol, the maximum deviation was 14 kcal/mol (for SO2), and the next five largest deviations were 7-8 kcal/mol. A comparison with results from other methods tested with the Gl set shows the FN-DQMC method to be competitive. For EMAD the list is: FN-DQMC, 2.9 kcal/mol; DFT (LDA), 40 kcal/mol; DFT (B3LYP or B3PW91), 2.9 kcal/mol; CCSD(T)/aug-cc-pVQZ, 2.8 kcal/mol; CCSD(T)/extrapolated, 1.3 kcal/mol; semi-empirical G1 theory, 1.6 kcal/mol. Additional calculations with different trial functions indicated that node location errors could be significantly reduced with multi determinant trial functions. Computational effort was not compared, but linear scaling for FN-DQMC was mentioned as an important advantage for QMC. 0J. A. Pople, M. Head-Gordon, D. J. Fox, K. Raghavachari, and L. A. Curtiss, Gaussian-I theory - A general procedure for prediction of molecular energies, J. Chem. Phys. 90, 5622 (1989).
We develop a formalism and present an algorithm for optimization of the trial wave function used in fixed-node diffusion quantum Monte Carlo (DMC) methods. The formalism is based on the DMC mixed estimator of the ground-state probability density. We take advantage of a basic property of the walker configuration distribution generated in a DMC calculation, to (i) project out a multideterminant expansion of the fixed-node ground-state wave function and (ii) to define a cost function that relates the fixed-node ground-state and the noninteracting trial wave functions. We show that (a) locally smoothing out the kink of the fixed-node ground-state wave function at the node generates a new trial wave function with better nodal structure and (b) we argue that the noise in the fixed-node wave function resulting from finite sampling plays a beneficial role, allowing the nodes to adjust toward the ones of the exact many-body ground state in a simulated annealing-like process. Based on these principles, we propose a method to improve both single determinant and multideterminant expansions of the trial wave function. The method can be generalized to other wave-function forms such as pfaffians. We test the method in a model system where benchmark configuration-interaction calculations can be performed and most components of the Hamiltonian are evaluated analytically. Comparing the DMC calculations with the exact solutions, we find that the trial wave function is systematically improved. The overlap of the optimized trial wave function and the exact ground state converges to 100% even starting from wave functions orthogonal to the exact ground state. Similarly, the DMC total energy and density converges to the exact solutions for the model. In the optimization process we find an optimal noninteracting nodal potential of density-functional-like form whose existence was predicted in a previous publication [Phys. Rev. B 77, 245110 (2008)]. Tests of the method are extended to a model system with a conventional Coulomb interaction where we show we can obtain the exact Kohn-Sham effective potential from the DMC data.