Related Experiment Videos
Chemical reaction rates from ring polymer molecular dynamics
Ian R Craig1, David E Manolopoulos
1Physical and Theoretical Chemistry Laboratory, Oxford University, South Parks Road, Oxford OX1 3QZ, United Kingdom.
The Journal of Chemical Physics
|April 20, 2005
Summary
Ring-polymer molecular dynamics (RPMD) can now calculate reaction rate coefficients for condensed phase reactions. This new method offers comparable accuracy to existing techniques but is simpler to implement for complex chemical systems.
Area of Science:
- Chemical Physics
- Computational Chemistry
- Reaction Dynamics
Background:
- Calculating reaction rate coefficients in condensed phases is crucial for understanding chemical processes.
- Existing methods like classical Wigner and centroid molecular dynamics have limitations in accuracy and applicability.
- Accurate modeling of quantum tunneling and system-bath interactions is essential.
Purpose of the Study:
- To adapt the ring-polymer molecular dynamics (RPMD) method for calculating Kubo-transformed flux-side correlation functions.
- To determine rate coefficients for condensed phase reactions using the adapted RPMD method.
- To assess the accuracy and ease of application of the new RPMD approach compared to existing methods.
Main Methods:
- Adaptation of the ring-polymer molecular dynamics (RPMD) method.
- Calculation of approximate Kubo-transformed flux-side correlation functions.
- Application to a quartic double-well potential model with a harmonic oscillator bath.
- Development of a logarithmic discretization scheme for the Ohmic bath.
- Analysis of transmission through an Eckart barrier.
Main Results:
- The adapted RPMD method yields results comparable in accuracy to classical Wigner and centroid molecular dynamics.
- The method simplifies calculations by avoiding Wigner transforms and separate path integral computations.
- A novel discretization scheme for the Ohmic bath requires only nine modes for converged results.
- Satisfactory description of deep quantum tunneling regimes was achieved for Eckart barrier transmission.
- Exact quantum-mechanical rate constant for parabolic barrier transmission was analytically demonstrated.
Conclusions:
- The adapted RPMD method provides an accurate and more accessible approach for calculating condensed phase reaction rate coefficients.
- The method shows promise for broader applications in complex chemical systems.
- The efficient bath discretization scheme significantly reduces computational cost.
- The method accurately captures quantum mechanical tunneling effects, crucial for many reactions.