Mass-Conserved Reaction-Diffusion Systems

• Mass-conserved reaction-diffusion systems are models where local reactions redistribute conserved mass, enabling the study of intracellular patterning and phase separation phenomena.
• They utilize mathematical constructs like reactive nullclines and flux-balance subspaces to determine pattern-forming instabilities and interface dynamics.
• These systems find applications in explaining cell polarity, self-assembly, and phase-separation in both biological and materials science contexts with robust numerical methods.

A mass-conserved reaction-diffusion (McRD) system is a class of reaction-diffusion equations in which interconversion of chemical species occurs exclusively via local reactions, but the total amount (mass) of all species is strictly conserved across the spatial domain. Such systems are central to the modeling of intracellular pattern formation, cell polarity, phase-separation dynamics, and more generally, to any context where global conservation laws tightly constrain local dynamics. These systems exhibit distinct dynamical phenomena—including mass-redistribution instabilities, geometric pattern selection via phase-plane mechanisms, interface motion governed by mesoscale laws, and coarsening dynamics analogous to phase separation—arising from the interplay of nonlinear kinetics and diffusive transport under conservation constraints (Brauns et al., 2018, Frey et al., 14 Dec 2025, Tateno et al., 2020).

1. Mathematical Formulation and Conservation Principles

The canonical example of a mass-conserved reaction-diffusion system is the two-component model on a bounded domain $\Omega\subset\mathbb{R}^n$ (typically with no-flux boundaries):

$$\begin{cases} \partial_t u(x,t) = D_u \nabla^2 u + f(u,v), \ \partial_t v(x,t) = D_v \nabla^2 v - f(u,v), \end{cases}$$

with the global conservation law

$$M := \int_\Omega [u(x,t) + v(x,t)]\,dx = \text{constant}.$$

Here, $u$ and $v$ represent densities of two chemical species (such as membrane-bound and cytosolic protein forms), $f(u,v)$ encodes local interconversion, and $D_u$, $D_v$ are their respective diffusion coefficients (Brauns et al., 2018).

General $m$-species systems allow for constitutive relations and stoichiometric weights, with the requirement that total mass is preserved: $$\frac{d}{dt} \sum_{i=1}^m \int_\Omega u_i(x, t)\, dx = \sum_{i=1}^m \int_\Omega f_i(u)\, dx = 0,$$ with reaction terms $f_i(u)$ satisfying $\sum_i f_i(u) = 0$ (Sun et al., 2021, Herberg et al., 2013).

In more complex geometries, such as coupled bulk-surface models and evolving domains, spatial integration includes both volume and surface components, and mass conservation is enforced through the accounting of all reservoirs and appropriate boundary fluxes (Mackenzie et al., 2019, Bao et al., 2014).

2. Phase-Space Structure and Pattern-Forming Instabilities

The geometric structure of mass-conserved systems is captured in phase space. The key objects are:

• Reactive Nullcline (NC): The locus $f(u,v) = 0$ in the $(u,v)$-plane defines the set of local reactive equilibria, parametrized by the conserved total density $n = u + v$.
• Flux-Balance Subspace (FBS): Steady states require local diffusive flux-balance, $D_u \nabla u + D_v \nabla v = 0$, which constrains steady state $(u,v)$ profiles to lie along lines (or hyperplanes for $N>2$) of the form $(D_u/D_v)u + v = \eta_0$ in phase space (Brauns et al., 2018, Frey et al., 14 Dec 2025).

Pattern formation is driven by interplay between these geometric structures. The classical "Turing" (mass-redistribution) instability in McRD systems occurs when the slope of the nullcline is more negative than the FBS: $$-\frac{f_u}{f_v} < -\frac{D_u}{D_v}\qquad\Longleftrightarrow\qquad \frac{f_u}{D_u} > \frac{f_v}{D_v},$$ where $f_u$ and $f_v$ are derivatives of $f$ with respect to $u$ and $v$ evaluated at homogeneous equilibrium. This geometric criterion characterizes the onset of patterning as a redistribution instability uniquely linked to the global conservation law (Brauns et al., 2018, Frey et al., 14 Dec 2025). Precise instability regions can be computed as subspaces in parameter space (e.g., diffusion coefficients, reaction rates).

3. Interface Dynamics, Wave-Pinning, and Coarsening

Upon nonlinear saturation, mass-conserved systems typically produce patterns consisting of nearly homogeneous "plateau" regions separated by narrow interfaces. The dynamics of these interfaces, including their motion, interaction, and stability, are strongly influenced by mass conservation. For two-plateau mesa states, the plateau concentrations are determined as intersections of FBS and NC in phase space; interfaces reside where these structures cross in the laterally unstable region of the NC (Brauns et al., 2018).

Key mesoscale features:

• Interface width:

$$\ell_{\text{interface}} \simeq \frac{\pi}{q_{\mathrm{max}}(n_0)},$$

with $q_{\mathrm{max}}^{2} = f_u/D_u - f_v/D_v$ at the inflection density $n_0$ (Brauns et al., 2018, Frey et al., 14 Dec 2025).

• Wave-pinning: Mass-conserved bistable systems generate stable stationary fronts (single transition layers) that pin at locations determined by global mass and Maxwell-type area constraints—unlike classical bistable RD systems, where fronts generically travel (Kuwamura et al., 2023, Miller et al., 2022).
• Coarsening: Differences in domain sizes drive mass flux through interfaces, leading to growth of larger domains at the expense of smaller ones, governed by laws analogous to Ostwald ripening; e.g., domain growth $\ell(t) \sim t^{1/3}$ in diffusion-limited regimes (Tateno et al., 2020, Frey et al., 14 Dec 2025).

In multiple spatial dimensions, interfaces acquire curvature, and interface dynamics involve effective surface tension and curvature-driven coarsening. The local "chemical potential" exhibits a jump proportional to interface curvature, reminiscent of the Young-Laplace law for classical phase separation (Tateno et al., 2020).

4. Analytical, Variational, and Dynamical Principles

The dynamical equations of McRD systems can be derived from energetic and variational principles for systems with reversibility/detailed balance. The Maxwell-Stefan framework with mass-action kinetics yields a strongly coupled quasilinear parabolic PDE, possessing a strict Lyapunov (free energy) functional: $$\mathcal{F}[c] = \int_\Omega \sum_{i} c_i(\ln(c_i/c_i^*) - 1)\, dx,$$ which decreases along solutions via a split into diffusion and reactive dissipation (Herberg et al., 2013, Wang et al., 2020).

Far-from-equilibrium systems, lacking detailed balance, may not possess classical free energy structure but can still admit gradient flow or generalized (entropy-dissipation) principles in limit regimes (Frey et al., 14 Dec 2025, Bao et al., 2014). In the singular limit of fast reactions and slow diffusion, McRD systems reduce to nonlinear diffusion equations for total density, and Maxwell constructions select admissible plateau values (Miller et al., 2022, Kuwamura et al., 2023).

Global existence and uniform bounds are ensured under quasi-positivity, mass control, and polynomial growth constraints. Techniques include $L^p$-energy methods, duality estimates, and the use of "intermediate-sum" and entropy-dissipation structures, for which detailed theorems cover critical and supercritical nonlinearities in both one and higher dimensions (Sun et al., 2021, Fellner et al., 2015, Lankeit et al., 2021, Morgan et al., 2021).

5. Applications: Pattern Formation in Biology and Materials

Mass-conserved reaction-diffusion systems are a canonical modeling framework for phenomena where total protein, ion, or particle content is conserved:

• Intracellular patterning: McRD models explain the formation and maintenance of protein patterns in the cytosol and on the membrane, as in the Min protein system of E. coli (pole-to-pole oscillations), Cdc42-mediated yeast polarity, and small GTPase-F-actin-based cell migration (Frey et al., 14 Dec 2025, Hughes et al., 10 Apr 2025).
• Cell polarization: Wave-pinning, bistability, and interface nucleation in McRD systems give paradigms for polarized domain formation, robust to cell size and geometry, and are validated by direct comparison to cellular imaging (Miller et al., 2022, Kuwamura et al., 2023).
• Coarsening and self-assembly: In synthetic and biomolecular settings, McRD systems recapitulate phase-separation dynamics, droplet coarsening, and interfacial phenomena in non-equilibrium contexts (Tateno et al., 2020).
• Bulk-surface and domain-coupled models: Coupling McRD equations on evolving bulk and membrane domains captures the interplay of diffusion, boundary kinetics, advection, and moving geometry in realistic cell and tissue morphodynamics (Mackenzie et al., 2019).

Well-posed numerical discretizations achieve robust mass conservation and accuracy for such systems on complex, evolving domains (Mackenzie et al., 2019).

6. Generalizations and Outlook

Recent theoretical advances extend McRD frameworks to:

• Multi-component and multi-conserved quantities: The FBS becomes a multidimensional hyperplane, and local equilibrium structure forms $(N-1)$-dimensional manifolds, with mass constraints inducing nontrivial high-dimensional pattern selection (Brauns et al., 2018).
• Weakly broken conservation: Addition of slow source/sink terms allows the core McRD dynamics to organize fast pattern formation before eventual arrest of coarsening at tunable scales (Brauns et al., 2018).
• Active matter and motility-induced phase separation: Phase-space and nullcline geometry generalize to cases with density-dependent diffusion and nontrivial mechanical feedback (Brauns et al., 2018).
• Irreversible and non-equilibrium networks: Analytical and energetic approaches accommodate generalized dissipation structures and nonequilibrium driving (Wang et al., 2020).

Open directions include sharp characterization of pattern-forming bifurcations, singular limit analysis in higher dimensions, and rigorous classification of metastable and multistable structures for generic nonlinearities.

---

References:

• "Phase-space geometry of mass-conserving reaction-diffusion dynamics" (Brauns et al., 2018)
• "Pattern Formation Beyond Turing: Physical Principles of Mass-Conserving Reaction--Diffusion Systems" (Frey et al., 14 Dec 2025)
• "Surface-tension-driven coarsening in mass-conserved reaction-diffusion systems" (Tateno et al., 2020)
• "Analysis of mass controlled reaction-diffusion systems with nonlinearities having critical growth rates" (Sun et al., 2021)
• "Reaction-diffusion systems of Maxwell-Stefan type with reversible mass-action kinetics" (Herberg et al., 2013)
• "Global classical solutions for mass-conserving, (super)-quadratic reaction-diffusion systems in three and higher space dimensions" (Fellner et al., 2015)
• "Global existence in reaction-diffusion systems with mass control under relaxed assumptions merely referring to cross-absorptive effects" (Lankeit et al., 2021)
• "A Conservative Finite Element ALE Scheme for Mass-Conserving Reaction-Diffusion Equations on Evolving Two-Dimensional Domains" (Mackenzie et al., 2019)
• "Existence of spiky stationary solutions to a mass-conserved reaction-diffusion model" (Morita et al., 2023)
• "A mass conserved reaction-diffusion system reveals switching between coexisting polar and oscillatory cell motility states" (Hughes et al., 10 Apr 2025)
• "Single Transition Layer in Mass-Conserving Reaction-Diffusion Systems with Bistable Nonlinearity" (Kuwamura et al., 2023)
• "Generation and motion of interfaces in a mass-conserving reaction-diffusion system" (Miller et al., 2022)
• "Well-posedness and exponential equilibration of a volume-surface reaction-diffusion system with nonlinear boundary coupling" (Bao et al., 2014)
• "Field Theory of Reaction-Diffusion: Mass Action with an Energetic Variational Approach" (Wang et al., 2020)
• "Global Existence of Solutions to Reaction Diffusion Systems with Mass Transport Type Boundary Conditions" (Sharma et al., 2015)
• "Mass conservative reaction diffusion systems describing cell polarity" (Latos et al., 2020)
• "Global well-posedness for volume-surface reaction-diffusion systems" (Morgan et al., 2021)
• "Global dynamics and spectrum comparison of a reaction-diffusion system with mass conservation" (Latos et al., 2015)