Binding energies of first row diatomics in the light of the interacting quantum atoms approach

Binding energies of first row diatomics are revisited within the interacting quantum atoms (IQA) approach. This is a formalism in chemical bonding theory based upon the quantum theory of atoms in molecules. It is characterized by the preservation of the energetic identity of atoms within molecules. Quantum mechanically computed binding energies are recovered in IQA as a sum of small atomic deformation energies and large pairwise interaction terms. We show how this partition responds faithfully to chemical intuition, and how the different evolution of deformations and interactions accounts in a unified manner for the subtle variations of the binding energy of these molecules.     

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.     

Dissociation energy of the water dimer from quantum Monte Carlo calculations

We report a study of the electronic dissociation energy of the water dimer using quantum Monte Carlo techniques. We have performed variational quantum Monte Carlo and diffusion quantum Monte Carlo (DMC) calculations of the electronic ground state of the water monomer and dimer using all-electron and pseudopotential approaches. We have used Slater-Jastrow trial wave functions with B3LYP type single-particle orbitals, into which we have incorporated backflow correlations. When backflow correlations are introduced, the total energy of the water monomer decreases by about 4-5 mhartree, yielding a DMC energy of -76.428 30(5) hartree, which is only 10 mhartree above the experimental value. In our pseudopotential DMC calculations, we have compared the total energies of the water monomer and dimer obtained using the locality approximation with those from the variational scheme recently proposed by Casula [Phys. Rev. B 74, 161102(R) (2006)]. The time step errors in the Casula scheme are larger, and the extrapolation of the energy to zero time step always lies above the result obtained with the locality approximation. However, the errors cancel when energy differences are taken, yielding electronic dissociation energies within error bars of each other. The dissociation energies obtained in our various all-electron and pseudopotential calculations range between 5.03(7) and 5.47(9) kcalmol and are in good agreement with experiment. Our calculations give monomer dipole moments which range between 1.897(2) and 1.909(4) D and dimer dipole moments which range between 2.628(6) and 2.672(5) D.     

DFT-LDA pseudopotentials in quantum Monte Carlo

We investigate the portability of standard norm-conserving pseudopotentials outside the density functional theory-local density approximation (DFT-LDA) framework, i.e., their use and interpretation as electron-ion effective potentials in valence-only diffusion Monte Carlo simulations. While first-principles many-body pseudopotentials are not available in the literature yet, the use of approximate pseudopotentials in quantum Monte Carlo simulations is becoming widespread. Here we attempt a systematic analysis of the portability of norm-conserving pseudopotentials generated within DFT-LDA, focusing on a model many-body system, the two-electron valence-only ion. Our results indicate that the portability is good in most cases, hence the use of pseudopotentials in quantum Monte Carlo simulations is in general a reasonable approximation but suggest that in some cases this approximation may be relevant. © 1997 John Wiley & Sons, Inc.     

Prospects for release-node quantum Monte Carlo

We perform release-node quantum Monte Carlo simulations on the first row diatomic molecules in order to assess how accurately their ground-state energies can be obtained. An analysis of the fermion-boson energy difference is shown to be strongly dependent on the nuclear charge, Z, which in turn determines the growth of variance of the release-node energy. It is possible to use maximum entropy analysis to extrapolate to ground-state energies only for the low Z elements. For the higher Z dimers beyond boron, the error growth is too large to allow accurate data for long enough imaginary times. Within the limit of our statistics we were able to estimate, in atomic units, the ground-state energy of Li(2) (-14.9947(1)), Be(2) (-29.3367(7)), and B(2)(-49.410(2)).     

Rydberg states with quantum Monte Carlo

Calculations on Rydberg states are performed using quantum Monte Carlo methods. Excitation energies and singlet-triplet splittings are calculated for two model systems, the carbon atom (3P and 1P) and carbon monoxide ((1Sigma and 3Sigma). Kohn-Sham wave functions constructed from open-shell localized Hartree-Fock orbitals are used as trial and guide functions. The fixed-node diffusion quantum Monte Carlo (FN-DMC) method depends strongly on the wave function's nodal hypersurface. Nodal artefacts are investigated for the ground state of the carbon atom. Their effect on the FN-DMC results can be analyzed quantitatively. FN-DMC leads to accurate excitation energies but to less accurate singlet-triplet splittings. Variational Monte Carlo calculations are able to reproduce the experimental results for both the excitation energies and the singlet-triplet splittings.     

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.     

Time-dependent quantum Monte Carlo and the stochastic quantization

We examine the relation between the recently proposed time-dependent quantum Monte Carlo (TDQMC) method and the principles of stochastic quantization. In both TDQMC and stochastic quantization, particle motion obeys stochastic guidance equations to preserve quantum equilibrium. In this way the probability density of the Monte Carlo particles corresponds to the modulus square of the many-body wave function at all times. However, in TDQMC, the motion of particles and guide waves occurs in physical space unlike in stochastic quantization where it occurs in configuration space. Hence, the practical calculation of time evolution of many-body fully correlated quantum systems becomes feasible within the TDQMC methodology. We illustrate the TDQMC technique by calculating the symmetric and antisymmetric ground state of a model one-dimensional helium atom, and the time evolution of the dipole moment when the atom is irradiated by a strong ultrashort laser pulse.     

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.     

Gaussian basis sets for the first and second row atoms

Gaussian basis sets consisting for first row atoms of 7 s-type and 3 p-type and for second row atoms of 10 s-type and 6 p-type functions with optimized exponents are reported. These basis sets consists of at least two functions per atomic orbital.

---

Zusammenfassung Es werden für die Atome der ersten und zweiten Reihe Basissätze aus Gaußfunktionen mitgeteilt, die aus 7 Funktionen vom s-Typ und 3 Funktionen vom p-Typ für die Elemente der ersten Reihe und 10 Funktionen vom s-Typ und 6 Funktionen vom p-Typ für die Elemente der zweiten Reihe mit optimierten Exponenten bestehen. Diese Basissätze bestehen aus wenigstens zwei Funktionen pro Atomorbital.

---

Résumé Une base de 7 gaussiens du type s et 3 du type p est presenté pour les éléments du premier rang et de 10 gaussiens du type s et 6 du type p pour des éléments du deuxième rang; les exposants sont optimisés. Les bases consistent au moins en deux fonctions par orbital atomique.

     

Dynamic correlations with time-dependent quantum Monte Carlo

In this paper, we solve quantum many-body problem by propagating ensembles of trajectories and guiding waves in physical space. We introduce the "effective potential" correction within the recently proposed time-dependent quantum Monte Carlo methodology to incorporate the nonlocal quantum correlation effects between the electrons. The associated correlation length is calculated by adaptive kernel density estimation over the walker distribution. The general formalism is developed and tested on one-dimensional helium atom in laser field of different intensities and carrier frequencies. Good agreement with exact results for the atomic ionization is obtained.     

Vanadium oxide compounds with quantum Monte Carlo

Calculations with the diffusion quantum Monte Carlo method are presented for vanadium oxide molecules VO0/+0(n) with n = 1-4 and for V2O5. Atomization and ionization energies are calculated as well as oxygen abstraction energies. The fixed-node approximation is compared for guide functions with orbitals from B3LYP and BP86 calculations and higher accuracy was obtained with the latter orbitals. Additionally, all-electron and pseudopotential calculations are compared for the oxygen atom. The overall accuracy is found to be comparable to CCSD(T) calculations where experimental data is available.     

Approaching chemical accuracy with quantum Monte Carlo

A quantum Monte Carlo study of the atomization energies for the G2 set of molecules is presented. Basis size dependence of diffusion Monte Carlo atomization energies is studied with a single determinant Slater-Jastrow trial wavefunction formed from Hartree-Fock orbitals. With the largest basis set, the mean absolute deviation from experimental atomization energies for the G2 set is 3.0 kcal/mol. Optimizing the orbitals within variational Monte Carlo improves the agreement between diffusion Monte Carlo and experiment, reducing the mean absolute deviation to 2.1 kcal/mol. Moving beyond a single determinant Slater-Jastrow trial wavefunction, diffusion Monte Carlo with a small complete active space Slater-Jastrow trial wavefunction results in near chemical accuracy. In this case, the mean absolute deviation from experimental atomization energies is 1.2 kcal/mol. It is shown from calculations on systems containing phosphorus that the accuracy can be further improved by employing a larger active space.     

Energy-consistent pseudopotentials for quantum Monte Carlo calculations

The authors present scalar-relativistic energy-consistent Hartree-Fock pseudopotentials for the main-group elements. The pseudopotentials do not exhibit a singularity at the nucleus and are therefore suitable for quantum Monte Carlo (QMC) calculations. They demonstrate their transferability through extensive benchmark calculations of atomic excitation spectra as well as molecular properties. In particular, they compute the vibrational frequencies and binding energies of 26 first- and second-row diatomic molecules using post-Hartree-Fock methods, finding excellent agreement with the corresponding all-electron values. They also show their pseudopotentials give superior accuracy than other existing pseudopotentials constructed specifically for QMC. Finally, valence basis sets of different sizes (VnZ with n=D,T,Q,5 for first and second rows, and n=D,T for third to fifth rows) optimized for our pseudopotentials are also presented.     

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     

Bond breaking with auxiliary-field quantum Monte Carlo

Bond stretching mimics different levels of electron correlation and provides a challenging test bed for approximate many-body computational methods. Using the recently developed phaseless auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well-studied molecules BH and N(2) and in the H(50) chain. To control the sign/phase problem, the phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. With single Slater determinants from unrestricted Hartree-Fock as trial wave function, the phaseless AF QMC method generally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) method in mapping the potential-energy curve. In both BH and N(2), we also study the use of multiple-determinant trial wave functions from multiconfiguration self-consistent-field calculations. The increase in computational cost versus the gain in statistical and systematic accuracy are examined. With such trial wave functions, excellent results are obtained across the entire region between equilibrium and the dissociation limit.     

Spectroscopic data for the LiH molecule from pseudopotential quantum Monte Carlo calculations

Quantum Monte Carlo and quantum chemistry techniques are used to investigate pseudopotential models of the lithium hydride (LiH) molecule. Interatomic potentials are calculated and tested by comparing with the experimental spectroscopic constants and well depth. Two recently developed pseudopotentials are tested, and the effects of introducing a Li core polarization potential are investigated. The calculations are sufficiently accurate to isolate the errors from the pseudopotentials and core polarization potential. Core-valence correlation and core relaxation are found to be important in determining the interatomic potential.     

Diffusion Monte Carlo Calculations of Zero-Point Energies of Methanol and Deuterated Methanol

Diffusion Monte Carlo (DMC) simulations have been used to obtain quantum zero-point energies of methanol and all its isotopologs and isotopomers, using a new, accurate semi-global potential energy surface. This potential energy surface is a precise, permutationally invariant fit to 6676 ab initio energies, obtained at the CCSD(T)-F12b/aug-cc-pVDZ level of theory. Quantum zero-point energies of deuterated methanol isotopomers are very close to each other and so a simple statistical argument can be used to estimate the populations of each isotopomer at very low-temperatures. The DMC simulations also indicate that there is virtually zero probability for H/D exchange in the zero-point state. © 2019 Wiley Periodicals, Inc.