Coarse-Grained Kinetic Computations for Rare Events: Application to Micelle Formation

Abstract

We discuss a coarse-grained approach to the computation of rare events in the context of grand canonical Monte Carlo (GCMC) simulations of self-assembly of surfactant molecules into micelles. The basic assumption is that the {\it computational} system dynamics can be decomposed into two parts -- fast (noise) and slow (reaction coordinates) dynamics, so that the system can be described by an effective, coarse grained Fokker-Planck (FP) equation. While such an assumption may be valid in many circumstances, an explicit form of FP equation is not always available. In our computations we bypass the analytic derivation of such an effective FP equation. The effective free energy gradient and the state-dependent magnitude of the random noise, which are necessary to formulate the effective Fokker-Planck equation, are obtained from ensembles of short bursts of microscopic simulations {\it with judiciously chosen initial conditions}. The reaction coordinate in our micelle formation problem is taken to be the size of a cluster of surfactant molecules. We test the validity of the effective FP description in this system and reconstruct a coarse-grained free energy surface in good agreement with full-scale GCMC simulations. We also show that, for very small clusters, the cluster size seizes to be a good reaction coordinate for a one-dimensional effective description. We discuss possible ways to improve the current model and to take higher-dimensional coarse-grained dynamics into account.

References in corpus (5)

• Escaping free-energy minima
• Coarse Molecular Dynamics of a Peptide Fragment: Free Energy, Kinetics, and Long-Time Dynamics Computations
• Coarse Stability and Bifurcation Analysis Using Stochastic Simulators: Kinetic Monte Carlo Examples
• Three eras of micellization
• Coarse Grained Computations for a Micellar System

Cited by in corpus (14)

• Gene regulatory networks: a coarse-grained, equation-free approach to multiscale computation
• Stochastic Physics-Informed Neural Ordinary Differential Equations
• Variable-free exploration of stochastic models: a gene regulatory network example
• Optimal Langevin modelling of out-of-equilibrium molecular dynamics simulations
• Statistically optimal continuous free energy surfaces from biased simulations and multistate reweighting
• Coarse Grained Computations for a Micellar System
• Learning the temporal evolution of multivariate densities via normalizing flows
• Fitting Effective Diffusion Models to Data Associated with a "Glassy Potential": Estimation, Classical Inference Procedures and Some Heuristics
• Extracting Stochastic Governing Laws by Nonlocal Kramers-Moyal Formulas
• Coarse molecular-dynamics analysis of stress-induced structural transitions in crystals
• Coarse Molecular-Dynamics Determination of the Onset of Structural Transitions: Melting of Crystalline Solids
• Stable developmental patterns of gene expression without morphogen gradients
• A Random Force is a Force, of Course, of Coarse: Decomposing Complex Enzyme Kinetics with Surrogate Models
• Quantifying biomolecular diffusion with a "spherical cow" model