Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias

Details Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias

2011-11-07

arxiv.org/abs/1111.1568v2

Physics

Condensed Matter

Statistical Mechanics

27 pages, 8 figures

Scientific paper

We propose an efficient method to compute reaction rate constants of thermally activated processes occurring in many-body systems at finite temperature. The method consists in two steps: first, paths are sampled using a transition path sampling (TPS) algorithm supplemented with a local Lyapunov bias favoring diverging trajectories. This enhances the probability of observing rare reactive trajectories between stable states during relatively short simulations. Secondly, reaction constants are eventually estimated from the unbiased fraction of reactive paths, yielded by an appropriate statistical data analysis tool, the multistate Bennett acceptance ratio (MBAR) package. In order to test our algorithm, we compute reaction constants for structural transitions in LJ38, a well studied Lennard-Jones cluster, comparing our results to values previously reported in the literature. Additionally, we apply our method to the calculation of reaction rates for vacancy migration in an {\alpha}-Iron crystal, for temperatures ranging from 300 K to 850 K. Vacancy diffusion rates associated with activation barriers at finite temperature are then evaluated, and compared to values obtained using the standard harmonic approximation.

Affiliated with

Also associated with

No associations

Rating

Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-691340