REVIEW 4 major objections 5 minor 27 references
Stochastic bubble dynamics in phase-separated scalar active matter
T0 review · 4 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Vapor bubbles in phase-separated active liquids mostly dissolve internally, and their areas obey a Langevin equation with constant negative drift and perimeter-proportional noise, yielding a generalized Gibbs distribution and lifetimes…
desk verdict A useful empirical characterization of MIPS bubbles as a one-variable Langevin process, but the central claim is under-tested: the fits use the same data and the Markov property is never checked. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the effective Langevin equation for bubble area, Eq. (2), written in Itô form. Its multiplicative noise amplitude, $\sqrt{2\Gamma A^{1/4}}$, is proportional to the square root of the bubble perimeter, so the diffusion coefficient $\Gamma\sqrt{A}$ scales with perimeter. From this equation, via the Fokker–Planck equation, the stationary distribution, first-passage time, and lifetime distribution follow. The lattice gas model—square blocks added or removed at boundary sites with rates set by the fitted drift and diffusion and a Boltzmann factor $e^{\epsilon m_i}$ favoring more-coordinated sites—supplies the morphological predictions, reproducing the decline of asphericity with area.
What would settle it
One concrete check is to measure the first and second moments of area increments at several lag times (e.g., 0.02, 0.05, 0.1, and 0.2 $\tau$) and test whether the inferred $\nu$ and $\Gamma$ stay fixed and the transition probabilities obey the Chapman–Kolmogorov equation; if the coefficients change with lag time or the increment distribution is non-Gaussian, the Markovian Itô description fails. A second check is to follow the bubbles that merge with the domain boundary and test whether their area statistics obey the same $\nu$ and $\Gamma$; if not, the restriction to internally dissolving bubbles biases the result.
Extended reading notes
Core claim
The central discovery is that the bubble area $A_t$ is a Markov process described by the Itô Langevin equation $dA_t = -\nu\,dt + \sqrt{2\Gamma A_t^{1/4}}\,dW_t$, with $\nu = 19.3\,\sigma^2/\tau$ and $\Gamma = 348\,\sigma^3/\tau$ fitted from particle simulations. This equation is equivalent to a Fokker–Planck equation whose stationary solution is the generalized Gibbs distribution $p(A) \propto A^{-1/2} \exp(-2\nu\sqrt{A}/\Gamma)$, which matches the simulated area distribution down to the resolution limit. It also yields the mean first-passage time $T(A) = A/\nu$ and an inverse Gaussian lifetime distribution that fit the simulated lifetimes. A companion lattice gas model with perimeter-proportional addition and removal rates reproduces the measured decrease of asphericity with area, supporting the interpretation that boundary fluctuations drive the process. The paper contrasts these results with Active Model B+, which predicts a peaked area distribution not seen in the particle simulations.
Load-bearing premise
The analysis assumes that the bubble area alone is a closed Markov process: the drift and noise depend only on the instantaneous area, with Gaussian white noise, and all other variables are either negligible or slaved to area; it also assumes that discarding the roughly 3% of bubbles that reach the boundary does not bias the fitted drift and diffusion.
Editorial extensions
If this is right
- Once $\nu$ and $\Gamma$ are fitted, the stationary area distribution and the mean lifetime $T(A) = A/\nu$ are predictions with no further free parameters.
- The inverse Gaussian form of the lifetime distribution means that measuring lifetimes in experiments or simulations is enough to extract the drift and noise strength.
- The effective line tension $\lambda = \nu/(\sqrt{\pi}\,\Gamma)$ turns the noise-dominated bubble statistic into a thermodynamic-like free energy, suggesting a consistent interfacial description even away from equilibrium.
- The lattice gas model implies that local boundary events, not global rearrangements, set the shape statistics, connecting area noise to the same perimeter mechanism.
- The general equation $dA_t = -\nu A^{\beta}\,dt + \sqrt{2\Gamma A^{\alpha}}\,dW_t$ with $\alpha = 1/4$, $\beta = 0$ defines a class of sub-demographic noise processes whose statistical properties are interesting in their own right.
Reading between the lines
- Beyond the paper: if the same perimeter-proportional noise holds in three dimensions, the analogous area (volume) exponent should change from $1/4$ to $1/3$, since the diffusion coefficient would scale with surface area, i.e., $V^{2/3}$.
- Beyond the paper: the measured discrepancy with Active Model B+ suggests continuum theories may need nonlocal terms; a testable extension is to compute the area distribution predicted by Active Model B+ in the same parameter regime and compare it quantitatively with Fig. 3(a).
- Beyond the paper: the Markov property could be tested by checking higher-order moments of area increments, for example the ratio of third to second moments, which would reveal non-Gaussian noise and break the Itô Langevin description.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies vapor bubbles inside the liquid domains of phase-separated active Brownian particles using large-scale particle-based simulations. It reports that 97% of bubbles nucleate, grow, and dissolve within the bulk of the liquid domain, and proposes that the bubble area A_t obeys an Itô Langevin equation dA_t = -ν dt + sqrt(2Γ A^{1/4}) dW_t (Eq. 2), with fitted constants ν = 19.3 σ²/τ and Γ = 348 σ³/τ. From this equation the authors derive a stationary area distribution, a linear mean first-passage time T(A) = A/ν, and lifetime distributions, and they compare these against the simulations. They also introduce a lattice gas model for bubble shape and show that the mean asphericity decreases with area. The paper contrasts its findings with Active Model B+ and highlights limitations of continuum theories.
Significance. The paper addresses a genuinely open question in motility-induced phase separation: the stochastic dynamics of vapor bubbles inside dense domains. The proposed effective Langevin equation is simple, physically interpretable (noise proportional to the bubble perimeter), and it makes explicit predictions that are partly confirmed by a large data set of 30,012 simulated bubbles. The connection to a broader class of sub-demographic multiplicative noise processes is appealing and could be of interest beyond this specific system. However, the validation presented is mostly in-sample consistency rather than independent prediction: the drift and diffusion coefficients are fitted from the same simulation increments used to test the stationary distribution, the MFPT is insensitive to the functional form of the diffusion coefficient, and the lifetime distribution requires an additional tuned parameter A_eff. The Markov and Gaussian-white-noise assumptions underlying Eq. (2) are not directly tested. These issues currently limit the strength of the claim that Eq. (2) 'fully characterizes' the bubble statistics.
major comments (4)
- [Eq. (4) and Fig. 3(a)] The stationary distribution p(A) ∝ A^{-1/2} exp(-2ν√A/Γ) is derived from the same fitted values of ν and Γ that were obtained from the simulation increments in Fig. 2(b). The agreement in Fig. 3(a) is therefore a consistency check, not an independent validation of the Langevin form. The same applies to the mean first-passage time T(A) = A/ν shown in Fig. 3(b): for any diffusion coefficient D(A) with D(0)=0, the equation -ν T' + D(A) T'' = -1 admits the solution T(A)=A/ν, so this result tests only the constant drift and not the perimeter-proportional noise. The authors should either explicitly label these as consistency tests or demonstrate prediction by, for example, fitting on one part of the data and predicting another part.
- [Eq. (2) and Fig. 2(b)] The Markov property and the Gaussian white-noise assumption are not tested. The drift and diffusion are estimated from moment fits at a single lag Δt = 0.1τ, and the paper does not check whether the inferred coefficients are independent of the sampling lag, whether the transition density P(A_t | A_0) satisfies the Chapman-Kolmogorov equation, or whether the conditional increment distribution is approximately Gaussian for small Δt. Without such checks, Eq. (2) is a Kramers-Moyal truncation whose validity as a closed stochastic evolution equation for A_t is unsupported. The authors should add at least one direct test of the Markov property, for example by comparing two-lag predictions with direct simulation or by showing that the inferred drift and diffusion do not change when Δt is varied.
- [Eq. (5) and Fig. 3(c,d)] The lifetime distribution in Eq. (5) uses a second fitted parameter A_eff = A0/2, and the population lifetime w(T) in Eq. (6) uses an additional cutoff Amin = 10σ², both chosen for best fit. The inverse-Gaussian curves in Fig. 3(c,d) are therefore not parameter-free predictions of Eq. (2). Moreover, the 'Langevin dynamics' curves in those panels are numerical solutions of the very equation that was fitted, so they cannot independently validate the model against the particle simulations. The authors should clarify the status of A_eff and Amin, ideally deriving them from the model or treating the lifetime comparison as a demonstration that the model can be tuned to match the data.
- [Fig. 2(b) and Fig. 3(a)] The small-area regime is explicitly identified as problematic: the first two drift data points deviate from the constant value, and the stationary distribution deviates for the smallest bubble sizes. Since bubble lifetimes are controlled by the dynamics near A=0, these deviations are not merely cosmetic. The paper should quantify the range of A over which Eq. (2) is validated and discuss whether the cutoff Amin (used for normalization and in Eq. (6)) biases the lifetime statistics. If the model is intended only for intermediate and large areas, the claim that it 'fully characterizes' lifetime statistics should be weakened accordingly.
minor comments (5)
- [Reference [19]] The supplementary material is referenced with a placeholder title and link ('link-provided-after-publication'); since many technical details and additional figures are deferred to [19], the manuscript is not self-contained for review. The authors should provide the supplemental material as part of the submission.
- [Fig. 2(a)] It would be helpful to show a Gaussian fit (or other comparison) overlaid on the empirical distribution of ΔA for a few representative values of A, to support the Gaussian white-noise assumption used in Eq. (2).
- [Conclusion] The text contains a typo: 'motility-induces phase separation' should read 'motility-induced phase separation'. Also, the Péclet number is written as 'P´eclet' with inconsistent accents.
- [Fig. 4(b)] The lattice gas parameter ε = 1.1 is fitted to the asphericity data, but no error bars or sensitivity analysis are shown. A brief statement of how ε was determined and how sensitive the curves are to this parameter would strengthen the morphological comparison.
- [Eq. (8)] The general family dA_t = -νA^β dt + sqrt(2ΓA^α) dW_t is introduced at the end without any further analysis. This is fine as outlook, but the notation 'sub-demographic noise' is used without definition; a short explanation would help readers unfamiliar with demographic noise.
Circularity Check
No significant circularity: the Langevin description is an empirical fit whose derived stationary/lifetime statistics are consistency checks, not independent predictions, and no load-bearing argument reduces to a self-citation.
full rationale
The paper does not claim a first-principles derivation; Eq. (2) is an effective Langevin model fitted to conditional moments of the ABP simulations (Fig. 2b). The stationary distribution Eq. (4) and MFPT T(A)=A/ν are mathematical consequences of that fitted SDE, and comparing them with the same simulation data is an in-sample consistency check rather than an independent prediction. That is a statistical limitation, and the Markov closure is an untested assumption, but it is not circular: the fitted drift and diffusion do not by construction determine the empirical stationary histogram, the mean first-passage times, or the lifetime distributions, and the numerical solution of Eq. (2) uses only parameters fitted to short-time increments to reproduce lifetime statistics. The lattice-gas model uses the same fitted coefficients but is tested against a separate morphological observable (asphericity), so it is not a re-derivation of the input. Self-citations (Refs. [4], [26]) are incidental and not load-bearing. No equation is defined in terms of the quantity it is used to predict, and no prediction is statistically forced by the fitting procedure.
Assumptions & free parameters
free parameters (5)
- ν (drift rate) =
19.3 σ²/τ
- Γ (noise amplitude coefficient) =
348 σ³/τ
- ε (lattice gas interaction) =
1.1
- A_eff (effective area in lifetime approximation) =
A0/2
- A_min (cutoff for normalization/lifetime integral) =
4.3752 σ² and 10 σ²
assumptions (4)
- domain assumption Bubble area evolves as a Markov process with drift and diffusion depending only on instantaneous area
- domain assumption Area increments are Gaussian white noise
- domain assumption Selection of interior-dissolving bubbles does not bias the statistics
- domain assumption Coarse-grained density field with mesh size ξ = 4.375σ faithfully represents bubbles
Cite this review
Pith. "Pith review of Stochastic bubble dynamics in phase-separated scalar active matter." pith.science (2026). https://pith.science/paper/546QWBEP
@misc{pith2026250111442,
author = {Pith},
title = {Pith review of: Stochastic bubble dynamics in phase-separated scalar active matter},
year = {2026},
howpublished = {\url{https://pith.science/paper/546QWBEP}},
note = {Machine review of arXiv:2501.11442}
}
read the original abstract
In ABP systems, phase separation is accompanied by the emergence of vapor bubbles within liquid domains. Using large-scale particle-based simulations, we study the stochastic dynamics of these bubbles and find that most nucleate, grow, and dissolve within liquid domains. We show that their area dynamics can be described by a Langevin equation with a constant negative drift and noise proportional to the perimeter, fully characterizing bubble area and lifetime statistics. Additionally, we develop a lattice gas model that captures the morphological properties, including the decrease in bubble asphericity with increasing area. These findings provide new insights into phase separation in active matter and highlight limitations in current continuum theories.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
Symmetry, thermodynamics, and topology in active matter,
Mark J. Bowick, Nikta Fakhri, M. Cristina Marchetti, and Sriram Ramaswamy, “Symmetry, thermodynamics, and topology in active matter,” Phys. Rev. X 12, 010501 (2022)
work page 2022
-
[2]
Control of protein-based pattern forma- tion via guiding cues,
Tom Burkart, Manon C. Wigbers, Laeschkir W¨ urthner, and Erwin Frey, “Control of protein-based pattern forma- tion via guiding cues,” Nature Reviews Physics 4, 511– 527 (2022)
work page 2022
-
[3]
Motility-induced phase separation,
Michael E. Cates and Julien Tailleur, “Motility-induced phase separation,” Annual Review of Condensed Matter Physics 6, 219–244 (2015)
work page 2015
-
[4]
Self-organization of Protein Patterns,
Erwin Frey and Fridtjof Brauns, “Self-organization of Protein Patterns,” in Active Matter and Nonequilibrium Statistical Physics: Lecture Notes of the Les Houches Summer School: Volume 112, September 2018 (Oxford University Press, 2022)
work page 2018
-
[5]
Athermal phase separation of self-propelled particles with no alignment,
Yaouen Fily and M. Cristina Marchetti, “Athermal phase separation of self-propelled particles with no alignment,” Phys. Rev. Lett. 108, 235702 (2012)
work page 2012
-
[6]
Structure and dynamics of a phase-separating active colloidal fluid,
Gabriel S. Redner, Michael F. Hagan, and Aparna Baskaran, “Structure and dynamics of a phase-separating active colloidal fluid,” Phys. Rev. Lett. 110, 055701 (2013)
work page 2013
-
[7]
Dy- namical clustering and phase separation in suspensions of self-propelled colloidal particles,
Ivo Buttinoni, Julian Bialk´ e, Felix K¨ ummel, Hartmut L¨ owen, Clemens Bechinger, and Thomas Speck, “Dy- namical clustering and phase separation in suspensions of self-propelled colloidal particles,” Phys. Rev. Lett. 110, 238301 (2013)
work page 2013
-
[8]
Continuum theory of phase separation kinetics for ac- tive brownian particles,
Joakim Stenhammar, Adriano Tiribocchi, Rosalind J. Allen, Davide Marenduzzo, and Michael E. Cates, “Continuum theory of phase separation kinetics for ac- tive brownian particles,” Phys. Rev. Lett. 111, 145702 (2013)
work page 2013
Show all 27 references
-
[9]
Free Energy of a Nonuniform System. I. Interfacial Free Energy,
John W. Cahn and John E. Hilliard, “Free Energy of a Nonuniform System. I. Interfacial Free Energy,” The Journal of Chemical Physics 28, 258–267 (1958)
1958
-
[10]
Theory of dy- namic critical phenomena,
P. C. Hohenberg and B. I. Halperin, “Theory of dy- namic critical phenomena,” Rev. Mod. Phys.49, 435–479 (1977)
1977
-
[11]
Phase behaviour of ac- tive brownian particles: the role of dimensionality,
Joakim Stenhammar, Davide Marenduzzo, Rosalind J. Allen, and Michael E. Cates, “Phase behaviour of ac- tive brownian particles: the role of dimensionality,” Soft Matter 10, 1489–1499 (2014)
2014
-
[12]
Curvature-dependent tension and tangential flows at the interface of motility-induced phases,
Adam Patch, Daniel M. Sussman, David Yllanes, and M. Cristina Marchetti, “Curvature-dependent tension and tangential flows at the interface of motility-induced phases,” Soft Matter 14, 7435–7445 (2018)
2018
-
[13]
Cluster phases and bubbly phase separation in active fluids: Reversal of the ostwald process,
Elsen Tjhung, Cesare Nardini, and Michael E. Cates, “Cluster phases and bubbly phase separation in active fluids: Reversal of the ostwald process,” Phys. Rev. X 8, 031080 (2018)
2018
-
[14]
Statistical properties of microphase and bubbly phase- separated active fluids,
Giordano Fausti, Michael E. Cates, and Cesare Nardini, “Statistical properties of microphase and bubbly phase- separated active fluids,” Phys. Rev. E 110, L042103 (2024)
2024
-
[15]
Negative interfacial tension in phase- separated active brownian particles,
Julian Bialk´ e, Jonathan T. Siebert, Hartmut L¨ owen, and Thomas Speck, “Negative interfacial tension in phase- separated active brownian particles,” Phys. Rev. Lett. 115, 098301 (2015)
2015
-
[16]
Role of Repulsive Forces in Determining the Equilib- rium Structure of Simple Liquids,
John D. Weeks, David Chandler, and Hans C. Andersen, “Role of Repulsive Forces in Determining the Equilib- rium Structure of Simple Liquids,” The Journal of Chem- ical Physics 54, 5237–5247 (1971)
1971
-
[17]
Morris, and Sylvie Pic, A Physical Introduction to Suspension Dynamics , Cam- bridge Texts in Applied Mathematics (Cambridge Uni- versity Press, 2011)
´Elisabeth Guazzelli, Jeffrey F. Morris, and Sylvie Pic, A Physical Introduction to Suspension Dynamics , Cam- bridge Texts in Applied Mathematics (Cambridge Uni- versity Press, 2011). 6
2011
-
[18]
Modeling active col- loids: From active brownian particles to hydrodynamic and chemical fields,
Andreas Z¨ ottl and Holger Stark, “Modeling active col- loids: From active brownian particles to hydrodynamic and chemical fields,” Annual Review of Condensed Mat- ter Physics 14, 109–127 (2023)
2023
-
[19]
Supplemental material for “title of the paper
“Supplemental material for “title of the paper”,” Available online (2024), see Supplemental Material at link-provided-after-publication for further details and an additional analysis of the model, which includes Refs. [???]
2024
-
[20]
¨Uber Brownsche Molekularbewegung unter Einwirkung ¨ ausserer Kr¨ afte und deren Zusam- menhang mit der veralgemeinerten Diffusionsgleichung,
M Smoluchowski, “ ¨Uber Brownsche Molekularbewegung unter Einwirkung ¨ ausserer Kr¨ afte und deren Zusam- menhang mit der veralgemeinerten Diffusionsgleichung,” Ann. d. Phys 48, 1103 (1916)
1916
-
[21]
Zur Theorie der Fall- und Steigversuche an Teilchen mit Brownscher Bewegung,
E. Schr¨ odinger, “Zur Theorie der Fall- und Steigversuche an Teilchen mit Brownscher Bewegung,” Physikalische Zeitschrift 16, 189–295 (1915)
1915
-
[22]
Notiz ¨ uber die Berechnung der Brownschen Molekularbewegung bei der Ehrenhaft- Millikanschen Versuchsanordnung,
M. v. Smoluchowski, “Notiz ¨ uber die Berechnung der Brownschen Molekularbewegung bei der Ehrenhaft- Millikanschen Versuchsanordnung,” Physikalische Zeitschrift 16, 318–321 (1915)
1915
-
[23]
The inverse gaussian distribution,
Albert W. Marshall and Ingram Olkin, “The inverse gaussian distribution,” in Life Distributions: Structure of Nonparametric, Semiparametric, and Parametric Fami- lies (Springer New York, New York, NY, 2007) pp. 451– 471
2007
-
[24]
The aspherity of random walks,
J Rudnick and G Gaspari, “The aspherity of random walks,” Journal of Physics A: Mathematical and General 19, L191 (1986)
1986
-
[25]
The Shapes of Random Walks,
Joseph Rudnick and George Gaspari, “The Shapes of Random Walks,” Science 237, 384–389 (1987)
1987
-
[26]
Shapes of semiflexible poly- mer rings,
Karen Alim and Erwin Frey, “Shapes of semiflexible poly- mer rings,” Phys. Rev. Lett. 99, 198102 (2007)
2007
-
[27]
Active phase separa- tion: new phenomenology from non-equilibrium physics,
M. E. Cates and C. Nardini, “Active phase separa- tion: new phenomenology from non-equilibrium physics,” (2024), arXiv:2412.02854 [cond-mat.soft]
2024 arXiv
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.