skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Should thermostatted ring polymer molecular dynamics be used to calculate thermal reaction rates?

Journal Article · · Journal of Chemical Physics
DOI:https://doi.org/10.1063/1.4928599· OSTI ID:22493532
 [1];  [2]
  1. Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW (United Kingdom)
  2. Computation-based Science and Technology Research Center, Cyprus Institute, 20 Kavafi St., Nicosia 2121 (Cyprus)

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 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{sub 2}, D + MuH, and F + H{sub 2}, and the prototypical polyatomic reaction H + CH{sub 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 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.

OSTI ID:
22493532
Journal Information:
Journal of Chemical Physics, Vol. 143, Issue 7; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
Country of Publication:
United States
Language:
English