Coarse-grained stochastic processes and kinetic Monte Carlo simulators for the diffusion of interacting particles

3D Off-Lattice Coarse-Grained Monte Carlo Simulations for Nucleation of Alkaline Aluminosilicate Gels

This work presents a 3D off-lattice coarse-grained Monte Carlo (CGMC) approach to simulate the nucleation of alkaline aluminosilicate gels, their nanostructure particle size, and their pore size distribution. In this model, four monomer species are coarse-grained with different particle sizes. The novelty is extending the previous on-lattice approach from White et al. (2012 and 2020) by implementing a full off-lattice numerical implementation to consider tetrahedral geometrical constraints when aggregating the particles into clusters. Aggregation of the dissolved silicate and aluminate monomers was simulated until reaching the equilibrium condition of 16.46% and 17.04% in particle number, respectively. The cluster size formation was analyzed as a function of iteration step evolution. The obtained equilibrated nano-structure was digitized to obtain the pore size distribution and this was compared with the on-lattice CGMC and measurement results from White et al. The observed difference highlighted the importance of the developed off-lattice CGMC approach to better describe the nanostructure of aluminosilicate gels.

Emergence of sector and spiral patterns from a two-species mutualistic cross-feeding model

The ubiquitous existence of microbial communities marks the importance of understanding how species interact within the community to coexist and their spatial organization. We study a two-species mutualistic cross-feeding model through a stochastic cellular automaton on a square lattice using kinetic Monte Carlo simulation. Our model encapsulates the essential dynamic processes such as cell growth, and nutrient excretion, diffusion and uptake. Focusing on the interplay among nutrient diffusion and individual cell division, we discover three general classes of colony morphology: co-existing sectors, co-existing spirals, and engulfment. When the cross-feeding nutrient is widely available, either through high excretion or fast diffusion, a stable circular colony with alternating species sector emerges. When the consumer cells rely on being spatially close to the producers, we observe a stable spiral. We also see one species being engulfed by the other when species interfaces merge due to stochastic fluctuation. By tuning the diffusion rate and the growth rate, we are able to gain quantitative insights into the structures of the sectors and the spirals.

Acceleration scheme for particle transport in kinetic Monte Carlo methods

Kinetic Monte Carlo (kMC) simulations are frequently used to study (electro-)chemical processes within science and engineering. kMC methods provide insight into the interplay of stochastic processes and can link atomistic material properties with macroscopic characteristics. Significant problems concerning the computational demand arise if processes with large time disparities are competing. Acceleration algorithms are required to make slow processes accessible. Especially, the accelerated superbasin kMC (AS-kMC) scheme has been frequently applied within chemical reaction networks. For larger systems, the computational overhead of the AS-kMC is significant as the computation of the superbasins is done during runtime and comes with the need for large databases. Here, we propose a novel acceleration scheme for diffusion and transport processes within kMC simulations. Critical superbasins are detected during the system initialization. Scaling factors for the critical rates within the superbasins, as well as a lower bound for the number of sightings, are derived. Our algorithm exceeds the AS-kMC in the required simulation time, which we demonstrate with a 1D-chain example. In addition, we apply the acceleration scheme to study the time-of-flight (TOF) of charge carriers within organic semiconductors. In this material class, time disparities arise due to a significant spread of transition rates. The acceleration scheme allows a significant acceleration up to a factor of 65 while keeping the error of the TOF values negligible. The computational overhead is negligible, as all superbasins only need to be computed once.

A fast species redistribution approach to accelerate the kinetic Monte Carlo simulation for heterogeneous catalysis

The first-principles kinetic Monte Carlo (kMC) simulation has been demonstrated as a reliable multiscale modeling approach in silico to disclose the interplay among all the elementary steps in a complex reaction network for heterogeneous catalysis. Heterogeneous catalytic systems frequently contain fast surface diffusion processes of some adsorbates while the elementary steps in it would be much slower than those in fast diffusion. Consequently, the kMC simulation for these systems is easily trapped in the sub-basins of a super basin on a potential energy surface due to the continuous and repeated sampling of these fast processes, which would significantly increase the total accessible simulation time and even make it impossible to get the reasonable simulation results using the kMC simulation. In this work, we present an improved fast species redistribution (FSR) method for the kMC simulation to overcome the stiffness problem resulting from the low-barrier surface diffusion to accelerate the heterogeneous catalytic kMC simulation. Taking CO oxidations on Pt(111) and Pt(100) as examples, we demonstrate that the FSR approach can properly reproduce the results of an equivalent first-principles microkinetic model simulation with more reasonable reaction rates. The improved kMC simulation based on the FSR method can accurately incorporate the effect of the fast diffusion of species on the surface and provide several orders of magnitude of acceleration compared to the standard kMC simulation.

A lattice kinetic Monte-Carlo method for simulating chromosomal dynamics and other (non-)equilibrium bio-assemblies

Biological assemblies in living cells such as chromosomes constitute large many-body systems that operate in a fluctuating, out-of-equilibrium environment. Since a brute-force simulation of that many degrees of freedom is currently computationally unfeasible, it is necessary to perform coarse-grained stochastic simulations. Here, we develop all tools necessary to write a lattice kinetic Monte-Carlo (LKMC) algorithm capable of performing such simulations. We discuss the validity and limits of this approach by testing the results of the simulation method in simple settings. Importantly, we illustrate how at large external forces Metropolis-Hastings kinetics violate the fluctuation-dissipation and steady-state fluctuation theorems and discuss better alternatives. Although this simulation framework is rather general, we demonstrate our approach using a DNA polymer with interacting SMC condensin loop-extruding enzymes. Specifically, we show that the scaling behavior of the loop-size distributions that we obtain in our LKMC simulations of this SMC-DNA system is consistent with that reported in other studies using Brownian dynamics simulations and analytic approaches. Moreover, we find that the irreversible dynamics of these enzymes under certain conditions result in frozen, sterically jammed polymer configurations, highlighting a potential pitfall of this approach.

Analysis of the lattice kinetic Monte Carlo method in systems with external fields

The lattice kinetic Monte Carlo (LKMC) method is studied in the context of Brownian particles subjected to drift forces, here principally represented by external fluid flow. LKMC rate expressions for particle hopping are derived that satisfy detailed balance at equilibrium while also providing correct dynamical trajectories in advective-diffusive situations. Error analyses are performed for systems in which collections of particles undergo Brownian motion while also being advected by plug and parabolic flows. We demonstrate how the flow intensity, and its associated drift force, as well as its gradient, each impact the accuracy of the method in relation to reference analytical solutions and Brownian dynamics simulations. Finally, we show how a non-uniform grid that everywhere retains full microscopic detail may be employed to increase the computational efficiency of lattice kinetic Monte Carlo simulations of particles subjected to drift forces arising from the presence of external fields.

Adaptive coarse-grained Monte Carlo simulation of reaction and diffusion dynamics in heterogeneous plasma membranes

Background

An adaptive coarse-grained (kinetic) Monte Carlo (ACGMC) simulation framework is applied to reaction and diffusion dynamics in inhomogeneous domains. The presented model is relevant to the diffusion and dimerization dynamics of epidermal growth factor receptor (EGFR) in the presence of plasma membrane heterogeneity and specifically receptor clustering. We perform simulations representing EGFR cluster dissipation in heterogeneous plasma membranes consisting of higher density clusters of receptors surrounded by low population areas using the ACGMC method. We further investigate the effect of key parameters on the cluster lifetime.

Results

Coarse-graining of dimerization, rather than of diffusion, may lead to computational error. It is shown that the ACGMC method is an effective technique to minimize error in diffusion-reaction processes and is superior to the microscopic kinetic Monte Carlo simulation in terms of computational cost while retaining accuracy. The low computational cost enables sensitivity analysis calculations. Sensitivity analysis indicates that it may be possible to retain clusters of receptors over the time scale of minutes under suitable conditions and the cluster lifetime may depend on both receptor density and cluster size.

Conclusions

The ACGMC method is an ideal platform to resolve large length and time scales in heterogeneous biological systems well beyond the plasma membrane and the EGFR system studied here. Our results demonstrate that cluster size must be considered in conjunction with receptor density, as they synergistically affect EGFR cluster lifetime. Further, the cluster lifetime being of the order of several seconds suggests that any mechanisms responsible for EGFR aggregation must operate on shorter timescales (at most a fraction of a second), to overcome dissipation and produce stable clusters observed experimentally.

Coarse-grained kinetic Monte Carlo models: Complex lattices, multicomponent systems, and homogenization at the stochastic level

On-lattice kinetic Monte Carlo (KMC) simulations have extensively been applied to numerous systems. However, their applicability is severely limited to relatively short time and length scales. Recently, the coarse-grained MC (CGMC) method was introduced to greatly expand the reach of the lattice KMC technique. Herein, we extend the previous spatial CGMC methods to multicomponent species and/or site types. The underlying theory is derived and numerical examples are presented to demonstrate the method. Furthermore, we introduce the concept of homogenization at the stochastic level over all site types of a spatially coarse-grained cell. Homogenization provides a novel coarsening of the number of processes, an important aspect for complex problems plagued by the existence of numerous microscopic processes (combinatorial complexity). As expected, the homogenized CGMC method outperforms the traditional KMC method on computational cost while retaining good accuracy.

Temporal coarse-graining of microscopic-lattice kinetic Monte Carlo simulations via tau leaping

A coarse-time-step method is presented that enables the execution of multiple events at each time increment of microscopic-lattice kinetic Monte Carlo simulations. The method employs the n-fold method to create groups of reactions in which the tau-leap algorithm of Gillespie, originally proposed for well-mixed systems, is applied. Creation of groups of reactions is an essential step to avoid violation of the leap condition that arises when the tau-leap algorithm is applied to a single site. The method is general, very easy to implement, and can result in substantial computational savings when global updating is employed. An illustrative example from crystal growth of a simple cubic lattice with the solid-on-solid approximation is presented.

An applied mathematics perspective on stochastic modelling for climate

Systematic strategies from applied mathematics for stochastic modelling in climate are reviewed here. One of the topics discussed is the stochastic modelling of mid-latitude low-frequency variability through a few teleconnection patterns, including the central role and physical mechanisms responsible for multiplicative noise. A new low-dimensional stochastic model is developed here, which mimics key features of atmospheric general circulation models, to test the fidelity of stochastic mode reduction procedures. The second topic discussed here is the systematic design of stochastic lattice models to capture irregular and highly intermittent features that are not resolved by a deterministic parametrization. A recent applied mathematics design principle for stochastic column modelling with intermittency is illustrated in an idealized setting for deep tropical convection; the practical effect of this stochastic model in both slowing down convectively coupled waves and increasing their fluctuations is presented here.

Coarse-grained lattice kinetic Monte Carlo simulation of systems of strongly interacting particles

A general approach is presented for spatially coarse-graining lattice kinetic Monte Carlo (LKMC) simulations of systems containing strongly interacting particles. While previous work has relied on approximations that are valid in the limit of weak interactions, here we show that it is possible to compute coarse-grained transition rates for strongly interacting systems without a large computational burden. A two-dimensional square lattice is employed on which a collection of (supersaturated) strongly interacting particles is allowed to reversibly evolve into clusters. A detailed analysis is presented of the various approximations applied in LKMC coarse graining, and a number of numerical closure rules are contrasted and compared. In each case, the overall cluster size distribution and individual cluster structures are used to assess the accuracy of the coarse-graining approach. The resulting closure approach is shown to provide an excellent coarse-grained representation of the systems considered in this study.

Acceleration of Monte Carlo simulations through spatial updating in the grand canonical ensemble

A new grand canonical Monte Carlo algorithm for continuum fluid models is proposed. The method is based on a generalization of sequential Monte Carlo algorithms for lattice gas systems. The elementary moves, particle insertions and removals, are constructed by analogy with those of a lattice gas. The updating is implemented by selecting points in space (spatial updating) either at random or in a definitive order (sequential). The type of move, insertion or removal, is deduced based on the local environment of the selected points. Results on two-dimensional square-well fluids indicate that the sequential version of the proposed algorithm converges faster than standard grand canonical algorithms for continuum fluids. Due to the nature of the updating, additional reduction of simulation time may be achieved by parallel implementation through domain decomposition.

Continuum mesoscopic framework for multiple interacting species and processes on multiple site types and/or crystallographic planes

While recently derived continuum mesoscopic equations successfully bridge the gap between microscopic and macroscopic physics, so far they have been derived only for simple lattice models. In this paper, general deterministic continuum mesoscopic equations are derived rigorously via nonequilibrium statistical mechanics to account for multiple interacting surface species and multiple processes on multiple site types and/or different crystallographic planes. Adsorption, desorption, reaction, and surface diffusion are modeled. It is demonstrated that contrary to conventional phenomenological continuum models, microscopic physics, such as the interaction potential, determines the final form of the mesoscopic equation. Models of single component diffusion and binary diffusion of interacting particles on single-type site lattice and of single component diffusion on complex microporous materials' lattices consisting of two types of sites are derived, as illustrations of the mesoscopic framework. Simplification of the diffusion mesoscopic model illustrates the relation to phenomenological models, such as the Fickian and Maxwell-Stefan transport models. It is demonstrated that the mesoscopic equations are in good agreement with lattice kinetic Monte Carlo simulations for several prototype examples studied.

Spectral method study of domain coarsening

The self-organization of particles in a two phase system in the coexistence region through a diffusive mechanism is known as Ostwald ripening. The late stage of Ostwald ripening is studied here through the use of a mesoscopic model in conjunction with a special configuration that allows for the direct measurement of system characteristics such as domain size. The mesoscopic model is a stochastic partial integrodifferential equation and is studied through the use of recently developed and benchmarked spectral schemes. The size of the growing region is not observed to increase as a power law in this model during the late stages of self-organization in contrast to the power law growth observed in simulations of the earlier stages of self-organization. The results included here also demonstrate the effect of adjusting the interparticle interaction on the morphological evolution of the system.

Multiscale spatial Monte Carlo simulations: Multigriding, computational singular perturbation, and hierarchical stochastic closures

Monte Carlo (MC) simulation of most spatially distributed systems is plagued by several problems, namely, execution of one process at a time, large separation of time scales of various processes, and large length scales. Recently, a coarse-grained Monte Carlo (CGMC) method was introduced that can capture large length scales at reasonable computational times. An inherent assumption in this CGMC method revolves around a mean-field closure invoked in each coarse cell that is inaccurate for short-ranged interactions. Two new approaches are explored to improve upon this closure. The first employs the local quasichemical approximation, which is applicable to first nearest-neighbor interactions. The second, termed multiscale CGMC method, employs singular perturbation ideas on multiple grids to capture the entire cluster probability distribution function via short microscopic MC simulations on small, fine-grid lattices by taking advantage of the time scale separation of multiple processes. Computational strategies for coupling the fast process at small length scales (fine grid) with the slow processes at large length scales (coarse grid) are discussed. Finally, the binomial tau-leap method is combined with the multiscale CGMC method to execute multiple processes over the entire lattice and provide additional computational acceleration. Numerical simulations demonstrate that in the presence of fast diffusion and slow adsorption and desorption processes the two new approaches provide more accurate solutions in comparison to the previously introduced CGMC method.

Acceleration of Markov chain Monte Carlo simulations through sequential updating

Strict detailed balance is not necessary for Markov chain Monte Carlo simulations to converge to the correct equilibrium distribution. In this work, we propose a new algorithm which only satisfies the weaker balance condition, and it is shown analytically to have better mobility over the phase space than the Metropolis algorithm satisfying strict detailed balance. The new algorithm employs sequential updating and yields better sampling statistics than the Metropolis algorithm with random updating. We illustrate the efficiency of the new algorithm on the two-dimensional Ising model. The algorithm is shown to identify the correct equilibrium distribution and to converge faster than the Metropolis algorithm with strict detailed balance. The main advantages of the new algorithm are its simplicity and the feasibility of parallel implementation through domain decomposition.

The significance of lipid composition for membrane activity: New concepts and ways of assessing function

In the last decade or so, it has been realised that membranes do not just have a lipid-bilayer structure in which proteins are embedded or with which they associate. Structures are dynamic and contain areas of heterogeneity which are vital for their formation. In this review, we discuss some of the ways in which these dynamic and heterogeneous structures have implications during stress and in relation to certain human diseases. A particular stress is that of temperature which may instigate adaptation in poikilotherms or appropriate defensive responses during fever in mammals. Recent data emphasise the role of membranes in sensing temperature changes and in controlling a regulatory loop with chaperone proteins. This loop seems to need the existence of specific membrane microdomains and also includes association of chaperone (heat stress) proteins with the membrane. The role of microdomains is then discussed further in relation to various human pathologies such as cardiovascular disease, cancer and neurodegenerative diseases. The concept of modifying membrane lipids (lipid therapy) as a means for treating such pathologies is then introduced. Examples are given when such methods have been shown to have benefit. In order to study membrane microheterogeneity in detail and to elucidate possible molecular mechanisms that account for alteration in membrane function, new methods are needed. In the second part of the review, we discuss ultra-sensitive and ultra-resolution imaging techniques. These include atomic force microscopy, single particle tracking, single particle tracing and various modern fluorescence methods. Finally, we deal with computing simulation of membrane systems. Such methods include coarse-grain techniques and Monte Carlo which offer further advances into molecular dynamics. As computational methods advance they will have more application by revealing the very subtle interactions that take place between the lipid and protein components of membranes - and which are so essential to their function.

Spatially adaptive grand canonical ensemble Monte Carlo simulations

A spatially adaptive Monte Carlo method is introduced directly from the underlying microscopic mechanisms, which satisfies detailed balance, gives the correct noise, and describes accurately dynamic and equilibrium states for adsorption-desorption (grand canonical ensemble) processes. It enables simulations of large scales while capturing sharp gradients with molecular resolution at significantly reduced computational cost. A posteriori estimates, in the sense used in finite-elements methods, are developed for assessing errors (information loss) in coarse-graining and guiding mesh generation.

Binomial distribution based τ-leap accelerated stochastic simulation

Recently, Gillespie introduced the tau-leap approximate, accelerated stochastic Monte Carlo method for well-mixed reacting systems [J. Chem. Phys. 115, 1716 (2001)]. In each time increment of that method, one executes a number of reaction events, selected randomly from a Poisson distribution, to enable simulation of long times. Here we introduce a binomial distribution tau-leap algorithm (abbreviated as BD-tau method). This method combines the bounded nature of the binomial distribution variable with the limiting reactant and constrained firing concepts to avoid negative populations encountered in the original tau-leap method of Gillespie for large time increments, and thus conserve mass. Simulations using prototype reaction networks show that the BD-tau method is more accurate than the original method for comparable coarse-graining in time.

Promoter-induced reactive phase separation in surface reactions

Promoters are adsorbed mobile species which do not directly participate in a catalytic surface reaction, but can influence its rate. Often, they are characterized by strong attractive interactions with one of the reactants. We show that these conditions lead to a Turing instability of the uniform state and to the formation of reaction-induced periodic concentration patterns. Experimentally such patterns are observed in catalytic water formation on a Rh(110) surface in the presence of coadsorbed potassium.