Correcting Rate Constants from Anharmonic Molecular Dynamics for Quantum Effects.

Area of Science:

• Computational Chemistry
• Chemical Kinetics
• Reaction Dynamics

Background:

• Standard rigid-rotor harmonic-oscillator models often yield inaccurate rate constants due to neglecting anharmonicity.
• Reactive molecular dynamics (MD) simulations offer a path to anharmonic calculations but are limited by classical approximations and force field accuracy.
• Accurate calculation of chemical reaction rates is crucial for understanding and predicting chemical processes.

Purpose of the Study:

• To develop and validate a hybrid computational approach for accurate calculation of anharmonic rate constants.
• To combine the strengths of classical molecular dynamics (MD) and quantum-mechanical (QM) electronic structure calculations.
• To assess the accuracy and applicability of the new method for hydrogen abstraction reactions.

Main Methods:

• Developed an approximation pairing anharmonic information from classical MD with high-accuracy QM energies and frequencies.
• Applied the scheme to hydrogen abstraction reactions in methane systems for benchmarking.
• Investigated specific reaction pathways, including CH3• + H• → CH22• + H2, to identify failure modes.

Main Results:

• The corrected rate constants showed a standard deviation of 2.6 compared to experimental values for methane hydrogen abstractions.
• Analysis identified failure archetypes for the CH3• + H• reaction, defining the method's application range.
• Within its application range, the method achieved a standard deviation of 2.1, demonstrating improved accuracy.

Conclusions:

• The proposed hybrid method effectively incorporates anharmonicity for more accurate rate constant calculations.
• The computational efficiency and scaling allow for application to larger, more complex systems, such as hydrogen abstraction from 2-butanone by HO2•.
• This approach bridges the gap between classical dynamics efficiency and QM accuracy for chemical kinetics.