Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDissecting the hydrogen bond: a Quantum Monte Carlo approach
174
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a variational MonteCarlo (VMC) and lattice regularized diffusion MonteCarlo (LRDMC) study of the binding energy and dispersion curve of the water dimer. As a variation ansatz we use the JAGP wave function, an implementation of the resonating valence bond (RVB) idea. Actually one the aim of the present work is to investigate how the bonding of two water molecules, as a prototype of the hydrogen-bonded complexes, could be described within an JAGP approach. Using a pseudopotential for the inert core of the Oxygen, with a full optimization of the variational parameters, we obtain at the VMC level a binding energy of -4.5(0.1) Kcal/mol, while LRDMC calculations gives -4.9(0.1) Kcal/mol (experiment 5 Kcal/Mol). The calculated dispersion curve reproduces both at the VMC and LRDMC level the miminum position and the curvature.The quality of the WF gives us the possibility to dissect the binding energy in different contributions by appropriately switching off determinantal and Jastrow terms in the JAGP: we estimate the dynamical contribution to the binding energy to be of the order of 1.4(0.2) Kcal/Mol whereas the covalent contribution about 1.0(0.2) Kcal/Mol. JAGP reveales thus a promising WF for describing systems where both dispersive and covalent forces play an important role.
rate research
Read More
Path integral Monte Carlo approach is used to study the coupled quantum dynamics of the electron and nuclei in hydrogen molecule ion. The coupling effects are demonstrated by comparing differences in adiabatic Born--Oppenheimer and non-adiabatic simulations, and inspecting projections of the full three-body dynamics onto adiabatic Born--Oppenheimer approximation. Coupling of electron and nuclear quantum dynamics is clearly seen. Nuclear pair correlation function is found to broaden by 0.040 a_0 and average bond length is larger by 0.056 a_0. Also, non-adiabatic correction to the binding energy is found. Electronic distribution is affected less, and therefore, we could say that the adiabatic approximation is better for the electron than for the nuclei.
We study the potential energy surface of the ozone molecule by means of Quantum Monte Carlo simulations based on the resonating valence bond concept. The trial wave function consists of an antisymmetrized geminal power arranged in a single-determinant that is multiplied by a Jastrow correlation factor. Whereas the determinantal part incorporates static correlation effects, the augmented real-space correlation factor accounts for the dynamics electron correlation. The accuracy of this approach is demonstrated by computing the potential energy surface for the ozone molecule in three vibrational states: symmetric, asymmetric and scissoring. We find that the employed wave function provides a detailed description of rather strongly-correlated multi-reference systems, which is in quantitative agreement with experiment.
We explore the extended Koopmans theorem (EKT) within the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method. The EKT allows for the direct calculation of electron addition and removal spectral functions using reduced density matrices of the $N$-particle system, and avoids the need for analytic continuation. The lowest level of EKT with AFQMC, called EKT1-AFQMC, is benchmarked using small molecules, 14-electron and 54-electron uniform electron gas supercells, and diamond at the $Gamma$-point. Via comparison with numerically exact results (when possible) and coupled-cluster methods, we find that EKT1-AFQMC can reproduce the qualitative features of spectral functions for Koopmans-like charge excitations with errors in peak locations of less than 0.25 eV in a finite basis. We also note the numerical difficulties that arise in the EKT1-AFQMC eigenvalue problem, especially when back-propagated quantities are very noisy. We show how a systematic higher order EKT approach can correct errors in EKT1-based theories with respect to the satellite region of the spectral function. Our work will be of use for the study of low-energy charge excitations and spectral functions in correlated molecules and solids where AFQMC can be reliably performed.
We calculate the linear and non-linear susceptibilities of periodic longitudinal chains of hydrogen dimers with different bond-length alternations using a diffusion quantum Monte Carlo approach. These quantities are derived from the changes in electronic polarization as a function of applied finite electric field - an approach we recently introduced and made possible by the use of a Berry-phase, many-body electric-enthalpy functional. Calculated susceptibilities and hyper-susceptibilities are found to be in excellent agreement with the best estimates available from quantum chemistry - usually extrapolations to the infinite-chain limit of calculations for chains of finite length. It is found that while exchange effects dominate the proper description of the susceptibilities, second hyper-susceptibilities are greatly affected by electronic correlations. We also assess how different approximations to the nodal surface of the many-body wavefunction affect the accuracy of the calculated susceptibilities.
Monte Carlo simulations using a hybrid quantum and classical mechanical potential were performed for crystal and amorphous-like HCl-water(n) clusters The subsystem composed by HCl and one water molecule was treated within Density Functional Theory, and a classical force field was used for the rest of the system. Simulations performed at 200 K suggest that the energetic feasibility of HCl dissociation strongly depends on its initial placement within the cluster. An important degree of ionization occurs only if HCl is incorporated into the surface. We observe that local melting does not play a crucial role in the ionization process.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
