TY - JOUR

T1 - Comparisons of classical chemical dynamics simulations of the unimolecular decomposition of classical and quantum microcanonical ensembles

AU - Manikandan, Paranjothy

AU - Hase, William L.

N1 - Funding Information:
This material is based upon work supported by the National Science Foundation under Grant No. CHE-0957521 and the Robert A. Welch Foundation under Grant No. D-0005. Support was also provided by the High Performance Computing Center (HPCC) at Texas Tech University, under the direction of Philip W. Smith.

PY - 2012/5/14

Y1 - 2012/5/14

N2 - Previous studies have shown that classical trajectory simulations often give accurate results for short-time intramolecular and unimolecular dynamics, particularly for initial non-random energy distributions. To obtain such agreement between experiment and simulation, the appropriate distributions must be sampled to choose initial coordinates and momenta for the ensemble of trajectories. If a molecules classical phase space is sampled randomly, its initial decomposition will give the classical anharmonic microcanonical (RRKM) unimolecular rate constant for its decomposition. For the work presented here, classical trajectory simulations of the unimolecular decomposition of quantum and classical microcanonical ensembles, at the same fixed total energy, are compared. In contrast to the classical microcanonical ensemble, the quantum microcanonical ensemble does not sample the phase space randomly. The simulations were performed for CH 4, C 2H 5, and Cl - - -CH 3Br using both analytic potential energy surfaces and direct dynamics methods. Previous studies identified intrinsic RRKM dynamics for CH 4 and C 2H 5, but intrinsic non-RRKM dynamics for Cl - - -CH 3Br. Rate constants calculated from trajectories obtained by the time propagation of the classical and quantum microcanonical ensembles are compared with the corresponding harmonic RRKM estimates to obtain anharmonic corrections to the RRKM rate constants. The relevance and accuracy of the classical trajectory simulation of the quantum microcanonical ensemble, for obtaining the quantum anharmonic RRKM rate constant, is discussed.

AB - Previous studies have shown that classical trajectory simulations often give accurate results for short-time intramolecular and unimolecular dynamics, particularly for initial non-random energy distributions. To obtain such agreement between experiment and simulation, the appropriate distributions must be sampled to choose initial coordinates and momenta for the ensemble of trajectories. If a molecules classical phase space is sampled randomly, its initial decomposition will give the classical anharmonic microcanonical (RRKM) unimolecular rate constant for its decomposition. For the work presented here, classical trajectory simulations of the unimolecular decomposition of quantum and classical microcanonical ensembles, at the same fixed total energy, are compared. In contrast to the classical microcanonical ensemble, the quantum microcanonical ensemble does not sample the phase space randomly. The simulations were performed for CH 4, C 2H 5, and Cl - - -CH 3Br using both analytic potential energy surfaces and direct dynamics methods. Previous studies identified intrinsic RRKM dynamics for CH 4 and C 2H 5, but intrinsic non-RRKM dynamics for Cl - - -CH 3Br. Rate constants calculated from trajectories obtained by the time propagation of the classical and quantum microcanonical ensembles are compared with the corresponding harmonic RRKM estimates to obtain anharmonic corrections to the RRKM rate constants. The relevance and accuracy of the classical trajectory simulation of the quantum microcanonical ensemble, for obtaining the quantum anharmonic RRKM rate constant, is discussed.

UR - http://www.scopus.com/inward/record.url?scp=84862903092&partnerID=8YFLogxK

U2 - 10.1063/1.4714219

DO - 10.1063/1.4714219

M3 - Article

C2 - 22583280

AN - SCOPUS:84862903092

VL - 136

JO - The Journal of Chemical Physics

JF - The Journal of Chemical Physics

SN - 0021-9606

IS - 18

M1 - 184110

ER -