• In this paper, we study the extreme statistics in the complex Ginibre ensemble of $N \times N$ random matrices with complex Gaussian entries, but with no other symmetries. All the $N$ eigenvalues are complex random variables and their joint distribution can be interpreted as a $2d$ Coulomb gas with a logarithmic repulsion between any pair of particles and in presence of a confining harmonic potential $v(r) \propto r^2$. We study the statistics of the eigenvalue with the largest modulus $r_{\max}$ in the complex plane. The typical and large fluctuations of $r_{\max}$ around its mean had been studied before, and they match smoothly to the right of the mean. However, it remained a puzzle to understand why the large and typical fluctuations to the left of the mean did not match. In this paper, we show that there is indeed an intermediate fluctuation regime that interpolates smoothly between the large and the typical fluctuations to the left of the mean. Moreover, we compute explicitly this "intermediate deviation function" (IDF) and show that it is universal, i.e. independent of the confining potential $v(r)$ as long as it is spherically symmetric and increases faster than $\ln r^2$ for large $r$ with an unbounded support. If the confining potential $v(r)$ has a finite support, i.e. becomes infinite beyond a finite radius, we show via explicit computation that the corresponding IDF is different. Interestingly, in the borderline case where the confining potential grows very slowly as $v(r) \sim \ln r^2$ for $r \gg 1$ with an unbounded support, the intermediate regime disappears and there is a smooth matching between the central part and the left large deviation regime.
  • We study the statistics of the kinetic (or equivalently potential) energy for $N$ non-interacting fermions in a $1d$ harmonic trap of frequency $\omega$, at finite temperature $T$. Remarkably, we find an exact solution for the full distribution of the kinetic energy, at any temperature $T$ and for any $N$, using a non-trivial mapping to an integrable Calogero-Moser-Sutherland model. As a function of temperature $T$, and for large $N$, we identify: (i) a quantum regime, for $T \sim \hbar \omega$, where quantum fluctuations dominate and (ii) a thermal regime, for $T \sim N \hbar \omega$, governed by thermal fluctuations. We show how the mean, the variance as well as the large deviation function associated with the distribution of the kinetic energy cross over from the quantum to the thermal regime as temperature increases.
  • We study the dynamics of a single active Brownian particle (ABP) in a two-dimensional harmonic trap. The active particle has an intrinsic time scale $D_R^{-1}$ set by the rotational diffusion with diffusion constant $D_R$. The harmonic trap also induces a relaxational time-scale $\mu^{-1}$. We show that the competition between these two time scales leads to a nontrivial time evolution for the ABP. At short times a strongly anisotropic motion emerges leading to anomalous persistence/first-passage properties. At long-times, the stationary position distribution in the trap exhibits two different behaviours: a Gaussian peak at the origin in the strongly passive limit ($D_R \to \infty$) and a delocalised ring away from the origin in the opposite strongly active limit ($D_R \to 0$). The predicted stationary behaviours in these limits are in agreement with recent experimental observations.
  • We consider a one-dimensional gas of $N$ charged particles confined by an external harmonic potential and interacting via the one-dimensional Coulomb potential. For this system we show that in equilibrium the charges settle, on an average, uniformly and symmetrically on a finite region centred around the origin. We study the statistics of the position of the rightmost particle $x_{\max}$ and show that the limiting distribution describing its typical fluctuations is different from the Tracy-Widom distribution found in the one-dimensional log-gas. We also compute the large deviation functions which characterise the atypical fluctuations of $x_{\max}$ far away from its mean value. In addition, we study the gap between the two rightmost particles as well as the index $N_+$, i.e., the number of particles on the positive semi-axis. We compute the limiting distributions associated to the typical fluctuations of these observables as well as the corresponding large deviation functions. We provide numerical supports to our analytical predictions. Part of these results were announced in a recent Letter, Phys. Rev. Lett. 119, 060601 (2017).
  • We compute the joint statistics of the momenta $p_i$ of $N$ non-interacting fermions in a trap, near the Fermi edge, with a particular focus on the largest one $p_{\max}$. For a $1d$ harmonic trap, momenta and positions play a symmetric role and hence, the joint statistics of momenta is identical to that of the positions. In particular, $p_{\max}$, as $x_{\max}$, is distributed according to the Tracy-Widom distribution. Here we show that novel "momentum edge statistics" emerge when the curvature of the potential vanishes, i.e. for "flat traps" near their minimum, with $V(x) \sim x^{2n}$ and $n>1$. These are based on generalisations of the Airy kernel that we obtain explicitly. The fluctuations of $p_{\max}$ are governed by new universal distributions determined from the $n$-th member of the second Painlev\'e hierarchy of non-linear differential equations, with connections to multicritical random matrix models. Finite temperature extensions and possible experimental signatures in cold atoms are discussed.
  • We study a system of 1D noninteracting spinless fermions in a confining trap at finite temperature. We first derive a useful and general relation for the fluctuations of the occupation numbers valid for arbitrary confining trap, as well as for both canonical and grand canonical ensembles. Using this relation, we obtain compact expressions, in the case of the harmonic trap, for the variance of certain observables of the form of sums of a function of the fermions' positions, $\mathcal{L}=\sum_n h(x_n)$. Such observables are also called linear statistics of the positions. As anticipated, we demonstrate explicitly that these fluctuations do depend on the ensemble in the thermodynamic limit, as opposed to averaged quantities, which are ensemble independent. We have applied our general formalism to compute the fluctuations of the number of fermions $\mathcal{N}_+$ on the positive axis at finite temperature. Our analytical results are compared to numerical simulations. We discuss the universality of the results with respect to the nature of the confinement.
  • The one-point distribution of the height for the continuum Kardar-Parisi-Zhang (KPZ) equation is determined numerically using the mapping to the directed polymer in a random potential at high temperature. Using an importance sampling approach, the distribution is obtained over a large range of values, down to a probability density as small as 10^{-1000} in the tails. Both short and long times are investigated and compared with recent analytical predictions for the large-deviation forms of the probability of rare fluctuations. At short times the agreement with the analytical expression is spectacular. We observe that the far left and right tails, with exponents 5/2 and 3/2 respectively, are preserved until large time. We present some evidence for the predicted non-trivial crossover in the left tail from the 5/2 tail exponent to the cubic tail of Tracy-Widom, although the details of the full scaling form remains beyond reach.
  • We establish an exact mapping between (i) the equilibrium (imaginary time) dynamics of non-interacting fermions trapped in a harmonic potential at temperature $T=1/\beta$ and (ii) non-intersecting Ornstein-Uhlenbeck (OU) particles constrained to return to their initial positions after time $\beta$. Exploiting the determinantal structure of the process we compute the universal correlation functions both in the bulk and at the edge of the trapped Fermi gas. The latter corresponds to the top path of the non-intersecting OU particles, and leads us to introduce and study the time-periodic Airy$_2$ process, ${\cal A}^b_2(u)$, depending on a single parameter, the period $b$. The standard Airy$_2$ process is recovered for $b=+\infty$. We discuss applications of our results to the real time quantum dynamics of trapped fermions.
  • We study analytically the Wigner function $W_N({\bf x},{\bf p})$ of $N$ noninteracting fermions trapped in a smooth confining potential $V({\bf x})$ in $d$ dimensions. At zero temperature, $W_N({\bf x},{\bf p})$ is constant over a finite support in the phase space $({\bf x},{\bf p})$ and vanishes outside. Near the edge of this support, we find a universal scaling behavior of $W_N({\bf x},{\bf p})$ for large $N$. The associated scaling function is independent of the precise shape of the potential as well as the spatial dimension $d$. We further generalize our results to finite temperature $T>0$. We show that there exists a low temperature regime $T \sim e_N/b$ where $e_N$ is an energy scale that depends on $N$ and the confining potential $V({\bf x})$, where the Wigner function at the edge again takes a universal scaling form with a $b$-dependent scaling function. This temperature dependent scaling function is also independent of the potential as well as the dimension $d$. Our results generalize to any $d\geq 1$ and $T \geq 0$ the $d=1$ and $T=0$ results obtained by Bettelheim and Wiegman [Phys. Rev. B ${\bf 84}$, 085102 (2011)].
  • We study $N$ noninteracting fermions in a domain bounded by a hard wall potential in $d \geq 1$ dimensions. We show that for large $N$, the correlations at the edge of the Fermi gas (near the wall) at zero temperature are described by a universal kernel, different from the universal edge kernel valid for smooth potentials. We compute this $d$ dimensional hard edge kernel exactly for a spherical domain and argue, using a generalized method of images, that it holds close to any sufficiently smooth boundary. As an application we compute the quantum statistics of the position of the fermion closest to the wall. Our results are then extended in several directions, including non-smooth boundaries such as a wedge, and also to finite temperature.
  • We consider a one-dimensional classical Coulomb gas of $N$ like-charges in a harmonic potential -- also known as the one-dimensional one-component plasma (1dOCP). We compute analytically the probability distribution of the position $x_{\max}$ of the rightmost charge in the limit of large $N$. We show that the typical fluctuations of $x_{\max}$ around its mean are described by a non-trivial scaling function, with asymmetric tails. This distribution is different from the Tracy-Widom distribution of $x_{\max}$ for the Dyson's log-gas. We also compute the large deviation functions of $x_{\max}$ explicitly and show that the system exhibits a third-order phase transition, as in the log-gas. Our theoretical predictions are verified numerically.
  • We consider the early time regime of the Kardar-Parisi-Zhang (KPZ) equation in $1+1$ dimensions in curved (or droplet) geometry. We show that for short time $t$, the probability distribution $P(H,t)$ of the height $H$ at a given point $x$ takes the scaling form $P(H,t) \sim \exp{\left(-\Phi_{\rm drop}(H)/\sqrt{t} \right)}$ where the rate function $\Phi_{\rm drop}(H)$ is computed exactly. While it is Gaussian in the center, i.e., for small $H$, the PDF has highly asymmetric non-Gaussian tails which we characterize in detail. This function $\Phi_{\rm drop}(H)$ is surprisingly reminiscent of the large deviation function describing the stationary fluctuations of finite size models belonging to the KPZ universality class. Thanks to a recently discovered connection between KPZ and free fermions, our results have interesting implications for the fluctuations of the rightmost fermion in a harmonic trap at high temperature and the full couting statistics at the edge.
  • We consider a one-dimensional random walk (RW) with a continuous and symmetric jump distribution, $f(\eta)$, characterized by a L\'evy index $\mu \in (0,2]$, which includes standard random walks ($\mu=2$) and L\'evy flights ($0<\mu<2$). We study the survival probability, $q(x_0,n)$, representing the probability that the RW stays non-negative up to step $n$, starting initially at $x_0 \geq 0$. Our main focus is on the $x_0$-dependence of $q(x_0,n)$ for large $n$. We show that $q(x_0,n)$ displays two distinct regimes as $x_0$ varies: (i) for $x_0= O(1)$ ("quantum regime"), the discreteness of the jump process significantly alters the standard scaling behavior of $q(x_0,n)$ and (ii) for $x_0 = O(n^{1/\mu})$ ("classical regime") the discrete-time nature of the process is irrelevant and one recovers the standard scaling behavior (for $\mu =2$ this corresponds to the standard Brownian scaling limit). The purpose of this paper is to study how precisely the crossover in $q(x_0,n)$ occurs between the quantum and the classical regime as one increases $x_0$.
  • We review recent advances on the record statistics of strongly correlated time series, whose entries denote the positions of a random walk or a L\'evy flight on a line. After a brief survey of the theory of records for independent and identically distributed random variables, we focus on random walks. During the last few years, it was indeed realized that random walks are a very useful "laboratory" to test the effects of correlations on the record statistics. We start with the simple one-dimensional random walk with symmetric jumps (both continuous and discrete) and discuss in detail the statistics of the number of records, as well as of the ages of the records, i.e., the lapses of time between two successive record breaking events. Then we review the results that were obtained for a wide variety of random walk models, including random walks with a linear drift, continuous time random walks, constrained random walks (like the random walk bridge) and the case of multiple independent random walkers. Finally, we discuss further observables related to records, like the record increments, as well as some questions raised by physical applications of record statistics, like the effects of measurement error and noise.
  • We consider $N$ non-interacting fermions in an isotropic $d$-dimensional harmonic trap. We compute analytically the cumulative distribution of the maximal radial distance of the fermions from the trap center at zero temperature. While in $d=1$ the limiting distribution (in the large $N$ limit), properly centered and scaled, converges to the squared Tracy-Widom distribution of the Gaussian Unitary Ensemble in Random Matrix Theory, we show that for all $d>1$, the limiting distribution converges to the Gumbel law. These limiting forms turn out to be universal, i.e., independent of the details of the trapping potential for a large class of isotropic trapping potentials. We also study the position of the right-most fermion in a given direction in $d$ dimensions and, in the case of a harmonic trap, the maximum momentum, and show that they obey similar Gumbel statistics. Finally, we generalize these results to low but finite temperature.
  • 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.
  • 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.
  • We study one-dimensional discrete as well as continuous time random walks, either with a fixed number of steps (for discrete time) $n$ or on a fixed time interval $T$ (for continuous time). In both cases, we focus on symmetric probability distribution functions (PDF) of jumps with a finite support $[-g_{max}, g_{max}]$. For continuous time random walks (CTRWs), the waiting time $\tau$ between two consecutive jumps is a random variable whose probability distribution (PDF) has a power law tail $\Psi(\tau) \propto \tau^{-1-\gamma}$, with $0<\gamma<1$. We obtain exact results for the joint statistics of the gap between the first two maximal positions of the random walk and the time elapsed between them. We show that for large $n$ (or large time $T$ for CTRW), this joint PDF reaches a stationary joint distribution which exhibits an interesting concentration effect in the sense that a gap close to its maximum possible value, $g\approx g_{max}$, is much more likely to be achieved by two successive jumps rather than by a long walk between the first two maxima. Our numerical simulations confirm this concentration effect.
  • The cytoskeleton is an inhomogeneous network of semi-flexible filaments, which are involved in a wide variety of active biological processes. Although the cytoskeletal filaments can be very stiff and embedded in a dense and cross-linked network, it has been shown that, in cells, they typically exhibit significant bending on all length scales. In this work we propose a model of a semi-flexible filament deformed by different types of cross-linkers for which one can compute and investigate the bending spectrum. Our model allows to couple the evolution of the deformation of the semi-flexible polymer with the stochastic dynamics of linkers which exert transversal forces onto the filament. We observe a $q^{-2}$ dependence of the bending spectrum for some biologically relevant parameters and in a certain range of wavenumbers $q$. However, generically, the spatially localized forcing and the non-thermal dynamics both introduce deviations from the thermal-like $q^{-2}$ spectrum.
  • 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\'evy 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{\pi} = 1.53362\ldots$.
  • 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.
  • 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.
  • Let $\mathcal{P}_{\beta}^{(V)} (N_{\cal I})$ be the probability that a $N\times N$ $\beta$-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}_{\beta}^{(V)} (N_{\cal I})$ for large $N$. We show that this probability scales for large $N$ as $\mathcal{P}_{\beta}^{(V)} (N_{\cal I})\approx \exp\left(-\beta N^2 \psi^{(V)}(N_{\cal I} /N)\right)$, where $\beta$ is the Dyson index of the ensemble. The rate function $\psi^{(V)}(k_{\cal I})$, independent of $\beta$, is computed in terms of single integrals that can be easily evaluated numerically. The general formalism is then applied to the classical $\beta$-Gaussian (${\cal I}=[-L,L]$), $\beta$-Wishart (${\cal I}=[1,L]$) and $\beta$-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.
  • 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.
  • 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<\mu\le 2$, known as the L\'evy 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.