REVIEW 2 major objections 5 minor 41 references
Influence of active breathing on rheology and jamming of amorphous solids: insights from microscopic and mesoscale analysis
T0 review · 2 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Periodic breathing of particles can fluidize a jammed amorphous solid once the oscillation amplitude exceeds a critical value.
desk verdict Solid qualitative case for activity-induced fluidization via breathing, but the headline threshold and exponent need stronger low-shear-rate convergence checks before I'd trust them. 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 central object is the dimensionless amplitude a of sinusoidal diameter oscillations, with period T = 820τ0, applied to each particle with a phase choice that keeps the packing fraction constant. The carrying mechanism is the combination of local strain injection—each breathing particle acts as a dilating or contracting elastic inclusion—with long-range stress redistribution. In the mesoscale active elasto-plastic model (AEPM), the total stress is the sum of external, active, and plastic contributions; the active contribution is computed by convolving the active deformation field with an elastic propagator F derived from the linear response to an isotropic inclusion, while plastic events relax stress through the Eshelby kernel G. The agreement between microscopic and mesoscale results identifies this elastic propagator response as the ingredient that carries the fluidization.
What would settle it
Extend the shear-rate-controlled simulations to rates below $10^{-7}$, or run long enough to verify convergence at the current lowest rates; if the apparent yield stress continues to decrease and a Newtonian plateau appears for amplitudes below 0.0506, the reported threshold and yield-stress picture would need revision.
Extended reading notes
Core claim
The central discovery is an activity-induced unjamming transition driven by particle size oscillations. In the microscopic model, the number of plastic events per breathing cycle stays low and nearly constant until a threshold amplitude, then rises sharply, while the von Mises stress passes through a cusp at the same amplitude. Under imposed shear, flow curves retain Herschel-Bulkley (yield-stress) form for low amplitudes, and the fitted yield stress decreases continuously with amplitude, extrapolating to zero at a* = 0.0506 with an exponent of about 0.39. Above the threshold, flow curves develop a constant viscosity at low shear rates, the signature of Newtonian flow. The mesoscale elasto-plastic model reproduces this softening and the yield-stress-to-Newtonian crossover by coupling local active dilations and contractions to stress redistribution through an elastic propagator.
Load-bearing premise
The claimed yield-stress regime below threshold rests on Herschel-Bulkley fits to flow curves measured only down to shear rates near $10^{-7}$, with no demonstrated convergence of the steady state at the lowest rates.
Editorial extensions
If this is right
- At low breathing amplitudes the jammed solid keeps a finite yield stress, and the yield stress falls continuously as amplitude increases until rigidity disappears at a* = 0.0506.
- Above the critical amplitude, the shear-rate-dependent viscosity develops a constant plateau at low rates, meaning the unsheared system is a Newtonian fluid.
- The fluidization is generated entirely by internal activity, with no external forcing or thermal fluctuations.
- The mesoscale model matches the microscopic trends, indicating that only local strain injection plus long-range elastic coupling are needed to capture the transition.
- Creep tests under constant stress confirm the picture: at low activity strain saturates, while at higher activity it flows steadily.
Reading between the lines
- If the mechanism is generic, experimental systems with controllable size oscillations—such as responsive colloids or cells—should show the same yield-stress drop and Newtonian plateau in oscillatory rheology.
- The fitted exponent of about 0.39 for the yield stress vanishing near the threshold suggests a possible scaling law that could be tested across different particle interactions, softness, or packing fractions.
- The same mesoscale framework could be applied to other internal activators such as swelling, contraction, or local contractility, which would link it to tissue remodeling and active gel theories.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript studies dense two-dimensional assemblies of soft particles whose diameters oscillate periodically in time ('active breathing'), using large-scale Langevin dynamics simulations (N = 16,000; packing fraction 0.94) combined with a mesoscale tensorial elasto-plastic model in which local dilations are coupled to long-range elastic stress redistribution. The main claims are: (i) for breathing amplitudes below a critical value a* ≈ 0.0506 the system is a yield-stress solid whose yield stress decreases with amplitude and vanishes as (a* − a)^ζ with ζ ≈ 0.39; (ii) above a* the low-shear-rate response becomes Newtonian, with a constant-viscosity plateau; and (iii) the mesoscale model reproduces the same qualitative softening and fluidization, identifying internally generated stresses as the fluidization mechanism. Supporting evidence includes a sharp rise in plastic events per breathing cycle near a ≈ 0.05, a cusp in the von Mises stress, a collapse of σ − σ_y versus shear rate at low amplitudes, and constant-stress creep curves consistent with reduced mechanical resistance.
Significance. The question is timely: whether purely internal, isotropic size oscillations can continuously soften and fluidize a jammed amorphous solid, and how this manifests in rheology. The paper has clear strengths: a well-defined microscopic model; a mesoscale model built from a derived linear-elastic propagator for dilating inclusions (Appendix B), rather than fitted to the MD flow curves; error-barred event statistics in Fig. 3; and a falsifiable structural prediction (collapse of σ − σ_y with γ^{1/2} for a < a*). The qualitative transition is supported by multiple independent observables (events per cycle, von Mises stress, flow-curve shape). However, the headline quantitative outputs (a* and ζ) and the Newtonian-plateau identification rest on Herschel-Bulkley fits to data at shear rates as low as 10^-7 without demonstrated steady state or reported uncertainties, so the precise exponents are not yet established. The mesoscale agreement is qualitative and depends on parameters that are not reported. With the requested convergence checks, error bars, and parameter reporting, the paper would make a solid contribution to the active-matter and jamming literature.
major comments (2)
- [Sec. III B 1, Fig. 4] The central quantitative results—the vanishing yield stress with a* = 0.0506 and ζ ≈ 0.39, and the identification of a Newtonian regime for a ≥ 0.054—are obtained from flow curves extending down to shear rates near 10^-7, but the manuscript gives no evidence that these lowest-rate points are in steady state and provides no error bars on σ_y in Fig. 4(c). Since T = 820 τ_0, a shear rate of 10^-7 accumulates a strain of only about 8 × 10^-5 per breathing cycle; for a < a*, where plastic events are rare (Fig. 3), the run time required to reach a statistically stationary stress is not established, and an under-sampled or transient low-rate point would bias the HB extrapolation of σ_y. That bias would propagate through the power-law fit (a* − a)^ζ, and the situation is delicate because the a = 0.050 point lies within about 1% of the fitted a* = 0.0506. The same caveat applies to the data collapse and the γ^{1/2} claim in Fig. 4(d). The creep measurements in Appendix A, at one applied stress per amplitude and for amplitudes up to only a = 0.048, neither bracket σ_y nor probe the near-threshold region, so they cannot corroborate the extrapolation. Please report the averaging protocol (number of cycles and independent runs), add error bars or confidence intervals to σ_y(a), provide convergence checks (e.g., block averages) at the lowest shear rates, and test the sensitivity of a* and ζ to the fit range and to the inclusion or exclusion of the slowest points.
- [Sec. II B, Fig. 5] The mesoscale results in Fig. 5 and the statement in Sec. III B 2 that the model 'captures this behavior' and 'validates the modeling approach' cannot be properly assessed as presented: the manuscript does not report the yield-threshold distribution (form, mean, width), the relaxation and elastic-recovery timescales τ and τ_el, the active-strain amplitude and its mapping to the MD amplitude a, the system size, or the number of disorder realizations. The paper itself acknowledges that the assumption of zero stress at an active site 'is not very realistic' (Sec. II B). Given the number of free parameters, the qualitative agreement in Fig. 5 is weak evidence unless the parameter values are specified and the trend is shown to be robust to reasonable parameter choices. Please include the full parameter list (reproducing the relevant values rather than only citing Ref. [31]) and a minimal robustness check of the softening-to-fluidization trend.
minor comments (5)
- [Sec. III B 1] The quantity σ_c in δσ_xy = σ_xy − σ_c (Fig. 4(d)) is never defined; it should be stated that σ_c is the Herschel-Bulkley yield stress σ_y.
- [Sec. II B, Eqs. (4)-(5)] The sentence concerning the zero-wavevector components, '−σ_xy(0,0) = σ_xy(0,0) = 0', is garbled and should be rewritten. In addition, the kernel 2 q_x q_y / q^2 has a direction-dependent limit as q → 0, so setting the active stress at q = 0 to zero is a specific regularization of the uniform mode that should be stated and justified.
- [Appendix B, Eq. (B1)] In Eq. (B1) the second unit vector is printed as e_x; it should be e_y.
- [Sec. IV] The concluding paragraph asserts that at low amplitudes plastic events 'remain synchronized with the breathing cycle' and at higher amplitudes 'span multiple cycles,' but no phase-resolved or multi-cycle analysis is presented in Sec. III; either add the supporting analysis or soften these claims.
- [References] Ref. [25] cites the arXiv preprint of the Bonn et al. yield-stress review; the published version (Rev. Mod. Phys. 89, 035005 (2017)) should be cited.
Circularity Check
No circular derivation: the fluidization threshold and rheology are outputs of MD simulation, the mesoscale model is independently constructed from linear elasticity, and the cited prior threshold is an external benchmark.
full rationale
The paper's central claims are not circular. The microscopic model prescribes only the pair potential (Eq. 1), the breathing rule d_i(t)=d_i0(1+a cos(ωt+ψ_i)) (Eq. 2), and athermal Langevin dynamics (Eq. 3). The plastic-event counts, von Mises stress, flow curves, and yield stresses are all measured outputs of the simulation, not inputs. The critical amplitude a*=0.0506 is obtained by fitting the simulated σ_y(a) to (a*−a)^ζ (Sec. III B 1, Fig. 4(c)), and the independent observation of a Newtonian viscosity plateau at large a corroborates the fluidization without being fed into the fit. The comparison with a_c=0.0498 uses Ref. [23] (Tjhung and Berthier), which is external to the present authors and is used as a benchmark, not as a premise. The mesoscale AEPM is built from the linear-elastic response to dilating inclusions derived in Appendix B, plus standard Eshelby and elasto-plastic ingredients (Refs. [31,32]); no MD flow-curve data are injected into the EPM, and the agreement is explicitly qualitative. Self-citations (Refs. [30], [31], [32]) concern numerical methods and model framework and are not load-bearing for the fluidization claim. The Herschel-Bulkley extrapolation at shear rates near 1e-7 raises a legitimate steady-state/error-bar concern, but that is a question of statistical convergence and fit robustness, not circular identification of inputs with outputs.
Assumptions & free parameters
free parameters (5)
- EPM yield threshold distribution (mean and width)
- EPM relaxation timescale tau and elastic recovery time tau_el
- EPM active-strain amplitude and forcing schedule
- Herschel-Bulkley parameters (sigma_y, A, beta) per amplitude a =
sigma_y values shown in Fig. 4(c)
- Rigidity threshold a* and exponent zeta =
a* = 0.0506, zeta = 0.39
assumptions (5)
- domain assumption Uniformly distributed breathing phases psi_i = 2 pi i/N keep the instantaneous packing fraction constant (Eq. 2).
- domain assumption Athermal Langevin dynamics with harmonic repulsive interactions (Eqs. 1 and 3) captures the mechanical response of dense amorphous solids.
- domain assumption A contact change between stroboscopic frames is a faithful definition of a plastic event.
- ad hoc to paper The mesoscale EPM assumes an incompressible medium and uses point-like active inclusions with the stress at the inclusion center set to zero.
- standard math Linear elasticity solutions (Eshelby kernel and the dilation-inclusion response) apply to the mesoscale plastic and active stress propagation.
Cite this review
Pith. "Pith review of Influence of active breathing on rheology and jamming of amorphous solids: insights from microscopic and mesoscale analysis." pith.science (2026). https://pith.science/paper/57LQD7M7
@misc{pith2026250514520,
author = {Pith},
title = {Pith review of: Influence of active breathing on rheology and jamming of amorphous solids: insights from microscopic and mesoscale analysis},
year = {2026},
howpublished = {\url{https://pith.science/paper/57LQD7M7}},
note = {Machine review of arXiv:2505.14520}
}
read the original abstract
We study the flow behavior and unjamming transition in dense assemblies of actively deforming particles that periodically change size, a process that we refer to as breathing. Using extensive molecular dynamics simulations and a complementary mesoscale elasto-plastic model, we explore how this internal activity influences plasticity and rheology. At low amplitudes of breathing, the system remains jammed and displays localized, reversible rearrangements. As the amplitude of the breathing increases beyond a critical threshold, the system undergoes an activity-induced fluidization marked by a surge in plastic events and a drop in yield stress. The flow curve analysis reveals a transition from yield-stress behavior to Newtonian flow at high activity. The mesoscale model captures these trends and provides insight into the role of stress redistribution due to local active deformations. Our findings highlight the potential of internal active driving to tune the mechanical state of amorphous materials without external forcing.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[31]
P. Jaiswal, I. Procaccia, C. Rainone, and M. Singh, Phys- ical Review Letters116, 085501 (2016)
work page 2016
-
[1]
Defining a plastic event Within the microscopic simulations, to develop the phenomenology similar to the elastoplastic model where local active sites are intrinsic descriptors, we need a strat- egy to define a plastic event at the particulate scale. A possible way to define such an event is by identifying whether any of the constituent particles changes n...
-
[2]
Measuring stresses Since the loss of rigidity of the material at some thresh- old breathing amplitude is akin to yielding, we would like to track the evolution of the stress that develops in the system under steady state conditions, for increasinga. Since a breathing event involves the development of dif- ferent stress components as discussed above, we th...
-
[3]
Microscopic model Next we probe the response of the athermal assembly of breathing particles to macroscopic deformation via an imposed shear-rate˙γ. Initially we prepare the steady state configurations corresponding to a breathing ampli- tudeaand subsequently impose the external shear. For every such{a,˙γ}, we obtain the steadily flowing state and record ...
work page 2000
-
[4]
Mesoscale model We show the mesoscale data for rheological response in Fig. 5. While the rheology remains largely unaffected by the activity at high shear rates, as evidenced by the col- lapse of the flow curves in this regime, we observe a pro- nounced decrease in shear stress with increasing activity amplitude at low to intermediate shear rates. This so...
-
[5]
Te Vrugt and R
M. Te Vrugt and R. Wittkowski, The European Physical Journal E48, 12 (2025)
2025
-
[6]
M. C. Marchetti, J.-F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Reviews of Modern Physics85, 1143 (2013)
work page 2013
- [7]
Show all 41 references
-
[8]
M. F. Copeland and D. B. Weibel, Soft Matter5, 1174 (2009)
2009
-
[9]
Theurkauff, C
I. Theurkauff, C. Cottin-Bizonne, J. Palacci, C. Ybert, and L. Bocquet, Physical Review Letters108, 268303 (2012)
2012
-
[10]
Buttinoni, J
I. Buttinoni, J. Bialké, F. Kümmel, H. Löwen, C. Bechinger, and T. Speck, Physical Review Letters 110, 238301 (2013)
2013
-
[11]
Poujade, E
M. Poujade, E. Grasland-Mongrain, A. Hertzog, J. Jouanneau, P. Chavrier, B. Ladoux, A. Buguin, and P. Silberzan, Proceedings of the National Academy of Sciences104, 15988 (2007)
2007
-
[12]
Mehes and T
E. Mehes and T. Vicsek, Integrative Biology6, 831 (2014)
2014
-
[13]
T. E. Angelini, E. Hannezo, X. Trepat, M. Marquez, J. J. Fredberg, and D. A. Weitz, Proceedings of the National Academy of Sciences108, 4714 (2011)
2011
-
[14]
S. M. Rafelski and J. A. Theriot, Annual Review of Bio- chemistry73, 209 (2004)
2004
-
[15]
Sugino, H
Y. Sugino, H. Ebata, Y. Sowa, A. Ikeda, and D. Mizuno, arXiv preprint arXiv:2401.15658 (2024)
2024 arXiv
-
[16]
S. M. Zehnder, M. Suaris, M. M. Bellaire, and T. E. Angelini, Biophysical Journal108, 247 (2015)
2015
-
[17]
Taloni, E
A. Taloni, E. Kardash, O. U. Salman, L. Truskinovsky, S. Zapperi, and C. A. M. La Porta, Physical Review Let- ters114, 208101 (2015)
2015
-
[18]
D. Bi, J. Lopez, J. M. Schwarz, and M. L. Manning, Nature Physics11, 1074 (2015)
2015
-
[19]
S.-Z. Lin, M. Merkel, and J.-F. m. c. Rupprecht, Phys. Rev. Lett.130, 058202 (2023)
2023
-
[20]
Lawson-Keister and M
E. Lawson-Keister and M. L. Manning, Current Opinion 9 in Cell Biology72, 146 (2021)
2021
-
[21]
Ranft, M
J. Ranft, M. Basan, J. Elgeti, J.-F. Joanny, J. Prost, and F. Jülicher, Proceedings of the National Academy of Sciences107, 20863 (2010)
2010
-
[22]
Matoz-Fernandez, E
D. Matoz-Fernandez, E. Agoritsas, J.-L. Barrat, E. Bertin, and K. Martens, Physical review letters118, 158105 (2017)
2017
-
[23]
D. A. Matoz-Fernandez, K. Martens, R. Sknepnek, J.-L. Barrat, and S. Henkes, Soft matter13, 3205 (2017)
2017
-
[24]
Vicsek, A
T. Vicsek, A. Czirók, E. Ben-Jacob, I. Cohen, and O. Shochet, Physical Review Letters75, 1226 (1995)
1995
-
[25]
Henkes, Y
S. Henkes, Y. Fily, and M. C. Marchetti, Physical Review E84, 040301 (2011)
2011
-
[26]
Wiese, K
R. Wiese, K. Kroy, and D. Levis, Physical Review Letters 131, 178302 (2023)
2023
-
[27]
Tjhung and L
E. Tjhung and L. Berthier, Physical Review E96, 050601 (2017)
2017
-
[28]
Li, Q.-L
Z.-Q. Li, Q.-L. Lei, and Y.-Q. Ma, Proceedings of the National Academy of Sciences122, e2421518122 (2025)
2025
-
[29]
Manneville, arXiv preprint arXiv:1502.05281 (2015)
D.Bonn, J.Paredes, M.M.Denn, L.Berthier, T.Divoux, and S. Manneville, arXiv preprint arXiv:1502.05281 (2015)
2015 arXiv
-
[30]
E.D.Knowlton, D.J.Pine, andL.Cipelletti, SoftMatter 10, 6931 (2014)
2014
-
[32]
Heussinger and J.-L
C. Heussinger and J.-L. Barrat, Physical review letters 102, 218303 (2009)
2009
-
[33]
A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, W. M. Brown, P. S. Crozier, P. J. In’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, et al., Com- puter Physics Communications271, 108171 (2022)
2022
-
[34]
Cabriolu, J
R. Cabriolu, J. Horbach, P. Chaudhuri, and K. Martens, Soft matter15, 415 (2019)
2019
-
[35]
Le Goff, Theses, Université Grenoble Alpes (2021), URLhttps://theses.hal.science/tel-03426944
M. Le Goff, Theses, Université Grenoble Alpes (2021), URLhttps://theses.hal.science/tel-03426944
2021
-
[36]
Nicolas, E
A. Nicolas, E. E. Ferrero, K. Martens, and J.-L. Barrat, Reviews of Modern Physics90, 045006 (2018)
2018
-
[37]
Ghosh, P
T. Ghosh, P. Sollich, and S. K. Nandi, Soft Matter21, 3047 (2025)
2025
-
[38]
Nicolas and J.-L
A. Nicolas and J.-L. Barrat, Faraday discussions167, 567 (2013)
2013
-
[39]
Fiocco, G
D. Fiocco, G. Foffi, and S. Sastry, Physical Review E88, 020301 (2013)
2013
-
[40]
Regev, T
I. Regev, T. Lookman, and C. Reichhardt, Physical Re- view E88, 062401 (2013)
2013
-
[41]
Lechner and C
W. Lechner and C. Dellago, Soft Matter5, 646 (2009)
2009
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.