Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDissociation energy of the water dimer from Quantum Monte Carlo calculations
101
0
0.0
(
0
)
Added by
Idoia G. de Gurtubay
Publication date
2007
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We report a study of the electronic dissociation energy of the water dimer using quantum Monte Carlo (QMC) techniques. We have performed variational quantum Monte Carlo (VMC) 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-like 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 mHa, yielding a DMC energy of -76.42830(5) Ha, which is only 10 mHa 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) kcal/mol and are in good agreement with experiment. Our calculations give monomer dipole moments which range between 1.897(2) and 1.909(4) Debye and dimer dipole moments which range between 2.628(6) and 2.672(5) Debye.
rate research
Read More
We have used diffusion Monte Carlo (DMC) simulations to calculate the energy barrier for H$_2$ dissociation on the Mg(0001) surface. The calculations employ pseudopotentials and systematically improvable B-spline basis sets to expand the single particle orbitals used to construct the trial wavefunctions. Extensive tests on system size, time step, and other sources of errors, performed on periodically repeated systems of up to 550 atoms, show that all these errors together can be reduced to $sim 0.03$ eV. The DMC dissociation barrier is calculated to be $1.18 pm 0.03$ eV, and is compared to those obtained with density functional theory using various exchange-correlation functionals, with values ranging between 0.44 and 1.07 eV.
An ab-initio method for determining the dynamical structure function of an interacting many--body quantum system has been devised by combining a generalized integral transform method with Quantum Monte Carlo methods. As a first application, the coherent and, separately, the incoherent excitation spectrum of bulk atomic 4He has been computed, both in the low and intermediate momentum range. The peculiar form of the kernel in the integral transform of the dynamical structure function allows to predict, without using any model, both position and width of the collective excitations in the maxon--roton region, as well as the second collective peak. A prediction of the dispersion of the single--particle modes described by the incoherent part is also presented.
Extended solids are frequently simulated as finite systems with periodic boundary conditions, which due to the long-range nature of the Coulomb interaction may lead to slowly decaying finite- size errors. In the case of Quantum-Monte-Carlo simulations, which are based on real space, both real-space and momentum-space solutions to this problem exist. Here, we describe a hybrid method which using real-space data models the spherically averaged structure factor in momentum space. We show that (i) by integration our hybrid method exactly maps onto the real-space model periodic Coulomb-interaction (MPC) method and (ii) therefore our method combines the best of both worlds (real-space and momentum-space). One can use known momentum-resolved behavior to improve convergence where MPC fails (e.g., at surface-like systems). In contrast to pure momentum-space methods, our method only deals with a simple single-valued function and, hence, better lends itself to interpolation with exact small-momentum data as no directional information is needed. By virtue of integration, the resulting finite-size corrections can be written as an addition to MPC.
We report results of both Diffusion Quantum Monte Carlo(DMC) method and Reptation Quantum Monte Carlo(RMC) method on the potential energy curve of the helium dimer. We show that it is possible to obtain a highly accurate description of the helium dimer. An improved stochastic reconfiguration technique is employed to optimize the many-body wave function, which is the starting point for highly accurate simulations based on the Diffusion Quantum Monte Carlo(DMC) and Reptation Quantum Monte Carlo (RMC) methods. We find that the results of these methods are in excellent agreement with the best theoretical results at short range, especially recently developed Reptation Quantum Monte Carlo(RMC) method, yield practically accurate results with reduced statistical error, which gives very excellent agreement across the whole potential. For the equilibrium internuclear distance of 5.6 bohr, the calculated electronic energy with Reptation Quantum Monte Carlo(RMC) method is 5.807483599$pm$0.000000015 hartrees and the corresponding well depth is -11.003$pm$0.005 K.
We report an accurate study of interactions between Benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory (DFT) using different van der Waals (vdW) functionals. In our QMC calculations, we use accurate correlated trial wave functions including three-body Jastrow factors, and backflow transformations. We consider two benzene molecules in the parallel displaced (PD) geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the CCSD(T)/CBS limit is -2.65(2) kcal/mol [E. Miliordos et al, J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
