Search

We consider a toy model for the study of monitored dynamics in many-body quantum systems. We study the stochastic Schrödinger equation resulting from continuous monitoring with a rate Γ of a random Hermitian operator, drawn from the Gaussian unitary ensemble (GUE) at every time t. Due to invariance by unitary transformations, the dynamics of the eigenvalues {λα}α=1n of the density matrix decouples from that of the eigenvectors, and is exactly described by stochastic equations that we derive. 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 size n→∞ it matches with the 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 short times tΓ=O(1), the spectrum is in a Coulomb gas regime, with a well-defined continuous spectral distribution in the n→∞ limit. In that case, all moments of the density matrix become self-averaging and it is possible to exactly characterize the entanglement spectrum. 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 ∼O(Γt/n) play an essential role. We are still able to characterize the asymptotic behaviors of the entanglement entropy in this regime.

We consider a system of $N$ spinless fermions, interacting with each other via a power-law interaction $\epsilon/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 $\epsilon$, $E_{N}=E_{N}^{≤ft(0)}+\epsilon E_{N}^{≤ft(1)}+O≤ft(\epsilon^{2})$. We calculate $E_{N}^{≤ft(1)}$ 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}^{≤ft(1)}$, 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}^{≤ft(1)}$ significantly simplifies in the case where $n$ is even.

We consider critical one-dimensional quantum systems initially prepared in their ground state 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 two-point chiral correlation function can be obtained from solving two coupled ordinary stochastic differential equations. In contrast to the general expectation of heating, we demonstrate that, over the noise realizations, the system reaches a nontrivial and universal stationary distribution of states characterized by broad tails of physical quantities. As an example, we analyze the entanglement entropy production associated to a given interval of size ℓ. The corresponding stationary distribution has a 3/2 right tail for all ℓ and converges to a one-sided Levy stable for large ℓ. We obtain a similar result for 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 show that this stationary distribution for the energy density emerges even if the initial state is prepared at finite temperature. Our results are benchmarked via analytical and numerical calculations for a chain of noninteracting spinless fermions with excellent agreement.

We present a systematic short time expansion for the generating function of the one point height probability distribution for the KPZ equation with droplet initial condition, which goes much beyond previous studies. The expansion is checked against a numerical evaluation of the known exact Fredholm determinant expression. We also obtain the next order term for the Brownian initial condition. Although initially devised for short time, a resummation of the series allows to obtain also the long time large deviation function, found to agree with previous works using completely different techniques. Unexpected similarities with stationary large deviations of TASEP with periodic and open boundaries are discussed. Two additional applications are given. (i) Our method is generalized to study the linear statistics of the Airy point process, i.e. of the GUE edge eigenvalues. We obtain the generating function of the cumulants of the empirical measure to a high order. The second cumulant is found to match the result in the bulk obtained from the Gaussian free field by Borodin and Ferrari [1,2], but we obtain systematic corrections to the Gaussian free field (higher cumulants, expansion towards the edge). This also extends a result of Basor and Widom [3] to a much higher order. We obtain large deviation functions for the Airy point process for a variety of linear statistics test functions. (ii) We obtain results for the counting statistics of trapped fermions at the edge of the Fermi gas in both the high and the low temperature limits.

We calculate the effective action for disordered elastic manifolds in the ground state (equilibrium) up to 3-loop order. This yields the renormalization-group β-function to third order in ε=4−d, in an expansion in the dimension d around the upper critical dimension d=4. The calculations are performed using exact RG, and several other techniques, which allow us to resolve consistently the problems associated with the cusp of the renormalized disorder.

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 $\beta$. In the fermion model $\beta$ controls the strength of the interaction, $\beta=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 $\beta$ the variance of ${\cal N}_{\cal D}$ grows as $A_{\beta} \log N + B_{\beta}$ for $N \gg 1$ and we obtain a formula for $A_\beta$ and $B_\beta$. This is based on an explicit calculation for $\beta\in\left\{ 1,2,4\right\} $ and on a conjecture that we formulate for general $\beta$. 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/\beta$, 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 $\beta$. In some cases the fermion density exhibits interesting transitions, for example we obtain a noninteracting fermion formulation of the Gross-Witten-Wadia model.

The O-acetylated form of GD2, almost exclusively expressed in cancerous tissues, is considered to be a promising therapeutic target for neuroectoderm-derived tumors, especially for breast cancer. Our recent data have shown that 9-O-acetylated GD2 (9-OAcGD2) is the major O-acetylated ganglioside species in breast cancer cells. In 2015, Baumann et al. proposed that Cas 1 domain containing 1 (CASD1), which is the only known human sialyl-O-acetyltransferase, plays a role in GD3 O-acetylation. However, the mechanisms of ganglioside O-acetylation remain poorly understood. The aim of this study was to determine the involvement of CASD1 in GD2 O-acetylation in breast cancer. The role of CASD1 in OAcGD2 synthesis was first demonstrated using wild type CHO and CHOΔCasd1 cells as cellular models. Overexpression using plasmid transfection and siRNA strategies was used to modulate CASD1 expression in SUM159PT breast cancer cell line. Our results showed that OAcGD2 expression was reduced in SUM159PT that was transiently depleted for CASD1 expression. Additionally, OAcGD2 expression was increased in SUM159PT cells transiently overexpressing CASD1. The modulation of CASD1 expression using transient transfection strategies provided interesting insights into the role of CASD1 in OAcGD2 and OAcGD3 biosynthesis, and it highlights the importance of further studies on O-acetylation mechanisms.

Abstract Background Magnesium sulphate (MgSO4) is the only treatment approved for fetal neuroprotection. No information on its use is available in the absence of a national registry of neonatal practices. The objective of our study was to evaluate the use of MgSO4 for fetal neuroprotection in French tertiary maternity hospitals (FTMH). Methods Online and phone survey of all FTMH between August 2014 and May 2015. A participation was expected from one senior obstetrician, one senior anaesthetist and one senior neonatologist from each FTMH. Information was obtained from 63/63 (100%) FTMH and 138/189 (73%) physicians. Use of MgSO4 for fetal neuroprotection, regimen and injection protocols, reasons for non-use were the main outcome measures. Results 60.3% of FTMH used MgSO4 for fetal neuroprotection. No significant difference was observed between university and non-university hospitals or according to the annual number of births. Protocols differed especially in terms of the maximum gestational age (3%

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\'ech\'e (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.

We consider the Kardar-Parisi-Zhang (KPZ) equation for the stochastic growth of an interface of height $h(x,t)$ on the positive half line with boundary condition $\partial_x h(x,t)|_{x=0}=A$. It is equivalent to a continuum directed polymer (DP) in a random potential in half-space with a wall at $x=0$ either repulsive $A>0$, or attractive $A - \frac{1}{2}$ the large time PDF is the GSE Tracy-Widom distribution. For $A= \frac{1}{2}$, the critical point at which the DP binds to the wall, we obtain the GOE Tracy-Widom distribution. In the critical region, $A+\frac{1}{2} = \epsilon t^{-1/3} \to 0$ with fixed $\epsilon = \mathcal{O}(1)$, we obtain a transition kernel continuously depending on $\epsilon$. Our work extends the results obtained previously for $A=+\infty$, $A=0$ and $A=- \frac{1}{2}$.

The expression and biological functions of oncofetal markers GD2 and GD3 were extensively studied in neuroectoderm-derived cancers in order to characterize their potential as therapeutic targets. Using immunological approaches, we previously identified GD3, GD2, and OAcGD2 expression in breast cancer (BC) cell lines. However, antibodies specific for O-acetylated gangliosides are not exempt of limitations, as they only provide information on the expression of a limited set of O-acetylated ganglioside species. Consequently, the aim of the present study was to use structural approaches in order to apprehend ganglioside diversity in melanoma, neuroblastoma, and breast cancer cells, focusing on O-acetylated species that are usually lost under alkaline conditions and require specific analytical procedures. We used purification and extraction methods that preserve the O-acetyl modification for the analysis of native gangliosides by MALDI-TOF. We identified the expression of GM1, GM2, GM3, GD2, GD3, GT2, and GT3 in SK-Mel28 (melanoma), LAN-1 (neuroblastoma), Hs 578T, SUM 159PT, MDA-MB-231, MCF-7 (BC), and BC cell lines over-expressing GD3 synthase. Among O-acetylated gangliosides, we characterized the expression of OAcGM1, OAcGD3, OAcGD2, OAcGT2, and OAcGT3. Furthermore, the experimental procedure allowed us to clearly identify the position of the sialic acid residue that carries the O-acetyl group on b- and c-series gangliosides by MS/MS fragmentation. These results show that ganglioside O-acetylation occurs on both inner and terminal sialic acid residue in a cell type-dependent manner, suggesting different O-acetylation pathways for gangliosides. They also highlight the limitation of immuno-detection for the complete identification of O-acetylated ganglioside profiles in cancer cells.

We investigate the short-time regime of the KPZ equation in $1+1$ dimensions and develop a unifying method to obtain the height distribution in this regime, valid whenever an exact solution exists in the form of a Fredholm Pfaffian or determinant. These include the droplet and stationary initial conditions in full space, previously obtained by a different method. The novel results concern the droplet initial condition in a half space for several Neumann boundary conditions: hard wall, symmetric, and critical. In all cases, the height probability distribution takes the large deviation form $P(H,t) \sim \exp( - \Phi(H)/\sqrt{t})$ for small time. We obtain the rate function $\Phi(H)$ analytically for the above cases. It has a Gaussian form in the center with asymmetric tails, $|H|^{5/2}$ on the negative side, and $H^{3/2}$ on the positive side. The amplitude of the left tail for the half-space is found to be half the one of the full space. As in the full space case, we find that these left tails remain valid at all times. In addition, we present here (i) a new Fredholm Pfaffian formula for the solution of the hard wall boundary condition and (ii) two Fredholm determinant representations for the solutions of the hard wall and the symmetric boundary respectively.

Building upon the one-step replica symmetry breaking formalism, duly understood and ramified, we show that the sequence of ordered extreme values of a general class of Euclidean-space logarithmically correlated random energy models (logREMs) behave in the thermodynamic limit as a randomly shifted decorated exponential Poisson point process. The distribution of the random shift is determined solely by the large-distance ("infra-red", IR) limit of the model, and is equal to the free energy distribution at the critical temperature up to a translation. the decoration process is determined solely by the small-distance ("ultraviolet", UV) limit, in terms of the biased minimal process. Our approach provides connections of the replica framework to results in the probability literature and sheds further light on the freezing/duality conjecture which was the source of many previous results for log-REMs. In this way we derive the general and explicit formulae for the joint probability density of depths of the first and second minima (as well its higher-order generalizations) in terms of model-specific contributions from UV as well as IR limits. In particular, we show that the second min statistics is largely independent of details of UV data, whose influence is seen only through the mean value of the gap. For a given log-correlated field this parameter can be evaluated numerically, and we provide several numerical tests of our theory using the circular model of $1/f$-noise.

We analyse by numerical simulations and scaling arguments the avalanche statistics of 1-dimensional elastic interfaces in random media driven at a single point. Both global and local avalanche sizes are power-law distributed, with universal exponents given by the depinning roughness exponent ζ and the interface dimension d, and distinct from their values in the uniformly driven case. A crossover appears between uniformly driven behaviour for small avalanches, and point-driven behaviour for large avalanches. The scale of the crossover is controlled by the ratio between the stiffness of the pulling spring and the elasticity of the interface; it is visible both in the global and local avalanche-size distributions, as in the average spatial avalanche shape. Our results are relevant to model experiments involving locally driven elastic manifolds at low temperatures, such as magnetic domain walls or vortex lines in superconductors. [ABSTRACT FROM AUTHOR]

We investigate the statistics of the local time T=∫&#95;{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(T) satisfies the large deviation principle P(T)∼e^{-TI(T/T)} in the large observation time limit T→∞. Remarkably, we find that in the presence of drift the rate function I(ρ) is nonanalytic: we interpret its singularity as a dynamical phase transition of first order. We then extend these results by studying the statistics of the amount of time 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(R=T) that the particle does not exit the interval up to time T. We find that the conditional end-point 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.

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)∼k^{μ-1} for k→0. 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 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 scaling predictions that partially disagree with previous works, but that are corroborated by numerical simulations. We also provide an interpretation of the anomalous scaling in terms of the abrupt jumps in the line's shape that dominate the average value of the observable.

We investigate and solve the weak noise theory for the semidiscrete O'Connell-Yor directed polymer. In the large deviation regime, the most probable evolution of the partition function obeys a classical nonlinear system which is a nonstandard discretization of the nonlinear Schrödinger 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 non-standard Fredholm determinant framework that we develop. This allows us 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.

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≫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.

We study the dynamics of the outliers for a large number of independent Brownian particles in one dimension. We derive the multitime 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 we study the counting statistics at the edge. Finally, we derive the multitime joint distribution of the running maximum, as well as the joint distribution of the arrival time of the first particle at several space points.

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 semicircle 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≫1, the covariance matrix takes scaling forms that we compute explicitly both in the bulk and at the edge of the support of the semicircle. 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 obtained by Agarwal et al. [J. Stat. Phys. 176, 1463 (2019)0022-471510.1007/s10955-019-02349-6]. 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 in Gorin and Kleptsyn [arXiv:2009.02006 (2023)] 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.

We study the nonequilibrium 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 long 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 us 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, which 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.

We study the crossover from the macroscopic fluctuation theory (MFT), which describes one-dimensional 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 nonlinear system which interpolates between the derivative and the standard nonlinear 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.

We consider an overdamped run-and-tumble particle in two dimensions, with self-propulsion in an orientation that stochastically rotates by 90^{∘} at a constant rate, clockwise or counterclockwise 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(x,y,t) of the particle's position, and in particular, the steady-state distribution P&#95;{st}(x,y) that is reached in the long-time limit. We also find P(x,y,t) 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.

In the 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 and the distribution of their sizes and of the distances between them. We also discuss applications to depinning avalanches with long-range elasticity.

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 nonlinear hydrodynamic equations of the WNT are solved analytically through a connection 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 nonvanishing 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.

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 one 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 one 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 one to study the final equilibrium at large N.

We study two interacting identical run-and-tumble particles (RTPs) 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 RTPs 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)∼|y|. In this case, for D=0 we find that P(y) contains a δ function part at y=0, signaling a strong clustering effect, together with a smooth exponential component. For D>0, the δ function part broadens, leading instead to weak clustering. The second example is the harmonic attraction V(y)∼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.

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 nonlinear 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.

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.

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.

We study the position distribution P(R[over ⃗],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(R[over ⃗],N)→P(R,N) where R=|R[over ⃗]|. 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&#95;{c}∼O(N) located in the large-deviation regime of P(R,N). For RR&#95;{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&#95;{0})^{α-1}/v&#95;{0}, parametrized by α>0, we show that, for large N, P(R,N)∼exp[-Nψ&#95;{d,α}(R/N)], and we compute exactly the rate function ψ&#95;{d,α}(z) for any d and α. We show that the transition manifests itself as a singularity of this rate function at R=R&#95;{c} and that its order depends continuously on d and α. We also compute the distribution of the condensate size for R>R&#95;{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 ∼10^{-100}.

We revisit the problem of an elastic line (such as 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 in the mathematics literature, and their connections to random matrix theory, we find that for a single line and a single strong defect this transition in the presence of point disorder coincides with the Baik-Ben Arous-Péché (BBP) transition for the appearance of outliers in the spectrum of a perturbed random matrix in the Gaussian unitary ensemble. 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&#95;{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 or (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, closely related to the kernel describing the largest real eigenvalues of the real Ginibre ensemble. The case of many columns and many nonintersecting 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. Connections with recent results on the generalized Rosenzweig-Porter model suggest that the localization of many polymers occurs gradually upon increasing their lengths.

We develop a first-principles 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 N&#95;{D} in a domain D of macroscopic size in the bulk of the support. We show that the variance of N&#95;{D} grows as N^{(d-1)/d}(A&#95;{d}logN+B&#95;{d}) for large N, and obtain the explicit dependence of A&#95;{d},B&#95;{d} on the potential and on the size of 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 D, in any dimension, supported by exact results for d=1.

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 Letter 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 number of clusters. Our work has concrete importance for experimental applications where the cluster statistics is a key probe of avalanche dynamics.

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 noninstantaneous. 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)10.1103/PhysRevLett.124.090603].