entropy-logo

Journal Browser

Journal Browser

Understanding Molecular Dynamics via Stochastic Processes

A special issue of Entropy (ISSN 1099-4300).

Deadline for manuscript submissions: closed (30 September 2017) | Viewed by 79350

Special Issue Editors


E-Mail Website
Guest Editor
Department of Physics, University of Roma “La Sapienza”, Piazzale A. Moro 2, 00185 Roma, Italy
Interests: methods in molecular dynamics simulation of systems of statistical mechanical interest; equilibrium and non equilibrium molecular dynamics; rare events; computer simulation of complex molecular systems
Special Issues, Collections and Topics in MDPI journals

E-Mail Website
Guest Editor
Department of Physics, Computer Science and Mathematics of The Former Physics Department, University of Modena and Reggio Emilia, Via G. Campi 213/A, 41100 Modena, Italy
Interests: molecular dynamics simulation of condensed matter systems; solvation in h-bonded liquids; friction at the nanoscale; non equilibrium molecular dynamics; free energy computation in solution phase; simulation of tribological phenomena; transport properties in molecular fluids
Special Issues, Collections and Topics in MDPI journals

E-Mail Website
Guest Editor
Freie Universität Berlin, ZIB, Takustraße 7, 14195 Berlin, Germany
Interests: rare events statistics in molecular dynamics; coarse graining in molecular dynamics; conformation dynamics
Special Issues, Collections and Topics in MDPI journals

Special Issue Information

Dear Colleagues,

Contrary to what the name seems to suggest, Molecular Dynamics (MD) is not only about generating a dynamical trajectory of a system of molecules, but also, and foremost, about understanding the statistical properties of this trajectory. As a result, to make sense of this solution is to extract its main statistical features, which must be done within the probabilistic framework of Statistical Mechanics. Pushing this viewpoint one step further, an MD trajectory can be thought of not so much as a small piece of the actual trajectory of a realistic system of molecules, but rather as a sampling device for the statistical mechanics properties of this system.

When speaking about equilibrium quantities like free energy, etc., the probabilistic interpretation of an MD trajectory as a sampling device has long be recognized, and so has its computational advantages. One is that the numerical accuracy of an MD trajectory should be analyzed in terms of the statistics it produces rather than by pathwise comparison with an actual solution, which simplifies matters. A second advantage is that the sampling performances can be improved by fiddling with MD dynamics as long as this does not affect the statistical properties one is aiming at. This is why, for example, Monte Carlo is used as a perfectly valid alternative to MD to compute equilibrium quantities: from the present perspective, the two approaches are not so different in spirit.

Contrary to a widespread opinion in the community of practitioners, there is no reason to restrict this type of approach to the computation of the standard equilibrium quantities of Statistical Mechanics: the same philosophy in which MD is fully integrated within a probabilistic perspective can be applied to understand dynamical properties such as correlation functions, transport coefficients, pathways and rates of rare events, etc. In this context, also, identifying the right statistical quantities first, then using MD or whichever modification thereof to sample them can prove valuable. However, developing the right probabilistic framework for the study of dynamical phenomena, such as rare events, is a formidable challenge. This part of Statistical Mechanics is still much less developed. It also requires more sophisticated tools from Stochastic Processes Theory, for dynamical properties are multiple- rather than single-time statistical properties of the systems, i.e., one must deal with a stochastic process rather than with random variables. As a result the probability distributions relevant to dynamical phenomena are more complicated objects, often not even readily available. However, establishing what these distributions are and how to use MD as a tool to sample them efficiently is the right way to go. The steady growth in computing power as well as the development of various computational tricks may permit to generate ever longer trajectories in ever bigger systems; however bare trajectories have very little use without the right probabilistic framework to under- stand their meaning and looking at them in their gory details may even be more confusing than helpful.

It is in this spirit that, in this Special Issue, we would like to collect papers focusing, with a pedagogical aim, on the importance of stochastic process modeling, to understand and put on solid basis classical statistical mechanics and MD simulations, showing that not brute force MD but intelligent use of probability will give us the possibility to solve the most challenging problems when modeling physical, chemical and biological processes on all space and time scales. To comply with this goal we have invited a group of highly qualified researchers working along these lines to produce some useful glimpses of the ongoing process.

Planned Papers*

  • Ron Elber, Juan M. Bello-Rivas, Piao Ma, Alfredo Cardenas and Arman Fathizadeh, Calculating isocommittor surfaces as optimal reaction coordinates with milestoning
  • Anastasia S. Georgiou, Juan M. Bello-Rivas, C. William Gear, Hau-Tieng Wu, Eliodoro Chiavazzo and Ioannis G. Kevrekidis, An exploration algorithm for stochastic simulators driven by energy gradients
  • Wei Zhang and Christof Schütte, Reliable approximation of long relaxation timescales in molecular dynamics.
  • Josh Fass, David A. Sivak, Gavin E. Crooks and John D. Chodera, Which integrator is best for biomolecular simulation? Comparing the efficiencies of Langevin integrators
  • Carsten Hartmann, Lorenz Richter, Christof Schütte and Wei Zhang, Variational characterisation of free energy: theory and algorithms
  • Mattias Sachs, Vincent Danos and Ben Leimkuhler, Variable-coefficient stochastic particle models and their stationary states.
  • Robert L Jack, Marcus Kaiser and Johannes Zimmer, Symmetries and geometrical properties of dynamical fluctuations in molecular dynamics"
  • Robert Skeel and Youhan Fang, Comparing Markov Chain Samplers for Molecular Simulation
  • Chloe Gao and David T. Limmer, Transport coefficients from path ensemble free energies
  • Carsten Hartmann, Susana Gomes and Grigorios A. Pavliotis, On the linear response of a system of infinitely many randomly perturbed oscillators

[*The above list represents only planned manuscripts. Some of these manuscripts have not been received by the Editorial Office yet. Papers submitted to MDPI journals are subject to peer-review]

Prof. Dr. Giovanni Ciccotti
Prof. Dr. Mauro Ferrario
Prof. Dr. Christof Schuette
Guest Editors

Manuscript Submission Information

Manuscripts should be submitted online at www.mdpi.com by registering and logging in to this website. Once you are registered, click here to go to the submission form. Manuscripts can be submitted until the deadline. All submissions that pass pre-check are peer-reviewed. Accepted papers will be published continuously in the journal (as soon as accepted) and will be listed together on the special issue website. Research articles, review articles as well as short communications are invited. For planned papers, a title and short abstract (about 100 words) can be sent to the Editorial Office for announcement on this website.

Submitted manuscripts should not have been published previously, nor be under consideration for publication elsewhere (except conference proceedings papers). All manuscripts are thoroughly refereed through a single-blind peer-review process. A guide for authors and other relevant information for submission of manuscripts is available on the Instructions for Authors page. Entropy is an international peer-reviewed open access monthly journal published by MDPI.

Please visit the Instructions for Authors page before submitting a manuscript. The Article Processing Charge (APC) for publication in this open access journal is 2600 CHF (Swiss Francs). Submitted papers should be well formatted and use good English. Authors may use MDPI's English editing service prior to publication or during author revisions.

Benefits of Publishing in a Special Issue

  • Ease of navigation: Grouping papers by topic helps scholars navigate broad scope journals more efficiently.
  • Greater discoverability: Special Issues support the reach and impact of scientific research. Articles in Special Issues are more discoverable and cited more frequently.
  • Expansion of research network: Special Issues facilitate connections among authors, fostering scientific collaborations.
  • External promotion: Articles in Special Issues are often promoted through the journal's social media, increasing their visibility.
  • e-Book format: Special Issues with more than 10 articles can be published as dedicated e-books, ensuring wide and rapid dissemination.

Further information on MDPI's Special Issue polices can be found here.

Published Papers (12 papers)

Order results
Result details
Select all
Export citation of selected articles as:

Editorial

Jump to: Research

3 pages, 188 KiB  
Editorial
Molecular Dynamics vs. Stochastic Processes: Are We Heading Anywhere?
by Giovanni Ciccotti, Mauro Ferrario and Christof Schütte
Entropy 2018, 20(5), 348; https://doi.org/10.3390/e20050348 - 7 May 2018
Cited by 5 | Viewed by 4069
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)

Research

Jump to: Editorial

28 pages, 2989 KiB  
Article
Quantifying Configuration-Sampling Error in Langevin Simulations of Complex Molecular Systems
by Josh Fass, David A. Sivak, Gavin E. Crooks, Kyle A. Beauchamp, Benedict Leimkuhler and John D. Chodera
Entropy 2018, 20(5), 318; https://doi.org/10.3390/e20050318 - 26 Apr 2018
Cited by 26 | Viewed by 7415
Abstract
While Langevin integrators are popular in the study of equilibrium properties of complex systems, it is challenging to estimate the timestep-induced discretization error: the degree to which the sampled phase-space or configuration-space probability density departs from the desired target density due to the [...] Read more.
While Langevin integrators are popular in the study of equilibrium properties of complex systems, it is challenging to estimate the timestep-induced discretization error: the degree to which the sampled phase-space or configuration-space probability density departs from the desired target density due to the use of a finite integration timestep. Sivak et al., introduced a convenient approach to approximating a natural measure of error between the sampled density and the target equilibrium density, the Kullback-Leibler (KL) divergence, in phase space, but did not specifically address the issue of configuration-space properties, which are much more commonly of interest in molecular simulations. Here, we introduce a variant of this near-equilibrium estimator capable of measuring the error in the configuration-space marginal density, validating it against a complex but exact nested Monte Carlo estimator to show that it reproduces the KL divergence with high fidelity. To illustrate its utility, we employ this new near-equilibrium estimator to assess a claim that a recently proposed Langevin integrator introduces extremely small configuration-space density errors up to the stability limit at no extra computational expense. Finally, we show how this approach to quantifying sampling bias can be applied to a wide variety of stochastic integrators by following a straightforward procedure to compute the appropriate shadow work, and describe how it can be extended to quantify the error in arbitrary marginal or conditional distributions of interest. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Graphical abstract

702 KiB  
Article
Langevin Dynamics with Variable Coefficients and Nonconservative Forces: From Stationary States to Numerical Methods
by Matthias Sachs, Benedict Leimkuhler and Vincent Danos
Entropy 2017, 19(12), 647; https://doi.org/10.3390/e19120647 - 29 Nov 2017
Cited by 29 | Viewed by 8792
Abstract
Langevin dynamics is a versatile stochastic model used in biology, chemistry, engineering, physics and computer science. Traditionally, in thermal equilibrium, one assumes (i) the forces are given as the gradient of a potential and (ii) a fluctuation-dissipation relation holds between stochastic and dissipative [...] Read more.
Langevin dynamics is a versatile stochastic model used in biology, chemistry, engineering, physics and computer science. Traditionally, in thermal equilibrium, one assumes (i) the forces are given as the gradient of a potential and (ii) a fluctuation-dissipation relation holds between stochastic and dissipative forces; these assumptions ensure that the system samples a prescribed invariant Gibbs-Boltzmann distribution for a specified target temperature. In this article, we relax these assumptions, incorporating variable friction and temperature parameters and allowing nonconservative force fields, for which the form of the stationary state is typically not known a priori. We examine theoretical issues such as stability of the steady state and ergodic properties, as well as practical aspects such as the design of numerical methods for stochastic particle models. Applications to nonequilibrium systems with thermal gradients and active particles are discussed. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Graphical abstract

3738 KiB  
Article
Statistical Measures to Quantify Similarity between Molecular Dynamics Simulation Trajectories
by Jenny Farmer, Fareeha Kanwal, Nikita Nikulsin, Matthew C. B. Tsilimigras and Donald J. Jacobs
Entropy 2017, 19(12), 646; https://doi.org/10.3390/e19120646 - 29 Nov 2017
Cited by 32 | Viewed by 9084
Abstract
Molecular dynamics simulation is commonly employed to explore protein dynamics. Despite the disparate timescales between functional mechanisms and molecular dynamics (MD) trajectories, functional differences are often inferred from differences in conformational ensembles between two proteins in structure-function studies that investigate the effect of [...] Read more.
Molecular dynamics simulation is commonly employed to explore protein dynamics. Despite the disparate timescales between functional mechanisms and molecular dynamics (MD) trajectories, functional differences are often inferred from differences in conformational ensembles between two proteins in structure-function studies that investigate the effect of mutations. A common measure to quantify differences in dynamics is the root mean square fluctuation (RMSF) about the average position of residues defined by C α -atoms. Using six MD trajectories describing three native/mutant pairs of beta-lactamase, we make comparisons with additional measures that include Jensen-Shannon, modifications of Kullback-Leibler divergence, and local p-values from 1-sample Kolmogorov-Smirnov tests. These additional measures require knowing a probability density function, which we estimate by using a nonparametric maximum entropy method that quantifies rare events well. The same measures are applied to distance fluctuations between C α -atom pairs. Results from several implementations for quantitative comparison of a pair of MD trajectories are made based on fluctuations for on-residue and residue-residue local dynamics. We conclude that there is almost always a statistically significant difference between pairs of 100 ns all-atom simulations on moderate-sized proteins as evident from extraordinarily low p-values. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

479 KiB  
Article
Variational Characterization of Free Energy: Theory and Algorithms
by Carsten Hartmann, Lorenz Richter, Christof Schütte and Wei Zhang
Entropy 2017, 19(11), 626; https://doi.org/10.3390/e19110626 - 20 Nov 2017
Cited by 30 | Viewed by 7585
Abstract
The article surveys and extends variational formulations of the thermodynamic free energy and discusses their information-theoretic content from the perspective of mathematical statistics. We revisit the well-known Jarzynski equality for nonequilibrium free energy sampling within the framework of importance sampling and Girsanov change-of-measure [...] Read more.
The article surveys and extends variational formulations of the thermodynamic free energy and discusses their information-theoretic content from the perspective of mathematical statistics. We revisit the well-known Jarzynski equality for nonequilibrium free energy sampling within the framework of importance sampling and Girsanov change-of-measure transformations. The implications of the different variational formulations for designing efficient stochastic optimization and nonequilibrium simulation algorithms for computing free energies are discussed and illustrated. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

883 KiB  
Article
Transport Coefficients from Large Deviation Functions
by Chloe Ya Gao and David T. Limmer
Entropy 2017, 19(11), 571; https://doi.org/10.3390/e19110571 - 25 Oct 2017
Cited by 17 | Viewed by 5383
Abstract
We describe a method for computing transport coefficients from the direct evaluation of large deviation functions. This method is general, relying on only equilibrium fluctuations, and is statistically efficient, employing trajectory based importance sampling. Equilibrium fluctuations of molecular currents are characterized by their [...] Read more.
We describe a method for computing transport coefficients from the direct evaluation of large deviation functions. This method is general, relying on only equilibrium fluctuations, and is statistically efficient, employing trajectory based importance sampling. Equilibrium fluctuations of molecular currents are characterized by their large deviation functions, which are scaled cumulant generating functions analogous to the free energies. A diffusion Monte Carlo algorithm is used to evaluate the large deviation functions, from which arbitrary transport coefficients are derivable. We find significant statistical improvement over traditional Green–Kubo based calculations. The systematic and statistical errors of this method are analyzed in the context of specific transport coefficient calculations, including the shear viscosity, interfacial friction coefficient, and thermal conductivity. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

506 KiB  
Article
Symmetries and Geometrical Properties of Dynamical Fluctuations in Molecular Dynamics
by Robert L. Jack, Marcus Kaiser and Johannes Zimmer
Entropy 2017, 19(10), 562; https://doi.org/10.3390/e19100562 - 22 Oct 2017
Cited by 5 | Viewed by 4655
Abstract
We describe some general results that constrain the dynamical fluctuations that can occur in non-equilibrium steady states, with a focus on molecular dynamics. That is, we consider Hamiltonian systems, coupled to external heat baths, and driven out of equilibrium by non-conservative forces. We [...] Read more.
We describe some general results that constrain the dynamical fluctuations that can occur in non-equilibrium steady states, with a focus on molecular dynamics. That is, we consider Hamiltonian systems, coupled to external heat baths, and driven out of equilibrium by non-conservative forces. We focus on the probabilities of rare events (large deviations). First, we discuss a PT (parity-time) symmetry that appears in ensembles of trajectories where a current is constrained to have a large (non-typical) value. We analyse the heat flow in such ensembles, and compare it with non-equilibrium steady states. Second, we consider pathwise large deviations that are defined by considering many copies of a system. We show how the probability currents in such systems can be decomposed into orthogonal contributions that are related to convergence to equilibrium and to dissipation. We discuss the implications of these results for modelling non-equilibrium steady states. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

414 KiB  
Article
Comparing Markov Chain Samplers for Molecular Simulation
by Robert D. Skeel and Youhan Fang
Entropy 2017, 19(10), 561; https://doi.org/10.3390/e19100561 - 21 Oct 2017
Cited by 6 | Viewed by 4539
Abstract
Markov chain Monte Carlo sampling propagators, including numerical integrators for stochastic dynamics, are central to the calculation of thermodynamic quantities and determination of structure for molecular systems. Efficiency is paramount, and to a great extent, this is determined by the integrated autocorrelation time [...] Read more.
Markov chain Monte Carlo sampling propagators, including numerical integrators for stochastic dynamics, are central to the calculation of thermodynamic quantities and determination of structure for molecular systems. Efficiency is paramount, and to a great extent, this is determined by the integrated autocorrelation time (IAcT). This quantity varies depending on the observable that is being estimated. It is suggested that it is the maximum of the IAcT over all observables that is the relevant metric. Reviewed here is a method for estimating this quantity. For reversible propagators (which are those that satisfy detailed balance), the maximum IAcT is determined by the spectral gap in the forward transfer operator, but for irreversible propagators, the maximum IAcT can be far less than or greater than what might be inferred from the spectral gap. This is consistent with recent theoretical results (not to mention past practical experience) suggesting that irreversible propagators generally perform better if not much better than reversible ones. Typical irreversible propagators have a parameter controlling the mix of ballistic and diffusive movement. To gain insight into the effect of the damping parameter for Langevin dynamics, its optimal value is obtained here for a multidimensional quadratic potential energy function. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

515 KiB  
Article
Reliable Approximation of Long Relaxation Timescales in Molecular Dynamics
by Wei Zhang and Christof Schütte
Entropy 2017, 19(7), 367; https://doi.org/10.3390/e19070367 - 18 Jul 2017
Cited by 10 | Viewed by 5693
Abstract
Many interesting rare events in molecular systems, like ligand association, protein folding or conformational changes, occur on timescales that often are not accessible by direct numerical simulation. Therefore, rare event approximation approaches like interface sampling, Markov state model building, or advanced reaction coordinate-based [...] Read more.
Many interesting rare events in molecular systems, like ligand association, protein folding or conformational changes, occur on timescales that often are not accessible by direct numerical simulation. Therefore, rare event approximation approaches like interface sampling, Markov state model building, or advanced reaction coordinate-based free energy estimation have attracted huge attention recently. In this article we analyze the reliability of such approaches. How precise is an estimate of long relaxation timescales of molecular systems resulting from various forms of rare event approximation methods? Our results give a theoretical answer to this question by relating it with the transfer operator approach to molecular dynamics. By doing so we also allow for understanding deep connections between the different approaches. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

7112 KiB  
Article
An Exploration Algorithm for Stochastic Simulators Driven by Energy Gradients
by Anastasia S. Georgiou, Juan M. Bello-Rivas, Charles William Gear, Hau-Tieng Wu, Eliodoro Chiavazzo and Ioannis G. Kevrekidis
Entropy 2017, 19(7), 294; https://doi.org/10.3390/e19070294 - 22 Jun 2017
Cited by 8 | Viewed by 5957
Abstract
In recent work, we have illustrated the construction of an exploration geometry on free energy surfaces: the adaptive computer-assisted discovery of an approximate low-dimensional manifold on which the effective dynamics of the system evolves. Constructing such an exploration geometry involves geometry-biased sampling (through [...] Read more.
In recent work, we have illustrated the construction of an exploration geometry on free energy surfaces: the adaptive computer-assisted discovery of an approximate low-dimensional manifold on which the effective dynamics of the system evolves. Constructing such an exploration geometry involves geometry-biased sampling (through both appropriately-initialized unbiased molecular dynamics and through restraining potentials) and, machine learning techniques to organize the intrinsic geometry of the data resulting from the sampling (in particular, diffusion maps, possibly enhanced through the appropriate Mahalanobis-type metric). In this contribution, we detail a method for exploring the conformational space of a stochastic gradient system whose effective free energy surface depends on a smaller number of degrees of freedom than the dimension of the phase space. Our approach comprises two steps. First, we study the local geometry of the free energy landscape using diffusion maps on samples computed through stochastic dynamics. This allows us to automatically identify the relevant coarse variables. Next, we use the information garnered in the previous step to construct a new set of initial conditions for subsequent trajectories. These initial conditions are computed so as to explore the accessible conformational space more efficiently than by continuing the previous, unbiased simulations. We showcase this method on a representative test system. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

2600 KiB  
Article
Ion Hopping and Constrained Li Diffusion Pathways in the Superionic State of Antifluorite Li2O
by Ajay Annamareddy and Jacob Eapen
Entropy 2017, 19(5), 227; https://doi.org/10.3390/e19050227 - 18 May 2017
Cited by 10 | Viewed by 7677
Abstract
Li2O belongs to the family of antifluorites that show superionic behavior at high temperatures. While some of the superionic characteristics of Li2O are well-known, the mechanistic details of ionic conduction processes are somewhat nebulous. In this work, we first [...] Read more.
Li2O belongs to the family of antifluorites that show superionic behavior at high temperatures. While some of the superionic characteristics of Li2O are well-known, the mechanistic details of ionic conduction processes are somewhat nebulous. In this work, we first establish an onset of superionic conduction that is emblematic of a gradual disordering process among the Li ions at a characteristic temperature Tα (~1000 K) using reported neutron diffraction data and atomistic simulations. In the superionic state, the Li ions are observed to portray dynamic disorder by hopping between the tetrahedral lattice sites. We then show that string-like ionic diffusion pathways are established among the Li ions in the superionic state. The diffusivity of these dynamical string-like structures, which have a finite lifetime, shows a remarkable correlation to the bulk diffusivity of the system. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

2983 KiB  
Article
Calculating Iso-Committor Surfaces as Optimal Reaction Coordinates with Milestoning
by Ron Elber, Juan M. Bello-Rivas, Piao Ma, Alfredo E. Cardenas and Arman Fathizadeh
Entropy 2017, 19(5), 219; https://doi.org/10.3390/e19050219 - 11 May 2017
Cited by 46 | Viewed by 6835
Abstract
Reaction coordinates are vital tools for qualitative and quantitative analysis of molecular processes. They provide a simple picture of reaction progress and essential input for calculations of free energies and rates. Iso-committor surfaces are considered the optimal reaction coordinate. We present an algorithm [...] Read more.
Reaction coordinates are vital tools for qualitative and quantitative analysis of molecular processes. They provide a simple picture of reaction progress and essential input for calculations of free energies and rates. Iso-committor surfaces are considered the optimal reaction coordinate. We present an algorithm to compute efficiently a sequence of isocommittor surfaces. These surfaces are considered an optimal reaction coordinate. The algorithm analyzes Milestoning results to determine the committor function. It requires only the transition probabilities between the milestones, and not transition times. We discuss the following numerical examples: (i) a transition in the Mueller potential; (ii) a conformational change of a solvated peptide; and (iii) cholesterol aggregation in membranes. Full article
(This article belongs to the Special Issue Understanding Molecular Dynamics via Stochastic Processes)
Show Figures

Figure 1

Back to TopTop