Direct evaluation of the temperature dependence of the rate constant based on the quantum instanton approximation
Abstract
A general method for the direct evaluation of the temperature dependence of the quantum-mechanical reaction rate constant in many-dimensional systems is described. The method is based on the quantum instanton approximation for the rate constant, thermodynamic integration with respect to the inverse temperature, and the path integral Monte Carlo evaluation. It can describe deviations from the Arrhenius law due to the coupling of rotations and vibrations, zero-point energy, tunneling, corner-cutting, and other nuclear quantum effects. The method is tested on the Eckart barrier and the full-dimensional H+H2→H2+H reaction. In the temperature range from 300 to 1500 K, the error of the present method remains within 13% despite the very large deviations from the Arrhenius law. The direct approach makes the calculations much more efficient, and the efficiency is increased even further (by up to two orders of magnitude in the studied reactions) by using optimal estimators for reactant and transition state thermal energies. Which of the estimators is optimal, however, depends on the system and the strength of constraint in a constrained simulation.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- May 2010
- DOI:
- arXiv:
- arXiv:1004.0201
- Bibcode:
- 2010JChPh.132s4106B
- Keywords:
-
- atom-molecule reactions;
- hydrogen neutral atoms;
- hydrogen neutral molecules;
- Monte Carlo methods;
- quantum chemistry;
- reaction rate constants;
- thermodynamic properties;
- 82.20.Pm;
- 82.20.Xr;
- 82.20.Db;
- 82.60.-s;
- 82.30.Cf;
- Rate constants reaction cross sections and activation energies;
- Quantum effects in rate constants;
- Transition state theory and statistical theories of rate constants;
- Chemical thermodynamics;
- Atom and radical reactions;
- chain reactions;
- molecule-molecule reactions;
- Physics - Chemical Physics;
- Physics - Computational Physics;
- Quantum Physics
- E-Print:
- 11 pages, 8 figures, submitted to J. Chem. Phys.