Stationary densities and delocalized domain walls in asymmetric exclusion processes competing for finite pools of resources

This paper demonstrates that a pair of antiparallel asymmetric exclusion processes competing for finite particle reservoirs exhibits an extended region of delocalized domain walls in its phase diagram, leading to large particle number fluctuations in the thermodynamic limit, a behavior distinct from standard TASEP models where such walls typically exist only along a specific line.

The Gist

The Big Picture: A Traffic Jam with a Twist

Imagine you are managing a busy highway system, but with a very specific rule: there is a limited number of cars in the entire universe. You can't just buy more cars, and you can't throw any away. The total number of cars is fixed.

In this paper, the scientists study a model of two parallel highways (let's call them Lane A and Lane B) where cars only move in one direction.

• Lane A goes from Left to Right.
• Lane B goes from Right to Left.

At the ends of these highways, there are two giant parking lots (Reservoir 1 and Reservoir 2). Cars enter the highways from these lots and exit back into them.

The Catch: The rate at which cars enter or leave the highway depends on how many cars are currently sitting in the parking lots.

• If a parking lot is full, it's easy for cars to jump onto the highway (high entry rate).
• If a parking lot is empty, it's hard for cars to leave the highway and get back in (low exit rate).

The researchers wanted to see: What happens to the traffic flow when the total number of cars is limited and shared between these two lanes?

---

The Main Discovery: The "Ghost" Traffic Jam

In standard traffic models (where you can have infinite cars), traffic jams usually happen in one of two ways:

1. Low Density: The road is empty; cars zoom along.
2. High Density: The road is packed; cars are stuck.
3. The Shockwave (Domain Wall): Sometimes, you get a sudden transition from "empty road" to "packed road." In normal models, this shockwave is usually pinned to a specific spot (like a construction zone) or exists only on a very thin, specific line of control.

The Surprise:
In this specific model with limited resources, the scientists found a strange new state called Delocalized Domain Walls (DDW).

The Analogy:
Imagine a traffic jam that isn't stuck in one place. Instead, the boundary between "free-flowing traffic" and "gridlock" is shaking back and forth wildly across the entire highway.

• One moment, the jam is at the start of the road.
• The next moment, it's at the end.
• Over time, if you take a photo, the jam looks like a fuzzy, diagonal blur stretching across the whole road.

Usually, in physics, these "fuzzy" states only happen if you tune your controls (like the speed limit or entry rate) to a perfect, razor-thin line. If you move even a tiny bit away from that line, the fuzziness disappears.

The Breakthrough:
In this paper, the scientists found that this "fuzzy, shaking jam" state exists over a huge, wide area of the control settings. It's not a razor-thin line; it's a wide, open field. This means that for a wide range of conditions, the traffic density on the road will fluctuate wildly and unpredictably, even if the system is huge and stable.

---

Why This Matters: The "Balancing Act"

Here is the most counter-intuitive part of the story:

Even though the traffic on the highways is chaotic, shaking, and fluctuating wildly, the parking lots remain perfectly calm.

• The Highways: The number of cars on the road jumps up and down significantly.
• The Parking Lots: The number of cars in the lots stays almost exactly the same.

The Metaphor:
Think of the highways as a sponge and the parking lots as a bucket of water.
If you squeeze the sponge (the highway), water squirts out and moves around chaotically inside the sponge. However, because the total amount of water is fixed, the bucket doesn't need to change its water level much to compensate. The chaos is contained within the sponge.

The researchers proved mathematically that even as the system gets infinitely large (the "thermodynamic limit"), the traffic on the road stays chaotic, but the parking lots become perfectly stable.

---

Real-World Connections

Why should a regular person care about mathy traffic jams? The paper suggests this model applies to biology:

• Protein Synthesis: Imagine a cell making proteins. The "highways" are strands of mRNA, and the "cars" are ribosomes (the machines that build proteins). The "parking lots" are the supply of free ribosomes floating in the cell.
• The Lesson: If the cell has a limited supply of ribosomes, the process of building proteins might not be smooth. Instead, the ribosomes might clump up and spread out in a chaotic, fluctuating wave along the mRNA strand. This could mean that protein production isn't steady; it might surge and dip wildly, even if the cell's overall resources are constant.

Summary in Three Points

1. The Setup: Two one-way roads sharing a fixed pool of cars.
2. The Surprise: Instead of a steady flow or a fixed jam, the system creates a "shaking" jam that moves back and forth across the whole road.
3. The Result: This shaking happens over a wide range of conditions (not just a tiny line), and while the road is chaotic, the supply of cars in the parking lots remains perfectly stable.

This discovery changes how we understand how limited resources affect flow in systems ranging from traffic on roads to the machinery of life inside our cells.

Technical Summary

1. Problem Statement and Motivation

The paper investigates the non-equilibrium steady states of a system composed of two Totally Asymmetric Simple Exclusion Processes (TASEP) lanes coupled antiparallelly to two particle reservoirs with finite capacity.

• Context: Standard TASEP models usually assume open boundaries with fixed entry/exit rates or infinite reservoirs. However, many real-world systems (e.g., ribosome translation in cells, traffic in closed networks) operate under finite resource constraints.
• Specific Gap: Previous studies on TASEP with finite resources often assumed a single reservoir or constant exit rates independent of reservoir population. This work explores a model (originally proposed in Haldar et al., Phys. Rev. E 111, 014154 (2025)) where both entry and exit rates are dynamically controlled by the instantaneous populations of the reservoirs.
• Key Question: How does the competition for a finite pool of particles between two coupled TASEP lanes affect stationary densities, phase diagrams, and the nature of domain walls (shocks), particularly when the system parameters are symmetric?

2. Model Description

The system consists of two 1D lattices ($T_1$ and $T_2$) of length $L$ and two point-like reservoirs ($R_1$ and $R_2$) with particle counts $N_1$ and $N_2$.

• Dynamics:
  • Particles hop unidirectionally ($T_1$: left-to-right; $T_2$: right-to-left) with rate 1.
  • Hard-core exclusion: Max one particle per site.
  • Finite Resources: The total particle number $N_0$ is strictly conserved: $N_0 = N_1 + N_2 + \sum n_j$.
• Rate Functions: The effective entry ($\alpha_{eff}$) and exit ($\beta_{eff}$) rates depend on reservoir populations:
  • $\alpha_{eff}^{(1)} = \alpha_1 (N_1/L)$, $\beta_{eff}^{(1)} = \beta_1 (1 - N_2/L)$
  • $\alpha_{eff}^{(2)} = \alpha_2 (N_2/L)$, $\beta_{eff}^{(2)} = \beta_2 (1 - N_1/L)$
  • Note: Entry rates increase with reservoir population; exit rates decrease (crowding effect).
• Parameter Space: The authors simplify the study by assuming symmetry ($\alpha_1 = \alpha_2 \equiv \alpha$ and $\beta_1 = \beta_2 \equiv \beta$). The system is governed by three parameters: $\alpha$, $\beta$, and the filling factor $\mu = N_0 / 2L$.

3. Methodology

The authors employ a dual approach:

1. Mean-Field Theory (MFT):
   • Neglects spatial correlations to derive equations of motion for bulk densities.
   • Utilizes Particle-Hole Symmetry transformations to map the phase diagram for $\mu > 1$ from the $\mu < 1$ regime.
   • Derives analytical expressions for phase boundaries and steady-state densities.
2. Monte Carlo Simulations (MCS):
   • Extensive numerical simulations using random sequential updates.
   • System sizes up to $L=1000$ and $2 \times 10^9$ time steps.
   • Used to validate MFT predictions, calculate density profiles, and analyze fluctuations (standard deviations) in the thermodynamic limit.

4. Key Contributions and Results

A. Symmetry of Stationary States

Contrary to previous studies (e.g., Ref. [53]) where asymmetric phases (e.g., Low-Density in one lane, High-Density in the other) were possible, this symmetric parameter space forces statistically identical stationary densities in both lanes ($\rho^{(1)} = \rho^{(2)}$).

• Allowed Phases: Only symmetric combinations exist:
  • LD-LD: Low density in both lanes.
  • HD-HD: High density in both lanes.
  • MC-MC: Maximal current in both lanes.
  • DW-DW: Domain walls in both lanes.
• Exclusion of Asymmetry: The authors prove analytically that asymmetric phases (LD-HD, LD-DW, etc.) are impossible under these symmetric boundary conditions.

B. Delocalized Domain Walls (DDW)

The most significant finding is the nature of the DW-DW phase.

• Extended Region: Unlike standard TASEP models where delocalized domain walls exist only along a specific line in parameter space (coexistence line), here the DDW phase occupies an extended region in the $(\alpha, \beta)$ plane.
• Mechanism: Due to global particle conservation and the freedom of the domain walls to shift relative to each other (one shifts left, the other shifts right to conserve $N_0$), the walls are not pinned.
• Fluctuations: The position of the domain walls fluctuates over a distance scaling with $L$. Consequently, the stationary densities exhibit large fluctuations that survive in the thermodynamic limit ($L \to \infty$). This is a rare feature in driven diffusive systems.

C. Phase Diagram Topology

The phase diagrams in the $(\alpha, \beta)$ plane vary with the filling factor $\mu$:

• $0 < \mu < 1/2$: Only LD-LD and DW-DW phases exist.
• $1/2 < \mu < 3/2$: All four phases (LD-LD, HD-HD, MC-MC, DW-DW) coexist. They meet at a multicritical point $(\alpha_c, \beta_c)$.
• $3/2 < \mu < 2$: Only HD-HD and DW-DW phases exist.
• Transitions: All phase transitions are continuous (second-order), characterized by the bulk density as the order parameter. The multicritical point traces a hyperbolic trajectory as $\mu$ varies.

D. Fluctuation Analysis

• TASEP Lanes: In the DW-DW phase, the relative fluctuations of the average density do not vanish as $L \to \infty$. This implies macroscopic instability in particle numbers within the lanes.
• Reservoirs: Surprisingly, despite the large fluctuations in the lanes, the relative fluctuations in the reservoir populations ($N_1, N_2$) vanish in the thermodynamic limit. The reservoirs act as a stabilizing buffer, maintaining equal populations ($N_1 = N_2$) on average.

5. Significance and Implications

• Theoretical Novelty: The discovery of an extended region for delocalized domain walls challenges the conventional understanding of TASEP phase diagrams, where such features are typically confined to lines. It demonstrates that finite resource constraints coupled with symmetric dynamics can induce robust macroscopic fluctuations.
• Biological Relevance: The model serves as a minimal description for ribosome translation on mRNA strands with a fixed pool of ribosomes. The results suggest that under certain entry/exit rate conditions, ribosome density along the mRNA could exhibit large, persistent fluctuations even in large systems, potentially affecting protein synthesis efficiency and regulation.
• Traffic Dynamics: The findings apply to closed traffic networks with a fixed number of vehicles, predicting regimes where traffic jams (domain walls) are not static but fluctuate significantly across the entire network.
• Methodological Insight: The study highlights that Mean-Field Theory can yield accurate results even in systems with global conservation laws, provided the reservoirs effectively decouple the conservation constraints of individual lanes, making them behave similarly to open TASEPs.

Conclusion

This paper establishes that a pair of TASEP lanes competing for finite, symmetric resources exhibits a unique phase behavior characterized by identical stationary densities in both lanes and an extended phase of delocalized domain walls. This leads to non-vanishing density fluctuations in the thermodynamic limit, a phenomenon with profound implications for understanding resource-limited transport in biological and physical systems.