Non-stationary current fluctuations in 1D boundary-driven diffusive systems via Macroscopic Fluctuation Theory

This study extends Macroscopic Fluctuation Theory to non-steady-state processes by deriving exact expressions for current variance and the cumulant generating function in one-dimensional boundary-driven diffusive systems, demonstrating that the framework can quantitatively describe current fluctuations during relaxation toward a steady state.

The Gist

Imagine a long, narrow hallway filled with people (particles) who are trying to move from one end to the other. At the left door, people are constantly entering and leaving based on how crowded the room is outside. The same happens at the right door. This is a "boundary-driven system."

Usually, scientists study what happens after everyone has settled into a steady rhythm—a "non-equilibrium steady state" (NESS). But this paper asks a different question: What happens while the system is still waking up? What are the chaotic fluctuations of people moving through the hallway before the steady rhythm is established?

The authors use a powerful mathematical toolkit called Macroscopic Fluctuation Theory (MFT). Think of MFT as a "weather forecast" for crowds. Instead of tracking every single person, it predicts the probability of different crowd patterns and flow rates. While MFT has been great at predicting steady weather, this paper applies it to the "stormy" period of relaxation.

Here is a breakdown of their findings using simple analogies:

1. The Two Types of "Starting Lines"

The researchers looked at two different ways the hallway could start out, which changes how the crowd behaves:

• The "Annealed" Start (The Party): Imagine the people are already in the hallway, but they are jittery and shifting around randomly due to thermal energy (like a party where everyone is dancing). The starting positions are fluid and fluctuating.
• The "Quenched" Start (The Frozen Line): Imagine the people are frozen in place at the start. Their positions are fixed and rigid, with no initial jitter.

The Finding: The paper proves that the "Party" start (Annealed) leads to more chaos (higher variance) in how many people cross a specific point than the "Frozen Line" start (Quenched). Because the people were already wiggling around at the start, the total number of people passing through fluctuates more wildly.

2. The "Traffic Jam" vs. "Free Flow" (Diffusion Models)

They tested their theory on two specific types of "crowds":

• The "Exclusion" Crowd (SEP): Imagine people in a hallway who cannot pass each other. If you are in front of someone, you are stuck. This is like a single-file line.
• The "Independent" Crowd (IRW/RBM): Imagine people in a hallway who can walk through each other like ghosts, or a crowd of non-interacting Brownian particles.

The Finding: In the "Exclusion" crowd, the movement is slower and less fluctuating because people block each other. In the "Independent" crowd, people move more freely, leading to larger fluctuations. The authors derived exact formulas showing exactly how much the "traffic jam" effect suppresses the noise compared to the "ghost" crowd.

3. The "Time Travel" of Fluctuations

One of the most interesting discoveries is how the "noise" (fluctuations) changes over time.

• Early Times (The Short Hop): If you watch for a very short time, the crowd hasn't felt the influence of the far end of the hallway yet. It acts like an infinite hallway with just one door. The fluctuations grow slowly (proportional to the square root of time, $\sqrt{T}$).
• Late Times (The Long Haul): If you watch for a long time, the crowd feels the pressure from both doors. The system settles into a steady flow. Now, the fluctuations grow linearly with time ($T$).

The Finding: The paper maps out the exact "crossover" moment where the system stops acting like a short hop and starts acting like a long, steady flow. They showed that the mathematical framework (MFT) can perfectly describe this transition, even when the initial conditions and the boundary doors are interacting in complex ways.

4. The "Magic Trick" of Math (RBM)

For a specific type of crowd called Reflective Brownian Motion (RBM)—which is like a crowd of non-interacting particles bouncing off walls—the authors performed a "magic trick." They used a mathematical transformation (Cole-Hopf) to turn a very messy, non-linear equation into a simple, linear one.

The Result: This allowed them to write down the exact formula for the probability of any specific flow rate. They didn't just guess; they solved it perfectly. They showed that the statistics of this crowd are essentially the difference between two simple "coin-flip" processes (Poisson processes), making the complex behavior surprisingly simple to describe.

Summary

In short, this paper takes a sophisticated theory used for steady states and successfully applies it to the messy, chaotic period of relaxation.

• They proved that how you start (frozen vs. jittery) changes how much the flow fluctuates.
• They showed that crowd rules (blocking vs. passing) change the size of those fluctuations.
• They mapped out exactly how the system transitions from a short-term chaotic state to a long-term steady state.

The paper concludes that Macroscopic Fluctuation Theory is not just for steady states; it is a robust, universal tool for understanding how physical systems relax and find their balance, even when they are far from equilibrium.

Technical Summary

Technical Summary: Non-stationary current fluctuations in 1D boundary-driven diffusive systems via Macroscopic Fluctuation Theory

Problem Statement
Macroscopic Fluctuation Theory (MFT) has established itself as a robust framework for analyzing large deviations in non-equilibrium steady states (NESS) and infinite non-stationary systems. However, its application to finite non-stationary systems—specifically the relaxation process of a system coupled to particle reservoirs at both ends—remains limited. In such finite systems, the statistical behavior of the integrated current $Q_T$ is governed by a complex interplay between initial conditions and boundary driving. As the system relaxes toward a steady state, the scaling of current fluctuations transitions from an initial diffusive regime ($O(\sqrt{T})$) to a steady-state linear regime ($O(T)$). Understanding the large deviation properties in this crossover region, where initial and boundary contributions are coupled, has been a significant gap in the theoretical understanding of non-equilibrium relaxation dynamics.

Methodology
The authors apply the MFT framework to one-dimensional boundary-driven diffusive systems. The methodology involves the following steps:

1. Formalism Derivation: Starting from stochastic lattice gases (e.g., Symmetric Simple Exclusion Process - SEP, Independent Random Walk - IRW), the authors derive the MFT action and the associated coupled non-linear partial differential equations (MFT equations) for the density field $\rho(x,t)$ and the auxiliary field $H(x,t)$. This derivation utilizes the Martin-Siggia-Rose (MSR) path integral formalism and the hydrodynamic limit ($\Lambda \to \infty$).
2. Boundary Conditions: The study explicitly formulates boundary conditions for three regimes: fast-boundary ($\theta < 1$), slow-boundary ($\theta > 1$), and marginal-boundary ($\theta = 1$). The marginal regime is treated as the general case, with fast and slow regimes emerging as asymptotic limits.
3. Perturbative Expansion: For systems with a constant diffusion coefficient $D$ and arbitrary mobility $\sigma(\rho)$, the authors perform a perturbative expansion with respect to the counting field $\lambda$. This allows for the analytical calculation of the first and second cumulants (mean and variance) of the integrated current without requiring a full solution to the non-linear MFT equations.
4. Exact Solution for RBM: For Reflective Brownian Motion (RBM), which corresponds to the Independent Random Walk (IRW) with transport coefficients $D(\rho)=1$ and $\sigma(\rho)=2\rho$, the authors utilize the Cole-Hopf transformation ($P=e^H, Q=\rho e^{-H}$). This transformation linearizes the coupled MFT equations into independent forward and backward diffusion equations, enabling an exact analytical derivation of the full Scaled Cumulant Generating Function (SCGF) for arbitrary boundary speeds and initial conditions.

Key Contributions and Results

• Current Variance for General $\sigma(\rho)$:
  The authors derive exact analytical expressions for the current variance in systems with constant $D$ and arbitrary $\sigma(\rho)$. They distinguish between two initial conditions:

  • Annealed: Initial particle positions fluctuate according to thermal equilibrium.
  • Quenched: The initial density profile is fixed.
    The results show that the current variance for the annealed case is strictly larger than that for the quenched case ($\langle Q_T^2 \rangle_{c,A} > \langle Q_T^2 \rangle_{c,Q}$), quantifying the contribution of initial thermal fluctuations to the total current fluctuations. The derived variance exhibits a crossover from $O(\sqrt{T})$ at short times ($T \ll L^2$) to $O(T)$ at long times ($T \gg L^2$), consistent with the transition from single-reservoir dominance to NESS behavior. These analytical results are validated against Monte Carlo simulations of the SEP, showing excellent quantitative agreement.

• Exact SCGF for RBM:
  The paper provides an exact derivation of the SCGF for the integrated current in RBM under both annealed and quenched initial conditions and for arbitrary boundary coupling strengths.

  • The solution reveals that the statistics of the integrated current follow a Skellam distribution (the difference of two independent Poisson processes).
  • The authors analyze the Large Deviation Function (LDF) across different time regimes. They observe that in the negative current regime, the LDF for the quenched case increases more rapidly than for the annealed case, indicating that negative current fluctuations are suppressed when the initial configuration is fixed.
  • The study confirms that the marginal-boundary regime captures the full physics, reproducing the fast-boundary Dirichlet limits ($\rho(0,t)=\rho_L$) as the boundary coupling strength tends to infinity.

• Consistency Checks:
  The derived results are shown to be consistent with:

  • Previous results for semi-infinite systems in the limit $L \to \infty$.
  • Known results for NESS in the long-time limit ($T \to \infty$), recovering the additivity principle for the steady-state current variance.

Significance and Claims
The paper claims to demonstrate that non-stationary current fluctuations during the relaxation process of finite systems can be systematically and quantitatively described within the MFT framework. By extending MFT beyond steady states and infinite systems, the authors provide a unified description of how initial conditions and boundary driving jointly determine the statistical properties of current fluctuations.

The work establishes MFT as a robust tool for analyzing the crossover between transient and steady-state dynamics. Specifically, the exact solvability of the RBM case serves as a benchmark, confirming that MFT can capture the microscopic dynamics of relaxation processes accurately. The authors note that while the exact SCGF derivation is specific to RBM (due to the linearity afforded by the Cole-Hopf transformation), the perturbative approach for current variance is applicable to a broader class of diffusive systems with arbitrary mobility.

The paper concludes by identifying the extension of these exact methods to systems with strong many-body interactions, such as the finite-system SEP, as a primary future challenge. This would require mapping the finite-system MFT equations to integrable systems, a task currently solved only for infinite systems.