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.     

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.     

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     

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.     

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.     

量子Monte Carlo处理激发态

提出了用于电子激发态的剩余函数变分量子ＭｏｎｔｅＣａｒｌｏ（ＳＦＶＭＣ）方法,已经证明：若激发态的初始波函数与基态的初始波函数属于对称性不同的不可约表示时,该激发态的ＳＦＶＭＣ方法与基态的ＳＦＶＭＣ方法完全相同;若激发态的初始波函数与基态的初始波函数有相同的对称性时,只要对激发态的初始波函数作正交性修正,则其态的ＳＦＶＭＣ方法亦可推到该激发态的情况。文章导出了这第二类激发态的ＳＦＶＭＣ方法的详细计     

One- and two-body densities of carbon isoelectronic series in their low-lying multiplet states from explicitly correlated wave functions

The (3)P ground state and both the (1)D and (1)S excited states arising from the low-lying 1s(2)2s(2)2p(2) configuration of the carbon isoelectronic series are studied starting from explicitly correlated multiconfigurational wave functions. One- and two-body densities in position space have been calculated and different one- and two-body expectation values have been obtained. The effects of electronic correlations have been systematically studied. All the calculations have been done by means of variational Monte Carlo.     

Quantum Monte Carlo for electronic excitations of free-base porphyrin

Accurate calculations of allowed and nonallowed transitions in porphyrin are reported. Using the quantum Monte Carlo method in the diffusion Monte Carlo variant, the vertical transition between the ground state singlet and the second excited state singlet as well as the adiabatic transition between the ground state and the lowest triplet state have been computed for this 162-electron system. The present theoretical results are compared to experiment and to results of other theoretical methods. The diffusion Monte Carlo energy differences are found to be in excellent agreement with experiment.     

Energies of the first row atoms from quantum Monte Carlo

All-electron variational and diffusion quantum Monte Carlo calculations of the ground state energies of the first row atoms (from Li to Ne) are reported. The authors use trial wave functions of four types: single-determinant Slater-Jastrow wave functions, multideterminant Slater-Jastrow wave functions, single-determinant Slater-Jastrow wave functions with backflow transformations, and multideterminant Slater-Jastrow wave functions with backflow transformations. At the diffusion quantum Monte Carlo level and using their multideterminant Slater-Jastrow wave functions with backflow transformations, they recover 99% or more of the correlation energies for Li, Be, B, C, N, and Ne, 97% for O, and 98% for F.     

Excited states of boron isoelectronic series from explicitly correlated wave functions

The ground state and some low-lying excited states arising from the 1s2 2s2p2 configuration of the boron isoelectronic series are studied starting from explicitly correlated multideterminant wave functions. One- and two-body densities in position space have been calculated and different expectation values such as , , , , , and , where r, r12, and R stand for the electron-nucleus, interelectronic, and two electron center of mass coordinates, respectively, have been obtained. The energetic ordering of the excited states and the fulfillment of the Hund's rules is analyzed systematically along the isoelectronic series in terms of the electron-electron and electron-nucleus potential energies. The effects of electronic correlations have been systematically studied by comparing the correlated results with the corresponding noncorrelated ones. All the calculations have been done by using the variational Monte Carlo method.     

Jastrow correlated and quantum Monte Carlo calculations for the low-lying states of the carbon atom

Different computational methods are employed to calculate excitation energies of the carbon atom. Explicitly correlated wave functions have been obtained in a Variational Monte Carlo calculation. Fixed node Diffusion Monte Carlo calculations for the lowest energy excited states of a given symmetry are reported. A systematic and quantitative analysis of the performance of the different schemes in the calculation of the excitation energy of up to 27 excited states of the carbon atom is carried out. The quality of the different methods have been studied in terms of the deviation with respect to the experimental excitation energies. A good agreement with the experimental values has been reached.     

Quantum Monte Carlo with Jastrow-valence-bond wave functions

We consider the use in quantum Monte Carlo calculations of two types of valence bond wave functions based on strictly localized active orbitals, namely valence bond self-consistent-field and breathing-orbital valence bond wave functions. Complemented by a Jastrow factor, these Jastrow-valence-bond wave functions are tested by computing the equilibrium well depths of the four diatomic molecules C(2), N(2), O(2), and F(2) in both variational Monte Carlo and diffusion Monte Carlo. We show that it is possible to design compact wave functions based on chemical grounds that are capable of describing both static and dynamic electron correlations. These wave functions can be systematically improved by inclusion of valence bond structures corresponding to additional bonding patterns.     

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.     

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.     

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.     

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

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

Path integral ground state study of finite-size systems: application to small (parahydrogen)N (N=2-20) clusters

We use the path integral ground state method to study the energetic and structural properties of small para-H2 clusters of sizes ranging from 2 to 20 molecules. A fourth order formula is used to approximate the short imaginary-time propagator and two interaction potentials are considered. Our results are compared to those of exact basis set calculations and other quantum Monte Carlo methods when available. We find that for all cluster sizes considered, our results show a lower ground state energy than literature values obtained by diffusion Monte Carlo and variational Monte Carlo. For the dimer and trimer, ground state energies are in good agreement with exact results obtained using the discrete variable representation. Structural properties are found to be insensitive to the choice of interaction potential. We explore the use of Pekeris coordinates to analyze the importance of linear arrangement in trimers and for trimers within clusters of larger size.     

Quantum Monte Carlo calculations for ground and excited states

A brief overview of the diffusion quantum Monte Carlo method is given. We illustrate the application to ground‐state calculations by a study of the relative stability of carbon clusters near the crossover to fullerene stability, thereby determining the smallest stable fullerene. The application to excited states is illustrated via a study of excitonic states in small hydrogenated silicon clusters. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

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     

A ground state potential energy surface for H2 using Monte Carlo methods

Using variational Monte Carlo and a simple explicitly correlated wave function we have computed the Born-Oppenheimer energy of the H2 ground state (X 1Sigmag+) at 24 internuclear distances. We have also calculated the diagonal correction to the Born-Oppenheimer approximation and the lowest-order relativistic corrections at each distance using variational Monte Carlo techniques. The nonadiabatic values are evaluated from numerical derivatives of the wave function with respect to the nuclear coordinates. With this potential energy surface we have computed several of the lowest vibrational-rotational energies for this system. Our results are in good agreement with the best values found in the literature.