
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$ noninteracting 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 nontrivial mapping to an integrable CalogeroMoserSutherland 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
twodimensional 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 timescale $\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/firstpassage properties. At longtimes, 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 onedimensional gas of $N$ charged particles confined by an
external harmonic potential and interacting via the onedimensional 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 TracyWidom distribution found in the
onedimensional loggas. 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
semiaxis. 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$ noninteracting
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 TracyWidom 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
nonlinear 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 onepoint distribution of the height for the continuum
KardarParisiZhang (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 largedeviation 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 nontrivial crossover in the left tail from the 5/2 tail
exponent to the cubic tail of TracyWidom, although the details of the full
scaling form remains beyond reach.

We establish an exact mapping between (i) the equilibrium (imaginary time)
dynamics of noninteracting fermions trapped in a harmonic potential at
temperature $T=1/\beta$ and (ii) nonintersecting OrnsteinUhlenbeck (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 nonintersecting OU
particles, and leads us to introduce and study the timeperiodic 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 nonsmooth
boundaries such as a wedge, and also to finite temperature.

We consider a onedimensional classical Coulomb gas of $N$ likecharges in a
harmonic potential  also known as the onedimensional onecomponent 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
nontrivial scaling function, with asymmetric tails. This distribution is
different from the TracyWidom distribution of $x_{\max}$ for the Dyson's
loggas. We also compute the large deviation functions of $x_{\max}$ explicitly
and show that the system exhibits a thirdorder phase transition, as in the
loggas. Our theoretical predictions are verified numerically.

We consider the early time regime of the KardarParisiZhang (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 nonGaussian 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 onedimensional 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 nonnegative 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 discretetime 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 onedimensional 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$ noninteracting 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 TracyWidom 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 rightmost 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 onedimensional 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$ noninteracting spinless 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 KardarParisiZhang equation at finite time.
Consequently our result provides a finite $T$ generalization of the TracyWidom
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
presentday cold atom experiments, most notably our detailed predictions for
the temperature dependence of the fluctuations near the edge.

We study onedimensional 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 semiflexible 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
crosslinked 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
semiflexible filament deformed by different types of crosslinkers for which
one can compute and investigate the bending spectrum. Our model allows to
couple the evolution of the deformation of the semiflexible 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 nonthermal
dynamics both introduce deviations from the thermallike $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
KardarParisiZhang (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 nonmonotonic 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 onedimensional 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) spinglass 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 nearextreme eigenvalues at the edge of the spectrum of the
coupling matrix $J$, and are thus non selfaveraging. 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.