Mode-coupling theory for aging in active glasses: relaxation dynamics and evolution towards steady state
Pith reviewed 2026-05-10 17:56 UTC · model grok-4.3
The pith
Activity modifies the mode-coupling critical point, so the distance to it sets how relaxation time grows with waiting time in active glasses.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The aging solutions of the mode-coupling theory for active glasses show that the two-point correlation function decays more slowly with growing waiting time tw and that the relaxation time tr increases. The activity-modified critical point λC determines the aging: the quench distance from λC governs the dynamics and fixes the exponent δ in the relation tr ∼ tw^δ. This δ decreases with increasing self-propulsion force f0, while the dependence on persistence time τp varies with the nature of the activity.
What carries the argument
The activity-shifted mode-coupling critical point λC, whose distance from the quench controls the aging exponent δ that relates relaxation time to waiting time.
If this is right
- Two-point correlation functions decay more slowly as waiting time grows.
- Relaxation time follows a power law with waiting time, tr ∼ tw^δ.
- The exponent δ decreases when the self-propulsion force f0 is increased.
- How δ changes with persistence time τp depends on the specific form of activity.
Where Pith is reading between the lines
- The same distance-to-critical-point logic could be checked in other non-equilibrium glassy models that add persistent driving.
- Varying f0 and τp independently in particle simulations would provide a direct test of the predicted trends in δ.
- The framework suggests a route to connect aging exponents in active biological matter to an effective distance from a modified glass transition.
Load-bearing premise
The usual mode-coupling approximations remain valid after activity is added, so the critical point can be located without further uncontrolled steps.
What would settle it
Simulations or experiments that measure the aging exponent δ at several values of self-propulsion force f0 while holding the quench distance fixed, to check whether δ decreases as f0 increases.
Figures
read the original abstract
Aging refers to the evolution of system properties with waiting time $t_w$. It is a key feature of glassy dynamics. Recent experiments have demonstrated aging in biological systems that are inherently active with a magnitude of self-propulsion force $f_0$ and a persistence time $\tau_p$. Thus, what governs the aging dynamics in these active systems has fundamental importance. We formulate a generic mode-coupling theory (MCT) of active glasses to address this question. The aging solutions of the theory show that the two-point correlation function decays more slowly with growing $t_w$, and the relaxation time $t_r$ increases. The activity-modification of the MCT critical point, $\lambda_\text{C}$, has profound significance for active aging: the quench distance from $\lambda_\text{C}$ governs aging and determines $\delta$, where $t_r\sim t_w^\delta$. $\delta$ decreases with increasing $f_0$, in agreement with existing simulations. However, the variation with $\tau_p$ depends on the nature of activity. Our work has fundamental theoretical implications for active glasses and paves the way for a deeper understanding of the aging dynamics in biological systems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper formulates a mode-coupling theory (MCT) for aging in active glasses driven by self-propulsion force f0 and persistence time τp. It derives aging solutions for the two-time density correlators showing slower decay with waiting time tw and increasing relaxation time tr. The central result is that an activity-modified critical point λC controls the dynamics: the quench distance from λC determines the aging exponent δ in the power-law relation tr ∼ tw^δ, with δ decreasing as f0 increases (in agreement with simulations) while the τp dependence varies with the nature of activity.
Significance. If the derivation holds, the work provides a first-principles MCT framework for aging in active glasses, extending passive MCT to systems with persistent forces and offering a predictive link between activity parameters and aging exponents. This has direct implications for biological glasses and could guide experiments on living systems. The stated agreement with simulations is a positive feature, though the lack of quantified error estimates limits the strength of that support.
major comments (2)
- [§3] §3 (formulation of the two-time MCT equations): the memory kernel is constructed by inserting the persistent active force term (with finite τp) into the standard mode-coupling vertex derived from the static structure factor. No explicit demonstration is given that non-factorizable correlations between active velocity and density fluctuations are absent or absorbed without affecting the closure; this assumption is load-bearing for the claim that the quench distance from λC alone governs δ(f0, τp).
- [§4.2] §4.2 (aging solutions and δ extraction): the exponent δ is reported to decrease with f0 and to depend on the nature of activity for τp, but the manuscript does not show whether λC is computed independently from the active structure factor or adjusted to match the distance that produces the observed δ; if the latter, the mapping from quench distance to δ becomes circular and the predictive power for varying τp is undermined.
minor comments (2)
- [Abstract] The abstract claims quantitative agreement with simulations but provides no details on the simulated systems, the range of f0 and τp, or error metrics on δ; this should be added to the main text or a supplementary table.
- [§2] Notation for the activity-modified critical point λC is introduced without an explicit equation contrasting it to the passive λC; a dedicated equation in §2 would improve clarity.
Simulated Author's Rebuttal
We thank the referee for their thorough reading and valuable comments, which have helped us improve the clarity of the manuscript. We address each major comment point by point below and have revised the text to strengthen the presentation of the derivation and the independence of the critical point.
read point-by-point responses
-
Referee: §3 (formulation of the two-time MCT equations): the memory kernel is constructed by inserting the persistent active force term (with finite τp) into the standard mode-coupling vertex derived from the static structure factor. No explicit demonstration is given that non-factorizable correlations between active velocity and density fluctuations are absent or absorbed without affecting the closure; this assumption is load-bearing for the claim that the quench distance from λC alone governs δ(f0, τp).
Authors: We agree that an explicit statement of the closure assumptions strengthens the presentation. In the revised manuscript we have added a dedicated paragraph in §3 that justifies the factorization: the active velocity is treated as an external persistent drive whose cross-correlations with density fluctuations are negligible in the overdamped, long-time regime relevant to glassy aging (consistent with the standard MCT treatment of driven systems). This closure is the same as that used in prior active MCT formulations and does not alter the structure of the memory kernel or the resulting dependence on the distance to λC. The added text also references the supporting literature on which the approximation rests. revision: yes
-
Referee: §4.2 (aging solutions and δ extraction): the exponent δ is reported to decrease with f0 and to depend on the nature of activity for τp, but the manuscript does not show whether λC is computed independently from the active structure factor or adjusted to match the distance that produces the observed δ; if the latter, the mapping from quench distance to δ becomes circular and the predictive power for varying τp is undermined.
Authors: λC is obtained entirely from the static structure factor of the active system (computed in §2 via the activity-modified pair correlations) and is inserted into the standard MCT equation for the non-ergodicity parameter; no dynamical fitting or adjustment is performed. The aging exponent δ is then predicted directly from the distance to this independently determined λC. The comparison with simulations is therefore a validation, not an input. To make this separation explicit we have added a new panel to Fig. 4 showing λC(f0, τp) computed from statics alone, together with a short clarifying sentence in §4.2. revision: yes
Circularity Check
No significant circularity detected
full rationale
The derivation inserts the active force term (with parameters f0 and τp) into the standard MCT memory kernel, solves the resulting equations for the non-ergodicity parameter to locate the activity-modified critical point λ_C, and then applies the established MCT scaling relations for the aging regime to obtain the exponent δ from the quench distance to λ_C. No step reduces by construction to a fitted parameter or self-citation; the equations are closed under the usual mode-coupling vertex, and the predicted δ(f0) trend is compared to independent simulations rather than fitted to them. The central claim therefore rests on the internal consistency of the extended MCT equations rather than on any self-referential loop.
Axiom & Free-Parameter Ledger
free parameters (1)
- activity strength f0 and persistence τp
axioms (1)
- domain assumption Mode-coupling closure approximations remain valid for active particles
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the activity-modification of the MCT critical point, λC, has profound significance for active aging: the quench distance from λC governs aging and determines δ, where tr∼tw^δ
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
D(t,s)=2λC²(t,s)+Δ(t−s) and Σ(t,s)=4λC(t,s)R(t,s) with Δ(t)=f0² exp(−t/τp) for ABP
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
L. F. Cugliandolo and J. Kurchan, Phys. Rev. Lett.71, 173 (1993)
work page 1993
-
[2]
S. K. Nandi and S. Ramaswamy, Phys. Rev. Lett.109, 115702 (2012)
work page 2012
-
[3]
S. K. Nandi and S. Ramaswamy, Phys. Rev. E94, 012607 (2016)
work page 2016
-
[4]
P. Lunkenheimer, R. Wehn, U. Schneider, and A. Loidl, Phys. Rev. Lett.95, 055702 (2005)
work page 2005
-
[5]
B. Abou, D. Bonn, and J. Meunier, Phys. Rev. E64, 021510 (2001)
work page 2001
-
[6]
I. M. Hodge, Science267, 1945 (1995)
work page 1945
- [7]
-
[8]
E. Zhou, X. Trepat, C. Park, G. Lenormand, M. Oliver, S. Mijailovich, C. Hardin, D. Weitz, J. Butler, and J. Fredberg, Proc. Nat. Acad. Sci. (USA)106, 10632 (2009)
work page 2009
- [9]
-
[10]
T. E. Angelini, E. Hannezo, X. Trepat, M. Marquez, J. J. Fredberg, and D. A. Weitz, Proc. Nat. Acad. Sci. (USA) 108, 4714 (2011)
work page 2011
-
[11]
J.-A. Park, J. H. Kim, D. Bi, J. A. Mitchel, N. T. Qazvini, K. Tantisira, C. Y. Park, M. McGill, S.-H. Kim, B. Gweon, et al., Nat. Mat.14, 1040 (2015)
work page 2015
-
[12]
C. Malinverno, S. Corallino, F. Giavazzi, M. Bergert, Q. Li, M. Leoni, A. Disanza, E. Frittoli, A. Oldani, E. Martini,et al., Nat. Mat.16, 587 (2017)
work page 2017
-
[13]
L. Atia, D. Bi, Y. Sharma, J. A. Mitchel, B. Gweon, S. A. Koehler, S. J. DeCamp, B. Lan, J. H. Kim, R. Hirsch, et al., Nat. phys.14, 613 (2018)
work page 2018
- [14]
-
[15]
F. Giavazzi, C. Malinverno, G. Scita, and R. Cerbino, Front. Phys.6, 120 (2018)
work page 2018
-
[16]
R. Cerbino, S. Villa, A. Palamidessi, E. Frittoli, G. Scita, and F. Giavazzi, Soft Matter17, 3550 (2021)
work page 2021
-
[17]
K. Nishizawa, K. Fujiwara, M. Ikenaga, N. Nakajo, M. Yanagisawa, and D. Mizuno, Sci. Rep.7, 15143 (2017)
work page 2017
-
[18]
L. Jawerth, E. Fischer-Friedrich, S. Saha, J. Wang, T. Franzmann, X. Zhang, J. Sachweh, M. Ruer, M. Ijavi, S. Saha,et al., Science370, 1317 (2020)
work page 2020
-
[19]
I. Alshareedah, M. M. Moosa, M. Pham, D. A. Potoyan, and P. R. Banerjee, Nat. Commun.12, 6620 (2021)
work page 2021
- [20]
-
[21]
L. Deng, X. Trepat, J. P. Butler, E. Millet, K. G. Morgan, D. A. Weitz, and J. J. Fredberg, Nat. Mat.5, 636 (2006)
work page 2006
- [22]
-
[23]
S. C. Takatori and K. K. Mandadapu, arXiv preprint arXiv:2003.05618 (2020), https://doi.org/10.48550/arXiv.2003.05618
-
[24]
H. Lama, M. J. Yamamoto, Y. Furuta, T. Shimaya, and K. A. Takeuchi, PNAS Nexus , pgae238 (2024)
work page 2024
-
[25]
N. Klongvessa, F. Ginot, C. Ybert, C. Cottin-Bizonne, and M. Leocmach, Phys. Rev. Lett.123, 248004 (2019)
work page 2019
-
[26]
N. Klongvessa, C. Ybert, C. Cottin-Bizonne, T. Kawasaki, and M. Leocmach, J. Chem. Phys. 156, 154509 (2022)
work page 2022
- [27]
- [28]
-
[29]
M. Poujade, E. Grasland-Mongrain, A. Hertzog, J. Jouanneau, P. Chavrier, B. Ladoux, A. Buguin, and P. Silberzan, Proc. Nat. Acad. Sci. (USA)104, 15988 (2007)
work page 2007
-
[30]
A. Brugu´ es, E. Anon, V. Conte, J. H. Veldhuis, M. Gupta, J. Colombelli, J. J. Mu˜ noz, G. W. Brodland, B. Ladoux, and X. Trepat, Nat. phys.10, 683 (2014)
work page 2014
- [31]
-
[32]
D. T. Tambe, C. Corey Hardin, T. E. Angelini, K. Ra- jendran, C. Y. Park, X. Serra-Picamal, E. H. Zhou, M. H. Zaman, J. P. Butler, D. A. Weitz,et al., Nat. Mat.10, 469 (2011)
work page 2011
-
[33]
A. N. Malmi-Kakkada, X. Li, H. S. Samanta, S. Sinha, and D. Thirumalai, Phys. Rev. X8, 021025 (2018)
work page 2018
-
[34]
S. Sadhukhan, S. Dey, S. Karmakar, and S. K. Nandi, Euro. Phys. J. Spe. Top.233, 3193 (2024)
work page 2024
-
[35]
P. Maiuri, J.-F. Rupprecht, S. Wieser, V. Ruprecht, O. B´ enichou, N. Carpi, M. Coppey, S. De Beco, N. Gov, C.-P. Heisenberg,et al., Cell161, 374 (2015)
work page 2015
-
[36]
I. M. Wortel, I. Niculescu, P. M. Kolijn, N. S. Gov, R. J. de Boer, and J. Textor, Biophys. J.120, 2609 (2021)
work page 2021
-
[37]
L. M. Janssen, J. Phys.: Condens. Matter31, 503002 (2019)
work page 2019
-
[38]
L. Berthier, E. Flenner, and G. Szamel, J. Chem. Phys. 150, 200901 (2019)
work page 2019
- [39]
- [40]
-
[41]
A. Liluashvili, J. ´Onody, and T. Voigtmann, Phys. Rev. E96, 062608 (2017)
work page 2017
- [42]
-
[43]
V. E. Debets and L. M. C. Janssen, J. Chem. Phys.157, 224902 (2022)
work page 2022
-
[44]
S. K. Nandi and N. S. Gov, Soft matter13, 7609 (2017)
work page 2017
-
[45]
S. K. Nandi, R. Mandal, P. J. Bhuyan, C. Dasgupta, M. Rao, and N. S. Gov, Proc. Nat. Acad. Sci. (USA) 115, 7688 (2018). 9
work page 2018
-
[46]
K. Paul, A. Mutneja, S. K. Nandi, and S. Karmakar, Proc. Nat. Acad. Sci. (USA120, e2217073120 (2023)
work page 2023
- [47]
- [48]
- [49]
- [50]
-
[51]
V. E. Debets, X. M. D. Wit, and L. M. Janssen, Phys. Rev. Lett.127, 278002 (2021)
work page 2021
- [52]
-
[53]
Y.-E. Keta, R. Mandal, P. Sollich, R. L. Jack, and L. Berthier, Soft Matter19, 3871 (2023)
work page 2023
-
[54]
Y.-E. Keta, R. L. Jack, and L. Berthier, Phys. Rev. Lett. 129, 048002 (2022)
work page 2022
- [55]
- [56]
- [57]
-
[58]
D. Bonn, S. Tanase, B. Abou, H. Tanaka, and J. Meu- nier, Phys. Rev. Lett.89, 015701 (2002)
work page 2002
-
[59]
X. Di, K. Win, G. B. McKenna, T. Narita, F. Lequeux, S. R. Pullela, and Z. Cheng, Phys. Rev. Lett.106, 095701 (2011)
work page 2011
- [60]
-
[61]
B. Riechers, L. A. Roed, S. Mehri, T. S. Ingebrigtsen, T. Hecksher, J. C. Dyre, and K. Niss, Science Advances 8, eabl9809 (2022)
work page 2022
- [62]
- [63]
- [64]
- [65]
- [66]
- [67]
-
[68]
T. Castellani and A. Cavagna, Journal of Statistical Mechanics: Theory and Experiment2005, P05012 (2005)
work page 2005
-
[69]
D. R. Reichman and P. Charbonneau, J. Stat. Mech.: Theor. Exp.2005, P05013 (2005)
work page 2005
- [70]
-
[71]
J.-P. Bouchaud, L. Cugliandolo, J. Kurchan, and M. M´ ezard, Physica A226, 243 (1996)
work page 1996
-
[72]
G¨ otze,Complex dynamics of glass-forming liquids: A mode-coupling theory, Vol
W. G¨ otze,Complex dynamics of glass-forming liquids: A mode-coupling theory, Vol. 143 (Oxford University Press, USA, 2009)
work page 2009
-
[73]
E. Zaccarelli, G. Foffi, F. Sciortino, P. Tartaglia, and K. A. Dawson, Europhys. Lett.55, 157 (2001)
work page 2001
- [74]
-
[75]
N. Klongvessa, C. Ybert, C. Cottin-Bizonne, T. Kawasaki, and M. Leocmach, J. Chem. Phys. 156(2022)
work page 2022
- [76]
- [77]
- [78]
-
[79]
R. Farhadifar, J.-C. R¨ oper, B. Aigouy, S. Eaton, and F. J¨ ulicher, Curr. Biol.17, 2095 (2007)
work page 2095
-
[80]
A. G. Fletcher, M. Osterfield, R. E. Baker, and S. Y. Shvartsman, Biophys. J.106, 2291 (2014)
work page 2014
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.