pith. machine review for the scientific record. sign in

arxiv: 2604.26346 · v1 · submitted 2026-04-29 · ❄️ cond-mat.soft · cond-mat.stat-mech· physics.bio-ph

Recognition: unknown

Coexistence of patterned phases in chemically active multicomponent mixtures

Chengjie Luo , Yicheng Qiang , Guido L. A. Kusters , David Zwicker
Authors on Pith no claims yet

Pith reviewed 2026-05-07 12:49 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mechphysics.bio-ph
keywords chemically active mixturespatterned phasesphase coexistenceLyapunov functionalgeneralized Gibbs phase rulemodular combinationphase separationreaction-diffusion
0
0 comments X

The pith

Minimizing a Lyapunov functional in chemically active mixtures yields a generalized Gibbs phase rule for coexisting patterns.

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

The paper studies mixtures where physical interactions drive phase separation while chemical reactions can generate patterns like those in reaction-diffusion systems. It identifies a Lyapunov functional for a class of such reactions, a quantity that decreases over time and reaches a minimum at stationary states. Minimizing this functional produces a generalized Gibbs phase rule that dictates the number of patterns that can coexist. The work also shows that intricate patterns can be formed by combining simpler, independent phases in a modular way. This provides a systematic way to understand and predict the complex stationary structures that arise in these non-equilibrium systems.

Core claim

For chemically active multicomponent mixtures, a Lyapunov functional exists for a class of chemical reactions. Minimization of this functional identifies a generalized Gibbs phase rule governing the number of coexisting patterns and demonstrates that complex patterns can be created by the modular combination of independent phases.

What carries the argument

A Lyapunov functional for the class of chemical reactions, which is minimized to determine the stationary patterned states and enforce the coexistence rules.

Load-bearing premise

That a Lyapunov functional exists and can be minimized for the specific class of chemical reactions and physical interactions considered.

What would settle it

A counterexample where a chemically active mixture reaches a stationary state violating the generalized phase rule, such as more patterns coexisting than predicted, or where the functional does not decrease during evolution.

Figures

Figures reproduced from arXiv: 2604.26346 by Chengjie Luo, David Zwicker, Guido L. A. Kusters, Yicheng Qiang.

Figure 1
Figure 1. Figure 1: FIG. 1. Complex, multi-scale pattern in a mixture of six components view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Schematic of the mapping of two linear reactions view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Reaction rate ratios view at source ↗
read the original abstract

Chemically active mixtures exhibit complex patterns that emerge from the interplay of physical interactions and reactions among components. Individually, these two processes are well-understood: Physical interactions can give rise to phase separation, whereas reactions can form reaction-diffusion patterns. To understand the combination of both processes, we identify a Lyapunov functional for a class of chemical reactions. By minimizing this functional, we identify a generalized Gibbs phase rule that governs the number of coexisting patterns, and we demonstrate that complex patterns can be created by the modular combination of independent phases. Our theory unveils complex stationary patterns in chemically active mixtures and provides a framework for analyzing more complex systems.

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 manuscript claims that for a class of chemically active multicomponent mixtures, a Lyapunov functional can be identified whose minimization yields a generalized Gibbs phase rule governing the number of coexisting patterned phases; complex patterns are further shown to arise via modular combination of independent phases.

Significance. If the Lyapunov functional is rigorously established and the minimization argument holds, the work supplies a variational principle that unifies phase separation with reaction-diffusion dynamics, offering a predictive framework for stationary pattern selection in active mixtures. This would be a notable advance for soft-matter theory, with implications for designing multicomponent systems.

major comments (2)
  1. [Section identifying and deriving the Lyapunov functional] The central claim that minimization of the identified functional produces a generalized phase rule rests on F being a true Lyapunov functional (dF/dt ≤ 0 with equality only at stationary states). The manuscript must explicitly state the functional F, the precise class of reaction terms for which dF/dt ≤ 0 holds identically, and the calculation of its time derivative along the coupled Cahn-Hilliard/reaction-diffusion dynamics (including any assumptions on interaction potentials). Without these steps, the counting argument for the phase rule cannot be verified.
  2. [Derivation of the generalized Gibbs phase rule] In the section deriving the generalized phase rule, the counting of independent constraints at the minima must be shown explicitly, including how the number of coexisting patterns is determined from the degrees of freedom after minimization. Any assumptions about boundedness of F or the form of the reaction kinetics should be stated, as these directly affect whether stationary states coincide with global minima.
minor comments (2)
  1. [Abstract] The abstract is concise but omits any indication of the functional's form or the reaction class; a single sentence or parenthetical example would improve accessibility without lengthening the text.
  2. [Notation and definitions] Notation for the functional and the phase-rule variables should be introduced with a clear table or list of definitions to avoid ambiguity when the counting argument is presented.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive comments that highlight areas where additional explicit detail will strengthen the presentation. We appreciate the recognition of the potential unifying role of the Lyapunov functional. Below we respond point by point to the major comments and indicate the revisions made.

read point-by-point responses

  1. Referee: [Section identifying and deriving the Lyapunov functional] The central claim that minimization of the identified functional produces a generalized phase rule rests on F being a true Lyapunov functional (dF/dt ≤ 0 with equality only at stationary states). The manuscript must explicitly state the functional F, the precise class of reaction terms for which dF/dt ≤ 0 holds identically, and the calculation of its time derivative along the coupled Cahn-Hilliard/reaction-diffusion dynamics (including any assumptions on interaction potentials). Without these steps, the counting argument for the phase rule cannot be verified.

    Authors: We agree that the explicit steps are necessary for independent verification. In the revised manuscript we have inserted a new subsection that (i) writes the functional F in full, (ii) states the precise class of reaction terms (those derivable from a scalar potential linear in the chemical potentials) for which dF/dt ≤ 0 holds identically, and (iii) provides the complete time-derivative calculation along the coupled dynamics, making explicit the assumptions of symmetric pairwise interaction potentials and mass-conserving reactions. These additions directly support the subsequent counting argument. revision: yes

  2. Referee: [Derivation of the generalized Gibbs phase rule] In the section deriving the generalized phase rule, the counting of independent constraints at the minima must be shown explicitly, including how the number of coexisting patterns is determined from the degrees of freedom after minimization. Any assumptions about boundedness of F or the form of the reaction kinetics should be stated, as these directly affect whether stationary states coincide with global minima.

    Authors: We have expanded the phase-rule derivation to include an explicit enumeration of degrees of freedom and constraints. After minimization, the conditions of equal chemical potentials and equal reaction affinities across interfaces reduce the number of independent variables by the expected amount, yielding the generalized rule for the maximum number of coexisting patterned phases. We now state the assumptions on boundedness (F is bounded from below by the stabilizing gradient terms) and on the reaction kinetics (they derive from a potential, guaranteeing that stationary states are global minima of F). revision: yes

Circularity Check

0 steps flagged

No circularity: generalized phase rule follows from explicit minimization of identified Lyapunov functional

full rationale

The paper constructs an explicit Lyapunov functional F for a stated class of reactions, demonstrates dF/dt ≤ 0 along the dynamics, and then applies standard variational minimization to count independent constraints at stationary points, yielding the generalized phase rule. No equation or claim reduces the output count of coexisting patterns to a fitted parameter, a self-referential definition, or a load-bearing self-citation whose validity is presupposed. The modular combination of phases is likewise obtained by direct superposition at the minima of F. The derivation chain is therefore self-contained and does not collapse to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the existence of a Lyapunov functional for the reaction class and the validity of its minimization to produce the phase rule; no free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption A Lyapunov functional exists for the considered class of chemical reactions in multicomponent mixtures
    Invoked as the starting point for minimization to obtain the generalized phase rule.

pith-pipeline@v0.9.0 · 5413 in / 1089 out tokens · 108604 ms · 2026-05-07T12:49:00.534782+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

48 extracted references · 5 canonical work pages · 1 internal anchor

  1. [1]

    C. A. Weber, D. Zwicker, F. J¨ ulicher, and C. F. Lee, Reports on Progress in Physics82, 064601 (2019)

    2019

  2. [2]

    M. C. Cross and P. C. Hohenberg, Reviews of modern physics 65, 851 (1993)

    1993

  3. [3]

    J. W. Gibbs, Trans. Conn. Acad. Arts Sci.3, 1 (1876)

  4. [4]

    S. Mao, M. S. Chakraverti-Wuerthwein, H. Gaudio, and A. Koˇsmrlj, Physical review letters125, 218003 (2020)

    2020

  5. [5]

    Qiang, C

    Y. Qiang, C. Luo, and D. Zwicker, arXiv preprint arXiv:2405.01138 (2024)

    work page arXiv 2024

  6. [6]

    S. C. Glotzer, D. Stauffer, and N. Jan, Physical review letters72, 4109 (1994)

    1994

  7. [7]

    J. J. Christensen, K. Elder, and H. C. Fogedby, Phys. Rev. E54, R2212 (1996)

    1996

  8. [8]

    Zwicker, A

    D. Zwicker, A. A. Hyman, and F. J¨ ulicher, Phys. Rev. E92, 012317 (2015)

    2015

  9. [9]

    Zwicker, Current Opinion in Colloid & Interface Science61, 101606 (2022)

    D. Zwicker, Current Opinion in Colloid & Interface Science61, 101606 (2022)

    2022

  10. [10]

    Bauermann, G

    J. Bauermann, G. Bartolucci, J. Boekhoven, F. J¨ ulicher, and C. A. Weber, Phys. Rev. X15, 041027 (2025)

    2025

  11. [11]

    Liu and N

    F. Liu and N. Goldenfeld, Phys. Rev. A39, 4805 (1989)

    1989

  12. [12]

    Ohta and K

    T. Ohta and K. Kawasaki, Macromolecules19, 2621 (1986)

    1986

  13. [13]

    C. B. Muratov, Phys. Rev. E66, 066108 (2002)

    2002

  14. [14]

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

    work page arXiv 2025

  15. [15]

    Turing, B Jack Copeland519(1952)

    A. Turing, B Jack Copeland519(1952)

    1952

  16. [16]

    Carati and R

    D. Carati and R. Lefever, Phys. Rev. E56, 3127 (1997)

    1997

  17. [17]

    Aslyamov, F

    T. Aslyamov, F. Avanzini, E. Fodor, and M. Esposito, Phys. Rev. Lett.131, 138301 (2023)

    2023

  18. [18]

    Avanzini, T

    F. Avanzini, T. Aslyamov, ´E. Fodor, and M. Esposito, J. Chem. Phys.161, 174108 (2024)

    2024

  19. [19]

    Onsager, Phys

    L. Onsager, Phys. Rev.37, 405 (1931)

    1931

  20. [20]

    J¨ ulicher, S

    F. J¨ ulicher, S. W. Grill, and G. Salbreux, Reports on Progress in Physics81, 076601 (2018)

    2018

  21. [21]

    Zwicker, O

    D. Zwicker, O. W. Paulin, and C. ter Burg, Rep. Prog. Phys.88, 116601 (2025)

    2025

  22. [22]

    Kirschbaum and D

    J. Kirschbaum and D. Zwicker, Journal of The Royal Society Interface18(2021)

    2021

  23. [23]

    Zwicker, Curr

    D. Zwicker, Curr. Opin. Colloid Interface Sci.61, 101606 (2022)

    2022

  24. [24]

    L. Qiao, P. Gispert, L. S. Stelzl, and F. Schmid, Active trans- port as a mechanism of microphase selection in biomolecular condensates (2026), arXiv:2604.08316 [physics.bio-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  25. [25]

    A. M. Rumyantsev and A. Johner, ACS Macro Letters14, 472 (2025)

    2025

  26. [26]

    P. I. Flory, J. Chem. Phys.10, 51 (1942)

    1942

  27. [27]

    M. L. Huggins, J. Chem. Phys.9, 440 (1941)

    1941

  28. [28]

    C. Luo, N. Hess, D. Aierken, Y. Qiang, J. A. Joseph, and D. Zwicker, ACS Macro Letters14, 1484 (2025)

    2025

  29. [29]

    J. Xu, C. Wang, A.-C. Shi, and P. Zhang, Communications in Computational Physics21, 1 (2017)

    2017

  30. [30]

    R. R. Netz, D. Andelman, and M. Schick, Physical review letters 79, 1058 (1997)

    1997

  31. [31]

    Weith, A

    V. Weith, A. Krekhov, and W. Zimmermann, The Journal of Chemical Physics139(2013)

    2013

  32. [32]

    Luo and D

    C. Luo and D. Zwicker, Phys. Rev. E108, 034206 (2023)

    2023

  33. [33]

    Archer, A

    A. Archer, A. M. Rucklidge, and E. Knobloch, Physical Review E92, 012324 (2015)

    2015

  34. [34]

    Archer, A

    A. Archer, A. Rucklidge, and E. Knobloch, Physical review letters111, 165501 (2013)

    2013

  35. [35]

    C. Duan, M. Zhao, Y. Qiang, L. Chen, W. Li, F. Qiu, and A.-C. Shi, Macromolecules51, 7713 (2018)

    2018

  36. [36]

    W. M. Jacobs, Journal of Chemical Theory and Computation 19, 3429 (2023)

    2023

  37. [37]

    Rossetto, M

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

    work page arXiv 2025

  38. [38]

    Greve, G

    D. Greve, G. Lovato, T. Frohoff-H¨ ulsmann, and U. Thiele, Phys- ical Review Letters134, 018303 (2025)

    2025

  39. [39]

    Agrawal, Y

    A. Agrawal, Y. N. Fang, S. Rizvi, N. F. Noor Azman, M. V. Tirrell, A. Karim, and A. E. Neitzel, Macromolecules58, 7776 (2025)

    2025

  40. [40]

    S. Y. Li, M. Kardar, Z. Feng, and W. Taylor, Phys. Rev. E113, 044306 (2026)

    2026

  41. [41]

    1 2 X i κi|∂xϕi(β, x)|2 − X s 1 2 |∂xψs(β, x)|2 dx =−J β 2 Lβ 1 Lβ Z Lβ 0

    A. Al-Hiyasat, D. W. Swartz, J. Gore, and M. Kardar, Spa- tiotemporal noise stabilizes unbounded diversity in strongly- competitive communities (2026), arXiv:2602.13423 [q-bio.PE]. 6 Appendix A: Dynamics of multicomponent mixtures with multiple long-range interactions We consider ad-dimensional incompressible, isothermal mixture composed ofNcomponents wit...

    work page arXiv 2026

  42. [42]

    Two components Specializing the mobility matrix to two componentsAandB, we obtain Λ =D 0 mA 0 0m B − D0 mA +m B m2 A mAmB mAmB m2 B = D0mAmB mA +m B 1−1 −1 1 .(F3) We assume that the only flavor of effective charges is given byqA andq B =−q Ar1. The charge matrixQ ij ≡P s q(s) i q(s) j then becomes Q= X s q(s) i q(s) j =q 2 A 1−r 1 −r1 r2 1 .(F4) The reac...

  43. [43]

    Four components In the case of four components, we find Λ =D 0   mA 0 0 0 0m B 0 0 0 0m C 0 0 0 0m D   − D0 mA +m B +m C +m D   m2 A mAmB mAmC mAmD mAmB m2 B mBmC mBmD mAmC mBmC m2 C mCmD mAmD mBmD mCmD m2 D   =− D0 mA +m B +m C +m D ×   −mA(mB +m C +m D)m AmB mAmC mAmD mAmB −mB(mA +m C +m D)m BmC mBmD mAmC mBmC −mC(mA +m B +m D)m CmD mAmD ...

  44. [44]

    With the choicem 2s+1 =r sm2s+2 for each pairs= 1,

    Generalization This construction extends directly to2nspecies comprisingnindependent reaction pairs. With the choicem 2s+1 =r sm2s+2 for each pairs= 1, . . . , n, the reaction matrixK ij is block-diagonal withn2×2blocks. Each block has the form K(s) =D 0m2s+1q2 2s+1 −1r s 1−r s ,(F14) where the component2s+ 1carries the positive charge of flavors. Thus, a...

  45. [45]

    For example, for a system system comprising four components with reactionA ⇀↽ BandA ⇀↽ C, still at most two phases can form

    Components involved in multiple reactions If a component is involved in more than one reaction, phases with more complicated patterns can form, but the generalized Gibbs phase rule still applies. For example, for a system system comprising four components with reactionA ⇀↽ BandA ⇀↽ C, still at most two phases can form. If we setk AB =k BA =k AC =k CA = 0....

  46. [46]

    2b Associated with the phase diagram shown in Fig

    Amplitude of patterns in Fig. 2b Associated with the phase diagram shown in Fig. 2b of the main text, we in Fig. S2 show the amplitude of the volume fraction ϕA (ϕC) in the patterned phaseLAM AB (LAMCD) as a function ofχandϕ C. The fact that the amplitude goes to zero at the spinodal line indicates that there is a continuous transition from homogeneous ph...

  47. [47]

    Here we show the coexistence of a patterned AB-rich phase and a C-rich homogeneous phase in Fig

    Coexistence of a patterned phase and homogeneous for three components, one reactionN=3,S=1 For three species A, B, and C with one reactionA ⇀↽ B, it is easy to imagine the coexistence of two homogeneous phases: one enriched in A and B, and the other enriched in C. Here we show the coexistence of a patterned AB-rich phase and a C-rich homogeneous phase in Fig. S3

  48. [48]

    To improve numerical stability, we add a solvent componentSwith zero interaction with other components and zero charge to the6-component system mentioned in the main text

    Patterns in a system with seven components We here show additional figures for some simulations of systems with seven components. To improve numerical stability, we add a solvent componentSwith zero interaction with other components and zero charge to the6-component system mentioned in the main text. The solvent component is passive and does not participa...