Extreme current fluctuations in lattice gases: beyond nonequilibrium steady states

Exact Solution of the Macroscopic Fluctuation Theory for the Symmetric Exclusion Process

We present the first exact solution for the time-dependent equations of the macroscopic fluctuation theory (MFT) for the symmetric simple exclusion process by combining a generalization of the canonical Cole-Hopf transformation with the inverse scattering method. For the step initial condition with two densities, we obtain exact and compact formulas for the optimal density profile and the response field that produce a required fluctuation, both at initial and final times. The large deviation function of the current is derived and coincides with the formula obtained previously by microscopic calculations. This provides the first analytic confirmation of the validity of the MFT for an interacting model in the time-dependent regime.

Inverse Scattering Method Solves the Problem of Full Statistics of Nonstationary Heat Transfer in the Kipnis-Marchioro-Presutti Model

We determine the full statistics of nonstationary heat transfer in the Kipnis-Marchioro-Presutti lattice gas model at long times by uncovering and exploiting complete integrability of the underlying equations of the macroscopic fluctuation theory. These equations are closely related to the derivative nonlinear Schrödinger equation (DNLS), and we solve them by the Zakharov-Shabat inverse scattering method (ISM) adapted by D. J. Kaup and A. C. Newell, J. Math. Phys. 19, 798 (1978)JMAPAQ0022-248810.1063/1.523737 for the DNLS. We obtain explicit results for the exact large deviation function of the transferred heat for an initially localized heat pulse, where we uncover a nontrivial symmetry relation.

Solitons in fluctuating hydrodynamics of diffusive processes

We demonstrate that fluid mechanical systems arising from large fluctuations of one-dimensional statistical processes generically exhibit solitons and nonlinear waves. We derive the explicit form of these solutions and examine their properties for the specific cases of the Kipnis-Marchioro-Presutti model (KMP) and the symmetric exclusion process (SEP). We show that the two fluid systems are related by a nonlinear transformation but still have markedly different properties. In particular, the KMP fluid has a nontrivial sound wave spectrum exhibiting birefringence, whereas sound waves for the SEP fluid are essentially trivial. The appearance of sound waves and soliton configurations in the KMP model is related to the onset of instabilities.

Structure of the optimal path to a fluctuation

Macroscopic fluctuations have become an essential tool to understand physics far from equilibrium due to the link between their statistics and nonequilibrium ensembles. The optimal path leading to a fluctuation encodes key information on this problem, shedding light on, e.g., the physics behind the enhanced probability of rare events out of equilibrium, the possibility of dynamic phase transitions, and new symmetries. This makes the understanding of the properties of these optimal paths a central issue. Here we derive a fundamental relation which strongly constrains the architecture of these optimal paths for general d-dimensional nonequilibrium diffusive systems, and implies a nontrivial structure for the dominant current vector fields. Interestingly, this general relation (which encompasses and explains previous results) makes manifest the spatiotemporal nonlocality of the current statistics and the associated optimal trajectories.

Dynamical Symmetry Breaking and Phase Transitions in Driven Diffusive Systems

We study the probability distribution of a current flowing through a diffusive system connected to a pair of reservoirs at its two ends. Sufficient conditions for the occurrence of a host of possible phase transitions both in and out of equilibrium are derived. These transitions manifest themselves as singularities in the large deviation function, resulting in enhanced current fluctuations. Microscopic models which implement each of the scenarios are presented, with possible experimental realizations. Depending on the model, the singularity is associated either with a particle-hole symmetry breaking, which leads to a continuous transition, or in the absence of the symmetry with a first-order phase transition. An exact Landau theory which captures the different singular behaviors is derived.

Short-time height distribution in the one-dimensional Kardar-Parisi-Zhang equation: Starting from a parabola

We study the probability distribution P(H,t,L) of the surface height h(x=0,t)=H in the Kardar-Parisi-Zhang (KPZ) equation in 1+1 dimension when starting from a parabolic interface, h(x,t=0)=x^{2}/L. The limits of L→∞ and L→0 have been recently solved exactly for any t>0. Here we address the early-time behavior of P(H,t,L) for general L. We employ the weak-noise theory-a variant of WKB approximation-which yields the optimal history of the interface, conditioned on reaching the given height H at the origin at time t. We find that at small HP(H,t,L) is Gaussian, but its tails are non-Gaussian and highly asymmetric. In the leading order and in a proper moving frame, the tails behave as -lnP=f_{+}|H|^{5/2}/t^{1/2} and f_{-}|H|^{3/2}/t^{1/2}. The factor f_{+}(L,t) monotonically increases as a function of L, interpolating between time-independent values at L=0 and L=∞ that were previously known. The factor f_{-} is independent of L and t, signaling universality of this tail for a whole class of deterministic initial conditions.

Weak additivity principle for current statistics in d dimensions

The additivity principle (AP) allows one to compute the current distribution in many one-dimensional nonequilibrium systems. Here we extend this conjecture to general d-dimensional driven diffusive systems, and validate its predictions against both numerical simulations of rare events and microscopic exact calculations of three paradigmatic models of diffusive transport in d=2. Crucially, the existence of a structured current vector field at the fluctuating level, coupled to the local mobility, turns out to be essential to understand current statistics in d>1. We prove that, when compared to the straightforward extension of the AP to high d, the so-called weak AP always yields a better minimizer of the macroscopic fluctuation theory action for current statistics.

Large Deviations of Surface Height in the Kardar-Parisi-Zhang Equation

Using the weak-noise theory, we evaluate the probability distribution P(H,t) of large deviations of height H of the evolving surface height h(x,t) in the Kardar-Parisi-Zhang equation in one dimension when starting from a flat interface. We also determine the optimal history of the interface, conditioned on reaching the height H at time t. We argue that the tails of P behave, at arbitrary time t>0, and in a proper moving frame, as -lnP∼|H|^{5/2} and ∼|H|^{3/2}. The 3/2 tail coincides with the asymptotic of the Gaussian orthogonal ensemble Tracy-Widom distribution, previously observed at long times.

Macroscopic fluctuation theory and first-passage properties of surface diffusion

We investigate nonequilibrium fluctuations of a solid surface governed by the stochastic Mullins-Herring equation with conserved noise. This equation describes surface diffusion of adatoms accompanied by their exchange between the surface and the bulk of the solid, when desorption of adatoms is negligible. Previous works dealt with dynamic scaling behavior of the fluctuating interface. Here we determine the probability that the interface first reaches a large given height at a specified time. We also find the optimal time history of the interface conditional on this nonequilibrium fluctuation. We obtain these results by developing a macroscopic fluctuation theory of surface diffusion.

Survival of interacting diffusing particles inside a domain with absorbing boundary

Suppose that a d-dimensional domain is filled with a gas of (in general, interacting) diffusive particles with density n_{0}. A particle is absorbed whenever it reaches the domain boundary. Employing macroscopic fluctuation theory, we evaluate the probability P that no particles are absorbed during a long time T. We argue that the most likely gas density profile, conditional on this event, is stationary throughout most of the time T. As a result, P decays exponentially with T for a whole class of interacting diffusive gases in any dimension. For d=1 the stationary gas density profile and P can be found analytically. In higher dimensions we focus on the simple symmetric exclusion process (SSEP) and show that -lnP≃D_{0}TL^{d-2}s(n_{0}), where D_{0} is the gas diffusivity, and L is the linear size of the system. We calculate the rescaled action s(n_{0}) for d=1, for rectangular domains in d=2, and for spherical domains. Near close packing of the SSEP s(n_{0}) can be found analytically for domains of any shape and in any dimension.

Excess white noise to probe transport mechanisms in a membrane channel

Current fluctuation analysis has been successfully used over the years to investigate the physical properties of different systems. Here, we perform single-channel time-resolved current experiments in a protein channel to evaluate the different transport mechanisms governing the channel function. Using different salts of monovalent and divalent cations in a wide range of concentrations and applied potentials, we analyze current fluctuations focusing on the voltage dependence of the additional white noise that appears in the low-frequency range of the spectra. We demonstrate that the channel displays two characteristic transport regimes: at low salt concentrations (10 mM to 1 M) ion permeation is controlled by the protein fixed charges that induce accumulation or exclusion of ions to preserve local electroneutrality. At high salt concentrations (>1M) adsorption processes associated to the binding of cations to the channel charges regulate the transport properties.

Survival of a static target in a gas of diffusing particles with exclusion

Let a lattice gas of constant density, described by the symmetric simple exclusion process, be brought in contact with a "target": a spherical absorber of radius R. Employing the macroscopic fluctuation theory (MFT), we evaluate the probability P(T) that no gas particle hits the target until a long but finite time T. We also find the most likely gas density history conditional on the nonhitting. The results depend on the dimension of space d and on the rescaled parameter ℓ=R/√[D(0)T], where D(0) is the gas diffusivity. For small ℓ and d>2, P(T) is determined by an exact stationary solution of the MFT equations that we find. For large ℓ, and for any ℓ in one dimension, the relevant MFT solutions are nonstationary. In this case, lnP(T) scales differently with relevant parameters, and it also depends on whether the initial condition is random or deterministic. The latter effects also occur if the lattice gas is composed of noninteracting random walkers. Finally, we extend the formalism to a whole class of diffusive gases of interacting particles.