Effect of Weak Non-Conservative Dynamics on Pattern Formation in Scalar Active Matter

Technical Summary: Effect of Weak Non-Conservative Dynamics on Pattern Formation in Scalar Active Matter

Problem Statement
Active matter systems, ranging from bacterial colonies to synthetic active colloids, are inherently driven far from thermodynamic equilibrium by continuous energy consumption. While standard theoretical frameworks like Active Model B (AMB) describe phase separation in scalar active matter under strict mass conservation, many biological and synthetic systems exhibit weak violations of mass conservation due to particle birth-death processes, growth, degradation, or chemical exchange with reservoirs. The specific influence of these weak non-conservative dynamics on the late-time morphology, pattern selection, and coarsening behavior of scalar active systems remains underexplored. This work addresses the gap in understanding how such non-conservative reactions modify the nonlinear evolution and steady-state structures of active phase separation.

Methodology
The author constructs a minimal continuum model by extending Active Model B with a weak non-conservative reaction term. The system is described by a coarse-grained scalar order parameter field, $\phi(\mathbf{r}, t)$, representing particle density. The dynamics are governed by the equation of motion:
$$\frac{\partial \phi}{\partial t} = -\Gamma_\phi \nabla \cdot \mathbf{J} + \gamma \phi \left(1 - \frac{\phi}{\phi_c}\right)\left(1 + \frac{\phi}{\phi_c}\right)$$
where $\mathbf{J} = -\nabla \mu + \sqrt{2\Gamma_\phi \phi}\eta$ is the current, and the chemical potential $\mu$ includes a standard Landau free energy contribution and an active term $\mu_{active} = \zeta |\nabla \phi|^2$ that breaks time-reversal symmetry. The parameter $\zeta$ controls activity strength, while $\gamma$ (with $\gamma \leq 0$) controls the strength of the non-conservative reaction.

The study employs two primary approaches:

1. Linear Stability Analysis: To determine the initial instability and the dispersion relation $\sigma(q)$, analyzing how the reaction term and activity affect the growth rate of perturbations around a homogeneous steady state.
2. Numerical Simulations: Extensive simulations on a $512 \times 512$ lattice with periodic boundary conditions to observe late-time morphologies. The system is initialized with a uniform density and evolved to steady states.
3. Structural Characterization: Quantitative analysis using equal-time two-point correlation functions $C(r)$, static structure factors $S(k)$, cluster statistics (largest cluster fraction $L_c$, number of isolated clusters $N_c$, circularity $C$), droplet size distributions, and the hexatic bond-orientational order parameter $\psi_6$.

Key Contributions and Results

• Arrest of Coarsening and Microphase Separation: The primary finding is that even weak non-conservative reactions ($\gamma < 0$) fundamentally arrest the coarsening process typical of conserved active systems. Instead of evolving toward bulk phase separation, the system stabilizes into nonequilibrium microphase-separated states with a well-defined intrinsic length scale.
• Reversal of Morphological Transitions: In the conserved limit ($\gamma = 0$), increasing the magnitude of negative activity ($\zeta$) drives a transition from droplets to worm-like structures and finally to interconnected labyrinthine patterns (spinodal decomposition). In contrast, the presence of a weak non-conservative reaction reverses this sequence. As activity increases, the system transitions from labyrinthine patterns to worm-like structures and eventually to isolated, polydisperse droplets.
• Distinct Roles of Parameters:
  • Activity ($\zeta$): Predominantly controls the global domain morphology and the characteristic wavelength of the patterns. Increasing $|\zeta|$ reduces the characteristic length scale ($\lambda^*$) and drives the transition from labyrinths to isolated droplets.
  • Reaction Coefficient ($\gamma$): Primarily promotes microphase separation and enhances local hexagonal ordering. Increasing the magnitude of $\gamma$ suppresses domain coalescence, increases the number of isolated clusters, and reduces the average droplet size, leading to a narrower size distribution.
• Linear Stability Insights: The analysis reveals that the reaction term enters the dispersion relation as a $q$-independent constant, shifting the growth rate curve vertically without altering the position of the most unstable wavenumber $q^*$. The active contribution ($\zeta |\nabla \phi|^2$) does not appear in the linear order, meaning neither term modifies the initial wavelength selection. The qualitative differences in morphology arise from the nonlinear evolution following the initial instability.
• Emergence of Hexatic Order: In the droplet phase, the system exhibits local hexagonal packing. The non-conservative reaction term enhances this local hexatic order ($\psi_6 \to 1$) by regulating the droplet size distribution and facilitating size homogenization, which allows for ordered packing. This is distinct from classical Ostwald ripening, as the non-conserved nature of the order parameter prevents unlimited growth.
• Scaling Behavior: Scaled correlation functions $C(r/\lambda^*)$ collapse onto a universal master curve for varying $\zeta$ and $\gamma$, indicating that while the parameters rescale the characteristic length, the underlying normalized spatial organization remains robust.

Significance
The paper demonstrates that weak violations of mass conservation are not merely perturbations but fundamentally alter the nonlinear coarsening dynamics of active phase separation. By introducing a minimal framework that couples active transport with weak non-conservative reactions, the study provides a mechanism to stabilize microphase-separated states that would otherwise coarsen into bulk phases. The results suggest that non-conservative dynamics serve as a simple yet effective mechanism for controlling nonequilibrium self-organization. This framework offers a basis for understanding pattern formation in biological systems (e.g., bacterial colonies with birth-death cycles) and for designing synthetic active materials with tunable microstructures, such as isolated droplets with specific ordering, without requiring complex external forcing.