Source author record

Baruch Meerson

Baruch Meerson 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

60works
18topics
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

60 published item(s)

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.

preprint2021arXiv

Statistics of first-passage Brownian functionals

We study the distribution of first-passage functionals ${\cal A}= \int_0^{t_f} x^n(t)\, dt$, where $x(t)$ is a Brownian motion (with or without drift) with diffusion constant $D$, starting at $x_0>0$, and $t_f$ is the first-passage time to the origin. In the driftless case, we compute exactly, for all $n>-2$, the probability density $P_n(A|x_0)=\text{Prob}.(\mathcal{A}=A)$. This probability density has an essential singular tail as $A\to 0$ and a power-law tail $\sim A^{-(n+3)/(n+2)}$ as $A\to \infty$. The former is reproduced by the optimal fluctuation method (OFM), which also predicts the optimal paths of the conditioned process for small $A$. For the case with a drift toward the origin, where no exact solution is known for general $n>-1$, the OFM predicts the distribution tails. For $A\to 0$ it predicts the same essential singular tail as in the driftless case. For $A\to \infty$ it predicts a stretched exponential tail $-\ln P_n(A|x_0)\sim A^{1/(n+1)}$ for all $n>0$. In the limit of large Péclet number $\text{Pe}= μx_0/(2D)\gg 1$, where $μ$ is the drift velocity, the OFM predicts a large-deviation scaling for all $A$: $-\ln P_n(A|x_0)\simeq\text{Pe}\, Φ_n\left(z= A/\bar{A}\right)$, where $\bar{A}=x_0^{n+1}/{μ(n+1)}$ is the mean value of $\mathcal{A}$. We compute the rate function $Φ_n(z)$ analytically for all $n>-1$. For $n>0$ $Φ_n(z)$ is analytic for all $z$, but for $-1<n<0$ it is non-analytic at $z=1$, implying a dynamical phase transition. The order of this transition is $2$ for $-1/2<n<0$, while for $-1<n<-1/2$ the order of transition changes continuously with $n$. Finally, we apply the OFM to the case of $μ<0$ (drift away from the origin). We show that, when the process is conditioned on reaching the origin, the distribution of $\mathcal{A}$ coincides with the distribution of $\mathcal{A}$ for $μ>0$ with the same $|μ|$.

preprint2020arXiv

The Airy distribution: experiment, large deviations and additional statistics

The Airy distribution (AD) describes the probability distribution of the area under a Brownian excursion. The AD is prominent in several areas of physics, mathematics and computer science. Here we use a dilute colloidal system to directly measure, for the first time, the AD in experiment. We also show how two different techniques of theory of large deviations - the Donsker-Varadhan formalism and the optimal fluctuation method - manifest themselves in the AD. We advance the theory of the AD by calculating, at large and small areas, the position distribution of a Brownian excursion conditioned on a given area, and measure its mean in the experiment. For large areas, we uncover two singularities in the large deviation function, which can be interpreted as dynamical phase transitions of third order. For small areas the position distribution coincides with the Ferrari-Spohn distribution, and we identify the reason for this coincidence.

preprint2020arXiv

Toward the full short-time statistics of an active Brownian particle on the plane

We study the position distribution of a single active Brownian particle (ABP) on the plane. We show that this distribution has a compact support, the boundary of which is an expanding circle. We focus on a short-time regime and employ the optimal fluctuation method (OFM) to study large deviations of the particle position coordinates $x$ and $y$. We determine the optimal paths of the ABP, conditioned on reaching specified values of $x$ and $y$, and the large deviation functions of the marginal distributions of $x$, and of $y$. These marginal distributions match continuously with "near tails" of the $x$ and $y$ distributions of typical fluctuations, studied earlier. We also calculate the large deviation function of the joint $x$ and $y$ distribution $P(x,y,t)$ in a vicinity of a special "zero-noise" point, and show that $\ln P(x,y,t)$ has a nontrivial self-similar structure as a function of $x$, $y$ and $t$. The joint distribution vanishes extremely fast at the expanding circle, exhibiting an essential singularity there. This singularity is inherited by the marginal $x$- and $y$-distributions. We argue that this fingerprint of the short-time dynamics remains there at all times.

preprint2020arXiv

Velocity fluctuations of stochastic reaction fronts propagating into an unstable state: strongly pushed fronts

The empirical velocity of a reaction-diffusion front, propagating into an unstable state, fluctuates because of the shot noises of the reactions and diffusion. Under certain conditions these fluctuations can be described as a diffusion process in the reference frame moving with the average velocity of the front. Here we address pushed fronts, where the front velocity in the deterministic limit is affected by higher-order reactions and is therefore larger than the linear spread velocity. For a subclass of these fronts -- strongly pushed fronts -- the effective diffusion constant $D_f\sim 1/N$ of the front can be calculated, in the leading order, via a perturbation theory in $1/N \ll 1$, where $N\gg 1$ is the typical number of particles in the transition region. This perturbation theory, however, overestimates the contribution of a few fast particles in the leading edge of the front. We suggest a more consistent calculation by introducing a spatial integration cutoff at a distance beyond which the average number of particles is of order 1. This leads to a non-perturbative correction to $D_f$ which even becomes dominant close to the transition point between the strongly and weakly pushed fronts. At the transition point we obtain a logarithmic correction to the $1/N$ scaling of $D_f$. We also uncover another, and quite surprising, effect of the fast particles in the leading edge of the front. Because of these particles, the position fluctuations of the front can be described as a diffusion process only on very long time intervals with a duration $Δt \gg τ_N$, where $τ_N$ scales as $N$. At intermediate times the position fluctuations of the front are anomalously large and non-diffusive. Our extensive Monte-Carlo simulations of a particular reacting lattice gas model support these conclusions.

preprint2019arXiv

Optimal paths of non-equilibrium stochastic fields: the Kardar-Parisi-Zhang interface as a test case

Atypically large fluctuations in macroscopic non-equilibrium systems continue to attract interest. Their probability can often be determined by the optimal fluctuation method (OFM). The OFM brings about a conditional variational problem, the solution of which describes the "optimal path" of the system which dominates the contribution of different stochastic paths to the desired statistics. The OFM proved efficient in evaluating the probabilities of rare events in a host of systems. However, theoretically predicted optimal paths were observed in stochastic simulations only in diffusive lattice gases, where the predicted optimal density patterns are either stationary, or travel with constant speed. Here we focus on the one-point height distribution of the paradigmatic Kardar-Parisi-Zhang interface. Here the optimal paths, corresponding to the distribution tails at short times, are intrinsically non-stationary and can be predicted analytically. Using the mapping to the directed polymer in a random potential at high temperature, we obtain "snapshots" of the optimal paths in Monte-Carlo simulations which probe the tails with an importance sampling algorithm. For each tail we observe a very narrow "tube" of height profiles around a single optimal path which agrees with the analytical prediction. The agreement holds even at long times, supporting earlier assertions of the validity of the OFM in the tails well beyond the weak-noise limit.

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

Dynamical phase transition in large deviation statistics of the Kardar-Parisi-Zhang equation

We study the short-time behavior of the probability distribution $\mathcal{P}(H,t)$ of the surface height $h(x=0,t)=H$ in the Kardar-Parisi-Zhang (KPZ) equation in $1+1$ dimension. The process starts from a stationary interface: $h(x,t=0)$ is given by a realization of two-sided Brownian motion constrained by $h(0,0)=0$. We find a singularity of the large deviation function of $H$ at a critical value $H=H_c$. The singularity has the character of a second-order phase transition. It reflects spontaneous breaking of the reflection symmetry $x \leftrightarrow -x$ of optimal paths $h(x,t)$ predicted by the weak-noise theory of the KPZ equation. At $|H|\gg |H_c|$ the corresponding tail of $\mathcal{P}(H)$ scales as $-\ln \mathcal{P} \sim |H|^{3/2}/t^{1/2}$ and agrees, at any $t>0$, with the proper tail of the Baik-Rains distribution, previously observed only at long times. The other tail of $\mathcal{P}$ scales as $-\ln \mathcal{P} \sim |H|^{5/2}/t^{1/2}$ and coincides with the corresponding tail for the sharp-wedge initial condition.

preprint2016arXiv

Large Deviations of Surface Height in the Kardar-Parisi-Zhang Equation

Using the weak-noise theory, we evaluate the probability distribution $\mathcal{P}(H,t)$ of large deviations of height $H$ of the evolving surface height $h(x,t)$ in the Kardar-Parisi-Zhang (KPZ) equation in one dimension when starting from a flat interface. We also determine the optimal history of the interface, conditioned on reaching the height $H$ at time $t$. We argue that the tails of $\mathcal{P}$ behave, at arbitrary time $t>0$, and in a proper moving frame, as $-\ln \mathcal{P}\sim |H|^{5/2}$ and $\sim |H|^{3/2}$. The $3/2$ tail coincides with the asymptotic of the Gaussian orthogonal ensemble Tracy-Widom distribution, previously observed at long times.

preprint2016arXiv

Large-displacement statistics of the rightmost particle of the one-dimensional branching Brownian motion

Consider a one-dimensional branching Brownian motion, and rescale the coordinate and time so that the rates of branching and diffusion are both equal to $1$. If $X_1(t)$ is the position of the rightmost particle of the branching Brownian motion at time $t$, the empirical velocity $c$ of this rightmost particle is defined as $c=X_1(t)/t$. Using the Fisher-KPP equation, we evaluate the probability distribution ${\mathcal P(c,t)}$ of this empirical velocity $c$ in the long time $t$ limit for $c > 2$. It was already known that, for a single seed particle, ${\mathcal P(c,t)} \sim \exp \,[-(c^2/4-1)t]$ up to a prefactor that can depend on $c$ and $t$. Here we show how to determine this prefactor. The result can be easily generalized to the case of multiple seed particles and to branching random walks associated to other traveling wave equations.

preprint2016arXiv

Short-time height distribution in 1d KPZ equation: starting from a parabola

We study the probability distribution $\mathcal{P}(H,t,L)$ of the surface height $h(x=0,t)=H$ in the Kardar-Parisi-Zhang (KPZ) equation in $1+1$ dimension when starting from a parabolic interface, $h(x,t=0)=x^2/L$. The limits of $L\to\infty$ and $L\to 0$ have been recently solved exactly for any $t>0$. Here we address the early-time behavior of $\mathcal{P}(H,t,L)$ for general $L$. We employ the weak-noise theory - a variant of WKB approximation -- which yields the optimal history of the interface, conditioned on reaching the given height $H$ at the origin at time $t$. We find that at small $H$ $\mathcal{P}(H,t,L)$ is Gaussian, but its tails are non-Gaussian and highly asymmetric. In the leading order and in a proper moving frame, the tails behave as $-\ln \mathcal{P}= f_{+}|H|^{5/2}/t^{1/2}$ and $f_{-}|H|^{3/2}/t^{1/2}$. The factor $f_{+}(L,t)$ monotonically increases as a function of $L$, interpolating between time-independent values at $L=0$ and $L=\infty$ that were previously known. The factor $f_{-}$ is independent of $L$ and $t$, signalling universality of this tail for a whole class of deterministic initial conditions.

preprint2016arXiv

Statistics of large currents in the Kipnis-Marchioro-Presutti model in a ring geometry

We use the macroscopic fluctuation theory to determine the statistics of large currents in the Kipnis-Marchioro-Presutti (KMP) model in a ring geometry. About 10 years ago this simple setting was instrumental in identifying a breakdown of the additivity principle in a class of lattice gases at currents exceeding a critical value. Building on earlier work, we assume that, for supercritical currents, the optimal density profile, conditioned on the given current, has the form of a traveling wave (TW). For the KMP model we find this TW analytically, in terms of elliptic functions, for any supercritical current $I$. Using this TW solution, we evaluate, up to a pre-exponential factor, the probability distribution $P(I)$. We obtain simple asymptotics of the TW and of $P(I)$ for currents close to the critical current, and for currents much larger than the critical current. In the latter case we show that $-\ln P (I) \sim I\ln I$, whereas the optimal density profile acquires a soliton-like shape. Our analytic results are in a very good agreement with Monte-Carlo simulations and numerical solutions of Hurtado and Garrido (2011).

preprint2015arXiv

Extinction of oscillating populations

Established populations often exhibit oscillations in their sizes. If a population is isolated, intrinsic stochasticity of elemental processes can ultimately bring it to extinction. Here we study extinction of oscillating populations in a stochastic version of the Rosenzweig-MacArthur predator-prey model. To this end we extend a WKB approximation (after Wentzel, Kramers and Brillouin) of solving the master equation to the case of extinction from a limit cycle in the space of population sizes. We evaluate the extinction rates and find the most probable paths to extinction by applying Floquet theory to the dynamics of an effective WKB Hamiltonian. We show that the entropic barriers to extinction change in a non-analytic way as the system passes through the Hopf bifurcation. We also study the subleading pre-exponential factors of the WKB approximation.

preprint2015arXiv

First Invader Dynamics in Diffusion-Controlled Absorption

We investigate the average time for the earliest particle to hit a spherical absorber when a homogeneous gas of freely diffusing particles with density $ρ$ and diffusivity $D$ is prepared in a deterministic state and is initially separated by a minimum distance $\ell$ from this absorber. In the high-density limit, this first absorption time scales as $\frac{\ell^2}{D}\frac{1}{\lnρ\ell}$ in one dimension; we also obtain the first absorption time in three dimensions. In one dimension, we determine the probability that the $k^{\rm th}$-closest particle is the first one to hit the absorber. At large $k$, this probability decays as $k^{1/3}\exp(-Ak^{2/3})$, with $A= 1.93299\ldots$ analytically calculable. As a corollary, the characteristic hitting time $T_k$ for the $k^{\rm th}$-closest particle scales as $k^{4/3}$; this corresponds to superdiffusive but still subballistic motion.

preprint2015arXiv

Full absorption statistics of diffusing particles with exclusion

Suppose that an infinite lattice gas of constant density $n_0$, whose dynamics are described by the symmetric simple exclusion process, is brought in contact with a spherical absorber of radius $R$. Employing the macroscopic fluctuation theory and assuming the additivity principle, we evaluate the probability distribution ${\mathcal P}(N)$ that $N$ particles are absorbed during a long time $T$. The limit of $N=0$ corresponds to the survival problem, whereas $N\gg \bar{N}$ describes the opposite extreme. Here $\bar{N}=4πR D_0 n_0 T$ is the \emph{average} number of absorbed particles (in three dimensions), and $D_0$ is the gas diffusivity. For $n_0\ll 1$ the exclusion effects are negligible, and ${\mathcal P}(N)$ can be approximated, for not too large $N$, by the Poisson distribution with mean $\bar{N}$. For finite $n_0$, ${\mathcal P}(N)$ is non-Poissonian. We show that $-\ln{\mathcal P}(N) \simeq n_0 N^2/\bar{N}$ at $N\gg \bar{N}$. At sufficiently large $N$ and $n_0<1/2$ the most likely density profile of the gas, conditional on the absorption of $N$ particles, is non-monotonic in space. We also establish a close connection between this problem and that of statistics of current in finite open systems.

preprint2015arXiv

Macroscopic fluctuation theory and first-passage properties of surface diffusion

We investigate non-equilibrium fluctuations of a solid surface governed by the stochastic Mullins-Herring equation with conserved noise. This equation describes surface diffusion of adatoms accompanied by their exchange between the surface and the bulk of the solid, when desorption of adatoms is negligible. Previous works dealt with dynamic scaling behavior of the fluctuating interface. Here we determine the probability that the interface first reaches a large given height at a specified time. We also find the optimal time history of the interface conditional on this non-equilibrium fluctuation. We obtain these results by developing a macroscopic fluctuation theory of surface diffusion.

preprint2015arXiv

Mortality, Redundancy, and Diversity in Stochastic Search

We investigate a stochastic search process in one dimension under the competing roles of mortality, redundancy, and diversity of the searchers. This picture represents a toy model for the fertilization of an oocyte by sperm. A population of $N$ independent and mortal diffusing searchers all start at $x=L$ and attempt to reach the target at $x=0$. When mortality is irrelevant, the search time scales as $τ_D/\ln N$ for $\ln N\gg 1$, where $τ_D\sim L^2/D$ is the diffusive time scale. Conversely, when the mortality rate $μ$ of the searchers is sufficiently large, the search time scales as $\sqrt{τ_D/μ}$, independent of $N$. When searchers have distinct and high mortalities, a subpopulation with a non-trivial optimal diffusivity are most likely to reach the target. We also discuss the effect of chemotaxis on the search time and its fluctuations.

preprint2015arXiv

Survival of interacting diffusing particles inside a domain with absorbing boundary

Suppose that a $d$-dimensional domain is filled with a gas of (in general, interacting) diffusive particles with density $n_0$. A particle is absorbed whenever it reaches the domain boundary. Employing macroscopic fluctuation theory, we evaluate the probability ${\mathcal P}$ that no particles are absorbed during a long time $T$. We argue that the most likely gas density profile, conditional on this event, is stationary throughout most of the time $T$. As a result, ${\mathcal P}$ decays exponentially with $T$ for a whole class of interacting diffusive gases in any dimension. For $d=1$ the stationary gas density profile and ${\mathcal P}$ can be found analytically. In higher dimensions we focus on the simple symmetric exclusion process (SSEP) and show that $-\ln {\mathcal P}\simeq D_0TL^{d-2} \,s(n_0)$, where $D_0$ is the gas diffusivity, and $L$ is the linear size of the system. We calculate the rescaled action $s(n_0)$ for $d=1$, for rectangular domains in $d=2$, and for spherical domains. Near close packing of the SSEP $s(n_0)$ can be found analytically for domains of any shape and in any dimension.

preprint2015arXiv

The number statistics and optimal history of non-equilibrium steady states of mortal diffusing particles

Suppose that a point-like steady source at $x=0$ injects particles into a half-infinite line. The particles diffuse and die. At long times a non-equilibrium steady state sets in, and we assume that it involves many particles. If the particles are non-interacting, their total number $N$ in the steady state is Poisson-distributed with mean $\bar{N}$ predicted from a deterministic reaction-diffusion equation. Here we determine the most likely density history of this driven system conditional on observing a given $N$. We also consider two prototypical examples of \emph{interacting} diffusing particles: (i) a family of mortal diffusive lattice gases with constant diffusivity (as illustrated by the simple symmetric exclusion process with mortal particles), and (ii) random walkers that can annihilate in pairs. In both examples we calculate the variances of the (non-Poissonian) stationary distributions of $N$.

preprint2014arXiv

Extreme Fluctuations of Current in the Symmetric Simple Exclusion Process: a Non-Stationary Setting

We use the macroscopic fluctuation theory (MFT) to evaluate the probability distribution P of extreme values of integrated current J at a specified time t=T in the symmetric simple exclusion process (SSEP) on an infinite line. As shown recently [Phys. Rev. E 89, 010101(R) (2014)], the SSEP belongs to the elliptic universality class. Here, for very large currents, the diffusion terms of the MFT equations can be neglected compared with the terms coming from the shot noise. Using the hodograph transformation and an additional change of variables, we reduce the "inviscid" MFT equations to Laplace's equation in an extended space. This opens the way to an exact solution. Here we solve the extreme-current problem for a flat deterministic initial density profile with an arbitrary density 0<n<1. The solution yields the most probable density history of the system conditional on the extreme current, and leads to a super-Gaussian extreme-current statistics, - ln P = F(n) J^3/T, in agreement with Derrida and Gerschenfeld [J. Stat. Phys. 137, 978 (2009)]. We calculate the function F(n) analytically. It is symmetric with respect to the half-filling density n=1/2, diverges as n approached 0 or 1, and exhibits a singularity F(n) |n-1/2| at the half-filling density n=1/2.

preprint2014arXiv

Large fluctuations in diffusion-controlled absorption

Suppose that $N_0$ independently diffusing particles, each with diffusivity $D$, are initially released at $x=\ell>0$ on the semi-infinite interval $0\leq x<\infty$ with an absorber at $x=0$. We determine the probability ${\cal P}(N)$ that $N$ particles survive until time $t=T$. We also employ macroscopic fluctuation theory to find the most likely history of the system, conditional on there being exactly $N$ survivors at time $t=T$. Depending on the basic parameter $\ell/\sqrt{4DT}$, very different histories can contribute to the extreme cases of $N=N_0$ (all particles survive) and $N=0$ (no survivors). For large values of $\ell/\sqrt{4DT}$, the leading contribution to ${\cal P}(N=0)$ comes from an effective point-like quasiparticle that contains all the $N_0$ particles and moves ballistically toward the absorber until absorption occurs.

preprint2014arXiv

Survival of a static target in a gas of diffusing particles with exclusion

Let a lattice gas of constant density, described by the symmetric simple exclusion process, be brought in contact with a "target": a spherical absorber of radius $R$. Employing the macroscopic fluctuation theory (MFT), we evaluate the probability ${\mathcal P}(T)$ that no gas particle hits the target until a long but finite time $T$. We also find the most likely gas density history conditional on the non-hitting. The results depend on the dimension of space $d$ and on the rescaled parameter $\ell=R/\sqrt{D_0T}$, where $D_0$ is the gas diffusivity. For small $\ell$ and $d>2$, ${\mathcal P}(T)$ is determined by an exact stationary solution of the MFT equations that we find. For large $\ell$, and for any $\ell$ in one dimension, the relevant MFT solutions are non-stationary. In this case $\ln {\mathcal P}(T)$ scales differently with relevant parameters, and it also depends on whether the initial condition is random or deterministic. The latter effects also occur if the lattice gas is composed of non-interacting random walkers. Finally, we extend the formalism to a whole class of diffusive gases of interacting particles.

preprint2014arXiv

Two-strain competition in quasi-neutral stochastic disease dynamics

We develop a new perturbation method for studying quasi-neutral competition in a broad class of stochastic competition models, and apply it to the analysis of fixation of competing strains in two epidemic models. The first model is a two-strain generalization of the stochastic Susceptible-Infected-Susceptible (SIS) model. Here we extend previous results due to Parsons and Quince (2007), Parsons et al (2008) and Lin, Kim and Doering (2012). The second model, a two-strain generalization of the stochastic Susceptible-Infected-Recovered (SIR) model with population turnover, has not been studied previously. In each of the two models, when the basic reproduction numbers of the two strains are identical, a system with an infinite population size approaches a point on the deterministic coexistence line (CL): a straight line of fixed points in the phase space of sub-population sizes. Shot noise drives one of the strain populations to fixation, and the other to extinction, on a time scale proportional to the total population size. Our perturbation method explicitly tracks the dynamics of the probability distribution of the sub-populations in the vicinity of the CL. We argue that, whereas the slow strain has a competitive advantage for mathematically "typical" initial conditions, it is the fast strain that is more likely to win in the important situation when a few infectives of both strains are introduced into a susceptible population.

preprint2013arXiv

Extreme Current Fluctuations in a Nonstationary Stochastic Heat Flow

We employ the Hamiltonian formalism of macroscopic fluctuation theory to study large deviations of integrated current in the Kipnis-Marchioro-Presutti (KMP) model of stochastic hear flow when starting from a step-like initial condition. The KMP model belongs to the hyperbolic universality class where diffusion remains relevant no matter how large the fluctuating current is. The extreme current statistics for the KMP model turns out to be sub-Gaussian, as distinguished from the super-Gaussian statistics found for the Symmetric Simple Exclusion Process and other models of the elliptic class. The most probable time history of the system, which dominates the extreme current statistics of the KMP model, involves two large-amplitude solitary pulses: of the energy density field and of the conjugate "momentum" field. The coupled pulses propagate with a constant speed, but their amplitudes slowly grow with time, as the energy density pulse collects most of the available energy on its way.

preprint2013arXiv

Extreme Current Fluctuations in Lattice Gases: Beyond Nonequilibrium Steady States

We use the macroscopic fluctuation theory (MFT) to study large current fluctuations in non-stationary diffusive lattice gases. We identify two universality classes of these fluctuations which we call elliptic and hyperbolic. They emerge in the limit when the deterministic mass flux is small compared to the mass flux due to the shot noise. The two classes are determined by the sign of compressibility of \emph{effective fluid}, obtained by mapping the MFT into an inviscid hydrodynamics. An example of the elliptic class is the Symmetric Simple Exclusion Process where, for some initial conditions, we can solve the effective hydrodynamics exactly. This leads to a super-Gaussian extreme current statistics conjectured by Derrida and Gerschenfeld (2009) and yields the optimal path of the system. For models of the hyperbolic class the deterministic mass flux cannot be neglected, leading to a different extreme current statistics.

preprint2013arXiv

Immigration-extinction dynamics of stochastic populations

How high should be the rate of immigration into a stochastic population in order to significantly reduce the probability of observing the population extinct? Is there any relation between the population size distributions with and without immigration? Under what conditions can one justify the simple patch occupancy models which ignore the population distribution and its dynamics in a patch, and treat a patch simply as either occupied or empty? We address these questions by exactly solving a simple stochastic model obtained by adding a steady immigration to a variant of the Verhulst model: a prototypical model of an isolated stochastic population.

preprint2013arXiv

Impact of Colored Environmental Noise on the Extinction of a Long-Lived Stochastic Population: Role of the Allee Effect

We study the combined impact of a colored environmental noise and demographic noise on the extinction risk of a long-lived and well-mixed isolated stochastic population which exhibits the Allee effect. The environmental noise modulates the population birth and death rates. Assuming that the Allee effect is strong, and the environmental noise is positively correlated and Gaussian, we derive a Fokker-Planck equation for the joint probability distribution of the population sizes and environmental fluctuations. In WKB approximation this equation reduces to an effective two-dimensional Hamiltonian mechanics, where the most likely path to extinction and the most likely environmental fluctuation are encoded in an instanton-like trajectory in the phase space. The mean time to extinction (MTE) is related to the mechanical action along this trajectory. We obtain new analytic results for short-correlated, long-correlated and relatively weak environmental noise. The population-size dependence of the MTE changes from exponential for weak environmental noise to no dependence for strong noise, implying a greatly increased extinction risk. The theory is readily extendable to population switches between different metastable states, and to stochastic population explosion, due to a combined action of demographic and environmental noise.

preprint2012arXiv

Emergence of fluctuating traveling front solutions in macroscopic theory of noisy invasion fronts

The position of an invasion front, propagating into an unstable state, fluctuates because of the shot noise coming from the discreteness of reacting particles and stochastic character of the reactions and diffusion. A recent macroscopic theory [Meerson and Sasorov, Phys. Rev. E 84, 030101(R) (2011)] yields the probability of observing, during a long time, an unusually slow front. The theory is formulated as an effective classical Hamiltonian field theory which operates with the density field and the conjugate "momentum" field. Further, the theory assumes that the most probable density field history of an unusually slow front represents, up to small corrections, a traveling front solution of the Hamilton equations. Here we verify this assumption by solving the Hamilton equations numerically for models belonging to the directed percolation universality class.

preprint2012arXiv

Fluctuations of Current in Non-Stationary Diffusive Lattice Gases

We employ the macroscopic fluctuation theory to study fluctuations of integrated current in one-dimensional lattice gases with a step-like initial density profile. We analytically determine the variance of the current fluctuations for a class of diffusive processes with a density-independent diffusion coefficient, but otherwise arbitrary. Our calculations rely on a perturbation theory around the noiseless hydrodynamic solution. We consider both quenched and annealed types of averaging (the initial condition is allowed to fluctuate in the latter situation). The general results for the variance are specialized to a few interesting models including the symmetric exclusion process and the Kipnis-Marchioro-Presutti model. We also probe large deviations of the current for the symmetric exclusion process. This is done by numerically solving the governing equations of the macroscopic fluctuation theory using an efficient iteration algorithm.

preprint2012arXiv

Minimizing the population extinction risk by migration

Many populations in nature are fragmented: they consist of local populations occupying separate patches. A local population is prone to extinction due to the shot noise of birth and death processes. A migrating population from another patch can dramatically delay the extinction. What is the optimal migration rate that minimizes the extinction risk of the whole population? Here we answer this question for a connected network of model habitat patches with different carrying capacities.

preprint2012arXiv

Multiple extinction routes in stochastic population models

Isolated populations ultimately go extinct because of the intrinsic noise of elementary processes. In multi-population systems extinction of a population may occur via more than one route. We investigate this generic situation in a simple predator-prey (or infected-susceptible) model. The predator and prey populations may coexist for a long time but ultimately both go extinct. In the first extinction route the predators go extinct first, whereas the prey thrive for a long time and then also go extinct. In the second route the prey go extinct first causing a rapid extinction of the predators. Assuming large sub-population sizes in the coexistence state, we compare the probabilities of each of the two extinction routes and predict the most likely path of the sub-populations to extinction. We also suggest an effective three-state master equation for the probabilities to observe the coexistence state, the predator-free state and the empty state.

preprint2012arXiv

Survival of the Scarcer

We investigate extinction dynamics in the paradigmatic model of two competing species A and B that reproduce (A-->2A, B-->2B), self-regulate by annihilation (2A-->0, 2B-->0), and compete (A+B-->A, A+B-->B). For a finite system that is in the well-mixed limit, a quasi-stationary state arises which describes coexistence of the two species. Because of discrete noise, both species eventually become extinct in time that is exponentially long in the quasi-stationary population size. For a sizable range of asymmetries in the growth and competition rates, the paradoxical situation arises in which the numerically disadvantaged species according to the deterministic rate equations survives much longer.

preprint2012arXiv

Velocity fluctuations of noisy reaction fronts propagating into a metastable state: testing theory in stochastic simulations

The position of a reaction front, propagating into a metastable state, fluctuates because of the shot noise of reactions and diffusion. A recent theory [B. Meerson, P.V. Sasorov, and Y. Kaplan, Phys. Rev. E 84, 011147 (2011)] gave a closed analytic expression for the front diffusion coefficient in the weak noise limit. Here we test this theory in stochastic simulations involving reacting and diffusing particles on a one-dimensional lattice. We also investigate a small noise-induced systematic shift of the front velocity compared to the prediction from the spatially continuous deterministic reaction-diffusion equation.

preprint2012arXiv

Void formation in diffusive lattice gases

What is the probability that a macroscopic void will spontaneously arise, at a specified time T, in an initially homogeneous gas? We address this question for diffusive lattice gases, and also determine the most probable density history leading to the void formation. We employ the macroscopic fluctuation theory by Bertini et al and consider both annealed and quenched averaging procedures (the initial condition is allowed to fluctuate in the annealed setting). We show that in the annealed case the void formation probability is given by the equilibrium Boltzmann-Gibbs formula, so the probability is independent of T (and also of the void shape, as only the volume matters). In the quenched case, which is intrinsically non-equilibrium, we evaluate the void formation probability analytically for non-interacting random walkers and probe it numerically for the simple symmetric exclusion process. For voids that are small compared with the diffusion length, the equilibrium result for the void formation probability is recovered. We also re-derive our main results for non-interacting random walkers from an exact microscopic analysis.

preprint2011arXiv

Extinction rates of established spatial populations

This paper deals with extinction of an isolated population caused by intrinsic noise. We model the population dynamics in a "refuge" as a Markov process which involves births and deaths on discrete lattice sites and random migrations between neighboring sites. In extinction scenario I the zero population size is a repelling fixed point of the on-site deterministic dynamics. In extinction scenario II the zero population size is an attracting fixed point, corresponding to what is known in ecology as Allee effect. Assuming a large population size, we develop WKB (Wentzel-Kramers-Brillouin) approximation to the master equation. The resulting Hamilton's equations encode the most probable path of the population toward extinction and the mean time to extinction. In the fast-migration limit these equations coincide, up to a canonical transformation, with those obtained, in a different way, by Elgart and Kamenev (2004). We classify possible regimes of population extinction with and without an Allee effect and for different types of refuge and solve several examples analytically and numerically. For a very strong Allee effect the extinction problem can be mapped into the over-damped limit of theory of homogeneous nucleation due to Langer (1969). In this regime, and for very long systems, we predict an optimal refuge size that maximizes the mean time to extinction.

preprint2011arXiv

Negative velocity fluctuations of pulled reaction fronts

The position of a reaction front, propagating into an unstable state, fluctuates because of the shot noise. What is the probability that the fluctuating front moves considerably slower than its deterministic counterpart? Can the noise arrest the front motion for some time, or even make it move in the wrong direction? We present a WKB theory that assumes many particles in the front region and answers these questions for the microscopic model A->2A, 2A->A and random walk.

preprint2011arXiv

Switching between phenotypes and population extinction

Many types of bacteria can survive under stress by switching stochastically between two different phenotypes: the "normals" who multiply fast, but are vulnerable to stress, and the "persisters" who hardly multiply, but are resilient to stress. Previous theoretical studies of such bacterial populations have focused on the \emph{fitness}: the asymptotic rate of unbounded growth of the population. Yet for an isolated population of established (and not very large) size, a more relevant measure may be the population \emph{extinction risk} due to the interplay of adverse extrinsic variations and intrinsic noise of birth, death and switching processes. Applying a WKB approximation to the pertinent master equation of such a two-population system, we quantify the extinction risk, and find the most likely path to extinction under both favorable and adverse conditions. Analytical results are obtained both in the biologically relevant regime when the switching is rare compared with the birth and death processes, and in the opposite regime of frequent switching. We show that rare switches are most beneficial in reducing the extinction risk.

preprint2011arXiv

Velocity fluctuations of population fronts propagating into metastable states

The position of propagating population fronts fluctuates because of the discreteness of the individuals and stochastic character of processes of birth, death and migration. Here we consider a Markov model of a population front propagating into a metastable state, and focus on the weak noise limit. For typical, small fluctuations the front motion is diffusive, and we calculate the front diffusion coefficient. We also determine the probability distribution of rare, large fluctuations of the front position and, for a given average front velocity, find the most likely population density profile of the front. Implications of the theory for population extinction risk are briefly considered.

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.

preprint2010arXiv

Navier-Stokes hydrodynamics of thermal collapse in a freely cooling granular gas

We employ Navier-Stokes granular hydrodynamics to investigate the long-time behavior of clustering instability in a freely cooling dilute granular gas in two dimensions. We find that, in circular containers, the homogeneous cooling state (HCS) of the gas loses its stability via a sub-critical pitchfork bifurcation. There are no time-independent solutions for the gas density in the supercritical region, and we present analytical and numerical evidence that the gas develops thermal collapse unarrested by heat diffusion. To get more insight, we switch to a simpler geometry of a narrow-sector-shaped container. Here the HCS loses its stability via a transcritical bifurcation. For some initial conditions a time-independent inhomogeneous density profile sets in, qualitatively similar to that previously found in a narrow-channel geometry. For other initial conditions, however, the dilute gas develops thermal collapse unarrested by heat diffusion. We determine the dynamic scalings of the flow close to collapse analytically and verify them in hydrodynamic simulations. The results of this work imply that, in dimension higher than one, Navier-Stokes hydrodynamics of a dilute granular gas is prone to finite-time density blowups. This provides a natural explanation to the formation of densely packed clusters of particles in a variety of initially dilute granular flows.

preprint2010arXiv

Stochastic models of population extinction

Theoretical ecologists have long sought to understand how the persistence of populations depends on biotic and abiotic factors. Classical work showed that demographic stochasticity causes the mean time to extinction to increase exponentially with population size, whereas variation in environmental conditions can lead to a power law scaling. Recent work has focused especially on the influence of the autocorrelation structure ("color") of environmental noise. In theoretical physics, there is a burst of research activity in analyzing large fluctuations in stochastic population dynamics. This research provides powerful tools for determining extinction times and characterizing the pathway to extinction. It yields, therefore, sharp insights into extinction processes and has great potential for further applications in theoretical biology.

preprint2010arXiv

Time-resolved extinction rates of stochastic populations

Extinction of a long-lived isolated stochastic population can be described as an exponentially slow decay of quasi-stationary probability distribution of the population size. We address extinction of a population in a two-population system in the case when the population turnover -- renewal and removal -- is much slower than all other processes. In this case there is a time scale separation in the system which enables one to introduce a short-time quasi-stationary extinction rate W_1 and a long-time quasi-stationary extinction rate W_2, and develop a time-dependent theory of the transition between the two rates. It is shown that W_1 and W_2 coincide with the extinction rates when the population turnover is absent, and present but very slow, respectively. The exponentially large disparity between the two rates reflects fragility of the extinction rate in the population dynamics without turnover.

preprint2009arXiv

Emergence of stability in a stochastically driven pendulum: beyond the Kapitsa effect

We consider a prototypical nonlinear system which can be stabilized by multiplicative noise: an underdamped non-linear pendulum with a stochastically vibrating pivot. A numerical solution of the pertinent Fokker-Planck equation shows that the upper equilibrium point of the pendulum can become stable even when the noise is white, and the "Kapitsa pendulum" effect is not at work. The stabilization occurs in a strong-noise regime where WKB approximation does not hold.

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

Extinction of an infectious disease: a large fluctuation in a non-equilibrium system

We develop a theory of first passage processes in stochastic non-equilibrium systems of birth-death type using two closely related epidemiological models as examples. Our method employs the probability generating function technique in conjunction with the eikonal approximation. In this way the problem is reduced to finding the optimal path to extinction: a heteroclinic trajectory of an effective multi-dimensional classical Hamiltonian system. We compute this trajectory and mean extinction time of the disease numerically and uncover a non-monotone, spiral path to extinction of a disease. We also obtain analytical results close to a bifurcation point, where the problem is described by a Hamiltonian previously identified in one-species population models.

preprint2008arXiv

Noise driven unlimited population growth

Demographic noise causes unlimited population growth in a broad class of models which, without noise, would predict a stable finite population. We study this effect on the example of a stochastic birth-death model which includes immigration, binary reproduction and death. The unlimited population growth proceeds as an exponentially slow decay of a metastable probability distribution (MPD) of the population. We develop a systematic WKB theory, complemented by the van Kampen system size expansion, for the MPD and for the decay time. Important signatures of the MPD is a power-law tail (such that all the distribution moments, except the zeroth one, diverge) and the presence in the solution of two different WKB modes.

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 fluctuating environment

Environmental noise can cause an exponential reduction in the mean time to extinction (MTE) of an isolated population. We study this effect on an example of a stochastic birth-death process with rates modulated by a colored Gaussian noise. A path integral formulation yields a transparent way of evaluating the MTE and finding the optimal realization of the environmental noise that determines the most probable path to extinction. The population-size dependence of the MTE changes from exponential in the absence of the environmental noise to a power law for a short-correlated noise and to no dependence for long-correlated noise. We also establish the validity domains of the limits of white noise and adiabatic noise.

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.

preprint2008arXiv

The Knudsen temperature jump and the Navier-Stokes hydrodynamics of granular gases driven by thermal walls

Thermal wall is a convenient idealization of a rapidly vibrating plate used for vibrofluidization of granular materials. The objective of this work is to incorporate the Knudsen temperature jump at thermal wall in the Navier-Stokes hydrodynamic modeling of dilute granular gases of monodisperse particles that collide nearly elastically. The Knudsen temperature jump manifests itself as an additional term, proportional to the temperature gradient, in the boundary condition for the temperature. Up to a numerical pre-factor of order unity, this term is known from kinetic theory of elastic gases. We determine the previously unknown numerical pre-factor by measuring, in a series of molecular dynamics (MD) simulations, steady-state temperature profiles of a gas of elastically colliding hard disks, confined between two thermal walls kept at different temperatures, and comparing the results with the predictions of a hydrodynamic calculation employing the modified boundary condition. The modified boundary condition is then applied, without any adjustable parameters, to a hydrodynamic calculation of the temperature profile of a gas of inelastic hard disks driven by a thermal wall. We find the hydrodynamic prediction to be in very good agreement with MD simulations of the same system. The results of this work pave the way to a more accurate hydrodynamic modeling of driven granular gases.

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.

preprint2007arXiv

van der Waals-like phase separation instability of a driven granular gas in three dimensions

We show that the van der Waals-like phase separation instability of a driven granular gas at zero gravity, previously investigated in two-dimensional settings, persists in three dimensions. We consider a monodisperse granular gas driven by a thermal wall of a three-dimensional rectangular container at zero gravity. The basic steady state of this system, as described by granular hydrodynamic equations, involves a denser and colder layer of granulate located at the wall opposite to the driving wall. When the inelastic energy loss is sufficiently high, the driven granular gas exhibits, in some range of average densities, negative compressibility in the directions parallel to the driving wall. When the lateral dimensions of the container are sufficiently large, the negative compressibility causes spontaneous symmetry breaking of the basic steady state and a phase separation instability. Event-driven molecular dynamics simulations confirm and complement our theoretical predictions.

preprint2006arXiv

Phase separation of a driven granular gas in annular geometry

This work investigates phase separation of a monodisperse gas of inelastically colliding hard disks confined in a two-dimensional annulus, the inner circle of which represents a "thermal wall". When described by granular hydrodynamic equations, the basic steady state of this system is an azimuthally symmetric state of increased particle density at the exterior circle of the annulus. When the inelastic energy loss is sufficiently large, hydrodynamics predicts spontaneous symmetry breaking of the annular state, analogous to the van der Waals-like phase separation phenomenon previously found in a driven granular gas in rectangular geometry. At a fixed aspect ratio of the annulus, the phase separation involves a "spinodal interval" of particle area fractions, where the gas has negative compressibility in the azimuthal direction. The heat conduction in the azimuthal direction tends to suppress the instability, as corroborated by a marginal stability analysis of the basic steady state with respect to small perturbations. To test and complement our theoretical predictions we performed event-driven molecular dynamics (MD) simulations of this system. We clearly identify the transition to phase separated states in the MD simulations, despite large fluctuations present, by measuring the probability distribution of the amplitude of the fundamental Fourier mode of the azimuthal spectrum of the particle density. We find that the instability region, predicted from hydrodynamics, is always located within the phase separation region observed in the MD simulations. This implies the presence of a binodal (coexistence) region, where the annular state is metastable. The phase separation persists when the driving and elastic walls are interchanged, and also when the elastic wall is replaced by weakly inelastic one.

preprint2005arXiv

Towards a continuum theory of clustering in a freely cooling inelastic gas

We performed molecular dynamics simulations to investigate the clustering instability of a freely cooling dilute gas of inelastically colliding disks in a quasi-one-dimensional setting. We observe that, as the gas cools, the shear stress becomes negligibly small, and the gas flows by inertia only. Finite-time singularities, intrinsic in such a flow, are arrested only when close-packed clusters are formed. We observe that the late-time dynamics of this system are describable by the Burgers equation with vanishing viscosity, and predict the long-time coarsening behavior.

preprint2003arXiv

Development of a density inversion in driven granular gases

Granular materials fluidized by a rapidly vibrating bottom plate often develop a fascinating density inversion: a heavier layer of granulate supported by a lower-density region. We employ the Navier-Stokes granular hydrodynamics to follow a density inversion as it develops in time. Assuming a dilute low-Mach-number flow, we derive a reduced time-dependent model of the late stage of the dynamics. The model looks especially simple in the Lagrangian coordinates. The time-dependent solution describes the growth of a density peak at an intermediate height. A transient temperature minimum is predicted to develop in the region of the density peak. The temperature minimum disappears at later times, as the system approaches the steady state. At late times, the predictions of the low-Mach-number model are in good agreement with a numerical solution of the full hydrodynamic equations. At an early stage of the dynamics, pressure oscillations are predicted.

preprint2003arXiv

Giant fluctuations at a granular phase separation threshold

We investigate a phase separation instability that occurs in a system of nearly elastically colliding hard spheres driven by a thermal wall. If the aspect ratio of the confining box exceeds a threshold value, granular hydrostatics predict phase separation: the formation of a high-density region coexisting with a low-density region along the wall that is opposite to the thermal wall. Event-driven molecular dynamic simulations confirm this prediction. The theoretical bifurcation curve agrees with the simulations quantitatively well below and well above the threshold. However, in a wide region of aspect ratios around the threshold, the system is dominated by fluctuations, and the hydrostatic theory breaks down. Two possible scenarios of the origin of the giant fluctuations are discussed.

preprint1997arXiv

Universality in Dynamic Coarsening of a Fractal Cluster

Dynamics of coarsening of a statistically homogeneous fractal cluster, created by a morphological instability of diffusion-controlled growth, is investigated theoretically. An exact mathematical setting of the problem is presented that obeys a global conservation law. A statistical mean field theory is developed that accounts for shadowing during the growth instability and assumes that the total mass and fractal dimension of the cluster remain constant. The coarsening dynamics are shown to be self-similar, and the dynamic scaling exponents are calculated for any Euclidean dimension.