Showing 1–50 of 275 results for author: Doussal, P L

Large deviations in statistics of the local time and occupation time for a run and tumble particle

Authors: Soheli Mukherjee, Pierre Le Doussal, Naftali R. Smith

Abstract: We investigate the statistics of the local time $\mathcal{T} = \int_0^T δ(x(t)) dt$ that a run and tumble particle (RTP) $x(t)$ in one dimension spends at the origin, with or without an external drift. By relating the local time to the number of times the RTP crosses the origin, we find that the local time distribution $P(\mathcal{T})$ satisfies the large deviation principle… ▽ More We investigate the statistics of the local time $\mathcal{T} = \int_0^T δ(x(t)) dt$ that a run and tumble particle (RTP) $x(t)$ in one dimension spends at the origin, with or without an external drift. By relating the local time to the number of times the RTP crosses the origin, we find that the local time distribution $P(\mathcal{T})$ satisfies the large deviation principle $P(\mathcal{T}) \sim \, e^{-T \, I(\mathcal{T} / T)} $ in the large observation time limit $T \to \infty$. Remarkably, we find that in presence of drift the rate function $I(ρ)$ is nonanalytic: We interpret its singularity as dynamical phase transitions of first order. We then extend these results by studying the statistics of the amount of time $\mathcal{R}$ that the RTP spends inside a finite interval (i.e., the occupation time), with qualitatively similar results. In particular, this yields the long-time decay rate of the probability $P(\mathcal{R} = T)$ that the particle does not exit the interval up to time $T$. We find that the conditional endpoint distribution exhibits an interesting change of behavior from unimodal to bimodal as a function of the size of the interval. To study the occupation time statistics, we extend the Donsker-Varadhan large-deviation formalism to the case of RTPs, for general dynamical observables and possibly in the presence of an external potential. △ Less

Submitted 11 May, 2024; originally announced May 2024.

Comments: 17 pages, 5 figures

A Dyson Brownian motion model for weak measurements in chaotic quantum systems

Authors: Federico Gerbino, Pierre Le Doussal, Guido Giachetti, Andrea De Luca

Abstract: We consider a toy model for the study of monitored dynamics in a many-body quantum systems. We study the stochastic Schrodinger equation resulting from the continuous monitoring with a rate $Γ$ of a random hermitian operator chosen at every time from the gaussian unitary ensemble (GUE). Due to invariance by unitary transformations, the dynamics of the eigenvalues $\{λ_α\}_{α=1}^n$ of the density m… ▽ More We consider a toy model for the study of monitored dynamics in a many-body quantum systems. We study the stochastic Schrodinger equation resulting from the continuous monitoring with a rate $Γ$ of a random hermitian operator chosen at every time from the gaussian unitary ensemble (GUE). Due to invariance by unitary transformations, the dynamics of the eigenvalues $\{λ_α\}_{α=1}^n$ of the density matrix can be decoupled from that of the eigenvectors. Thus, stochastic equations are derived that exactly describe the dynamics of $λ$'s. We consider two regimes: in the presence of an extra dephasing term, which can be generated by imperfect quantum measurements, the density matrix has a stationary distribution, and we show that in the limit of large sizes the distribution of $λ$'s is described by an inverse Marchenko Pastur distribution. In the case of perfect measurements instead, purification eventually occurs and we focus on finite-time dynamics. In this case, remarkably, we find an exact solution for the joint probability distribution of $λ$'s at each time $t$ and for each size $n$. Two relevant regimes emerge: at small times $tΓ= O(1)$, the spectrum is in a Coulomb gas regime, with a well-defined continuous spectral distribution in the limit of $n\to\infty$. In that case, all moments of the density matrix become self-averaging and it is possible to characterize the entanglement spectrum exactly. In the limit of large times $t Γ= O(n)$ one enters instead a regime in which the eigenvalues are exponentially separated $\log(λ_α/λ_β) = O(Γt/n)$, but fluctuations $\sim O(\sqrt{Γt/n})$ play an essential role. We are still able to characterize the asymptotic behaviors of entanglement entropy in this regime. △ Less

Submitted 29 June, 2024; v1 submitted 1 January, 2024; originally announced January 2024.

Comments: 29 pages, 3 figures

Journal ref: Quantum Reports, 6(2), 200-230 (2024)

Anomalous scaling of heterogeneous elastic lines: a new picture from sample to sample fluctuations

Authors: Maximilien Bernard, Pierre Le Doussal, Alberto Rosso, Christophe Texier

Abstract: We study a discrete model of an heterogeneous elastic line with internal disorder, submitted to thermal fluctuations. The monomers are connected through random springs with independent and identically distributed elastic constants drawn from $p(k)\sim k^{μ-1}$ for $k\to0$. When $μ>1$, the scaling of the standard Edwards-Wilkinson model is recovered. When $μ<1$, the elastic line exhibits an anomalo… ▽ More We study a discrete model of an heterogeneous elastic line with internal disorder, submitted to thermal fluctuations. The monomers are connected through random springs with independent and identically distributed elastic constants drawn from $p(k)\sim k^{μ-1}$ for $k\to0$. When $μ>1$, the scaling of the standard Edwards-Wilkinson model is recovered. When $μ<1$, the elastic line exhibits an anomalous scaling of the type observed in many growth models and experiments. Here we derive and use the exact expression for the exact probability distribution of the line shape at equilibrium, as well as the spectral properties of the matrix containing the random couplings, to fully characterize the sample to sample fluctuations. Our results lead to novel scaling predictions that partially disagree with previous works, but which are corroborated by numerical simulations. We also provide a novel interpretation of the anomalous scaling in terms of the abrupt jumps in the line's shape that dominate the average value of the observable. △ Less

Submitted 14 March, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

Exact first-order effect of interactions on the ground-state energy of harmonically-confined fermions

Authors: Pierre Le Doussal, Naftali R. Smith, Nathan Argaman

Abstract: We consider a system of $N$ spinless fermions, interacting with each other via a power-law interaction $ε/r^n$, and trapped in an external harmonic potential $V(r) = r^2/2$, in $d=1,2,3$ dimensions. For any $0 < n < d+2$, we obtain the ground-state energy $E_N$ of the system perturbatively in $ε$, $E_{N}=E_{N}^{\left(0\right)}+εE_{N}^{\left(1\right)}+O\left(ε^{2}\right)$. We calculate… ▽ More We consider a system of $N$ spinless fermions, interacting with each other via a power-law interaction $ε/r^n$, and trapped in an external harmonic potential $V(r) = r^2/2$, in $d=1,2,3$ dimensions. For any $0 < n < d+2$, we obtain the ground-state energy $E_N$ of the system perturbatively in $ε$, $E_{N}=E_{N}^{\left(0\right)}+εE_{N}^{\left(1\right)}+O\left(ε^{2}\right)$. We calculate $E_{N}^{\left(1\right)}$ exactly, assuming that $N$ is such that the "outer shell" is filled. For the case of $n=1$ (corresponding to a Coulomb interaction for $d=3$), we extract the $N \gg 1$ behavior of $E_{N}^{\left(1\right)}$, focusing on the corrections to the exchange term with respect to the leading-order term that is predicted from the local density approximation applied to the Thomas-Fermi approximate density distribution. The leading correction contains a logarithmic divergence, and is of particular importance in the context of density functional theory. We also study the effect of the interactions on the fermions' spatial density. Finally, we find that our result for $E_{N}^{\left(1\right)}$ significantly simplifies in the case where $n$ is even. △ Less

Submitted 18 April, 2024; v1 submitted 15 November, 2023; originally announced November 2023.

Comments: 40 pages, 5 figures

"Tattered" membrane

Authors: Pierre Le Doussal, Leo Radzihovsky

Abstract: Ideal crystalline membranes, realized by graphene and other atomic monolayers, exhibit rich physics - a universal anomalous elasticity of the critical "flat" phase characterized by a negative Poisson ratio, universally singular elastic moduli, order-from-disorder and a crumpling transition. We formulate a generalized $D$-dimensional field theory, parameterized by an $O(d)\times O(D)$ tensor field… ▽ More Ideal crystalline membranes, realized by graphene and other atomic monolayers, exhibit rich physics - a universal anomalous elasticity of the critical "flat" phase characterized by a negative Poisson ratio, universally singular elastic moduli, order-from-disorder and a crumpling transition. We formulate a generalized $D$-dimensional field theory, parameterized by an $O(d)\times O(D)$ tensor field with an {\it energetic} longitudinal constraint. For a soft constraint the resulting field theory describes a new class of a fluctuating "tattered" membranes, exhibiting a nonzero density of topological connectivity defects - slits, cracks and faults at an effective medium level. For hard, infinite-coupling constraint, the model reproduces the conventional crystalline membrane and its crumpling transition, and thereby demonstrates the essence of the difference between an elastic membrane and conventional field theories. Two additional fixed points emerge within the critical manifold, (i) globally attractive, "isotropic" $O(d)\times O(D)$, and (ii) "transverse", which in $D=2$ is the exact "dual" of the elastic membrane. Their properties are obtained in general $D,d$ from the renormalization group and the self-consistent screening analyses. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 6 pages + 25 pages of Supplemental Material, 2 figures

Linear statistics for Coulomb gases: higher order cumulants

Authors: Benjamin De Bruyne, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We consider $N$ classical particles interacting via the Coulomb potential in spatial dimension $d$ and in the presence of an external trap, at equilibrium at inverse temperature $β$. In the large $N$ limit, the particles are confined within a droplet of finite size. We study smooth linear statistics, i.e. the fluctuations of sums of the form ${\cal L}_N = \sum_{i=1}^N f({\bf x}_i)$, where… ▽ More We consider $N$ classical particles interacting via the Coulomb potential in spatial dimension $d$ and in the presence of an external trap, at equilibrium at inverse temperature $β$. In the large $N$ limit, the particles are confined within a droplet of finite size. We study smooth linear statistics, i.e. the fluctuations of sums of the form ${\cal L}_N = \sum_{i=1}^N f({\bf x}_i)$, where ${\bf x}_i$'s are the positions of the particles and where $f({\bf x}_i)$ is a sufficiently regular function. There exists at present standard results for the first and second moments of ${\cal L}_N$ in the large $N$ limit, as well as associated Central Limit Theorems in general dimension and for a wide class of confining potentials. Here we obtain explicit expressions for the higher order cumulants of ${\cal L}_N$ at large $N$, when the function $f({\bf x})=f(|{\bf x}|)$ and the confining potential are both rotationnally invariant. A remarkable feature of our results is that these higher cumulants depend only on the value of $f'(|{\bf x}|)$ and its higher order derivatives evaluated exactly at the boundary of the droplet, which in this case is a $d$-dimensional sphere. In the particular two-dimensional case $d=2$ at the special value $β=2$, a connection to the Ginibre ensemble allows us to derive these results in an alternative way using the tools of determinantal point processes. Finally we also obtain the large deviation form of the full probability distribution function of ${\cal L}_N$. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 19 pages

Journal ref: J. Phys. A: Math. Theor. 57, 155002 (2024)

Timeliness criticality in complex systems

Authors: José Moran, Matthijs Romeijnders, Pierre Le Doussal, Frank P. Pijpers, Utz Weitzel, Debabrata Panja, Jean-Philippe Bouchaud

Abstract: In complex systems, external parameters often determine the phase in which the system operates, i.e., its macroscopic behavior. For nearly a century, statistical physics has extensively studied systems' transitions across phases, (universal) critical exponents, and related dynamical properties. Here we consider the functionality of systems, notably operations in socio-technical ones, production in… ▽ More In complex systems, external parameters often determine the phase in which the system operates, i.e., its macroscopic behavior. For nearly a century, statistical physics has extensively studied systems' transitions across phases, (universal) critical exponents, and related dynamical properties. Here we consider the functionality of systems, notably operations in socio-technical ones, production in economic ones and, more generally, any schedule-based system, where timing is of crucial importance. We introduce a stylized model of delay propagation on temporal networks, where the magnitude of delay-mitigating buffer acts as a control parameter. The model exhibits {\it timeliness criticality}, a novel form of critical behavior. We characterize fluctuations near criticality, commonly referred to as ``avalanches'', and identify the corresponding critical exponents. The model exhibits timeliness criticality also when run on real-world temporal systems such as production networks. Additionally, we explore potential connections with the Mode-Coupling Theory of glasses, the depinning transition and the directed polymer problem. △ Less

Submitted 18 June, 2024; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: 12 pages of main paper, 12 figures, 9 pages of supplementary information

Journal ref: Nature Physics 2024, advanced online publication

Dynamics at the edge for independent diffusing particles

Authors: Pierre Le Doussal

Abstract: We study the dynamics of the outliers for a large number of independent Brownian particles in one dimension. We derive the multi-time joint distribution of the position of the rightmost particle, by two different methods. We obtain the two time joint distribution of the maximum and second maximum positions, and study the counting statistics at the edge. Finally we derive the multi-time joint distr… ▽ More We study the dynamics of the outliers for a large number of independent Brownian particles in one dimension. We derive the multi-time joint distribution of the position of the rightmost particle, by two different methods. We obtain the two time joint distribution of the maximum and second maximum positions, and study the counting statistics at the edge. Finally we derive the multi-time joint distribution of the running maximum, as well as the joint distribution of the arrival time of the first particle at several space points. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: 44 pages, 8 figures

Non-equilibrium phase transitions in active rank diffusions

Authors: Léo Touzo, Pierre Le Doussal

Abstract: We consider N run and tumble particles in one dimension interacting via a linear 1D Coulomb potential, an active version of the rank diffusion problem. It was solved previously for N = 2 leading to a stationary bound state in the attractive case. Here the evolution of the density fields is obtained in the large N limit in terms of two coupled Burger's type equations. In the attractive case the exa… ▽ More We consider N run and tumble particles in one dimension interacting via a linear 1D Coulomb potential, an active version of the rank diffusion problem. It was solved previously for N = 2 leading to a stationary bound state in the attractive case. Here the evolution of the density fields is obtained in the large N limit in terms of two coupled Burger's type equations. In the attractive case the exact stationary solution describes a non-trivial N-particle bound state, which exhibits transitions between a phase where the density is smooth with infinite support, a phase where the density has finite support and exhibits "shocks", i.e. clusters of particles, at the edges, and a fully clustered phase. In presence of an additional linear potential, the phase diagram, obtained for either sign of the interaction, is even richer, with additional partially expanding phases, with or without shocks. Finally, a general self-consistent method is introduced to treat more general interactions. The predictions are tested through extensive numerical simulations. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: 8 pages, 6 Figures. Supp. Mat.: 18 pages, 2 Figures

Probing the large deviations for the Beta random walk in random medium

Authors: Alexander K. Hartmann, Alexandre Krajenbrink, Pierre Le Doussal

Abstract: We consider a discrete-time random walk on a one-dimensional lattice with space and time-dependent random jump probabilities, known as the Beta random walk. We are interested in the probability that, for a given realization of the jump probabilities (a sample), a walker starting at the origin at time $t=0$ is at position beyond $ξ\sqrt{T/2}$ at time $T$. This probability fluctuates from sample to… ▽ More We consider a discrete-time random walk on a one-dimensional lattice with space and time-dependent random jump probabilities, known as the Beta random walk. We are interested in the probability that, for a given realization of the jump probabilities (a sample), a walker starting at the origin at time $t=0$ is at position beyond $ξ\sqrt{T/2}$ at time $T$. This probability fluctuates from sample to sample and we study the large-deviation rate function which characterizes the tails of its distribution at large time $T \gg 1$. It is argued that, up to a simple rescaling, this rate function is identical to the one recently obtained exactly by two of the authors for the continuum version of the model. That continuum model also appears in the macroscopic fluctuation theory of a class of lattice gases, e.g. in the so-called KMP model of heat transfer. An extensive numerical simulation of the Beta random walk, based on an importance sampling algorithm, is found in good agreement with the detailed analytical predictions. A first-order transition in the tilted measure, predicted to occur in the continuum model, is also observed in the numerics. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 17 pages

Fluctuations in the active Dyson Brownian motion and the overdamped Calogero-Moser model

Authors: Leo Touzo, Pierre Le Doussal, Gregory Schehr

Abstract: Recently, we introduced the active Dyson Brownian motion model (DBM), in which $N$ run-and-tumble particles interact via a logarithmic repulsive potential in the presence of a harmonic well. We found that in a broad range of parameters the density of particles converges at large $N$ to the Wigner semi-circle law, as in the passive case. In this paper, we provide an analytical support for this nume… ▽ More Recently, we introduced the active Dyson Brownian motion model (DBM), in which $N$ run-and-tumble particles interact via a logarithmic repulsive potential in the presence of a harmonic well. We found that in a broad range of parameters the density of particles converges at large $N$ to the Wigner semi-circle law, as in the passive case. In this paper, we provide an analytical support for this numerical observation, by studying the fluctuations of the positions of the particles in the nonequilibrium stationary state of the active DBM, in the regime of weak noise and large persistence time. In this limit, we obtain an analytical expression for the covariance between the particle positions for any $N$ from the exact inversion of the Hessian matrix of the system. We show that, when the number of particles is large $N \gg 1$, the covariance matrix takes scaling forms that we compute explicitly both in the bulk and at the edge of the support of the semi-circle. In the bulk, the covariance scales as $N^{-1}$, while at the edge, it scales as $N^{-2/3}$. Remarkably, we find that these results can be transposed directly to an equilibrium model, the overdamped Calogero-Moser model in the low temperature limit, providing an analytical confirmation of the numerical results by Agarwal, Kulkarni and Dhar. For this model, our method also allows us to obtain the equilibrium two-time correlations and their dynamical scaling forms both in the bulk and at the edge. Our predictions at the edge are reminiscent of a recent result in the mathematics literature by Gorin and Kleptsyn on the (passive) DBM. That result can be recovered by the present methods, and also, as we show, using the stochastic Airy operator. Finally, our analytical predictions are confirmed by precise numerical simulations, in a wide range of parameters. △ Less

Submitted 21 December, 2023; v1 submitted 26 July, 2023; originally announced July 2023.

Comments: 36 pages, 13 figures. New results added (e.g. anomalous diffusion and linear statistics) not included in published version. Accepted for publication in Phys. Rev. E

The weak noise theory of the O'Connell-Yor polymer as an integrable discretisation of the nonlinear Schrodinger equation

Authors: Alexandre Krajenbrink, Pierre Le Doussal

Abstract: We investigate and solve the weak noise theory for the semi-discrete O'Connell-Yor directed polymer. In the large deviation regime, the most probable evolution of the partition function obeys a classical non-linear system which is a non-standard discretisation of the nonlinear Schrodinger equation with mixed initial-final conditions. We show that this system is integrable and find its general solu… ▽ More We investigate and solve the weak noise theory for the semi-discrete O'Connell-Yor directed polymer. In the large deviation regime, the most probable evolution of the partition function obeys a classical non-linear system which is a non-standard discretisation of the nonlinear Schrodinger equation with mixed initial-final conditions. We show that this system is integrable and find its general solution through an inverse scattering method and a novel Fredholm determinant framework that we develop. This allows to obtain the large deviation rate function of the free energy of the polymer model from its conserved quantities and to study its convergence to the large deviations of the Kardar-Parisi-Zhang equation. Our model also degenerates to the classical Toda chain, which further substantiates the applicability of our Fredholm framework. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: 62 pages, 10 figures

Replica-symmetry breaking transitions in the large deviations of the ground-state of a spherical spin-glass

Authors: Bertrand Lacroix-A-Chez-Toine, Yan V. Fyodorov, Pierre Le Doussal

Abstract: We derive, within the replica formalism, a generalisation of the Crisanti-Sommers formula to describe the large deviation function (LDF) ${\cal L}(e)$ for the speed-$N$ atypical fluctuations of the intensive ground-state energy $e$ of a generic spherical spin-glass in the presence of a random external magnetic field of variance $Γ$. We then analyse our exact formula for the LDF in much detail for… ▽ More We derive, within the replica formalism, a generalisation of the Crisanti-Sommers formula to describe the large deviation function (LDF) ${\cal L}(e)$ for the speed-$N$ atypical fluctuations of the intensive ground-state energy $e$ of a generic spherical spin-glass in the presence of a random external magnetic field of variance $Γ$. We then analyse our exact formula for the LDF in much detail for the Replica symmetric, single step Replica Symmetry Breaking (1-RSB) and Full Replica Symmetry Breaking (FRSB) situations. Our main qualitative conclusion is that the level of RSB governing the LDF may be different from that for the typical ground state. We find that while the deepest ground-states are always controlled by a LDF of replica symmetric form, beyond a finite threshold $e\geq e_{t}$ a replica-symmetry breaking starts to be operative. These findings resolve the puzzling discrepancy between our earlier replica calculations for the $p=2$ spherical spin-glass and the rigorous results by Dembo and Zeitouni which we are able to reproduce invoking an 1-RSB pattern. Finally at an even larger critical energy $e_{c}\geq e_{t}$, acting as a "wall", the LDF diverges logarithmically, which we interpret as a change in the large deviation speed from $N$ to a faster growth. In addition, we show that in the limit $Γ\to 0$ the LDF takes non-trivial scaling forms (i) ${\cal L}(e) \sim G((e-e_c)/Γ)$ in the vicinity of the wall (ii) ${\cal L}(e) \sim Γ^{ην} F((e-e_{\rm typ})/Γ^ν)$ in the vicinity of the typical energy, characterised by two new exponents $η\geq 1$ and $ν$ characterising universality classes. Via matching the latter allows us to formulate several conjectures concerning the regime of {\it typical fluctuations}, identified as $e-e_{\rm typ} \sim N^{-1/η}$ and $Γ\sim N^{-1/(ην)}$. △ Less

Submitted 6 July, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: 58 pages, 14 figures

Interacting, running and tumbling: the active Dyson Brownian motion

Authors: Leo Touzo, Pierre Le Doussal, Gregory Schehr

Abstract: We introduce and study a model in one dimension of $N$ run-and-tumble particles (RTP) which repel each other logarithmically in the presence of an external quadratic potential. This is an "active'' version of the well-known Dyson Brownian motion (DBM) where the particles are subjected to a telegraphic noise, with two possible states $\pm$ with velocity $\pm v_0$. We study analytically and numerica… ▽ More We introduce and study a model in one dimension of $N$ run-and-tumble particles (RTP) which repel each other logarithmically in the presence of an external quadratic potential. This is an "active'' version of the well-known Dyson Brownian motion (DBM) where the particles are subjected to a telegraphic noise, with two possible states $\pm$ with velocity $\pm v_0$. We study analytically and numerically two different versions of this model. In model I a particle only interacts with particles in the same state, while in model II all the particles interact with each other. In the large time limit, both models converge to a steady state where the stationary density has a finite support. For finite $N$, the stationary density exhibits singularities, which disappear when $N \to +\infty$. In that limit, for model I, using a Dean-Kawasaki approach, we show that the stationary density of $+$ (respectively $-$) particles deviates from the DBM Wigner semi-circular shape, and vanishes with an exponent $3/2$ at one of the edges. In model II, the Dean-Kawasaki approach fails but we obtain strong evidence that the density in the large $N$ limit retains a Wigner semi-circular shape. △ Less

Submitted 23 November, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: Main text: 8 pages, 6 Figures. Supp. Mat.: 28 pages, 15 Figures. Typos corrected

Journal ref: EPL 142, 61004 (2023)

Out of equilibrium dynamics of repulsive ranked diffusions: the expanding crystal

Authors: Ana Flack, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We study the non-equilibrium Langevin dynamics of $N$ particles in one dimension with Coulomb repulsive linear interactions. This is a dynamical version of the so-called jellium model (without confinement) also known as ranked diffusion. Using a mapping to the Lieb-Liniger model of quantum bosons, we obtain an exact formula for the joint distribution of the positions of the $N$ particles at time… ▽ More We study the non-equilibrium Langevin dynamics of $N$ particles in one dimension with Coulomb repulsive linear interactions. This is a dynamical version of the so-called jellium model (without confinement) also known as ranked diffusion. Using a mapping to the Lieb-Liniger model of quantum bosons, we obtain an exact formula for the joint distribution of the positions of the $N$ particles at time $t$, all starting from the origin. A saddle point analysis shows that the system converges at large time to a linearly expanding crystal. Properly rescaled, this dynamical state resembles the equilibrium crystal in a time dependent effective quadratic potential. This analogy allows to study the fluctuations around the perfect crystal, which, to leading order, are Gaussian. There are however deviations from this Gaussian behavior, which embody long-range correlations of purely dynamical origin, characterized by the higher order cumulants of, e.g., the gaps between the particles, that we calculate exactly. We complement these results using a recent approach by one of us in terms of a noisy Burgers equation. In the large $N$ limit, the mean density of the gas can be obtained at any time from the solution of a deterministic viscous Burgers equation. This approach provides a quantitative description of the dense regime at shorter times. Our predictions are in good agreement with numerical simulations for finite and large $N$. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 26 pages, 8 figures

Journal ref: Phys. Rev. E 107, 064105 (2023)

Stationary time correlations for fermions after a quench in the presence of an impurity

Authors: Gabriel Gouraud, Pierre Le Doussal, Gregory Schehr

Abstract: We consider the quench dynamics of non-interacting fermions in one dimension in the presence of a finite-size impurity at the origin. This impurity is characterized by general momentum-dependent reflection and transmission coefficients which are changed from ${\sf r}_0(k), {\sf t}_0(k)$ to ${\sf r}(k), {\sf t}(k)$ at time $t=0$. The initial state is at equilibrium with ${\sf t}_0(k)=0$ such that t… ▽ More We consider the quench dynamics of non-interacting fermions in one dimension in the presence of a finite-size impurity at the origin. This impurity is characterized by general momentum-dependent reflection and transmission coefficients which are changed from ${\sf r}_0(k), {\sf t}_0(k)$ to ${\sf r}(k), {\sf t}(k)$ at time $t=0$. The initial state is at equilibrium with ${\sf t}_0(k)=0$ such that the system is cut in two independent halves with ${\sf r}_0^R(k)$, ${\sf r}_0^L(k)$ respectively to the right and to the left of the impurity. We obtain the exact large time limit of the multi-time correlations. These correlations become time translationally invariant, and are non-zero in two different regimes: (i) for $x=O(1)$ where the system reaches a non-equilibrium steady state (NESS) (ii) for $x \sim t$, i.e., the ray-regime. For a repulsive impurity these correlations are independent of ${\sf r}_0^R(k)$, ${\sf r}_0^L(k)$, while in the presence of bound states they oscillate and memory effects persist. We show that these nontrivial relaxational properties can be retrieved in a simple manner from the large time behaviour of the single particle wave functions. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: Main text: 7 pages, 3 figures. Supp. Mat.: 31 pages, 4 figures

Journal ref: EPL, 142 (2023) 41001

Clusters in the critical branching Brownian motion

Authors: Benoît Ferté, Pierre Le Doussal, Alberto Rosso, Xiangyu Cao

Abstract: Brownian particles that are replicated and annihilated at equal rate have strongly correlated positions, forming a few compact clusters separated by large gaps. We characterize the distribution of the particles at a given time, using a definition of clusters in terms a coarse-graining length recently introduced by some of us. We show that, in a non-extinct realization, the average number of cluste… ▽ More Brownian particles that are replicated and annihilated at equal rate have strongly correlated positions, forming a few compact clusters separated by large gaps. We characterize the distribution of the particles at a given time, using a definition of clusters in terms a coarse-graining length recently introduced by some of us. We show that, in a non-extinct realization, the average number of clusters grows as $\sim t^{D_{\mathrm{f}}/2}$ where $D_{\mathrm{f}} \approx 0.22$ is the Haussdoff dimension of the boundary of the super-Brownian motion, found by Mueller, Mytnik, and Perkins. We also compute the distribution of gaps between consecutive particles. We find two regimes separated by the characteristic length scale $\ell = \sqrt{D/β}$ where $D$ is the diffusion constant and $β$ the branching rate. The average number of gaps greater than $g$ decays as $\sim g^{D_{\mathrm{f}}-2}$ for $g\ll \ell$ and $\sim g^{-D_{\mathrm{f}}}$ for $g \gg \ell$. Finally, conditioned on the number of particles $n$, the above distributions are valid for $g \ll \sqrt{n}$; the average number of gaps greater than $g \gg \sqrt{n}$ is much less than one, and decays as $\simeq 4 (g/\sqrt{n})^{-2}$, in agreement with the universal gap distribution predicted by Ramola, Majumdar, and Schehr. Our results interpolate between a dense super-Brownian motion regime and a large-gap regime, unifying two previously independent approaches. △ Less

Submitted 22 February, 2023; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 20 pages, 9 figure; v2, minor revisions, accepted version

Journal ref: J. Phys. A: Math. Theor. 56 115002 (2023)

Stationary measures of the KPZ equation on an interval from Enaud-Derrida's matrix product ansatz representation

Authors: Guillaume Barraquand, Pierre Le Doussal

Abstract: The stationary measures of the Kardar-Parisi-Zhang equation on an interval have been computed recently. We present a rather direct derivation of this result by taking the weak asymmetry limit of the matrix product ansatz for the asymmetric simple exclusion process. We rely on the matrix product ansatz representation of Enaud and Derrida, which allows to express the steady-state in terms of re-weig… ▽ More The stationary measures of the Kardar-Parisi-Zhang equation on an interval have been computed recently. We present a rather direct derivation of this result by taking the weak asymmetry limit of the matrix product ansatz for the asymmetric simple exclusion process. We rely on the matrix product ansatz representation of Enaud and Derrida, which allows to express the steady-state in terms of re-weighted simple random walks. In the continuum limit, its measure becomes a path integral (or re-weighted Brownian motion) of the form encountered in Liouville quantum mechanics, recovering the recent formula. △ Less

Submitted 18 November, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: v2: Added Section 2.4 about summation over the zero mode. 14 pages

Density profile of noninteracting fermions in a rotating $2d$ trap at finite temperature

Authors: Manas Kulkarni, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We study the average density of $N$ spinless noninteracting fermions in a $2d$ harmonic trap rotating with a constant frequency $Ω$ and in the presence of an additional repulsive central potential $γ/r^2$. The average density at zero temperature was recently studied in Phys. Rev. A $\textbf{103}$, 033321 (2021) and an interesting multi-layered "wedding cake" structure with a "hole" at the center w… ▽ More We study the average density of $N$ spinless noninteracting fermions in a $2d$ harmonic trap rotating with a constant frequency $Ω$ and in the presence of an additional repulsive central potential $γ/r^2$. The average density at zero temperature was recently studied in Phys. Rev. A $\textbf{103}$, 033321 (2021) and an interesting multi-layered "wedding cake" structure with a "hole" at the center was found for the density in the large $N$ limit. In this paper, we study the average density at finite temperature. We demonstrate how this "wedding-cake" structure is modified at finite temperature. These large $N$ results warrant going much beyond the standard Local Density Approximation. We also generalize our results to a wide variety of trapping potentials and demonstrate the universality of the associated scaling functions both in the bulk and at the edges of the "wedding-cake". △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 18 pages, 7 figures, 1 table

Journal ref: Phys. Rev. A 107, 023302 (2023)

Exact position distribution of a harmonically-confined run-and-tumble particle in two dimensions

Authors: Naftali R. Smith, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We consider an overdamped run-and-tumble particle in two dimensions, with self propulsion in an orientation that stochastically rotates by 90 degrees at a constant rate, clockwise or counter-clockwise with equal probabilities. In addition, the particle is confined by an external harmonic potential of stiffness $μ$, and possibly diffuses. We find the exact time-dependent distribution… ▽ More We consider an overdamped run-and-tumble particle in two dimensions, with self propulsion in an orientation that stochastically rotates by 90 degrees at a constant rate, clockwise or counter-clockwise with equal probabilities. In addition, the particle is confined by an external harmonic potential of stiffness $μ$, and possibly diffuses. We find the exact time-dependent distribution $P\left(x,y,t\right)$ of the particle's position, and in particular, the steady-state distribution $P_{\text{st}}\left(x,y\right)$ that is reached in the long-time limit. We also find $P\left(x,y,t\right)$ for a "free" particle, $μ=0$. We achieve this by showing that, under a proper change of coordinates, the problem decomposes into two statistically-independent one-dimensional problems, whose exact solution has recently been obtained. We then extend these results in several directions, to two such run-and-tumble particles with a harmonic interaction, to analogous systems of dimension three or higher, and by allowing stochastic resetting. △ Less

Submitted 16 November, 2022; v1 submitted 21 July, 2022; originally announced July 2022.

Comments: 15 pages, 5 figures

Journal ref: Phys. Rev. E 106, 054133 (2022)

A stationary model of non-intersecting directed polymers

Authors: Guillaume Barraquand, Pierre Le Doussal

Abstract: We consider the partition function $Z_{\ell}(\vec x,0\vert \vec y,t)$ of $\ell$ non-intersecting continuous directed polymers of length $t$ in dimension $1+1$, in a white noise environment, starting from positions $\vec x$ and terminating at positions $\vec y$. When $\ell=1$, it is well known that for fixed $x$, the field $\log Z_1(x,0\vert y,t)$ solves the Kardar-Parisi-Zhang equation and admits… ▽ More We consider the partition function $Z_{\ell}(\vec x,0\vert \vec y,t)$ of $\ell$ non-intersecting continuous directed polymers of length $t$ in dimension $1+1$, in a white noise environment, starting from positions $\vec x$ and terminating at positions $\vec y$. When $\ell=1$, it is well known that for fixed $x$, the field $\log Z_1(x,0\vert y,t)$ solves the Kardar-Parisi-Zhang equation and admits the Brownian motion as a stationary measure. In particular, as $t$ goes to infinity, $Z_1(x,0\vert y,t)/Z_1(x,0\vert 0,t) $ converges to the exponential of a Brownian motion $B(y)$. In this article, we show an analogue of this result for any $\ell$. We show that $Z_{\ell}(\vec x,0\vert \vec y,t)/Z_{\ell}(\vec x,0\vert \vec 0,t) $ converges as $t$ goes to infinity to an explicit functional $Z_{\ell}^{\rm stat}(\vec y)$ of $\ell$ independent Brownian motions. This functional $Z_{\ell}^{\rm stat}(\vec y)$ admits a simple description as the partition sum for $\ell$ non-intersecting semi-discrete polymers on $\ell$ lines. We discuss applications to the endpoints and midpoints distribution for long non-crossing polymers and derive explicit formulas in the case of two polymers. To obtain these results, we show that the stationary measure of the O'Connell-Warren multilayer stochastic heat equation is given by a collection of independent Brownian motions. This in turn is shown via analogous results in a discrete setup for the so-called log-gamma polymer and exploit the connection between non-intersecting log-gamma polymers and the geometric RSK correspondence found in arXiv:1110.3489. . △ Less

Submitted 20 February, 2023; v1 submitted 16 May, 2022; originally announced May 2022.

Comments: 30 pages, 8 figures

Journal ref: Journal of Physics A: Mathematical and Theoretical, Volume 56, Number 4, 2023

The crossover from the Macroscopic Fluctuation Theory to the Kardar-Parisi-Zhang equation controls the large deviations beyond Einstein's diffusion

Authors: Alexandre Krajenbrink, Pierre Le Doussal

Abstract: We study the crossover from the macroscopic fluctuation theory (MFT) which describes 1D stochastic diffusive systems at late times, to the weak noise theory (WNT) which describes the Kardar-Parisi-Zhang (KPZ) equation at early times. We focus on the example of the diffusion in a time-dependent random field, observed in an atypical direction which induces an asymmetry. The crossover is described by… ▽ More We study the crossover from the macroscopic fluctuation theory (MFT) which describes 1D stochastic diffusive systems at late times, to the weak noise theory (WNT) which describes the Kardar-Parisi-Zhang (KPZ) equation at early times. We focus on the example of the diffusion in a time-dependent random field, observed in an atypical direction which induces an asymmetry. The crossover is described by a non-linear system which interpolates between the derivative and the standard non-linear Schrodinger equations in imaginary time. We solve this system using the inverse scattering method for mixed-time boundary conditions introduced by us to solve the WNT. We obtain the rate function which describes the large deviations of the sample-to-sample fluctuations of the cumulative distribution of the tracer position. It exhibits a crossover as the asymmetry is varied, recovering both MFT and KPZ limits. We sketch how it is consistent with extracting the asymptotics of a Fredholm determinant formula, recently derived for sticky Brownian motions. The crossover mechanism studied here should generalize to a larger class of models described by the MFT. Our results apply to study extremal diffusion beyond Einstein's theory. △ Less

Submitted 10 April, 2022; originally announced April 2022.

Comments: 46 pages

Clusters in an epidemic model with long-range dispersal

Authors: Xiangyu Cao, Pierre Le Doussal, Alberto Rosso

Abstract: In presence of long range dispersal, epidemics spread in spatially disconnected regions known as clusters. Here, we characterize exactly their statistical properties in a solvable model, in both the supercritical (outbreak) and critical regimes. We identify two diverging length scales, corresponding to the bulk and the outskirt of the epidemic. We reveal a nontrivial critical exponent that governs… ▽ More In presence of long range dispersal, epidemics spread in spatially disconnected regions known as clusters. Here, we characterize exactly their statistical properties in a solvable model, in both the supercritical (outbreak) and critical regimes. We identify two diverging length scales, corresponding to the bulk and the outskirt of the epidemic. We reveal a nontrivial critical exponent that governs the cluster number, the distribution of their sizes and of the distances between them. We also discuss applications to depinning avalanches with long range elasticity. △ Less

Submitted 2 September, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 19 pages, 8 figures; v2: added an application to the Covid-19 outbreak in the US; v3: minor changes, published version

Journal ref: Phys. Rev. Lett. 129, 108301 (2022)

More on the Brownian force model: avalanche shapes, tip driven, higher $d$

Authors: Pierre Le Doussal

Abstract: The Brownian force model (BFM) is the mean-field model for the avalanches of an elastic interface slowly driven in a random medium. It describes the spatio-temporal statistics of the velocity field, and, to some extent is analytically tractable. We extend our previous studies to obtain several observables for the BFM with short range elasticity, related to the local jump sizes $S(x)$ and to the av… ▽ More The Brownian force model (BFM) is the mean-field model for the avalanches of an elastic interface slowly driven in a random medium. It describes the spatio-temporal statistics of the velocity field, and, to some extent is analytically tractable. We extend our previous studies to obtain several observables for the BFM with short range elasticity, related to the local jump sizes $S(x)$ and to the avalanche spatial extension in $d=1$, or the avalanche span in $d>1$. In $d=1$ we consider both driving (i) by an imposed force (ii) by an imposed displacement "at the tip" and obtain in each case the mean spatial shape $\langle S(x) \rangle$ at fixed extension, or at fixed seed to edge distance. We find that near an edge $x_e$, $S(x) \simeq σ|x-x_e|^3$ where $σ$ has a universal distribution that we obtain. We also obtain the spatiotemporal shape near the edge. In $d>1$ we obtain (i) the mean shape $\langle S(x_1,x_\perp) \rangle$ for a fixed span, which exhibits a non-trivial dependence in the transverse distance to the seed $x_\perp$ (ii) the mean shape around a point which has not moved, $\langle S(x) \rangle_{S(0)=0}$, which vanishes at the center as $|x|^{b_d}$ with non trivial exponents, $b_1=4$, $b_2=2 \sqrt{2}$ and $b_3=\frac{1}{2} (\sqrt{17}-1)$. We obtain the probability distributions in any $d$ for the maximal radius of an avalanche and the minimal distance of approach to a given point, as well as the probability of not hitting a cone in $d=2$. These results equivalently apply to the continuum limit of some branching diffusions, as detailed in a companion paper. △ Less

Submitted 22 March, 2022; v1 submitted 20 March, 2022; originally announced March 2022.

Comments: 75 pages, 8 figures, Companion paper arXiv:2203:10512

Equivalence of mean-field avalanches and branching diffusions: From the Brownian force model to the super-Brownian motion

Authors: Pierre Le Doussal

Abstract: We point out that the mean-field theory of avalanches in the dynamics of elastic interfaces, the so-called Brownian force model (BFM) developed recently in non-equilibrium statistical physics, is equivalent to the so-called super-Brownian motion (SBM) developed in probability theory, a continuum limit of branching processes related to {\it space-embedded} Galton-Watson trees. In particular the exa… ▽ More We point out that the mean-field theory of avalanches in the dynamics of elastic interfaces, the so-called Brownian force model (BFM) developed recently in non-equilibrium statistical physics, is equivalent to the so-called super-Brownian motion (SBM) developed in probability theory, a continuum limit of branching processes related to {\it space-embedded} Galton-Watson trees. In particular the exact solvability property recently (re-)discovered from the field theory in mean-field avalanches (the "instanton equation") maps onto the so-called Dawson-Watanabe 1968 duality property. In the light of this correspondence we compare the results obtained independently in the two fields, and transport some of them from one field to the other. In particular, we discuss a scaling limit of the branching Brownian motion which maps onto the continuum field theory of mean-field avalanches △ Less

Submitted 22 March, 2022; v1 submitted 20 March, 2022; originally announced March 2022.

Comments: Main text 16 pages, total with appendices 43 pages, companion paper arXiv:2203.10544

Half-space stationary Kardar-Parisi-Zhang equation beyond the Brownian case

Authors: Guillaume Barraquand, Alexandre Krajenbrink, Pierre Le Doussal

Abstract: We study the Kardar-Parisi-Zhang equation on the half-line $x \geqslant 0$ with Neumann type boundary condition. Stationary measures of the KPZ dynamics were characterized in recent work: they depend on two parameters, the boundary parameter $u$ of the dynamics, and the drift $-v$ of the initial condition at infinity. We consider the fluctuations of the height field when the initial condition is g… ▽ More We study the Kardar-Parisi-Zhang equation on the half-line $x \geqslant 0$ with Neumann type boundary condition. Stationary measures of the KPZ dynamics were characterized in recent work: they depend on two parameters, the boundary parameter $u$ of the dynamics, and the drift $-v$ of the initial condition at infinity. We consider the fluctuations of the height field when the initial condition is given by one of these stationary processes. At large time $t$, it is natural to rescale parameters as $(u,v)=t^{-1/3}(a,b)$ to study the critical region. In the special case $a+b=0$, treated in previous works, the stationary process is simply Brownian. However, these Brownian stationary measures are particularly relevant in the bound phase ($a<0$) but not in the unbound phase. For instance, starting from the flat or droplet initial data, the height field near the boundary converges to the stationary process with $a>0$ and $b=0$, which is not Brownian. For $a+b\geqslant 0$, we determine exactly the large time distribution $F_{a,b}^{\rm stat}$ of the height function $h(0,t)$. As an application, we obtain the exact covariance of the height field in a half-line at two times $1\ll t_1\ll t_2$ starting from stationary initial data, as well as estimates, when starting from droplet initial data, in the limit $t_1/t_2\to 1$. △ Less

Submitted 11 July, 2022; v1 submitted 21 February, 2022; originally announced February 2022.

Comments: 37 pages. v2: minor edits (This is the accepted version)

Journal ref: J. Phys. A: Math. Theor., Vol. 55, Num. 27 (2022)

Quench dynamics of noninteracting fermions with a delta impurity

Authors: Gabriel Gouraud, Pierre Le Doussal, Gregory Schehr

Abstract: 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. I… ▽ More 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. △ Less

Submitted 9 March, 2022; v1 submitted 10 February, 2022; originally announced February 2022.

Comments: 57 pages, 12 figures. References as well as a discussion on GGE added

Journal ref: J. Phys. A: Math. Theor. 55, 395001 (2022)

Counting statistics for non-interacting fermions in a rotating trap

Authors: Naftali R. Smith, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: 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… ▽ More 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. △ Less

Submitted 24 April, 2022; v1 submitted 26 December, 2021; originally announced December 2021.

Comments: 21 pages, 7 figures

Journal ref: Phys. Rev. A 105, 043315 (2022)

Universal out-of-equilibrium dynamics of 1D critical quantum systems perturbed by noise coupled to energy

Authors: Alexios Christopoulos, Pierre Le Doussal, Denis Bernard, Andrea De Luca

Abstract: We consider critical one dimensional quantum systems initially prepared in their groundstate and perturbed by a smooth noise coupled to the energy density. By using conformal field theory, we deduce a universal description of the out-of-equilibrium dynamics. In particular, the full time-dependent distribution of any $2$--pt chiral correlation function can be obtained from solving two coupled ordin… ▽ More We consider critical one dimensional quantum systems initially prepared in their groundstate and perturbed by a smooth noise coupled to the energy density. By using conformal field theory, we deduce a universal description of the out-of-equilibrium dynamics. In particular, the full time-dependent distribution of any $2$--pt chiral correlation function can be obtained from solving two coupled ordinary stochastic differential equations. In contrast with the general expectation of heating, we demonstrate that the system reaches a non-trivial and universal stationary state characterized by broad distributions. As an example, we analyse the local energy density: while its first moment diverges exponentially fast in time, the stationary distribution, which we derive analytically, is symmetric around a negative median and exhibits a fat tail with $3/2$ decay exponent. We obtain a similar result for the entanglement entropy production associated to a given interval of size $\ell$. The corresponding stationary distribution has a $3/2$ right tail for all $\ell$, and converges to a one-sided Levy stable for large $\ell$. Our results are benchmarked via analytical and numerical calculations for a chain of non-interacting spinless fermions with excellent agreement. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: 9 pages, 5 figures, Supplemental material (16 pages)

Journal ref: Phys. Rev. X 13, 011043, 2023

Ranked diffusion, delta Bose gas and Burgers equation

Authors: Pierre Le Doussal

Abstract: We study the diffusion of $N$ particles in one dimension interacting via a drift proportional to their rank. In the attractive case (self-gravitating gas) a mapping to the Lieb Liniger quantum model allows to obtain stationary time correlations, return probabilities and the decay rate to the stationary state. The rank field obeys a Burgers equation, which we analyze. It allows to obtain the statio… ▽ More We study the diffusion of $N$ particles in one dimension interacting via a drift proportional to their rank. In the attractive case (self-gravitating gas) a mapping to the Lieb Liniger quantum model allows to obtain stationary time correlations, return probabilities and the decay rate to the stationary state. The rank field obeys a Burgers equation, which we analyze. It allows to obtain the stationary density at large $N$ in an external potential $V(x)$ (in the repulsive case). In the attractive case the decay rate to the steady state is found to depend on the initial condition if its spatial decay is slow enough. Coulomb gas methods allow to study the final equilibrium at large $N$. △ Less

Submitted 21 August, 2021; originally announced August 2021.

Comments: Main text 6 pages, Supp Mat. 25 pages. 4 figures

Inverse scattering solution of the weak noise theory of the Kardar-Parisi-Zhang equation with flat and Brownian initial conditions

Authors: Alexandre Krajenbrink, Pierre Le Doussal

Abstract: We present the solution of the weak noise theory (WNT) for the Kardar-Parisi-Zhang equation in one dimension at short time for flat initial condition (IC). The non-linear hydrodynamic equations of the WNT are solved analytically through a connexion to the Zakharov-Shabat (ZS) system using its classical integrability. This approach is based on a recently developed Fredholm determinant framework pre… ▽ More We present the solution of the weak noise theory (WNT) for the Kardar-Parisi-Zhang equation in one dimension at short time for flat initial condition (IC). The non-linear hydrodynamic equations of the WNT are solved analytically through a connexion to the Zakharov-Shabat (ZS) system using its classical integrability. This approach is based on a recently developed Fredholm determinant framework previously applied to the droplet IC. The flat IC provides the case for a non-vanishing boundary condition of the ZS system and yields a richer solitonic structure comprising the appearance of multiple branches of the Lambert function. As a byproduct, we obtain the explicit solution of the WNT for the Brownian IC, which undergoes a dynamical phase transition. We elucidate its mechanism by showing that the related spontaneous breaking of the spatial symmetry arises from the interplay between two solitons with different rapidities. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 20 pages

Full counting statistics for interacting trapped fermions

Authors: Naftali R. Smith, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We study $N$ spinless fermions in their ground state confined by an external potential in one dimension with long range interactions of the general Calogero-Sutherland type. For some choices of the potential this system maps to standard random matrix ensembles for general values of the Dyson index $β$. In the fermion model $β$ controls the strength of the interaction, $β=2$ corresponding to the no… ▽ More We study $N$ spinless fermions in their ground state confined by an external potential in one dimension with long range interactions of the general Calogero-Sutherland type. For some choices of the potential this system maps to standard random matrix ensembles for general values of the Dyson index $β$. In the fermion model $β$ controls the strength of the interaction, $β=2$ corresponding to the noninteracting case. We study the quantum fluctuations of the number of fermions ${\cal N}_{\cal D}$ in a domain $\cal{D}$ of macroscopic size in the bulk of the Fermi gas. We predict that for general $β$ the variance of ${\cal N}_{\cal D}$ grows as $A_β \log N + B_β$ for $N \gg 1$ and we obtain a formula for $A_β$ and $B_β$. This is based on an explicit calculation for $β\in\left\{ 1,2,4\right\} $ and on a conjecture that we formulate for general $β$. This conjecture further allows us to obtain a universal formula for the higher cumulants of ${\cal N}_{\cal D}$. Our results for the variance in the microscopic regime are found to be consistent with the predictions of the Luttinger liquid theory with parameter $K = 2/β$, and allow to go beyond. In addition we present families of interacting fermion models in one dimension which, in their ground states, can be mapped onto random matrix models. We obtain the mean fermion density for these models for general interaction parameter $β$. In some cases the fermion density exhibits interesting transitions, for example we obtain a noninteracting fermion formulation of the Gross-Witten-Wadia model. △ Less

Submitted 26 November, 2021; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: 61 pages, 6 figures

Journal ref: SciPost Phys. 11, 110 (2021)

Steady state of the KPZ equation on an interval and Liouville quantum mechanics

Authors: Guillaume Barraquand, Pierre Le Doussal

Abstract: We obtain a simple formula for the stationary measure of the height field evolving according to the Kardar-Parisi-Zhang equation on the interval $[0,L]$ with general Neumann type boundary conditions and any interval size. This is achieved using the recent results of Corwin and Knizel (arXiv:2103.12253) together with Liouville quantum mechanics. Our formula allows to easily determine the stationary… ▽ More We obtain a simple formula for the stationary measure of the height field evolving according to the Kardar-Parisi-Zhang equation on the interval $[0,L]$ with general Neumann type boundary conditions and any interval size. This is achieved using the recent results of Corwin and Knizel (arXiv:2103.12253) together with Liouville quantum mechanics. Our formula allows to easily determine the stationary measure in various limits: KPZ fixed point on an interval, half-line KPZ equation, KPZ fixed point on a half-line, as well as the Edwards-Wilkinson equation on an interval. △ Less

Submitted 9 September, 2021; v1 submitted 31 May, 2021; originally announced May 2021.

Comments: Letter + supplementary material. v4: typos corrected, final version

Stationary nonequilibrium bound state of a pair of run and tumble particles

Authors: Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We study two interacting identical run and tumble particles (RTP's) in one dimension. Each particle is driven by a telegraphic noise, and in some cases, also subjected to a thermal white noise with a corresponding diffusion constant $D$. We are interested in the stationary bound state formed by the two RTP's in the presence of a mutual attractive interaction. The distribution of the relative coord… ▽ More We study two interacting identical run and tumble particles (RTP's) in one dimension. Each particle is driven by a telegraphic noise, and in some cases, also subjected to a thermal white noise with a corresponding diffusion constant $D$. We are interested in the stationary bound state formed by the two RTP's in the presence of a mutual attractive interaction. The distribution of the relative coordinate $y$ indeed reaches a steady state that we characterize in terms of the solution of a second-order differential equation. We obtain the explicit formula for the stationary probability $P(y)$ of $y$ for two examples of interaction potential $V(y)$. The first one corresponds to $V(y) \sim |y|$. In this case, for $D=0$ we find that $P(y)$ contains a delta function part at $y=0$, signaling a strong clustering effect, together with a smooth exponential component. For $D>0$, the delta function part broadens, leading instead to weak clustering. The second example is the harmonic attraction $V(y) \sim y^2$ in which case, for $D=0$, $P(y)$ is supported on a finite interval. We unveil an interesting relation between this two-RTP model with harmonic attraction and a three-state single RTP model in one dimension, as well as with a four-state single RTP model in two dimensions. We also provide a general discussion of the stationary bound state, including examples where it is not unique, e.g., when the particles cannot cross due to an additional short-range repulsion. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: 21 pages, 8 figures

Journal ref: Phys. Rev. E 104, 044103 (2021)

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

Authors: Gabriel Gouraud, Pierre Le Doussal, Gregory Schehr

Abstract: 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… ▽ More 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. △ Less

Submitted 17 April, 2021; originally announced April 2021.

Comments: Main text: 6 pages, 1 figure Supp mat: 15 pages, 6 figures

Journal ref: EPL 137, 50003 (2022)

The KPZ equation in a half space with flat initial condition and the unbinding of a directed polymer from an attractive wall

Authors: Guillaume Barraquand, Pierre Le Doussal

Abstract: We present an exact solution for the height distribution of the KPZ equation at any time $t$ in a half space with flat initial condition. This is equivalent to obtaining the free energy distribution of a polymer of length $t$ pinned at a wall at a single point. In the large $t$ limit a binding transition takes place upon increasing the attractiveness of the wall. Around the critical point we find… ▽ More We present an exact solution for the height distribution of the KPZ equation at any time $t$ in a half space with flat initial condition. This is equivalent to obtaining the free energy distribution of a polymer of length $t$ pinned at a wall at a single point. In the large $t$ limit a binding transition takes place upon increasing the attractiveness of the wall. Around the critical point we find the same statistics as in the Baik-Ben--Arous-Péché transition for outlier eigenvalues in random matrix theory. In the bound phase, we obtain the exact measure for the endpoint and the midpoint of the polymer at large time. We also unveil curious identities in distribution between partition functions in half-space and certain partition functions in full space for Brownian type initial condition. △ Less

Submitted 23 February, 2023; v1 submitted 16 April, 2021; originally announced April 2021.

Comments: 27 pages. v4: typos corrected. v3: references added, appendix reorganized (some sections numbers have changed). Accepted version

Journal ref: Phys. Rev. E 104, 024502 (2021)

Wigner function for noninteracting fermions in hard wall potentials

Authors: Benjamin De Bruyne, David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: The Wigner function $W_N({\bf x}, {\bf p})$ is a useful quantity to characterize the quantum fluctuations of an $N$-body system in its phase space. Here we study $W_N({\bf x}, {\bf p})$ for $N$ noninteracting spinless fermions in a $d$-dimensional spherical hard box of radius $R$ at temperature $T=0$. In the large $N$ limit, the local density approximation (LDA) predicts that… ▽ More The Wigner function $W_N({\bf x}, {\bf p})$ is a useful quantity to characterize the quantum fluctuations of an $N$-body system in its phase space. Here we study $W_N({\bf x}, {\bf p})$ for $N$ noninteracting spinless fermions in a $d$-dimensional spherical hard box of radius $R$ at temperature $T=0$. In the large $N$ limit, the local density approximation (LDA) predicts that $W_N({\bf x}, {\bf p}) \approx 1/(2 π\hbar)^d$ inside a finite region of the $({\bf x}, {\bf p})$ plane, namely for $|{\bf x}| < R$ and $|{\bf p}| < k_F$ where $k_F$ is the Fermi momentum, while $W_N({\bf x}, {\bf p})$ vanishes outside this region, or "droplet", on a scale determined by quantum fluctuations. In this paper we investigate systematically, in this quantum region, the structure of the Wigner function along the edge of this droplet, called the Fermi surf. In one dimension, we find that there are three distinct edge regions along the Fermi surf and we compute exactly the associated nontrivial scaling functions in each regime. We also study the momentum distribution $\hat ρ_N(p)$ and find a striking algebraic tail for very large momenta $\hat ρ_N(p) \propto 1/p^4$, well beyond $k_F$, reminiscent of a similar tail found in interacting quantum systems (discussed in the context of Tan's relation). We then generalize these results to higher $d$ and find, remarkably, that the scaling function close to the edge of the box is universal, i.e., independent of the dimension~$d$. △ Less

Submitted 30 April, 2021; v1 submitted 11 April, 2021; originally announced April 2021.

Comments: 31 pages, 14 figures

Journal ref: Phys. Rev. A 104, 013314 (2021)

The inverse scattering of the Zakharov-Shabat system solves the weak noise theory of the Kardar-Parisi-Zhang equation

Authors: Alexandre Krajenbrink, Pierre Le Doussal

Abstract: We solve the large deviations of the Kardar-Parisi-Zhang (KPZ) equation in one dimension at short time by introducing an approach which combines field theoretical, probabilistic and integrable techniques. We expand the program of the weak noise theory, which maps the large deviations onto a non-linear hydrodynamic problem, and unveil its complete solvability through a connection to the integrabili… ▽ More We solve the large deviations of the Kardar-Parisi-Zhang (KPZ) equation in one dimension at short time by introducing an approach which combines field theoretical, probabilistic and integrable techniques. We expand the program of the weak noise theory, which maps the large deviations onto a non-linear hydrodynamic problem, and unveil its complete solvability through a connection to the integrability of the Zakharov-Shabat system. Exact solutions, depending on the initial condition of the KPZ equation, are obtained using the inverse scattering method and a Fredholm determinant framework recently developed. These results, explicit in the case of the droplet geometry, open the path to obtain the complete large deviations for general initial conditions. △ Less

Submitted 25 September, 2022; v1 submitted 31 March, 2021; originally announced March 2021.

Comments: 35 pages - V3 Some details of the derivation of the general solution to the P,Q system have been added in Section S-E

Journal ref: Phys. Rev. Lett. 127, 064101 (2021)

Condensation transition in the late-time position of a Run-and-Tumble particle

Authors: Francesco Mori, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We study the position distribution $P(\vec{R},N)$ of a run-and-tumble particle (RTP) in arbitrary dimension $d$, after $N$ runs. We assume that the constant speed $v>0$ of the particle during each running phase is independently drawn from a probability distribution $W(v)$ and that the direction of the particle is chosen isotropically after each tumbling. The position distribution is clearly isotro… ▽ More We study the position distribution $P(\vec{R},N)$ of a run-and-tumble particle (RTP) in arbitrary dimension $d$, after $N$ runs. We assume that the constant speed $v>0$ of the particle during each running phase is independently drawn from a probability distribution $W(v)$ and that the direction of the particle is chosen isotropically after each tumbling. The position distribution is clearly isotropic, $P(\vec{R},N)\to P(R,N)$ where $R=|\vec{R}|$. We show that, under certain conditions on $d$ and $W(v)$ and for large $N$, a condensation transition occurs at some critical value of $R=R_c\sim O(N)$ located in the large deviation regime of $P(R,N)$. For $R<R_c$ (subcritical fluid phase), all runs are roughly of the same size in a typical trajectory. In contrast, an RTP trajectory with $R>R_c$ is typically dominated by a `condensate', i.e., a large single run that subsumes a finite fraction of the total displacement (supercritical condensed phase). Focusing on the family of speed distributions $W(v)=α(1-v/v_0)^{α-1}/v_0$, parametrized by $α>0$, we show that, for large $N$, $P(R,N)\sim \exp\left[-Nψ_{d,α}(R/N)\right]$ and we compute exactly the rate function $ψ_{d,α}(z)$ for any $d$ and $α$. We show that the transition manifests itself as a singularity of this rate function at $R=R_c$ and that its order depends continuously on $d$ and $α$. We also compute the distribution of the condensate size for $R>R_c$. Finally, we study the model when the total duration $T$ of the RTP, instead of the total number of runs, is fixed. Our analytical predictions are confirmed by numerical simulations, performed using a constrained Markov chain Monte Carlo technique, with precision $\sim 10^{-100}$. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 44 pages, 14 figures

Journal ref: Phys. Rev. E 103, 062134 (2021)

Thermal buckling transition of crystalline membranes in a field

Authors: Pierre Le Doussal, Leo Radzihovsky

Abstract: Two dimensional crystalline membranes in isotropic embedding space exhibit a flat phase with anomalous elasticity, relevant e.g., for graphene. Here we study their thermal fluctuations in the absence of exact rotational invariance in the embedding space. An example is provided by a membrane in an orientational field, tuned to a critical buckling point by application of in-plane stresses. Through a… ▽ More Two dimensional crystalline membranes in isotropic embedding space exhibit a flat phase with anomalous elasticity, relevant e.g., for graphene. Here we study their thermal fluctuations in the absence of exact rotational invariance in the embedding space. An example is provided by a membrane in an orientational field, tuned to a critical buckling point by application of in-plane stresses. Through a detailed analysis, we show that the transition is in a new universality class. The self-consistent screening method predicts a second order transition, with modified anomalous elasticity exponents at criticality, while the RG suggests a weakly first order transition. △ Less

Submitted 21 June, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: 5 pages (main text) + 25 pages (supplementary material), 3 figures; a published version for PRL, with minor edits and added references

Journal ref: Phys. Rev. Lett. 127, 015702 (2021)

Matrix Kesten Recursion, Inverse-Wishart Ensemble and Fermions in a Morse Potential

Authors: Tristan Gautié, Jean-Philippe Bouchaud, Pierre Le Doussal

Abstract: The random variable $1+z_1+z_1z_2+\dots$ appears in many contexts and was shown by Kesten to exhibit a heavy tail distribution. We consider natural extensions of this variable and its associated recursion to $N \times N$ matrices either real symmetric $β=1$ or complex Hermitian $β=2$. In the continuum limit of this recursion, we show that the matrix distribution converges to the inverse-Wishart en… ▽ More The random variable $1+z_1+z_1z_2+\dots$ appears in many contexts and was shown by Kesten to exhibit a heavy tail distribution. We consider natural extensions of this variable and its associated recursion to $N \times N$ matrices either real symmetric $β=1$ or complex Hermitian $β=2$. In the continuum limit of this recursion, we show that the matrix distribution converges to the inverse-Wishart ensemble of random matrices. The full dynamics is solved using a mapping to $N$ fermions in a Morse potential, which are non-interacting for $β=2$. At finite $N$ the distribution of eigenvalues exhibits heavy tails, generalizing Kesten's results in the scalar case. The density of fermions in this potential is studied for large $N$, and the power-law tail of the eigenvalue distribution is related to the properties of the so-called determinantal Bessel process which describes the hard edge universality of random matrices. For the discrete matrix recursion, using free probability in the large $N$ limit, we obtain a self-consistent equation for the stationary distribution. The relation of our results to recent works of Rider and Valkó, Grabsch and Texier, as well as Ossipov, is discussed. △ Less

Submitted 31 July, 2021; v1 submitted 20 January, 2021; originally announced January 2021.

Comments: 44 pages, 5 figures, 6 appendices

Journal ref: J. Phys. A: Math. Theor. 54 255201 (2021)

Impurities in systems of noninteracting trapped fermions

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We study the properties of spin-less non-interacting fermions trapped in a confining potential in one dimension but in the presence of one or more impurities which are modelled by delta function potentials. We use a method based on the single particle Green's function. For a single impurity placed in the bulk, we compute the density of the Fermi gas near the impurity. Our results, in addition to r… ▽ More We study the properties of spin-less non-interacting fermions trapped in a confining potential in one dimension but in the presence of one or more impurities which are modelled by delta function potentials. We use a method based on the single particle Green's function. For a single impurity placed in the bulk, we compute the density of the Fermi gas near the impurity. Our results, in addition to recovering the Friedel oscillations at large distance from the impurity, allow the exact computation of the density at short distances. We also show how the density of the Fermi gas is modified when the impurity is placed near the edge of the trap in the region where the unperturbed system is described by the Airy gas. Our method also allows us to compute the effective potential felt by the impurity both in the bulk and at the edge. In the bulk this effective potential is shown to be a universal function only of the local Fermi wave vector, or equivalently of the local fermion density. When the impurity is placed near the edge of the Fermi gas, the effective potential can be expressed in terms of Airy functions. For an attractive impurity placed far outside the support of the fermion density, we show that an interesting transition occurs where a single fermion is pulled out of the Fermi sea and forms a bound state with the impurity. This is a quantum analogue of the well-known Baik-Ben Arous-Péché (BBP) transition, known in the theory of spiked random matrices. The density at the location of the impurity plays the role of an order parameter. We also consider the case of two impurities in the bulk and compute exactly the effective force between them mediated by the background Fermi gas. △ Less

Submitted 24 February, 2021; v1 submitted 27 December, 2020; originally announced December 2020.

Comments: 35 pages, 8 figures. Revised version: references added and typos corrected

Journal ref: SciPost Phys. 10, 082 (2021)

Spherical spin glass model with external field

Authors: Jinho Baik, Elizabeth Collins-Woodfin, Pierre Le Doussal, Hao Wu

Abstract: We analyze the free energy and the overlaps in the 2-spin spherical Sherrington Kirkpatrick spin glass model with an external field for the purpose of understanding the transition between this model and the one without an external field. We compute the limiting values and fluctuations of the free energy as well as three types of overlaps in the setting where the strength of the external field goes… ▽ More We analyze the free energy and the overlaps in the 2-spin spherical Sherrington Kirkpatrick spin glass model with an external field for the purpose of understanding the transition between this model and the one without an external field. We compute the limiting values and fluctuations of the free energy as well as three types of overlaps in the setting where the strength of the external field goes to zero as the dimension of the spin variable grows. In particular, we consider overlaps with the external field, the ground state, and a replica. Our methods involve a contour integral representation of the partition function along with random matrix techniques. We also provide computations for the matching between different scaling regimes. Finally, we discuss the implications of our results for susceptibility and for the geometry of the Gibbs measure. Some of the findings of this paper are confirmed rigorously by Landon and Sosoe in their recent paper which came out independently and simultaneously. △ Less

Submitted 5 April, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 75 pages, 9 figures, 5 tables. Version 2: a few small typos are corrected. One author's last name (Collins-Woodfin) is changed from version 1. Version 3: corrected a few small typos and added one reference [11]

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

Authors: David S. Dean, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr, Naftali R. Smith

Abstract: 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… ▽ More 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. △ Less

Submitted 27 September, 2020; originally announced September 2020.

Comments: 35 pages, 9 figures

Journal ref: J. Phys. A: Math. Theor. 54 084001 (2021)

Tilted elastic lines with columnar and point disorder, non-Hermitian quantum mechanics and spiked random matrices: pinning and localization

Authors: Alexandre Krajenbrink, Pierre Le Doussal, Neil O'Connell

Abstract: We revisit the problem of an elastic line (e.g. a vortex line in a superconductor) subject to both columnar disorder and point disorder in dimension $d=1+1$. Upon applying a transverse field, a delocalization transition is expected, beyond which the line is tilted macroscopically. We investigate this transition in the fixed tilt angle ensemble and within a one-way model where backward jumps are ne… ▽ More We revisit the problem of an elastic line (e.g. a vortex line in a superconductor) subject to both columnar disorder and point disorder in dimension $d=1+1$. Upon applying a transverse field, a delocalization transition is expected, beyond which the line is tilted macroscopically. We investigate this transition in the fixed tilt angle ensemble and within a one-way model where backward jumps are neglected. From recent results about directed polymers and their connections to random matrix theory, we find that for a single line and a single strong defect this transition in presence of point disorder coincides with the Baik-Ben Arous-Peche (BBP) transition for the appearance of outliers in the spectrum of a perturbed random matrix in the GUE. This transition is conveniently described in the polymer picture by a variational calculation. In the delocalized phase, the ground state energy exhibits Tracy-Widom fluctuations. In the localized phase we show, using the variational calculation, that the fluctuations of the occupation length along the columnar defect are described by $f_{KPZ}$, a distribution which appears ubiquitously in the Kardar-Parisi-Zhang universality class. We then consider a smooth density of columnar defect energies. Depending on how this density vanishes at its lower edge we find either (i) a delocalized phase only (ii) a localized phase with a delocalization transition. We analyze this transition which is an infinite-rank extension of the BBP transition. The fluctuations of the ground state energy of a single elastic line in the localized phase (for fixed columnar defect energies) are described by a Fredholm determinant based on a new kernel. The case of many columns and many non-intersecting lines, relevant for the study of the Bose glass phase, is also analyzed. The ground state energy is obtained using free probability and the Burgers equation. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: 45 pages

Journal ref: Phys. Rev. E 103, 042120 (2021)

Spatial clustering of depinning avalanches in presence of long-range interactions

Authors: Clément Le Priol, Pierre Le Doussal, Alberto Rosso

Abstract: Disordered elastic interfaces display avalanche dynamics at the depinning transition. For short-range interactions, avalanches correspond to compact reorganizations of the interface well described by the depinning theory. For long-range elasticity, an avalanche is a collection of spatially disconnected clusters. In this paper we determine the scaling properties of the clusters and relate them to t… ▽ More Disordered elastic interfaces display avalanche dynamics at the depinning transition. For short-range interactions, avalanches correspond to compact reorganizations of the interface well described by the depinning theory. For long-range elasticity, an avalanche is a collection of spatially disconnected clusters. In this paper we determine the scaling properties of the clusters and relate them to the roughness exponent of the interface. The key observation of our analysis is the identification of a Bienaym{é}-Galton-Watson process describing the statistics of the cluster number. Our work has a concrete importance for experimental applications where the cluster statistics is a key probe of avalanche dynamics. △ Less

Submitted 23 August, 2020; originally announced August 2020.

Comments: 5pages + 6 pages SupMat, 4 + 7 figures

Journal ref: Phys. Rev. Lett. 126, 025702 (2021)

Counting statistics for non-interacting fermions in a $d$-dimensional potential

Authors: Naftali R. Smith, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We develop a first-principle approach to compute the counting statistics in the ground-state of $N$ noninteracting spinless fermions in a general potential in arbitrary dimensions $d$ (central for $d>1$). In a confining potential, the Fermi gas is supported over a bounded domain. In $d=1$, for specific potentials, this system is related to standard random matrix ensembles. We study the quantum flu… ▽ More We develop a first-principle approach to compute the counting statistics in the ground-state of $N$ noninteracting spinless fermions in a general potential in arbitrary dimensions $d$ (central for $d>1$). In a confining potential, the Fermi gas is supported over a bounded domain. In $d=1$, for specific potentials, this system is related to standard random matrix ensembles. We study the quantum fluctuations of the number of fermions ${\cal N}_{\cal D}$ in a domain $\cal{D}$ of macroscopic size in the bulk of the support. We show that the variance of ${\cal N}_{\cal D}$ grows as $N^{(d-1)/d} (A_d \log N + B_d)$ for large $N$, and obtain the explicit dependence of $A_d, B_d$ on the potential and on the size of ${\cal D}$ (for a spherical domain in $d>1$). This generalizes the free-fermion results for microscopic domains, given in $d=1$ by the Dyson-Mehta asymptotics from random matrix theory. This leads us to conjecture similar asymptotics for the entanglement entropy of the subsystem $\cal{D}$, in any dimension, supported by exact results for $d=1$. △ Less

Submitted 18 March, 2021; v1 submitted 3 August, 2020; originally announced August 2020.

Comments: Main text: 8 pages, 1 figure. Supplemental material: 22 pages, 6 figures

Journal ref: Phys. Rev. E 103, 030105 (2021)

Universal Properties of a Run-and-Tumble Particle in Arbitrary Dimension

Authors: Francesco Mori, Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: We consider an active run-and-tumble particle (RTP) in $d$ dimensions, starting from the origin and evolving over a time interval $[0,t]$. We examine three different models for the dynamics of the RTP: the standard RTP model with instantaneous tumblings, a variant with instantaneous runs and a general model in which both the tumblings and the runs are non-instantaneous. For each of these models, w… ▽ More We consider an active run-and-tumble particle (RTP) in $d$ dimensions, starting from the origin and evolving over a time interval $[0,t]$. We examine three different models for the dynamics of the RTP: the standard RTP model with instantaneous tumblings, a variant with instantaneous runs and a general model in which both the tumblings and the runs are non-instantaneous. For each of these models, we use the Sparre Andersen theorem for discrete-time random walks to compute exactly the probability that the $x$ component does not change sign up to time $t$, showing that it does not depend on $d$. As a consequence of this result, we compute exactly other $x$-component properties, namely the distribution of the time of the maximum and the record statistics, showing that they are universal, i.e. they do not depend on $d$. Moreover, we show that these universal results hold also if the speed $v$ of the particle after each tumbling is random, drawn from a generic probability distribution. Our findings are confirmed by numerical simulations. Some of these results have been announced in a recent Letter [Phys. Rev. Lett. 124, 090603 (2020)]. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: 34 pages, 18 figures. This is an extended version of arXiv:2001.01492, published in Physical Review Letters

Journal ref: Phys. Rev. E 102, 042133 (2020)

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

Authors: Pierre Le Doussal, Satya N. Majumdar, Gregory Schehr

Abstract: 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 under… ▽ More 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. △ Less

Submitted 18 March, 2020; originally announced March 2020.

Comments: 7 + 15 pages, 8 figures

Journal ref: EPL, 130 (2020) 40002

Half-space stationary Kardar-Parisi-Zhang equation

Authors: Guillaume Barraquand, Alexandre Krajenbrink, Pierre Le Doussal

Abstract: We study the solution of the Kardar-Parisi-Zhang (KPZ) equation for the stochastic growth of an interface of height $h(x,t)$ on the positive half line, equivalently the free energy of the continuum directed polymer in a half space with a wall at $x=0$. The boundary condition $\partial_x h(x,t)|_{x=0}=A$ corresponds to an attractive wall for $A<0$, and leads to the binding of the polymer to the wal… ▽ More We study the solution of the Kardar-Parisi-Zhang (KPZ) equation for the stochastic growth of an interface of height $h(x,t)$ on the positive half line, equivalently the free energy of the continuum directed polymer in a half space with a wall at $x=0$. The boundary condition $\partial_x h(x,t)|_{x=0}=A$ corresponds to an attractive wall for $A<0$, and leads to the binding of the polymer to the wall below the critical value $A=-1/2$. Here we choose the initial condition $h(x,0)$ to be a Brownian motion in $x>0$ with drift $-(B+1/2)$. When $A+B \to -1$, the solution is stationary, i.e. $h(\cdot,t)$ remains at all times a Brownian motion with the same drift, up to a global height shift $h(0,t)$. We show that the distribution of this height shift is invariant under the exchange of parameters $A$ and $B$. For any $A,B > - 1/2$, we provide an exact formula characterizing the distribution of $h(0,t)$ at any time $t$, using two methods: the replica Bethe ansatz and a discretization called the log-gamma polymer, for which moment formulae were obtained. We analyze its large time asymptotics for various ranges of parameters $A,B$. In particular, when $(A, B) \to (-1/2, -1/2)$, the critical stationary case, the fluctuations of the interface are governed by a universal distribution akin to the Baik-Rains distribution arising in stationary growth on the full-line. It can be expressed in terms of a simple Fredholm determinant, or equivalently in terms of the Painlevé II transcendent. This provides an analog for the KPZ equation, of some of the results recently obtained by Betea-Ferrari-Occelli in the context of stationary half-space last-passage-percolation. From universality, we expect that limiting distributions found in both models can be shown to coincide. △ Less

Submitted 13 January, 2021; v1 submitted 8 March, 2020; originally announced March 2020.

Comments: 53 pages

Journal ref: Journal of Statistical Physics (2020) 181:1149-1203