Fixed-node Diffusion Monte Carlo Energies Research Articles

Accurate Atomic Correlation and Total Energies for Correlation Consistent Effective Core Potentials.

Very recently, we introduced a set of correlation consistent effective core potentials (ccECPs) constructed within full many-body approaches. By employing significantly more accurate correlated approaches, we were able to reach a new level of accuracy for the resulting effective core Hamiltonians. We also strived for simplicity of use and easy transferability into a variety of electronic structure methods in quantum chemistry and condensed matter physics. Here, as a reference for future use, we present exact or nearly exact total energy calculations for these ccECPs. The calculations cover H-Kr elements and are based on the state-of-the-art configuration interaction (CI), coupled-cluster (CC), and quantum Monte Carlo (QMC) calculations with systematically eliminated/improved errors. In particular, we carry out full CI/CCSD(T)/CCSDT(Q) calculations with cc-pVnZ with up to n = 6 basis sets and we estimate the complete basis set limits. Using combinations of these approaches, we achieved an accuracy of ≈1-10 mHa for K-Zn atoms and ≈0.1-0.3 mHa for all other elements-within about 1% or better of the ccECP total correlation energies. We also estimate the corresponding kinetic energies within the feasible limit of full CI calculations. In order to provide data for QMC calculations, we include fixed-node diffusion Monte Carlo energies for each element that give quantitative insights into the fixed-node biases for single-reference trial wave functions. The results offer a clear benchmark for future high-accuracy calculations in a broad variety of correlated wave function methods such as CI and CC as well is in stochastic approaches such as real space sampling QMC.

Non-orthogonal determinants in multi-Slater-Jastrow trial wave functions for fixed-node diffusion Monte Carlo.

The accuracy and efficiency of ab initio Quantum Monte Carlo (QMC) algorithms benefit greatly from compact variational trial wave functions that accurately reproduce ground state properties of a system. We investigate the possibility of using multi-Slater-Jastrow trial wave functions with non-orthogonal determinants by optimizing identical single particle orbitals independently in separate determinants. As a test case, we compute variational and fixed-node diffusion Monte Carlo (FN-DMC) energies of a C2 molecule. For a given multi-determinant expansion, we find that this non-orthogonal orbital optimization results in a consistent improvement in the variational energy and the FN-DMC energy on the order of a few tenths of an eV. In some cases, fewer non-orthogonal determinants are required compared to orthogonal ones in order to achieve similar accuracy in FN-DMC. Our calculations indicate that trial wave functions with non-orthogonal determinants can improve computed energies in a QMC calculation when compared to their orthogonal counterparts.

Noncovalent Interactions by Fixed-Node Diffusion Monte Carlo: Convergence of Nodes and Energy Differences vs Gaussian Basis-Set Size

Convergence of fixed-node (FN) shape and FN diffusion Monte Carlo (FNDMC) interaction energies is studied vs the Gaussian basis set saturation level in HF and CH4 dimers and one-determinant Slater-Jastrow trial wave functions (ΨT). The tested 25 distinct basis sets obtained by stepwise trimming of aug-VDZ and aug-VTZ bases suggest minimum basis set requirements to achieve reasonable results. A single selected trimmed basis set, about 2 times smaller in size than aug-VTZ, is extensively tested on a set of 12 noncovalent complexes including formic acid dimer, benzene-methane, or coronene-H2. The results indicate that equivalent noncovalent FNDMC energy differences are available at costs lower than assumed before. Additional insights from electron density differences and comparison of dimer vs monomer ΨT nodes explain this observation.

Bias cancellation in one-determinant fixed-node diffusion Monte Carlo: Insights from fermionic occupation numbers.

The accuracy of the fixed-node diffusion Monte Carlo (FNDMC) depends on the node location of the supplied trial state Ψ_{T}. The practical FNDMC approaches available for large systems rely on compact yet effective Ψ_{T}, most often containing an explicitly correlated single Slater determinant (SD). However, SD nodes may be better suited to one system than to another, which may possibly lead to inaccurate FNDMC energy differences. It remains a challenge how to estimate nonequivalence or appropriateness of SDs. Here we use the differences of a measure based on the Euclidean distance between the natural orbital occupation number (NOON) vector of the SD and the exact solution in the NOON vector space, which can be viewed as a measure of SD nonequivalence and as a qualitative measure of the expected degree of nondynamic-correlation-related bias in FNDMC energy differences. This is explored on a set of small noncovalent complexes and covalent bond breaking of Si_{2} vs N_{2}. It turns out that NOON-based measures well reflect the magnitude and sign of the bias present in the data available, thus providing insights into the nature of bias cancellation in SD FNDMC energy differences.

Fixed-node diffusion Monte Carlo potential energy curve of the fluorine molecule F2 using selected configuration interaction trial wavefunctions.

The potential energy curve of the F2 molecule is calculated with Fixed-Node Diffusion Monte Carlo (FN-DMC) using Configuration Interaction (CI)-type trial wavefunctions. To keep the number of determinants reasonable and thus make FN-DMC calculations feasible in practice, the CI expansion is restricted to those determinants that contribute the most to the total energy. The selection of the determinants is made using the CIPSI approach (Configuration Interaction using a Perturbative Selection made Iteratively). The trial wavefunction used in FN-DMC is directly issued from the deterministic CI program; no Jastrow factor is used and no preliminary multi-parameter stochastic optimization of the trial wavefunction is performed. The nodes of CIPSI wavefunctions are found to reduce significantly the fixed-node error and to be systematically improved upon increasing the number of selected determinants. To reduce the non-parallelism error of the potential energy curve, a scheme based on the use of a R-dependent number of determinants is introduced. Using Dunning's cc-pVDZ basis set, the FN-DMC energy curve of F2 is found to be of a quality similar to that obtained with full configuration interaction/cc-pVQZ.

Diffusion Monte Carlo method with symmetry-constrained variational principle for degenerate excited states

Fixed-node diffusion Monte Carlo (FN-DMC) is an accurate and useful method for estimating the wave function and the energy of the quantum ground state of a many-fermion system. However, it has been shown that difficulties with the method may occur when it is applied to a degenerate excited state because the nodal surface of the degenerate trial function is generally insufficient to impose the complete symmetry properties of the trial function on the FN-DMC wave function. As a result, the tiling theorem and the symmetry-constrained variational principle may be violated by FN-DMC when the excited state is degenerate. There are two practical consequences for the study of degenerate excited states: The FN-DMC energy may lie below the energy of the lowest stationary state that transforms according to the same degenerate irreducible representation as the trial function; and the convergence of the FN-DMC energy with improvements in the trial function may not be quadratic. In this paper a diffusion Monte Carlo method for degenerate excited states is presented. It provides a direct generalization of the FN-DMC method, and when applied to the study of degenerate excited states, it has the support of the tiling theorem and the symmetry-constrained variational principle. The method is applied to the lowest degenerate state of a simple test problem in which FN-DMC has been shown to violate both the tiling theorem and the symmetry-constrained variational principle. The numerical results support the assertion that this method for degenerate excited states satisfies both the tiling theorem and the symmetry-constrained variational principle.

Wave functions for quantum Monte Carlo calculations in solids: Orbitals from density functional theory with hybrid exchange-correlation functionals

We investigate how the fixed-node diffusion Monte Carlo energy of solids depends on single-particle orbitals used in Slater--Jastrow wave functions. We demonstrate that the dependence can be significant, in particular in the case of 3d transition-metal compounds, which we adopt as examples. We illustrate how exchange-correlation functionals with variable exact-exchange component can be exploited to reduce the fixed-node errors. On the basis of these results we argue that the fixed-node quantum Monte Carlo provides a variational approach for optimization of effective hamiltonians with parameters.

Accuracy of electronic wave functions in quantum Monte Carlo: The effect of high-order correlations

Compact and accurate wave functions can be constructed by quantum Monte Carlo methods. Typically, these wave functions consist of a sum of a small number of Slater determinants multiplied by a Jastrow factor. In this paper we study the importance of including high-order, nucleus-three-electron correlations in the Jastrow factor. An efficient algorithm based on the theory of invariants is used to compute the high-body correlations. We observe significant improvements in the variational Monte Carlo energy and in the fluctuations of the local energies but not in the fixed-node diffusion Monte Carlo energies. Improvements for the ground states of physical, fermionic atoms are found to be smaller than those for the ground states of fictitious, bosonic atoms, indicating that errors in the nodal surfaces of the fermionic wave functions are a limiting factor.