First-passage fingerprints of water diffusion near glutamine surfaces

Estimating Position-Dependent and Anisotropic Diffusivity Tensors from Molecular Dynamics Trajectories: Existing Methods and Future Outlook

Confinement can substantially alter the physicochemical properties of materials by breaking translational isotropy and rendering all physical properties position-dependent. Molecular dynamics (MD) simulations have proven instrumental in characterizing such spatial heterogeneities and probing the impact of confinement on materials' properties. For static properties, this is a straightforward task and can be achieved via simple spatial binning. Such an approach, however, cannot be readily applied to transport coefficients due to lack of natural extensions of autocorrelations used for their calculation in the bulk. The prime example of this challenge is diffusivity, which, in the bulk, can be readily estimated from the particles' mobility statistics, which satisfy the Fokker-Planck equation. Under confinement, however, such statistics will follow the Smoluchowski equation, which lacks a closed-form analytical solution. This brief review explores the rich history of estimating profiles of the diffusivity tensor from MD simulations and discusses various approximate methods and algorithms developed for this purpose. Besides discussing heuristic extensions of bulk methods, we overview more rigorous algorithms, including kernel-based methods, Bayesian approaches, and operator discretization techniques. Additionally, we outline methods based on applying biasing potentials or imposing constraints on tracer particles. Finally, we discuss approaches that estimate diffusivity from mean first passage time or committor probability profiles, a conceptual framework originally developed in the context of collective variable spaces describing rare events in computational chemistry and biology. In summary, this paper offers a concise survey of diverse approaches for estimating diffusivity from MD trajectories, highlighting challenges and opportunities in this area.

Controlling Uncertainty of Empirical First-Passage Times in the Small-Sample Regime

We derive general bounds on the probability that the empirical first-passage time τ[over ¯]_{n}≡∑_{i=1}^{n}τ_{i}/n of a reversible ergodic Markov process inferred from a sample of n independent realizations deviates from the true mean first-passage time by more than any given amount in either direction. We construct nonasymptotic confidence intervals that hold in the elusive small-sample regime and thus fill the gap between asymptotic methods and the Bayesian approach that is known to be sensitive to prior belief and tends to underestimate uncertainty in the small-sample setting. We prove sharp bounds on extreme first-passage times that control uncertainty even in cases where the mean alone does not sufficiently characterize the statistics. Our concentration-of-measure-based results allow for model-free error control and reliable error estimation in kinetic inference, and are thus important for the analysis of experimental and simulation data in the presence of limited sampling.

Statistical physics of inhomogeneous transport: Unification of diffusion laws and inference from first-passage statistics

Characterization of composite materials, whose properties vary in space over microscopic scales, has become a problem of broad interdisciplinary interest. In particular, estimation of the inhomogeneous transport coefficients, e.g., the diffusion coefficient or the heat conductivity, which shape important processes in biology and engineering, is a challenging task. The analysis of such systems is further complicated because two alternative formulations of the inhomogeneous transport equations exist in the literature-the Smoluchowski and Fokker-Planck equations, which are also related to the so-called Ito-Stratonovich dilemma. Using the theory of statistical physics, we show that the two formulations, usually regarded as distinct models, are physically equivalent. From this result we develop efficient estimates for the transverse space-dependent diffusion coefficient in fluids near a phase boundary. Our method requires only measurements of escape probabilities and mean exit times of molecules leaving a narrow spatial region. We test our estimates in three case studies: (i) a Langevin model of a Büttikker-Landauer ratchet; atomistic molecular-dynamics simulations of liquid-water molecules in contact with (ii) vapor, and (iii) soap (surfactant) film which has promising applications in physical chemistry. Our analysis reveals that near the surfactant monolayer the mobility of water molecules is slowed down almost twice with respect to the bulk liquid. Moreover, the diffusion coefficient of water correlates with the transition from hydrophilic to hydrophobic parts of the film.

Capture of a diffusive prey by multiple predators in confined space

The first passage search of a diffusing target (prey) by multiple searchers (predators) in confinement is an important problem in the stochastic process literature. While the analogous problem in open space has been studied in some detail, a systematic study in confined space is still lacking. In this paper, we study the first passage times for this problem in one, two, and three dimensions. Due to confinement, the survival probability of the target takes a form ∼e^{-t/τ} at large times t. The characteristic capture timescale τ associated with the rare capture events are rather challenging to measure. We use a computational algorithm that allows us to estimate τ with high accuracy. We study in detail the behavior of τ as a function of the system parameters, namely, the number of searchers N, the relative diffusivity r of the target with respect to the searcher, and the system size. We find that τ deviates from the ∼1/N scaling seen in the case of a static target, and this deviation varies continuously with r and the spatial dimensions.

Mean exit time and escape probability for the Ornstein-Uhlenbeck process

This paper studies the dynamics of the Ornstein-Uhlenbeck process by the deterministic quantities such as the mean exit time and escape probability. By solving the elliptic partial equations, we obtain explicit solutions to both mentioned problems using the special functions. We find that the mean exit time is longer for smaller noise, and the maximum depends on the middle of the given interval. Moreover, the likelihood that the solution orbits exiting the interval from left or right relies on the middle of the interval. The Monte Carlo simulations are carried out to support the obtained results.