Breakdown of hydrodynamics in a simple one-dimensional fluid

Blast in a One-Dimensional Cold Gas: From Newtonian Dynamics to Hydrodynamics

A gas composed of a large number of atoms evolving according to Newtonian dynamics is often described by continuum hydrodynamics. Proving this rigorously is an outstanding open problem, and precise numerical demonstrations of the equivalence of the hydrodynamic and microscopic descriptions are rare. We test this equivalence in the context of the evolution of a blast wave, a problem that is expected to be at the limit where hydrodynamics could work. We study a one-dimensional gas at rest with instantaneous localized release of energy for which the hydrodynamic Euler equations admit a self-similar scaling solution. Our microscopic model consists of hard point particles with alternating masses, which is a nonintegrable system with strong mixing dynamics. Our extensive microscopic simulations find a remarkable agreement with Euler hydrodynamics, with deviations in a small core region that are understood as arising due to heat conduction.

Coexistence of Ballistic and Fourier Regimes in the β Fermi-Pasta-Ulam-Tsingou Lattice

Commonly, thermal transport properties of one-dimensional systems are found to be anomalous. Here, we perform a numerical and theoretical study of the β-Fermi-Pasta-Ulam-Tsingou chain, considered a prototypical model for one-dimensional anharmonic crystals, in contact with thermostats at different temperatures. We give evidence that, in steady state conditions, the local wave energy spectrum can be naturally split into modes that are essentially ballistic (noninteracting or scarcely interacting) and kinetic modes (interacting enough to relax to local thermodynamic equilibrium). We show numerically that the well-known divergence of the energy conductivity is related to how the transition region between these two sets of modes shifts in k space with the system size L, due to properties of the collision integral of the system. Moreover, we show that the kinetic modes are responsible for a macroscopic behavior compatible with Fourier's law. Our work sheds light on the long-standing problem of the applicability of standard thermodynamics in one-dimensional nonlinear chains, testbed for understanding the thermal properties of nanotubes and nanowires.

Heat flux in one-dimensional systems

Understanding heat transport in one-dimensional systems remains a major challenge in theoretical physics, both from the quantum as well as from the classical point of view. In fact, steady states of one-dimensional systems are commonly characterized by macroscopic inhomogeneities, and by long-range correlations, as well as large fluctuations that are typically absent in standard three-dimensional thermodynamic systems. These effects violate locality-material properties in the bulk may be strongly affected by the boundaries, leading to anomalous energy transport-and they make more problematic the interpretation of mechanical microscopic quantities in terms of thermodynamic observables. Here, we revisit the problem of heat conduction in chains of classical nonlinear oscillators, following a Lagrangian and a Eulerian approach. The Eulerian definition of the flux is composed of a convective and a conductive component. The former component tends to prevail at large temperatures where the system behavior is increasingly gaslike. Finally, we find that the convective component tends to be negative in the presence of a negative pressure.

A violation of universality in anomalous Fourier's law

Since the discovery of long-time tails, it has been clear that Fourier's law in low dimensions is typically anomalous, with a size-dependent heat conductivity, though the nature of the anomaly remains puzzling. The conventional wisdom, supported by renormalization-group arguments and mode-coupling approximations within fluctuating hydrodynamics, is that the anomaly is universal in 1d momentum-conserving systems and belongs in the Lévy/Kardar-Parisi-Zhang universality class. Here we challenge this picture by using a novel scaling method to show unambiguously that universality breaks down in the paradigmatic 1d diatomic hard-point fluid. Hydrodynamic profiles for a broad set of gradients, densities and sizes all collapse onto an universal master curve, showing that (anomalous) Fourier's law holds even deep into the nonlinear regime. This allows to solve the macroscopic transport problem for this model, a solution which compares flawlessly with data and, interestingly, implies the existence of a bound on the heat current in terms of pressure. These results question the renormalization-group and mode-coupling universality predictions for anomalous Fourier's law in 1d, offering a new perspective on transport in low dimensions.

Probing local equilibrium in nonequilibrium fluids

We use extensive computer simulations to probe local thermodynamic equilibrium (LTE) in a quintessential model fluid, the two-dimensional hard-disks system. We show that macroscopic LTE is a property much stronger than previously anticipated, even in the presence of important finite-size effects, revealing a remarkable bulk-boundary decoupling phenomenon in fluids out of equilibrium. This allows us to measure the fluid's equation of state in simulations far from equilibrium, with an excellent accuracy comparable to the best equilibrium simulations. Subtle corrections to LTE are found in the fluctuations of the total energy which strongly point to the nonlocality of the nonequilibrium potential governing the fluid's macroscopic behavior out of equilibrium.

Scaling laws and bulk-boundary decoupling in heat flow

When driven out of equilibrium by a temperature gradient, fluids respond by developing a nontrivial, inhomogeneous structure according to the governing macroscopic laws. Here we show that such structure obeys strikingly simple scaling laws arbitrarily far from equilibrium, provided that both macroscopic local equilibrium and Fourier's law hold. Extensive simulations of hard disk fluids confirm the scaling laws even under strong temperature gradients, implying that Fourier's law remains valid in this highly nonlinear regime, with putative corrections absorbed into a nonlinear conductivity functional. In addition, our results show that the scaling laws are robust in the presence of strong finite-size effects, hinting at a subtle bulk-boundary decoupling mechanism which enforces the macroscopic laws on the bulk of the finite-sized fluid. This allows one to measure the marginal anomaly of the heat conductivity predicted for hard disks.

Symmetries in fluctuations far from equilibrium

Fluctuations arise universally in nature as a reflection of the discrete microscopic world at the macroscopic level. Despite their apparent noisy origin, fluctuations encode fundamental aspects of the physics of the system at hand, crucial to understand irreversibility and nonequilibrium behavior. To sustain a given fluctuation, a system traverses a precise optimal path in phase space. Here we show that by demanding invariance of optimal paths under symmetry transformations, new and general fluctuation relations valid arbitrarily far from equilibrium are unveiled. This opens an unexplored route toward a deeper understanding of nonequilibrium physics by bringing symmetry principles to the realm of fluctuations. We illustrate this concept studying symmetries of the current distribution out of equilibrium. In particular we derive an isometric fluctuation relation that links in a strikingly simple manner the probabilities of any pair of isometric current fluctuations. This relation, which results from the time-reversibility of the dynamics, includes as a particular instance the Gallavotti-Cohen fluctuation theorem in this context but adds a completely new perspective on the high level of symmetry imposed by time-reversibility on the statistics of nonequilibrium fluctuations. The new symmetry implies remarkable hierarchies of equations for the current cumulants and the nonlinear response coefficients, going far beyond Onsager's reciprocity relations and Green-Kubo formulas. We confirm the validity of the new symmetry relation in extensive numerical simulations, and suggest that the idea of symmetry in fluctuations as invariance of optimal paths has far-reaching consequences in diverse fields.

Large fluctuations of the macroscopic current in diffusive systems: a numerical test of the additivity principle

Most systems, when pushed out of equilibrium, respond by building up currents of locally conserved observables. Understanding how microscopic dynamics determines the averages and fluctuations of these currents is one of the main open problems in nonequilibrium statistical physics. The additivity principle is a theoretical proposal that allows to compute the current distribution in many one-dimensional nonequilibrium systems. Using simulations, we validate this conjecture in a simple and general model of energy transport, both in the presence of a temperature gradient and in canonical equilibrium. In particular, we show that the current distribution displays a Gaussian regime for small current fluctuations, as prescribed by the central limit theorem, and non-Gaussian (exponential) tails for large current deviations, obeying in all cases the Gallavotti-Cohen fluctuation theorem. In order to facilitate a given current fluctuation, the system adopts a well-defined temperature profile different from that of the steady state and in accordance with the additivity hypothesis predictions. System statistics during a large current fluctuation is independent of the sign of the current, which implies that the optimal profile (as well as higher-order profiles and spatial correlations) are invariant upon current inversion. We also demonstrate that finite-time joint fluctuations of the current and the profile are well described by the additivity functional. These results suggest the additivity hypothesis as a general and powerful tool to compute current distributions in many nonequilibrium systems.

Test of the additivity principle for current fluctuations in a model of heat conduction

The additivity principle allows to compute the current distribution in many one-dimensional (1D) nonequilibrium systems. Using simulations, we confirm this conjecture in the 1D Kipnis-Marchioro-Presutti model of heat conduction for a wide current interval. The current distribution shows both Gaussian and non-Gaussian regimes, and obeys the Gallavotti-Cohen fluctuation theorem. We verify the existence of a well-defined temperature profile associated to a given current fluctuation. This profile is independent of the sign of the current, and this symmetry extends to higher-order profiles and spatial correlations. We also show that finite-time joint fluctuations of the current and the profile are described by the additivity functional. These results suggest the additivity hypothesis as a general and powerful tool to compute current distributions in many nonequilibrium systems.

Dry and wet granular shock waves

The formation of a shock wave in one-dimensional granular gases is considered, for both the dry and the wet cases, and the results are compared with the analytical shock wave solution in a sticky gas. Numerical simulations show that the behavior of the shock wave in both cases tends asymptotically to the sticky limit. In the inelastic gas (dry case) there is a very close correspondence to the sticky gas, with one big cluster growing in the center of the shock wave, and a step-like stationary velocity profile. In the wet case, the shock wave has a nonzero width which is marked by two symmetric heavy clusters performing breathing oscillations with slowly increasing amplitude. All three models have the same asymptotic energy dissipation law, which is important in the context of the free cooling scenario. For the early stage of the shock formation and asymptotic oscillations we provide analytical results as well.

Simplest piston problem. I. Elastic collisions

We study the dynamics of three elastic particles in a finite interval where two light particles are separated by a heavy "piston." The piston undergoes surprisingly complex motion that is oscillatory at short time scales but seemingly chaotic at longer scales. The piston also makes long-duration excursions close to the ends of the interval that stem from the breakdown of energy equipartition. Many of these dynamical features can be understood by mapping the motion of three particles on the line onto the trajectory of an elastic billiard in a highly skewed tetrahedral region. We exploit this picture to construct a qualitative random walk argument that predicts a power-law tail, with exponent -3/2, for the distribution of time intervals between successive piston crossings of the interval midpoint. These predictions are verified by numerical simulations.