BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gregory Schehr

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

ResearcherUnclaimed source record
cond-mat.stat-mechmath.PRcond-mat.dis-nnmath-phmath.MPcond-mat.quant-gasq-fin.STcond-mat.softphysics.data-anhep-thnlin.SIq-fin.GNq-fin.PR

Catalog footprint

What is connected

81works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

81 published item(s)

preprint2026arXiv

The distribution of the maximum of independent resetting Brownian motions

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

The Wishart--Rosenzweig--Porter random matrix ensemble

In recent years the Rosenzweig--Porter (RP) ensemble, obtained by adding a diagonal matrix with independent and identically distributed elements to a Gaussian random matrix, has been widely used as a minimal model for the emergence of fractal eigenstates in complex many-body systems. A key open question concerns the robustness of its phase diagram when the assumption of independent and uncorrelated entries is relaxed -- an assumption that simplifies its analysis, but is generally violated in realistic quantum systems. In this work, we take a first step in this direction by considering a deformed Wishart (rather than Gaussian) random matrix, which we dub the ``Wishart--RP'' ensemble. Using perturbation theory, as well as the cavity and replica methods and the Dyson Brownian motion approach, we characterize its phase diagram and localization properties. Remarkably, we show that the level compressibility, which quantifies spectral correlations in the fractal phase, coincides with that of the Gaussian RP model, thereby extending the universality conjectured in [SciPost Phys. 14, 110 (2023)] beyond the fully uncorrelated setting. We confirm our results with numerical tests.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quench dynamics of noninteracting fermions with a delta impurity

We study the out-of-equilibrium dynamics of noninteracting fermions in one dimension and in continuum space, in the presence of a delta impurity potential at the origin whose strength $g$ is varied at time $t=0$. The system is prepared in its ground state with $g=g_0=+\infty$, with two different densities and Fermi wave-vectors $k_L$ and $k_R$ on the two half-spaces $x>0$ and $x<0$ respectively. It then evolves for $t>0$ as an isolated system, with a finite impurity strength $g$. We compute exactly the time dependent density and current. For a fixed position $x$ and in the large time limit $t \to \infty$, the system reaches a non-equilibrium stationary state (NESS). We obtain analytically the correlation kernel, density, particle current, and energy current in the NESS, and characterize their relaxation, which is algebraic in time. In particular, in the NESS, we show that, away from the impurity, the particle density displays oscillations which are the non-equilibrium analog of the Friedel oscillations. In the regime of "rays", $x/t=ξ$ fixed with $x, t \to \infty$, we compute the same quantities and observe the emergence of two light cones, associated to the Fermi velocities $k_L$ and $k_R$ in the initial state. Interestingly, we find non trivial quantum correlations between two opposite rays with velocities $ξ$ and $-ξ$ which we compute explicitly. We extend to a continuum setting and to a correlated initial state the analytical methods developed in a recent work of Ljubotina, Sotiriadis and Prosen, in the context of a discrete fermionic chain with an impurity. We also generalize our results to an initial state at finite temperature, recovering, via explicit calculations, some predictions of conformal field theory in the low energy limit.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hole probability for noninteracting fermions in a $d$-dimensional trap

The hole probability, i.e., the probability that a region is void of particles, is a benchmark of correlations in many body systems. We compute analytically this probability $P(R)$ for a spherical region of radius $R$ in the case of $N$ noninteracting fermions in their ground state in a $d$-dimensional trapping potential. Using a connection to the Laguerre-Wishart ensembles of random matrices, we show that, for large $N$ and in the bulk of the Fermi gas, $P(R)$ is described by a universal scaling function of $k_F R$, for which we obtain an exact formula ($k_F$ being the local Fermi wave-vector). It exhibits a super exponential tail $P(R)\propto e^{- κ_d (k_F R)^{d+1}}$ where $κ_d$ is a universal amplitude, in good agreement with existing numerical simulations. When $R$ is of the order of the radius of the Fermi gas, the hole probability is described by a large deviation form which is not universal and which we compute exactly for the harmonic potential. Similar results also hold in momentum space.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)].

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tracy-Widom distributions for the Gaussian orthogonal and symplectic ensembles revisited: a skew-orthogonal polynomials approach

We study the distribution of the largest eigenvalue in the "Pfaffian" classical ensembles of random matrix theory, namely in the Gaussian orthogonal (GOE) and Gaussian symplectic (GSE) ensembles, using semi-classical skew-orthogonal polynomials, in analogue to the approach of Nadal and Majumdar (NM) for the Gaussian unitary ensemble (GUE). Generalizing the techniques of Adler, Forrester, Nagao and van Moerbeke, and using "overlapping Pfaffian" identities due to Knuth, we explicitly construct these semi-classical skew-orthogonal polynomials in terms of the semi-classical orthogonal polynomials studied by NM in the case of the GUE. With these polynomials we obtain expressions for the cumulative distribution functions of the largest eigenvalue in the GOE and the GSE. Further, by performing asymptotic analysis of these skew-orthogonal polynomials in the limit of large matrix size, we obtain an alternative derivation of the Tracy-Widom distributions for GOE and GSE. This asymptotic analysis relies on a certain Pfaffian identity, the proof of which employs the characterization of Pfaffians in terms of perfect matchings and link diagrams.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large time zero temperature dynamics of the spherical p=2-spin glass model of finite size

We revisit the long time dynamics of the spherical fully connected $p = 2$-spin glass model when the number of spins $N$ is large but {\it finite}. At $T=0$ where the system is in a (trivial) spin-glass phase, and on long time scale $t \gtrsim {\cal O}{(N^{2/3})}$ we show that the behavior of physical observables, like the energy, correlation and response functions, is controlled by the density of near-extreme eigenvalues at the edge of the spectrum of the coupling matrix $J$, and are thus non self-averaging. We show that the late time decay of these observables, once averaged over the disorder, is controlled by new universal exponents which we compute exactly.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite N corrections to the limiting distribution of the smallest eigenvalue of Wishart complex matrices

We study the probability distribution function (PDF) of the smallest eigenvalue of Laguerre-Wishart matrices $W = X^\dagger X$ where $X$ is a random $M \times N$ ($M \geq N$) matrix, with complex Gaussian independent entries. We compute this PDF in terms of semi-classical orthogonal polynomials, which are deformations of Laguerre polynomials. By analyzing these polynomials, and their associated recurrence relations, in the limit of large $N$, large $M$ with $M/N \to 1$ -- i.e. for quasi-square large matrices $X$ -- we show that this PDF, in the hard edge limit, can be expressed in terms of the solution of a Painlevé III equation, as found by Tracy and Widom, using Fredholm operators techniques. Furthermore, our method allows us to compute explicitly the first $1/N$ corrections to this limiting distribution at the hard edge. Our computations confirm a recent conjecture by Edelman, Guionnet and Péché. We also study the soft edge limit, when $M-N \sim {\cal O}(N)$, for which we conjecture the form of the first correction to the limiting distribution of the smallest eigenvalue.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The density of eigenvalues seen from the soft edge of random matrices in the Gaussian beta-ensembles

We characterize the phenomenon of "crowding" near the largest eigenvalue $λ_{\max}$ of random $N \times N$ matrices belonging to the Gaussian $β$-ensemble of random matrix theory, including in particular the Gaussian orthogonal ($β=1$), unitary ($β=2$) and symplectic ($β= 4$) ensembles. We focus on two distinct quantities: (i) the density of states (DOS) near $λ_{\max}$, $ρ_{\rm DOS}(r,N)$, which is the average density of eigenvalues located at a distance $r$ from $λ_{\max}$ (or the density of eigenvalues seen from $λ_{\max}$) and (ii) the probability density function of the gap between the first two largest eigenvalues, $p_{\rm GAP}(r,N)$. Using heuristic arguments as well as well numerical simulations, we generalize our recent exact analytical study of the Hermitian case (corresponding to $β= 2$). We also discuss some applications of these two quantities to statistical physics models.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Branching Brownian Motion Conditioned on Particle Numbers

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Diffusion in periodic, correlated random forcing landscapes

We study the dynamics of a Brownian particle in a strongly correlated quenched random potential defined as a periodically-extended (with period $L$) finite trajectory of a fractional Brownian motion with arbitrary Hurst exponent $H \in (0,1)$. While the periodicity ensures that the ultimate long-time behavior is diffusive, the generalised Sinai potential considered here leads to a strong logarithmic confinement of particle trajectories at intermediate times. These two competing trends lead to dynamical frustration and result in a rich statistical behavior of the diffusion coefficient $D_L$: Although one has the typical value $D^{\rm typ}_L \sim \exp(-βL^H)$, we show via an exact analytical approach that the positive moments ($k>0$) scale like $\langle D^k_L \rangle \sim \exp{[-c' (k βL^{H})^{1/(1+H)}]}$, and the negative ones as $\langle D^{-k}_L \rangle \sim \exp(a' (k βL^{H})^2)$, $c'$ and $a'$ being numerical constants and $β$ the inverse temperature. These results demonstrate that $D_L$ is strongly non-self-averaging. We further show that the probability distribution of $D_L$ has a log-normal left tail and a highly singular, one-sided log-stable right tail reminiscent of a Lifshitz singularity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universal Order and Gap Statistics of Critical Branching Brownian Motion

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalous fluctuations of currents in Sinai-type random chains with strongly correlated disorder

We study properties of a random walk in a generalized Sinai model, in which a quenched random potential is a trajectory of a fractional Brownian motion with arbitrary Hurst parameter H, 0< H <1, so that the random force field displays strong spatial correlations. In this case, the disorder-average mean-square displacement grows in proportion to log^{2/H}(n), n being time. We prove that moments of arbitrary order k of the steady-state current J_L through a finite segment of length L of such a chain decay as L^{-(1-H)}, independently of k, which suggests that despite a logarithmic confinement the average current is much higher than its Fickian counterpart in homogeneous systems. Our results reveal a paradoxical behavior such that, for fixed n and L, the mean square displacement decreases when one varies H from 0 to 1, while the average current increases. This counter-intuitive behavior is explained via an analysis of representative realizations of disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Asymmetric Lévy flights in the presence of absorbing boundaries

We consider a one dimensional asymmetric random walk whose jumps are identical, independent and drawn from a distribution ϕ(η) displaying asymmetric power law tails (i.e. ϕ(η) \sim c/η^{α+1} for large positive jumps and ϕ(η) \sim c/(γ|η|^{α+1}) for large negative jumps, with 0 < α< 2). In absence of boundaries and after a large number of steps n, the probability density function (PDF) of the walker position, x_n, converges to an asymmetric Lévy stable law of stability index αand skewness parameter β=(γ-1)/(γ+1). In particular the right tail of this PDF decays as c n/x_n^{1+α}. Much less is known when the walker is confined, or partially confined, in a region of the space. In this paper we first study the case of a walker constrained to move on the positive semi-axis and absorbed once it changes sign. In this case, the persistence exponent θ_+, which characterizes the algebraic large time decay of the survival probability, can be computed exactly and we show that the tail of the PDF of the walker position decays as c \, n/[(1-θ_+) \, x_n^{1+α}]. This last result can be generalized in higher dimensions such as a planar Lévy walker confined in a wedge with absorbing walls. Our results are corroborated by precise numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fluctuating interfaces subject to stochastic resetting

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Near-extreme eigenvalues and the first gap of Hermitian random matrices

We study the phenomenon of "crowding" near the largest eigenvalue $λ_{\max}$ of random $N \times N$ matrices belonging to the Gaussian Unitary Ensemble (GUE) of random matrix theory. We focus on two distinct quantities: (i) the density of states (DOS) near $λ_{\max}$, $ρ_{\rm DOS}(r,N)$, which is the average density of eigenvalues located at a distance $r$ from $λ_{\max}$ and (ii) the probability density function of the gap between the first two largest eigenvalues, $p_{\rm GAP}(r,N)$. In the edge scaling limit where $r = {\cal O}(N^{-1/6})$, which is described by a double scaling limit of a system of unconventional orthogonal polynomials, we show that $ρ_{\rm DOS}(r,N)$ and $p_{\rm GAP}(r,N)$ are characterized by scaling functions which can be expressed in terms of the solution of a Lax pair associated to the Painlevé XXXIV equation. This provides an alternative and simpler expression for the gap distribution, which was recently studied by Witte, Bornemann and Forrester in Nonlinearity 26, 1799 (2013). Our expressions allow to obtain precise asymptotic behaviors of these scaling functions both for small and large arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Near-extreme statistics of Brownian motion

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Extremes of N vicious walkers for large N: application to the directed polymer and KPZ interfaces

We compute the joint probability density function (jpdf) P_N(M, τ_M) of the maximum M and its position τ_M for N non-intersecting Brownian excursions, on the unit time interval, in the large N limit. For N \to \infty, this jpdf is peaked around M = \sqrt{2N} and τ_M = 1/2, while the typical fluctuations behave for large N like M - \sqrt{2N} \propto s N^{-1/6} and τ_M - 1/2 \propto w N^{-1/3} where s and w are correlated random variables. One obtains an explicit expression of the limiting jpdf P(s,w) in terms of the Tracy-Widom distribution for the Gaussian Orthogonal Ensemble (GOE) of Random Matrix Theory and a psi-function for the Hastings-McLeod solution to the Painlevé II equation. Our result yields, up to a rescaling of the random variables s and w, an expression for the jpdf of the maximum and its position for the Airy_2 process minus a parabola. This latter describes the fluctuations in many different physical systems belonging to the Kardar-Parisi-Zhang (KPZ) universality class in 1+1 dimensions. In particular, the marginal probability density function (pdf) P(w) yields, up to a model dependent length scale, the distribution of the endpoint of the directed polymer in a random medium with one free end, at zero temperature. In the large w limit one shows the asymptotic behavior \log P(w) \sim - w^3/12.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite-temperature and finite-time scaling of the directed polymer free-energy with respect to its geometrical fluctuations

We study the fluctuations of the directed polymer in 1+1 dimensions in a Gaussian random environment with a finite correlation length ξ and at finite temperature. We address the correspondence between the geometrical transverse fluctuations of the directed polymer, described by its roughness, and the fluctuations of its free-energy, characterized by its two-point correlator. Analytical arguments are provided in favor of a generic scaling law between those quantities, at finite time, non-vanishing ξ and explicit temperature dependence. Numerical results are in good agreement both for simulations on the discrete directed polymer and on a continuous directed polymer (with short-range correlated disorder). Applications to recent experiments on liquid crystals are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Large deviations of the top eigenvalue of large Cauchy random matrices

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lévy flights on the half line

We study the probability distribution function (pdf) of the position of a Lévy flight of index 0<α<2 in presence of an absorbing wall at the origin. The solution of the associated fractional Fokker-Planck equation can be constructed using a perturbation scheme around the Brownian solution (corresponding to α= 2), as an expansion in ε= 2 - α. We obtain an explicit analytical solution, exact at the first order in ε, which allows us to conjecture the precise asymptotic behavior of this pdf, including the first subleading corrections, for any α. Careful numerical simulations, as well as an exact computation for α= 1, confirm our conjecture.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the joint distribution of the maximum and its position of the Airy2 process minus a parabola

The maximal point of the Airy2 process minus a parabola is believed to describe the scaling limit of the end-point of the directed polymer in a random medium, which was proved to be true for a few specific cases. Recently two different formulas for the joint distribution of the location and the height of this maximal point were obtained, one by Moreno Flores, Quastel and Remenik, and the other by Schehr. The first formula is given in terms of the Airy function and an associated operator, and the second formula is expressed in terms of the Lax pair equations of the Painleve II equation. We give a direct proof that these two formulas are the same.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Record statistics and persistence for a random walk with a drift

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Record Statistics for Multiple Random Walks

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Super-Rough Glassy Phase of the Random Field XY Model in Two Dimensions

We study both analytically, using the renormalization group (RG) to two loop order, and numerically, using an exact polynomial algorithm, the disorder-induced glass phase of the two-dimensional XY model with quenched random symmetry-breaking fields and without vortices. In the super-rough glassy phase, i.e. below the critical temperature $T_c$, the disorder and thermally averaged correlation function $B(r)$ of the phase field $θ(x)$, $B(r) = \bar{<[θ(x) - θ(x+ r) ]^2>}$ behaves, for $r \gg a$, as $B(r) \simeq A(τ) \ln^2 (r/a)$ where $r = |r|$ and $a$ is a microscopic length scale. We derive the RG equations up to cubic order in $τ= (T_c-T)/T_c$ and predict the universal amplitude ${A}(τ) = 2τ^2-2τ^3 + {\cal O}(τ^4)$. The universality of $A(τ)$ results from nontrivial cancellations between nonuniversal constants of RG equations. Using an exact polynomial algorithm on an equivalent dimer version of the model we compute ${A}(τ)$ numerically and obtain a remarkable agreement with our analytical prediction, up to $τ\approx 0.5$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Distribution of the time at which N vicious walkers reach their maximal height

We study the extreme statistics of N non-intersecting Brownian motions (vicious walkers) over a unit time interval in one dimension. Using path-integral techniques we compute exactly the joint distribution of the maximum M and of the time τ_M at which this maximum is reached. We focus in particular on non-intersecting Brownian bridges ("watermelons without wall") and non-intersecting Brownian excursions ("watermelons with a wall"). We discuss in detail the relationships between such vicious walkers models in watermelons configurations and stochastic growth models in curved geometry on the one hand and the directed polymer in a disordered medium (DPRM) with one free end-point on the other hand. We also check our results using numerical simulations of Dyson's Brownian motion and confront them with numerical simulations of the Polynuclear Growth Model (PNG) and of a model of DPRM on a discrete lattice. Some of the results presented here were announced in a recent letter [J. Rambeau and G. Schehr, Europhys. Lett. 91, 60006 (2010)].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Maximum relative height of elastic interfaces in random media

The distribution of the maximal relative height (MRH) of self-affine one-dimensional elastic interfaces in a random potential is studied. We analyze the ground state configuration at zero driving force, and the critical configuration exactly at the depinning threshold, both for the random-manifold and random-periodic universality classes. These configurations are sampled by exact numerical methods, and their MRH distributions are compared with those with the same roughness exponent and boundary conditions, but produced by independent Fourier modes with normally distributed amplitudes. Using Pickands' theorem we derive an exact analytical description for the right tail of the latter. After properly rescaling the MRH distributions we find that corrections from the Gaussian independent modes approximation are in general small, as previously found for the average width distribution of depinning configurations. In the large size limit all corrections are finite except for the ground-state in the random-periodic class whose MRH distribution becomes, for periodic boundary conditions, indistinguishable from the Airy distribution. We find that the MRH distributions are, in general, sensitive to changes of boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Two stock options at the races: Black-Scholes forecasts

Suppose one buys two very similar stocks and is curious about how much, after some time T, one of them will contribute to the overall asset, expecting, of course, that it should be around 1/2 of the sum. Here we examine this question within the classical Black and Scholes (BS) model, focusing on the evolution of the probability density function P(w) of a random variable w = a_T^{(1)}/(a_T^{(1)} + a_T^{(2)}) where a_T^{(1)} and a_T^{(2)} are the values of two (either European- or the Asian-style) options produced by two absolutely identical BS stochastic equations. We show that within the realm of the BS model the behavior of P(w) is surprisingly different from common-sense-based expectations. For the European-style options P(w) always undergoes a transition, (when T approaches a certain threshold value), from a unimodal to a bimodal form with the most probable values being close to 0 and 1, and, strikingly, w =1/2 being the least probable value. This signifies that the symmetry between two options spontaneously breaks and just one of them completely dominates the sum. For path-dependent Asian-style options we observe the same anomalous behavior, but only for a certain range of parameters. Outside of this range, P(w) is always a bell-shaped function with a maximum at w = 1/2.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Universal Order Statistics of Random Walks

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extremal statistics of curved growing interfaces in 1+1 dimensions

We study the joint probability distribution function (pdf) of the maximum M of the height and its position X_M of a curved growing interface belonging to the universality class described by the Kardar-Parisi-Zhang equation in 1+1 dimensions. We obtain exact results for the closely related problem of p non-intersecting Brownian bridges where we compute the joint pdf P_p(M,τ_M) where τ_M is there the time at which the maximal height M is reached. Our analytical results, in the limit p \to \infty, become exact for the interface problem in the growth regime. We show that our results, for moderate values of p \sim 10 describe accurately our numerical data of a prototype of these systems, the polynuclear growth model in droplet geometry. We also discuss applications of our results to the ground state configuration of the directed polymer in a random potential with one fixed endpoint.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Extreme value statistics from the Real Space Renormalization Group: Brownian Motion, Bessel Processes and Continuous Time Random Walks

We use the Real Space Renormalization Group (RSRG) method to study extreme value statistics for a variety of Brownian motions, free or constrained such as the Brownian bridge, excursion, meander and reflected bridge, recovering some standard results, and extending others. We apply the same method to compute the distribution of extrema of Bessel processes. We briefly show how the continuous time random walk (CTRW) corresponds to a non standard fixed point of the RSRG transformation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The longest excursion of fractional Brownian motion : numerical evidence of non-Markovian effects

We study, using exact numerical simulations, the statistics of the longest excursion l_{\max}(t) up to time t for the fractional Brownian motion with Hurst exponent 0<H<1. We show that in the large t limit, < l_{\max}(t) > \propto Q_\infty t where Q_\infty \equiv Q_\infty(H) depends continuously on H, and in a non trivial way. These results are compared with exact analytical results obtained recently for a renewal process with an associated persistence exponent θ= 1-H. This comparison shows that Q_\infty(H) carries the clear signature of non-Markovian effects for H\neq 1/2. The pre-asymptotic behavior of < l_{\max}(t)> is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The longest excursion of stochastic processes in nonequilibrium systems

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

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Super-Aging in two-dimensional random ferromagnets

We study the aging properties, in particular the two-time autocorrelations, of the two-dimensional randomly diluted Ising ferromagnet below the critical temperature via Monte-Carlo simulations. We find that the autocorrelation function displays additive aging $C(t,t_w)=C_{st}(t)+C_{ag}(t,t_w)$, where the stationary part $C_{st}$ decays algebraically. The aging part shows anomalous scaling $C_{ag}(t,t_w)={\cal C}(h(t)/h(t_w))$, where $h(u)$ is a non-homogeneous function excluding a $t/t_w$ scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Non Markovian persistence in the diluted Ising model at criticality

We investigate global persistence properties for the non-equilibrium critical dynamics of the randomly diluted Ising model. The disorder averaged persistence probability $\bar{{P}_c}(t)$ of the global magnetization is found to decay algebraically with an exponent $θ_c$ that we compute analytically in a dimensional expansion in $d=4-ε$. Corrections to Markov process are found to occur already at one loop order and $θ_c$ is thus a novel exponent characterizing this disordered critical point. Our result is thoroughly compared with Monte Carlo simulations in $d=3$, which also include a measurement of the initial slip exponent. Taking carefully into account corrections to scaling, $θ_c$ is found to be a universal exponent, independent of the dilution factor $p$ along the critical line at $T_c(p)$, and in good agreement with our one loop calculation.

0 citations0 reviewsNo code linkedOpen access