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

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.

Large deviations of the interface height in the Golubović-Bruinsma model of stochastic growth

We study large deviations of the one-point height H of a stochastic interface, governed by the Golubović-Bruinsma equation, ∂_{t}h=-ν∂_{x}^{4}h+(λ/2)(∂_{x}h)^{2}+sqrt[D]ξ(x,t), where h(x,t) is the interface height at point x and time t and ξ(x,t) is the Gaussian white noise. The interface is initially flat, and H is defined by the relation h(x=0,t=T)=H. We focus on the short-time limit, T≪T_{NL}, where T_{NL}=ν^{5/7}(Dλ^{2})^{-4/7} is the characteristic nonlinear time of the system. In this limit typical, small fluctuations of H are unaffected by the nonlinear term, and they are Gaussian. However, the large-deviation tails of the probability distribution P(H,T) "feel" the nonlinearity already at short times, and they are non-Gaussian and asymmetric. We evaluate these tails using the optimal fluctuation method (OFM). The lower tail scales as -lnP(H,T)∼H^{5/2}/T^{1/2}. It coincides with its analog for the Kardar-Parisi-Zhang (KPZ) equation, and we point out to the mechanism of this universality. The upper tail scales as -lnP(H,T)∼H^{11/6}/T^{5/6}, it is different from the upper tail of the KPZ equation. We also compute the large deviation function of H numerically and verify our asymptotic results for the tails.

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.

Imitating nonequilibrium steady states using time-varying equilibrium force in many-body diffusive systems

An equivalence between nonequilibrium steady states (NESS) driven by a time-independent force and stochastic pumps (SP) stirred by a time-varying conservative force is studied for general many-body diffusive systems. When the particle density and current of NESS are imitated by SP time-averaged counterparts, we prove that the entropy production rate in the SP is always greater than that of the NESS, provided that the conductivity of the particle current is concave as a function of the particle density. Searching for a SP protocol that saturates the entropy production bound reveals an unexpected connection with traffic waves, where a high density region propagates against the direction of the particle current.

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.