Generalized Active Model B+ (AMB+)

• Generalized Active Model B+ is a framework for non-equilibrium phase separation in active matter, incorporating conserved order parameter dynamics.
• It breaks time-reversal symmetry by introducing scalar and rotational active currents that modify phase behavior and interface dynamics.
• Extensions of AMB+ include reaction-driven terms and off-critical effects, resulting in tunable hierarchical microphase structures.

Generalized Active Model B+ (AMB+), introduced as an extension of the passive Model B for conserved order parameter dynamics, is a paradigmatic framework for studying non-equilibrium phase separation in active matter. AMB+ systematically incorporates time-reversal symmetry (TRS) breaking through two distinct active current contributions: a scalar, curl-free current proportional to λ, and a rotational, nonintegrable term proportional to ξ. The model describes the evolution of a conserved scalar field φ(r, t), relevant for binary mixtures and active Brownian particle suspensions, and captures both macroscale (bulk) and microscale (arrested) phase separation regimes depending on these active couplings. Recent generalizations have extended AMB+ to include reactions and off-critical effects, yielding a rich phenomenology including tunable pattern formation, amplitude equation bifurcations, and multiscale hierarchical structures (Yadav et al., 17 Jun 2025, Mondal et al., 12 Jan 2026, Li et al., 2021).

1. Formal Structure of Generalized AMB⁺

The dynamics of AMB+ are governed by the continuity equation for a conserved order parameter φ(r, t): $\partial_t\phi = -\nabla\cdot \boldsymbol{J}$ where the total current is decomposed into: $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$ The equilibrium chemical potential $\mu_E$ derives from a Ginzburg–Landau functional: $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$

$$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$

The λ-term ($J_\lambda = -\lambda\nabla|\nabla\phi|^2$) is the lowest-order scalar correction that breaks TRS yet remains gradient-driven, while the ξ-term ($J_\xi = \xi (\nabla^2\phi) \nabla\phi$) fundamentally breaks this structure by introducing a rotational, nonintegrable contribution. In the presence of chemical reactions or off-criticality, further extensions incorporate a quadratic term $g\phi^2$ and a reaction sink $-\Gamma\phi$, modifying the chemical potential and mass conservation, respectively (Mondal et al., 12 Jan 2026, Li et al., 2021).

2. Physical Origin and Interpretation of Active Currents

The two active terms in AMB+ originate from different physical mechanisms of TRS breaking:

• λ-term (rotation-free): This is a scalar correction, $O(\nabla^4\phi^2)$, representing the lowest-order TRS-violating contribution that maintains curl-free structure. Physically, it modifies the effective chemical potential, shifting coexistence densities and giving rise to non-equilibrium steady states. It can be absorbed into a nonequilibrium chemical potential $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$0.
• ξ-term (rotational current): This term, proportional to $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$1, cannot be written as the gradient of a scalar. It introduces interfacial tangential flows and generates circulating surface currents, leading to fundamentally nonequilibrium behaviors, especially at interfaces.

The interplay between these two terms enables the model to capture regimes where conventional Ostwald ripening dynamics are either altered (forward Ostwald, λ–ξ/2 > 0) or reversed (reverse Ostwald, λ–ξ/2 < 0), resulting in either coarsening to bulk phase separation or arrested microphase separation, respectively (Yadav et al., 17 Jun 2025).

3. Macroscale and Microscale Phase Separation Kinetics

Phase separation outcomes in AMB+ depend sensitively on the relative magnitudes and signs of $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$2 and $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$3:

• Macroscale Phase Separation (MPS): For $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$4, standard coarsening occurs with two distinct growth regimes:
  • Early time: Lifshitz–Slyozov scaling $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$5 characteristic of bulk diffusion.
  • Late time: Surface-diffusion dominated regime with $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$6 due to strong interfacial currents. Crossover time $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$7 for $$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$8.
  • Correlation functions exhibit dynamical scaling with superuniversal morphology at early times, though late-time morphologies reveal parameter dependence due to off-criticality.
• Microscale Phase Separation ($$\boldsymbol{J} = \underbrace{ -\nabla\mu_E }_{\text{Model B}} \underbrace{ -\lambda\nabla|\nabla\phi|^2 }_{\text{rotation-free active}} \underbrace{ +\xi\,(\nabla^2\phi)\nabla\phi }_{\text{rotational active}}$$9PS): For $\mu_E$0, the system undergoes reverse Ostwald ripening: small droplets grow at the expense of larger ones, leading to arrested coarsening and a steady state with finite length $\mu_E$1. For fixed $\mu_E$2, $\mu_E$3, and the late-time steady state forms a hexagonal droplet array.

These kinetic regimes and morphologies have been established via large-scale simulations using forward-Euler discretization and periodic domains, with systematic averaging over many realizations (Yadav et al., 17 Jun 2025).

4. Reaction-Driven Pattern Formation and Amplitude Equation

Introducing a reversible reaction $\mu_E$4 at rate $\mu_E$5 adds a nonconservative sink $\mu_E$6, breaking the conservation of $\mu_E$7. This induces several new phenomena:

• Linear stability analysis yields a dispersion relation $\mu_E$8, with a critical rate $\mu_E$9 and corresponding preferred wavenumber $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$0.
• Pattern formation below threshold ($$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$1) is characterized by spatial modulations at $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$2.

A multiscale analysis leads to a complex amplitude equation for the envelope $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$3 of roll patterns: $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$4 where the sign and magnitude of $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$5 determine the bifurcation nature:

• Always supercritical for $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$6.
• For $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$7, the transition can be subcritical, with the boundary given by $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$8.
• The Eckhaus instability band is determined independently of $$F[\phi] = \int d^d r\, \left[ -\tfrac{1}{2}\phi^2 + \tfrac{1}{4}\phi^4 + \tfrac{1}{2}|\nabla\phi|^2 \right ]$$9 (Mondal et al., 12 Jan 2026).

This amplitude equation recovers several well-known models in appropriate limits (passive Cahn–Hilliard, asymmetric Cahn–Hilliard, AMB+, etc.).

5. Model AB+ and Hierarchical Microphase Separation

The further generalization to "Model AB+" (editor's term: AMB+ with both conservative and nonconservative TRS breaking) incorporates both diffusive/active and reaction/chemical mechanisms. Its governing equation includes: $$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$0 where the current $$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$1 contains the AMB+ terms, and the reaction field $$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$2 enforces nonconserved dynamics.

This setting yields hierarchical microphase separation:

• Small-scale droplets ("1-in-2") stabilized by AMB+–mediated reversed Ostwald, with their radius fixed by curvature corrections from $$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$3.
• Large-scale bubbles ("2-in-1") set by the competition between Model B transport and Model A conversion, scaling as $$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$4.
• The resulting steady state, termed "bubbly microphase separation," presents a dual hierarchy of scales, absent in purely conservative or purely reactive models (Li et al., 2021).

6. Regimes, Phase Diagram, and Special Limits

The generalized AMB+ parameter space can be summarized as follows:

Model Limit	Dominant Terms	Regime / Bifurcation
$$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$5	Passive (Model B)	Lifshitz–Slyozov, bulk coarsening
$$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$6 only	Rotation-free active	Shifts coexistence, modifies binodal
$$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$7 only	Rotational active	Interfacial flows, microscale PS
$$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$8	Reaction, nonconserving	Preferred pattern $$\mu_E = \delta F / \delta\phi = -\phi + \phi^3 - \nabla^2\phi$$9, roll states
$J_\lambda = -\lambda\nabla|\nabla\phi|^2$0	Quadratic, off-critical	Tunable super/subcritical transition
Both conservative + reaction (AB+)	All above	Bubbly hierarchical microphase

Special cases include:

• $J_\lambda = -\lambda\nabla|\nabla\phi|^2$1, only supercritical transitions possible (original AMB+ with reaction).
• $J_\lambda = -\lambda\nabla|\nabla\phi|^2$2: recovers Cahn–Hilliard with reaction.
• $J_\lambda = -\lambda\nabla|\nabla\phi|^2$3: returns to purely conserved AMB+.
• For sufficiently slow reactions ($J_\lambda = -\lambda\nabla|\nabla\phi|^2$4), hierarchical microphases develop with small droplets embedded in large bubbles.

7. Significance in Active Matter and Outlook

Generalized Active Model B+ provides a rigorous mesoscopic foundation for understanding phase behavior in active matter, capturing phenomena such as cluster phases, motility-induced phase separation (MIPS), and active microemulsions beyond equilibrium frameworks. Its systematic parameterization allows precise control and prediction of transitions between bulk, microphase, and hierarchically organized steady states using well-characterized deterministic and stochastic PDEs. Current studies focus on nonlinear pattern selection, multistability, and the role of noise. A plausible implication is the wider relevance of AMB+ to biological pattern formation and synthetic active materials, particularly in interpreting nonequilibrium selection mechanisms dictated by interfacial and bulk activity (Yadav et al., 17 Jun 2025, Mondal et al., 12 Jan 2026, Li et al., 2021).