Extracting the diffusion tensor from molecular dynamics simulation with Milestoning
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712 (United States)
We propose an algorithm to extract the diffusion tensor from Molecular Dynamics simulations with Milestoning. A Kramers-Moyal expansion of a discrete master equation, which is the Markovian limit of the Milestoning theory, determines the diffusion tensor. To test the algorithm, we analyze overdamped Langevin trajectories and recover a multidimensional Fokker-Planck equation. The recovery process determines the flux through a mesh and estimates local kinetic parameters. Rate coefficients are converted to the derivatives of the potential of mean force and to coordinate dependent diffusion tensor. We illustrate the computation on simple models and on an atomically detailed system—the diffusion along the backbone torsions of a solvated alanine dipeptide.
- OSTI ID:
- 22415451
- Journal Information:
- Journal of Chemical Physics, Vol. 142, Issue 1; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
The probability equation for the cosmological comoving curvature perturbation
Recovering position-dependent diffusion from biased molecular dynamics simulations