Emergent Hydrodynamics in an Exclusion Process with Long-Range Interactions

This paper investigates the symmetric Dyson exclusion process, a lattice gas with long-range Coulomb-type interactions, by mapping its stochastic generator to a free-fermion quantum chain to derive and validate a conjectured non-local hydrodynamic equation that accurately describes the system's ballistic scaling and limit shapes.

The Gist

Imagine a crowded dance floor where everyone is trying to move around, but they have very strict rules. This is the world of stochastic particle systems, a field of physics that studies how groups of tiny particles move and interact.

Usually, when we look at a crowd of particles (like gas molecules in a room), we expect them to behave like a fluid. If you push them from one side, they flow smoothly to the other. In most cases, the movement of a particle depends only on its immediate neighbors—like a person in a line only caring about the person directly in front of them. This is called local interaction.

However, the paper you shared introduces a very special, unusual dance floor called the Symmetric Dyson Exclusion Process (SDEP). Here, the rules are twisted in a fascinating way.

The "Telepathic" Dance Floor

In this specific model, particles don't just care about their neighbors. They are "telepathic." Every particle feels a gentle, long-range pull or push from every other particle on the entire floor, no matter how far away they are.

Think of it like this:

• Normal Crowd: If you are in a line at a coffee shop, you only care about the person directly in front of you. If they move, you move.
• The SDEP Crowd: If you are in this line, you can "feel" the presence of the person at the very back of the line. If they move, it subtly changes your desire to move, even though you are miles apart.

This "long-range" connection is mathematically similar to how electric charges repel each other (Coulomb interactions) or how eigenvalues of random matrices behave (Dyson's Brownian motion).

The Big Discovery: A "Ghost" Current

The authors wanted to know: If we zoom out and look at this crowd from a distance (like a satellite view), how does the density of people change over time?

Usually, physics tells us that the flow of people (current) depends only on how crowded the spot is right now.

• Normal Physics: "If it's crowded here, people flow out." (Local)
• This Paper's Discovery: "If it's crowded here, people flow out based on the average crowding of the entire room." (Non-local)

They found that the flow of particles is governed by a mathematical tool called the Hilbert Transform.

• The Analogy: Imagine you are trying to predict the traffic flow on a highway. In a normal highway, you look at the cars right next to you. In this "SDEP highway," the speed of the cars at mile marker 10 depends on a "ghostly average" of the traffic density at mile markers 1, 50, and 100 simultaneously. The current is a non-local functional—it's a recipe that requires knowing the state of the whole system to determine what happens at one spot.

The "Melting Iceberg" and the Arctic Curve

To prove this, the authors simulated what happens when you start with a solid block of particles (like a frozen iceberg of people) and let them spread out.

In normal systems, the edge of the block melts away smoothly and diffusively, like butter melting on toast.
In this SDEP system, something magical happens:

1. Frozen Zones: The center of the block stays perfectly packed (density = 1), and the empty space outside stays perfectly empty (density = 0). These areas are "frozen."
2. The Fluctuating Zone: Between the frozen center and the empty outside, there is a chaotic, melting region where the density is somewhere between 0 and 1.
3. The Arctic Curve: The boundary between the "frozen" and "fluctuating" zones forms a sharp, distinct curve. The authors call this an Arctic Curve.

The Metaphor: Imagine an iceberg melting in the ocean. Usually, the edges are fuzzy and irregular. But in this model, the iceberg melts in a way that creates a perfect, sharp geometric shape separating the solid ice from the water. The "Arctic Curve" is the sharp line where the ice stops being solid and starts becoming liquid chaos.

Why Does This Matter?

This paper is a breakthrough because it breaks a long-held rule in physics. For decades, scientists believed that if you zoomed out far enough, all particle systems would behave like simple, local fluids (where only neighbors matter).

This paper shows that long-range interactions can create a new kind of fluid dynamics where the whole system talks to itself instantly.

• Real-world connection: While this is a theoretical model, it helps us understand complex systems where long-range forces matter, such as:
  • Quantum systems: How electrons behave in certain materials.
  • Random Matrices: Used in everything from nuclear physics to the stock market.
  • Biological systems: How proteins or cells might coordinate over long distances.

Summary

The authors took a complex mathematical model of particles that "feel" each other across long distances. They discovered that when these particles move, they don't just follow local rules. Instead, they create a global, non-local flow that can be described by a beautiful, sharp boundary (the Arctic Curve) separating order from chaos. It's like finding a new law of physics where the crowd knows what everyone else is doing, not just who is standing next to them.

Technical Summary

1. Problem Statement

The paper addresses a fundamental gap in the theory of stochastic interacting particle systems: the emergence of macroscopic hydrodynamics in systems with long-range interactions.

• Context: In standard short-range systems (e.g., the Symmetric Simple Exclusion Process, SSEP), macroscopic behavior is well-understood. The coarse-grained density $\rho(x,t)$ evolves according to a local continuity equation $\partial_t \rho + \partial_x j(\rho) = 0$, where the current $j$ is a local function of the density.
• The Challenge: For systems with long-range interactions, it is unclear whether the current remains a local functional of the density or becomes non-local. Exact results for such systems are scarce.
• The Model: The authors focus on the Symmetric Dyson Exclusion Process (SDEP). This is a lattice gas with exclusion (at most one particle per site) and long-range, Coulomb-type interactions.
  • It arises as the maximal-activity limit of the SSEP (conditioned on a large number of jumps).
  • It is a discrete version of Dyson's Brownian motion on the unitary group.
  • It maps to the ground state of the spin-1/2 XX quantum chain.

2. Methodology

The authors employ a multi-faceted approach combining exact microscopic mappings, heuristic derivations, and large-scale numerical simulations.

• Exact Mapping (Doob Transform):
  The SDEP is defined by an intensity matrix (stochastic generator) $H_N$. The authors utilize an exact Doob transform (ground-state transformation) to map this stochastic generator onto the Hamiltonian of the spin-1/2 XX quantum chain ($H_{XX}$).
  $$H_N = \hat{\Psi}_N (H_{XX}^N - E_N) \hat{\Psi}_N^{-1}$$
  Since the XX chain is a system of free fermions, this mapping allows the authors to leverage exact results from quantum many-body theory, specifically the free-fermion representation.

• Heuristic Derivation of Hydrodynamics:
  Starting from the exact microscopic expression for the instantaneous current, the authors perform a coarse-graining analysis. They separate the interaction into short-range (local) and long-range (non-local) components.

  • The long-range part sums to a Hilbert transform of the density profile, denoted as $H\rho$.
  • The short-range part is assumed to depend only on the local density $\rho$ under local equilibrium.
  • This leads to an ansatz for the current involving hyperbolic functions of the Hilbert transform.

• Complex Hydrodynamic Formulation:
  The authors demonstrate that the derived non-local single-field equation is mathematically equivalent to a local two-field system. By introducing a complex field $P(z,t) = \pi\rho + i\pi H\rho$, the evolution equation transforms into a complex inviscid Burgers (or complex Hopf) equation:
  $$\partial_t P - i \sin(P) \partial_x P = 0$$
  This equivalence allows for the use of analytic function theory (specifically, the properties of analytic functions in the upper half-plane) to solve the hydrodynamic equations.

• Numerical Verification:
  The theoretical predictions are tested against Kinetic Monte Carlo (KMC) simulations of the microscopic SDEP dynamics. The authors simulate the "melting" of initial block states (single and double blocks of particles) and compare the resulting density profiles and space-time trajectories with the analytical solutions of the hydrodynamic equations.

3. Key Contributions and Results

A. Emergence of Non-Local Hydrodynamics

The primary theoretical result is the derivation of the exact macroscopic hydrodynamic equation for the SDEP. Unlike standard diffusive systems, the SDEP exhibits ballistic (Eulerian) scaling ($z=1$) governed by a non-local current:
$$\partial_t \rho(x,t) + \partial_x j[\rho](x,t) = 0$$
$$j[\rho](x,t) = \frac{1}{\pi} \sin(\pi \rho(x,t)) \sinh(\pi H\rho(x,t))$$
where $H$ is the Hilbert transform. This proves that long-range interactions can destroy the locality of the current, making it a functional of the global density profile.

B. Equivalence to Local Two-Field System

The paper establishes that the non-local single-field description is equivalent to a local system of two coupled conservation laws for densities $\rho$ and $\tilde{\rho}$ (where $\tilde{\rho} = H\rho$):
$$\begin{cases} \partial_t (\pi\rho) + \partial_x [\frac{1}{\pi} \sin(\pi\rho) \sinh(\pi\tilde{\rho})] = 0 \\ \partial_t (\pi\tilde{\rho}) + \partial_x [\frac{1}{\pi} \cos(\pi\rho) \cosh(\pi\tilde{\rho})] = 0 \end{cases}$$
This equivalence simplifies the analysis and connects the problem to the complex Hopf equation, a known integrable system.

C. Limit Shapes and Arctic Curves

The authors solve the hydrodynamic equations for specific initial conditions (single and double blocks of particles).

• Limit Shape Phenomenon: They observe the emergence of sharp boundaries separating "frozen" regions (where density is strictly 0 or 1) from "fluctuating" regions (where $0 < \rho < 1$).
• Arctic Curves: The boundary between these regions is the "arctic curve." The authors derive explicit formulas for these curves using the discriminant of the polynomial equations arising from the method of characteristics.
• Agreement: The analytically derived arctic curves and density profiles show excellent agreement with Monte Carlo simulations, validating the hydrodynamic conjecture.

D. Connection to Continuous Dyson Gas

In the low-density limit ($\rho \ll 1$), the derived equation reduces to the known hydrodynamic equation for the continuous Dyson gas ($\partial_t \rho + \partial_x [\pi \rho H\rho] = 0$). However, the lattice model exhibits unique features (like the semi-circle expansion at large times) due to the hard-core exclusion constraint, which are absent in the continuous case.

4. Significance

• Paradigm for Non-Local Hydrodynamics: The SDEP serves as a rare, exactly solvable example of a system where long-range interactions lead to genuinely non-local hydrodynamic equations. It challenges the assumption that macroscopic currents are always local functions of density.
• Bridge Between Fields: The work creates a strong link between stochastic particle systems, random matrix theory (Dyson Brownian motion), and quantum integrable systems (XX chain, free fermions).
• Analytical Tools: The demonstration that non-local equations can be mapped to local complex PDEs (complex Hopf) provides a powerful new tool for analyzing other long-range interacting systems.
• Future Directions: The paper opens avenues for studying non-local Kardar-Parisi-Zhang (KPZ) regimes, macroscopic fluctuation theory with hard-core constraints, and boundary-driven phase transitions in long-range systems.

In summary, the paper successfully derives and validates a closed, non-local hydrodynamic equation for a long-range interacting lattice gas, revealing rich emergent phenomena like arctic curves and establishing a deep connection between stochastic exclusion processes and integrable quantum chains.