Pith. sign in

REVIEW 3 major objections 5 minor 30 references

The physics and mathematics of living and dying matter

T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Finite-fuel active matter settles into passive equilibrium; the active phase is a transient.

desk verdict A clean finite-fuel active matter template whose equilibrium endpoint is built in by construction; the useful physics is there, but the life-as-transient claim owes more to the setup than to emergent dynamics. read the letter →

arxiv 2507.14294 v1 pith:5X7TUXMY submitted 2025-07-18 cond-mat.soft physics.bio-ph

Tanniemola B. Liverpool , Kristian K. Müller-Nedebock , Xichen Chao , Chang Yuan This is my paper · ORCID
classification cond-mat.softphysics.bio-ph
keywords activematterfinitefuelLangevindynamicsVicsekmodelflockingtransitionnon-equilibriumtransientthermodynamicconsistencydepletion
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper introduces active matter models that explicitly track stored energy—the 'fuel' that drives self-propulsion—instead of assuming an infinite energy reservoir. The authors aim to show that such finite-fuel systems are thermodynamically consistent and that their true steady state is a passive equilibrium; the ordered active phase is a possibly long-lived transient. Using a two-dimensional flocking model related to the Vicsek model, they show that a disordered passive system becomes an aligned flock while fuel lasts, then loses order and relaxes to equilibrium as fuel runs out. If correct, the result gives a general picture: activity in a closed finite-energy system is an episode, not a permanent state.

What carries the argument

The central object is the non-symmetric mobility matrix $G_i$ that couples fuel consumption to self-propulsion in the Langevin equations. Choosing the noise covariance as $A_i = 2\Theta\Gamma_i$, where $\Gamma_i$ is the symmetric part of $G_i$, makes the Fokker-Planck equation have a steady state $P \propto e^{-H/\Theta}$ with $H = U + \Phi$, absorbing the fuel potential into an effective Hamiltonian. This choice is what gives the model thermodynamic consistency and gives the finite-fuel system a passive equilibrium endpoint.

What would settle it

Simulate the same model with a fuel consumption rate that increases with local alignment; if the system still reaches the passive Gibbs distribution after the fuel is gone and shows no lingering order, the framework survives, whereas persisting order or a steady state that differs from $e^{-(U+\Phi)/\Theta}$ would falsify it.

Watch

Extended reading notes

Core claim

The authors claim that keeping track of stored energy turns the non-equilibrium active state into a transient: the full system's steady state is the equilibrium Gibbs-Boltzmann distribution $P_{\mathrm{eq}} = Z^{-1} e^{-H/\Theta}$ with $H = U + \Phi$, where $U$ is the mechanical alignment potential and $\Phi$ is the fuel potential. In their two-dimensional Vicsek-type flocking model, simulations show the order parameter rising from zero to a plateau while fuel is consumed, then decaying to zero when the fuel runs out, returning the system to a disordered passive equilibrium. The central discovery is that the active, ordered phase is a possibly long-lived transient between two equilibrium states.

Load-bearing premise

The framework assumes fuel is consumed at a rate that depends only on how much fuel remains, with particle motion and alignment never feeding back into the rate of fuel use; if that coupling is reversed, the clean return to passive equilibrium is no longer automatic.

Editorial extensions

If this is right

  • Any active model expressible as Langevin equations with explicit finite stored fuel should show the same disorder-order-disorder lifecycle as its fuel is consumed.
  • In a closed system, the ordered active state cannot be a true steady state; it lasts only while fuel is being converted into mechanical work.
  • The effective Hamiltonian $H = U + \Phi$ provides a bookkeeping tool for tracking mechanical potential, fuel potential, and total energy through the active phase and its decay.
  • The fuel potential $\phi(n)$ controls the duration and shape of the active transient, so different fuels yield different flocking lifetimes without changing the equilibrium endpoint.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper leaves open the reverse coupling in which mechanical motion changes fuel consumption; if that coupling is included, the monotonic fuel depletion and exact Gibbs-Boltzmann endpoint may fail, so the framework's boundary is a concrete next test.
  • Translated to the lab, a suspension of synthetic swimmers powered by a finite chemical reservoir should show an order parameter that rises, plateaus, and decays to zero on a timescale set by the fuel potential; measuring that curve would test the transient picture directly.
  • Applied to organisms that replenish their energy stores, the framework would need a term for fuel intake, which would likely replace the passive equilibrium endpoint with a sustained non-equilibrium cycle rather than a transient.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript introduces a class of overdamped Langevin active-matter models with an explicit fuel variable n_i per particle. The mechanical degrees of freedom (positions and orientations) evolve under an aligning potential U, while the active self-propulsion term is proportional to dn_i/dt. Fuel dynamics is one-way: Eq. (3) is independent of positions and orientations. With the multiplicative noise choice A_i = 2ΘΓ_i, the authors write a Fokker–Planck equation and claim a Gibbs–Boltzmann steady state e^{-H/Θ}, where H=U+Φ. Simulations of a two-dimensional Vicsek-like flocking model show a disorder→order→disorder sequence as the fuel is consumed. The paper concludes that equilibrium (death) is the true steady state and that the active non-equilibrium state (life) is a transient between two equilibrium states, and claims the framework applies to any active model described by Langevin equations.

Significance. If correct, the paper would give a finite-fuel, thermodynamically consistent extension of standard active models and a concrete way to follow the onset and dissipation of activity. The Fokker–Planck steady-state calculation is standard and the simulations illustrate the intended disorder–order–disorder behavior. However, the main conceptual conclusion is substantially built into the model construction: the one-way fuel coupling and the choice of the active force make the passive endpoint a design feature rather than an emergent prediction. The claimed universality is not established. The paper is best viewed as a well-defined minimal model of finite-fuel activity, and its value would be clearer if the claims were narrowed accordingly.

major comments (3)
  1. [Model, Eq. (3)] The one-way coupling is load-bearing. Eq. (3) contains no dependence on r_i or θ_i, so the fuel variables evolve as an autonomous stochastic process and the active term -v0 u_i dn_i/dt in Eq. (1) is a pre-scheduled, decaying external drive. The conclusion that the active state is a transient between two equilibrium states is therefore a direct consequence of the model construction; the text itself says 'by definition' at the point where this is stated. The concluding claim that the framework can be applied to any active model described by Langevin equations is not supported: if fuel consumption depends on the mechanical variables, or if motion replenishes fuel, the monotonic depletion and the clean return to the passive equilibrium of Eq. (8) are not guaranteed. The caveat before Eq. (3) and footnote [21] acknowledge the restriction, but the conclusion does not.
  2. [Eq. (1), model section] The active self-propulsion term is proportional to the total time derivative dn_i/dt, which includes the noise in Eq. (3). Even when n_i reaches the minimum of ϕ at n=0, the noise term in dn_i/dt is generically nonzero, so the active force does not strictly vanish when the fuel runs out. The statement that the active terms 'vanish when the fuel runs out' is only valid if the deterministic consumption rate -∂ϕ/∂n is used instead of the total derivative, or if n is absorbed at 0 with a reflecting boundary. Please clarify this point, as the passivity of the final state is central to the main conclusion.
  3. [Noise choice, Eq. (7)] The equilibrium end-state is put in by hand through the noise choice A=2ΘΓ and the fuel potential's minimum at n=0. This is a legitimate modeling choice, but it means that the identification 'equilibrium (death) is the true steady state' is not a nontrivial prediction of the dynamics. The paper should present these choices as the design mechanism responsible for the transient life/death behavior and should not present that behavior as an emergent property of active matter.
minor comments (5)
  1. [Eq. (10)] The definition is written as v_a(t) ≡ ⟨|v_i(t)|⟩, but the right-hand side is |(1/N)∑ v_i|; the angle brackets and absolute value are in the wrong order. Please correct the notation.
  2. [Figure 3 caption] The caption says 'Rescaled by the initial fuel potential Φ(n0)', but the panels show H, U and Φ; please state explicitly which normalization is used for each curve.
  3. [Figure 4] The main text does not report the number of independent realizations or error bars for the ensemble averages in Fig. 4; please provide this information or state where it appears in the Supplementary Information.
  4. [Simulation parameters] The main text gives ζθ, ζ, ζn, R and J0, but omits the temperature Θ and the self-propulsion speed v0; these should be stated for reproducibility.
  5. [References] Reference [9] and reference [19] are the same Vicsek et al. paper; please cite it only once.

Circularity Check

3 steps flagged · score 8.0 of 10

The headline 'life as transient, death as steady state' is imposed by the model's design requirements and one-way fuel equation, and the paper itself calls the transient 'by definition'.

  1. self definitional [Model setup before Eq. (1); Conclusion]
    "We want to construct models in which (1) the active terms are proportional to the rate of consumption of fuel (stored energy) and vanish when the fuel runs out and (2) when the fuel runs out the system behaves as a passive system with potential energy U that evolves to a state of equilibrium at the same temperature as its environment with a detailed balance condition. ... we identify equilibrium (death) as the true steady state and that the active non-equilibrium state (life) is a possibly long-lived transient between two equilibrium states (death)."

    The conclusion is the construction requirement restated: 'passive equilibrium when fuel runs out' is exactly 'death as true steady state', and 'active terms vanish when the fuel runs out' is exactly 'life as transient'. The endpoint is an input design goal, not a result derived from independent active-matter dynamics.

  2. self definitional [Equations (1)-(3) and text before Eq. (3)]
    "at first we will only consider the case in which the consumption of stored energy affects the mechanical variables, generating the active terms, but not vice-versa [21]. ... dn_i/dt = - (1/zeta_n) partial Phi/partial n_i + eta^n_i ... The active state is than by definition a, possibly long-lived but, transient state between two passive states."

    Eq. (3) contains no positional or orientational variables, so n_i(t) is an autonomous stochastic process; the active force in Eq. (1) is a pre-scheduled function of this independent fuel clock. Fuel necessarily depletes and activity necessarily ceases regardless of the mechanical state. The paper's own 'by definition' sentence concedes that the transient is an imposed property, not an emergent prediction.

1 more flagged steps
  1. self definitional [Equations (6)-(8)]
    "a useful choice is A_i = 2 Theta Gamma_i(X_i) ... So our choice for the noise fluctuations leads us to a simple, explicit expression for the steady-state distribution. ... P_eq({X_i}) = 1/Z e^{-H({X_i})/Theta}."

    The Boltzmann steady state for the full H = U + Phi is engineered by choosing the noise matrix as 2 Theta Gamma. Since Phi is minimized at n = 0, the designed equilibrium is the depleted-fuel passive state. The later observation that the system 'returns to equilibrium' after fuel runs out is the chosen noise rule plus the chosen fuel potential, not an independent consequence of active-matter physics.

full rationale

The paper is transparent that its models are 'by construction' thermodynamically consistent, and the headline conclusion 'death as true steady state, life as transient' is exactly the pair of design requirements stated before Eq. (1): active terms vanish when fuel runs out, and the fuel-depleted system is passive with a detailed-balance equilibrium. The fuel equation (3) is autonomous, so the depletion clock runs independently of the mechanics; the paper itself calls the active state 'by definition' a transient. The noise choice A = 2 Theta Gamma engineers the Boltzmann steady state for H = U + Phi, and since Phi is minimized at n = 0, the engineered equilibrium is the dead state. Thus the central interpretive claim is an input assumption, not a derived prediction. The flocking simulations, the pseudo-steady state, and the density-dependent disorder-order-disorder transition are genuine numerical outputs not encoded in the inputs, and the FP steady-state existence result is standard; the self-citations [27,28] are not load-bearing here. The concluding universality claim for 'any active model which can be described in terms of Langevin equations' outruns the one-way-coupling model, as footnote [21] concedes. Overall score 8: the headline result is forced by definition, although the paper also contains independent simulation content.

Assumptions & free parameters 5 free parameters · 5 assumptions · 1 invented entities

The central claim rests on explicit modeling choices: one-way fuel-mechanics coupling, a positive-definite mobility condition, and the noise construction that forces a Gibbs equilibrium. The fuel potential coefficients and friction/alignment parameters are hand-set, not fitted to experiments. The only invented entity is the fuel degree of freedom, which has no independent empirical handle in this preprint.

free parameters (5)
  • fuel potential coefficients k_i = k1=2, k2=0.0026, k3=3, k4=2
    Chosen so all four fuel potentials run out at similar times; comparisons across potentials depend on this calibration.
  • friction and alignment parameters = zeta_theta=200, zeta_r=8, zeta_n=8, J0=8.75, R=1
    Hand-set simulation parameters used for all runs; the observed flocking window depends on these values.
  • initial fuel amount n0 = finite, n0 > 0; exact value deferred to SI
    Initial condition for all particles; the duration of the active transient depends on it.
  • self-propulsion speed v0 = not stated in main text; constrained by v0^2 < 4 zeta_r zeta_n
    Controls the active force and the positive-definiteness condition for the noise; exact value appears in SI.
  • thermostat temperature Theta = not stated in main text; varied in SI
    Sets the noise magnitude and hence the equilibrium fluctuations; needed to reproduce the reported transitions.
assumptions (5)
  • standard math The Ito discretization and divergence-free conditions d/dX_i * Gamma_i = d/dX_i * Omega_i = 0 are used to derive the Fokker-Planck equation.
    Invoked after Eq. (4)-(8); this is standard stochastic calculus, but the divergence-free condition is stated rather than proved in the main text.
  • domain assumption Noise matrix A_i = 2*Theta*Gamma_i, with Gamma_i the symmetric part of G_i, gives a Gibbs steady state at temperature Theta.
    Eq. (7) and the resulting steady state in Eq. (8). This is the thermodynamic consistency choice and assumes the equilibrium temperature is the environment temperature.
  • domain assumption The matrix Gamma_i is strictly positive definite, requiring v0^2 < 4 zeta_r zeta_n.
    Required to define the noise through the square root of Gamma; if violated, the multiplicative-noise construction breaks down.
  • ad hoc to paper Fuel consumption affects mechanical variables but not vice versa, corresponding to a very anisotropic reaction landscape.
    Stated after Eq. (3); this restriction makes fuel depletion exogenous and guarantees return to a passive equilibrium, so the central life-death arc rests on it.
  • standard math The Fokker-Planck equation has a unique stable steady state given by the Gibbs distribution.
    Attributed to refs [25-28]; the main text does not prove uniqueness or convergence.
invented entities (1)
  • Fuel variable n_i(t), a stored-energy degree of freedom per particle
    purpose: Keeps track of remaining fuel; active self-propulsion is proportional to its consumption rate and vanishes at n=0.
    No experimental calibration or falsifiable prediction tied to n_i is given; it is a modeling device, and the 'dead' equilibrium at n=0 is assumed by construction.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The physics and mathematics of living and dying matter." pith.science (2026). https://pith.science/paper/5X7TUXMY

@misc{pith2026250714294,
  author       = {Pith},
  title        = {Pith review of: The physics and mathematics of living and dying matter},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5X7TUXMY}},
  note         = {Machine review of arXiv:2507.14294}
}
read the original abstract

We introduce and study a class of active matter models in which we keep track of fuel (stored energy) consumption. They are by construction, thermodynamically consistent. Using these models it is possible for us to observe and follow how active behaviour develops and also how it dissipates as the energy runs out. It is also straightforward to define, calculate and keep track of macroscopic thermodynamic quantities.

Figures

Figures reproduced from arXiv: 2507.14294 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

30 extracted references · 23 canonical work pages

  1. [21]

    It is not difficult to generalise the model to do so

    work page

  2. [1]

    J. R. Baylis, J. H. Yeon, M. H. Thomson, A. Kazerooni, X. Wang, A. E. St. John, E. B. Lim, D. Chien, A. Lee, J. Q. Zhang,et al., Science advances1, e1500379 (2015)

    work page 2015

  3. [2]

    Di Leonardo, L

    R. Di Leonardo, L. Angelani, D. Dell’Arciprete, G. Ruocco, V. Iebba, S. Schippa, M. P. Conte, F. Mecarini, F. De Angelis, and E. Di Fabrizio, Pro- ceedings of the National Academy of Sciences107, 9541 (2010)

    work page 2010

  4. [3]

    Krishnamurthy, S

    S. Krishnamurthy, S. Ghosh, D. Chatterji, R. Ganapathy, and A. Sood, Nature Physics12, 1134 (2016)

    work page 2016

  5. [4]

    Stenhammar, R

    J. Stenhammar, R. Wittkowski, D. Marenduzzo, and M. E. Cates, Science advances2, e1501850 (2016)

    work page 2016

  6. [5]

    Frangipane, D

    G. Frangipane, D. Dell’Arciprete, S. Petracchini, C. Maggi, F. Saglimbeni, S. Bianchi, G. Vizsnyiczai, M. L. Bernardini, and R. Di Leonardo, Elife7, e36608 (2018)

    work page 2018

  7. [6]

    J. Arlt, V. A. Martinez, A. Dawson, T. Pilizota, and W. C. Poon, Nature communications9, 768 (2018)

    work page 2018

  8. [7]

    Y. Chen, X. Li, C. Scheibner, V. Vitelli, and G. Huang, Nature communications12, 5935 (2021)

    work page 2021

Show all 30 references
  1. [8]

    Prost and R

    J. Prost and R. Bruinsma, Europhysics Letters (EPL) 33, 321 (1996)

    1996

  2. [10]

    Toner and Y

    J. Toner and Y. Tu, Phys. Rev. Lett.75, 4326 (1995)

    1995

  3. [11]

    Toner, Y

    J. Toner, Y. Tu, and S. Ramaswamy, Annals of Physics 318, 170 (2005), special Issue

    2005

  4. [12]

    J. R. Howse, R. a. L. Jones, A. J. Ryan, T. Gough, R. Vafabakhsh, and R. Golestanian, Phys. Rev. Lett. 99, 8 (2007), arXiv:0706.4406

    2007 arXiv

  5. [13]

    Bricard, J.-B

    A. Bricard, J.-B. Caussin, N. Desreumaux, O. Dauchot, and D. Bartolo, Nature503, 95 (2013)

    2013

  6. [14]

    Julicher, K

    F. Julicher, K. Kruse, J. Prost, and J. Joanny, Phys. Rep.449, 3 (2007)

    2007

  7. [15]

    Ramaswamy, Annu

    S. Ramaswamy, Annu. Rev. Condens. Matter Phys.1, 323 (2010)

    2010

  8. [16]

    M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Rev. Mod. Phys.85, 1143 (2013)

    2013

  9. [17]

    Bechinger, R

    C. Bechinger, R. Di Leonardo, H. L¨ owen, C. Reichhardt, G. Volpe, and G. Volpe, Rev. Mod. Phys.88(2016), 10.1103/RevModPhys.88.045006, arXiv:1602.00081

    2016 arXiv

  10. [18]

    Gompper, R

    G. Gompper, R. G. Winkler, T. Speck, A. Solon, C. Nar- dini, F. Peruani, H. L¨ owen, R. Golestanian, U. B. Kaupp, L. Alvarez, T. Kiørboe, E. Lauga, W. C. K. Poon, A. DeSimone, S. Mui˜ nos-Landin, A. Fischer, N. A. S¨ oker, F. Cichos, R. Kapral, P. Gaspard, M. Ripoll, F. Sagues...

    2020

  11. [19]

    Vicsek, A

    T. Vicsek, A. Czir´ ok, E. Ben-Jacob, I. Cohen, and O. Shochet, Phys. Rev. Lett.75, 1226 (1995)

    1995

  12. [20]

    Chat´ e, F

    H. Chat´ e, F. Ginelli, and F. Raynaud, Phys. Rev. E77, 046113 (2008)

    2008

  13. [22]

    N. G. V. Kampen,Stochastic processes in physics and chemistry(Elsevier, Amsterdam, 2001)

    2001

  14. [23]

    Risken,The Fokker-Planck Equation(Springer: Berlin, 1984)

    H. Risken,The Fokker-Planck Equation(Springer: Berlin, 1984)

    1984

  15. [24]

    For modelling the dynamics of physical problems at equi- librium, one often chooses the discretisation that gives rise to the correct steady state equilibrium probability distribution

  16. [25]

    P. A. Markowich and C. Villani, Mat. Contemp.19, 1 (2000)

    2000

  17. [26]

    Guionnet and B

    A. Guionnet and B. Zegarlinski, S´ eminaire Probab. Strasbg.36, 1 (2002)

    2002

  18. [27]

    T. B. Liverpool, Phys. Rev. E101, 042107 (2020)

    2020

  19. [28]

    Cameron, M

    S. Cameron, M. Mosayebi, R. Bennett, and T. B. Liver- pool, Phys. Rev. E108, 1 (2023), arXiv:2201.10813

    2023 arXiv

  20. [29]

    D. J. Higham, SIAM Rev.43, 525 (2001)

    2001

  21. [30]

    Derrida, J

    B. Derrida, J. Stat. Mech. Theory Exp. (2007), 10.1088/1742-5468/2007/07/P07023, arXiv:0703762 [cond-mat]

    2007 doi

  22. [31]

    Jona-Lasinio, Nonlinear Process

    G. Jona-Lasinio, Nonlinear Process. Geophys.30, 253 (2023)

    2023

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.