Quantum diffusion in liquid para-hydrogen from ring-polymer molecular dynamics

A path integral centroid molecular dynamics method for Bose and Fermi statistics

We present the aspects of the path integral molecular dynamics (PIMD) method relative to (a) its theoretical fundamental principles, (b) its applicability to model quantum systems and (c) its implementation as a simulation tool. The PIMD method is based on the discretized path integral representation of quantum mechanics. In this representation, a quantum particle is isomorphic to a closed polymer chain. The problem of the indistinguishability of quantum particles is tackled with a non-local exchange potential. When the exact density matrix of the quantum particles is used, the exchange potential is exact. We use a high temperature approximation to the density matrix leading to an approximate form of the exchange potential. This quantum molecular dynamics method allows the simulation of collections of quantum particles at finite temperature. Our algorithm can be made to scale linearly with the number of quantum states on which the density matrix is projected. Therefore, it can be optimized to run efficiently on parallel computers. We apply the PIMD method to the electron plasma in 3-dimension. The kinetic and potential energies are calculated and compared with results for similar systems simulated with a variational Monte Carlo method. Both results show good agreements with each other at all the densities studied. The method is then use to model the thermodynamic behavior of a simple alkali metal. In these simulations, ions and valence electrons are treated as classical and quantum particles, respectively. The simple metal undergoes a phase transformation upon heating. Furthermore, to demonstrate the richness of behaviors that can be studied with the PIMD method, we also report on the metal to insulator transition in a hydrogenoid lattice. Finally, in previous studies of alkali metals, electrons interacted with ions via local pseudo-potentials, an extension of the method to modeling electrons in non-local pseudo-potentials is also presented with applications.

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange-correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. In order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

We have theoretically investigated the finite system size effect of hadronic matter on its transport coefficients such as shear viscosity, bulk viscosity, and electrical conductivity. We have used a hadron resonance gas model to calculate the thermodynamical quantities such as entropy density, speed of sound and the above transport coefficients. All these quantities are found to be sensitive to finite system size effects of hadronic matter. The effect of finite system size is found to be more when the system is at low temperatures and is reduced at high temperatures. Owing to the intimate linking between system size and centrality, we have presented the centrality dependence of transport coefficients. We have also explored to link of our results with the macroscopic picture of hydrodynamical evolution.

Dynamical quantum phase transitions reveal singularities in quench dynamics, characterized by the emergence of Loschmidt echo zeros at critical times, which usually exist only in the thermodynamical limit but are absent in finite size quantum systems. In this work, we propose a theoretical scheme to probe Loschmidt echo zeros and observe dynamical quantum phase transitions in finite size systems by applying a two-step quenching protocol, which offers an experimentally feasible approach to study Loschmidt echo zeros. Using the transverse Ising model as a testbed, we identify that the exact Loschmidt echo zeros can be always accessed by tuning the quench rate, when the quench is across the phase transition point. The associated rate function displays divergence at critical times, accompanying with the change of the dynamical topological order parameter. The critical times are influenced by the quench rate, system size, and momentum modes, embodying the interplay between finite-size effects and critical dynamics.

The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)(2) is about 30% weaker when quantum nuclear effects are included, while that in (H(2)O)(2) is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation.

We examine the short-time accuracy of a class of approximate quantum dynamical techniques that includes the centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD) methods. Both of these methods are based on the path integral molecular dynamics (PIMD) technique for calculating the exact static equilibrium properties of quantum mechanical systems. For Kubo-transformed real-time correlation functions involving operators that are linear functions of positions or momenta, the RPMD and (adiabatic) CMD approximations differ only in the choice of the artificial mass matrix of the system of ring polymer beads that is employed in PIMD. The obvious ansatz for a general method of this type is therefore to regard the elements of the PIMD (or Parrinello-Rahman) mass matrix as an adjustable set of parameters that can be chosen to improve the accuracy of the resulting approximation. We show here that this ansatz leads uniquely to the RPMD approximation when the criterion that is used to select the mass matrix is the short-time accuracy of the Kubo-transformed correlation function. In particular, we show that the leading error in the RPMD position autocorrelation function is O(t(8)) and the error in the velocity autocorrelation function is O(t(6)), for a general anharmonic potential. The corresponding errors in the CMD approximation are O(t(6)) and O(t(4)), respectively.

Effects of temperature and isotopic substitution on electron attachment dynamics of guanine–cytosine base pair: Ring-polymer and classical molecular dynamics simulations

In a previous joint experimental and theoretical study of the barrierless chemical reaction C(1D) + H2 at low temperatures (300-50 K) [K. M. Hickson, J.-C. Loison, H. Guo and Y. V. Suleimanov, J. Phys. Chem. Lett., 2015, 6, 4194], excellent agreement was found between experimental thermal rate constants and theoretical estimates based on ring polymer molecular dynamics (RPMD) over the two lowest singlet potential energy surfaces (PESs). Here, we extend this work to one of its isotopologues, C(1D) + D2, over the same temperature range. Experimental and RPMD results are in good agreement when contributions from both PESs to this chemical reaction are included in the RPMD simulations. The deviation between experiment and the RPMD calculations does not exceed 25% and both results exhibit a slight negative temperature dependence. The first excited 1A'' PES plays a more important role than the ground 1A' PES as the temperature is decreased, similar to our previous studies of the C(1D) + H2 reaction but with a more pronounced effect. The small differences in temperature dependence between the earlier and present experimental studies of C(1D) + H2/D2 reactions are discussed in terms of the use of non-equilibrium populations of ortho/para-H2/D2. We argue that RPMD provides a very convenient and reliable tool to study low-temperature chemical reactions.

To fast and accurately compute rate coefficients of the H/D + CH4 → H2/HD + CH3reactions, we propose a segmented strategy for fitting suitable potential energy surface (PES), on which ring-polymer molecular dynamics (RPMD) simulations are performed. On the basis of recently developed permutation invariant polynomial neural-network approach [J. Li et al., J. Chem. Phys. 142, 204302 (2015)], PESs in local configuration spaces are constructed. In this strategy, global PES is divided into three parts, including asymptotic, intermediate, and interaction parts, along the reaction coordinate. Since less fitting parameters are involved in the local PESs, the computational efficiency for operating the PES routine is largely enhanced by a factor of ∼20, comparing with that for global PES. On interaction part, the RPMD computational time for the transmission coefficient can be further efficiently reduced by cutting off the redundant part of the child trajectories. For H + CH4, good agreements among the present RPMD rates and those from previous simulations as well as experimental results are found. For D + CH4, on the other hand, qualitative agreement between present RPMD and experimental results is predicted.

ABSTRACTThermal quantum time-correlation functions are of fundamental importance in quantum dynamics, allowing experimentally measurable properties such as reaction rates, diffusion constants and vibrational spectra to be computed from first principles. Since the exact quantum solution scales exponentially with system size, there has been considerable effort in formulating reliable linear-scaling methods involving exact quantum statistics and approximate quantum dynamics modelled with classical-like trajectories. Here, we review recent progress in the field with the development of methods including centroid molecular dynamics , ring polymer molecular dynamics (RPMD) and thermostatted RPMD (TRPMD). We show how these methods have recently been obtained from ‘Matsubara dynamics’, a form of semiclassical dynamics which conserves the quantum Boltzmann distribution. We also apply the Matsubara formalism to reaction rate theory, rederiving t → 0+ quantum transition-state theory (QTST) and showing that Matsubara-TST, like RPMD-TST, is equivalent to QTST. We end by surveying areas for future progress.

i-PI: A Python interface for ab initio path integral molecular dynamics simulations

Equilibrium magnesium isotope fractionation between aqueous Mg2+ and carbonate minerals: Insights from path integral molecular dynamics

We have developed a unified quantum molecular dynamics method, in which the centroid molecular dynamics including the quantum effect of atomic nuclei is unified with the semiempirical molecular orbital calculation for the atomic configuration generated in the molecular dynamics. This unified method has successfully been applied to a couple of diatomic molecular systems for examining its performance. The method has yielded excellent agreement with exact quantum results and taken sufficiently short CPU time for practical use. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003

We present calculations of tunneling splittings in selected small water clusters, based on a recently developed path integral molecular dynamics (PIMD) method. The ground-rotational-state tunneling motions associated with the largest splittings in the water dimer, trimer, and hexamer are considered, and we show that the PIMD predictions are in very good agreement with benchmark quantum and experimental results. As the tunneling spectra are highly sensitive to both the details of the quantum dynamics and the potential energy surface, our calculations are a validation of the MB-Pol surface as well as the accuracy of PIMD. The favorable scaling of PIMD with system size paves the way for calculations of tunneling splittings in large, nonrigid molecular systems with motions that cannot be treated accurately by other methods, such as the semiclassical instanton.