Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHow to obtain thermostatted ring polymer molecular dynamics from exact quantum dynamics and when to use it
90
0
0.0
(
0
)
Added by
Timothy Hele MChem PhD
Publication date
2015
fields
Physics
and research's language is
English
Authors
Timothy J. H. Hele
Ask ChatGPT about the research
No Arabic abstract
We obtain thermostatted ring polymer molecular dynamics (TRPMD) from exact quantum dynamics via Matsubara dynamics, a recently-derived form of linearization which conserves the quantum Boltzmann distribution. Performing a contour integral in the complex quantum Boltzmann distribution of Matsubara dynamics, replacement of the imaginary Liouvillian which results with a Fokker-Planck term gives TRPMD. We thereby provide error terms between TRPMD and quantum dynamics and predict the systems in which they are likely to be small. Using a harmonic analysis we show that careful addition of friction causes the correct oscillation frequency of the higher ring-polymer normal modes in a harmonic well, which we illustrate with calculation of the position-squared autocorrelation function. However, no physical friction parameter will produce the correct fluctuation dynamics for a parabolic barrier. The results in this paper are consistent with previous numerical studies and advise the use of TRPMD for the computation of spectra.
rate research
Read More
We apply Thermostatted Ring Polymer Molecular Dynamics (TRPMD), a recently-proposed approximate quantum dynamics method, to the computation of thermal reaction rates. Its short-time Transition-State Theory (TST) limit is identical to rigorous Quantum Transition-State Theory, and we find that its long-time limit is independent of the location of the dividing surface. TRPMD rate theory is then applied to one-dimensional model systems, the atom-diatom bimolecular reactions H+H$_2$, D+MuH and F+H$_2$, and the prototypical polyatomic reaction H+CH$_4$. Above the crossover temperature, the TRPMD rate is virtually invariant to the strength of the friction applied to the internal ring-polymer normal modes, and beneath the crossover temperature the TRPMD rate generally decreases with increasing friction, in agreement with the predictions of Kramers theory. We therefore find that TRPMD is approximately equal to, or less accurate than, Ring Polymer Molecular Dynamics (RPMD) for symmetric reactions, and for certain asymmetric systems and friction parameters closer to the quantum result, providing a basis for further assessment of the accuracy of this method.
Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostat to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the desirable features of RPMD: thermostatted RPMD (TRPMD) is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.
We present a new non-adiabatic ring polymer molecular dynamics (NRPMD) method based on the spin mapping formalism, which we refer to as the spin-mapping NRPMD (SM-NRPMD) approach. We derive the path-integral partition function expression using the spin coherent state basis for the electronic states and the ring polymer formalism for the nuclear degrees of freedom (DOFs). This partition function provides an efficient sampling of the quantum statistics. Using the basic property of the Stratonovich-Weyl transformation, we derive a Hamiltonian which we propose for the dynamical propagation of the coupled spin mapping variables and the nuclear ring polymer. The accuracy of the SM-NRPMD method is numerically demonstrated by computing nuclear position and population auto-correlation functions of non-adiabatic model systems. The results from SM-NRPMD agree very well with the numerically exact results. The main advantage of using the spin mapping variables over the harmonic oscillator mapping variables is numerically demonstrated, where the former provides nearly time-independent expectation values of physical observables for systems under thermal equilibrium, the latter can not preserve the initial quantum Boltzmann distribution. We also explicitly demonstrate that SM-NRPMD provides invariant dynamics upon various ways of partitioning the state-dependent and state-independent potentials.
In a previous article [J. Chem. Phys. 138, 084108 (2013)], we showed that the $tto 0_+$ limit of ring-polymer molecular dynamics (RPMD) rate-theory is also the $tto 0_+$ limit of a new type of quantum flux-side time-correlation function, in which the dividing surfaces are invariant to imaginary-time translation; in other words, that RPMD transition-state theory (RPMD-TST) is a $tto 0_+$ quantum transition-state theory (QTST). Recently, Jang and Voth [J. Chem. Phys. 144, 084110 (2016)] rederived this quantum $tto 0_+$ limit, and claimed that it gives instead the centroid-density approximation. Here we show that the $tto 0_+$ limit derived by Jang and Voth is in fact RPMD-TST.
Path integral-based simulation methodologies play a crucial role for the investigation of nuclear quantum effects by means of computer simulations. However, these techniques are significantly more demanding than corresponding classical simulations. To reduce this numerical effort, we recently proposed a method, based on a rigorous Hamiltonian formulation, which restricts the quantum modeling to a small but relevant spatial region within a larger reservoir where particles are treated classically. In this work, we extend this idea and show how it can be implemented along with state-of-the-art path integral simulation techniques, such as ring polymer and centroid molecular dynamics, which allow the approximate calculation of both quantum statistical and quantum dynamical properties. To this end, we derive a new integration algorithm which also makes use of multiple time-stepping. The scheme is validated via adaptive classical--path-integral simulations of liquid water. Potential applications of the proposed multiresolution method are diverse and include efficient quantum simulations of interfaces as well as complex biomolecular systems such as membranes and proteins.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
