Source author record

Michael Assaf

Michael Assaf appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

30 published item(s)

preprint2022arXiv

Classification of anomalous diffusion in animal movement data using power spectral analysis

The field of movement ecology has seen a rapid increase in high-resolution data in recent years, leading to the development of numerous statistical and numerical methods to analyse relocation trajectories. Data are often collected at the level of the individual and for long periods that may encompass a range of behaviours. Here, we use the power spectral density (PSD) to characterise the random movement patterns of a black-winged kite (Elanus caeruleus) and a white stork (Ciconia ciconia). The tracks are first segmented and clustered into different behaviours (movement modes), and for each mode we measure the PSD and the ageing properties of the process. For the foraging kite we find $1/f$ noise, previously reported in ecological systems mainly in the context of population dynamics, but not for movement data. We further suggest plausible models for each of the behavioural modes by comparing both the measured PSD exponents and the distribution of the single-trajectory PSD to known theoretical results and simulations.

preprint2022arXiv

Ergodicity breaking in area-restricted search of avian predators

Quantifying and comparing patterns of dynamical ecological systems require averaging over measurable quantities. For example, to infer variation in movement and behavior, metrics such as step length and velocity are averaged over large ensembles. Yet, in nonergodic systems such averaging is inconsistent; thus, identifying ergodicity breaking is essential in ecology. Using rich high-resolution movement datasets ($>\! 7 \times 10^7$ localizations) from 70 individuals and continuous-time random walk modeling, we find subdiffusive behavior and ergodicity breaking in the localized movement of three species of avian predators. Small-scale, within-patch movement was found to be qualitatively different, not inferrable and separated from large-scale inter-patch movement. Local search is characterized by long power-law-distributed waiting times with diverging mean, giving rise to ergodicity breaking in the form of considerable variability uniquely observed at this scale. This implies that wild animal movement is scale specific with no typical waiting time at the local scale.

preprint2022arXiv

Extreme outbreak dynamics in epidemic models

Motivated by recent epidemic outbreaks, including those of COVID-19, we solve the canonical problem of calculating the dynamics and likelihood of extensive outbreaks in a population within a large class of stochastic epidemic models with demographic noise, including the Susceptible-Infected-Recovered (SIR) model and its general extensions. In the limit of large populations, we compute the probability distribution for all extensive outbreaks, including those that entail unusually large or small (extreme) proportions of the population infected. Our approach reveals that, unlike other well-known examples of rare events occurring in discrete-state stochastic systems, the statistics of extreme outbreaks emanate from a full continuum of Hamiltonian paths, each satisfying unique boundary conditions with a conserved probability flux.

preprint2022arXiv

Fluctuations and first-passage properties of systems of Brownian particles with reset

We study stationary fluctuations in two models involving $N$ Brownian particles undergoing stochastic resetting to the origin in 1d. We start with the basic reset model where the particles reset independently (model A). Then we introduce nonlocal interparticle correlations by postulating that only the particle farthest from the origin is reset (model B). At long times both models approach nonequilibrium steady states. In the limit of $N\to \infty$, the steady-state particle density in model A has an infinite support, whereas in model B it has a compact support. A finite system radius, which scales at large $N$ as $\ln N$, appears in model A when $N$ is finite. In both models we study stationary fluctuations of the center of mass of the system and of the system's radius due to the random character of the Brownian motion and of the resetting events. In model A we determine exact distributions of these two quantities. The variance of the center of mass for both models scales as $1/N$. The variance of the radius is independent of $N$ in model A and exhibits an unusual scaling $(\ln N)/N$ in model B. The latter scaling is intimately related to the $1/f$ noise in the radius autocorrelations. Finally, we evaluate the mean first-passage time (MFPT) to a distant target in model A, model B, and the BBM. For model A we obtain an exact asymptotic expression for the MFPT which scales as $1/N$. For model B, and for the "Brownian bees" model, we propose a sharp upper bound for the MFPT. The bound assumes an ``evaporation" scenario, where the first passage requires multiple attempts of a single particle, which breaks away from the rest of the particles, to reach the target. The resulting MFPT for model B and the Brownian bees model scales exponentially with $\sqrt{N}$. We verify this bound by performing highly efficient weighted-ensemble simulations of the first passage in model B.

preprint2022arXiv

Unravelling the origins of anomalous diffusion: from molecules to migrating storks

Anomalous diffusion or, more generally, anomalous transport, with nonlinear dependence of the mean-squared displacement on the measurement time, is ubiquitous in nature. It has been observed in processes ranging from microscopic movement of molecules to macroscopic, large-scale paths of migrating birds. Using data from multiple empirical systems, spanning 12 orders of magnitude in length and 8 orders of magnitude in time, we employ a method to detect the individual underlying origins of anomalous diffusion and transport in the data. This method decomposes anomalous transport into three primary effects: long-range correlations ("Joseph effect"), fat-tailed probability density of increments ("Noah effect"), and non-stationarity ("Moses effect"). We show that such a decomposition of real-life data allows to infer nontrivial behavioral predictions, and to resolve open questions in the fields of single particle tracking in living cells and movement ecology.

preprint2021arXiv

Phase transition in non-Markovian animal exploration model with preferential returns

We study a non-Markovian and nonstationary model of animal mobility incorporating both exploration and memory in the form of preferential returns. We derive exact results for the probability of visiting a given number of sites and develop a practical WKB approximation to treat the nonstationary problem. We further show that this model adequately describes empirical movement data of Egyptian fruit bats (Rousettus aegyptiacus) when accounting for inter-individual variation in the population. Finally, we study the probability of visiting any site a given number of times and derive the corresponding mean-field equation. Here, we find a remarkable phase transition occurring at preferential returns which scale linearly with past visits. Following empirical evidence, we suggest that this phase transition reflects a trade-off between extensive and intensive foraging modes.

preprint2020arXiv

Population Dynamics in a Changing Environment: Random versus Periodic Switching

Environmental changes greatly influence the evolution of populations. Here, we study the dynamics of a population of two strains, one growing slightly faster than the other, competing for resources in a time-varying binary environment modeled by a carrying capacity switching either randomly or periodically between states of abundance and scarcity. The population dynamics is characterized by demographic noise (birth and death events) coupled to a varying environment. We elucidate the similarities and differences of the evolution subject to a stochastically- and periodically-varying environment. Importantly, the population size distribution is generally found to be broader under intermediate and fast random switching than under periodic variations, which results in markedly different asymptotic behaviors between the fixation probability of random and periodic switching. We also determine the detailed conditions under which the fixation probability of the slow strain is maximal.

preprint2019arXiv

Extinction risk of a Metapopulation under the Allee Effect

We study the extinction risk of a fragmented population residing on a network of patches coupled by migration, where the local patch dynamics include the Allee effect. We show that mixing between patches dramatically influences the population's viability. Slow migration is shown to always increase the population's global extinction risk compared to the isolated case. At fast migration, we demonstrate that synchrony between patches minimizes the population's extinction risk. Moreover, we discover a critical migration rate that maximizes the extinction risk of the population, and identify an early-warning signal when approaching this state. Our theoretical results are confirmed via the highly-efficient weighted ensemble method. Notably, our analysis can also be applied to studying switching in gene regulatory networks with multiple transcriptional states.

preprint2019arXiv

Population switching under a time-varying environment

We examine the switching dynamics of a stochastic population subjected to a deterministically time-varying environment. Our approach is demonstrated in the realm of ecology on a problem of population establishment. Here, by assuming a constant immigration pressure along with a strong Allee effect, at the deterministic level one obtains a critical population size beyond which the system experiences establishment. Notably the latter has been shown to be strongly influenced by the interplay between demographic and environmental noise. We consider two prototypical examples for environmental variations: a temporary environmental change, and a periodically-varying environment. By employing a semi-classical approximation we compute, within exponential accuracy, the change in the establishment probability and mean establishment time of the population, due to the environmental variability. Our analytical results are verified by using a modified Gillespie algorithm which accounts for explicitly time-dependent reaction rates. Finally, our theoretical approach can also be useful in studying switching dynamics in gene regulatory networks under extrinsic variations.

preprint2016arXiv

Colonization of a territory by a stochastic population under a strong Allee effect and a low immigration pressure

We study the dynamics of colonization of a territory by a stochastic population at low immigration pressure. We assume a sufficiently strong Allee effect that introduces, in deterministic theory, a large critical population size for colonization. At low immigration rates, the average pre-colonization population size is small thus invalidating the WKB approximation to the master equation. We circumvent this difficulty by deriving an exact zero-flux solution of the master equation and matching it with an approximate non-zero-flux solution of the pertinent Fokker-Planck equation in a small region around the critical population size. This procedure provides an accurate evaluation of the quasi-stationary probability distribution of population sizes in the pre-colonization state, and of the mean time to colonization, for a wide range of immigration rates. At sufficiently high immigration rates our results agree with WKB results obtained previously. At low immigration rates the results can be very different.

preprint2016arXiv

Effect of reaction step-size noise on the switching dynamics of stochastic populations

In genetic circuits, when the mRNA lifetime is short compared to the cell cycle, proteins are produced in geometrically-distributed bursts, which greatly affects the cellular switching dynamics between different metastable phenotypic states. Motivated by this scenario, we study a general problem of switching or escape in stochastic populations, where influx of particles occurs in groups or bursts, sampled from an arbitrary distribution. The fact that the step size of the influx reaction is a-priori unknown, and in general, may fluctuate in time with a given correlation time and statistics, introduces an additional non-demographic step-size noise into the system. Employing the probability generating function technique in conjunction with Hamiltonian formulation, we are able to map the problem in the leading order onto solving a stationary Hamilton-Jacobi equation. We show that bursty influx exponentially decreases the mean escape time compared to the "usual case" of single-step influx. In particular, close to bifurcation we find a simple analytical expression for the mean escape time, which solely depends on the mean and variance of the burst-size distribution. Our results are demonstrated on several realistic distributions and compare well with numerical Monte-Carlo simulations.

preprint2016arXiv

Rare events in stochastic populations under bursty reproduction

Recently, a first step was made by the authors towards a systematic investigation of the effect of reaction-step-size noise - uncertainty in the step size of the reaction - on the dynamics of stochastic populations. This was done by investigating the effect of bursty influx on the switching dynamics of stochastic populations. Here we extend this formalism to account for bursty reproduction processes, and improve the accuracy of the formalism to include subleading-order corrections. Bursty reproduction appears in various contexts, where notable examples include bursty viral production from infected cells, and reproduction of mammals involving varying number of offspring. The main question we quantitatively address is how bursty reproduction affects the overall fate of the population. We consider two complementary scenarios: population extinction and population survival; in the former a population gets extinct after maintaining a long-lived metastable state, whereas in the latter a population proliferates despite undergoing a deterministic drift towards extinction. In both models reproduction occurs in bursts, sampled from an arbitrary distribution. In the extinction problem, we show that bursty reproduction broadens the quasi-stationary distribution of population sizes in the metastable state, which results in an exponential decrease of the mean time to extinction. In the survival problem, bursty reproduction yields an exponential increase in survival probability of the population. Close to the bifurcation limit our analytical results simplify considerably and are shown to depend solely on the mean and variance of the burst-size distribution. Our formalism is demonstrated on several realistic distributions which all compare well with numerical Monte-Carlo simulations.

preprint2016arXiv

The dynamics of simple gene-network motifs subject to extrinsic fluctuations

Cellular processes do not follow deterministic rules; even in identical environments genetically identical cells can make random choices leading to different phenotypes. This randomness originates from fluctuations present in the biomolecular interaction networks. Most previous work has been focused on the intrinsic noise (IN) of these networks. Yet, especially for high-copy-number biomolecules, extrinsic or environmental noise (EN) has been experimentally shown to dominate the variation. Here, we develop an analytical formalism that allows for calculation of the effect of EN on gene-expression motifs. We introduce a method for modeling bounded EN as an auxiliary species in the master equation. The method is fully generic and is not limited to systems with small EN magnitudes. We focus our study on motifs that can be viewed as the building blocks of genetic switches: a nonregulated gene, a self-inhibiting gene, and a self-promoting gene. The role of the EN properties (magnitude, correlation time, and distribution) on the statistics of interest are systematically investigated, and the effect of fluctuations in different reaction rates is compared. Due to its analytical nature, our formalism can be used to quantify the effect of EN on the dynamics of biochemical networks and can also be used to improve the interpretation of data from single-cell gene-expression experiments.

preprint2015arXiv

Genetic Toggle Switch in the Absence of Cooperative Binding: Exact Results

We present an analytical treatment of a genetic switch model consisting of two mutually inhibiting genes operating without cooperative binding of the corresponding transcription factors. Previous studies have numerically shown that these systems can exhibit bimodal dynamics without possessing two stable fixed points at the deterministic level. We analytically show that bimodality is induced by the noise and find the critical repression strength that controls a transition between the bimodal and non-bimodal regimes. We also identify characteristic polynomial scaling laws of the mean switching time between bimodal states. These results, independent of the model under study, reveal essential differences between these systems and systems with cooperative binding, where there is no critical threshold for bimodality and the mean switching time scales exponentially with the system size.

preprint2015arXiv

Stochastic dynamics and logistic population growth

The Verhulst model is probably the best known macroscopic rate equation in population ecology. It depends on two parameters, the intrinsic growth rate and the carrying capacity. These parameters can be estimated for different populations and are related to the reproductive fitness and the competition for limited resources, respectively. We investigate analytically and numerically the simplest possible microscopic scenarios that give rise to the logistic equation in the deterministic mean-field limit. We provide a definition of the two parameters of the Verhulst equation in terms of microscopic parameters. In addition, we derive the conditions for extinction or persistence of the population by employing either the "momentum-space" spectral theory or the "real-space" Wentzel-Kramers-Brillouin (WKB) approximation to determine the probability distribution function and the mean time to extinction of the population. Our analytical results agree well with numerical simulations.

preprint2013arXiv

Cooperation dilemma in finite populations under fluctuating environments

We present a novel approach allowing the study of rare events like fixation under fluctuating environments, modeled as extrinsic noise, in evolutionary processes characterized by the dominance of one species. Our treatment consists of mapping the system onto an auxiliary model, exhibiting metastable species coexistence, that can be analyzed semiclassically. This approach enables us to study the interplay between extrinsic and demographic noise on the statistics of interest. We illustrate our theory by considering the paradigmatic prisoner's dilemma game whose evolution is described by the probability that cooperators fixate the population and replace all defectors. We analytically and numerically demonstrate that extrinsic noise may drastically enhance the cooperation fixation probability and even change its functional dependence on the population size. These results, which generalize earlier works in population genetics, indicate that extrinsic noise may help sustain and promote a much higher level of cooperation than static settings.

preprint2013arXiv

Extrinsic noise driven phenotype switching in a self-regulating gene

Due to inherent noise in intracellular networks cellular decisions can be random, so genetically identical cells can display different phenotypic behavior even in identical environments. Most previous work in understanding the decision-making process has focused on the role of intrinsic noise in these systems. Yet, especially in the high copy-number regime, extrinsic noise has been shown to be much more significant. Here, using a prototypical example of a bistable self-regulating gene model, we develop a theoretical framework describing the combined effect of intrinsic and extrinsic noise on the dynamics of stochastic genetic switches. Employing our theory and Monte Carlo simulations, we show that extrinsic noise not only significantly alters the lifetimes of the phenotypic states, but can induce bistability in unexpected regions of parameter space, and may fundamentally change the escape mechanism. These results have implications for interpreting experimentally observed heterogeneity in cellular populations and for stochastic modeling of cellular decision processes.

preprint2012arXiv

Effect of weak rotation on large-scale circulation cessations in turbulent convection

We investigate the effect of weak rotation on the large-scale circulation (LSC) of turbulent Rayleigh-Bénard convection, using the theory for cessations in a low-dimensional stochastic model of the flow previously studied. We determine the cessation frequency of the LSC as a function of rotation, and calculate the statistics of the amplitude and azimuthal velocity fluctuations of the LSC as a function of the rotation rate for different Rayleigh numbers. Furthermore, we show that the tails of the reorientation probability distribution function remain unchanged for rotating systems, while the distribution of the LSC amplitude and correspondingly the cessation frequency are strongly affected by rotation. Our results are in close agreement with experimental observations.

preprint2012arXiv

Metastability and anomalous fixation in evolutionary games on scale-free networks

We study the influence of complex graphs on the metastability and fixation properties of a set of evolutionary processes. In the framework of evolutionary game theory, where the fitness and selection are frequency-dependent and vary with the population composition, we analyze the dynamics of snowdrift games (characterized by a metastable coexistence state) on scale-free networks. Using an effective diffusion theory in the weak selection limit, we demonstrate how the scale-free structure affects the system's metastable state and leads to anomalous fixation. In particular, we analytically and numerically show that the probability and mean time of fixation are characterized by stretched exponential behaviors with exponents depending on the network's degree distribution.

preprint2011arXiv

Determining the stability of genetic switches: explicitly accounting for mRNA noise

Cells use genetic switches to shift between alternate gene expression states, e.g., to adapt to new environments or to follow a developmental pathway. Here, we study the dynamics of switching in a generic-feedback on/off switch. Unlike protein-only models, we explicitly account for stochastic fluctuations of mRNA, which have a dramatic impact on switch dynamics. Employing the WKB theory to treat the underlying chemical master equations, we obtain accurate results for the quasi-stationary distributions of mRNA and protein copy numbers and for the mean switching time, starting from either state. Our analytical results agree well with Monte Carlo simulations. Importantly, one can use the approach to study the effect of varying biological parameters on switch stability.

preprint2011arXiv

Fixation of a Deleterious Allele under Mutation Pressure and Finite Selection Intensity

The mean fixation time of a deleterious mutant allele is studied beyond the diffusion approximation. As in Kimura's classical work [M. Kimura, Proc. Natl. Acad. Sci. U.S.A. Vol.77, 522 (1980)], that was motivated by the problem of fixation in the presence of amorphic or hypermorphic mutations, we consider a diallelic model at a single locus comprising a wild-type A and a mutant allele A' produced irreversibly from A at small uniform rate v. The relative fitnesses of the mutant homozygotes A'A', mutant heterozygotes A'A and wild-type homozygotes AA are 1-s, 1-h and 1, respectively, where it is assumed that v<< s. Here, we adopt an approach based on the direct treatment of the underlying Markov chain (birth-death process) obeyed by the allele frequency (whose dynamics is prescribed by the Moran model), which allows to accurately account for the effects of large fluctuations. After a general description of the theory, we focus on the case of a deleterious mutant allele (i.e. s>0) and discuss three situations: when the mutant is (i) completely dominant (s=h); (ii) completely recessive (h=0), and (iii) semi-dominant (h=s/2). Our theoretical predictions for the mean fixation time and the quasi-stationary distribution of the mutant population in the coexistence state, are shown to be in excellent agreement with numerical simulations. Furthermore, when s is finite, we demonstrate that our results are superior to those of the diffusion theory that is shown to be an accurate approximation only when N_e s^2 << 1, where N_e is the effective population size.

preprint2011arXiv

Rare fluctuations and large-scale circulation cessations in turbulent convection

In turbulent Rayleigh-Benard convection, a large-scale circulation (LSC) develops in a nearly vertical plane, and is maintained by rising and falling plumes detaching from the unstable thermal boundary layers. Rare but large fluctuations in the LSC amplitude can lead to extinction of the LSC (a cessation event), followed by the re-emergence of another LSC with a different (random) azimuthal orientation. We extend previous models of the LSC dynamics to include momentum and thermal diffusion in the azimuthal plane, and calculate the tails of the probability distributions of both the amplitude and azimuthal angle. Our analytical results are in very good agreement with experimental data.

preprint2010arXiv

Fixation in Evolutionary Games under Non-Vanishing Selection

One of the most striking effect of fluctuations in evolutionary game theory is the possibility for mutants to fixate (take over) an entire population. Here, we generalize a recent WKB-based theory to study fixation in evolutionary games under non-vanishing selection, and investigate the relation between selection intensity w and demographic (random) fluctuations. This allows the accurate treatment of large fluctuations and yields the probability and mean times of fixation beyond the weak selection limit. The power of the theory is demonstrated on prototypical models of cooperation dilemmas with multiple absorbing states. Our predictions compare excellently with numerical simulations and, for finite w, significantly improve over those of the Fokker-Planck approximation.

preprint2010arXiv

Large Fluctuations and Fixation in Evolutionary Games

We study large fluctuations in evolutionary games belonging to the coordination and anti-coordination classes. The dynamics of these games, modeling cooperation dilemmas, is characterized by a coexistence fixed point separating two absorbing states. We are particularly interested in the problem of fixation that refers to the possibility that a few mutants take over the entire population. Here, the fixation phenomenon is induced by large fluctuations and is investigated by a semi-classical WKB (Wentzel-Kramers-Brillouin) theory generalized to treat stochastic systems possessing multiple absorbing states. Importantly, this method allows us to analyze the combined influence of selection and random fluctuations on the evolutionary dynamics \textit{beyond} the weak selection limit often considered in previous works. We accurately compute, including pre-exponential factors, the probability distribution function in the long-lived coexistence state and the mean fixation time necessary for a few mutants to take over the entire population in anti-coordination games, and also the fixation probability in the coordination class. Our analytical results compare excellently with extensive numerical simulations. Furthermore, we demonstrate that our treatment is superior to the Fokker-Planck approximation when the selection intensity is finite.

preprint2010arXiv

Large fluctuations in stochastic population dynamics: momentum space calculations

Momentum-space representation renders an interesting perspective to theory of large fluctuations in populations undergoing Markovian stochastic gain-loss processes. This representation is obtained when the master equation for the probability distribution of the population size is transformed into an evolution equation for the probability generating function. Spectral decomposition then brings about an eigenvalue problem for a non-Hermitian linear differential operator. The ground-state eigenmode encodes the stationary distribution of the population size. For long-lived metastable populations which exhibit extinction or escape to another metastable state, the quasi-stationary distribution and the mean time to extinction or escape are encoded by the eigenmode and eigenvalue of the lowest excited state. If the average population size in the stationary or quasi-stationary state is large, the corresponding eigenvalue problem can be solved via WKB approximation amended by other asymptotic methods. We illustrate these ideas in several model examples.

preprint2009arXiv

Extinction of metastable stochastic populations

We investigate extinction of a long-lived self-regulating stochastic population, caused by intrinsic (demographic) noise. Extinction typically occurs via one of two scenarios depending on whether the absorbing state n=0 is a repelling (scenario A) or attracting (scenario B) point of the deterministic rate equation. In scenario A the metastable stochastic population resides in the vicinity of an attracting fixed point next to the repelling point n=0. In scenario B there is an intermediate repelling point n=n_1 between the attracting point n=0 and another attracting point n=n_2 in the vicinity of which the metastable population resides. The crux of the theory is WKB method which assumes that the typical population size in the metastable state is large. Starting from the master equation, we calculate the quasi-stationary probability distribution of the population sizes and the (exponentially long) mean time to extinction for each of the two scenarios. When necessary, the WKB approximation is complemented (i) by a recursive solution of the quasi-stationary master equation at small n and (ii) by the van Kampen system-size expansion, valid near the fixed points of the deterministic rate equation. The theory yields both entropic barriers to extinction and pre-exponential factors, and holds for a general set of multi-step processes when detailed balance is broken. The results simplify considerably for single-step processes and near the characteristic bifurcations of scenarios A and B.

preprint2008arXiv

On population extinction risk in the aftermath of a catastrophic event

We investigate how a catastrophic event (modeled as a temporary fall of the reproduction rate) increases the extinction probability of an isolated self-regulated stochastic population. Using a variant of the Verhulst logistic model as an example, we combine the probability generating function technique with an eikonal approximation to evaluate the exponentially large increase in the extinction probability caused by the catastrophe. This quantity is given by the eikonal action computed over "the optimal path" (instanton) of an effective classical Hamiltonian system with a time-dependent Hamiltonian. For a general catastrophe the eikonal equations can be solved numerically. For simple models of catastrophic events analytic solutions can be obtained. One such solution becomes quite simple close to the bifurcation point of the Verhulst model. The eikonal results for the increase in the extinction probability caused by a catastrophe agree well with numerical solutions of the master equation.

preprint2008arXiv

Population extinction in a time-modulated environment

The extinction time of an isolated population can be exponentially reduced by a periodic modulation of its environment. We investigate this effect using, as an example, a stochastic branching-annihilation process with a time-dependent branching rate. The population extinction is treated in eikonal approximation, where it is described as an instanton trajectory of a proper reaction Hamiltonian. The modulation of the environment perturbs this trajectory and synchronizes it with the modulation phase. We calculate the corresponding change in the action along the instanton using perturbation techniques supported by numerical calculations. The techniques include a first-order theory with respect to the modulation amplitude, a second-order theory in the spirit of the Kapitsa pendulum effect, and adiabatic theory valid for low modulation frequencies.

preprint2007arXiv

Noise enhanced persistence in a biochemical regulatory network with feedback control

We find that discrete noise of inhibiting (signal) molecules can greatly delay the extinction of plasmids in a plasmid replication system: a prototypical biochemical regulatory network. We calculate the probability distribution of the metastable state of the plasmids and show on this example that the reaction rate equations may fail in predicting the average number of regulated molecules even when this number is large, and the time is much shorter than the mean extinction time.

preprint2007arXiv

Spectral theory of metastability and extinction in a branching-annihilation reaction

We apply the spectral method, recently developed by the authors, to calculate the statistics of a reaction-limited multi-step birth-death process, or chemical reaction, that includes as elementary steps branching A->2A and annihilation 2A->0. The spectral method employs the generating function technique in conjunction with the Sturm-Liouville theory of linear differential operators. We focus on the limit when the branching rate is much higher than the annihilation rate, and obtain accurate analytical results for the complete probability distribution (including large deviations) of the metastable long-lived state, and for the extinction time statistics. The analytical results are in very good agreement with numerical calculations. Furthermore, we use this example to settle the issue of the "lacking" boundary condition in the spectral formulation.