We investigate the performance of a class of compact and systematically improvable Jastrow-Slater wave functions for the efficient and accurate computation of structural properties, where the determinantal component is expanded with a perturbatively selected configuration interaction scheme (CIPSI). We concurrently optimize the molecular ground-state geometry and full wave function-Jastrow factor, orbitals, and configuration interaction coefficients-in variational Monte Carlo (VMC) for the prototypical case of 1,3- trans-butadiene, a small yet theoretically challenging π-conjugated system. We find that the CIPSI selection outperforms the conventional scheme of correlating orbitals within active spaces chosen by chemical intuition: it gives significantly better variational and diffusion Monte Carlo energies for all but the smallest expansions, and much smoother convergence of the geometry with the number of determinants. In particular, the optimal bond lengths and bond-length alternation of butadiene are converged to better than 1 mÅ with just a few thousand determinants, to values very close to the corresponding CCSD(T) results. The combination of CIPSI expansion and VMC optimization represents an affordable tool for the determination of accurate ground-state geometries in quantum Monte Carlo.
We investigate the performance of a class of compact and systematically improvable Jastrow-Slater wave functions for the efficient and accurate computation of structural properties, where the determinantal component is expanded with a perturbatively selected configuration interaction scheme (CIPSI). We concurrently optimize the molecular ground-state geometry and full wave function-Jastrow factor, orbitals, and configuration interaction coefficients-in variational Monte Carlo (VMC) for the prototypical case of 1,3- trans-butadiene, a small yet theoretically challenging π-conjugated system. We find that the CIPSI selection outperforms the conventional scheme of correlating orbitals within active spaces chosen by chemical intuition: it gives significantly better variational and diffusion Monte Carlo energies for all but the smallest expansions, and much smoother convergence of the geometry with the number of determinants. In particular, the optimal bond lengths and bond-length alternation of butadiene are converged to better than 1 mÅ with just a few thousand determinants, to values very close to the corresponding CCSD(T) results. The combination of CIPSI expansion and VMC optimization represents an affordable tool for the determination of accurate ground-state geometries in quantum Monte Carlo.
is investigated using truncated CIPSI-Jastrow CSF wave functions with Hartree-Fock orbitals within the framework of variational and diffusion quantum Monte Carlo. The truncation is performed based on the absolute value of the CI coefficients, and the Jastrow, molecular orbitals, and CI parameters are either partially or fully reoptimized with respect to the variational energy. Excellent absolute as well as bond dissociation energies are obtained at DMC level with very compact, fully optimized wave functions. By studying the expansions in more detail, we observe a change in the CI picture when reoptimizing the antisymmetric part of the CIPSI-Jastrow wave functions. Furthermore, we demonstrate that a decrease in the VMC energy as well as an improvement of the nodal surface quality can be achieved─with the same expansion size─if the CSFs are selected in the presence of a Jastrow correlation function, laying the foundation for a Jastrow selected CI scheme with quantum Monte Carlo.
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Þ
Modern quantum Monte Carlo (QMC) methods often capture electron correlation through both explicitly correlating Jastrow factors and small to mid-sized configuration interaction (CI) expansions. Here, we study the additional optimization difficulty created by including increasing numbers of CI parameters. We find evidence that the quality of Variational Mone Carlo (VMC) optimization can be limited by the ability to statistically resolve the CI parameters in the presence of a Jastrow factor. Although using larger statistical samples can mitigate this issue and bring an optimization closer to its true minimum, this approach is computationally intensive. We present evidence that similar gains to optimization quality can be had without increasing the sample size by avoiding CI parameters that are statistically the most difficult to resolve. Our findings suggest that, in addition to the value of using traditional selected configuration interaction (sCI) methods to prepare VMC wave functions, sCI-like methods can play an important role within VMC by improving the effectiveness of stochastic energy minimization.
We present a study of fixed- and partial-node approximations in Slater determinant basis sets, using full configuration interaction quantum Monte Carlo (FCIQMC) to perform sampling. Walker annihilation in the FCIQMC method allows partial-node simulations to be performed, relaxing the nodal constraint to converge to the FCI solution. This is applied to ab initio molecular systems, using symmetry-projected Jastrow mean-field wave functions for complete active space (CAS) problems. Convergence and the sign problem within the partial-node approximation are studied, which is shown to eventually be limited in its use due to the large walker populations required. However, the fixed-node approximation results in an accurate and practical method. We apply these approaches to various molecular systems and active spaces, including ferrocene and acenes. This also provides a test of symmetry-projected Jastrow mean-field wave functions in variational Monte Carlo for a new set of problems. For trans -polyacetylene molecules and acenes, we find that the time to perform a constant number of fixed-node FCIQMC iterations scales as and, respectively, resulting in an efficient method for CAS-based problems that can be applied accurately to large active spaces.
Multideterminant wavefunctions, while having a long history in quantum chemistry, are increasingly being used in highly accurate quantum Monte Carlo calculations. Since the accuracy of QMC is ultimately limited by the quality of the trial wavefunction, multi-Slater determinants wavefunctions offer an attractive alternative to Slater-Jastrow and more sophisticated wavefunction ansatz for several reasons. They can be efficiently calculated, straightforwardly optimized, and systematically improved by increasing the number of included determinants. In spite of their potential, however, the convergence properties of multi-Slater determinant wavefunctions with respect to orbital set choice and excited determinant selection are poorly understood, which hinders the application of these wavefunctions to large systems and solids. In this paper, by performing QMC calculations on the equilibrium and stretched carbon dimer, we find that convergence of the recovered correlation energy with respect to number of determinants can depend quite strongly on basis set and determinant selection methods, especially where there is strong correlation. We demonstrate that properly chosen orbital sets and determinant selection techniques from quantum chemistry methods can dramatically reduce the required number of determinants (and thus the computational cost) to reach a given accuracy, which we argue shows clear need for an automatic QMC-only method for selecting determinants and generating optimal orbital sets.
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.