Ring-Polymer Molecular Dynamics and Kinetics for the H + CH → H + CH Reaction Using the Full-Dimensional Potential Energy Surface.

The H + CH → H + CH reaction is important in understanding the production mechanisms of anionic molecules in interstellar environments. Herein, the rate coefficients for the H + CH → H + CH reaction were calculated using ring-polymer molecular dynamics (RPMD), classical molecular dynamics (MD), and quasi-classical trajectory (QCT) approaches on a newly developed ab initio potential energy surface (PES) in full dimensions. PES was constructed by fitting a large number of ab initio energy points and their gradients using the permutationally invariant polynomial basis set method. There was no barrier in the reaction coordinates, which was a collinear-dominated reaction, and the reaction proceeded exothermically. It is found that the fitted PES provides the appropriate thermal rate coefficients based on all RPMD, classical MD, and QCT simulations at higher temperatures. The evaluation of the rate coefficients at lower temperatures should be conducted carefully because the fitting of the PES associated with the long-range interaction should be further improved. The spatial distribution of the nucleus allows a more effective attraction between the reactants.