Exact milestoning
Abstract
A new theory and an exact computer algorithm for calculating kinetics and thermodynamic properties of a particle system are described. The algorithm avoids trapping in metastable states, which are typical challenges for Molecular Dynamics (MD) simulations on rough energy landscapes. It is based on the division of the full space into Voronoi cells. Prior knowledge or coarse sampling of space points provides the centers of the Voronoi cells. Short time trajectories are computed between the boundaries of the cells that we call milestones and are used to determine fluxes at the milestones. The flux function, an essential component of the new theory, provides a complete description of the statistical mechanics of the system at the resolution of the milestones. We illustrate the accuracy and efficiency of the exact Milestoning approach by comparing numerical results obtained on a model system using exact Milestoning with the results of long trajectories and with a solution of the corresponding FokkerPlanck equation. The theory uses an equation that resembles the approximate Milestoning method that was introduced in 2004 [A. K. Faradjian and R. Elber, J. Chem. Phys. 120(23), 1088010889 (2004)]. However, the current formulation is exact and is still significantly more efficient than straightforward MDmore »
 Authors:
 Department of Chemistry, Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, Texas 78712 (United States)
 (United States)
 Publication Date:
 OSTI Identifier:
 22416200
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Chemical Physics; Journal Volume: 142; Journal Issue: 9; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCURACY; ALGORITHMS; APPROXIMATIONS; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; FOKKERPLANCK EQUATION; MATHEMATICAL SOLUTIONS; METASTABLE STATES; MOLECULAR DYNAMICS METHOD; STATISTICAL MECHANICS; THERMODYNAMIC PROPERTIES; TRAJECTORIES; TRAPPING
Citation Formats
BelloRivas, Juan M., Elber, Ron, and Department of Chemistry, University of Texas at Austin, Austin, Texas 78712. Exact milestoning. United States: N. p., 2015.
Web. doi:10.1063/1.4913399.
BelloRivas, Juan M., Elber, Ron, & Department of Chemistry, University of Texas at Austin, Austin, Texas 78712. Exact milestoning. United States. doi:10.1063/1.4913399.
BelloRivas, Juan M., Elber, Ron, and Department of Chemistry, University of Texas at Austin, Austin, Texas 78712. 2015.
"Exact milestoning". United States.
doi:10.1063/1.4913399.
@article{osti_22416200,
title = {Exact milestoning},
author = {BelloRivas, Juan M. and Elber, Ron and Department of Chemistry, University of Texas at Austin, Austin, Texas 78712},
abstractNote = {A new theory and an exact computer algorithm for calculating kinetics and thermodynamic properties of a particle system are described. The algorithm avoids trapping in metastable states, which are typical challenges for Molecular Dynamics (MD) simulations on rough energy landscapes. It is based on the division of the full space into Voronoi cells. Prior knowledge or coarse sampling of space points provides the centers of the Voronoi cells. Short time trajectories are computed between the boundaries of the cells that we call milestones and are used to determine fluxes at the milestones. The flux function, an essential component of the new theory, provides a complete description of the statistical mechanics of the system at the resolution of the milestones. We illustrate the accuracy and efficiency of the exact Milestoning approach by comparing numerical results obtained on a model system using exact Milestoning with the results of long trajectories and with a solution of the corresponding FokkerPlanck equation. The theory uses an equation that resembles the approximate Milestoning method that was introduced in 2004 [A. K. Faradjian and R. Elber, J. Chem. Phys. 120(23), 1088010889 (2004)]. However, the current formulation is exact and is still significantly more efficient than straightforward MD simulations on the system studied.},
doi = {10.1063/1.4913399},
journal = {Journal of Chemical Physics},
number = 9,
volume = 142,
place = {United States},
year = 2015,
month = 3
}

We propose an algorithm to extract the diffusion tensor from Molecular Dynamics simulations with Milestoning. A KramersMoyal 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 FokkerPlanck 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 themore »

Modeling kinetics and equilibrium of membranes with fields: Milestoning analysis and implication to permeation
Coarse graining of membrane simulations by translating atomistic dynamics to densities and fields with Milestoning is discussed. The space of the membrane system is divided into cells and the different cells are characterized by order parameters presenting the number densities. The dynamics of the order parameters are probed with Milestoning. The methodology is illustrated here for a phospholipid membrane system (a hydrated bilayer of DOPC (1,2dioleoylsnglycero3phosphocholine) lipid molecules). Significant inhomogeneity in membrane internal number density leads to complex free energy landscape and local maps of transition times. Dynamics and distributions of cavities within the membrane assist the permeation of nonpolarmore » 
Exact and almost exact solutions to the VlasovMaxwell system
Exact and almost exact solutions to the VlasovMaxwell systems describing a variety of plasma configurations with density, temperature, and current gradients are presented. Possible consequences of these solutions are also discussed. 
Exact and almost exact VlasovMaxwell equilibria
The author has attempted here to obtain a class of confined exact nonlinear solutions to the MaxwellVlasov system in cylindrical geometry. The problem becomes more and more difficult as a larger and larger set of plasma parameters (density = n, drift speed = u, temperature = T) is allowed to have radial dependence. The solution to the problem with only density gradients was presented by Harris, and constitutes a kinetic derivation of Bennett profiles for a cylindrical pinch.