From Near-Integrable to Far-from-Integrable: A Unified Picture of Thermalization and Heat Transport

This paper presents a unified theoretical framework for thermalization and heat transport in a one-dimensional diatomic hard-point gas by mapping a phase diagram that reveals three universal dynamical regimes—kinetic, hydrodynamic, and intermediate—across the spectrum from near-integrable to far-from-integrable systems, thereby challenging conventional views on system-size dependence and establishing a consistent description of relaxation dynamics.

The Gist

The Big Question: How Does Chaos Take Over Order?

Imagine you have a room full of people (particles) who are all standing perfectly still in a line. Suddenly, you tell them to start walking.

• Scenario A: Everyone walks in perfect lockstep, never bumping into each other, just swapping places smoothly. The system stays "ordered" forever.
• Scenario B: Everyone starts running wildly, bumping into each other, shouting, and eventually, the whole room becomes a chaotic, noisy mess where everyone is moving randomly. This is thermalization (reaching a state of equilibrium).

For decades, physicists have known how systems reach this chaotic state when they are slightly disturbed (like a room where people bump into each other a little bit). But they didn't have a complete map for what happens when the disturbance is huge, or when the system is somewhere in the middle.

This paper provides that complete map. It studies a specific model (a 1D gas of bouncing balls with different weights) to show exactly how a system moves from "perfect order" to "total chaos" and how heat moves through it.

---

The Cast of Characters: The "Hard-Point Gas"

The scientists used a model called a Diatomic Hard-Point (DHP) Gas.

• The Setup: Imagine a 1D hallway with balls bouncing back and forth.
• The Twist: The balls have two different weights (light and heavy).
• The Control Knob ($\delta$): This is the "non-integrability strength."
  • If the weights are identical ($\delta = 0$), the balls just pass through each other like ghosts. Nothing changes; the system never gets chaotic.
  • If the weights are different ($\delta > 0$), they collide and bounce. The bigger the difference in weight, the more chaotic the collisions become.

The Three Zones of Chaos (The Phase Diagram)

The paper reveals that depending on how different the weights are and how many balls you have, the system behaves in three distinct ways. Think of this like driving a car on different types of roads:

1. The "Highway" Zone (Near-Integrable)

• When: The weight difference is very small.
• What happens: The collisions are gentle. The balls mostly keep their original energy patterns.
• The Analogy: Imagine a highway where cars are driving fast but rarely change lanes. If one car speeds up, it takes a long time for that speed change to spread to the rest of the traffic.
• The Physics: The system relaxes (calms down) very slowly. The time it takes to reach equilibrium depends heavily on how small the weight difference is (specifically, it scales with $1/\delta^2$). It's like waiting for a whisper to travel across a huge stadium.

2. The "Traffic Jam" Zone (Far-from-Integrable)

• When: The weight difference is huge (e.g., a ping-pong ball hitting a bowling ball).
• What happens: The collisions are violent and frequent. The system forgets its initial state almost instantly.
• The Analogy: Imagine a mosh pit. As soon as one person moves, everyone else is jostled immediately. The "memory" of where you started is lost in seconds.
• The Physics: The system relaxes very fast, but the pattern of relaxation changes. Instead of a smooth decay, it follows a "power law" (a specific mathematical curve). Crucially, even a small group of people in this mosh pit will behave chaotically. You don't need a stadium to see the chaos; the chaos is so intense it happens everywhere.

3. The "Transition" Zone (The Bogoliubov Phase)

• When: The weight difference is in the middle.
• What happens: This is the messy middle ground. The system starts behaving like the "Highway" (slow, orderly) but then suddenly switches to the "Traffic Jam" (fast, chaotic).
• The Analogy: It's like a calm river that suddenly hits a waterfall. You see the smooth flow first, then the turbulence.
• The Physics: This is where the famous "three-stage" theory by physicist Bogoliubov happens: First, the balls fly freely (ballistic); then they collide locally (kinetic); finally, they move as a collective wave (hydrodynamic).

The Great Surprise: Size Doesn't Always Matter

For a long time, physicists thought that "hydrodynamic" effects (where the whole system moves like a fluid wave) only happened in huge systems (like the atmosphere or a large ocean).

This paper flips that idea on its head.
They found that if the "chaos" (the weight difference) is strong enough, even a tiny system (just a few dozen particles) will act like a fluid wave. You don't need a massive crowd to create a mosh pit; you just need enough energy and difference to make the collisions intense.

The "Order of Operations" Paradox

The paper also solves a philosophical puzzle about limits. Imagine you are trying to predict the future of this gas. You have two choices:

1. Make the system infinitely big first, then make the weights identical.
2. Make the weights identical first, then make the system infinitely big.

The result is totally different!

• If you make the system big first, it acts like a calm highway (Kinetic dominance).
• If you make the weights identical first, it acts like a fluid wave (Hydrodynamic dominance).

This is like saying: "If you build a bigger city before you fix the traffic lights, you get gridlock. But if you fix the traffic lights before you build the city, you get smooth flow." The order in which you apply the changes matters immensely.

Heat Transport: The Same Story Twice

The paper also looked at Heat Transport (how heat moves through the gas).

• The Discovery: The way heat moves is exactly the same as the way the system relaxes to equilibrium.
• The Analogy: If the system is in the "Highway" zone, heat moves normally (like a steady stream). If it's in the "Traffic Jam" zone, heat moves strangely (anomalous), getting stuck or moving too fast depending on the size of the system.
• Why it matters: This proves that "Thermalization" (getting hot/equilibrium) and "Heat Transport" (moving heat) are two sides of the same coin. You can't understand one without the other.

The Bottom Line

This paper draws a complete map of how order turns into chaos.

1. Small disturbances = Slow, orderly relaxation (Kinetic).
2. Big disturbances = Fast, fluid-like chaos (Hydrodynamic), even in small systems.
3. The Middle = A complex transition where both happen.

This unified picture helps scientists understand everything from how heat moves in tiny computer chips to how energy behaves in quantum computers. It tells us that the "size" of the system isn't the only thing that matters; the "strength" of the chaos is just as important.

Technical Summary

1. Problem Statement

The paper addresses a fundamental gap in nonequilibrium statistical physics: the lack of a unified theoretical framework describing relaxation dynamics and heat transport across the entire spectrum of non-integrability.

• Current State: Previous research has largely focused on near-integrable systems (e.g., dilute gases, weakly perturbed lattices), where thermalization is dominated by kinetic processes and follows a universal scaling law ($\tau \propto \delta^{-2}$, where $\delta$ is the non-integrability strength). Conversely, far-from-integrable regimes are often treated separately, where hydrodynamic effects dominate.
• The Gap: There is no comprehensive understanding of the transition between these regimes, nor a unified description linking thermalization (local energy relaxation in non-equilibrium states) with heat transport (decay of heat current fluctuations in equilibrium). Additionally, the interplay between system size ($N$) and non-integrability strength ($\delta$) in the thermodynamic limit remains ambiguous.

2. Methodology

The authors employ a combination of analytical theory and high-precision numerical simulations using the one-dimensional diatomic hard-point (DHP) gas model.

• Model: A 1D system of $N$ particles with alternating masses $m_1 = 1 - \delta/2$ and $m_2 = 1 + \delta/2$.
  • $\delta = 0$: Integrable limit (elastic collisions simply exchange velocities).
  • $\delta \in (0, 2)$: Non-integrable regime. The parameter $\delta$ serves as the sole intrinsic control for non-integrability strength.
• Analytical Approach:
  • Kinetic Theory: Utilizes the molecular chaos hypothesis (Stosszahlansatz) to derive evolution equations for mean energy and heat current.
  • Transition Matrix Analysis: Analyzes the eigenvalues of the collision matrix to determine the decay rates of energy fluctuations.
  • Timescale Hierarchy: Identifies three characteristic timescales:
    1. $\tau_f$: Mean free time (ballistic regime).
    2. $\tau_k$: Kinetic relaxation time (local collisions).
    3. $\tau_h \sim N/c_s$: Hydrodynamic relaxation time (collective modes).
• Numerical Simulations:
  • Event-Driven Molecular Dynamics: Simulates particle collisions to track the evolution of local energy and heat current.
  • Metrics:
    • Inverse Participation Ratio (IPR), $\Phi(t)$: Used to measure local energy relaxation (thermalization).
    • Heat Current Autocorrelation Function (HCAF), $C(t)$: Used to measure heat transport properties via the Green-Kubo formula.

3. Key Contributions

The paper constructs a comprehensive phase diagram that unifies thermalization and heat transport across all regimes of non-integrability. The primary contributions are:

1. Identification of Three Universal Dynamical Regimes:

   • Near-Integrable Regime (Kinetic Dominance): Local collisions dominate. Relaxation is exponential.
   • Far-from-Integrable Regime (Hydrodynamic Dominance): Collective modes dominate. Relaxation follows power-law decay.
   • Intermediate Regime (Bogoliubov Phase): A crossover region where ballistic, kinetic, and hydrodynamic effects coexist, characterized by a transition from exponential to power-law decay.

2. Unified Description of Thermalization and Transport:

   • Demonstrates that the relaxation mechanisms for local energy (non-equilibrium) and heat current (equilibrium) are governed by the same underlying dynamics and phase boundaries.
   • Resolves long-standing ambiguities in the DHP model regarding the divergence exponents of heat conductivity.

3. Resolution of the Thermodynamic Limit Paradox:

   • Clarifies that the relaxation behavior in the thermodynamic limit depends on the order of limits taken:
     • $N \to \infty$ first, then $\delta \to 0$: Kinetic dominance (normal transport).
     • $\delta \to 0$ first, then $N \to \infty$: Hydrodynamic dominance (anomalous transport).

4. Key Results

A. Thermalization Dynamics (Local Energy Relaxation)

• Scaling Laws:
  • Near-Integrable ($\delta \to 0$): The thermalization time $\tau$ scales as $\tau \propto \delta^{-2}$. The IPR decays exponentially. This is independent of system size $N$ for small $N$.
  • Far-from-Integrable ($\delta \to 2$): The thermalization time scales linearly with system size, $\tau \propto N$. The IPR decays as a power law ($\Phi(t) \sim t^{-2/3}$).
  • Critical Point: At $\delta = \sqrt{2}$, the system achieves equipartition in a single collision ($\lambda=0$), but strong velocity correlations ($g=1$) introduce significant hydrodynamic memory effects.
• Finite-Size Effects:
  • Hydrodynamic effects (power-law decay) can emerge in small systems if the system is sufficiently far from the integrable regime ($\delta$ is large). This challenges the conventional view that hydrodynamics only appear in the large-$N$ limit.
  • A critical system size $N_c \propto \delta^{-2}$ exists; below $N_c$, kinetics dominate; above $N_c$, hydrodynamics dominate.

B. Heat Transport (Equilibrium Fluctuations)

• HCAF Behavior: The Heat Current Autocorrelation Function mirrors the thermalization results.
  • Kinetic Regime: Exponential decay of HCAF, leading to normal heat conduction (Fourier's law, $\kappa$ is size-independent).
  • Hydrodynamic Regime: Power-law decay of HCAF ($C(t) \sim t^{-2/3}$), leading to anomalous heat conduction ($\kappa \propto N^{1/3}$).
• Unified Framework: The transition from normal to anomalous conduction is directly linked to the transition from kinetic to hydrodynamic dominance in the phase diagram.

C. The Phase Diagram

The authors map the relaxation behavior as a function of $\delta$ (non-integrability) and $N$ (system size).

• Region I (Left): Kinetic dominance. Exponential decay. Normal transport.
• Region II (Middle): Bogoliubov phase. Crossover from exponential to power-law.
• Region III (Right): Hydrodynamic dominance. Power-law decay. Anomalous transport.

5. Significance

• Theoretical Unification: This work provides the first quantitative framework linking integrability breaking, thermalization, and heat transport across the full parameter space. It bridges the gap between kinetic theory (Boltzmann) and hydrodynamic theory (Navier-Stokes/Fluctuating Hydrodynamics).
• Clarification of Scaling Laws: It explains why different studies on the DHP model have reported conflicting results (e.g., varying divergence exponents for heat conductivity); these differences arise from probing different regions of the phase diagram (different $\delta$ and $N$ combinations).
• Implications for Quantum Systems: The authors suggest that these classical findings offer critical insights for quantum thermalization, particularly in understanding prethermalization and the breakdown of integrability in quantum many-body systems.
• Material Design: The identification of a "kinetic transport" regime that is size-independent offers potential pathways for designing high-efficiency thermoelectric materials where anomalous heat conduction is undesirable.

In conclusion, the paper establishes that the relaxation dynamics of many-body systems are not monolithic but are determined by a competition between local collision kinetics and collective hydrodynamic modes, governed by the system's distance from integrability and its size.