Rare transitions to thin-layer turbulent condensates

Turbulent flows in a thin layer can develop an inverse energy cascade leading to spectral condensation of energy when the layer height is smaller than a certain threshold. These spectral condensates take the form of large-scale vortices in physical space. Recently, evidence for bistability was found in this system close to the critical height: depending on the initial conditions, the flow is either in a condensate state with most of the energy in the two-dimensional (2-D) large-scale modes, or in a three-dimensional (3-D) flow state with most of the energy in the small-scale modes. This bistable regime is characterised by the statistical properties of random and rare transitions between these two locally stable states. Here, we examine these statistical properties in thin-layer turbulent flows, where the energy is injected by either stochastic or deterministic forcing. To this end, by using a large number of direct numerical simulations (DNS), we measure the decay time $\tau_d$ of the 2-D condensate to 3-D flow state and the build-up time $\tau_b$ of the 2-D condensate. We show that both of these times $\tau_d,\tau_b$ follow an exponential distribution with mean values increasing faster than exponentially as the layer height approaches the threshold. We further show that the dynamics of large-scale kinetic energy may be modeled by a stochastic Langevin equation. From time-series analysis of DNS data, we determine the effective potential that shows two minima corresponding to the 2-D and 3-D states when the layer height is close to the threshold

A multiple replica approach to simulate reactive trajectories

A method to generate reactive trajectories, namely equilibrium trajectories leaving a metastable state and ending in another one is proposed. The algorithm is based on simulating in parallel many copies of the system, and selecting the replicas which have reached the highest values along a chosen one-dimensional reaction coordinate. This reaction coordinate does not need to precisely describe all the metastabilities of the system for the method to give reliable results. An extension of the algorithm to compute transition times from one metastable state to another one is also presented. We demonstrate the interest of the method on two simple cases: a one-dimensional two-well potential and a two-dimensional potential exhibiting two channels to pass from one metastable state to another one

Hybrid Pathwise Sensitivity Methods for Discrete Stochastic Models of Chemical Reaction Systems

Stochastic models are often used to help understand the behavior of intracellular biochemical processes. The most common such models are continuous time Markov chains (CTMCs). Parametric sensitivities, which are derivatives of expectations of model output quantities with respect to model parameters, are useful in this setting for a variety of applications. In this paper, we introduce a class of hybrid pathwise differentiation methods for the numerical estimation of parametric sensitivities. The new hybrid methods combine elements from the three main classes of procedures for sensitivity estimation, and have a number of desirable qualities. First, the new methods are unbiased for a broad class of problems. Second, the methods are applicable to nearly any physically relevant biochemical CTMC model. Third, and as we demonstrate on several numerical examples, the new methods are quite efficient, particularly if one wishes to estimate the full gradient of parametric sensitivities. The methods are rather intuitive and utilize the multilevel Monte Carlo philosophy of splitting an expectation into separate parts and handling each in an efficient manner.Comment: 30 pages. The numerical example section has been extensively rewritte

Kinetic approaches to lactose operon induction and bimodality

The quasi-equilibrium approximation is acceptable when molecular interactions are fast enough compared to circuit dynamics, but is no longer allowed when cellular activities are governed by rare events. A typical example is the lactose operon (lac), one of the most famous paradigms of transcription regulation, for which several theories still coexist to describe its behaviors. The lac system is generally analyzed by using equilibrium constants, contradicting single-event hypotheses long suggested by Novick and Weiner (1957). Enzyme induction as an all-or-none phenomenon. Proc. Natl. Acad. Sci. USA 43, 553-566) and recently refined in the study of (Choi et al., 2008. A stochastic single-molecule event triggers phenotype switching of a bacterial cell. Science 322, 442-446). In the present report, a lac repressor (LacI)-mediated DNA immunoprecipitation experiment reveals that the natural LacI-lac DNA complex built in vivo is extremely tight and long-lived compared to the time scale of lac expression dynamics, which could functionally disconnect the abortive expression bursts and forbid using the standard modes of lac bistability. As alternatives, purely kinetic mechanisms are examined for their capacity to restrict induction through: (i) widely scattered derepression related to the arrival time variance of a predominantly backward asymmetric random walk and (ii) an induction threshold arising in a single window of derepression without recourse to nonlinear multimeric binding and Hill functions. Considering the complete disengagement of the lac repressor from the lac promoter as the probabilistic consequence of a transient stepwise mechanism, is sufficient to explain the sigmoidal lac responses as functions of time and of inducer concentration. This sigmoidal shape can be misleadingly interpreted as a phenomenon of equilibrium cooperativity classically used to explain bistability, but which has been reported to be weak in this system

Weighted ensemble: Recent mathematical developments

The weighted ensemble (WE) method, an enhanced sampling approach based on periodically replicating and pruning trajectories in a set of parallel simulations, has grown increasingly popular for computational biochemistry problems, due in part to improved hardware and the availability of modern software. Algorithmic and analytical improvements have also played an important role, and progress has accelerated in recent years. Here, we discuss and elaborate on the WE method from a mathematical perspective, highlighting recent results which have begun to yield greater computational efficiency. Notable among these innovations are variance reduction approaches that optimize trajectory management for systems of arbitrary dimensionality.Comment: 12 pages, 10 figure

Monte Carlo Algorithm for Simulating Reversible Aggregation of Multisite Particles

We present an efficient and exact Monte Carlo algorithm to simulate reversible aggregation of particles with dedicated binding sites. This method introduces a novel data structure of dynamic bond tree to record clusters and sequences of bond formations. The algorithm achieves a constant time cost for processing cluster association and a cost between $\mathcal{O}(\log M)$ and $\mathcal{O}(M)$ for processing bond dissociation in clusters with $M$ bonds. The algorithm is statistically exact and can reproduce results obtained by the standard method. We applied the method to simulate a trivalent ligand and a bivalent receptor clustering system and obtained an average scaling of $\mathcal{O}(M^{0.45})$ for processing bond dissociation in acyclic aggregation, compared to a linear scaling with the cluster size in standard methods. The algorithm also demands substantially less memory than the conventional method.Comment: 8 pages, 3 figure