Assessing Parameters for Ring Polymer Molecular Dynamics Simulations at Low Temperatures: DH+H Chemical Reaction

Abstract

Ring polymer molecular dynamics (RPMD) is an accurate method for calculating thermal chemical reaction rates. It has recently been discovered that low-temperature calculations are strongly affected by the simulation parameters. Here, for the thermally activated reaction DH + H -> D + H2, we calculate the RPMD rate constants at T = 50, 100, and 300 K and demonstrate that for T >= 100 K the standard input parameters yield accurate results, but at low temperatures (e.g., 50 K) one must increase the asymptotic distance and force constant, and decrease the umbrella integration step.

References in corpus (7)

• Nuclear quantum effects in solids using a colored-noise thermostat
• Colored-noise thermostats à la carte
• Langevin equation with colored noise for constant-temperature molecular dynamics simulations
• The Impact of HD Cooling on the Formation of the First Stars
• An experimental and theoretical investigation of the C(1D) + D2 reaction
• Thermal Rate Coefficients for the Astrochemical Process C + CH${}^{+}$ $\to$ C${}_2^{+}$ + H by Ring Polymer Molecular Dynamics
• The low temperature D${}^{+}$ + H${}_2$ $\to$ HD + H${}^{+}$ reaction rate coefficient: a ring polymer molecular dynamics and quasi-classical trajectory study