Assessing Position-Dependent Diffusion from Biased Simulations and Markov State Model Analysis

TL;DR

This paper introduces a new automated method to accurately estimate position-dependent diffusion coefficients from complex molecular simulation data, improving understanding of kinetic properties in high-energy barrier systems.

Contribution

The authors develop a novel, general approach for assessing multidimensional, position-dependent diffusion within Markovian frameworks, applicable to complex molecular and porous media systems.

Findings

Method accurately estimates diffusion coefficients in analytic potentials.

Algorithm outperforms standard methods in heterogeneous porous media.

Application to molecular dynamics data demonstrates practical utility.

Abstract

A variety of enhanced statistical and numerical methods are now routinely used to extract comprehensible and relevant thermodynamic information from the vast amount of complex, high-dimensional data obtained from intensive molecular simulations. The characterization of kinetic properties, such as diffusion coefficients, of molecular systems with significantly high energy barriers, on the other hand, has received less attention. Among others, Markov state models, in which the long-time statistical dynamics of a system is approximated by a Markov chain on a discrete partition of configuration space, have seen widespread use in recent years, with the aim of tackling these fundamental issues. Here we propose a general, automatic method to assess multidimensional position-dependent diffusion coefficients within the framework of Markovian stochastic processes and Kramers-Moyal expansion. We apply the formalism to one- and two-dimensional analytic potentials and data from explicit solvent molecular dynamics simulations, including the water-mediated conformations of alanine dipeptide. Importantly, the developed algortihm presents significant improvement compared to standard methods when the transport of solute across three-dimensional heterogeneous porous media is studied, for example, the prediction of membrane permeation of drug molecules.