pith. sign in

arxiv: 2605.15903 · v1 · pith:K6TFAIIOnew · submitted 2026-05-15 · ❄️ cond-mat.soft · cond-mat.stat-mech

Active Model B^- from Mass-Conserving Reaction-Diffusion Systems

Davide Toffenetti , Beatrice Nettuno , Henrik Weyer , Erwin Frey This is my paper

Pith reviewed 2026-05-19 18:58 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mech
keywords mass-conserving reaction-diffusionactive model Bmicrophase separationfinite-wavelength instabilityinterfacial coefficientpattern formationconserved active matter
0
0 comments X

The pith

The late-time dynamics of a minimal three-component mass-conserving reaction-diffusion system reduce to Active Model B^- with negative high-density interfacial coefficient.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that a simple three-component system with conserved mass through reactions and diffusion evolves at late times according to an effective single-component active field theory called Active Model B^-. In this theory, the coefficient that controls the cost of interfaces depends on density and becomes negative when density is high. This change causes the system to form patterns with a preferred finite size rather than letting domains grow forever. The effective theory keeps the chemical potential as a state function thanks to the conservation, but does not give a pressure equation of state. This provides a mechanism for stable microphase separation in mass-conserving active systems.

Core claim

We show that the late-time dynamics of a minimal three-component mass-conserving reaction-diffusion system reduce to a scalar active field theory, Active Model B^-, in which a density-dependent interfacial coefficient κ(φ) turns negative at high density. This drives a finite-wavelength instability and stabilises microphase-separated patterns, in contrast to the unbounded coarsening of two-component mass-conserving systems. Unlike Active Model B^+, AMB^- retains a chemical potential that remains a state function, inherited from the underlying conservation law, but admits no equation of state for the pressure.

What carries the argument

The asymptotic reduction of the three-component reaction-diffusion equations to the Active Model B^- equation featuring a sign-changing density-dependent interfacial coefficient κ(φ).

If this is right

  • Microphase-separated patterns with a characteristic wavelength are stabilized instead of unbounded coarsening.
  • The chemical potential remains a state function due to the underlying mass conservation.
  • No equation of state for the pressure exists in this effective description.
  • The finite-wavelength instability arises specifically from the negative interfacial coefficient at high densities.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This approach may extend to other minimal mass-conserving reaction networks to derive effective active models.
  • Pattern formation in such systems could be controlled by adjusting the reaction rates to tune the density threshold for negative κ.
  • Direct comparison between full multi-component simulations and the reduced AMB^- model at late times would validate the mapping.

Load-bearing premise

The late-time reduction to the specific form of Active Model B^- with the sign change in the interfacial coefficient holds for the minimal three-component reaction network.

What would settle it

A simulation of the three-component system that exhibits continued unbounded coarsening at very late times without stabilizing to finite patterns would contradict the reduction to AMB^-.

Figures

Figures reproduced from arXiv: 2605.15903 by Beatrice Nettuno, Davide Toffenetti, Erwin Frey, Henrik Weyer.

Figure 1
Figure 1. Figure 1: (a) Reaction scheme of the three-component model [Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Phase diagram of phenomenological AMB−, obtained from finite-element simulations [34] in the (ϕ, κ ¯ 1) plane. Increasing κ1 at fixed mean density traverses a se￾quence of morphologies from the PFC limit (κ1 ≲ 0, periodic patterns) to the Model B limit (κ1 ≳ 1, macrophase separa￾tion), passing through coexistence of large and small droplets, microphase-separated droplets, and stripes/foams in the in￾termed… view at source ↗
read the original abstract

We show that the late-time dynamics of a minimal three-component mass-conserving reaction--diffusion system reduce to a scalar active field theory, Active Model B$^-$ (AMB$^-$), in which a density-dependent interfacial coefficient $\kappa(\phi)$ turns negative at high density. This drives a finite-wavelength instability and stabilises microphase-separated patterns, in contrast to the unbounded coarsening of two-component mass-conserving systems. Unlike Active Model B$^+$, AMB$^-$ retains a chemical potential that remains a state function, inherited from the underlying conservation law, but admits no equation of state for the pressure.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper claims that the late-time dynamics of a minimal three-component mass-conserving reaction-diffusion system reduce to the scalar active field theory Active Model B^-, in which a density-dependent interfacial coefficient κ(φ) turns negative at high density. This drives a finite-wavelength instability that stabilizes microphase-separated patterns, in contrast to the unbounded coarsening seen in two-component mass-conserving systems. The chemical potential remains a state function inherited from the underlying conservation law, but there is no equation of state for the pressure.

Significance. If the reduction is valid, the work supplies a microscopic derivation of AMB^- from an explicit mass-conserving reaction network, explaining the origin of the sign change in κ(φ) and the resulting stabilization of finite-wavelength patterns. This bridges reaction-diffusion models with active field theories and provides a route to understand microphase separation in systems where mass conservation is strictly enforced.

major comments (2)
  1. [§3.2] The adiabatic elimination of the two auxiliary fields to obtain the effective AMB^- dynamics for φ assumes a clear separation of timescales. No explicit computation of the eigenvalues of the reaction Jacobian (or equivalent linear stability analysis around the homogeneous state) is supplied to confirm that the fast modes remain stable and decoupled when κ(φ) < 0 drives the finite-k instability.
  2. [§4] The central claim that the late-time dynamics are captured by AMB^- lacks direct numerical validation: there are no quantitative comparisons (e.g., dispersion relations or pattern statistics) between simulations of the full three-component RD system and the reduced AMB^- equation in the regime where the microphase instability occurs.
minor comments (2)
  1. [Abstract] The abstract would be strengthened by naming the specific three-component reaction network and the parameter values at which κ(φ) changes sign.
  2. [§2] Notation for the auxiliary concentrations and the reaction rates could be introduced more explicitly in the model definition to aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment in turn below and have revised the manuscript to incorporate additional analysis and comparisons where appropriate.

read point-by-point responses

  1. Referee: [§3.2] The adiabatic elimination of the two auxiliary fields to obtain the effective AMB^- dynamics for φ assumes a clear separation of timescales. No explicit computation of the eigenvalues of the reaction Jacobian (or equivalent linear stability analysis around the homogeneous state) is supplied to confirm that the fast modes remain stable and decoupled when κ(φ) < 0 drives the finite-k instability.

    Authors: We agree that an explicit verification of the timescale separation strengthens the derivation. In the revised manuscript we have added a linear stability analysis of the homogeneous state for the full three-component system. We compute the eigenvalues of the reaction Jacobian (and the full reaction-diffusion operator) and show that the two fast modes retain negative real parts in the parameter regime where the slow mode develops a finite-wavenumber instability due to negative κ(φ). This analysis is now included in §3.2 together with the relevant parameter ranges, confirming that the adiabatic elimination remains valid. revision: yes

  2. Referee: [§4] The central claim that the late-time dynamics are captured by AMB^- lacks direct numerical validation: there are no quantitative comparisons (e.g., dispersion relations or pattern statistics) between simulations of the full three-component RD system and the reduced AMB^- equation in the regime where the microphase instability occurs.

    Authors: We acknowledge the value of direct numerical validation. We have performed additional simulations of both the original three-component mass-conserving reaction-diffusion system and the reduced Active Model B^- equation. In the revised manuscript we include quantitative comparisons of the dispersion relations extracted from linearised dynamics around the homogeneous state, as well as statistics of the emerging patterns (selected wavelength and modulation amplitude). These results are presented in a new figure in §4 and demonstrate close quantitative agreement between the full system and the effective theory throughout the microphase-separation regime. revision: yes

Circularity Check

0 steps flagged

No circularity: reduction derived from microscopic equations without self-referential inputs

full rationale

The paper derives the reduction of late-time dynamics from a three-component mass-conserving reaction-diffusion system to the scalar Active Model B^- theory, including the density-dependent interfacial coefficient κ(φ) that changes sign. This mapping is obtained by systematic elimination of auxiliary fields from the underlying conservation laws and reaction terms, without fitting the target AMB^- form to data or invoking self-citations as load-bearing justification for the central result. The derivation remains self-contained against the microscopic equations, with the sign change in κ(φ) emerging from the structure of the minimal reaction network rather than being presupposed.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The derivation rests on standard mass conservation and diffusive dynamics plus the assumption that higher-order modes can be adiabatically eliminated at late times; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption Total mass of each component is strictly conserved by the reaction terms.
    Stated in abstract as mass-conserving reaction-diffusion system.
  • domain assumption Late-time dynamics are captured by a single scalar density field after elimination of fast modes.
    Central to the reduction claim.

pith-pipeline@v0.9.0 · 5636 in / 1367 out tokens · 53311 ms · 2026-05-19T18:58:35.359533+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

46 extracted references · 46 canonical work pages

  1. [1]

    Frey and H

    E. Frey and H. Weyer, Annu. Rev. Biophys.55, 10.1146/annurev-biophys-030822-031638 (2026), (in press)

    work page doi:10.1146/annurev-biophys-030822-031638 2026

  2. [2]

    A. B. Goryachev and M. Leda, Molecular Biology of the Cell28, 370 (2017)

    work page 2017

  3. [3]

    Edelstein-Keshet, W

    L. Edelstein-Keshet, W. R. Holmes, M. Zajac, and M. Dutot, Philosophical Transactions of the Royal So- ciety B368, 20130003 (2013)

    work page 2013

  4. [4]

    Lutkenhaus, Annu

    J. Lutkenhaus, Annu. Rev. Biochem.76, 539 (2007)

    work page 2007

  5. [5]

    Loose, E

    M. Loose, E. Fischer-Friedrich, J. Ries, K. Kruse, and P. Schwille, Science320, 789 (2008)

    work page 2008

  6. [6]

    A. G. Vecchiarelli, M. Li, M. Mizuuchi, L. C. Hwang, Y. Seol, K. C. Neuman, and K. Mizuuchi, Proceedings of the National Academy of Sciences113, E1479 (2016)

    work page 2016

  7. [7]

    Glock, B

    P. Glock, B. Ramm, T. Heermann, S. Kretschmer, J. Schweizer, J. M¨ ucksch, G. Alag¨ oz, and P. Schwille, ACS Synthetic Biology8, 148 (2018)

    work page 2018

  8. [8]

    Otsuji, S

    M. Otsuji, S. Ishihara, C. Co, K. Kaibuchi, A. Mochizuki, and S. Kuroda, PLoS Comput. Biol.3, e108 (2007)

    work page 2007

  9. [9]

    A. B. Goryachev and A. V. Pokhilko, FEBS Letters582, 1437 (2008)

    work page 2008

  10. [10]

    Y. Mori, A. Jilkine, and L. Edelstein-Keshet, Biophys. J. 94, 3684 (2008)

    work page 2008

  11. [11]

    Halatek and E

    J. Halatek and E. Frey, Nat. Phys.14, 507 (2018)

    work page 2018

  12. [12]

    Brauns, J

    F. Brauns, J. Halatek, and E. Frey, Phys. Rev. X10, 041036 (2020)

    work page 2020

  13. [13]

    Ishihara, M

    S. Ishihara, M. Otsuji, and A. Mochizuki, Physical Re- view E75, 015203(R) (2007)

    work page 2007

  14. [14]

    Tateno and S

    M. Tateno and S. Ishihara, Physical Review Research3, 023198 (2021)

    work page 2021

  15. [15]

    Brauns, H

    F. Brauns, H. Weyer, J. Halatek, J. Yoon, and E. Frey, Phys. Rev. Lett.126, 104101 (2021)

    work page 2021

  16. [16]

    Weyer, F

    H. Weyer, F. Brauns, and E. Frey, Phys. Rev. E108, 064202 (2023)

    work page 2023

  17. [17]

    Weyer, C.-Y

    H. Weyer, C.-Y. Leung, and E. Frey (2026), manuscript in preparation

    work page 2026

  18. [18]

    Kolokolnikov, T

    T. Kolokolnikov, T. Erneux, and J. Wei, Physica D: Non- linear Phenomena214, 63 (2006)

    work page 2006

  19. [19]

    Kolokolnikov, M

    T. Kolokolnikov, M. J. Ward, and J. Wei, Physica D: Nonlinear Phenomena236, 104 (2007)

    work page 2007

  20. [20]

    Jacobs, J

    B. Jacobs, J. Molenaar, and E. E. Deinum, PLoS ONE 14, e0213188 (2019)

    work page 2019

  21. [21]

    Chiou, K

    J.-G. Chiou, K. D. Moran, and D. J. Lew, eLife10, e58768 (2021)

    work page 2021

  22. [22]

    M. C. Wigbers, T. H. Tan, F. Brauns, J. Liu, S. Z. 6 Swartz, E. Frey, and N. Fakhri, Nat. Phys.17, 578 (2021)

    work page 2021

  23. [23]

    M. E. Cates and C. Nardini, Rep. Prog. Phys.88, 056601 (2025)

    work page 2025

  24. [24]

    Wittkowski, A

    R. Wittkowski, A. Tiribocchi, J. Stenhammar, R. J. Allen, D. Marenduzzo, and M. E. Cates, Nat. Commun. 5, 4351 (2014)

    work page 2014

  25. [25]

    S. C. Glotzer, E. A. Di Marzio, and M. Muthukumar, Physical Review Letters74, 2034 (1995)

    work page 2034

  26. [26]

    Zwicker, A

    D. Zwicker, A. A. Hyman, and F. J¨ ulicher, Physical Re- view E92, 012317 (2015)

    work page 2015

  27. [27]

    Tjhung, C

    E. Tjhung, C. Nardini, and M. E. Cates, Phys. Rev. X 8, 031080 (2018)

    work page 2018

  28. [28]

    Stenhammar, A

    J. Stenhammar, A. Tiribocchi, R. J. Allen, D. Maren- duzzo, and M. E. Cates, Physical Review Letters111, 145702 (2013)

    work page 2013

  29. [29]

    Nardini, ´E

    C. Nardini, ´E. Fodor, E. Tjhung, F. van Wijland, J. Tailleur, and M. E. Cates, Phys. Rev. X7, 021007 (2017)

    work page 2017

  30. [30]

    Caballero, C

    F. Caballero, C. Nardini, and M. E. Cates, J. Stat. Mech. 2018, 123208 (2018)

    work page 2018

  31. [31]

    J. D. Murray,Mathematical Biology II: Spatial Mod- els and Biomedical Applications, 3rd ed. (Springer, New York, 2003)

    work page 2003

  32. [32]

    Gierer and H

    A. Gierer and H. Meinhardt, Kybernetik12, 30 (1972)

    work page 1972

  33. [33]

    Meinhardt, inCurrent Topics in Developmental Biol- ogy, Vol

    H. Meinhardt, inCurrent Topics in Developmental Biol- ogy, Vol. 81 (Academic Press, 2008) pp. 1–63

    work page 2008

  34. [34]

    (2026), see Supplemental Material at [URL] for a discus- sion of the slow-reactivation limitν→0 and for simula- tion details

    work page 2026

  35. [35]

    J. F. Robinson, T. Machon, and T. Speck, Physical Re- view E 10.1103/1bdc-1bjb (2025)

    work page doi:10.1103/1bdc-1bjb 2025

  36. [36]

    K. R. Elder, M. Katakowski, M. Haataja, and M. Grant, Phys. Rev. Lett.88, 245701 (2002)

    work page 2002

  37. [37]

    Weyer, T

    H. Weyer, T. A. Roth, and E. Frey, Nat. Phys.22, 94 (2026)

    work page 2026

  38. [38]

    Zhou and E

    D. Zhou and E. Frey (2026), in preparation

    work page 2026

  39. [39]

    P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys.49, 435 (1977)

    work page 1977

  40. [40]

    Burekovi´ c, F

    S. Burekovi´ c, F. D. Luca, and M. E. Cates, arXiv preprint arXiv:2601.16539 10.48550/arXiv.2601.16539 (2026), arXiv:2601.16539 [cond-mat.soft]

    work page doi:10.48550/arxiv.2601.16539 2026

  41. [41]

    COMSOL AB, Comsol multiphysics (2024),https:// www.comsol.com

    work page 2024

  42. [42]

    Active Model B − from Mass-Conserving Reaction-Diffusion Systems

    D. Toffenetti, B. Nettuno, H. Weyer, and E. Frey, Simu- lation files for “Active Model B − from Mass-Conserving Reaction-Diffusion Systems” (2026)

    work page 2026

  43. [43]

    Rossetto, M

    R. Rossetto, M. Ernst, and D. Zwicker, arXiv preprint 10.48550/arXiv.2511.03619 (2025)

    work page doi:10.48550/arxiv.2511.03619 2025

  44. [44]

    F. C. Thewes, Y. Qiang, O. W. Paulin, and D. Zwicker, arXiv preprint 10.48550/arXiv.2511.05214 (2025)

    work page doi:10.48550/arxiv.2511.05214 2025

  45. [45]

    M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys.65, 851 (1993)

    work page 1993

  46. [46]

    Thiele, A

    U. Thiele, A. J. Archer, and M. J. Robbins, Physical Review E87, 042915 (2013). End Matter Reduction of the three-component model toAMB − We give the algebraic steps underlying the reduction of Eq. (1) to Eq. (4) in the regimeν −1, dm ≪1. Mass conservation gives the continuity equation ∂tϕ=D c∇2η,(7) so the task is to expressηas a functional ofϕ. We use (...

    work page 2013