Do you want to publish a course? Click here

Hybrid Quantum and Classical Mechanical Monte Carlo Simulations of the Interaction of Hydrogen Chloride with Solid Water Clusters

57   0   0.0 ( 0 )
 Added by Kohanoff Jorge Jose
 Publication date 1997
  fields Physics
and research's language is English




Ask ChatGPT about the research

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.



rate research

Read More

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.
184 - S. Azadi , W. M. C. Foulkes , 2013
We use the diffusion quantum Monte Carlo (DMC) method to calculate the ground state phase diagram of solid molecular hydrogen and examine the stability of the most important insulating phases relative to metallic crystalline molecular hydrogen. We develop a new method to account for finite-size errors by combining the use of twist-averaged boundary conditions with corrections obtained using the Kwee-Zhang-Krakauer (KZK) functional in density functional theory. To study band-gap closure and find the metallization pressure, we perform accurate quasi-particle many-body calculations using the $GW$ method. In the static approximation, our DMC simulations indicate a transition from the insulating Cmca-12 structure to the metallic Cmca structure at around 375 GPa. The $GW$ band gap of Cmca-12 closes at roughly the same pressure. In the dynamic DMC phase diagram, which includes the effects of zero-point energy, the Cmca-12 structure remains stable up to 430 GPa, well above the pressure at which the $GW$ band gap closes. Our results predict that the semimetallic state observed experimentally at around 360 GPa [Phys. Rev. Lett. {bf 108}, 146402 (2012)] may correspond to the Cmca-12 structure near the pressure at which the band gap closes. The dynamic DMC phase diagram indicates that the hexagonal close packed $P6_3/m$ structure, which has the largest band gap of the insulating structures considered, is stable up to 220 GPa. This is consistent with recent X-ray data taken at pressures up to 183 GPa [Phys. Rev. B {bf 82}, 060101(R) (2010)], which also reported a hexagonal close packed arrangement of hydrogen molecules.
We propose the use of preconditioning in FCIQMC which, in combination with perturbative estimators, greatly increases the efficiency of the algorithm. The use of preconditioning allows a time step close to unity to be used (without time-step errors), provided that multiple spawning attempts are made per walker. We show that this approach substantially reduces statistical noise on perturbative corrections to initiator error, which improve the accuracy of FCIQMC but which can suffer from significant noise in the original scheme. Therefore, the use of preconditioning and perturbatively-corrected estimators in combination leads to a significantly more efficient algorithm. In addition, a simpler approach to sampling variational and perturbative estimators in FCIQMC is presented, which also allows the variance of the energy to be calculated. These developments are investigated and applied to benzene (30e,108o), an example where accurate treatment is not possible with the original method.
We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely-used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.
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.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا