ﻻ يوجد ملخص باللغة العربية
We introduce a novel approach for a fully quantum description of coupled electron-ion systems from first principles. It combines the variational quantum Monte Carlo (QMC) solution of the electronic part with the path integral (PI) formalism for the quantum nuclear dynamics. On the one hand, the PI molecular dynamics includes nuclear quantum effects by adding a set of fictitious classical particles (beads) aimed at reproducing nuclear quantum fluctuations via a harmonic kinetic term. On the other hand, variational QMC can provide Born-Oppenheimer (BO) potential energy surfaces with a precision comparable to the most advanced post Hartree-Fock approaches, and with a favorable scaling with the system size. To deal with the intrinsic QMC noise, we generalize the PI molecular dynamics using a Langevin thermostat correlated according to the covariance matrix of QMC nuclear forces. The variational parameters of the QMC wave function are evolved during the nuclear dynamics, such that the BO potential energy surface is unbiased. Statistical errors on the wave function parameters are reduced by resorting to bead grouping average, which we show to be accurate and well controlled. Our general algorithm relies on a Trotter breakup between the dynamics driven by ionic forces and the one set by the harmonic interbead couplings. The latter is exactly integrated even in presence of the Langevin thermostat, thanks to the mapping onto an Ornstein-Uhlenbeck process. This framework turns out to be very efficient also in the case of deterministic ionic forces. The new implementation is validated on the Zundel ion by direct comparison with standard PI Langevin dynamics calculations made with a coupled cluster potential energy surface. Nuclear quantum effects are confirmed to be dominant over thermal effects well beyond room temperature giving the excess proton an increased mobility by quantum tunneling.
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 simu
The properties of molecules and materials containing light nuclei are affected by their quantum mechanical nature. Modelling these quantum nuclear effects accurately requires computationally demanding path integral techniques. Considerable success ha
We present an overview of the variational and diffusion quantum Monte Carlo methods as implemented in the CASINO program. We particularly focus on developments made in the last decade, describing state-of-the-art quantum Monte Carlo algorithms and so
We investigate the use of optimized correlation consistent gaussian basis sets for the study of insulating solids with auxiliary-field quantum Monte Carlo (AFQMC). The exponents of the basis set are optimized through the minimization of the second or
We investigate the continuum limit that the number of beads goes to infinity in the ring polymer representation of thermal averages. Studying the continuum limit of the trajectory sampling equation sheds light on possible preconditioning techniques f