Estimating time-correlation functions by sampling and unbiasing dynamically activated events

Abstract

Transition path sampling is a rare-event method that estimates state-to-state timecorrelation functions in many-body systems from samples of short trajectories. In this framework, it is proposed to bias the importance function using the lowest Jacobian eigenvalue moduli along the dynamical trajectory. A lowest eigenvalue modulus is related to the lowest eigenvalue of the Hessian matrix and is evaluated here using the Lanczos algorithm as in activation-relaxation techniques. This results in favoring the sampling of activated trajectories and enhancing the occurrence of the rare reactive trajectories of interest, those corresponding to transitions between locally stable states. Estimating the time-correlation functions involves unbiasing the sample of simulated trajectories which is done using the multi-state Bennett acceptance ratio (MBAR) method. To assess the performance of our procedure, we compute the time-correlation function associated with the migration of a vacancy in α-iron. The derivative of the estimated time-correlation function yields a migration rate in agreement with the one given by transition state theory. Besides, we show that the information relative to rejected trajectories can be recycled within MBAR, resulting in a substantial speed-up. Unlike original transition path-sampling, our approach does not require computing the reversible work to confine the trajectory endpoints to a reactive state.

34 pages, 8 figures

References in corpus (9)

• Rational Construction of Stochastic Numerical Methods for Molecular Sampling
• Probing rare physical trajectories with Lyapunov weighted dynamics
• Some improvements of the ART method for finding transition pathways on potential energy surfaces
• Computing stationary distributions in equilibrium and non-equilibrium systems with Forward Flux Sampling
• Finite-temperature critical point of a glass transition
• Simulating structural transitions by direct transition current sampling: the example of LJ38
• Kramers equation and supersymmetry
• Path sampling with Stochastic Dynamics: Some new Algorithms
• Waste-recycling Monte Carlo with optimal estimates: application to free energy calculations in alloys

Cited by in corpus (4)

• Solute diffusion by self-interstitial defects and radiation-induced segregation in ferritic Fe-X (X=Cr, Cu, Mn, Ni, P, Si) dilute alloys
• Free energy calculations from adaptive molecular dynamics simulations with adiabatic reweighting
• Estimating thermodynamic expectations and free energies in expanded ensemble simulations: systematic variance reduction through conditioning
• Using Bayes formula to estimate rates of rare events in transition path sampling simulations