BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Rosalind J. Allen

Rosalind J. Allen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechBiological Physicscond-mat.softPopulations and Evolutionphysics.comp-phQuantitative MethodsCell BehaviorMolecular NetworksBiomoleculescond-mat.mtrl-sciphysics.chem-phphysics.flu-dynphysics.optics

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
13topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2015arXiv

Shaping the Growth Behaviour of Biofilms Initiated from Bacterial Aggregates

Bacterial biofilms are usually assumed to originate from individual cells deposited on a surface. However, many biofilm-forming bacteria tend to aggregate in the planktonic phase so that it is possible that many natural and infectious biofilms originate wholly or partially from pre-formed cell aggregates. Here, we use agent-based computer simulations to investigate the role of pre-formed aggregates in biofilm development. Focusing on the initial shape the aggregate forms on the surface, we find that the degree of spreading of an aggregate on a surface can play an important role in determining its eventual fate during biofilm development. Specifically, initially spread aggregates perform better when competition with surrounding unaggregated bacterial cells is low, while initially rounded aggregates perform better when competition with surrounding unaggregated cells is high. These contrasting outcomes are governed by a trade-off between aggregate surface area and height. Our results provide new insight into biofilm formation and development, and reveal new factors that may be at play in the social evolution of biofilm communities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Growth-dependent bacterial susceptibility to ribosome-targeting antibiotics

Bacterial growth environment strongly influences the efficacy of antibiotic treatment, with slow growth often being associated with decreased susceptibility. Yet in many cases the connection between antibiotic susceptibility and pathogen physiology remains unclear. We show that for ribosome-targeting antibiotics acting on Escherichia coli, a complex interplay exists between physiology and antibiotic action; for some antibiotics within this class faster growth indeed increases susceptibility, but for other antibiotics the opposite is true. Remarkably, these observations can be explained by a simple mathematical model that combines drug transport and binding with physiological constraints. Our model reveals that growth-dependent susceptibility is controlled by a single parameter characterizing the `reversibility' of antibiotic transport and binding. This parameter provides a spectrum-classification of antibiotic growth-dependent efficacy that appears to correspond at its extremes to existing binary classification schemes. In these limits the model predicts universal, parameter-free limiting forms for growth inhibition curves. The model also leads to non-trivial predictions for the drug susceptibility of a translation-mutant strain of E. coli, which we verify experimentally. Drug action and bacterial metabolism are mechanistically complex; nevertheless this study illustrates how coarse-grained models can be used to integrate pathogen physiology into drug design and treatment strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inherent variability in the kinetics of autocatalytic protein self-assembly

In small volumes, the kinetics of filamentous protein self-assembly is expected to show significant variability, arising from intrinsic molecular noise. This is not accounted for in existing deterministic models. We introduce a simple stochastic model including nucleation and autocatalytic growth via elongation and fragmentation, which allows us to predict the effects of molecular noise on the kinetics of autocatalytic self-assembly. We derive an analytic expression for the lag-time distribution, which agrees well with experimental results for the fibrillation of bovine insulin. Our expression decomposes the lag time variability into contributions from primary nucleation and autocatalytic growth and reveals how each of these scales with the key kinetic parameters. Our analysis shows that significant lag-time variability can arise from both primary nucleation and from autocatalytic growth, and should provide a way to extract mechanistic information on early-stage aggregation from small-volume experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lower glycolysis carries a higher flux than any biochemically possible alternative

The universality of many pathways of core metabolism suggests a strong role for evolutionary selection, but it remains unclear whether existing pathways have been selected from a large or small set of biochemical possibilities. To address this question, we construct "in silico" all possible biochemically feasible alternatives to the trunk pathway of glycolysis and gluconeogenesis, one of the most highly conserved pathways in metabolism. We show that, even though a large number of alternative pathways exist, the alternatives carry lower flux than the real pathway under typical physiological conditions. Alternative pathways that could potentially carry higher flux often lead to infeasible intermediate metabolite concentrations. We also find that if physiological conditions were different, different pathways could outperform those found in nature. Our results demonstrate how the rules of biochemistry restrict the alternatives that are open to evolution, and suggest that the existing trunk pathway of glycolysis and gluconeogenesis represents a maximal flux solution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Osmosis with active solutes

Despite much current interest in active matter, little is known about osmosis in active systems. Using molecular dynamics simulations, we investigate how active solutes perturb osmotic steady states. We find that solute activity increases the osmotic pressure, and can also expel solvent from the solution - i.e. cause reverse osmosis. The latter effect cannot be described by an effective temperature, but can be reproduced by mapping the active solution onto a passive one with the same degree of local structuring as the passive solvent component. Our results provide a basic framework for understanding active osmosis, and suggest that activity-induced structuring of the passive component may play a key role in the physics of active-passive mixtures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase behaviour of active Brownian particles: The role of dimensionality

Recently, there has been much interest in activity-induced phase separations in concentrated suspensions of "active Brownian particles" (ABPs), self-propelled spherical particles whose direction of motion relaxes through thermal rotational diffusion. To date, almost all these studies have been restricted to 2 dimensions. In this work we study activity-induced phase separation in 3D and compare the results with previous and new 2D simulations. To this end, we performed state-of-the-art Brownian dynamics simulations of up to 40 million ABPs -- such very large system sizes are unavoidable to evade finite size effects in 3D. Our results confirm the picture established for 2D systems in which an activity-induced phase separation occurs, with strong analogies to equilibrium gas-liquid spinodal decomposition, in spite of the purely non-equilibrium nature of the driving force behind the phase separation. However, we also find important differences between the 2D and 3D cases. Firstly, the shape and position of the phase boundaries is markedly different for the two cases. Secondly, for the 3D coarsening kinetics we find that the domain size grows in time according to the classical diffusive $t^{1/3}$ law, in contrast to the nonstandard subdiffusive exponent observed in 2D.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scalar $ϕ^4$ field theory for active-particle phase separation

Recent theories predict phase separation among orientationally disordered active particles whose propulsion speed decreases rapidly enough with density. Coarse-grained models of this process show time-reversal symmetry (detailed balance) to be restored for uniform states, but broken by gradient terms; hence detailed-balance violation is strongly coupled to interfacial phenomena. To explore the subtle generic physics resulting from such coupling we here introduce `Active Model B'. This is a scalar $ϕ^4$ field theory (or phase-field model) that minimally violates detailed balance via a leading-order square-gradient term. We find that this additional term has modest effects on coarsening dynamics, but alters the static phase diagram by creating a jump in (thermodynamic) pressure across flat interfaces. Both results are surprising, since interfacial phenomena are always strongly implicated in coarsening dynamics but are, in detailed-balance systems, irrelevant for phase equilibria.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The role of mechanical forces in the planar-to-bulk transition in growing Escherichia coli microcolonies

Mechanical forces are obviously important in the assembly of three-dimensional multicellular structures, but their detailed role is often unclear. We have used growing microcolonies of the bacterium \emph{Escherichia coli} to investigate the role of mechanical forces in the transition from two-dimensional growth (on the interface between a hard surface and a soft agarose pad) to three-dimensional growth (invasion of the agarose). We measure the position within the colony where the invasion transition happens, the cell density within the colony, and the colony size at the transition as functions of the concentration of the agarose. We use a phenomenological theory, combined with individual-based computer simulations, to show how mechanical forces acting between the bacterial cells, and between the bacteria and the surrounding matrix, lead to the complex phenomena observed in our experiments - in particular a non-trivial dependence of the colony size at the transition on the agarose concentration. Matching these approaches leads to a prediction for how the friction coefficient between the bacteria and the agarose should vary with agarose concentration. Our experimental conditions mimic numerous clinical and environmental scenarios in which bacteria invade soft matrices, as well as shedding more general light on the transition between two- and three-dimensional growth in multicellular assemblies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A continuum theory of phase separation kinetics for active Brownian particles

Active Brownian particles (ABPs), when subject to purely repulsive interactions, are known to undergo activity-induced phase separation broadly resembling an equilibrium (attraction-induced) gas-liquid coexistence. Here we present an accurate continuum theory for the dynamics of phase-separating ABPs, derived by direct coarse-graining, capturing leading-order density gradient terms alongside an effective bulk free energy. Such gradient terms do not obey detailed balance; yet we find coarsening dynamics closely resembling that of equilibrium phase separation. Our continuum theory is numerically compared to large-scale direct simulations of ABPs and accurately accounts for domain growth kinetics, domain topologies and coexistence densities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Automatic, optimized interface placement in forward flux sampling simulations

Forward flux sampling (FFS) provides a convenient and efficient way to simulate rare events in equilibrium or non-equilibrium systems. FFS ratchets the system from an initial state to a final state via a series of interfaces in phase space. The efficiency of FFS depends sensitively on the positions of the interfaces. We present two alternative methods for placing interfaces automatically and adaptively in their optimal locations, on-the-fly as an FFS simulation progresses, without prior knowledge or user intervention. These methods allow the FFS simulation to advance efficiently through bottlenecks in phase space by placing more interfaces where the probability of advancement is lower. The methods are demonstrated both for a single-particle test problem and for the crystallization of Yukawa particles. By removing the need for manual interface placement, our methods both facilitate the setting up of FFS simulations and improve their performance, especially for rare events which involve complex trajectories through phase space, with many bottlenecks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Malliavin weight sampling: a practical guide

Malliavin weight sampling (MWS) is a stochastic calculus technique for computing the derivatives of averaged system properties with respect to parameters in stochastic simulations, without perturbing the system's dynamics. It applies to systems in or out of equilibrium, in steady state or time-dependent situations, and has applications in the calculation of response coefficients, parameter sensitivities and Jacobian matrices for gradient-based parameter optimisation algorithms. The implementation of MWS has been described in the specific contexts of kinetic Monte Carlo and Brownian dynamics simulation algorithms. Here, we present a general theoretical framework for deriving the appropriate MWS update rule for any stochastic simulation algorithm. We also provide pedagogical information on its practical implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Osmosis in a minimal model system

Osmosis plays a central role in the function of living and soft matter systems. While the thermodynamics of osmosis is well understood, the underlying microscopic dynamical mechanisms remain the subject of discussion. Unraveling these mechanisms is a crucial prerequisite for eventually understanding osmosis in non-equilibrium systems. Here, we investigate the microscopic basis of osmosis, in a system at equilibrium, using molecular dynamics simulations of a minimal model in which repulsive solute and solvent particles differ only in their interactions with an external potential. For this system, we can derive a simple virial-like relation for the osmotic pressure. Our simulations support an intuitive picture in which the solvent concentration gradient, at osmotic equilibrium, arises from the balance between an outward force, caused by the increased total density in the solution, and an inward diffusive flux caused by the decreased solvent density in the solution. While more complex effects may occur in other osmotic systems, they are not required for a description of the basic physics of osmosis in this minimal model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Parasites on parasites: coupled fluctuations in stacked contact processes

We present a model for host-parasite dynamics which incorporates both vertical and horizontal transmission as well as spatial structure. Our model consists of stacked contact processes (CP), where the dynamics of the host is a simple CP on a lattice while the dynamics of the parasite is a secondary CP which sits on top of the host-occupied sites. In the simplest case, where infection does not incur any cost, we uncover a novel effect: a nonmonotonic dependence of parasite prevalence on host turnover. Inspired by natural examples of hyperparasitism, we extend our model to multiple levels of parasites and identify a transition between the maintenance of a finite and infinite number of levels, which we conjecture is connected to a roughening transition in models of surface-growth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Computing sensitivity coefficients in Brownian dynamics simulations by Malliavin weight sampling

We present a method for computing parameter sensitivities and response coefficients in Brownian dynamics simulations. The method involves tracking auxiliary variables (Malliavin weights) in addition to the usual particle positions, in an unperturbed simulation. The Malliavin weights sample the derivatives of the probability density with respect to the parameters of interest and are also interesting dynamical objects in themselves. Malliavin weight sampling is simple to implement, applies to equilibrium or nonequilibrium, steady state or time-dependent systems, and scales more efficiently than standard finite difference methods

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Computing the local pressure in molecular dynamics simulations

Computer simulations of inhomogeneous soft matter systems often require accurate methods for computing the local pressure. We present a simple derivation, based on the virial relation, of two equivalent expressions for the local (atomistic) pressure in a molecular dynamics simulation. One of these expressions, previously derived by other authors via a different route, involves summation over interactions between particles within the region of interest; the other involves summation over interactions across the boundary of the region of interest. We illustrate our derivation using simulations of a simple osmotic system; both expressions produce accurate results even when the region of interest over which the pressure is measured is very small.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mixed population of competing TASEPs with a shared reservoir of particles

We introduce a mean-field theoretical framework to describe multiple totally asymmetric simple exclusion processes (TASEPs) with different lattice lengths, entry and exit rates, competing for a finite reservoir of particles. We present relations for the partitioning of particles between the reservoir and the lattices: these relations allow us to show that competition for particles can have non-trivial effects on the phase behavior of individual lattices. For a system with non-identical lattices, we find that when a subset of lattices undergoes a phase transition from low to high density, the entire set of lattice currents becomes independent of total particle number. We generalize our approach to systems with a continuous distribution of lattice parameters, for which we demonstrate that measurements of the current carried by a single lattice type can be used to extract the entire distribution of lattice parameters. Our approach applies to populations of TASEPs with any distribution of lattice parameters, and could easily be extended beyond the mean-field case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mutational pathway determines whether drug gradients accelerate evolution of drug-resistant cells

Drug gradients are believed to play an important role in the evolution of bacteria resistant to antibiotics and tumors resistant to anti-cancer drugs. We use a statistical physics model to study the evolution of a population of malignant cells exposed to drug gradients, where drug resistance emerges via a mutational pathway involving multiple mutations. We show that a non-uniform drug distribution has the potential to accelerate the emergence of resistance when the mutational pathway involves a long sequence of mutants with increasing resistance, but if the pathway is short or crosses a fitness valley, the evolution of resistance may actually be slowed down by drug gradients. These predictions can be verified experimentally, and may help to improve strategies for combatting the emergence of resistance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-Stationary Forward Flux Sampling

We present a new method, Non-Stationary Forward Flux Sampling, that allows efficient simulation of rare events in both stationary and non-stationary stochastic systems. The method uses stochastic branching and pruning to achieve uniform sampling of trajectories in phase space and time, leading to accurate estimates for time-dependent switching propensities and time-dependent phase space probability densities. The method is suitable for equilibrium or non-equilibrium systems, in or out of stationary state, including non-Markovian or externally driven systems. We demonstrate the validity of the technique by applying it to a one-dimensional barrier crossing problem that can be solved exactly, and show its usefulness by applying it to the time-dependent switching of a genetic toggle switch.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Steady-state parameter sensitivity in stochastic modeling via trajectory reweighting

Parameter sensitivity analysis is a powerful tool in the building and analysis of biochemical network models. For stochastic simulations, parameter sensitivity analysis can be computationally expensive, requiring multiple simulations for perturbed values of the parameters. Here, we use trajectory reweighting to derive a method for computing sensitivity coefficients in stochastic simulations without explicitly perturbing the parameter values, avoiding the need for repeated simulations. The method allows the simultaneous computation of multiple sensitivity coefficients. Our approach recovers results originally obtained by application of the Girsanov measure transform in the general theory of stochastic processes [A. Plyasunov and A. P. Arkin, J. Comp. Phys. 221, 724 (2007)]. We build on these results to show how the method can be used to compute steady-state sensitivity coefficients from a single simulation run, and we present various efficiency improvements. For models of biochemical signaling networks the method has a particularly simple implementation. We demonstrate its application to a signaling network showing stochastic focussing and to a bistable genetic switch, and present exact results for models with linear propensity functions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Noise-driven oscillations in microbial population dynamics

Microbial populations in the natural environment are likely to experience growth conditions very different from those of a typical laboratory xperiment. In particular, removal rates of biomass and substrate are unlikely to be balanced under realistic environmental conditions. Here, we consider a single population growing on a substrate under conditions where the removal rates of substrate and biomass are not necessarily equal. For a large population, with deterministic growth dynamics, our model predicts that this system can show transient (damped) oscillations. For a small population, demographic noise causes these oscillations to be sustained indefinitely. These oscillations arise when the dynamics of changes in biomass are faster than the dynamics of the substrate, for example, due to a high microbial death rate and/or low substrate flow rates. We show that the same mechanism can produce sustained stochastic oscillations in a two-species, nutrient-cycling microbial ecosystem. Our results suggest that oscillatory population dynamics may be a common feature of small microbial populations in the natural environment, even in the absence of complex interspecies interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nucleation in a sheared Ising model: effects of external field

Simulations using the Forward Flux Sampling method have shown a nonmonotonic de- pendence of the homogeneous nucleation rate on the shear rate for a sheared two dimensional Ising model [R. J. Allen et al, arXiv cond-mat/0805.3029]. For quasi-equilibrium systems (i.e. in the absence of shear), Classical Nucleation Theory (CNT) predicts the dependence of the critical cluster size and the nucleation rate on the external magnetic field. We investigate the behaviour of the sheared Ising model as a function of the external field. At low exter- nal field strength, the same nonmonotonic behaviour holds and the peak in the nucleation rate is remarkably insensitive to the field strength. This suggests that the same external field-dependence holds for the enhancement of nucleation by shear at low shear rates and the suppression of shear at high shear rates. At high field strength, the nucleation behaviour is qualitatively different. We also analyse the size and shape of the largest cluster in the transition state configurations, as a function of the external field. In the sheared system, the transition state cluster becomes larger and more elongated as the field strength decreases. We compare our results for the sheared system to the predictions of the CNT for the quasi- equilibrium case, and find that the CNT cannot easily be used to describe nucleation in the system under shear.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A dynamical phase transition in a model for evolution with migration

Migration between different habitats is ubiquitous among biological populations. In this Letter, we study a simple quasispecies model for evolution in two different habitats, with different fitness landscapes, coupled through one-way migration. Our model applies to asexual, rapidly evolving organisms such as microbes. Our key finding is a dynamical phase transition at a critical value of the migration rate. The time to reach steady state diverges at this critical migration rate. Above the transition, the population is dominated by immigrants from the primary habitat. Below the transition, the genetic composition of the population is highly non-trivial, with multiple coexisting quasispecies which are not native to either habitat. Using results from localization theory, we show that the critical migration rate may be very small --- demonstrating that evolutionary outcomes can be very sensitive to even a small amount of migration.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Switching and growth for microbial populations in catastrophic responsive environments

Phase variation, or stochastic switching between alternative states of gene expression, is common among microbes, and may be important in coping with changing environments. We use a theoretical model to assess whether such switching is a good strategy for growth in environments with occasional catastrophic events. We find that switching can be advantageous, but only when the environment is responsive to the microbial population. In our model, microbes switch randomly between two phenotypic states, with different growth rates. The environment undergoes sudden "catastrophes", the probability of which depends on the composition of the population. We derive a simple analytical result for the population growth rate. For a responsive environment, two alternative strategies emerge. In the "no switching" strategy, the population maximises its instantaneous growth rate, regardless of catastrophes. In the "switching" strategy, the microbial switching rate is tuned to minimise the environmental response. Which of these strategies is most favourable depends on the parameters of the model. Previous studies have shown that microbial switching can be favourable when the environment changes in an unresponsive fashion between several states. Here, we demonstrate an alternative role for phase variation in allowing microbes to maximise their growth in catastrophic responsive environments.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A multi-purpose modular system for high-resolution microscopy at high hydrostatic pressure

We have developed a modular system for high-resolution microscopy at high hydrostatic pressure. The system consists of a pressurised cell of volume ~100 microlitres, a temperature controlled holder, a ram and a piston. We have made each of these components in several versions which can be interchanged to allow a wide range of applications. Here, we report two pressure cells with pressure ranges 0.1-700MPa and 0.1-100MPa, which can be combined with hollow or solid rams and pistons. Our system is designed to work with fluorescent samples (using a confocal or epifluorescence microscope), but also allows for transmitted light microscopy via the hollow ram and piston. The system allows precise control of pressure and temperature [-20-70C], as well as rapid pressure quenching. We demonstrate its performance and versatility with two applications: time-resolved imaging of colloidal phase transitions caused by pressure changes between 0.1MPa and 101MPa, and imaging the growth of Escherichia coli bacteria at 50MPa. We also show that the isotropic-nematic phase transition of pentyl-cyanobiphenyl (5CB) liquid crystal provides a simple, convenient and accurate method for calibrating pressure in the range 0.1-200MPa.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

DNA looping provides stability and robustness to the bacteriophage lambda switch

The bistable gene regulatory switch controlling the transition from lysogeny to lysis in bacteriophage lambda presents a unique challenge to quantitative modeling. Despite extensive characterization of this regulatory network, the origin of the extreme stability of the lysogenic state remains unclear. We have constructed a stochastic model for this switch. Using Forward Flux Sampling simulations, we show that this model predicts an extremely low rate of spontaneous prophage induction in a recA mutant, in agreement with experimental observations. In our model, the DNA loop formed by octamerization of CI bound to the O_L and O_R operator regions is crucial for stability, allowing the lysogenic state to remain stable even when a large fraction of the total CI is depleted by nonspecific binding to genomic DNA. DNA looping also ensures that the switch is robust to mutations in the order of the O_R binding sites. Our results suggest that DNA looping can provide a mechanism to maintain a stable lysogenic state in the face of a range of challenges including noisy gene expression, nonspecific DNA binding and operator site mutations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Forward Flux Sampling for rare event simulations

Rare events are ubiquitous in many different fields, yet they are notoriously difficult to simulate because few, if any, events are observed in a conventiona l simulation run. Over the past several decades, specialised simulation methods have been developed to overcome this problem. We review one recently-developed class of such methods, known as Forward Flux Sampling. Forward Flux Sampling uses a series of interfaces between the initial and final states to calculate rate constants and generate transition paths, for rare events in equilibrium or nonequilibrium systems with stochastic dynamics. This review draws together a number of recent advances, summarizes several applications of the method and highlights challenges that remain to be overcome.

0 citations0 reviewsNo code linkedOpen access