🔢 mathematics

Nonisospectral Integrability and Exact Current Fluctuations in the Two-Dimensional SSEP

This paper establishes that annealed current fluctuations in the two-dimensional symmetric simple exclusion process with a circular boundary can be reduced to a radially symmetric nonisospectral problem, enabling the derivation of a closed-form expression for the scaled cumulant generating function via scattering construction and scalar factorization.

Original authors: Tingfei Li

Published 2026-08-11

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Tingfei Li

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a crowded dance floor where everyone is trying to move around without bumping into each other. This is the world of "statistical physics," a branch of science that tries to predict how huge groups of tiny particles behave. Usually, we know the average outcome: if you wait long enough, the crowd spreads out evenly. But what about the rare, wild moments? What if, by pure chance, a huge chunk of the crowd suddenly rushes in one direction, leaving a gap behind? These rare events are called "large deviations." Scientists care about them because they reveal the hidden rules of how energy and matter move when things aren't in balance, like traffic jams or heat flowing through a wall. To study these, researchers use a model called the "Symmetric Simple Exclusion Process" (SSEP). Think of it as a game where particles hop randomly on a grid, but with one strict rule: no two particles can ever occupy the same spot at the same time. It's a simple game, but when you play it with millions of particles, the math gets incredibly tricky, especially when you try to count how many particles cross a specific line over time.

Now, picture a giant, invisible circular fence drawn on this dance floor. Inside the circle, the floor is packed with dancers (high density); outside, it's nearly empty (low density). The dancers can hop over the fence in any direction, but we are watching to see how many more hop out than in over a set period. This is the puzzle tackled in a new paper by Tingfei Li. The author asks: If we wait a long time, what is the exact probability of seeing a specific number of dancers cross that fence? While scientists have solved this for a straight line (a one-dimensional fence), the math for a curved, circular fence in two dimensions was a mystery. The paper reveals that this problem is actually "integrable," meaning it has a neat, exact mathematical solution, rather than just an approximation.

The author's main discovery is a precise formula that predicts the "scaled cumulant generating function." In plain English, this is a master key that unlocks the full story of the current fluctuations. It tells us not just the average number of dancers crossing, but the odds of every possible outcome, from the boring average to the extremely rare, chaotic surges. The paper proves that even though the dance floor is two-dimensional and the fence is round, the most likely path the particles take to create a rare fluctuation is actually perfectly symmetrical, like ripples spreading out from a stone dropped in a pond. This allows the complex 2D problem to be squashed down into a simpler 1D problem.

To solve this, the author uses a clever trick involving "scattering," a concept usually reserved for how light bounces off objects. Here, the "light" is a mathematical wave that bounces off the invisible fence of the problem. By tracking how this wave scatters, the author derives a closed-form expression (a neat, finished equation) for the statistics. This isn't just a guess or a computer simulation; it is a rigorous mathematical proof. The result shows that the statistics depend on a specific combination of the densities inside and outside the circle, and the size of the circle relative to the time observed.

The paper also uncovers a beautiful symmetry in the results. It turns out that the probability of a certain number of particles leaving the circle is mathematically linked to the probability of the same number entering, just shifted by a specific factor related to the density difference. This confirms that the system follows deep, underlying laws of thermodynamics, even in these rare, non-equilibrium moments. Furthermore, the author shows that the complex behavior of the whole crowd can be thought of as a continuous stream of tiny, independent "channels," each contributing a small piece to the total count. This "channel representation" offers a new way to visualize how the geometry of the circle shapes the flow of particles.

In short, this paper takes a notoriously difficult problem in two-dimensional physics and solves it exactly. It proves that for a circular counting region, the messy, complex dance of particles can be described by a clean, elegant formula. This gives scientists a new benchmark for understanding how currents fluctuate in curved spaces, moving beyond the simple straight-line cases we've known for decades. The work doesn't just fill a gap in the math; it opens a door to understanding how geometry itself influences the flow of matter in the quantum and classical worlds.

Technical Summary: Nonisospectral Integrability and Exact Current Fluctuations in the Two-Dimensional SSEP

Problem Statement
This work addresses the exact calculation of annealed current fluctuations in the two-dimensional Symmetric Simple Exclusion Process (SSEP). Specifically, it investigates the large-deviation statistics of the net particle number decrease ( $Q_T$ ) within a finite disk of radius $R$ over a finite time $T$ . The system is initialized with a step-like density profile: a Bernoulli density $\rho_1$ inside the disk and $\rho_2$ outside. The counting boundary is "passive," meaning particles cross it freely in both directions, and the observable records only the algebraic net transfer.

The author notes that while exact results exist for one-dimensional SSEP systems (including step currents and tagged particles), extending these to higher dimensions with finite, curved counting regions is non-trivial. In two dimensions, the counting boundary is curved, and the optimal transport history is coupled to the spatial geometry, preventing a direct reduction to standard planar one-interface settings.

Methodology
The study employs Macroscopic Fluctuation Theory (MFT) combined with nonisospectral integrability techniques. The approach proceeds through the following steps:

Radial Reduction via Convexity: The author utilizes the convexity of the MFT action functional and the relative entropy of the initial state. By applying rotational averaging over $SO(2)$, they demonstrate that the global minimizer of the full two-dimensional variational problem can be chosen to be radially symmetric. This reduces the problem to a one-dimensional radial formulation.
Coordinate Transformation: For the radial problem in two dimensions, a specific change of variables to the area coordinate $x = r^2/2$ is introduced. This transformation flattens the radial measure and converts the MFT equations into a linear-coefficient divergence form, a property unique to the two-dimensional radial case.
Nonisospectral Lax Representation: The resulting radial MFT equations are shown to admit a nonisospectral Lax pair. Unlike standard one-dimensional scattering problems where the spectral parameter is constant, here the spectral parameter $k$ evolves in time ($k(t; z) = 1/(z+4it)$) to compensate for the explicit spatial factor in the radial equation.
Scattering and Factorization: The problem is mapped to a half-line scattering problem. By introducing a diagonalizing gauge, the author derives a scalar factorization problem for the transfer matrix. The boundary conditions at the temporal endpoints ( $t=0$ and $t=1$ ) lead to a factorization of the form $\Phi_-(s)\Phi_+(s) = 1 + \omega K(s)$ , where $\omega$ encodes the densities and counting field, and $K(s)$ encodes the geometric scale.
Spectral Reconstruction: Using Cauchy integral formulas on the scalar factors, the author reconstructs the scaled cumulant generating function (SCGF) and the optimal density profiles.

Key Results
The paper derives a closed-form expression for the scaled cumulant generating function (SCGF), $\mu_2^{ann}(\lambda; r_0)$ , where $r_0 = R/\sqrt{T}$ is the scaled radius.

Exact SCGF: The result is given by:
$\mu_2^{ann}(\lambda; r_0) = F_2(\omega(\lambda); r_0) = 4r_0^2 \int_0^1 \sqrt{1-y^2} \log\left(1 + \omega e^{-r_0^2 y^2}\right) dy$
where $\omega(\lambda) = \rho_1(1-\rho_2)(e^\lambda - 1) + \rho_2(1-\rho_1)(e^{-\lambda} - 1)$ .
Large Deviation Rate Function: The rate function $I(q)$ is obtained via the Legendre transform of the SCGF. The author shows that the statistics are non-Gaussian for finite $r_0$ , deviating from the Gaussian approximation even for moderate system sizes.
Fluctuation Symmetry: The SCGF satisfies a Gallavotti-Cohen-type fluctuation symmetry $\mu_2^{ann}(\lambda) = \mu_2^{ann}(A - \lambda)$ , where $A = \log[\rho_2(1-\rho_1) / \rho_1(1-\rho_2)]$ .
Channel Representation: The solution admits a "spectral channel" representation, expressing the free energy as an integral over a continuum of effective Bernoulli channels, with the measure $\nu_2$ encoding the radial geometry.
Optimal Profiles: The method allows for the reconstruction of the optimal initial and final density profiles. For positive net outflow, the optimal initial profile is deformed near the counting circle, while the final profile exhibits a depleted interior and enriched exterior.

Significance and Claims
The paper claims to provide the first exact solution for time-dependent full counting statistics in a two-dimensional diffusive system with a finite, curved counting boundary.

Integrability Structure: The work identifies a distinct nonisospectral integrable structure specific to the two-dimensional radial geometry. This structure relies on the unique properties of the radial MFT equations in $d=2$ and does not follow from standard one-dimensional SSEP solutions or Bethe-ansatz methods (which are not assumed to exist in 2D).
Geometric Dependence: The results explicitly demonstrate how current fluctuations depend on the geometry and codimension of the counting set, bridging the gap between 1D step-current results and higher-dimensional transport.
Methodological Advance: The paper establishes that convexity arguments can reduce complex multidimensional variational problems to radial ones, which can then be solved via nonisospectral scattering techniques.

The author notes that while the scalar factorization relies on the specific 2D radial reduction, the framework suggests potential extensions to other diffusive lattice gases and controlled coarse-graining scenarios, though rigorous analysis of higher-dimensional nonisospectral closures remains an open question.

Technical Summary: Nonisospectral Integrability and Exact Current Fluctuations in the Two-Dimensional SSEP

More like this