Source author record

Satya N. Majumdar

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

140works

24topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

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.

preprint2022arXiv

Some considerations about reviewing and open-access in scientific publishing

Scientific research changed profoundly over the last 30 years, in all its aspects. Scientific publishing has changed as well, mainly because of the strong increased number of submitted papers and because of the appearance of Open Access journals and publishers. We propose some reflections on these issues.

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.

Satya N. Majumdar

Research overview

Is this your work?

140 publications

The distribution of the maximum of independent resetting Brownian motions

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

Metropolis Monte Carlo sampling: convergence, localization transition and optimality

An exactly solvable predator prey model with resetting

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

Counting statistics for non-interacting fermions in a rotating trap

Freezing transitions of Brownian particles in confining potentials

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

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

Number of distinct sites visited by a resetting random walker

Optimal Resetting Brownian Bridges

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

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

Some considerations about reviewing and open-access in scientific publishing

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

Generating stochastic trajectories with global dynamical constraints

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

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

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

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

Stability of large complex systems with heterogeneous relaxation dynamics

Statistics of first-passage Brownian functionals

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

Active Brownian Motion in Two Dimensions

Anderson-like localization transition of random walks with resetting

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

Distribution of the Time Between Maximum and Minimum of Random Walks

Ising model with stochastic resetting

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

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

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

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

Position distribution in a generalised run and tumble process

Quantum dynamics with stochastic reset

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

Stochastic Resetting and Applications

Time between the maximum and the minimum of a stochastic process

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

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

Universal survival probability for a ddd-dimensional run-and-tumble particle

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

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

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

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

Conditioned random walks and interaction-driven condensation

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

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

Exact extreme value statistics at mixed order transitions

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

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

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

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

Record statistics for random walk bridges

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

Convex Hulls of Random Walks: Large-Deviation Properties

Dynamical transition in the temporal relaxation of stochastic processes under resetting

Effective Langevin equations for constrained stochastic processes

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

Extreme value statistics of correlated random variables

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

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

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

Random walk with random resetting to the maximum

Spatial Extent of Branching Brownian Motion

Statistics of the longest interval in renewal processes

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

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

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

Branching Brownian Motion Conditioned on Particle Numbers

Condensation transition in joint large deviations of linear statistics

Constraint driven condensation in large fluctuations of linear statistics

Diffusion with resetting in arbitrary spatial dimension

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

Instanton Approach to Large NNN Harish-Chandra-Itzykson-Zuber Integrals

Long-range correlations in a locally driven exclusion process

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

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

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

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

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

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

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

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

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