Source author record

Satya N. Majumdar

Satya N. Majumdar 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

140works
24topics
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

140 published item(s)

preprint2026arXiv

The distribution of the maximum of independent resetting Brownian motions

The probability distribution of the maximum $M_t$ of a single resetting Brownian motion (RBM) of duration $t$ and resetting rate $r$, properly centred and scaled, is known to converge to the standard Gumbel distribution of the classical extreme value theory. This Gumbel law describes the typical fluctuations of $M_t$ around its average $\sim \ln (r t)$ for large $t$ on a scale of $O(1)$. Here we compute the large-deviation tails of this distribution when $M_t = O(t)$ and show that the large-deviation function has a singularity where the second derivative is discontinuous, signalling a dynamical phase transition. Then we consider a collection of independent RBMs with initial (and resetting) positions uniformly distributed with a density $ρ$ over the negative half-line. We show that the fluctuations in the initial positions of the particles modify the distribution of $M_t$. The average over the initial conditions can be performed in two different ways, in analogy with disordered systems: (i) the annealed case where one averages over all possible initial conditions and (ii) the quenched case where one considers only the contributions coming from typical initial configurations. We show that in the annealed case, the limiting distribution of the maximum is characterized by a new scaling function, different from the Gumbel law but the large-deviation function remains the same as in the single particle case. In contrast, for the quenched case, the limiting (typical) distribution remains Gumbel but the large-deviation behaviors are new and nontrivial. Our analytical results, both for the typical as well as for the large-deviation regime of $M_t$, are verified numerically with extremely high precision, down to $10^{-250}$ for the probability density of $M_t$.

preprint2024arXiv

Exact extreme, order and sum statistics in a class of strongly correlated system

Even though strongly correlated systems are abundant, only a few exceptional cases admit analytical solutions. In this paper we present a large class of solvable systems with strong correlations.. We consider a set of $N$ independent and identically distributed (i.i.d) random variables $\{X_1,\, X_2,\ldots, X_N\}$ whose common distribution has a parameter $Y$ (or a set of parameters) which itself is random with its own distribution. For a fixed value of this parameter $Y$, the $X_i$ variables are independent and we call them conditionally independent and identically distributed (c.i.i.d). However, once integrated over the distribution of the parameter $Y$, the $X_i$ variables get strongly correlated, yet retaining a solvable structure for various observables, such as for the sum and the extremes of $X_i$'s. This provides a simple procedure to generate a class of solvable strongly correlated systems. We illustrate how this procedure works via three physical examples where $N$ particles on a line perform independent (i) Brownian motions, (ii) ballistic motions with random initial velocities, and (iii) Lévy flights, but they get strongly correlated via {\it simultaneous resetting} to the origin. Our results are verified in numerical simulations. This procedure can be used to generate an endless variety of solvable strongly correlated systems.

preprint2023arXiv

Metropolis Monte Carlo sampling: convergence, localization transition and optimality

Among random sampling methods, Markov Chain Monte Carlo algorithms are foremost. Using a combination of analytical and numerical approaches, we study their convergence properties towards the steady state, within a random walk Metropolis scheme. Analysing the relaxation properties of some model algorithms sufficiently simple to enable analytic progress, we show that the deviations from the target steady-state distribution can feature a localization transition as a function of the characteristic length of the attempted jumps defining the random walk. While the iteration of the Monte Carlo algorithm converges to equilibrium for all choices of jump parameters, the localization transition changes drastically the asymptotic shape of the difference between the probability distribution reached after a finite number of steps of the algorithm and the target equilibrium distribution. We argue that the relaxation before and after the localisation transition is respectively limited by diffusion and rejection rates.

preprint2022arXiv

An exactly solvable predator prey model with resetting

We study a simple model of a diffusing particle (the prey) that on encounter with one of a swarm of diffusing predators can either perish or be reset to its original position at the origin. We show that the survival probability of the prey up to time $t$ decays algebraically as $\sim t^{-θ(p, γ)}$ where the exponent $θ$ depends continuously on two parameters of the model, with $p$ denoting the probability that a prey survives upon encounter with a predator and $γ= D_A/(D_A+D_B)$ where $D_A$ and $D_B$ are the diffusion constants of the prey and the predator respectively. We also compute exactly the probability distribution $P(N|t_c)$ of the total number of encounters till the capture time $t_c$ and show that it exhibits an anomalous large deviation form $P(N|t_c)\sim t_c^{- Φ\left(\frac{N}{\ln t_c}=z\right)}$ for large $t_c$. The rate function $Φ(z)$ is computed explicitly. Numerical simulations are in excellent agreement with our analytical results.

preprint2022arXiv

Condensation transition in large deviations of self-similar Gaussian processes with stochastic resetting

We study the fluctuations of the area $A(t)= \int_0^t x(τ)\, dτ$ under a self-similar Gaussian process (SGP) $x(τ)$ with Hurst exponent $H>0$ (e.g., standard or fractional Brownian motion, or the random acceleration process) that stochastically resets to the origin at rate $r$. Typical fluctuations of $A(t)$ scale as $\sim \sqrt{t}$ for large $t$ and on this scale the distribution is Gaussian, as one would expect from the central limit theorem. Here our main focus is on atypically large fluctuations of $A(t)$. In the long-time limit $t\to\infty$, we find that the full distribution of the area takes the form $P_{r}\left(A|t\right)\sim\exp\left[-t^αΦ\left(A/t^β\right)\right]$ with anomalous exponents $α=1/(2H+2)$ and $β= (2H+3)/(4H+4)$ in the regime of moderately large fluctuations, and a different anomalous scaling form $P_{r}\left(A|t\right)\sim\exp\left[-tΨ\left(A/t^{\left(2H+3\right)/2}\right)\right]$ in the regime of very large fluctuations. The associated rate functions $Φ(y)$ and $Ψ(w)$ depend on $H$ and are found exactly. Remarkably, $Φ(y)$ has a singularity that we interpret as a first-order dynamical condensation transition, while $Ψ(w)$ exhibits a second-order dynamical phase transition above which the number of resetting events ceases to be extensive. The parabolic behavior of $Φ(y)$ around the origin $y=0$ correctly describes the typical, Gaussian fluctuations of $A(t)$. Despite these anomalous scalings, we find that all of the cumulants of the distribution $P_{r}\left(A|t\right)$ grow linearly in time, $\langle A^n\rangle_c\approx c_n \, t$, in the long-time limit. For the case of reset Brownian motion (corresponding to $H=1/2$), we develop a recursive scheme to calculate the coefficients $c_n$ exactly and use it to calculate the first 6 nonvanishing cumulants.

preprint2022arXiv

Counting statistics for non-interacting fermions in a rotating trap

We study the ground state of $N \gg 1$ noninteracting fermions in a two-dimensional harmonic trap rotating at angular frequency $Ω>0$. The support of the density of the Fermi gas is a disk of radius $R_e$. We calculate the variance of the number of fermions ${\cal N}_R$ inside a disk of radius $R$ centered at the origin for $R$ in the bulk of the Fermi gas. We find rich and interesting behaviours in two different scaling regimes: (i) $Ω/ ω<1 $ and (ii) $1 - Ω/ ω= O(1/N)$, where $ω$ is the angular frequency of the oscillator. In the first regime (i) we find that ${\rm Var}\,{\cal N}_{R}\simeq\left(A\log N+B\right)\sqrt{N}$ and we calculate $A$ and $B$ as functions of $R/R_e$, $Ω$ and $ω$. We also predict the higher cumulants of ${\cal N}_{R}$ and the bipartite entanglement entropy of the disk with the rest of the system. In the second regime (ii), the mean fermion density exhibits a staircase form, with discrete plateaus corresponding to filling $k$ successive Landau levels, as found in previous studies. Here, we show that ${\rm Var}\,{\cal N}_{R}$ is a discontinuous piecewise linear function of $\sim (R/R_e) \sqrt{N}$ within each plateau, with coefficients that we calculate exactly, and with steps whose precise shape we obtain for any $k$. We argue that a similar piecewise linear behavior extends to all the cumulants of ${\cal N}_{R}$ and to the entanglement entropy. We show that these results match smoothly at large $k$ with the above results for $Ω/ω=O(1)$. These findings are nicely confirmed by numerical simulations. Finally, we uncover a universal behavior of ${\rm Var}\,{\cal N}_{R}$ near the fermionic edge. We extend our results to a three-dimensional geometry, where an additional confining potential is applied in the $z$ direction.

preprint2022arXiv

Freezing transitions of Brownian particles in confining potentials

We study the mean first passage time (MFPT) to an absorbing target of a one-dimensional Brownian particle subject to an external potential $v(x)$ in a finite domain. We focus on the cases in which the external potential is confining, of the form $v(x)=k|x-x_0|^n/n$, and where the particle's initial position coincides with $x_0$. We first consider a particle between an absorbing target at $x=0$ and a reflective wall at $x=c$. At fixed $x_0$, we show that when the target distance $c$ exceeds a critical value, there exists a nonzero optimal stiffness $k_{\rm opt}$ that minimizes the MFPT to the target. However, when $c$ lies below the critical value, the optimal stiffness $k_{\rm opt}$ vanishes. Hence, for any value of $n$, the optimal potential stiffness undergoes a continuous "freezing" transition as the domain size is varied. On the other hand, when the reflective wall is replaced by a second absorbing target, the freezing transition in $k_{\rm opt}$ becomes discontinuous. The phase diagram in the $(x_0,n)$-plane then exhibits three dynamical phases and metastability, with a "triple" point at $(x_0/c\simeq 0.17185$, $n\simeq 0.39539)$. For harmonic or higher order potentials $(n\ge 2)$, the MFPT always increases with $k$ at small $k$, for any $x_0$ or domain size. These results are contrasted with problems of diffusion under optimal resetting in bounded domains.

preprint2022arXiv

Gap probability and full counting statistics in the one dimensional one-component plasma

We consider the $1d$ one-component plasma (OCP) in thermal equilibrium, consisting of $N$ equally charged particles on a line, with pairwise Coulomb repulsion and confined by an external harmonic potential. We study two observables: (i) the distribution of the gap between two consecutive particles in the bulk and (ii) the distribution of the number of particles $N_I$ in a fixed interval $I=[-L,+L]$ inside the bulk, the so-called full-counting-statistics (FCS). For both observables, we compute, for large $N$, the distribution of the typical as well as atypical large fluctuations. We show that the distribution of the typical fluctuations of the gap are described by the scaling form ${\cal P}_{\rm gap, bulk}(g,N) \sim N H_α(g\,N)$, where $α$ is the interaction coupling and the scaling function $H_α(z)$ is computed explicitly. It has a faster than Gaussian tail for large $z$: $H_α(z) \sim e^{-z^3/(96 α)}$ as $z \to \infty$. Similarly, for the FCS, we show that the distribution of the typical fluctuations of $N_I$ is described by the scaling form ${\cal P}_{\rm FCS}(N_I,N) \sim 2α\, U_α[2 α(N_I - \bar{N}_I)]$, where $\bar{N}_I = L\,N/(2 α)$ is the average value of $N_I$ and the scaling function $U_α(z)$ is obtained explicitly. For both observables, we show that the probability of large fluctuations are described by large deviations forms with respective rate functions that we compute explicitly. Our numerical Monte-Carlo simulations are in good agreement with our analytical predictions.

preprint2022arXiv

Integer partitions and exclusion statistics: Limit shapes and the largest part of Young diagrams

We compute the limit shapes of the Young diagrams of the minimal difference $p$ partitions and provide a simple physical interpretation for the limit shapes. We also calculate the asymptotic distribution of the largest part of the Young diagram and show that the scaled distribution has a Gumbel form for all $p$. This Gumbel statistics for the largest part remains unchanged even for general partitions of the form $E=\sum_i n_i i^{1/ν}$ with $ν>0$ where $n_i$ is the number of times the part $i$ appears.

preprint2022arXiv

Number of distinct sites visited by a resetting random walker

We investigate the number $V_p(n)$ of distinct sites visited by an $n$-step resetting random walker on a $d$-dimensional hypercubic lattice with resetting probability $p$. In the case $p=0$, we recover the well-known result that the average number of distinct sites grows for large $n$ as $\langle V_0(n)\rangle\sim n^{d/2}$ for $d<2$ and as $\langle V_0(n)\rangle\sim n$ for $d>2$. For $p>0$, we show that $\langle V_p(n)\rangle$ grows extremely slowly as $\sim \left[\log(n)\right]^d$. We observe that the recurrence-transience transition at $d=2$ for standard random walks (without resetting) disappears in the presence of resetting. In the limit $p\to 0$, we compute the exact crossover scaling function between the two regimes. In the one-dimensional case, we derive analytically the full distribution of $V_p(n)$ in the limit of large $n$. Moreover, for a one-dimensional random walker, we introduce a new observable, which we call imbalance, that measures how much the visited region is symmetric around the starting position. We analytically compute the full distribution of the imbalance both for $p=0$ and for $p>0$. Our theoretical results are verified by extensive numerical simulations.

preprint2022arXiv

Optimal Resetting Brownian Bridges

We introduce a resetting Brownian bridge as a simple model to study search processes where the total search time $t_f$ is finite and the searcher returns to its starting point at $t_f$. This is simply a Brownian motion with a Poissonian resetting rate $r$ to the origin which is constrained to start and end at the origin at time $t_f$. We first provide a rejection-free algorithm to generate such resetting bridges in all dimensions by deriving an effective Langevin equation with an explicit space-time dependent drift $\tilde μ({\bf x},t)$ and resetting rate $\tilde r({\bf x}, t)$. We also study the efficiency of the search process in one-dimension by computing exactly various observables such as the mean-square displacement, the hitting probability of a fixed target and the expected maximum. Surprisingly, we find that there exists an optimal resetting rate $r^*$ that maximizes the search efficiency, even in the presence of a bridge constraint. We show however that the physical mechanism responsible for this optimal resetting rate for bridges is entirely different from resetting Brownian motions without the bridge constraint.

preprint2022arXiv

Reducing mean first passage times with intermittent confining potentials: a realization of resetting processes

During a random search, resetting the searcher's position from time to time to the starting point often reduces the mean completion time of the process. Although many different resetting models have been studied over the past ten years, only a few can be physically implemented. Here we study theoretically a protocol that can be realised experimentally and which exhibits unusual optimization properties. A Brownian particle is subject to an arbitrary confining potential $v(x)$ which is switched on and off intermittently at fixed rates. Motion is constrained between an absorbing wall located at the origin and a reflective wall. When the walls are sufficiently far apart, the interplay between free diffusion during the "off" phases and attraction toward the potential minimum during the "on" phases gives rise to rich behaviours, not observed in ideal resetting models. For potentials of the form $v(x)=k|x-x_0|^n/n$, with $n>0$, the switch-on and switch-off rates that minimise the mean first passage time (MFPT) to the origin undergo a continuous phase transition as the potential stiffness $k$ is varied. When $k$ is above a critical value $k_c$, potential intermittency enhances target encounter: the minimal MFPT is lower than the Kramer's time and is attained for a non-vanishing pair of switching rates. We focus on the harmonic case $n=2$, extending previous results for the piecewise linear potential ($n=1$) in unbounded domains. We also study the non-equilibrium stationary states emerging in this process.

preprint2022arXiv

Run-and-tumble particle in one-dimensional confining potential: Steady state, relaxation and first passage properties

We study the dynamics of a one-dimensional run and tumble particle subjected to confining potentials of the type $V(x) = α\, |x|^p$, with $p>0$. The noise that drives the particle dynamics is telegraphic and alternates between $\pm 1$ values. We show that the stationary probability density $P(x)$ has a rich behavior in the $(p, α)$-plane. For $p>1$, the distribution has a finite support in $[x_-,x_+]$ and there is a critical line $α_c(p)$ that separates an active-like phase for $α> α_c(p)$ where $P(x)$ diverges at $x_\pm$, from a passive-like phase for $α< α_c(p)$ where $P(x)$ vanishes at $x_\pm$. For $p<1$, the stationary density $P(x)$ collapses to a delta function at the origin, $P(x) = δ(x)$. In the marginal case $p=1$, we show that, for $α< α_c$, the stationary density $P(x)$ is a symmetric exponential, while for $α> α_c$, it again is a delta function $P(x) = δ(x)$. For the special cases $p=2$ and $p=1$, we obtain exactly the full time-dependent distribution $P(x,t)$, that allows us to study how the system relaxes to its stationary state. In addition, in these two cases, we also study analytically the full distribution of the first-passage time to the origin. Numerical simulations are in complete agreement with our analytical predictions.

preprint2021arXiv

Edge fluctuations and third-order phase transition in harmonically confined long-range systems

We study the distribution of the position of the rightmost particle $x_{\max}$ in a $N$-particle Riesz gas in one dimension confined in a harmonic trap. The particles interact via long-range repulsive potential, of the form $r^{-k}$ with $-2<k<\infty$ where $r$ is the inter-particle distance. In equilibrium at temperature $O(1)$, the gas settles on a finite length scale $L_N$ that depends on $N$ and $k$. We numerically observe that the typical fluctuation of $y_{\max} = x_{\max}/L_N$ around its mean is of $O(N^{-η_k})$. Over this length scale, the distribution of the typical fluctuations has a $N$ independent scaling form. We show that the exponent $η_k$ obtained from the Hessian theory predicts the scale of typical fluctuations remarkably well. The distribution of atypical fluctuations to the left and right of the mean $\langle y_{\max} \rangle$ are governed by the left and right large deviation functions, respectively. We compute these large deviation functions explicitly $\forall k>-2$. We also find that these large deviation functions describe a pulled to pushed type phase transition as observed in Dyson's log-gas ($k\to 0$) and $1d$ one component plasma ($k=-1$). Remarkably, we find that the phase transition remains $3^{\rm rd}$ order for the entire regime. Our results demonstrate the striking universality of the $3^{\rm rd}$ order transition even in models that fall outside the paradigm of Coulomb systems and the random matrix theory. We numerically verify our analytical expressions of the large deviation functions via Monte Carlo simulation using an importance sampling algorithm.

preprint2021arXiv

Generating stochastic trajectories with global dynamical constraints

We propose a method to exactly generate Brownian paths $x_c(t)$ that are constrained to return to the origin at some future time $t_f$, with a given fixed area $A_f = \int_0^{t_f}dt\, x_c(t)$ under their trajectory. We derive an exact effective Langevin equation with an effective force that accounts for the constraint. In addition, we develop the corresponding approach for discrete-time random walks, with arbitrary jump distributions including Lévy flights, for which we obtain an effective jump distribution that encodes the constraint. Finally, we generalise our method to other types of dynamical constraints such as a fixed occupation time on the positive axis $T_f=\int_0^{t_f}dt\, Θ\left[x_c(t)\right]$ or a fixed generalised quadratic area $\mathcal{A}_f=\int_0^{t_f}dt \,x_c^2(t)$.

preprint2021arXiv

Localization transition in the Discrete Non-Linear Schrödinger Equation: ensembles inequivalence and negative temperatures

We present a detailed account of a first-order localization transition in the Discrete Nonlinear Schrödinger Equation, where the localized phase is associated to the high energy region in parameter space. We show that, due to ensemble inequivalence, this phase is thermodynamically stable only in the microcanonical ensemble. In particular, we obtain an explicit expression of the microcanonical entropy close to the transition line, located at infinite temperature. This task is accomplished making use of large-deviation techniques, that allow us to compute, in the limit of large system size, also the subleading corrections to the microcanonical entropy. These subleading terms are crucial ingredients to account for the first-order mechanism of the transition, to compute its order parameter and to predict the existence of negative temperatures in the localized phase. All of these features can be viewed as signatures of a thermodynamic phase where the translational symmetry is broken spontaneously due to a condensation mechanism yielding energy fluctuations far away from equipartition: actually they prefer to participate in the formation of nonlinear localized excitations (breathers), typically containing a macroscopic fraction of the total energy.

preprint2021arXiv

Mean area of the convex hull of a run and tumble particle in two dimensions

We investigate the statistics of the convex hull for a single run-and-tumble particle in two dimensions. Run-and-tumble particle, also known as persistent random walker, has gained significant interest in the recent years due to its biological application in modelling the motion of bacteria. We consider two different statistical ensembles depending on whether (i) the total number of tumbles $n$ or (ii) the total observation time $t$ is kept fixed. Benchmarking the results on perimeter, we study the statistical properties of the area of the convex hull for RTP. Exploiting the connections to extreme value statistics, we obtain exact analytical expressions for the mean area for both ensembles. For fixed-$t$ ensemble, we show that the mean possesses a scaling form in $γt$ (with $γ$ being the tumbling rate) and the corresponding scaling function is exactly computed. Interestingly, we find that it exhibits crossover from $\sim t^3$ scaling at small times $\left( t \ll γ^{-1} \right)$ to $\sim t$ scaling at large times $\left( t \gg γ^{-1} \right)$. On the other hand, for fixed-$n$ ensemble, the mean expectedly grows linearly with $n$ for $n \gg 1$. All our analytical findings are supported with numerical simulations.

preprint2021arXiv

Optimal mean first-passage time of a Brownian searcher with resetting in one and two dimensions: Experiments, theory and numerical tests

We study experimentally, numerically and theoretically the optimal mean time needed by a Brownian particle, freely diffusing either in one or two dimensions, to reach, within a tolerance radius $R_{\text tol}$, a target at a distance $L$ from an initial position in the presence of resetting. The reset position is Gaussian distributed with width $σ$. We derived and tested two resetting protocols, one with a periodic and one with random (Poissonian) resetting times. We computed and measured the full first-passage probability distribution that displays spectacular spikes immediately after each resetting time for close targets. We study the optimal mean first-passage time as a function of the resetting period/rate for different target distances (values of the ratios $b=L/σ$) and target size ($a=R_\text{tol}/L$). We find an interesting phase transition at a critical value of $b$, both in one and two dimensions. The details of the calculations as well as experimental setup and limitations are discussed.

preprint2021arXiv

Record statistics for random walks and Lévy flights with resetting

We compute exactly the mean number of records $\langle R_N \rangle$ for a time-series of size $N$ whose entries represent the positions of a discrete time random walker on the line. At each time step, the walker jumps by a length $η$ drawn independently from a symmetric and continuous distribution $f(η)$ with probability $1-r$ (with $0\leq r < 1$) and with the complementary probability $r$ it resets to its starting point $x=0$. This is an exactly solvable example of a weakly correlated time-series that interpolates between a strongly correlated random walk series (for $r=0$) and an uncorrelated time-series (for $(1-r) \ll 1$). Remarkably, we found that for every fixed $r \in [0,1[$ and any $N$, the mean number of records $\langle R_N \rangle$ is completely universal, i.e., independent of the jump distribution $f(η)$. In particular, for large $N$, we show that $\langle R_N \rangle$ grows very slowly with increasing $N$ as $\langle R_N \rangle \approx (1/\sqrt{r})\, \ln N$ for $0<r <1$. We also computed the exact universal crossover scaling functions for $\langle R_N \rangle$ in the two limits $r \to 0$ and $r \to 1$. Our analytical predictions are in excellent agreement with numerical simulations.

preprint2021arXiv

Stability of large complex systems with heterogeneous relaxation dynamics

We study the probability of stability of a large complex system of size $N$ within the framework of a generalized May model, which assumes a linear dynamics of each population size $n_i$ (with respect to its equilibrium value): $ \frac{\mathrm{d}\, n_i}{\mathrm{d}t} = - a_i n_i - \sqrt{T} \sum_{j} J_{ij} n_j $. The $a_i>0$'s are the intrinsic decay rates, $J_{ij}$ is a real symmetric $(N\times N)$ Gaussian random matrix and $\sqrt{T}$ measures the strength of pairwise interaction between different species. Unlike in May's original homogeneous model, each species has now an intrinsic damping $a_i$ that may differ from one another. As the interaction strength $T$ increases, the system undergoes a phase transition from a stable phase to an unstable phase at a critical value $T=T_c$. We reinterpret the probability of stability in terms of the hitting time of the level $b=0$ of an associated Dyson Brownian Motion (DBM), starting at the initial position $a_i$ and evolving in `time' $T$. In the large $N \to \infty$ limit, using this DBM picture, we are able to completely characterize $T_c$ for arbitrary density $μ(a)$ of the $a_i$'s. For a specific flat configuration $a_i = 1 + σ\frac{i-1}{N}$, we obtain an explicit parametric solution for the limiting (as $N\to \infty$) spectral density for arbitrary $T$ and $σ$. For finite but large $N$, we also compute the large deviation properties of the probability of stability on the stable side $T < T_c$ using a Coulomb gas representation.

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 $|μ|$.

preprint2021arXiv

Statistics of the maximum and the convex hull of a Brownian motion in confined geometries

We consider a Brownian particle with diffusion coefficient $D$ in a $d$-dimensional ball of radius $R$ with reflecting boundaries. We study the maximum $M_x(t)$ of the trajectory of the particle along the $x$-direction at time $t$. In the long time limit, the maximum converges to the radius of the ball $M_x(t) \to R$ for $t\to \infty$. We investigate how this limit is approached and obtain an exact analytical expression for the distribution of the fluctuations $Δ(t) = [R-M_x(t)]/R$ in the limit of large $t$ in all dimensions. We find that the distribution of $Δ(t)$ exhibits a rich variety of behaviors depending on the dimension $d$. These results are obtained by establishing a connection between this problem and the narrow escape time problem. We apply our results in $d=2$ to study the convex hull of the trajectory of the particle in a disk of radius $R$ with reflecting boundaries. We find that the mean perimeter $\langle L(t)\rangle$ of the convex hull exhibits a slow convergence towards the perimeter of the circle $2πR$ with a stretched exponential decay $2πR-\langle L(t)\rangle \propto \sqrt{R}(Dt)^{1/4} \,e^{-2\sqrt{2Dt}/R}$. Finally, we generalise our results to other confining geometries, such as the ellipse with reflecting boundaries. Our results are corroborated by thorough numerical simulations.

preprint2020arXiv

Active Brownian Motion in Two Dimensions

We study the dynamics of a single active Brownian particle (ABP) in two spatial dimensions. The ABP has an intrinsic time scale $D_R^{-1}$ set by the rotational diffusion constant $D_R$. We show that, at short-times $t \ll D_R^{-1}$, the presence of `activness' results in a strongly anisotropic and non-diffusive dynamics in the $(xy)$ plane. We compute exactly the marginal distributions of the $x$ and $y$ position coordinates along with the radial distribution, which are all shown to be non-Brownian. In addition, we show that, at early times, the ABP has anomalous first-passage properties, characterized by non-Brownian exponents.

preprint2020arXiv

Anderson-like localization transition of random walks with resetting

We study several lattice random walk models with stochastic resetting to previously visited sites which exhibit a phase transition between an anomalous diffusive regime and a localization regime where diffusion is suppressed. The localized phase settles above a critical resetting rate, or rate of memory use, and the probability density asymptotically adopts in this regime a non-equilibrium steady state similar to that of the well known problem of diffusion with resetting to the origin. The transition occurs because of the presence of a single impurity site where the resetting rate is lower than on other sites, and around which the walker spontaneously localizes. Near criticality, the localization length diverges with a critical exponent that falls in the same class as the self-consistent theory of Anderson localization of waves in random media. The critical dimensions are also the same in both problems. Our study provides analytically tractable examples of localization transitions in path-dependent, reinforced stochastic processes, which can be also useful for understanding spatial learning by living organisms.

preprint2020arXiv

Current fluctuations in non-interacting run-and-tumble particles in one-dimension

We present a general framework to study the distribution of the flux through the origin up to time $t$, in a non-interacting one-dimensional system of particles with a step initial condition with a fixed density $ρ$ of particles to the left of the origin. We focus principally on two cases: (i) when the particles undergo diffusive dynamics (passive case) and (ii) run-and-tumble dynamics for each particle (active case). In analogy with disordered systems, we consider the flux distribution both for the annealed and the quenched initial conditions, for the passive and active particles. In the annealed case, we show that, for arbitrary particle dynamics, the flux distribution is a Poissonian with a mean $μ(t)$ that we compute exactly in terms of the Green's function of the single particle dynamics. For the quenched case, we show that, for the run-and-tumble dynamics, the quenched flux distribution takes an anomalous large deviation form at large times $P_{\rm qu}(Q,t) \sim \exp\left[-ρ\, v_0\, γ\, t^2 ψ_{\rm RTP}\left(\frac{Q}{ρv_0\,t} \right) \right]$, where $γ$ is the rate of tumbling and $v_0$ is the ballistic speed between two successive tumblings. In this paper, we compute the rate function $ψ_{\rm RTP}(q)$ and show that it is nontrivial. Our method also gives access to the probability of the rare event that, at time $t$, there is no particle to the right of the origin. For diffusive and run-and-tumble dynamics, we find that this probability decays with time as a stretched exponential, $\sim \exp(-c\, \sqrt{t})$ where the constant $c$ can be computed exactly. We verify our results for these large deviations by using an importance sampling Monte-Carlo method.

preprint2020arXiv

Distribution of the Time Between Maximum and Minimum of Random Walks

We consider a one-dimensional Brownian motion of fixed duration $T$. Using a path-integral technique, we compute exactly the probability distribution of the difference $τ=t_{\min}-t_{\max}$ between the time $t_{\min}$ of the global minimum and the time $t_{\max}$ of the global maximum. We extend this result to a Brownian bridge, i.e. a periodic Brownian motion of period $T$. In both cases, we compute analytically the first few moments of $τ$, as well as the covariance of $t_{\max}$ and $t_{\min}$, showing that these times are anti-correlated. We demonstrate that the distribution of $τ$ for Brownian motion is valid for discrete-time random walks with $n$ steps and with a finite jump variance, in the limit $n\to \infty$. In the case of Lévy flights, which have a divergent jump variance, we numerically verify that the distribution of $τ$ differs from the Brownian case. For random walks with continuous and symmetric jumps we numerically verify that the probability of the event "$τ= n$" is exactly $1/(2n)$ for any finite $n$, independently of the jump distribution. Our results can be also applied to describe the distance between the maximal and minimal height of $(1+1)$-dimensional stationary-state Kardar-Parisi-Zhang interfaces growing over a substrate of finite size $L$. Our findings are confirmed by numerical simulations. Some of these results have been announced in a recent Letter [Phys. Rev. Lett. 123, 200201 (2019)].

preprint2020arXiv

Ising model with stochastic resetting

We study the stationary properties of the Ising model that, while evolving towards its equilibrium state at temperature $T$ according to the Glauber dynamics, is stochastically reset to its fixed initial configuration with magnetisation $m_0$ at a constant rate $r$. Resetting breaks detailed balance and drives the system to a non-equilibrium stationary state where the magnetisation acquires a nontrivial distribution, leading to a rich phase diagram in the $(T,r)$ plane. We establish these results exactly in one-dimension and present scaling arguments supported by numerical simulations in two-dimensions. We show that resetting gives rise to a novel "pseudo-ferro" phase in the $(T,r)$ plane for $r > r^*(T)$ and $T>T_c$ where $r^*(T)$ is a crossover line separating the pseudo-ferro phase from a paramagnetic phase. This pseudo-ferro phase is characterised by a non-zero typical magnetisation and a vanishing gap near $m=0$ of the magnetisation distribution.

preprint2020arXiv

Kernels for noninteracting fermions via a Green's function approach with applications to step potentials

The quantum correlations of $N$ noninteracting spinless fermions in their ground state can be expressed in terms of a two-point function called the kernel. Here we develop a general and compact method for computing the kernel in a general trapping potential in terms of the Green's function for the corresponding single particle Schrödinger equation. For smooth potentials the method allows a simple alternative derivation of the local density approximation for the density and of the sine kernel in the bulk part of the trap in the large $N$ limit. It also recovers the density and the kernel of the so-called {\em Airy gas} at the edge. This method allows to analyse the quantum correlations in the ground state when the potential has a singular part with a fast variation in space. For the square step barrier of height $V_0$, we derive explicit expressions for the density and for the kernel. For large Fermi energy $μ>V_0$ it describes the interpolation between two regions of different densities in a Fermi gas, each described by a different sine kernel. Of particular interest is the {\em critical point} of the square well potential when $μ=V_0$. In this critical case, while there is a macroscopic number of fermions in the lower part of the step potential, there is only a finite $O(1)$ number of fermions on the shoulder, and moreover this number is independent of $μ$. In particular, the density exhibits an algebraic decay $\sim 1/x^2$, where $x$ is the distance from the jump. Furthermore, we show that the critical behaviour around $μ= V_0$ exhibits universality with respect with the shape of the barrier. This is established (i) by an exact solution for a smooth barrier (the Woods-Saxon potential) and (ii) by establishing a general relation between the large distance behavior of the kernel and the scattering amplitudes of the single-particle wave-function.

preprint2020arXiv

Mean perimeter and area of the convex hull of a planar Brownian motion in the presence of resetting

We compute exactly the mean perimeter and the mean area of the convex hull of a $2$-d Brownian motion of duration $t$ and diffusion constant $D$, in the presence of resetting to the origin at a constant rate $r$. We show that for any $t$, the mean perimeter is given by $\langle L(t)\rangle= 2 π\sqrt{\frac{D}{r}}\, f_1(rt)$ and the mean area is given by $\langle A(t) \rangle= 2π\frac{D}{r}\, f_2(rt)$ where the scaling functions $f_1(z)$ and $f_2(z)$ are computed explicitly. For large $t\gg 1/r$, the mean perimeter grows extremely slowly as $\langle L(t)\rangle \propto \ln (rt)$ with time. Likewise, the mean area also grows slowly as $\langle A(t)\rangle \propto \ln^2(rt)$ for $t\gg 1/r$. Our exact results indicate that the convex hull, in the presence of resetting, approaches a circular shape at late times. Numerical simulations are in perfect agreement with our analytical predictions.

preprint2020arXiv

Noninteracting trapped Fermions in double-well potentials: inverted parabola kernel

We study a system of $N$ noninteracting spinless fermions in a confining, double-well potential in one dimension. When the Fermi energy is close to the value of the potential at its local maximum we show that physical properties, such as the average density and the fermion position correlation functions, display a universal behavior that depends only on the local properties of the potential near its maximum. This behavior describes the merging of two Fermi gases, which are disjoint at sufficiently low Fermi energies. We describe this behavior in terms of a new correlation kernel that we compute analytically and we call it the inverted parabola kernel". As an application, we calculate the mean and variance of the number of particles in an interval of size $2L$ centered around the position of the local maximum, for sufficiently small $L$. Finally, we discuss the possibility of observing our results in experiments, as well as the extensions to nonzero temperature and to higher space dimensions.

preprint2020arXiv

Optimal mean first-passage time for a Brownian searcher subjected to resetting: experimental and theoretical results

We study experimentally and theoretically the optimal mean time needed by a free diffusing Brownian particle to reach a target at a distance L from an initial position in the presence of resetting. Both the initial position and the resetting position are Gaussian distributed with width $σ$. We derived and tested two resetting protocols, one with a periodic and one with random (Poissonian) resetting times. We computed and measured the full first-passage probability distribution that displays spectacular spikes immediately after each resetting time for close targets. We study the optimal mean first-passage time as a function of the resetting period/rate for different values of the ratio b = L/$σ$ and find an interesting phase transtion at a critical value b = bc. For bc < b < $\infty$, there is a metastable optimum time which disappears for b < bc. The intrinsic diffculties in implementing these protocols in experiments are also discussed.

preprint2020arXiv

Position distribution in a generalised run and tumble process

We study a class of stochastic processes of the type $\frac{d^n x}{dt^n}= v_0\, σ(t)$ where $n>0$ is a positive integer and $σ(t)=\pm 1$ represents an `active' telegraphic noise that flips from one state to the other with a constant rate $γ$. For $n=1$, it reduces to the standard run and tumble process for active particles in one dimension. This process can be analytically continued to any $n>0$ including non-integer values. We compute exactly the mean squared displacement at time $t$ for all $n>0$ and show that at late times while it grows as $\sim t^{2n-1}$ for $n>1/2$, it approaches a constant for $n<1/2$. In the marginal case $n=1/2$, it grows very slowly with time as $\sim \ln t$. Thus the process undergoes a {\em localisation} transition at $n=1/2$. We also show that the position distribution $p_n(x,t)$ remains time-dependent even at late times for $n\ge 1/2$, but approaches a stationary time-independent form for $n<1/2$. The tails of the position distribution at late times exhibit a large deviation form, $p_n(x,t)\sim \exp\left[-γ\, t\, Φ_n\left(\frac{x}{x^*(t)}\right)\right]$, where $x^*(t)= v_0\, t^n/Γ(n+1)$. We compute the rate function $Φ_n(z)$ analytically for all $n>0$ and also numerically using importance sampling methods, finding excellent agreement between them. For three special values $n=1$, $n=2$ and $n=1/2$ we compute the exact cumulant generating function of the position distribution at all times $t$.

preprint2020arXiv

Quantum dynamics with stochastic reset

We study non-equilibrium dynamics of integrable and non-integrable closed quantum systems whose unitary evolution is interrupted with stochastic resets, characterized by a reset rate $r$, that project the system to its initial state. We show that the steady state density matrix of a non-integrable system, averaged over the reset distribution, retains its off-diagonal elements for any finite $r$. Consequently a generic observable $\hat O$, whose expectation value receives contribution from these off-diagonal elements, never thermalizes under such dynamics for any finite $r$. We demonstrate this phenomenon by exact numerical studies of experimentally realizable models of ultracold bosonic atoms in a tilted optical lattice. For integrable Dirac-like fermionic models driven periodically between such resets, the reset-averaged steady state is found to be described by a family of generalized Gibbs ensembles (GGE s) characterized by $r$. We also study the spread of particle density of a non-interacting one-dimensional fermionic chain, starting from an initial state where all fermions occupy the left half of the sample, while the right half is empty. When driven by resetting dynamics, the density profile approaches at long times to a nonequilibrium stationary profile that we compute exactly. We suggest concrete experiments that can possibly test our theory.

preprint2020arXiv

Statistics of the Number of Records for Random Walks and Lévy Flights on a ${1D}$ Lattice

We study the statistics of the number of records $R_n$ for a symmetric, $n$-step, discrete jump process on a $1D$ lattice. At a given step, the walker can jump by arbitrary lattice units drawn from a given symmetric probability distribution. This process includes, as a special case, the standard nearest neighbor lattice random walk. We derive explicitly the generating function of the distribution $P(R_n)$ of the number of records, valid for arbitrary discrete jump distributions. As a byproduct, we provide a relatively simple proof of the generalized Sparre Andersen theorem for the survival probability of a random walk on a line, with discrete or continuous jump distributions. For the discrete jump process, we then derive the asymptotic large $n$ behavior of $P(R_n)$ as well as of the average number of records $E(R_n)$. We show that unlike the case of random walks with symmetric and continuous jump distributions where the record statistics is strongly universal (i.e., independent of the jump distribution for all $n$), the record statistics for lattice walks depends on the jump distribution for any fixed $n$. However, in the large $n$ limit, we show that the distribution of the scaled record number $R_n/E(R_n)$ approaches a universal, half-Gaussian form for any discrete jump process. The dependence on the jump distribution enters only through the scale factor $E(R_n)$, which we also compute in the large $n$ limit for arbitrary jump distributions. We present explicit results for a few examples and provide numerical checks of our analytical predictions.

preprint2020arXiv

Stochastic Resetting and Applications

In this Topical Review we consider stochastic processes under resetting, which have attracted a lot of attention in recent years. We begin with the simple example of a diffusive particle whose position is reset randomly in time with a constant rate $r$, which corresponds to Poissonian resetting, to some fixed point (e.g. its initial position). This simple system already exhibits the main features of interest induced by resetting: (i) the system reaches a nontrivial nonequilibrium stationary state (ii) the mean time for the particle to reach a target is finite and has a minimum, optimal, value as a function of the resetting rate $r$. We then generalise to an arbitrary stochastic process (e.g. Lévy flights or fractional Brownian motion) and non-Poissonian resetting (e.g. power-law waiting time distribution for intervals between resetting events). We go on to discuss multiparticle systems as well as extended systems, such as fluctuating interfaces, under resetting. We also consider resetting with memory which implies resetting the process to some randomly selected previous time. Finally we give an overview of recent developments and applications in the field.

preprint2020arXiv

Time between the maximum and the minimum of a stochastic process

We present an exact solution for the probability density function $P(τ=t_{\min}-t_{\max}|T)$ of the time-difference between the minimum and the maximum of a one-dimensional Brownian motion of duration $T$. We then generalise our results to a Brownian bridge, i.e. a periodic Brownian motion of period $T$. We demonstrate that these results can be directly applied to study the position-difference between the minimal and the maximal height of a fluctuating $(1+1)$-dimensional Kardar-Parisi-Zhang interface on a substrate of size $L$, in its stationary state. We show that the Brownian motion result is universal and, asymptotically, holds for any discrete-time random walk with a finite jump variance. We also compute this distribution numerically for Lévy flights and find that it differs from the Brownian motion result.

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

Universal gap statistics for random walks for a class of jump densities

We study the order statistics of a random walk (RW) of $n$ steps whose jumps are distributed according to symmetric Erlang densities $f_p(η)\sim |η|^p \,e^{-|η|}$, parametrized by a non-negative integer $p$. Our main focus is on the statistics of the gaps $d_{k,n}$ between two successive maxima $d_{k,n}=M_{k,n}-M_{k+1,n}$ where $M_{k,n}$ is the $k$-th maximum of the RW between step 1 and step $n$. In the limit of large $n$, we show that the probability density function of the gaps $P_{k,n}(Δ) = \Pr(d_{k,n} = Δ)$ reaches a stationary density $P_{k,n}(Δ) \to p_k(Δ)$. For large $k$, we demonstrate that the typical fluctuations of the gap, for $d_{k,n}= O(1/\sqrt{k})$ (and $n \to \infty$), are described by a non-trivial scaling function that is independent of $k$ and of the jump probability density function $f_p(η)$, thus corroborating our conjecture about the universality of the regime of typical fluctuations (see G. Schehr, S. N. Majumdar, Phys. Rev. Lett. 108, 040601 (2012)). We also investigate the large fluctuations of the gap, for $d_{k,n} = O(1)$ (and $n \to \infty$), and show that these two regimes of typical and large fluctuations of the gaps match smoothly.

preprint2020arXiv

Universal survival probability for a $d$-dimensional run-and-tumble particle

We consider an active run-and-tumble particle (RTP) in $d$ dimensions and compute exactly the probability $S(t)$ that the $x$-component of the position of the RTP does not change sign up to time $t$. When the tumblings occur at a constant rate, we show that $S(t)$ is independent of $d$ for any finite time $t$ (and not just for large $t$), as a consequence of the celebrated Sparre Andersen theorem for discrete-time random walks in one dimension. Moreover, we show that this universal result holds for a much wider class of RTP models in which the speed $v$ of the particle after each tumbling is random, drawn from an arbitrary probability distribution. We further demonstrate, as a consequence, the universality of the record statistics in the RTP problem.

preprint2020arXiv

Velocity and diffusion constant of an active particle in a one dimensional force field

We consider a run an tumble particle with two velocity states $\pm v_0$, in an inhomogeneous force field $f(x)$ in one dimension. We obtain exact formulae for its velocity $V_L$ and diffusion constant $D_L$ for arbitrary periodic $f(x)$ of period $L$. They involve the "active potential" which allows to define a global bias. Upon varying parameters, such as an external force $F$, the dynamics undergoes transitions from non-ergodic trapped states, to various moving states, some with non analyticities in the $V_L$ versus $F$ curve. A random landscape in the presence of a bias leads, for large $L$, to anomalous diffusion $x \sim t^μ$, $μ<1$, or to a phase with a finite velocity that we calculate.

preprint2019arXiv

A First-Order Dynamical Transition in the displacement distribution of a Driven Run-and-Tumble Particle

We study the probability distribution $P(X_N=X,N)$ of the total displacement $X_N$ of an $N$-step run and tumble particle on a line, in presence of a constant nonzero drive $E$. While the central limit theorem predicts a standard Gaussian form for $P(X,N)$ near its peak, we show that for large positive and negative $X$, the distribution exhibits anomalous large deviation forms. For large positive $X$, the associated rate function is nonanalytic at a critical value of the scaled distance from the peak where its first derivative is discontinuous. This signals a first-order dynamical phase transition from a homogeneous `fluid' phase to a `condensed' phase that is dominated by a single large run. A similar first-order transition occurs for negative large fluctuations as well. Numerical simulations are in excellent agreement with our analytical predictions.

preprint2019arXiv

Exact stationary state of a run-and-tumble particle with three internal states in a harmonic trap

We study the motion of a one-dimensional run-and-tumble particle with three discrete internal states in the presence of a harmonic trap of stiffness $μ.$ The three internal states, corresponding to positive, negative and zero velocities respectively, evolve following a jump process with rate $γ$. We compute the stationary position distribution exactly for arbitrary values of $μ$ and $γ$ which turns out to have a finite support on the real line. We show that the distribution undergoes a shape-transition as $β=γ/μ$ is changed. For $β<1,$ the distribution has a double-concave shape and shows algebraic divergences with an exponent $(β-1)$ both at the origin and at the boundaries. For $β>1,$ the position distribution becomes convex, vanishing at the boundaries and with a single, finite, peak at the origin. We also show that for the special case $β=1,$ the distribution shows a logarithmic divergence near the origin while saturating to a constant value at the boundaries.

preprint2017arXiv

Mean perimeter and mean area of the convex hull over planar random walks

We investigate the geometric properties of the convex hull over $n$ successive positions of a planar random walk, with a symmetric continuous jump distribution. We derive the large $n$ asymptotic behavior of the mean perimeter. In addition, we compute the mean area for the particular case of isotropic Gaussian jumps. While the leading terms of these asymptotics are universal, the subleading (correction) terms depend on finer details of the jump distribution and describe a "finite size effect" of discrete-time jump processes, allowing one to accurately compute the mean perimeter and the mean area even for small $n$, as verified by Monte Carlo simulations. This is particularly valuable for applications dealing with discrete-time jumps processes and ranging from the statistical analysis of single-particle tracking experiments in microbiology to home range estimations in ecology.

preprint2016arXiv

Conditioned random walks and interaction-driven condensation

We consider a discrete-time continuous-space random walk under the constraints that the number of returns to the origin (local time) and the total area under the walk are fixed. We first compute the joint probability of an excursion having area $a$ and returning to the origin for the first time after time $τ$. We then show how condensation occurs when the total area constraint is increased: an excursion containing a finite fraction of the area emerges. Finally we show how the phenomena generalises previously studied cases of condensation induced by several constraints and how it is related to interaction-driven condensation which allows us to explain the phenomenon in the framework of large deviation theory.

preprint2016arXiv

Convex Hulls of Multiple Random Walks: A Large-Deviation Study

We study the polygons governing the convex hull of a point set created by the steps of $n$ independent two-dimensional random walkers. Each such walk consists of $T$ discrete time steps, where $x$ and $y$ increments are i.i.d. Gaussian. We analyze area $A$ and perimeter $L$ of the convex hulls. We obtain probability densities for these two quantities over a large range of the support by using a large-deviation approach allowing us to study densities below $10^{-900}$. We find that the densities exhibit a universal scaling behavior as a function of $A/T$ and $L/\sqrt{T}$, respectively. As in the case of one walker ($n=1$), the densities follow Gaussian distributions for $L$ and $\sqrt{A}$, respectively. We also obtained the rate functions for the area and perimeter, rescaled with the scaling behavior of their maximum possible values, and found limiting functions for $T \rightarrow \infty$, revealing that the densities follow the large-deviation principle. These rate functions can be described by a power law for $n \rightarrow \infty$ as found in the $n=1$ case. We also investigated the behavior of the averages as a function of the number of walks $n$ and found good agreement with the predicted behavior.

preprint2016arXiv

Exact distributions of cover times for $N$ independent random walkers in one dimension

We study the probability density function (PDF) of the cover time $t_c$ of a finite interval of size $L$, by $N$ independent one-dimensional Brownian motions, each with diffusion constant $D$. The cover time $t_c$ is the minimum time needed such that each point of the entire interval is visited by at least one of the $N$ walkers. We derive exact results for the full PDF of $t_c$ for arbitrary $N \geq 1$, for both reflecting and periodic boundary conditions. The PDFs depend explicitly on $N$ and on the boundary conditions. In the limit of large $N$, we show that $t_c$ approaches its average value $\langle t_c \rangle \approx L^2/(16\, D \, \ln N)$, with fluctuations vanishing as $1/(\ln N)^2$. We also compute the centered and scaled limiting distributions for large $N$ for both boundary conditions and show that they are given by nontrivial $N$-independent scaling functions.

preprint2016arXiv

Exact extreme value statistics at mixed order transitions

We study extreme value statistics (EVS) for spatially extended models exhibiting mixed order phase transitions (MOT). These are phase transitions which exhibit features common to both first order (discontinuity of the order parameter) and second order (diverging correlation length) transitions. We consider here the truncated inverse distance squared Ising (TIDSI) model which is a prototypical model exhibiting MOT, and study analytically the extreme value statistics of the domain lengths. The lengths of the domains are identically distributed random variables except for the global constraint that their sum equals the total system size $L$. In addition, the number of such domains is also a fluctuating variable, and not fixed. In the paramagnetic phase, we show that the distribution of the largest domain length $l_{\max}$ converges, in the large $L$ limit, to a Gumbel distribution. However, at the critical point (for a certain range of parameters) and in the ferromagnetic phase, we show that the fluctuations of $l_{\max}$ are governed by novel distributions which we compute exactly. Our main analytical results are verified by numerical simulations.

preprint2016arXiv

Exact statistics of record increments of random walks and Lévy flights

We study the statistics of increments in record values in a time series $\{x_0=0,x_1, x_2, \ldots, x_n\}$ generated by the positions of a random walk (discrete time, continuous space) of duration $n$ steps. For arbitrary jump length distribution, including Lévy flights, we show that the distribution of the record increment becomes stationary, i.e., independent of $n$ for large $n$, and compute it explicitly for a wide class of jump distributions. In addition, we compute exactly the probability $Q(n)$ that the record increments decrease monotonically up to step $n$. Remarkably, $Q(n)$ is universal (i..e., independent of the jump distribution) for each $n$, decaying as $Q(n) \sim {\cal A}/\sqrt{n}$ for large $n$, with a universal amplitude ${\cal A} = e/\sqrtπ = 1.53362\ldots$.

preprint2016arXiv

Large deviations for the height in 1D Kardar-Parisi-Zhang growth at late times

We study the atypically large deviations of the height $H \sim {\cal O}(t)$ at the origin at late times in $1+1$-dimensional growth models belonging to the Kardar-Parisi-Zhang (KPZ) universality class. We present exact results for the rate functions for the discrete single step growth model, as well as for the continuum KPZ equation in a droplet geometry. Based on our exact calculation of the rate functions we argue that models in the KPZ class undergo a third order phase transition from a strong coupling to a weak coupling phase, at late times.

preprint2016arXiv

Non-interacting fermions at finite temperature in a $d$-dimensional trap: universal correlations

We study a system of $N$ non-interacting spin-less fermions trapped in a confining potential, in arbitrary dimensions $d$ and arbitrary temperature $T$. The presence of the trap introduces an edge where the average density of fermions vanishes. Far from the edge, near the center of the trap (the so called "bulk regime"), physical properties of the fermions have traditionally been understood using the Local Density Approximation. However, this approximation drastically fails near the edge where the density vanishes. In this paper we show that, even near the edge, novel universal properties emerge, independently of the details of the confining potential. We show that for large $N$, these fermions in a confining trap, in arbitrary dimensions and at finite temperature, form a determinantal point process. As a result, any $n$-point correlation function can be expressed as an $n \times n$ determinant whose entry is called the kernel. Near the edge, we derive the large $N$ scaling form of the kernels. In $d=1$ and $T=0$, this reduces to the so called Airy kernel, that appears in the Gaussian Unitary Ensemble (GUE) of random matrix theory. In $d=1$ and $T>0$ we show a remarkable connection between our kernel and the one appearing in the $1+1$-dimensional Kardar-Parisi-Zhang equation at finite time. Consequently our result provides a finite $T$ generalization of the Tracy-Widom distribution, that describes the fluctuations of the rightmost fermion at $T=0$. In $d>1$ and $T \geq 0$, while the connection to GUE no longer holds, the process is still determinantal whose analysis provides a new class of kernels, generalizing the $1d$ Airy kernel at $T=0$ obtained in random matrix theory. Some of our finite temperature results should be testable in present-day cold atom experiments, most notably our detailed predictions for the temperature dependence of the fluctuations near the edge.

preprint2016arXiv

Number statistics for $β$-ensembles of random matrices: applications to trapped fermions at zero temperature

Let $\mathcal{P}_β^{(V)} (N_{\cal I})$ be the probability that a $N\times N$ $β$-ensemble of random matrices with confining potential $V(x)$ has $N_{\cal I}$ eigenvalues inside an interval ${\cal I}=[a,b]$ of the real line. We introduce a general formalism, based on the Coulomb gas technique and the resolvent method, to compute analytically $\mathcal{P}_β^{(V)} (N_{\cal I})$ for large $N$. We show that this probability scales for large $N$ as $\mathcal{P}_β^{(V)} (N_{\cal I})\approx \exp\left(-βN^2 ψ^{(V)}(N_{\cal I} /N)\right)$, where $β$ is the Dyson index of the ensemble. The rate function $ψ^{(V)}(k_{\cal I})$, independent of $β$, is computed in terms of single integrals that can be easily evaluated numerically. The general formalism is then applied to the classical $β$-Gaussian (${\cal I}=[-L,L]$), $β$-Wishart (${\cal I}=[1,L]$) and $β$-Cauchy (${\cal I}=[-L,L]$) ensembles. Expanding the rate function around its minimum, we find that generically the number variance ${\rm Var}(N_{\cal I})$ exhibits a non-monotonic behavior as a function of the size of the interval, with a maximum that can be precisely characterized. These analytical results, corroborated by numerical simulations, provide the full counting statistics of many systems where random matrix models apply. In particular, we present results for the full counting statistics of zero temperature one-dimensional spinless fermions in a harmonic trap.

preprint2016arXiv

Record statistics for random walk bridges

We investigate the statistics of records in a random sequence $\{x_B(0)=0,x_B(1),\cdots, x_B(n)=x_B(0)=0\}$ of $n$ time steps. The sequence $x_B(k)$'s represents the position at step $k$ of a random walk `bridge' of $n$ steps that starts and ends at the origin. At each step, the increment of the position is a random jump drawn from a specified symmetric distribution. We study the statistics of records and record ages for such a bridge sequence, for different jump distributions. In absence of the bridge condition, i.e., for a free random walk sequence, the statistics of the number and ages of records exhibits a `strong' universality for all $n$, i.e., they are completely independent of the jump distribution as long as the distribution is continuous. We show that the presence of the bridge constraint destroys this strong `all $n$' universality. Nevertheless a `weaker' universality still remains for large $n$, where we show that the record statistics depends on the jump distributions only through a single parameter $0<μ\le 2$, known as the Lévy index of the walk, but are insensitive to the other details of the jump distribution. We derive the most general results (for arbitrary jump distributions) wherever possible and also present two exactly solvable cases. We present numerical simulations that verify our analytical results.

preprint2016arXiv

Universal ground state properties of free fermions in a $d$-dimensional trap

The ground state properties of $N$ spinless free fermions in a $d$-dimensional confining potential are studied. We find that any $n$-point correlation function has a simple determinantal structure that allows us to compute several properties exactly for large $N$. We show that the average density has a finite support with an edge, and near this edge the density exhibits a universal (valid for a wide class of potentials) scaling behavior for large $N$. The associated edge scaling function is computed exactly and generalizes to any $d$ the edge electron gas result of Kohn and Mattsson in $d=3$ [Phys. Rev. Lett. 81, 3487 (1998)]. In addition, we calculate the kernel (that characterizes any $n$-point correlation function) for large $N$ and show that, when appropriately scaled, it depends only on dimension $d$, but has otherwise universal scaling forms, at the edges. The edge kernel, for higher $d$, generalizes the Airy kernel in one dimension, well known from random matrix theory.

preprint2015arXiv

Convex Hulls of Random Walks: Large-Deviation Properties

We study the convex hull of the set of points visited by a two-dimensional random walker of T discrete time steps. Two natural observables that characterize the convex hull in two dimensions are its perimeter L and area A. While the mean perimeter <L> and the mean area <A> have been studied before, analytically and numerically, and exact results are known for large T (Brownian motion limit), little is known about the full distributions P(A) and P(L). In this paper, we provide numerical results for these distributions. We use a sophisticated large-deviation approach that allows us to study the distributions over a larger range of the support, where the probabilities P(A) and P(L) are as small as 10^{-300}. We analyze (open) random walks as well as (closed) Brownian bridges on the two-dimensional discrete grid as well as in the two-dimensional plane. The resulting distributions exhibit, for large T, a universal scaling behavior (independent of the details of the jump distributions) as a function of A/T and L/\sqrt{T}, respectively. We are also able to obtain the rate function, describing rare events at the tails of these distributions, via a numerical extrapolation scheme and find a linear and square dependence as a function of the rescaled perimeter and the rescaled area, respectively.

preprint2015arXiv

Dynamical transition in the temporal relaxation of stochastic processes under resetting

A stochastic process, when subject to resetting to its initial condition at a constant rate, generically reaches a non-equilibrium steady state. We study analytically how the steady state is approached in time and find an unusual relaxation mechanism in these systems. We show that as time progresses, an inner core region around the resetting point reaches the steady state, while the region outside the core is still transient. The boundaries of the core region grow with time as power laws at late times. Alternatively, at a fixed spatial point, the system undergoes a dynamical transition from the transient to the steady state at a characteristic space dependent timescale $t^*(x)$. We calculate analytically in several examples the large deviation function associated with this spatio-temporal fluctuation and show that generically it has a second order discontinuity at a pair of critical points characterizing the edges of the inner core. Our results are verified in the numerical simulations of several models, such as simple diffusion and fluctuating one-dimensional interfaces.

preprint2015arXiv

Effective Langevin equations for constrained stochastic processes

We propose a novel stochastic method to exactly generate Brownian paths conditioned to start at an initial point and end at a given final point during a fixed time $t_{f}$. These paths are weighted with a probability given by the overdamped Langevin dynamics. We show how these paths can be exactly generated by a local stochastic differential equation. The method is illustrated on the generation of Brownian bridges, Brownian meanders, Brownian excursions and constrained Ornstein-Uehlenbeck processes. In addition, we show how to solve this equation in the case of a general force acting on the particle. As an example, we show how to generate constrained path joining the two minima of a double-well. Our method allows to generate statistically independent paths, and is computationally very efficient.

preprint2015arXiv

Exact gap statistics for the random average process on a ring with a tracer

We study statistics of the gaps in Random Average Process (RAP) on a ring with particles hopping symmetrically, except one tracer particle which could be driven. These particles hop either to the left or to the right by a random fraction $η$ of the space available till next particle in the respective directions. The random fraction $η\in [0,~1)$ is chosen from a distribution $R(η)$. For non-driven tracer, when $R(η)$ satisfies a necessary and sufficient condition, the stationary joint distribution of the gaps between successive particles takes an universal form that is factorized except for a global constraint. Some interesting explicit forms of $R(η)$ are found which satisfy this condition. In case of driven tracer, the system reaches a current-carrying steady state where such factorization does not hold. Analytical progress has been made in the thermodynamic limit, where we computed the single site mass distribution inside the bulk. We have also computed the two point gap-gap correlation exactly in that limit. Numerical simulations support our analytical results.

preprint2015arXiv

Extreme value statistics of correlated random variables

Extreme value statistics (EVS) concerns the study of the statistics of the maximum or the minimum of a set of random variables. This is an important problem for any time-series and has applications in climate, finance, sports, all the way to physics of disordered systems where one is interested in the statistics of the ground state energy. While the EVS of uncorrelated variables are well understood, little is known for strongly correlated random variables. Only recently this subject has gained much importance both in statistical physics and in probability theory. In this note, we will first review the classical EVS for uncorrelated variables and discuss few examples of correlated variables where analytical progress can be made.

preprint2015arXiv

Finite temperature free fermions and the Kardar-Parisi-Zhang equation at finite time

We consider the system of $N$ one-dimensional free fermions confined by a harmonic well $V(x) = mω^2 {x^2}/{2}$ at finite inverse temperature $β= 1/T$. The average density of fermions $ρ_N(x,T)$ at position $x$ is derived. For $N \gg 1$ and $β\sim {\cal O}(1/N)$, $ρ_N(x,T)$ is given by a scaling function interpolating between a Gaussian at high temperature, for $β\ll 1/N$, and the Wigner semi-circle law at low temperature, for $β\gg N^{-1}$. In the latter regime, we unveil a scaling limit, for $β{\hbar ω}= b N^{-1/3}$, where the fluctuations close to the edge of the support, at $x \sim \pm \sqrt{2\hbar N/(mω)}$, are described by a limiting kernel $K^{\rm ff}_b(s,s')$ that depends continuously on $b$ and is a generalization of the Airy kernel, found in the Gaussian Unitary Ensemble of random matrices. Remarkably, exactly the same kernel $K^{\rm ff}_b(s,s')$ arises in the exact solution of the Kardar-Parisi-Zhang (KPZ) equation in 1+1 dimensions at finite time $t$, with the correspondence $t= b^3$.

preprint2015arXiv

On certain functionals of the maximum of Brownian motion and their applications

We consider a Brownian motion (BM) $x(τ)$ and its maximal value $x_{\max} = \max_{0 \leq τ\leq t} x(τ)$ on a fixed time interval $[0,t]$. We study functionals of the maximum of the BM, of the form ${\cal O}_{\max}(t)=\int_0^t\, V(x_{\max} - x(τ)) {\rm d} τ$ where $V(x)$ can be any arbitrary function and develop various analytical tools to compute their statistical properties. These tools rely in particular on (i) a "counting paths" method and (ii) a path-integral approach. In particular, we focus on the case where $V(x) = δ(x-r)$, with $r$ a real parameter, which is relevant to study the density of near-extreme values of the BM (the so called density of states), $ρ(r,t)$, which is the local time of the BM spent at given distance $r$ from the maximum. We also provide a thorough analysis of the family of functionals ${T}_α(t)=\int_0^t (x_{\max} - x(τ))^α\, {\rm d}τ$, corresponding to $V(x) = x^α$, with $α$ real. As $α$ is varied, $T_α(t)$ interpolates between different interesting observables. For instance, for $α=1$, $T_{α= 1}(t)$ is a random variable of the "area", or "Airy", type while for $α=-1/2$ it corresponds to the maximum time spent by a ballistic particle through a Brownian random potential. On the other hand, for $α= -1$, it corresponds to the cost of the optimal algorithm to find the maximum of a discrete random walk, proposed by Odlyzko. We revisit here, using tools of theoretical physics, the statistical properties of this algorithm which had been studied before using probabilistic methods. Finally, we extend our methods to constrained BM, including in particular the Brownian bridge, i.e., the Brownian motion starting and ending at the origin.

preprint2015arXiv

On the Gap and Time Interval between the First Two Maxima of Long Continuous Time Random Walks

We consider a one-dimensional continuous time random walk (CTRW) on a fixed time interval $T$ where at each time step the walker waits a random time $τ$, before performing a jump drawn from a symmetric continuous probability distribution function (PDF) $f(η)$, of Lévy index $0 < μ\leq 2$. Our study includes the case where the waiting time PDF $Ψ(τ)$ has a power law tail, $Ψ(τ) \propto τ^{-1 - γ}$, with $0< γ< 1$, such that the average time between two consecutive jumps is infinite. The random motion is sub-diffusive if $γ< μ/2$ (and super-diffusive if $γ> μ/2$). We investigate the joint PDF of the gap $g$ between the first two highest positions of the CTRW and the time $t$ separating these two maxima. We show that this PDF reaches a stationary limiting joint distribution $p(g,t)$ in the limit of long CTRW, $T \to \infty$. Our exact analytical results show a very rich behavior of this joint PDF in the $(γ, μ)$ plane, which we study in great detail. Our main results are verified by numerical simulations. This work provides a non trivial extension to CTRWs of the recent study in the discrete time setting by Majumdar et al. (J. Stat. Mech. P09013, 2014).

preprint2015arXiv

Random walk with random resetting to the maximum

We study analytically a simple random walk model on a one-dimensional lattice, where at each time step the walker resets to the maximum of the already visited positions (to the rightmost visited site) with a probability $r$, and with probability $(1-r)$, it undergoes symmetric random walk, i.e., it hops to one of its neighboring sites, with equal probability $(1-r)/2$. For $r=0$, it reduces to a standard random walk whose typical distance grows as $\sqrt{n}$ for large $n$. In presence of a nonzero resetting rate $0<r\le 1$, we find that both the average maximum and the average position grow ballistically for large $n$, with a common speed $v(r)$. Moreover, the fluctuations around their respective averages grow diffusively, again with the same diffusion coefficient $D(r)$. We compute $v(r)$ and $D(r)$ explicitly. We also show that the probability distribution of the difference between the maximum and the location of the walker, becomes stationary as $n\to \infty$. However, the approach to this stationary distribution is accompanied by a dynamical phase transition, characterized by a weakly singular large deviation function. We also show that $r=0$ is a special `critical' point, for which the growth laws are different from the $r\to 0$ case and we calculate the exact crossover functions that interpolate between the critical $(r=0)$ and the off-critical $(r\to 0)$ behavior for finite but large $n$.

preprint2015arXiv

Spatial Extent of Branching Brownian Motion

We study the one dimensional branching Brownian motion starting at the origin and investigate the correlation between the rightmost ($X_{\max}\geq 0$) and leftmost ($X_{\min} \leq 0$) visited sites up to time $t$. At each time step the existing particles in the system either diffuse (with diffusion constant $D$), die (with rate $a$) or split into two particles (with rate $b$). We focus on the regime $b \leq a$ where these two extreme values $X_{\max}$ and $X_{\min}$ are strongly correlated. We show that at large time $t$, the joint probability distribution function (PDF) of the two extreme points becomes stationary $P(X,Y,t \to \infty) \to p(X,Y)$. Our exact results for $p(X,Y)$ demonstrate that the correlation between $X_{\max}$ and $X_{\min}$ is nonzero, even in the stationary state. From this joint PDF, we compute exactly the stationary PDF $p(ζ)$ of the (dimensionless) span $ζ= {(X_{\max} - X_{\min})}/{\sqrt{D/b}}$, which is the distance between the rightmost and leftmost visited sites. This span distribution is characterized by a linear behavior ${p}(ζ) \sim \frac{1}{2} \left(1 + Δ\right) ζ$ for small spans, with $Δ= \left(\frac{a}{b} -1\right)$. In the critical case ($Δ= 0$) this distribution has a non-trivial power law tail ${p}(ζ) \sim 8 π\sqrt{3} /ζ^3$ for large spans. On the other hand, in the subcritical case ($Δ> 0$), we show that the span distribution decays exponentially as ${p}(ζ) \sim (A^2/2) ζ\exp \left(- \sqrtΔ~ζ\right)$ for large spans, where $A$ is a non-trivial function of $Δ$ which we compute exactly. We show that these asymptotic behaviors carry the signatures of the correlation between $X_{\max}$ and $X_{\min}$. Finally we verify our results via direct Monte Carlo simulations.

preprint2015arXiv

Statistics of the longest interval in renewal processes

We consider renewal processes where events, which can for instance be the zero crossings of a stochastic process, occur at random epochs of time. The intervals of time between events, $τ_{1},τ_{2},...$, are independent and identically distributed (i.i.d.) random variables with a common density $ρ(τ)$. Fixing the total observation time to $t$ induces a global constraint on the sum of these random intervals, which accordingly become interdependent. Here we focus on the largest interval among such a sequence on the fixed time interval $(0,t)$. Depending on how the last interval is treated, we consider three different situations, indexed by $α=$ I, II and III. We investigate the distribution of the longest interval $\ell^α_{\max}(t)$ and the probability $Q^α(t)$ that the last interval is the longest one. We show that if $ρ(τ)$ decays faster than $1/τ^2$ for large $τ$, then the full statistics of $\ell^α_{\max}(t)$ is given, in the large $t$ limit, by the standard theory of extreme value statistics for i.i.d. random variables, showing in particular that the global constraint on the intervals $τ_i$ does not play any role at large times in this case. However, if $ρ(τ)$ exhibits heavy tails, $ρ(τ)\simτ^{-1-θ}$ for large $τ$, with index $0 <θ<1$, we show that the fluctuations of $\ell^α_{\max}(t)/t$ are governed, in the large $t$ limit, by a stationary universal distribution which depends on both $θ$ and $α$, which we compute exactly. On the other hand, $Q^α(t)$ is generically different from its counterpart for i.i.d. variables (both for narrow or heavy tailed distributions $ρ(τ)$). In particular, in the case $0<θ<1$, the large $t$ behaviour of $Q^α(t)$ gives rise to universal constants (depending also on both $θ$ and $α$) which we compute exactly.

preprint2015arXiv

Survival probability of a Brownian motion in a planar wedge of arbitrary angle

We study the survival probability and the first-passage time distribution for a Brownian motion in a planar wedge with infinite absorbing edges. We generalize existing results obtained for wedge angles of the form $π/n$ with $n$ a positive integer to arbitrary angles, which in particular cover the case of obtuse angles. We give explicit and simple expressions of the survival probability and the first-passage time distribution in which the difference between an arbitrary angle and a submultiple of $π$ is contained in three additional terms. As an application, we obtain the short time development of the survival probability in a wedge of arbitrary angle.

preprint2015arXiv

The average number of distinct sites visited by a random walker on random graphs

We study the linear large $n$ behavior of the average number of distinct sites $S(n)$ visited by a random walker after $n$ steps on a large random graph. An expression for the graph topology dependent prefactor $B$ in $S(n) = Bn$ is proposed. We use generating function techniques to relate this prefactor to the graph adjacency matrix and then devise message-passing equations to calculate its value. Numerical simulations are performed to evaluate the agreement between the message passing predictions and random walk simulations on random graphs. Scaling with system size and average graph connectivity are also analysed.

preprint2015arXiv

Top eigenvalue of a random matrix: large deviations and third order phase transition

We study the fluctuations of the largest eigenvalue $λ_{\max}$ of $N \times N$ random matrices in the limit of large $N$. The main focus is on Gaussian $β$-ensembles, including in particular the Gaussian orthogonal ($β=1$), unitary ($β=2$) and symplectic ($β= 4$) ensembles. The probability density function (PDF) of $λ_{\max}$ consists, for large $N$, of a central part described by Tracy-Widom distributions flanked, on both sides, by two large deviations tails. While the central part characterizes the typical fluctuations of $λ_{\max}$ -- of order ${\cal O}(N^{-2/3})$ --, the large deviations tails are instead associated to extremely rare fluctuations -- of order ${\cal O}(1)$. Here we review some recent developments in the theory of these extremely rare events using a Coulomb gas approach. We discuss in particular the third-order phase transition which separates the left tail from the right tail, a transition akin to the so-called Gross-Witten-Wadia phase transition found in 2-d lattice quantum chromodynamics. We also discuss the occurrence of similar third-order transitions in various physical problems, including non-intersecting Brownian motions, conductance fluctuations in mesoscopic physics and entanglement in a bipartite system.

preprint2014arXiv

Branching Brownian Motion Conditioned on Particle Numbers

We study analytically the order and gap statistics of particles at time $t$ for the one dimensional branching Brownian motion, conditioned to have a fixed number of particles at $t$. The dynamics of the process proceeds in continuous time where at each time step, every particle in the system either diffuses (with diffusion constant $D$), dies (with rate $d$) or splits into two independent particles (with rate $b$). We derive exact results for the probability distribution function of $g_k(t) = x_k(t) - x_{k+1}(t)$, the distance between successive particles, conditioned on the event that there are exactly $n$ particles in the system at a given time $t$. We show that at large times these conditional distributions become stationary $P(g_k, t \to \infty|n) = p(g_k|n)$. We show that they are characterised by an exponential tail $p(g_k|n) \sim \exp[-\sqrt{\frac{|b - d|}{2 D}} ~g_k]$ for large gaps in the subcritical ($b < d$) and supercritical ($b > d$) phases, and a power law tail $p(g_k) \sim 8\left(\frac{D}{b}\right){g_k}^{-3}$ at the critical point ($b = d$), independently of $n$ and $k$. Some of these results for the critical case were announced in a recent letter [K. Ramola, S. N. Majumdar and G. Schehr, Phys. Rev. Lett. 112, 210602 (2014)].

preprint2014arXiv

Condensation transition in joint large deviations of linear statistics

Real space condensation is known to occur in stochastic models of mass transport in the regime in which the globally conserved mass density is greater than a critical value. It has been shown within models with factorised stationary states that the condensation can be understood in terms of sums of independent and identically distributed random variables: these exhibit condensation when they are conditioned to a large deviation of their sum. It is well understood that the condensation, whereby one of the random variables contributes a finite fraction to the sum, occurs only if the underlying probability distribution (modulo exponential) is heavy-tailed, i.e. decaying slower than exponential. Here we study a similar phenomenon in which condensation is exhibited for non-heavy-tailed distributions, provided random variables are additionally conditioned on a large deviation of certain linear statistics. We provide a detailed theoretical analysis explaining the phenomenon, which is supported by Monte Carlo simulations (for the case where the additional constraint is the sample variance) and demonstrated in several physical systems. Our results suggest that the condensation is a generic phenomenon that pertains to both typical and rare events.

preprint2014arXiv

Constraint driven condensation in large fluctuations of linear statistics

Condensation is the phenomenon whereby one of a sum of random variables contributes a finite fraction to the sum. It is manifested as an aggregation phenomenon in diverse physical systems such as coalescence in granular media, jamming in traffic and gelation in networks. We show here that the same condensation scenario, which normally happens only if the underlying probability distribution has tails heavier than exponential, can occur for light-tailed distributions in the presence of additional constraints. We demonstrate this phenomenon on the sample variance, whose probability distribution conditioned on the particular value of the sample mean undergoes a phase transition. The transition is manifested by a change in behavior of the large deviation rate function.

preprint2014arXiv

Diffusion with resetting in arbitrary spatial dimension

We consider diffusion in arbitrary spatial dimension d with the addition of a resetting process wherein the diffusive particle stochastically resets to a fixed position at a constant rate $r$. We compute the non-equilibrium stationary state which exhibits non-Gaussian behaviour. We then consider the presence of an absorbing target centred at the origin and compute the survival probability and mean time to absorption of the diffusive particle by the target. The mean absorption time is finite and has a minimum value at an optimal resetting rate $r^*$ which depends on dimension. Finally we consider the problem of a finite density of diffusive particles, each resetting to its own initial position. While the typical survival probability of the target at the origin decays exponentially with time regardless of spatial dimension, the average survival probability decays asymptotically as $\exp -A (\log t)^d$ where $A$ is a constant. We explain these findings using an interpretation as a renewal process and arguments invoking extreme value statistics.

preprint2014arXiv

First order transition for the optimal search time of Lévy flights with resetting

We study analytically an intermittent search process in one dimension. There is an immobile target at the origin and a searcher undergoes a discrete time jump process starting at $x_0\geq0$, where successive jumps are drawn independently from an arbitrary jump distribution $f(η)$. In addition, with a probability $0\leq r \leq1$ the position of the searcher is reset to its initial position $x_0$. The efficiency of the search strategy is characterized by the mean time to find the target, i.e., the mean first passage time (MFPT) to the origin. For arbitrary jump distribution $f(η)$, initial position $x_0$ and resetting probability $r$, we compute analytically the MFPT. For the heavy-tailed Lévy stable jump distribution characterized by the Lévy index $0<μ< 2$, we show that, for any given $x_0$, the MFPT has a global minimum in the $(μ,r)$ plane at $(μ^*(x_0),r^*(x_0))$. We find a remarkable first-order phase transition as $x_0$ crosses a critical value $x_0^*$ at which the optimal parameters change discontinuously. Our analytical results are in good agreement with numerical simulations.

preprint2014arXiv

Instanton Approach to Large $N$ Harish-Chandra-Itzykson-Zuber Integrals

We reconsider the large $N$ asymptotics of Harish-Chandra-Itzykson-Zuber integrals. We provide, using Dyson's Brownian motion and the method of instantons, an alternative, transparent derivation of the Matytsin formalism for the unitary case. Our method is easily generalized to the orthogonal and symplectic ensembles. We obtain an explicit solution of Matytsin's equations in the case of Wigner matrices, as well as a general expansion method in the dilute limit, when the spectrum of eigenvalues spreads over very wide regions.

preprint2014arXiv

Long-range correlations in a locally driven exclusion process

We show that the presence of a driven bond in an otherwise diffusive lattice gas with simple exclusion interaction results in long-range density-density correlation in its stationary state. In dimensions $d>1$ we show that in the thermodynamic limit this correlation decays as $C(r,s)\sim (r^2+s^2)^{-d}$ at large distances $r$ and $s$ away from the drive with $|r-s|>>1$. This is derived using an electrostatic analogy whereby $C(r,s)$ is expressed as the potential due to a configuration of electrostatic charges distributed in $2d$-dimension. At bulk density $ρ=1/2$ we show that the potential is that of a localized quadrupolar charge. At other densities the same is correct in leading order in the strength of the drive and is argued numerically to be valid at higher orders.

preprint2014arXiv

Non-local response in a lattice gas under a shear drive

In equilibrium, the effect of a spatially localised perturbation is typically confined around the perturbed region. Quite contrary to this, in a non-equilibrium stationary state often the entire system is affected. This appears to be a generic feature of non-equilibrium. We study such non-local response in the stationary state of a lattice gas with a shear drive at the boundary which keeps the system out of equilibrium. We show that a perturbation in the form of a localised blockage at the boundary, induces algebraically decaying density and current profile. In two examples, non-interacting particles and particles with simple exclusion, we analytically derive the power-law tail of the profiles.

preprint2014arXiv

On the Gap and Time Interval between the First Two Maxima of Long Random Walks

In the context of order statistics of discrete time random walks (RW), we investigate the statistics of the gap, $G_n$, and the number of time steps, $L_n$, between the two highest positions of a Markovian one-dimensional random walker, starting from $x_0 = 0$, after $n$ time steps (taking the $x$-axis vertical). The jumps $η_i = x_i - x_{i-1}$ are independent and identically distributed random variables drawn from a symmetric probability distribution function (PDF), $f(η)$, the Fourier transform of which has the small $k$ behavior $1 - \hat f(k) \propto |k|^μ$, with $0 < μ\leq 2$. For $μ=2$, the variance of the jump distribution is finite and the RW (properly scaled) converges to a Brownian motion. For $0<μ<2$, the RW is a Lévy flight of index $μ$. We show that the joint PDF of $G_n$ and $L_n$ converges to a well defined stationary bi-variate distribution $p(g,l)$ as the RW duration $n$ goes to infinity. We present a thorough analytical study of the limiting joint distribution $p(g,l)$, as well as of its associated marginals $p_{\rm gap}(g)$ and $p_{\rm time}(l)$, revealing a rich variety of behaviors depending on the tail of $f(η)$ (from slow decreasing algebraic tail to fast decreasing super-exponential tail). We also address the problem for a random bridge where the RW starts and ends at the origin after $n$ time steps. We show that in the large $n$ limit, the PDF of $G_n$ and $L_n$ converges to the {\it same} stationary distribution $p(g,l)$ as in the case of the free-end RW. Finally, we present a numerical check of our analytical predictions. Some of these results were announced in a recent letter [S. N. Majumdar, Ph. Mounaix, G. Schehr, Phys. Rev. Lett. {\bf 111}, 070601 (2013)].

preprint2014arXiv

Phase transitions and edge scaling of number variance in Gaussian random matrices

We consider $N\times N$ Gaussian random matrices, whose average density of eigenvalues has the Wigner semi-circle form over $[-\sqrt{2},\sqrt{2}]$. For such matrices, using a Coulomb gas technique, we compute the large $N$ behavior of the probability $\mathcal{P}_{\scriptscriptstyle N,L}(N_L)$ that $N_L$ eigenvalues lie within the box $[-L,L]$. This probability scales as $\mathcal{P}_{\scriptscriptstyle N,L}(N_L=κ_L N)\approx\exp\left(-β N^2 ψ_L(κ_L)\right)$, where $β$ is the Dyson index of the ensemble and $ψ_L(κ_L)$ is a $β$-independent rate function that we compute exactly. We identify three regimes as $L$ is varied: (i) $\, N^{-1}\ll L<\sqrt{2}$ (bulk), (ii) $\ L\sim\sqrt{2}$ on a scale of $\mathcal{O}(N^{-{2}/{3}})$ (edge) and (iii) $\ L > \sqrt{2}$ (tail). We find a dramatic non-monotonic behavior of the number variance $V_N(L)$ as a function of $L$: after a logarithmic growth $\propto \ln (N L)$ in the bulk (when $L \sim {\cal O}(1/N)$), $V_N(L)$ decreases abruptly as $L$ approaches the edge of the semi-circle before it decays as a stretched exponential for $L > \sqrt{2}$. This "drop-off" of $V_N(L)$ at the edge is described by a scaling function $\tilde V_β$ which smoothly interpolates between the bulk (i) and the tail (iii). For $β= 2$ we compute $\tilde V_2$ explicitly in terms of the Airy kernel. These analytical results, verified by numerical simulations, directly provide for $β=2$ the full statistics of particle-number fluctuations at zero temperature of 1d spinless fermions in a harmonic trap.

preprint2014arXiv

Random matrices and entanglement entropy of trapped Fermi gases

We exploit and clarify the use of random matrix theory for the calculation of the entanglement entropy of free Fermi gases. We apply this method to obtain analytic predictions for Renyi entanglement entropies of a one-dimensional gas trapped by a harmonic potential in all the relevant scaling regimes. We confirm our findings with accurate numerical calculations obtained by means of an ingenious discretisation of the reduced correlation matrix.

preprint2014arXiv

Statistics of the first passage time of Brownian motion conditioned by maximum value or area

We derive the moments of the first passage time for Brownian motion conditioned by either the maximum value or the area swept out by the motion. These quantities are the natural counterparts to the moments of the maximum value and area of Brownian excursions of fixed duration, which we also derive for completeness within the same mathematical framework. Various applications are indicated.

preprint2014arXiv

Universal Order and Gap Statistics of Critical Branching Brownian Motion

We study the order statistics of one dimensional branching Brownian motion in which particles either diffuse (with diffusion constant $D$), die (with rate $d$) or split into two particles (with rate $b$). At the critical point $b=d$ which we focus on, we show that, at large time $t$, the particles are collectively bunched together. We find indeed that there are two length scales in the system: (i) the diffusive length scale $\sim \sqrt{Dt}$ which controls the collective fluctuations of the whole bunch and (ii) the length scale of the gap between the bunched particles $\sim \sqrt{D/b}$. We compute the probability distribution function $P(g_k,t|n)$ of the $k$th gap $g_k = x_k - x_{k+1}$ between the $k$th and $(k+1)$th particles given that the system contains exactly $n>k$ particles at time $t$. We show that at large $t$, it converges to a stationary distribution $P(g_k,t\to \infty|n) = p(g_k|n)$ with an algebraic tail $p(g_k|n) \sim 8(D/b) g_k^{-3}$, for $g_k \gg 1$, independent of $k$ and $n$. We verify our predictions with Monte Carlo simulations.

preprint2014arXiv

Universal statistics of longest lasting records of random walks and Lévy flights

We study the record statistics of random walks after $n$ steps, $x_0, x_1,\ldots, x_n$, with arbitrary symmetric and continuous distribution $p(η)$ of the jumps $η_i = x_i - x_{i-1}$. We consider the age of the records, i.e. the time up to which a record survives. Depending on how the age of the current last record is defined, we propose three distinct sequences of ages (indexed by $α$ = I, II, III) associated to a given sequence of records. We then focus on the longest lasting record, which is the longest element among this sequence of ages. To characterize the statistics of these longest lasting records, we compute: (i) the probability that the record of the longest age is broken at step $n$, denoted by $Q^α(n)$, which we call the probability of record breaking and: (ii) the duration of the longest lasting record, $\ell_{\max}^α(n)$. We show that both $Q^α(n)$ and the full statistics of $\ell_{\max}^α(n)$ are universal, i.e. independent of the jump distribution $p(η)$. We compute exactly the large $n$ asymptotic behaviors of $Q^α(n)$ as well as $\langle \ell_{\max}^α(n)\rangle$ (when it exists) and show that each case gives rise to a different universal constant associated to random walks (including Lévy flights). While two of them appeared before in the excursion theory of Brownian motion, for which we provide here a simpler derivation, the third case gives rise to a non-trivial new constant $C^{\rm III} = 0.241749 \ldots$ associated to the records of random walks. Other observables characterizing the ages of the records, exhibiting an interesting universal behavior, are also discussed.

preprint2014arXiv

Wigner time-delay distribution in chaotic cavities and freezing transition

Using the joint distribution for proper time-delays of a chaotic cavity derived by Brouwer, Frahm & Beenakker [Phys. Rev. Lett. {\bf 78}, 4737 (1997)], we obtain, in the limit of large number of channels $N$, the large deviation function for the distribution of the Wigner time-delay (the sum of proper times) by a Coulomb gas method. We show that the existence of a power law tail originates from narrow resonance contributions, related to a (second order) freezing transition in the Coulomb gas.

preprint2014arXiv

Winding statistics of a Brownian particle on a ring

We consider a Brownian particle moving on a ring. We study the probability distributions of the total number of turns and the net number of counter-clockwise turns the particle makes till time t. Using a method based on the renewal properties of Brownian walker, we find exact analytical expressions of these distributions. This method serves as an alternative to the standard path integral techniques which are not always easily adaptable for certain observables. For large t, we show that these distributions have Gaussian scaling forms. We also compute large deviation functions associated to these distributions characterizing atypically large fluctuations. We provide numerical simulations in support of our analytical results.

preprint2013arXiv

Effect of Partial Absorption on Diffusion with Resetting

The effect of partial absorption on a diffusive particle which stochastically resets its position with a finite rate $r$ is considered. The particle is absorbed by a target at the origin with absorption `velocity' $a$; as the velocity $a$ approaches $\infty$ the absorption property of the target approaches that of a perfectly-absorbing target. The effect of partial absorption on first-passage time problems is studied, in particular, it is shown that the mean time to absorption (MTA) is increased by an additive term proportional to $1/a$. The results are extended to multiparticle systems where independent searchers, initially uniformly distributed with a given density, look for a single immobile target. It is found that the average survival probability $P^{av}$ is modified by a multiplicative factor which is a function of $1/a$, whereas the decay rate of the typical survival probability $P^{typ}$ is decreased by an additive term proportional to $1/a$.

preprint2013arXiv

Exact distributions of the number of distinct and common sites visited by N independent random walkers

We study the number of distinct sites S_N(t) and common sites W_N(t) visited by N independent one dimensional random walkers, all starting at the origin, after t time steps. We show that these two random variables can be mapped onto extreme value quantities associated to N independent random walkers. Using this mapping, we compute exactly their probability distributions P_N^d(S,t) and P_N^d(W,t) for any value of N in the limit of large time t, where the random walkers can be described by Brownian motions. In the large N limit one finds that S_N(t)/\sqrt{t} \propto 2 \sqrt{\log N} + \widetilde{s}/(2 \sqrt{\log N}) and W_N(t)/\sqrt{t} \propto \widetilde{w}/N where \widetilde{s} and \widetilde{w} are random variables whose probability density functions (pdfs) are computed exactly and are found to be non trivial. We verify our results through direct numerical simulations.

preprint2013arXiv

Exact record and order statistics of random walks via first-passage ideas

While records and order statistics of independent and identically distributed (i.i.d.) random variables X_1, ..., X_N are fully understood, much less is known for strongly correlated random variables, which is often the situation encountered in statistical physics. Recently, it was shown, in a series of works, that one-dimensional random walk (RW) is an interesting laboratory where the influence of strong correlations on records and order statistics can be studied in detail. We review here recent exact results which have been obtained for these questions about RW, using techniques borrowed from the study of first-passage problems. We also present a brief review of the well known (and not so well known) results for records and order statistics of i.i.d. variables.

preprint2013arXiv

Exact Statistics of the Gap and Time Interval Between the First Two Maxima of Random Walks

We investigate the statistics of the gap, G_n, between the two rightmost positions of a Markovian one-dimensional random walker (RW) after n time steps and of the duration, L_n, which separates the occurrence of these two extremal positions. The distribution of the jumps η_i's of the RW, f(η), is symmetric and its Fourier transform has the small k behavior 1-\hat{f}(k)\sim| k|^μwith 0 < μ\leq 2. We compute the joint probability density function (pdf) P_n(g,l) of G_n and L_n and show that, when n \to \infty, it approaches a limiting pdf p(g,l). The corresponding marginal pdf of the gap, p_{\rm gap}(g), is found to behave like p_{\rm gap}(g) \sim g^{-1 - μ} for g \gg 1 and 0<μ< 2. We show that the limiting marginal distribution of L_n, p_{\rm time}(l), has an algebraic tail p_{\rm time}(l) \sim l^{-γ(μ)} for l \gg 1 with γ(1<μ\leq 2) = 1 + 1/μ, and γ(0<μ<1) = 2. For l, g \gg 1 with fixed l g^{-μ}, p(g,l) takes the scaling form p(g,l) \sim g^{-1-2μ} \tilde p_μ(l g^{-μ}) where \tilde p_μ(y) is a (μ-dependent) scaling function. We also present numerical simulations which verify our analytic results.

preprint2013arXiv

Fluctuating interfaces subject to stochastic resetting

We study one-dimensional fluctuating interfaces of length $L$ where the interface stochastically resets to a fixed initial profile at a constant rate $r$. For finite $r$ in the limit $L \to \infty$, the system settles into a nonequilibrium stationary state with non-Gaussian interface fluctuations, which we characterize analytically for the Kardar-Parisi-Zhang and Edwards-Wilkinson universality class. Our results are corroborated by numerical simulations. We also discuss the generality of our results for a fluctuating interface in a generic universality class.

preprint2013arXiv

Index Distribution of Cauchy Random Matrices

Using a Coulomb gas technique, we compute analytically the probability $\mathcal{P}_β^{(C)}(N_+,N)$ that a large $N\times N$ Cauchy random matrix has $N_+$ positive eigenvalues, where $N_+$ is called the index of the ensemble. We show that this probability scales for large $N$ as $\mathcal{P}_β^{(C)}(N_+,N)\approx \exp\left[-βN^2 ψ_C(N_+/N)\right]$, where $β$ is the Dyson index of the ensemble. The rate function $ψ_C(κ)$ is computed in terms of single integrals that are easily evaluated numerically and amenable to an asymptotic analysis. We find that the rate function, around its minimum at $κ=1/2$, has a quadratic behavior modulated by a logarithmic singularity. As a consequence, the variance of the index scales for large $N$ as $\mathrm{Var}(N_+)\sim σ_C\ln N$, where $σ_C=2/(βπ^2)$ is twice as large as the corresponding prefactor in the Gaussian and Wishart cases. The analytical results are checked by numerical simulations and against an exact finite $N$ formula which, for $β=2$, can be derived using orthogonal polynomials.

preprint2013arXiv

Near-extreme statistics of Brownian motion

We study the statistics of near-extreme events of Brownian motion (BM) on the time interval [0,t]. We focus on the density of states (DOS) near the maximum ρ(r,t) which is the amount of time spent by the process at a distance r from the maximum. We develop a path integral approach to study functionals of the maximum of BM, which allows us to study the full probability density function (PDF) of ρ(r,t) and obtain an explicit expression for the moments, \langle [ρ(r,t)]^k \rangle, for arbitrary integer k. We also study near-extremes of constrained BM, like the Brownian bridge. Finally we also present numerical simulations to check our analytical results.

preprint2013arXiv

Persistence and First-Passage Properties in Non-equilibrium Systems

In this review we discuss the persistence and the related first-passage properties in extended many-body nonequilibrium systems. Starting with simple systems with one or few degrees of freedom, such as random walk and random acceleration problems, we progressively discuss the persistence properties in systems with many degrees of freedom. These systems include spins models undergoing phase ordering dynamics, diffusion equation, fluctuating interfaces etc. Persistence properties are nontrivial in these systems as the effective underlying stochastic process is non-Markovian. Several exact and approximate methods have been developed to compute the persistence of such non-Markov processes over the last two decades, as reviewed in this article. We also discuss various generalisations of the local site persistence probability. Persistence in systems with quenched disorder is discussed briefly. Although the main emphasis of this review is on the theoretical developments on persistence, we briefly touch upon various experimental systems as well.

preprint2013arXiv

Record-breaking statistics for random walks in the presence of measurement error and noise

We address the question of distance record-setting by a random walker in the presence of measurement error, $δ$, and additive noise, $γ$ and show that the mean number of (upper) records up to $n$ steps still grows universally as $< R_n> \sim n^{1/2}$ for large $n$ for all jump distributions, including Lévy flights, and for all $δ$ and $γ$. In contrast to the universal growth exponent of 1/2, the pace of record setting, measured by the pre-factor of $n^{1/2}$, depends on $δ$ and $γ$. In the absence of noise ($γ=0$), the pre-factor $S(δ)$ is evaluated explicitly for arbitrary jump distributions and it decreases monotonically with increasing $δ$ whereas, in case of perfect measurement $(δ=0)$, the corresponding pre-factor $T(γ)$ increases with $γ$. Our analytical results are supported by extensive numerical simulations and qualitatively similar results are found in two and three dimensions.

preprint2013arXiv

Sampling fractional Brownian motion in presence of absorption: a Markov Chain method

We study fractional Brownian motion (fBm) characterized by the Hurst exponent H. Using a Monte Carlo sampling technique, we are able to numerically generate fBm processes with an absorbing boundary at the origin at discrete times for a large number of 10^7 time steps even for small values like H=1/4. The results are compatible with previous analytical results that the distribution of (rescaled) endpoints y follow a power law P(y) y^ϕwith ϕ=(1-H)/H, even for small values of H. Furthermore, for the case H=0.5 we also study analytically the finite-length corrections to the first order, namely a plateau of P(y) for y->0 which decreases with increasing process length. These corrections are compatible with the numerical results.

preprint2013arXiv

Spatial extent of an outbreak in animal epidemics

Characterizing the spatial extent of epidemics at the outbreak stage is key to controlling the evolution of the disease. At the outbreak, the number of infected individuals is typically small, so that fluctuations around their average are important: then, it is commonly assumed that the susceptible-infected-recovered (SIR) mechanism can be described by a stochastic birth-death process of Galton-Watson type. The displacements of the infected individuals can be modelled by resorting to Brownian motion, which is applicable when long-range movements and complex network interactions can be safely neglected, as in case of animal epidemics. In this context, the spatial extent of an epidemic can be assessed by computing the convex hull enclosing the infected individuals at a given time. We derive the exact evolution equations for the mean perimeter and the mean area of the convex hull, and compare them with Monte Carlo simulations.

preprint2012arXiv

Bose-Einstein Condensation of a Gaussian Random Field in the Thermodynamic Limit

We derive the criterion for the Bose-Einstein condensation (BEC) of a Gaussian field $ϕ$ (real or complex) in the thermodynamic limit. The field is characterized by its covariance function and the control parameter is the intensity $u=\|ϕ\|_2^2/V$, where $V$ is the volume of the box containing the field. We show that for any dimension $d$ (including $d=1$), there is a class of covariance functions for which $ϕ$ exhibits a BEC as $u$ is increased through a critical value $u_c$. In this case, we investigate the probability distribution of the part of $u$ contained in the condensate. We show that depending on the parameters characterizing the covariance function and the dimension $d$, there can be two distinct types of condensate: a Gaussian distributed "normal" condensate with fluctuations scaling as $1/\sqrt{V}$, and a non Gaussian distributed "anomalous" condensate. A detailed analysis of the anomalous condensate is performed for a one-dimensional system ($d=1$). Extending this one-dimensional analysis to exactly the point of transition between normal and anomalous condensations, we find that the condensate at the transition point is still Gaussian distributed but with anomalously large fluctuations scaling as $\sqrt{\ln(L)/L}$, where $L$ is the system length. The conditional spectral density of $ϕ$, knowing $u$, is given for all the regimes (with and without BEC).

preprint2012arXiv

Can a Lamb Reach a Haven Before Being Eaten by Diffusing Lions?

We study the survival of a single diffusing lamb on the positive half line in the presence of N diffusing lions that all start at the same position L to the right of the lamb and a haven at x=0. If the lamb reaches this haven before meeting any lion, the lamb survives. We investigate the survival probability of the lamb, S_N(x,L), as a function of N and the respective initial positions of the lamb and the lions, x and L. We determine S_N(x,L) analytically for the special cases of N=1 and N--->oo. For large but finite N, we determine the unusual asymptotic form whose leading behavior is S_N(z)\simN^{-z^2}, with z=x/L. Simulations of the capture process very slowly converge to this asymptotic prediction as N reaches 10^{500}.

preprint2012arXiv

Invariant $β$-Wishart ensembles, crossover densities and asymptotic corrections to the Marchenko-Pastur law

We construct a diffusive matrix model for the $β$-Wishart (or Laguerre) ensemble for general continuous $β\in [0,2]$, which preserves invariance under the orthogonal/unitary group transformation. Scaling the Dyson index $β$ with the largest size $M$ of the data matrix as $β=2c/M$ (with $c$ a fixed positive constant), we obtain a family of spectral densities parametrized by $c$. As $c$ is varied, this density interpolates continuously between the Mar\vcenko-Pastur ($c\to \infty$ limit) and the Gamma law ($c\to 0$ limit). Analyzing the full Stieltjes transform (resolvent) equation, we obtain as a byproduct the correction to the Mar\vcenko-Pastur density in the bulk up to order 1/M for all $β$ and up to order $1/M^2$ for the particular cases $β=1,2$.

preprint2012arXiv

Large deviations of the top eigenvalue of large Cauchy random matrices

We compute analytically the probability density function (pdf) of the largest eigenvalue $λ_{\max}$ in rotationally invariant Cauchy ensembles of $N\times N$ matrices. We consider unitary ($β= 2$), orthogonal ($β=1$) and symplectic ($β=4$) ensembles of such heavy-tailed random matrices. We show that a central non-Gaussian regime for $λ_{\max} \sim \mathcal{O}(N)$ is flanked by large deviation tails on both sides which we compute here exactly for any value of $β$. By matching these tails with the central regime, we obtain the exact leading asymptotic behaviors of the pdf in the central regime, which generalizes the Tracy-Widom distribution known for Gaussian ensembles, both at small and large arguments and for any $β$. Our analytical results are confirmed by numerical simulations.

preprint2012arXiv

Number of Common Sites Visited by N Random Walkers

We compute analytically the mean number of common sites, W_N(t), visited by N independent random walkers each of length t and all starting at the origin at t=0 in d dimensions. We show that in the (N-d) plane, there are three distinct regimes for the asymptotic large t growth of W_N(t). These three regimes are separated by two critical lines d=2 and d=d_c(N)=2N/(N-1) in the (N-d) plane. For d<2, W_N(t)\sim t^{d/2} for large t (the N dependence is only in the prefactor). For 2<d<d_c(N), W_N(t)\sim t^ν where the exponent ν= N-d(N-1)/2 varies with N and d. For d>d_c(N), W_N(t) approaches a constant as t\to \infty. Exactly at the critical dimensions there are logaritmic corrections: for d=2, we get W_N(t)\sim t/[\ln t]^N, while for d=d_c(N), W_N(t)\sim \ln t for large t. Our analytical predictions are verified in numerical simulations.

preprint2012arXiv

Optimal diffusive search: nonequilibrium resetting versus equilibrium dynamics

We study first-passage time problems for a diffusive particle with stochastic resetting with a finite rate $r$. The optimal search time is compared quantitatively with that of an effective equilibrium Langevin process with the same stationary distribution. It is shown that the intermittent, nonequilibrium strategy with non-vanishing resetting rate is more efficient than the equilibrium dynamics. Our results are extended to multiparticle systems where a team of independent searchers, initially uniformly distributed with a given density, looks for a single immobile target. Both the average and the typical survival probability of the target are smaller in the case of nonequilibrium dynamics.

preprint2012arXiv

Record statistics and persistence for a random walk with a drift

We study the statistics of records of a one-dimensional random walk of n steps, starting from the origin, and in presence of a constant bias c. At each time-step the walker makes a random jump of length ηdrawn from a continuous distribution f(η) which is symmetric around a constant drift c. We focus in particular on the case were f(η) is a symmetric stable law with a Lévy index 0 < μ\leq 2. The record statistics depends crucially on the persistence probability which, as we show here, exhibits different behaviors depending on the sign of c and the value of the parameter μ. Hence, in the limit of a large number of steps n, the record statistics is sensitive to these parameters (c and μ) of the jump distribution. We compute the asymptotic mean record number <R_n> after n steps as well as its full distribution P(R,n). We also compute the statistics of the ages of the longest and the shortest lasting record. Our exact computations show the existence of five distinct regions in the (c, 0 < μ\leq 2) strip where these quantities display qualitatively different behaviors. We also present numerical simulation results that verify our analytical predictions.

preprint2012arXiv

Record Statistics for Multiple Random Walks

We study the statistics of the number of records R_{n,N} for N identical and independent symmetric discrete-time random walks of n steps in one dimension, all starting at the origin at step 0. At each time step, each walker jumps by a random length drawn independently from a symmetric and continuous distribution. We consider two cases: (I) when the variance σ^2 of the jump distribution is finite and (II) when σ^2 is divergent as in the case of Lévy flights with index 0 < μ< 2. In both cases we find that the mean record number <R_{n,N}> grows universally as \sim α_N \sqrt{n} for large n, but with a very different behavior of the amplitude α_N for N > 1 in the two cases. We find that for large N, α_N \approx 2 \sqrt{\log N} independently of σ^2 in case I. In contrast, in case II, the amplitude approaches to an N-independent constant for large N, α_N \approx 4/\sqrtπ, independently of 0<μ<2. For finite σ^2 we argue, and this is confirmed by our numerical simulations, that the full distribution of (R_{n,N}/\sqrt{n} - 2 \sqrt{\log N}) \sqrt{\log N} converges to a Gumbel law as n \to \infty and N \to \infty. In case II, our numerical simulations indicate that the distribution of R_{n,N}/\sqrt{n} converges, for n \to \infty and N \to \infty, to a universal nontrivial distribution, independently of μ. We discuss the applications of our results to the study of the record statistics of 366 daily stock prices from the Standard & Poors 500 index.

preprint2012arXiv

Reunion probability of N vicious walkers: typical and large fluctuations for large N

We consider three different models of N non-intersecting Brownian motions on a line segment [0,L] with absorbing (model A), periodic (model B) and reflecting (model C) boundary conditions. In these three cases we study a properly normalized reunion probability, which, in model A, can also be interpreted as the maximal height of N non-intersecting Brownian excursions on the unit time interval. We provide a detailed derivation of the exact formula for these reunion probabilities for finite N using a Fermionic path integral technique. We then analyse the asymptotic behavior of this reunion probability for large N using two complementary techniques: (i) a saddle point analysis of the underlying Coulomb gas and (ii) orthogonal polynomial method. These two methods are complementary in the sense that they work in two different regimes, respectively for L\ll O(\sqrt{N}) and L\geq O(\sqrt{N}). A striking feature of the large N limit of the reunion probability in the three models is that it exhibits a third-order phase transition when the system size L crosses a critical value L=L_c(N)\sim \sqrt{N}. This transition is akin to the Douglas-Kazakov transition in two-dimensional continuum Yang-Mills theory. While the central part of the reunion probability, for L \sim L_c(N), is described in terms of the Tracy-Widom distributions (associated to GOE and GUE depending on the model), the emphasis of the present study is on the large deviations of these reunion probabilities, both in the right [L \gg L_c(N)] and the left [L \ll L_c(N)] tails. In particular, for model B, we find that the matching between the different regimes corresponding to typical L \sim L_c(N) and atypical fluctuations in the right tail L \gg L_c(N) is rather unconventional, compared to the usual behavior found for the distribution of the largest eigenvalue of GUE random matrices.

preprint2011arXiv

A simple derivation of the Tracy-Widom distribution of the maximal eigenvalue of a Gaussian unitary random matrix

In this paper, we first briefly review some recent results on the distribution of the maximal eigenvalue of a $(N\times N)$ random matrix drawn from Gaussian ensembles. Next we focus on the Gaussian Unitary Ensemble (GUE) and by suitably adapting a method of orthogonal polynomials developed by Gross and Matytsin in the context of Yang-Mills theory in two dimensions, we provide a rather simple derivation of the Tracy-Widom law for GUE. Our derivation is based on the elementary asymptotic scaling analysis of a pair of coupled nonlinear recursion relations. As an added bonus, this method also allows us to compute the precise subleading terms describing the right large deviation tail of the maximal eigenvalue distribution. In the Yang-Mills language, these subleading terms correspond to non-perturbative (in $1/N$ expansion) corrections to the two-dimensional partition function in the so called `weak' coupling regime.

preprint2011arXiv

Diffusion with Optimal Resetting

We consider the mean time to absorption by an absorbing target of a diffusive particle with the addition of a process whereby the particle is reset to its initial position with rate $r$. We consider several generalisations of the model of M. R. Evans and S. N. Majumdar (2011), Diffusion with stochastic resetting, Phys. Rev. Lett. 106, 160601: (i) a space dependent resetting rate $r(x)$ ii) resetting to a random position $z$ drawn from a resetting distribution ${\cal P}(z)$ iii) a spatial distribution for the absorbing target $P_T(x)$. As an example of (i) we show that the introduction of a non-resetting window around the initial position can reduce the mean time to absorption provided that the initial position is sufficiently far from the target. We address the problem of optimal resetting, that is, minimising the mean time to absorption for a given target distribution. For an exponentially decaying target distribution centred at the origin we show that a transition in the optimal resetting distribution occurs as the target distribution narrows.

preprint2011arXiv

Diffusion with Stochastic Resetting

We study simple diffusion where a particle stochastically resets to its initial position at a constant rate r. A finite resetting rate leads to a nonequilibrium stationary state with non-Gaussian fluctuations for the particle position. We also show that the mean time to find a stationary target by a diffusive searcher is finite and has a minimum value at an optimal resetting rate r^*. Resetting also alters fundamentally the late time decay of the survival probability of a stationary target when there are multiple searchers: while the typical survival probability decays exponentially with time, the average decays as a power law with an exponent depending continuously on the density of searchers.

preprint2011arXiv

Large deviations of the maximal eigenvalue of random matrices

We present detailed computations of the 'at least finite' terms (three dominant orders) of the free energy in a one-cut matrix model with a hard edge a, in beta-ensembles, with any polynomial potential. beta is a positive number, so not restricted to the standard values beta = 1 (hermitian matrices), beta = 1/2 (symmetric matrices), beta = 2 (quaternionic self-dual matrices). This model allows to study the statistic of the maximum eigenvalue of random matrices. We compute the large deviation function to the left of the expected maximum. We specialize our results to the gaussian beta-ensembles and check them numerically. Our method is based on general results and procedures already developed in the literature to solve the Pastur equations (also called "loop equations"). It allows to compute the left tail of the analog of Tracy-Widom laws for any beta, including the constant term.

preprint2011arXiv

Long-range steady state density profiles induced by localized drive

We show that the presence of a localized drive in an otherwise diffusive system results in steady-state density and current profiles that decay algebraically to their global average value, away from the drive in two or higher dimensions. An analogy to an electrostatic problem is established, whereby the density profile induced by a driving bond maps onto the electrostatic potential due to an electric dipole located along the bond. The dipole strength is proportional to the drive, and is determined self-consistently by solving the electrostatic problem. The profile resulting from a localized configuration of more than one driving bond can be straightforwardly determined by the superposition principle of electrostatics. This picture is shown to hold even in the presence of exclusion interaction between particles.

preprint2011arXiv

Non-intersecting Brownian walkers and Yang-Mills theory on the sphere

We study a system of N non-intersecting Brownian motions on a line segment [0,L] with periodic, absorbing and reflecting boundary conditions. We show that the normalized reunion probabilities of these Brownian motions in the three models can be mapped to the partition function of two-dimensional continuum Yang-Mills theory on a sphere respectively with gauge groups U(N), Sp(2N) and SO(2N). Consequently, we show that in each of these Brownian motion models, as one varies the system size L, a third order phase transition occurs at a critical value L=L_c(N)\sim \sqrt{N} in the large N limit. Close to the critical point, the reunion probability, properly centered and scaled, is identical to the Tracy-Widom distribution describing the probability distribution of the largest eigenvalue of a random matrix. For the periodic case we obtain the Tracy-Widom distribution corresponding to the GUE random matrices, while for the absorbing and reflecting cases we get the Tracy-Widom distribution corresponding to GOE random matrices. In the absorbing case, the reunion probability is also identified as the maximal height of N non-intersecting Brownian excursions ("watermelons" with a wall) whose distribution in the asymptotic scaling limit is then described by GOE Tracy-Widom law. In addition, large deviation formulas for the maximum height are also computed.

preprint2011arXiv

Number of relevant directions in Principal Component Analysis and Wishart random matrices

We compute analytically, for large $N$, the probability $\mathcal{P}(N_+,N)$ that a $N\times N$ Wishart random matrix has $N_+$ eigenvalues exceeding a threshold $Nζ$, including its large deviation tails. This probability plays a benchmark role when performing the Principal Component Analysis of a large empirical dataset. We find that $\mathcal{P}(N_+,N)\approx\exp(-βN^2 ψ_ζ(N_+/N))$, where $β$ is the Dyson index of the ensemble and $ψ_ζ(κ)$ is a rate function that we compute explicitly in the full range $0\leq κ\leq 1$ and for any $ζ$. The rate function $ψ_ζ(κ)$ displays a quadratic behavior modulated by a logarithmic singularity close to its minimum $κ^\star(ζ)$. This is shown to be a consequence of a phase transition in an associated Coulomb gas problem. The variance $Δ(N)$ of the number of relevant components is also shown to grow universally (independent of $ζ)$ as $Δ(N)\sim (βπ^2)^{-1}\ln N$ for large $N$.

preprint2011arXiv

Perturbation Theory for Fractional Brownian Motion in Presence of Absorbing Boundaries

Fractional Brownian motion is a Gaussian process x(t) with zero mean and two-time correlations <x(t)x(s)> ~ t^{2H} + s^{2H} - |t-s|^{2H}, where H, with 0<H<1 is called the Hurst exponent. For H = 1/2, x(t) is a Brownian motion, while for H unequal 1/2, x(t) is a non-Markovian process. Here we study x(t) in presence of an absorbing boundary at the origin and focus on the probability density P(x,t) for the process to arrive at x at time t, starting near the origin at time 0, given that it has never crossed the origin. It has a scaling form P(x,t) ~ R(x/t^H)/t^H. Our objective is to compute the scaling function R(y), which up to now was only known for the Markov case H=1/2. We develop a systematic perturbation theory around this limit, setting H = 1/2 + epsilon, to calculate the scaling function R(y) to first order in epsilon. We find that R(y) behaves as R(y) ~ y^phi as y -> 0 (near the absorbing boundary), while R(y) ~ y^gamma exp(-y^2/2) as y -> oo, with phi = 1 - 4 epsilon + O(epsilon^2) and gamma = 1 - 2 epsilon + O(epsilon^2). Our epsilon-expansion result confirms the scaling relation phi = (1-H)/H proposed in Ref. [28]. We verify our findings via numerical simulations for H = 2/3. The tools developed here are versatile, powerful, and adaptable to different situations.

preprint2011arXiv

Phase transitions in the distribution of the Andreev conductance of superconductor-metal junctions with multiple transverse modes

We compute analytically the full distribution of Andreev conductance $G_{\mathrm{NS}}$ of a metal-superconductor interface with a large number $N_c$ of transverse modes, using a random matrix approach. The probability distribution $\mathcal{P}(G_{\mathrm{NS}},N_c)$ in the limit of large $N_c$ displays a Gaussian behavior near the average value $<G_{\mathrm{NS}}>= (2-\sqrt{2}) N_c$ and asymmetric power-law tails in the two limits of very small and very large $G_{\mathrm{NS}}$. In addition, we find a novel third regime sandwiched between the central Gaussian peak and the power law tail for large $G_{\mathrm{NS}}$. Weakly non-analytic points separate these four regimes---these are shown to be consequences of three phase transitions in an associated Coulomb gas problem.

preprint2011arXiv

The convex hull for a random acceleration process in two dimensions

We compute exactly the mean perimeter <L(T)> and the mean area <A(T)> of the convex hull of a random acceleration process of duration T in two dimensions. We use an exact mapping that relates, via Cauchy's formulae, the computation of the perimeter and the area of the convex hull of an arbitrary two dimensional stochastic process [x(t); y(t)] to the computation of the extreme value statistics of the associated one dimensional component process x(t). The latter can be computed exactly for the one dimensional random acceleration process even though the process in non-Markovian. Physically, our results are relevant in describing theaverage shape of a semi-flexible ideal polymer chain in two dimensions.

preprint2011arXiv

Universal Order Statistics of Random Walks

We study analytically the order statistics of a time series generated by the successive positions of a symmetric random walk of n steps with step lengths of finite variance σ^2. We show that the statistics of the gap d_{k,n}=M_{k,n} -M_{k+1,n} between the k-th and the (k+1)-th maximum of the time series becomes stationary, i.e, independent of n as n\to \infty and exhibits a rich, universal behavior. The mean stationary gap (in units of σ) exhibits a universal algebraic decay for large k, <d_{k,\infty}>/σ\sim 1/\sqrt{2πk}, independent of the details of the jump distribution. Moreover, the probability density (pdf) of the stationary gap exhibits scaling, Proba.(d_{k,\infty}=δ)\simeq (\sqrt{k}/σ) P(δ\sqrt{k}/σ), in the scaling regime when δ\sim <d_{k,\infty}>\simeq σ/\sqrt{2πk}. The scaling function P(x) is universal and has an unexpected power law tail, P(x) \sim x^{-4} for large x. For δ\gg <d_{k,\infty}> the scaling breaks down and the pdf gets cut-off in a nonuniversal way. Consequently, the moments of the gap exhibit an unusual multi-scaling behavior.

preprint2010arXiv

Area distribution and the average shape of a Lévy bridge

We consider a one dimensional Lévy bridge x_B of length n and index 0 < α< 2, i.e. a Lévy random walk constrained to start and end at the origin after n time steps, x_B(0) = x_B(n)=0. We compute the distribution P_B(A,n) of the area A = \sum_{m=1}^n x_B(m) under such a Lévy bridge and show that, for large n, it has the scaling form P_B(A,n) \sim n^{-1-1/α} F_α(A/n^{1+1/α}), with the asymptotic behavior F_α(Y) \sim Y^{-2(1+α)} for large Y. For α=1, we obtain an explicit expression of F_1(Y) in terms of elementary functions. We also compute the average profile < \tilde x_B (m) > at time m of a Lévy bridge with fixed area A. For large n and large m and A, one finds the scaling form < \tilde x_B(m) > = n^{1/α} H_α({m}/{n},{A}/{n^{1+1/α}}), where at variance with Brownian bridge, H_α(X,Y) is a non trivial function of the rescaled time m/n and rescaled area Y = A/n^{1+1/α}. Our analytical results are verified by numerical simulations.

preprint2010arXiv

Extreme Eigenvalues of Wishart Matrices: Application to Entangled Bipartite System

We discuss an application of the random matrix theory in the context of estimating the bipartite entanglement of a quantum system. We discuss how the Wishart ensemble (the earliest studied random matrix ensemble) appears in this quantum problem. The eigenvalues of the reduced density matrix of one of the subsystems have similar statistical properties as those of the Wishart matrices, except that their {\em trace is constrained to be unity}. We focus here on the smallest eigenvalue which serves as an important measure of entanglement between the two subsystems. In the hard edge case (when the two subsystems have equal sizes) one can fully characterize the probability distribution of the minimum eigenvalue for real, complex and quaternion matrices of all sizes. In particular, we discuss the important finite size effect due to the {\em fixed trace constraint}.

preprint2010arXiv

How many eigenvalues of a Gaussian random matrix are positive?

We study the probability distribution of the index ${\mathcal N}_+$, i.e., the number of positive eigenvalues of an $N\times N$ Gaussian random matrix. We show analytically that, for large $N$ and large $\mathcal{N}_+$ with the fraction $0\le c=\mathcal{N}_+/N\le 1$ of positive eigenvalues fixed, the index distribution $\mathcal{P}({\mathcal N}_+=cN,N)\sim\exp[-βN^2 Φ(c)]$ where $β$ is the Dyson index characterizing the Gaussian ensemble. The associated large deviation rate function $Φ(c)$ is computed explicitly for all $0\leq c \leq 1$. It is independent of $β$ and displays a quadratic form modulated by a logarithmic singularity around $c=1/2$. As a consequence, the distribution of the index has a Gaussian form near the peak, but with a variance $Δ(N)$ of index fluctuations growing as $Δ(N)\sim \log N/βπ^2$ for large $N$. For $β=2$, this result is independently confirmed against an exact finite $N$ formula, yielding $Δ(N)= \log N/2π^2 +C+\mathcal{O}(N^{-1})$ for large $N$, where the constant $C$ has the nontrivial value $C=(γ+1+3\log 2)/2π^2\simeq 0.185248...$ and $γ=0.5772...$ is the Euler constant. We also determine for large $N$ the probability that the interval $[ζ_1,ζ_2]$ is free of eigenvalues. Part of these results have been announced in a recent letter [\textit{Phys. Rev. Lett.} {\bf 103}, 220603 (2009)].

preprint2010arXiv

Maximum Distance Between the Leader and the Laggard for Three Brownian Walkers

We consider three independent Brownian walkers moving on a line. The process terminates when the left-most walker (the `Leader') meets either of the other two walkers. For arbitrary values of the diffusion constants D_1 (the Leader), D_2 and D_3 of the three walkers, we compute the probability distribution P(m|y_2,y_3) of the maximum distance m between the Leader and the current right-most particle (the `Laggard') during the process, where y_2 and y_3 are the initial distances between the leader and the other two walkers. The result has, for large m, the form P(m|y_2,y_3) \sim A(y_2,y_3) m^{-δ}, where δ= (2π-θ)/(π-θ) and θ= cos^{-1}(D_1/\sqrt{(D_1+D_2)(D_1+D_3)}. The amplitude A(y_2,y_3) is also determined exactly.

preprint2010arXiv

Maximum of N Independent Brownian Walkers till the First Exit From the Half Space

We consider the one-dimensional target search process that involves an immobile target located at the origin and $N$ searchers performing independent Brownian motions starting at the initial positions $\vec x = (x_1,x_2,..., x_N)$ all on the positive half space. The process stops when the target is first found by one of the searchers. We compute the probability distribution of the maximum distance $m$ visited by the searchers till the stopping time and show that it has a power law tail: $P_N(m|\vec x)\sim B_N (x_1x_2... x_N)/m^{N+1}$ for large $m$. Thus all moments of $m$ up to the order $(N-1)$ are finite, while the higher moments diverge. The prefactor $B_N$ increases with $N$ faster than exponentially. Our solution gives the exit probability of a set of $N$ particles from a box $[0,L]$ through the left boundary. Incidentally, it also provides an exact solution of the Laplace's equation in an $N$-dimensional hypercube with some prescribed boundary conditions. The analytical results are in excellent agreement with Monte Carlo simulations.

preprint2010arXiv

Probability distributions of Linear Statistics in Chaotic Cavities and associated phase transitions

We establish large deviation formulas for linear statistics on the $N$ transmission eigenvalues $\{T_i\}$ of a chaotic cavity, in the framework of Random Matrix Theory. Given any linear statistics of interest $A=\sum_{i=1}^N a(T_i)$, the probability distribution $\mathcal{P}_A(A,N)$ of $A$ generically satisfies the large deviation formula $\lim_{N\to\infty}[-2\log\mathcal{P}_A(Nx,N)/βN^2]=Ψ_A(x)$, where $Ψ_A(x)$ is a rate function that we compute explicitly in many cases (conductance, shot noise, moments) and $β$ corresponds to different symmetry classes. Using these large deviation expressions, it is possible to recover easily known results and to produce new formulas, such as a closed form expression for $v(n)=\lim_{N\to\infty}\mathrm{var}(\mathcal{T}_n)$ (where $\mathcal{T}_n=\sum_{i}T_i^n$) for arbitrary integer $n$. The universal limit $v^\star=\lim_{n\to\infty} v(n)=1/2πβ$ is also computed exactly. The distributions display a central Gaussian region flanked on both sides by non-Gaussian tails. At the junction of the two regimes, weakly non-analytical points appear, a direct consequence of phase transitions in an associated Coulomb gas problem. Numerical checks are also provided, which are in full agreement with our asymptotic results in both real and Laplace space even for moderately small $N$. Part of the results have been announced in [P. Vivo, S.N. Majumdar and O. Bohigas, {\it Phys. Rev. Lett.} {\bf 101}, 216809 (2008)].

preprint2010arXiv

Statistics of layered zigzags: a two-dimensional generalization of TASEP

A novel discrete growth model in 2+1 dimensions is presented in three equivalent formulations: i) directed motion of zigzags on a cylinder, ii) interacting interlaced TASEP layers, and iii) growing heap over 2D substrate with a restricted minimal local height gradient. We demonstrate that the coarse-grained behavior of this model is described by the two-dimensional Kardar-Parisi-Zhang equation. The coefficients of different terms in this hydrodynamic equation can be derived from the steady state flow-density curve, the so called `fundamental' diagram. A conjecture concerning the analytical form of this flow-density curve is presented and is verified numerically.

preprint2010arXiv

Switching and growth for microbial populations in catastrophic responsive environments

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

preprint2010arXiv

Universal First-passage Properties of Discrete-time Random Walks and Levy Flights on a Line: Statistics of the Global Maximum and Records

In these lecture notes I will discuss the universal first-passage properties of a simple correlated discrete-time sequence {x_0=0, x_1,x_2.... x_n} up to n steps where x_i represents the position at step i of a random walker hopping on a continuous line by drawing independently, at each time step, a random jump length from an arbitrary symmetric and continuous distribution (it includes, e.g., the Levy flights). I will focus on the statistics of two extreme observables associated with the sequence: (i) its global maximum and the time step at which the maximum occurs and (ii) the number of records in the sequence and their ages. I will demonstrate how the universal statistics of these observables emerge as a consequence of Pollaczek-Spitzer formula and the associated Sparre Andersen theorem.

preprint2009arXiv

Capture of particles undergoing discrete random walks

It is shown that particles undergoing discrete-time jumps in 3D, starting at a distance r0 from the center of an adsorbing sphere of radius R, are captured with probability (R - c sigma)/r0 for r0 much greater than R, where c is related to the Fourier transform of the scaled jump distribution and sigma is the distribution's root-mean square jump length. For particles starting on the surface of the sphere, the asymptotic survival probability is non-zero (in contrast to the case of Brownian diffusion) and has a universal behavior sigma/(R sqrt(6)) depending only upon sigma/R. These results have applications to computer simulations of reaction and aggregation.

preprint2009arXiv

Condensation of the roots of real random polynomials on the real axis

We introduce a family of real random polynomials of degree n whose coefficients a_k are symmetric independent Gaussian variables with variance <a_k^2> = e^{-k^α}, indexed by a real α\geq 0. We compute exactly the mean number of real roots <N_n> for large n. As αis varied, one finds three different phases. First, for 0 \leq α< 1, one finds that <N_n> \sim (\frac{2}π) \log{n}. For 1 < α< 2, there is an intermediate phase where < N_n > grows algebraically with a continuously varying exponent, < N_n > \sim \frac{2}π \sqrt{\frac{α-1}α} n^{α/2}. And finally for α> 2, one finds a third phase where <N_n> \sim n. This family of real random polynomials thus exhibits a condensation of their roots on the real line in the sense that, for large n, a finite fraction of their roots <N_n>/n are real. This condensation occurs via a localization of the real roots around the values \pm \exp{[\fracα{2}(k+{1/2})^{α-1} ]}, 1 \ll k \leq n.

preprint2009arXiv

Convex Hull of N Planar Brownian Motions: Exact Results and an Application to Ecology

We compute exactly the mean perimeter and area of the convex hull of N independent planar Brownian paths each of duration T, both for open and closed paths. We show that the mean perimeter < L_N > = α_N, \sqrt{T} and the mean area <A_N> = β_N T for all T. The prefactors α_N and β_N, computed exactly for all N, increase very slowly (logarithmically) with increasing N. This slow growth is a consequence of extreme value statistics and has interesting implication in ecological context in estimating the home range of a herd of animals with population size N.

preprint2009arXiv

Hitting probability for anomalous diffusion processes

We present the universal features of the hitting probability $Q(x,L)$, the probability that a generic stochastic process starting at $x$ and evolving in a box $[0,L]$ hits the upper boundary $L$ before hitting the lower boundary at 0. For a generic self-affine process (describing, for instance, the polymer translocation through a nanopore) we show that $Q(x,L)=Q(x/L)$ and the scaling function $Q(z)\sim z^ϕ$ as $z\to 0$ with $ϕ=θ/H$ where $H$ and $θ$ are respectively the Hurst exponent and the persistence exponent of the process. This result is verified in several exact calculations including when the process represents the position of a particle diffusing in a disordered potential. We also provide numerical supports for our analytical results.

preprint2009arXiv

Non-intersecting Brownian Interfaces and Wishart Random Matrices

We study a system of $N$ non-intersecting $(1+1)$-dimensional fluctuating elastic interfaces (`vicious bridges') at thermal equilibrium, each subject to periodic boundary condition in the longitudinal direction and in presence of a substrate that induces an external confining potential for each interface. We show that, for a large system and with an appropriate choice of the external confining potential, the joint distribution of the heights of the $N$ non-intersecting interfaces at a fixed point on the substrate can be mapped to the joint distribution of the eigenvalues of a Wishart matrix of size $N$ with complex entries (Dyson index $β=2$), thus providing a physical realization of the Wishart matrix. Exploiting this analogy to random matrix, we calculate analytically (i) the average density of states of the interfaces (ii) the height distribution of the uppermost and lowermost interfaces (extrema) and (iii) the asymptotic (large $N$) distribution of the center of mass of the interfaces. In the last case, we show that the probability density of the center of mass has an essential singularity around its peak which is shown to be a direct consequence of a phase transition in an associated Coulomb gas problem.

preprint2009arXiv

Phase Transitions in the Distribution of Bipartite Entanglement of a Random Pure State

Using a Coulomb gas method, we compute analytically the probability distribution of the Renyi entropies (a standard measure of entanglement) for a random pure state of a large bipartite quantum system. We show that, for any order q>1 of the Renyi entropy, there are two critical values at which the entropy's probability distribution changes shape. These critical points correspond to two different transitions in the corresponding charge density of the Coulomb gas: the disappearance of an integrable singularity at the origin and the detachement of a single-charge drop from the continuum sea of all the other charges. These transitions respectively control the left and right tails of the entropy's probability distribution, as verified also by Monte Carlo numerical simulations of the Coulomb gas equilibrium dynamics.

preprint2009arXiv

Random Convex Hulls and Extreme Value Statistics

In this paper we study the statistical properties of convex hulls of $N$ random points in a plane chosen according to a given distribution. The points may be chosen independently or they may be correlated. After a non-exhaustive survey of the somewhat sporadic literature and diverse methods used in the random convex hull problem, we present a unifying approach, based on the notion of support function of a closed curve and the associated Cauchy's formulae, that allows us to compute exactly the mean perimeter and the mean area enclosed by the convex polygon both in case of independent as well as correlated points. Our method demonstrates a beautiful link between the random convex hull problem and the subject of extreme value statistics. As an example of correlated points, we study here in detail the case when the points represent the vertices of $n$ independent random walks. In the continuum time limit this reduces to $n$ independent planar Brownian trajectories for which we compute exactly, for all $n$, the mean perimeter and the mean area of their global convex hull. Our results have relevant applications in ecology in estimating the home range of a herd of animals. Some of these results were announced recently in a short communication [Phys. Rev. Lett. {\bf 103}, 140602 (2009)].

preprint2009arXiv

The Index Distribution of Gaussian Random Matrices

We compute analytically, for large N, the probability distribution of the number of positive eigenvalues (the index N_{+}) of a random NxN matrix belonging to Gaussian orthogonal (β=1), unitary (β=2) or symplectic (β=4) ensembles. The distribution of the fraction of positive eigenvalues c=N_{+}/N scales, for large N, as Prob(c,N)\simeq\exp[-βN^2 Φ(c)] where the rate function Φ(c), symmetric around c=1/2 and universal (independent of $β$), is calculated exactly. The distribution has non-Gaussian tails, but even near its peak at c=1/2 it is not strictly Gaussian due to an unusual logarithmic singularity in the rate function.

preprint2009arXiv

The longest excursion of stochastic processes in nonequilibrium systems

We consider the excursions, i.e. the intervals between consecutive zeros, of stochastic processes that arise in a variety of nonequilibrium systems and study the temporal growth of the longest one l_{\max}(t) up to time t. For smooth processes, we find a universal linear growth < l_{\max}(t) > \simeq Q_{\infty} t with a model dependent amplitude Q_\infty. In contrast, for non-smooth processes with a persistence exponent θ, we show that < l_{\max}(t) > has a linear growth if θ< θ_c while < l_{\max}(t) > \sim t^{1-ψ} if θ> θ_c. The amplitude Q_{\infty} and the exponent ψare novel quantities associated to nonequilibrium dynamics. These behaviors are obtained by exact analytical calculations for renewal and multiplicative processes and numerical simulations for other systems such as the coarsening dynamics in Ising model as well as the diffusion equation with random initial conditions.

preprint2007arXiv

Extreme statistics of complex random and quantum chaotic states

An exact analytical description of extreme intensity statistics in complex random states is derived. These states have the statistical properties of the Gaussian and Circular Unitary Ensemble eigenstates of random matrix theory. Although the components are correlated by the normalization constraint, it is still possible to derive compact formulae for all values of the dimensionality N. The maximum intensity result slowly approaches the Gumbel distribution even though the variables are bounded, whereas the minimum intensity result rapidly approaches the Weibull distribution. Since random matrix theory is conjectured to be applicable to chaotic quantum systems, we calculate the extreme eigenfunction statistics for the standard map with parameters at which its classical map is fully chaotic. The statistical behaviors are consistent with the finite-N formulae.

preprint2005arXiv

Strong clustering of non-interacting, passive sliders driven by a Kardar-Parisi-Zhang surface

We study the clustering of passive, non-interacting particles moving under the influence of a fluctuating field and random noise, in one dimension. The fluctuating field in our case is provided by a surface governed by the Kardar-Parisi-Zhang (KPZ) equation and the sliding particles follow the local surface slope. As the KPZ equation can be mapped to the noisy Burgers equation, the problem translates to that of passive scalars in a Burgers fluid. We study the case of particles moving in the same direction as the surface, equivalent to advection in fluid language. Monte-Carlo simulations on a discrete lattice model reveal extreme clustering of the passive particles. The resulting Strong Clustering State is defined using the scaling properties of the two point density-density correlation function. Our simulations show that the state is robust against changing the ratio of update speeds of the surface and particles. In the equilibrium limit of a stationary surface and finite noise, one obtains the Sinai model for random walkers on a random landscape. In this limit, we obtain analytic results which allow closed form expressions to be found for the quantities of interest. Surprisingly, these results for the equilibrium problem show good agreement with the results in the non-equilibrium regime.

preprint2004arXiv

Equilibrium statistics of a slave estimator in Langevin processes

We analyze the statistics of an estimator, denoted by xi_t and referred to as the slave, for the equilibrium susceptibility of a one dimensional Langevin process x_t in a potential phi(x). The susceptibility can be measured by evolving the slave equation in conjunction with the original Langevin process. This procedure yields a direct estimate of the susceptibility and avoids the need, when performing numerical simulations, to include applied external fields explicitly. The success of the method however depends on the statistical properties of the slave estimator. The joint probability density function for x_t and xi_t is analyzed. In the case where the potential of the system has a concave component the probability density function of the slave acquires a power law tail characterized by a temperature dependent exponent. Thus we show that while the average value of the slave, in the equilibrium state, is always finite and given by the fluctuation dissipation relation, higher moments and indeed the variance may show divergences. The behavior of the power law exponent is analyzed in a general context and it is calculated explicitly in some specific examples. Our results are confirmed by numerical simulations and we discuss possible measurement discrepancies in the fluctuation dissipation relation which could arise due to this behavior.

preprint2000arXiv

Slow Logarithmic Decay of Magnetization in the Zero Temperature Dynamics of an Ising Spin Chain: Analogy to Granular Compaction

We study the zero temperature coarsening dynamics in an Ising chain in presence of a dynamically induced field that favors locally the `-' phase compared to the `+' phase. At late times, while the `+' domains still coarsen as $t^{1/2}$, the `-' domains coarsen slightly faster as $t^{1/2}\log (t)$. As a result, at late times, the magnetization decays slowly as, $m(t)=-1 +{\rm const.}/{\log (t)}$. We establish this behavior both analytically within an independent interval approximation (IIA) and numerically. In the zero volume fraction limit of the `+' phase, we argue that the IIA becomes asymptotically exact. Our model can be alternately viewed as a simple Ising model for granular compaction. At late times in our model, the system decays into a fully compact state (where all spins are `-') in a slow logarithmic manner $\sim 1/{\log (t)}$, a fact that has been observed in recent experiments on granular systems.

preprint1999arXiv

Nonequilibrium Phase Transition in a Model of Diffusion, Aggregation and Fragmentation

We study the nonequilibrium phase transition in a model of aggregation of masses allowing for diffusion, aggregation on contact and fragmentation. The model undergoes a dynamical phase transition in all dimensions. The steady state mass distribution decays exponentially for large mass in one phase. On the contrary, in the other phase it has a power law tail and in addition an infinite aggregate. The model is solved exactly within a mean field approximation which keeps track of the distribution of masses. In one dimension, by mapping to an equivalent lattice gas model, exact steady states are obtained in two extreme limits of the parameter space. Critical exponents and the phase diagram are obtained numerically in one dimension. We also study the time dependent fluctuations in an equivalent interface model in (1+1) dimension and compute the roughness exponent $χ$ and the dynamical exponent z analytically in some limits and numerically otherwise. Two new fixed points of interface fluctuations in (1+1) dimension are identified. We also generalize our model to include arbitrary fragmentation kernels and solve the steady states exactly for some special choices of these kernels via mappings to other solvable models of statistical mechanics.

preprint1997arXiv

Far from equilibrium dynamics of the Bose gas

We study the approach to equilibrium of a Bose gas to a superfluid state. We point out that dynamic scaling, characteristic of far from equilibrium phase-ordering systems, should hold. We stress the importance of a non-dissipative Josephson precession term in driving the system to a new universality class. A model of coarsening in d=2, involving a quench between two temperatures below the equilibrium superfluid transition temperature (T_c), is exactly solved and demonstrates the relevance of the Josephson term. Numerical results on quenches from above T_c in d=2,3 provide evidence for the scaling picture postulated.

preprint1993arXiv

Annihilation of Immobile Reactants on the Bethe Lattice

Two-particle annihilation reaction, A+A -> inert, for immobile reactants on the Bethe lattice is solved exactly for the initially random distribution. The process reaches an absorbing state in which no nearest-neighbor reactants are left. The approach of the concentration to the limiting value is exponential. The solution reproduces the known one-dimensional result which is further extended to the reaction A+B -> inert.