Self-organized hyperuniformity in a minimal model of population dynamics

Technical Summary: Self-Organized Hyperuniformity in a Minimal Model of Population Dynamics

Problem Statement
Hyperuniformity describes a state where a system is statistically disordered at short length scales yet exhibits vanishing density fluctuations at large scales, characterized by a structure factor $S(q) \to 0$ as $q \to 0$. While established in synthetic systems (e.g., jammed spheres, active rollers) and biological contexts (e.g., avian photoreceptors, leaf veins), the generic mechanisms driving hyperuniformity in non-equilibrium systems remain an open challenge. Previous models typically rely on fine-tuning control parameters to critical points of phase transitions (e.g., absorbing state transitions) or the emergence of effective center-of-mass conservation laws. This paper investigates whether hyperuniformity can arise in a biologically realistic population dynamics model without such conservation laws or parameter fine-tuning, driven instead by resource-mediated feedback.

Methodology
The authors introduce a minimal spatially extended model of population dynamics consisting of immobile agents competing for a shared, diffusive resource. The model comprises three core ingredients:

1. Stochastic Birth: Agents are introduced at a mean rate $\lambda$ per unit volume.
2. Viability and Death: An agent's internal "viability" coordinate $\nu_i(t)$ evolves according to a saddle-node (SN) bifurcation normal form, $\dot{\nu}_i = \nu_i^2 + \mu_i$. The bifurcation parameter $\mu_i$ is coupled linearly to the local resource concentration $c(x,t)$ via $\mu_i = c_{crit} - c(x,t)$. Agents are removed (die) when $\nu_i$ transitions from a stable negative fixed point (viable) to the positive half-line (non-viable).
3. Resource Dynamics: The resource field $c(x,t)$ is produced at a constant rate, consumed locally by viable agents, and diffuses with coefficient $D$. This creates a negative feedback loop: high resource levels allow agent accumulation, which increases consumption, thereby lowering resource levels and eventually triggering death.

The study employs a combination of stochastic simulations of the microscopic dynamics and a systematic coarse-graining procedure to derive a fluctuating hydrodynamic theory. This involves linearizing the dynamics around a mean-field steady state and analyzing the resulting system of delayed partial differential equations for density ($\delta\rho$) and resource ($\delta c$) fluctuations.

Key Contributions and Results

• Mechanism of Self-Organization: The model demonstrates that the system self-organizes into a steady state where the resource level asymptotically approaches the critical threshold $c_{crit}$ as the consumption rate decreases (or production increases). This occurs without fine-tuning an external parameter to a critical value; rather, the feedback loop drives the system to the critical point.
• Hyperuniformity without Conservation: Unlike previous non-equilibrium hyperuniform models, this system does not possess conservation laws (neither total mass nor center-of-mass conservation). Instead, hyperuniformity arises from the divergence of the interaction range mediated by the resource field. As the system approaches the critical threshold, the individual agent lifetime $\tau^*$ diverges ($\tau^* \sim \mu^*{}^{-1/2}$), allowing the resource field to adiabatically relax over a diverging screening length $\ell_D \rho^* \sim \mu^*{}^{-1/2}$.
• Hydrodynamic Theory: The authors derive a hydrodynamic description (Eqs. 10-11) that captures the system's behavior. The theory reveals that the resource field mediates long-range negative feedback interactions. In Fourier space, this results in a restoring kernel that diverges as $q^{-2}$, suppressing long-wavelength fluctuations.
• Scaling Behavior:
  • The structure factor exhibits Class I hyperuniformity, scaling as $S(q) \sim q^2$ for small wavenumbers $q$.
  • The number variance $\text{Var}[N(\ell)]$ scales sub-extensively as $\ell^{-1}$ (in 1D) over an intermediate range of length scales, before crossing over to a constant value at very large scales determined by the finite system size or the asymptotic limit.
  • The hyperuniform scaling is robust and insensitive to non-linearities in the model dynamics.
• Role of Time Delays: The analysis highlights that time-delay terms, arising naturally from the finite lifetime of agents, are essential for damping small-wavelength fluctuations and ensuring a bounded number variance. Without these delays (or an imposed singular death rate in a time-local alternative model), the system would not maintain bounded fluctuations at small scales.

Significance
The paper claims to identify a novel mechanism for hyperuniformity that is distinct from previously studied classes of non-equilibrium critical phenomena. Specifically:

1. No Phase Transition: The many-body system does not undergo a phase transition between distinct collective states; rather, it asymptotically approaches a critical resource threshold where interaction ranges diverge.
2. No Conservation Laws: Hyperuniformity is achieved without the effective emergence of center-of-mass conservation, challenging the notion that such conservation is a prerequisite for this state in non-equilibrium systems.
3. Biological Relevance: The mechanism is inspired by biological systems (e.g., plasma cell memory, cell cycle control) where survival depends on resource competition near a bifurcation. The findings suggest that hyperuniformity may be a generic feature of population dynamics where death processes are coupled through indirect resource-mediated interactions.
4. Theoretical Framework: The work provides a minimal, analytically tractable framework that bridges microscopic stochastic dynamics with macroscopic hydrodynamic behavior, offering a new avenue for understanding the origins of hyperuniformity in natural systems.

The authors conclude that their results open a new perspective on self-organized hyperuniformity, suggesting it can emerge from non-conservative population dynamics driven by resource competition and critical tuning of individual lifetimes, rather than requiring fine-tuned global control parameters or conservation laws.