Quantum Monte Carlo calculations of the potential energy curve of the helium dimer

We report results of two quantum Monte Carlo methods -- variational Monte Carlo and diffusion Monte Carlo -- on the potential energy curve of the helium dimer. In contrast to previous quantum Monte Carlo calculations on this system, we have employed trial wave functions of the Slater-Jastrow form and used the fixed node approximation for the fermion nodal surface. We find both methods to be in excellent agreement with the best theoretical results at short range. In addition, the diffusion Monte Carlo results give very good agreement across the whole potential energy curve, while the Slater-Jastrow wave function fails to bind the dimer at all.     

Computational methods in coupled electron-ion Monte Carlo simulations.

In the last few years, we have been developing a Monte Carlo simulation method to cope with systems of many electrons and ions in the Born-Oppenheimer approximation: the coupled electron-ion Monte Carlo method (CEIMC). Electronic properties in CEIMC are computed by quantum Monte Carlo rather than by density functional theory (DFT) based techniques. CEIMC can, in principle, overcome some of the limitations of the present DFT-based ab initio dynamical methods. The new method has recently been applied to high-pressure metallic hydrogen. Herein, we present a new sampling algorithm that we have developed in the framework of the reptation quantum Monte Carlo method chosen to sample the electronic degrees of freedom, thereby improving its efficiency. Moreover, we show herein that, at least for the case of metallic hydrogen, variational estimates of the electronic energies lead to an accurate sampling of the proton degrees of freedom.     

Computing the energy of a water molecule using multideterminants: a simple, efficient algorithm

Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.     

Variational Monte Carlo calculations for some cations and anions of the first‐row atoms using explicitly correlated wave functions

The variational Monte Carlo method is applied to calculate ground‐state energies of some cations and anions of the first‐row atoms. Accurate values providing between 80 and 90% of the correlation energy are obtained. Explicitly correlated wave functions including up to 42 variational parameters are used. The nondynamic correlation due to the 2s ? 2p near degeneracy effect is included by using a multideterminant wave function. The variational free parameters have been fixed by minimizing the energy that has shown to be a more convenient functional than the variance of the local energy, which is the most commonly employed method in variational Monte Carlo calculations. The energies obtained improve previous works using similar wave functions. © 2002 Wiley Periodicals, Inc.; DOI 10.1002/qua.10125     

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.     

Quantum Monte Carlo study of first-row atoms using transcorrelated variational Monte Carlo trial functions

The effect of using the transcorrelated variational Monte Carlo (TC-VMC) approach to construct a trial function for fixed node diffusion Monte Carlo (DMC) energy calculations has been investigated for the first-row atoms, Li to Ne. The computed energies are compared with fixed node DMC energies obtained using trial functions constructed from Hartree-Fock and density functional levels of theory. Despite major VMC energy improvement with TC-VMC trial functions, no improvement in DMC energy was observed using these trial functions for the first-row atoms studied. The implications of these results on the nodes of the trial wave functions are discussed.     

Quantum Monte Carlo study of the first-row atoms and ions

Quantum Monte Carlo calculations of the first-row atoms Li-Ne and their singly positively charged ions are reported. Multideterminant-Jastrow-backflow trial wave functions are used which recover more than 98% of the correlation energy at the variational Monte Carlo level and more than 99% of the correlation energy at the diffusion Monte Carlo level for both the atoms and ions. We obtain the first ionization potentials to chemical accuracy. We also report scalar relativistic corrections to the energies, mass-polarization terms, and one- and two-electron expectation values.     

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.     

Oscillator strengths of helium computed using Monte Carlo methods

We have optimized trial wave functions for the three lowest states of the helium atom with symmetry 1S, 1P, 1D, 3S, 3P, and 3D using variational Monte Carlo methods. With these wave functions we then computed dipole oscillator strengths for the 1S-1P, 1P-1D, 3S-3P, and 3P-3D transitions using the length, velocity, and acceleration forms. Our values are in good agreement with the best results found in the literature.     

应用Monte Carlo方法计算He原子包含电子相关作用的基态能量

应用Monte Carlo方法计算He原子包含电子相关波函数的基态能量,获得了与精确值非常接近的结果．实践表明,应用Monte Carlo方法有可能在多电子体系中直接采用包含任意2个电子间距离ry的函数作为变分函数来考虑电子相关作用．     

Electron pair localization function: a practical tool to visualize electron localization in molecules from quantum Monte Carlo data

In this work we introduce an electron localization function describing the pairing of electrons in a molecular system. This function, called "electron pair localization function," is constructed to be particularly simple to evaluate within a quantum Monte Carlo framework. Two major advantages of this function are the following: (i) the simplicity and generality of its definition; and (ii) the possibility of calculating it with quantum Monte Carlo at various levels of accuracy (Hartree-Fock, multiconfigurational wave functions, valence bond, density functional theory, variational Monte Carlo with explicitly correlated trial wave functions, fixed-node diffusion Monte Carlo, etc). A number of applications of the electron pair localization function to simple atomic and molecular systems are presented and systematic comparisons with the more standard electron localization function of Becke and Edgecombe are done. Results illustrate that the electron pair localization function is a simple and practical tool for visualizing electronic localization in molecular systems.     

Vibrational–rotational energies of all H2 isotopomers using Monte Carlo methods

Using variational Monte Carlo techniques, we have computed several of the lowest rotational–vibrational energies of all the hydrogen molecule isotopomers (H₂, HD, HT, D₂, DT, and T₂). These calculations do not require the excited states to be explicitly orthogonalized. We have examined both the usual Gaussian wave function form as well as a rapidly convergent Padé form. The high‐quality potential energy surfaces used in these calculations are taken from our earlier work and include the Born–Oppenheimer energy, the diagonal correction to the Born–Oppenheimer approximation, and the lowest‐order relativistic corrections at 24 internuclear points. Our energies are in good agreement with those determined by other methods. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Full optimization of Jastrow-Slater wave functions with application to the first-row atoms and homonuclear diatomic molecules

We pursue the development and application of the recently introduced linear optimization method for determining the optimal linear and nonlinear parameters of Jastrow-Slater wave functions in a variational Monte Carlo framework. In this approach, the optimal parameters are found iteratively by diagonalizing the Hamiltonian matrix in the space spanned by the wave function and its first-order derivatives, making use of a strong zero-variance principle. We extend the method to optimize the exponents of the basis functions, simultaneously with all the other parameters, namely, the Jastrow, configuration state function, and orbital parameters. We show that the linear optimization method can be thought of as a so-called augmented Hessian approach, which helps explain the robustness of the method and permits us to extend it to minimize a linear combination of the energy and the energy variance. We apply the linear optimization method to obtain the complete ground-state potential energy curve of the C(2) molecule up to the dissociation limit and discuss size consistency and broken spin-symmetry issues in quantum Monte Carlo calculations. We perform calculations for the first-row atoms and homonuclear diatomic molecules with fully optimized Jastrow-Slater wave functions, and we demonstrate that molecular well depths can be obtained with near chemical accuracy quite systematically at the diffusion Monte Carlo level for these systems.     

Rovibrationally averaged properties of H2 using Monte Carlo methods

Using variational Monte Carlo and a simple explicitly correlated wave function, we have computed 18 molecular properties of the hydrogen molecule (X¹∑) at 24 internuclear distances. These properties have been combined with rapidly convergent rovibrational wave functions to produce rovibrationally averaged properties for several of the lowest rotational and vibrational levels of this system. Our results are in very good agreement with previous values found in the literature. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Optimum and efficient sampling for variational quantum Monte Carlo

Quantum mechanics for many-body systems may be reduced to the evaluation of integrals in 3N dimensions using Monte Carlo, providing the Quantum Monte Carlo ab initio methods. Here we limit ourselves to expectation values for trial wave functions, that is to variational quantum Monte Carlo. Almost all previous implementations employ samples distributed as the physical probability density of the trial wave function, and assume the central limit theorem to be valid. In this paper we provide an analysis of random error in estimation and optimization that leads naturally to new sampling strategies with improved computational and statistical properties. A rigorous lower limit to the random error is derived, and an efficient sampling strategy presented that significantly increases computational efficiency. In addition the infinite variance heavy tailed random errors of optimum parameters in conventional methods are replaced with a Normal random error, strengthening the theoretical basis of optimization. The method is applied to a number of first row systems and compared with previously published results.     

Quantum Monte Carlo for 3d transition-metal atoms

The domain Green's function Monte Carlo method has been used to calculate the ground-state energy of the atoms Sc through Zn. The fixed node approximation with single-configuration explicitly correlated wave functions is used. A comparison with variational Monte Carlo energies is carried out. The quality of the ground-state energies reported here is similar to that achieved for few-electron atoms using similar techniques.     

Excited electronic state calculations by the transcorrelated variational Monte Carlo method: application to a helium atom

We have implemented the excited electronic state calculations for a helium atom by the transcorrelated variational Monte Carlo (TC-VMC) method. In this method, Jastrow-Slater-type wave function is efficiently optimized not only for the Jastrow factor but also for the Slater determinant. Since the formalism for the TC-VMC method is based on the variance minimization, excited states as well as the ground state calculations are feasible. It is found that both the first and the second excitation energies given by TC-VMC are much closer to the experimental data than those given by the variational Monte Carlo method with using the Hartree-Fock orbitals. The successful results in the TC-VMC method are considered to be due to the nodal optimization of the wave functions.     

Insulator to metal transition in fluid deuterium

We have investigated the insulator to metal transition in fluid deuterium using first principles simulations. Both density functional and quantum Monte Carlo calculations indicate that the electronic energy gap of the liquid vanishes at about ninefold compression and 3000 K. At these conditions the computed conductivity values are characteristic of a poor metal. These findings are consistent with those of recent shock wave experiments but the computed conductivity is larger than the measured value. From our ab initio results we conclude that the transition is driven by molecular dissociation rather than disorder and that both temperature and pressure play a key role in determining structural changes in the fluid.     

Properties of selected diatomics using variational Monte Carlo methods

Using variational Monte Carlo and highly accurate trial wave functions optimized by Filippi and Umrigar, we calculate a number of molecular properties for the ground state of Li2, Be2, B2, C2, N2, O2, and F2. This is the first time that many of these properties have been computed.     

A practical treatment for the three-body interactions in the transcorrelated variational Monte Carlo method: application to atoms from lithium to neon

We suggest a practical solution to dealing with the three-body interactions in the transcorrelated variational Monte Carlo method (TC-VMC). In the TC-VMC method, which was suggested in our previous paper [N. Umezawa and S. Tsuneyuki, J. Chem. Phys. 119, 10015 (2003)], the Jastrow-Slater-type wave function is efficiently optimized through a self-consistent procedure by minimizing the variance of the local energy. The three-body terms in the transcorrelated self-consistent-field equation, which have been simply ignored in our previous works, are efficiently calculated by the Monte Carlo numerical integration. We found that our treatment for the three-body interactions is successful for atoms from Li to Ne.