Exact short-time height distribution in the one-dimensional Kardar-Parisi-Zhang equation with Brownian initial condition

Weak noise theory of the O'Connell-Yor polymer as an integrable discretization of the nonlinear Schrödinger equation

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.

Crossover from the macroscopic fluctuation theory to the Kardar-Parisi-Zhang equation controls the large deviations beyond Einstein's diffusion

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.

Large deviations in chaotic systems: Exact results and dynamical phase transition

Large deviations in chaotic dynamics have potentially significant and dramatic consequences. We study large deviations of series of finite lengths N generated by chaotic maps. The distributions generally display an exponential decay with N, associated with large-deviation (rate) functions. We obtain the exact rate functions analytically for the doubling, tent, and logistic maps. For the latter two, the solution is given as a power series whose coefficients can be systematically calculated to any order. We also obtain the rate function for the cat map numerically, uncovering strong evidence for the existence of a remarkable singularity of it that we interpret as a second-order dynamical phase transition. Furthermore, we develop a numerical tool for efficiently simulating atypical realizations of sequences if the chaotic map is not invertible, and we apply it to the tent and logistic maps.

Exact short-time height distribution and dynamical phase transition in the relaxation of a Kardar-Parisi-Zhang interface with random initial condition

We consider the relaxation (noise-free) statistics of the one-point height H=h(x=0,t), where h(x,t) is the evolving height of a one-dimensional Kardar-Parisi-Zhang (KPZ) interface, starting from a Brownian (random) initial condition. We find that, at short times, the distribution of H takes the same scaling form -lnP(H,t)=S(H)/sqrt[t] as the distribution of H for the KPZ interface driven by noise, and we find the exact large-deviation function S(H) analytically. At a critical value H=H_{c}, the second derivative of S(H) jumps, signaling a dynamical phase transition (DPT). Furthermore, we calculate exactly the most likely history of the interface that leads to a given H, and show that the DPT is associated with spontaneous breaking of the mirror symmetry x↔-x of the interface. In turn, we find that this symmetry breaking is a consequence of the nonconvexity of a large-deviation function that is closely related to S(H), and describes a similar problem but in half space. Moreover, the critical point H_{c} is related to the inflection point of the large-deviation function of the half-space problem.

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

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.

Observing symmetry-broken optimal paths of the stationary Kardar-Parisi-Zhang interface via a large-deviation sampling of directed polymers in random media

Consider the short-time probability distribution P(H,t) of the one-point interface height difference h(x=0,τ=t)-h(x=0,τ=0)=H of the stationary interface h(x,τ) described by the Kardar-Parisi-Zhang (KPZ) equation. It was previously shown that the optimal path, the most probable history of the interface h(x,τ) which dominates the upper tail of P(H,t), is described by any of two ramplike structures of h(x,τ) traveling either to the left, or to the right. These two solutions emerge, at a critical value of H, via a spontaneous breaking of the mirror symmetry x↔-x of the optimal path, and this symmetry breaking is responsible for a second-order dynamical phase transition in the system. We simulate the interface configurations numerically by employing a large-deviation Monte Carlo sampling algorithm in conjunction with the mapping between the KPZ interface and the directed polymer in a random potential at high temperature. This allows us to observe the optimal paths, which determine each of the two tails of P(H,t), down to probability densities as small as 10^{-500}. At short times we observe mirror-symmetry-broken traveling optimal paths for the upper tail, and a single mirror-symmetric path for the lower tail, in good quantitative agreement with analytical predictions. At long times, even at moderate values of H, where the optimal fluctuation method is not supposed to apply, we still observe two well-defined dominating paths. Each of them violates the mirror symmetry x↔-x and is a mirror image of the other.

Inverse Scattering of the Zakharov-Shabat System Solves the Weak Noise Theory of the Kardar-Parisi-Zhang Equation

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.

Probing large deviations of the Kardar-Parisi-Zhang equation at short times with an importance sampling of directed polymers in random media

The one-point distribution of the height for the continuum Kardar-Parisi-Zhang equation is determined numerically using the mapping to the directed polymer in a random potential at high temperature. Using an importance sampling approach, the distribution is obtained over a large range of values, down to a probability density as small as 10^{-1000} in the tails. The short-time behavior is investigated and compared with recent analytical predictions for the large-deviation forms of the probability of rare fluctuations, showing a spectacular agreement with the analytical expressions. The flat and stationary initial conditions are studied in the full space, together with the droplet initial condition in the half-space.

Large fluctuations of a Kardar-Parisi-Zhang interface on a half line: The height statistics at a shifted point

We consider a stochastic interface h(x,t), described by the 1+1 Kardar-Parisi-Zhang (KPZ) equation on the half line x≥0 with the reflecting boundary at x=0. The interface is initially flat, h(x,t=0)=0. We focus on the short-time probability distribution P(H,L,t) of the height H of the interface at point x=L. Using the optimal fluctuation method, we determine the (Gaussian) body of the distribution and the strongly asymmetric non-Gaussian tails. We find that the slower-decaying tail scales as -sqrt[t]lnP≃|H|^{3/2}f_{-}(L/sqrt[|H|t]) and calculate the function f_{-} analytically. Remarkably, this tail exhibits a first-order dynamical phase transition at a critical value of L, L_{c}=0.60223⋯sqrt[|H|t]. The transition results from a competition between two different fluctuation paths of the system. The faster-decaying tail scales as -sqrt[t]lnP≃|H|^{5/2}f_{+}(L/sqrt[|H|t]). We evaluate the function f_{+} using a specially developed numerical method which involves solving a nonlinear second-order elliptic equation in Lagrangian coordinates. The faster-decaying tail also involves a sharp transition which occurs at a critical value L_{c}≃2sqrt[2|H|t]/π. This transition is similar to the one recently found for the KPZ equation on a ring, and we believe that it has the same fractional order, 5/2. It is smoothed, however, by small diffusion effects.

Coulomb-Gas Electrostatics Controls Large Fluctuations of the Kardar-Parisi-Zhang Equation

We establish a large deviation principle for the Kardar-Parisi-Zhang (KPZ) equation, providing precise control over the left tail of the height distribution for narrow wedge initial condition. Our analysis exploits an exact connection between the KPZ one-point distribution and the Airy point process-an infinite particle Coulomb gas that arises at the spectral edge in random matrix theory. We develop the large deviation principle for the Airy point process and use it to compute, in a straightforward and assumption-free manner, the KPZ large deviation rate function in terms of an electrostatic problem (whose solution we evaluate). This method also applies to the half-space KPZ equation, showing that its rate function is half of the full-space rate function. In addition to these long-time estimates, we provide rigorous proof of finite-time tail bounds on the KPZ distribution, which demonstrate a crossover between exponential decay with exponent 3 (in the shallow left tail) to exponent 5/2 (in the deep left tail). The full-space KPZ rate function agrees with the one computed in Sasorov et al. [J. Stat. Mech. (2017) 063203JSMTC61742-546810.1088/1742-5468/aa73f8] via a WKB approximation analysis of a nonlocal, nonlinear integrodifferential equation generalizing Painlevé II which Amir et al. [Commun. Pure Appl. Math. 64, 466 (2011)CPMAMV0010-364010.1002/cpa.20347] related to the KPZ one-point distribution.

Exact short-time height distribution for the flat Kardar-Parisi-Zhang interface

We determine the exact short-time distribution -lnP_{f}(H,t)=S_{f}(H)/sqrt[t] of the one-point height H=h(x=0,t) of an evolving 1+1 Kardar-Parisi-Zhang (KPZ) interface for flat initial condition. This is achieved by combining (i) the optimal fluctuation method, (ii) a time-reversal symmetry of the KPZ equation in 1+1 dimension, and (iii) the recently determined exact short-time height distribution -lnP_{st}(H,t)=S_{st}(H)/sqrt[t] for stationary initial condition. In studying the large-deviation function S_{st}(H) of the latter, one encounters two branches: an analytic and a nonanalytic. The analytic branch is nonphysical beyond a critical value of H where a second-order dynamical phase transition occurs. Here we show that, remarkably, it is the analytic branch of S_{st}(H) which determines the large-deviation function S_{f}(H) of the flat interface via a simple mapping S_{f}(H)=2^{-3/2}S_{st}(2H).

Landau theory of the short-time dynamical phase transitions of the Kardar-Parisi-Zhang interface

We study the short-time distribution P(H,L,t) of the two-point two-time height difference H=h(L,t)-h(0,0) of a stationary Kardar-Parisi-Zhang interface in 1+1 dimension. Employing the optimal-fluctuation method, we develop an effective Landau theory for the second-order dynamical phase transition found previously for L=0 at a critical value H=H_{c}. We show that |H| and L play the roles of inverse temperature and external magnetic field, respectively. In particular, we find a first-order dynamical phase transition when L changes sign, at supercritical H. We also determine analytically P(H,L,t) in several limits away from the second-order transition. Typical fluctuations of H are Gaussian, but the distribution tails are highly asymmetric. The tails -lnP∼|H|^{3/2}/sqrt[t] and -lnP∼|H|^{5/2}/sqrt[t], previously found for L=0, are enhanced for L≠0. At very large |L| the whole height-difference distribution P(H,L,t) is time-independent and Gaussian in H, -lnP∼|H|^{2}/|L|, describing the probability of creating a ramplike height profile at t=0.