Source author record

Pavel V. Sasorov

Pavel V. Sasorov 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

16works
9topics
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

16 published item(s)

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.

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.

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.

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

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

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

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.

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

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.

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.