Dynamical Symmetry Breaking and Phase Transitions in Driven Diffusive Systems

Variational Path Sampling of Rare Dynamical Events

This article reviews the concepts and methods of variational path sampling. These methods allow computational studies of rare events in systems driven arbitrarily far from equilibrium. Based upon a statistical mechanics of trajectory space and leveraging the theory of large deviations, they provide a perspective from which dynamical phenomena can be studied with the same types of ensemble reweighting ideas that have been used for static equilibrium properties. Applications to chemical, material, and biophysical systems are highlighted.

Programmable time crystals from higher-order packing fields

Time crystals are many-body systems that spontaneously break time-translation symmetry, and thus exhibit long-range spatiotemporal order and robust periodic motion. Recent results have shown that coupling an external packing field to density fluctuations in driven diffusive fluids can trigger a transition to a time-crystal phase. Here, we exploit this mechanism to engineer and control on-demand programmable continuous time crystals characterized by an arbitrary number of rotating condensates, which can be further enhanced with higher-order modes. We elucidate the underlying critical point, as well as general properties of the condensates' density profiles and velocities, demonstrating a scaling property of higher-order traveling condensates in terms of first-order ones. We illustrate our findings by solving the hydrodynamic equations for various paradigmatic driven diffusive systems, obtaining along the way a number of remarkable results, e.g., the possibility of explosive time crystal phases characterized by an abrupt, first-order-type transition. Overall, these results demonstrate the versatility and broad possibilities of this promising route to time crystals.

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

We investigate the statistics of the local time T=∫_{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.

Emerging universality classes in thermally assisted activation of interacting diffusive systems: A perturbative hydrodynamic approach

Thermal activation of a particle from a deep potential trap follows the Arrhenius law. Recently, this result has been generalized for interacting diffusive particles in the trap, revealing two universality classes-the Arrhenius class and the excluded volume class. The result was demonstrated with the aid of numerical analysis. Here, we present a perturbative hydrodynamic approach to analytically validate the existence and range of validity for the two universality classes.

Large deviations in statistics of the convex hull of passive and active particles: A theoretical study

We investigate analytically the distribution tails of the area A and perimeter L of a convex hull for different types of planar random walks. For N noninteracting Brownian motions of duration T we find that the large-L and -A tails behave as P(L)∼e^{-b_{N}L^{2}/DT} and P(A)∼e^{-c_{N}A/DT}, while the small-L and -A tails behave as P(L)∼e^{-d_{N}DT/L^{2}} and P(A)∼e^{-e_{N}DT/A}, where D is the diffusion coefficient. We calculated all of the coefficients (b_{N},c_{N},d_{N},e_{N}) exactly. Strikingly, we find that b_{N} and c_{N} are independent of N for N≥3 and N≥4, respectively. We find that the large-L (A) tails are dominated by a single, most probable realization that attains the desired L (A). The left tails are dominated by the survival probability of the particles inside a circle of appropriate size. For active particles and at long times, we find that large-L and -A tails are given by P(L)∼e^{-TΨ_{N}^{per}(L/T)} and P(A)∼e^{-TΨ_{N}^{area}(sqrt[A]/T)}, respectively. We calculate the rate functions Ψ_{N} exactly and find that they exhibit multiple singularities. We interpret these as DPTs of first order. We extended several of these results to dimensions d>2. Our analytic predictions display excellent agreement with existing results that were obtained from extensive numerical simulations.

Arrhenius law for interacting diffusive systems

Finding the mean time it takes for a particle to escape from a metastable state due to thermal fluctuations is a fundamental problem in physics, chemistry, and biology. Here, we consider the escape rate of interacting diffusive particles, from a deep potential trap within the framework of the macroscopic fluctuation theory-a nonequilibrium hydrodynamic theory. For systems without excluded volume, our investigation reveals adherence to the well-established Arrhenius law. However, in the presence of excluded volume, a universality class emerges, fundamentally altering the escape rate. Remarkably, the modified escape rate within this universality class is independent of the interactions at play. The universality class, demonstrating the importance of excluded volume effects, may bring insights to the interpretation of escape processes in the realm of chemical physics.

Finding the Effective Dynamics to Make Rare Events Typical in Chaotic Maps

Dynamical fluctuations or rare events associated with atypical trajectories in chaotic maps due to specific initial conditions can crucially determine their fate, as the may lead to stability islands or regions in phase space otherwise displaying unusual behavior. Yet, finding such initial conditions is a daunting task precisely because of the chaotic nature of the system. In this Letter, we circumvent this problem by proposing a framework for finding an effective topologically conjugate map whose typical trajectories correspond to atypical ones of the original map. This is illustrated by means of examples which focus on counterbalancing the instability of fixed points and periodic orbits, as well as on the characterization of a dynamical phase transition involving the finite-time Lyapunov exponent. The procedure parallels that of the application of the generalized Doob transform in the stochastic dynamics of Markov chains, diffusive processes, and open quantum systems, which in each case results in a new process having the prescribed statistics in its stationary state. This Letter thus brings chaotic maps into the growing family of systems whose rare fluctuations-sustaining prescribed statistics of dynamical observables-can be characterized and controlled by means of a large-deviation formalism.

Tricritical Behavior in Dynamical Phase Transitions

We identify a new scenario for dynamical phase transitions associated with time-integrated observables occurring in diffusive systems described by the macroscopic fluctuation theory. It is characterized by the pairwise meeting of first- and second-order bias-induced phase transition curves at two tricritical points. We formulate a simple, general criterion for its appearance and derive an exact Landau theory for the tricritical behavior. The scenario is demonstrated in three examples: the simple symmetric exclusion process biased by an activity-related structural observable; the Katz-Lebowitz-Spohn lattice gas model biased by its current; and in an active lattice gas biased by its entropy production.

Spectral signatures of symmetry-breaking dynamical phase transitions

Large deviation theory provides the framework to study the probability of rare fluctuations of time-averaged observables, opening new avenues of research in nonequilibrium physics. Some of the most appealing results within this context are dynamical phase transitions (DPTs), which might occur at the level of trajectories in order to maximize the probability of sustaining a rare event. While macroscopic fluctuation theory has underpinned much recent progress on the understanding of symmetry-breaking DPTs in driven diffusive systems, their microscopic characterization is still challenging. In this work we shed light on the general spectral mechanism giving rise to continuous DPTs not only for driven diffusive systems, but for any jump process in which a discrete Z_{n} symmetry is broken. By means of a symmetry-aided spectral analysis of the Doob-transformed dynamics, we provide the conditions whereby symmetry-breaking DPTs might emerge and how the different dynamical phases arise from the specific structure of the degenerate eigenvectors. In particular, we show explicitly how all symmetry-breaking features are encoded in the subleading eigenvectors of the degenerate subspace. Moreover, by partitioning configuration space into equivalence classes according to a proper order parameter, we achieve a substantial dimensional reduction which allows for the quantitative characterization of the spectral fingerprints of DPTs. We illustrate our predictions in several paradigmatic many-body systems, including (1) the one-dimensional boundary-driven weakly asymmetric exclusion process (WASEP), which exhibits a particle-hole symmetry-breaking DPT for current fluctuations, (2) the three- and four-state Potts model for spin dynamics, which displays discrete rotational symmetry-breaking DPTs for energy fluctuations, and (3) the closed WASEP which presents a continuous symmetry-breaking DPT into a time-crystal phase characterized by a rotating condensate.

Dynamical phase transition in the occupation fraction statistics for noncrossing Brownian particles

We consider a system of N noncrossing Brownian particles in one dimension. We find the exact rate function that describes the long-time large deviation statistics of their occupation fraction in a finite interval in space. Remarkably, we find that, for any general N≥2, the system undergoes N-1 dynamical phase transitions of second order. The N-1 transitions are the boundaries of N phases that correspond to different numbers of particles which are in the vicinity of the interval throughout the dynamics. We achieve this by mapping the problem to that of finding the ground-state energy for N noninteracting spinless fermions in a square-well potential. The phases correspond to different numbers of single-body bound states for the quantum problem. We also study the process conditioned on a given occupation fraction and the large-N limiting behavior.

Exact spatial correlations in single-file diffusion

Single-file diffusion refers to the motion of diffusive particles in narrow channels, so that they cannot bypass each other. This constraint leads to the subdiffusion of a tagged particle, called the tracer. This anomalous behavior results from the strong correlations that arise in this geometry between the tracer and the surrounding bath particles. Despite their importance, these bath-tracer correlations have long remained elusive, because their determination is a complex many-body problem. Recently, we have shown that, for several paradigmatic models of single-file diffusion such as the simple exclusion process, these bath-tracer correlations obey a simple exact closed equation. In this paper, we provide the full derivation of this equation, as well as an extension to another model of single-file transport: the double exclusion process. We also make the connection between our results and the ones obtained very recently by several other groups and which rely on the exact solution of different models obtained by the inverse scattering method.

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.

Anomalous fluctuations of renewal-reward processes with heavy-tailed distributions

For renewal-reward processes with a power-law decaying waiting time distribution, anomalously large probabilities are assigned to atypical values of the asymptotic processes. Previous works have revealed that this anomalous scaling causes a singularity in the corresponding large deviation function. In order to further understand this problem, we study in this article the scaling of variance in several renewal-reward processes: counting processes with two different power-law decaying waiting time distributions and a Knudsen gas (a heat conduction model). Through analytical and numerical analyses of these models, we find that the variances show an anomalous scaling when the exponent of the power law is -3. For a counting process with the power-law exponent smaller than -3, this anomalous scaling does not take place: this indicates that if we only consider the standard deviation from the expectation, any anomalous behavior will not be detected. In this case, we argue that anomalous scaling appears in higher order cumulants. Finally, many-body particles interacting through soft-core interactions with the boundary conditions employed in the Knudsen gas are studied using numerical simulations. We observe that the variance scaling becomes normal even though the power-law exponent in the boundary conditions is -3.

Dynamical phase transition to localized states in the two-dimensional random walk conditioned on partial currents

The study of dynamical large deviations allows for a characterization of stationary states of lattice gas models out of equilibrium conditioned on averages of dynamical observables. The application of this framework to the two-dimensional random walk conditioned on partial currents reveals the existence of a dynamical phase transition between delocalized band dynamics and localized vortex dynamics. We present a numerical microscopic characterization of the phases involved and provide analytical insight based on the macroscopic fluctuation theory. A spectral analysis of the microscopic generator shows that the continuous phase transition is accompanied by spontaneous Z_{2}-symmetry breaking whereby the stationary solution loses the reflection symmetry of the generator. Dynamical phase transitions similar to this one, which do not rely on exclusion effects or interactions, are likely to be observed in more complex nonequilibrium physics models.

Long-ranged correlations in large deviations of local clustering

In systems of diffusing particles, we investigate large deviations of a time-averaged measure of clustering around one particle. We focus on biased ensembles of trajectories, which realize large-deviation events. The bias acts on a single particle, but elicits a response that spans the whole system. We analyze this effect through the lens of macroscopic fluctuation theory, focusing on the coupling of the bias to hydrodynamic modes. This explains that the dynamical free energy has nontrivial scaling relationships with the system size, in 1 and 2 spatial dimensions. We show that the long-ranged response to a bias on one particle also has consequences when biasing two particles.

Building Continuous Time Crystals from Rare Events

Symmetry-breaking dynamical phase transitions (DPTs) abound in the fluctuations of nonequilibrium systems. Here, we show that the spectral features of a particular class of DPTs exhibit the fingerprints of the recently discovered time-crystal phase of matter. Using Doob's transform as a tool, we provide a mechanism to build classical time-crystal generators from the rare event statistics of some driven diffusive systems. An analysis of the Doob's smart field in terms of the order parameter of the transition then leads to the time-crystal lattice gas (TCLG), a model of driven fluid subject to an external packing field, which presents a clear-cut steady-state phase transition to a time-crystalline phase characterized by a matter density wave, which breaks continuous time-translation symmetry and displays rigidity and long-range spatiotemporal order, as required for a time crystal. A hydrodynamic analysis of the TCLG transition uncovers striking similarities, but also key differences, with the Kuramoto synchronization transition. Possible experimental realizations of the TCLG in colloidal fluids are also discussed.

Using Matrix Product States to Study the Dynamical Large Deviations of Kinetically Constrained Models

Here we demonstrate that tensor network techniques-originally devised for the analysis of quantum many-body problems-are well suited for the detailed study of rare event statistics in kinetically constrained models (KCMs). As concrete examples, we consider the Fredrickson-Andersen and East models, two paradigmatic KCMs relevant to the modeling of glasses. We show how variational matrix product states allow us to numerically approximate-systematically and with high accuracy-the leading eigenstates of the tilted dynamical generators, which encode the large deviation statistics of the dynamics. Via this approach, we can study system sizes beyond what is possible with other methods, allowing us to characterize in detail the finite size scaling of the trajectory-space phase transition of these models, the behavior of spectral gaps, and the spatial structure and "entanglement" properties of dynamical phases. We discuss the broader implications of our results.

Single-file transport in periodic potentials: The Brownian asymmetric simple exclusion process

Single-file Brownian motion in periodic structures is an important process in nature and technology, which becomes increasingly amenable for experimental investigation under controlled conditions. To explore and understand generic features of this motion, the Brownian asymmetric simple exclusion process (BASEP) was recently introduced. The BASEP refers to diffusion models where hard spheres are driven by a constant drag force through a periodic potential. Here we derive general properties of the rich collective dynamics in the BASEP. Average currents in the steady state change dramatically with the particle size and density. For an open system coupled to particle reservoirs, extremal current principles predict various nonequilibrium phases, which we verify by Brownian dynamics simulations. For general pair interactions we discuss connections to single-file transport by traveling-wave potentials and prove the impossibility of current reversals in systems driven by a constant drag and by traveling waves.

Exact large deviation statistics and trajectory phase transition of a deterministic boundary driven cellular automaton

We study the statistical properties of the long-time dynamics of the rule 54 reversible cellular automaton (CA), driven stochastically at its boundaries. This CA can be considered as a discrete-time and deterministic version of the Fredrickson-Andersen kinetically constrained model (KCM). By means of a matrix product ansatz, we compute the exact large deviation cumulant generating functions for a wide range of time-extensive observables of the dynamics, together with their associated rate functions and conditioned long-time distributions over configurations. We show that for all instances of boundary driving the CA dynamics occurs at the point of phase coexistence between competing active and inactive dynamical phases, similar to what happens in more standard KCMs. We also find the exact finite size scaling behavior of these trajectory transitions, and provide the explicit "Doob-transformed" dynamics that optimally realizes rare dynamical events.

Sampling rare events across dynamical phase transitions

Interacting particle systems with many degrees of freedom may undergo phase transitions to sustain atypical fluctuations of dynamical observables such as the current or the activity. In some cases, this leads to symmetry-broken space-time trajectories which enhance the probability of such events due to the emergence of ordered structures. Despite their conceptual and practical importance, these dynamical phase transitions (DPTs) at the trajectory level are difficult to characterize due to the low probability of their occurrence. However, during the last decade, advanced computational techniques have been developed to measure rare events in simulations of many-particle systems that allow the direct observation and characterization of these DPTs. Here we review the application of a particular rare-event simulation technique, based on cloning Monte Carlo methods, to characterize DPTs in paradigmatic stochastic lattice gases. In particular, we describe in detail some tricks and tips of the trade, paying special attention to the measurement of order parameters capturing the physics of the different DPTs, as well as to the finite-size effects (both in the system size and in the number of clones) that affect the measurements. Overall, we provide a consistent picture of the phenomenology associated with DPTs and their measurement.

Occupation-time statistics of a gas of interacting diffusing particles

The time that a diffusing particle spends in a certain region of space is known as the occupation time, or the residence time. Recently, the joint occupation-time statistics of an ensemble of noninteracting particles was addressed using the single-particle statistics. Here we employ the macroscopic fluctuation theory (MFT) to study the occupation-time statistics of many interacting particles. We find that interactions can significantly change the statistics and, in some models, even cause a singularity of the large-deviation function describing these statistics. This singularity can be interpreted as a dynamical phase transition. We also point out to a close relation between the MFT description of the occupation-time statistics of noninteracting particles and the level 2 large deviation formalism which describes the occupation-time statistics of a single particle.

Infinite family of universal profiles for heat current statistics in Fourier's law

Using tools from large deviation theory, we study fluctuations of the heat current in a model of d-dimensional incompressible fluid driven out of equilibrium by a temperature gradient. We find that the most probable temperature fields sustaining atypical values of the global current can be naturally classified in an infinite set of curves, allowing us to exhaustively analyze their topological properties and to define universal profiles onto which all optimal fields collapse. We also compute the statistics of empirical heat current, where we find remarkable logarithmic tails for large current fluctuations orthogonal to the thermal gradient. Finally, we determine explicitly a number of cumulants of the current distribution, finding interesting relations between them.

Tracer diffusion in crowded narrow channels

We summarise different results on the diffusion of a tracer particle in lattice gases of hard-core particles with stochastic dynamics, which are confined to narrow channels-single-files, comb-like structures and quasi-one-dimensional channels with the width equal to several particle diameters. We show that in such geometries a surprisingly rich, sometimes even counter-intuitive, behaviour emerges, which is absent in unbounded systems. This is well-documented for the anomalous diffusion in single-files. Less known is the anomalous dynamics of a tracer particle in crowded branching single-files-comb-like structures, where several kinds of anomalous regimes take place. In narrow channels, which are broader than single-files, one encounters a wealth of anomalous behaviours in the case where the tracer particle is subject to a regular external bias: here, one observes an anomaly in the temporal evolution of the tracer particle velocity, super-diffusive at transient stages, and ultimately a giant diffusive broadening of fluctuations in the position of the tracer particle, as well as spectacular multi-tracer effects of self-clogging of narrow channels. Interactions between a biased tracer particle and a confined crowded environment also produce peculiar patterns in the out-of-equilibrium distribution of the environment particles, very different from the ones appearing in unbounded systems. For moderately dense systems, a surprising effect of a negative differential mobility takes place, such that the velocity of a biased tracer particle can be a non-monotonic function of the force. In some parameter ranges, both the velocity and the diffusion coefficient of a biased tracer particle can be non-monotonic functions of the density. We also survey different results obtained for a tracer particle diffusion in unbounded systems, which will permit a reader to have an exhaustively broad picture of the tracer diffusion in crowded environments.

Brownian Asymmetric Simple Exclusion Process

We study the driven Brownian motion of hard rods in a one-dimensional cosine potential with a large amplitude compared to the thermal energy. In a closed system, we find surprising features of the steady-state current in dependence of the particle density. The form of the current-density relation changes greatly with the particle size and can exhibit both a local maximum and minimum. The changes are caused by an interplay of a barrier reduction, blocking, and exchange symmetry effect. The latter leads to a current equal to that of noninteracting particles for a particle size commensurate with the period length of the cosine potential. For an open system coupled to particle reservoirs, we predict five different phases of nonequilibrium steady states to occur. Our results show that the particle size can be of crucial importance for nonequilibrium phase transitions in driven systems. Possible experiments for demonstrating our findings are pointed out.

Anomalous Scaling of Dynamical Large Deviations

The typical values and fluctuations of time-integrated observables of nonequilibrium processes driven in steady states are known to be characterized by large deviation functions, generalizing the entropy and free energy to nonequilibrium systems. The definition of these functions involves a scaling limit, similar to the thermodynamic limit, in which the integration time τ appears linearly, unless the process considered has long-range correlations, in which case τ is generally replaced by τ^{ξ} with ξ≠1. Here, we show that such an anomalous power-law scaling in time of large deviations can also arise without long-range correlations in Markovian processes as simple as the Langevin equation. We describe the mechanism underlying this scaling using path integrals and discuss its physical consequences for more general processes.

Solvation in Space-time: Pretransition Effects in Trajectory Space

We demonstrate pretransition effects in space-time in trajectories of systems in which the dynamics displays a first-order phase transition between distinct dynamical phases. These effects are analogous to those observed for thermodynamic first-order phase transitions, most notably the hydrophobic effect in water. Considering the (infinite temperature) East model as an elementary example, we study the properties of "space-time solvation" by examining trajectories where finite space-time regions are conditioned to be inactive in an otherwise active phase. We find that solvating an inactive region of space-time within an active trajectory shows two regimes in the dynamical equivalent of solvation free energy: an "entropic" small solute regime in which uncorrelated fluctuations are sufficient to evacuate activity from the solute, and an "energetic" large solute regime which involves the formation of a solute-induced inactive domain with an associated active-inactive interface bearing a dynamical interfacial tension. We also show that as a result of this dynamical interfacial tension there is a dynamical analog of the hydrophobic collapse that drives the assembly of large hydrophobes in water. We discuss the general relevance of these results to the properties of dynamical fluctuations in systems with slow collective relaxation such as glass formers.

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.

Geometrical interpretation of dynamical phase transitions in boundary-driven systems

Dynamical phase transitions are defined as nonanalytic points of the large deviation function of current fluctuations. We show that for boundary-driven systems, many dynamical phase transitions can be identified using the geometrical structure of an effective potential of a Hamiltonian, recovered from the macroscopic fluctuation theory description. Using this method we identify new dynamical phase transitions that could not be recovered using existing perturbative methods. Moreover, using the Hamiltonian picture, an experimental scheme is suggested to demonstrate an analog of dynamical phase transitions in linear, rather than exponential, time.

Fluctuating hydrodynamics, current fluctuations, and hyperuniformity in boundary-driven open quantum chains

We consider a class of either fermionic or bosonic noninteracting open quantum chains driven by dissipative interactions at the boundaries and study the interplay of coherent transport and dissipative processes, such as bulk dephasing and diffusion. Starting from the microscopic formulation, we show that the dynamics on large scales can be described in terms of fluctuating hydrodynamics. This is an important simplification as it allows us to apply the methods of macroscopic fluctuation theory to compute the large deviation (LD) statistics of time-integrated currents. In particular, this permits us to show that fermionic open chains display a third-order dynamical phase transition in LD functions. We show that this transition is manifested in a singular change in the structure of trajectories: while typical trajectories are diffusive, rare trajectories associated with atypical currents are ballistic and hyperuniform in their spatial structure. We confirm these results by numerically simulating ensembles of rare trajectories via the cloning method, and by exact numerical diagonalization of the microscopic quantum generator.

Order and Symmetry Breaking in the Fluctuations of Driven Systems

Dynamical phase transitions (DPTs) in the space of trajectories are one of the most intriguing phenomena of nonequilibrium physics, but their nature in realistic high-dimensional systems remains puzzling. Here we observe for the first time a DPT in the current vector statistics of an archetypal two-dimensional (2D) driven diffusive system and characterize its properties using the macroscopic fluctuation theory. The complex interplay among the external field, anisotropy, and vector currents in 2D leads to a rich phase diagram, with different symmetry-broken fluctuation phases separated by lines of first- and second-order DPTs. Remarkably, different types of 1D order in the form of jammed density waves emerge to hinder transport for low-current fluctuations, revealing a connection between rare events and self-organized structures which enhance their probability.

Variational calculation of transport coefficients in diffusive lattice gases

A diffusive lattice gas is characterized by the diffusion coefficient depending only on the density. The Green-Kubo formula for diffusivity can be represented as a variational formula, but even when the equilibrium properties of a lattice gas are analytically known, the diffusion coefficient can be computed only in the exceptional situation when the lattice gas is gradient. In the general case, minimization over an infinite-dimensional space is required. We propose an approximation scheme based on minimizing over finite-dimensional subspaces of functions. The procedure is demonstrated for one-dimensional generalized exclusion processes in which each site can accommodate at most two particles. Our analytical predictions provide upper bounds for the diffusivity that are very close to simulation results throughout the entire density range. We also analyze nonequilibrium density profiles for finite chains coupled to reservoirs. The predictions for the profiles are in excellent agreement with simulations.

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.

Numerical study of continuous and discontinuous dynamical phase transitions for boundary-driven systems

The existence and search for thermodynamic phase transitions is of unfading interest. In this paper, we present numerical evidence of dynamical phase transitions occurring in boundary-driven systems with a constrained integrated current. It is shown that certain models exhibit a discontinuous transition between two different density profiles and a continuous transition between a time-independent and a time-dependent profile. We also verified that the Kipnis-Marchioro-Presutti model exhibits no phase transitions in a range much larger than previously explored.

Simple bounds on fluctuations and uncertainty relations for first-passage times of counting observables

Recent large deviation results have provided general lower bounds for the fluctuations of time-integrated currents in the steady state of stochastic systems. A corollary are so-called thermodynamic uncertainty relations connecting precision of estimation to average dissipation. Here we consider this problem but for counting observables, i.e., trajectory observables which, in contrast to currents, are non-negative and nondecreasing in time (and possibly symmetric under time reversal). In the steady state, their fluctuations to all orders are bound from below by a Conway-Maxwell-Poisson distribution dependent only on the averages of the observable and of the dynamical activity. We show how to obtain the corresponding bounds for first-passage times (times when a certain value of the counting variable is first reached) and their uncertainty relations. Just like entropy production does for currents, dynamical activity controls the bounds on fluctuations of counting observables.