Large temperature-up-jump simulations of a binary Lennard-Jones system
Pith reviewed 2026-05-16 11:23 UTC · model grok-4.3
The pith
The Tool-Narayanaswamy material-time prediction collapses autocorrelation functions better after moderate temperature up-jumps than after large ones in a binary Lennard-Jones liquid.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The triangular relation of the potential energy is well obeyed, allowing definition of a potential-energy-based material time ξ. The TN material-time prediction that all time-autocorrelation functions collapse to depend only on the material-time difference ξ₂ − ξ₁ is found to work better for the 0.43 to 0.48 temperature jump than for the 0.37 to 0.48 jump.
What carries the argument
The potential-energy-based material time ξ defined from the triangular relation of the potential energy, used to rescale time in aging curves.
If this is right
- The triangular relation holds after both moderate and large jumps, enabling a global material time.
- All five monitored quantities follow the material-time scaling for the smaller temperature jump.
- The formalism performs best when the system is never very far from equilibrium.
- The results are consistent with the general understanding that TN aging works well only for moderate deviations from equilibrium.
Where Pith is reading between the lines
- Allowing each observable its own material time might restore collapse for larger jumps.
- A locally defined material time could better capture dynamic heterogeneity after large jumps.
- The same protocol could be tested on other glass formers to check whether the moderate-jump preference is universal.
Load-bearing premise
The triangular relation for potential energy remains accurate enough to define a single global material time even after large temperature jumps that take the system far from equilibrium.
What would settle it
Failure of the self-intermediate scattering function or mean-square displacement to collapse onto a single master curve when plotted against material-time difference for the 0.43 to 0.48 jump would disprove the central claim.
Figures
read the original abstract
This paper presents simulations of the physical aging of a binary Kob-Andersen-type Lennard-Jones liquid following large temperature up-jumps from equilibrated states of high relaxation time. The purpose is to investigate how well the Tool-Narayanaswamy (TN) material-time concept works for this rather extreme case of aging. First the triangular relation of the potential energy is investigated. This is found to be well obeyed, making it possible to define a potential-energy-based material time $\xi$. We proceed to study aging toward equilibrium at the final temperature 0.48 for jumps from the two temperatures 0.43 and 0.37 (primarily), monitoring the following five quantities: the potential energy, the self-intermediate scattering function, the mean-square displacement, the dynamic susceptibility $\chi_4$, and the non-Gaussian parameter $\alpha_2$. The TN material-time prediction is that all time-autocorrelation functions should collapse to only depend on the material-time difference $\xi_2-\xi_1$. This is found to work better for the $0.43\to 0.48$ temperature jump than for the $0.37\to 0.48$ jump. Our findings thus confirm the general understanding that the TN aging formalism works best for systems that are never very far from equilibrium. This raises two questions for future work: Is the collapse significantly improved if each aging quantity is allowed its own material time? Can better collapse be obtained if the material-time is generalized to be locally defined (in order to reflect dynamic heterogeneity)?
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports MD simulations of physical aging in a binary Kob-Andersen Lennard-Jones liquid after large temperature up-jumps from equilibrated initial states at T=0.43 and T=0.37 to a final T=0.48. It first verifies that the potential energy obeys the triangular relation to within simulation accuracy, permitting definition of a single global material time ξ based on potential energy. It then monitors five observables (potential energy, self-intermediate scattering function, mean-square displacement, χ4, and α2) during aging and tests the TN prediction that all time-autocorrelation functions collapse onto curves depending only on the material-time difference ξ2−ξ1. The collapse is reported to be visibly better for the smaller (0.43→0.48) jump than for the larger (0.37→0.48) jump, consistent with the expectation that the TN formalism works best near equilibrium.
Significance. If the reported trends hold, the work supplies concrete numerical evidence on the range of validity of the Tool-Narayanaswamy material-time concept under large temperature jumps. The explicit comparison of collapse quality between two jump sizes quantifies how far from equilibrium the formalism remains useful and identifies two concrete directions for improvement (quantity-specific material times and locally defined material times). The independent verification of the triangular relation for potential energy provides a reproducible foundation for the material-time construction.
minor comments (2)
- [Results on collapse of dynamical quantities] The quality of the data collapse onto ξ is assessed visually; no quantitative metric (e.g., mean-squared deviation from a master curve or χ² per degree of freedom) or error bars on the collapse residuals are reported, which would allow a statistically grounded comparison between the two jumps.
- [Section on triangular relation] The manuscript does not state whether alternative choices for the material-time variable (e.g., based on pressure or on a different observable) were tested before settling on potential energy; an explicit statement that the triangular relation was checked for other quantities would strengthen the justification for the chosen definition of ξ.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the positive assessment. The referee summary correctly describes the simulations of physical aging in the binary Lennard-Jones liquid after large temperature up-jumps and the testing of the Tool-Narayanaswamy material-time concept. We appreciate the recognition of the significance of our findings in quantifying the range of validity of the TN formalism. As no specific major comments requiring changes were provided, we do not propose any revisions.
Circularity Check
No significant circularity detected
full rationale
The paper first verifies the triangular relation for potential energy directly in the simulations and uses the verified relation to define the material time ξ. The central test—whether autocorrelation functions collapse as a function of ξ differences—is then performed as an independent empirical comparison against this externally constructed time coordinate. No step reduces by construction to a fitted parameter, self-definition, or load-bearing self-citation; the collapse quality is reported as a simulation outcome that can be falsified by the data. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (2)
- Final temperature 0.48
- Initial temperatures 0.43 and 0.37
axioms (2)
- domain assumption The triangular relation for potential energy holds for the chosen temperature jumps.
- standard math Standard molecular-dynamics assumptions (Newtonian dynamics, periodic boundaries, Nose-Hoover thermostat).
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The TN material-time prediction is that all time-autocorrelation functions should collapse to only depend on the material-time difference ξ2−ξ1. This is found to work better for the 0.43→0.48 temperature jump than for the 0.37→0.48 jump.
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]
The (normalized) potential-energy time-autocorrelation functionC uu(t1, t2) defined in Eq. (4)
-
[2]
The incoherent intermediate scattering functionF s(t1, t2) calculated from the particle displacement from time t1 to timet 2, evaluated at the wavevector corresponding to the first maximum of the radial distribution function determined at the final temperature 0.48 (the peak position was virtually constant throughout all simulations), Fs(t1, t2) = * 1 NA ...
-
[3]
The mean-square displacement between timest 1 andt 2, ∆r2 (t1, t2) = * 1 NA NAX i=1 (ri(t2)−r i(t1))2 + .(6)
-
[4]
The dynamic susceptibilityχ 4(t1, t2) quantifying the degree of dynamic heterogeneity in terms of the variance of the incoherent intermediate scattering function [33], χ4(t1, t2, q) =N A "* 1 NA NAX i=1 cos(q·(r i(t2)−r i(t1) !2+ −F 2 s (t1, t2) # ; (7)
-
[5]
The non-Gaussian parameterα 2(t1, t2) (which is zero for Gaussian diffusion), α2(t1, t2) = 3 (∆r)4(t1, t2) 5 ⟨(∆r)2(t1, t2)⟩ 2 −1.(8) 8 0.00 0.25 0.50 Cuu(t1, t2) (a) T0.43→0.48 (b) T0.37→0.48 0.0 0.5 Fs(t1, t2) (c) (d) 10□2 10□1 100 101 ⟨(∆r)2⟩(t1, t2) (e) (f) 100 102 χ4(t1, t2) (g) (h) 100 102 10410□2 10□1 100 101 α2(t1, t2) (i) 100 102 104 (j) 100 102 ...
-
[6]
Mazurin, Relaxation phenomena in glass, J
O. Mazurin, Relaxation phenomena in glass, J. Non-Cryst. Solids25, 129 (1977)
work page 1977
-
[7]
G. W. Scherer,Relaxation in Glass and Composites(Wiley, New York, 1986)
work page 1986
-
[8]
T. Hecksher, N. B. Olsen, K. Niss, and J. C. Dyre, Physical aging of molecular glasses studied by a device allowing for rapid thermal equilibration, J. Chem. Phys.133, 174514 (2010)
work page 2010
-
[9]
T. Hecksher, N. B. Olsen, and J. C. Dyre, Communication: Direct tests of single-parameter aging, J. Chem. Phys.142, 241103 (2015)
work page 2015
- [10]
-
[11]
B. Riechers, L. A. Roed, S. Mehri, T. S. Ingebrigtsen, T. Hecksher, J. C. Dyre, and K. Niss, Predicting nonlinear physical aging of glasses from equilibrium relaxation via the material time, Sci. Adv.8, eabl9809 (2022)
work page 2022
-
[12]
K. Moch, C. Gainaru, and R. B¨ ohmer, Nonlinear susceptibilities and higher-order responses related to physical aging: Wiener–Volterra approach and extended Tool–Narayanaswamy–Moynihan models, J. Chem. Phys.161, 014502 (2024)
work page 2024
-
[13]
O. S. Narayanaswamy, A model of structural relaxation in glass, J. Amer. Ceram. Soc.54, 491 (1971). 10
work page 1971
-
[14]
Micoulaut, Relaxation and physical aging in network glasses: a review, Rep
M. Micoulaut, Relaxation and physical aging in network glasses: a review, Rep. Prog. Phys.79, 066504 (2016)
work page 2016
-
[15]
I. M. Hodge, Physical aging in polymer glasses, Science267, 1945 (1995)
work page 1945
-
[16]
J. M. Hutchinson, Physical aging of polymers, Prog. Polym. Sci.20, 703 (1995)
work page 1995
-
[17]
D. Cangialosi, V. M. Boucher, A. Alegria, and J. Colmenero, Physical aging in polymers and polymer nanocomposites: recent results and open questions, Soft Matter9, 8619 (2013)
work page 2013
-
[18]
G. B. McKenna and S. L. Simon, 50th anniversary perspective: Challenges in the dynamics and kinetics of glass-forming polymers, Macromolecules50, 6333 (2017)
work page 2017
-
[19]
B. Ruta, E. Pineda, and Z. Evenson, Relaxation processes and physical aging in metallic glasses, J. Phys.: Condens. Matter 29, 503002 (2017)
work page 2017
-
[20]
L. Song, W. Xu, J. Huo, F. Li, L.-M. Wang, M. D. Ediger, and J.-Q. Wang, Activation entropy as a key factor controlling the memory effect in glasses, Phys. Rev. Lett.125, 135501 (2020)
work page 2020
-
[21]
J. C. Mauro,Materials Kinetics: Transport and Rate Phenomena(Elsevier, Amsterdan, Netherlands, 2021)
work page 2021
-
[22]
A. Vila-Costa, M. Gonzalez-Silveira, C. Rodriguez-Tinoco, M. Rodriguez-Lopez, and J. Rodriguez-Viejo, Emergence of equilibrated liquid regions within the glass, Nat. Phys.19, 114 (2023)
work page 2023
- [23]
-
[24]
R. F. Lancelotti, L. Pereira, R. B. Rosante, M. Andersen, K.-U. Hess, D. B. Dingwell, E. D. Zanotto, and S. Sen, Relaxation kinetics in temperature jump experiments above calorimetric glass transition via fast scanning calorimetry, J. Am. Ceram. Soc.109, e70618 (2026)
work page 2026
-
[25]
A. J. Kovacs, Transition vitreuse dans les polymeres amorphes. Etude phenomenologique, Fortschr. Hochpolym.-Forsch. 3, 394 (1963)
work page 1963
-
[26]
J. Malek, How to distinguish nonexponentiality and nonlinearity in isothermal structural relaxation of glass-forming ma- terials, J. Phys. Chem. B128, 8074 (2024)
work page 2024
-
[27]
K. Niss, J. C. Dyre, and T. Hecksher, Long-time structural relaxation of glass-forming liquids: Simple or stretched expo- nential?, J. Chem. Phys.152, 041103 (2020)
work page 2020
-
[28]
A. Q. Tool, Relation between inelastic deformability and thermal expansion of glass in its annealing range, J. Amer. Ceram. Soc.29, 240 (1946)
work page 1946
-
[29]
I. M. Douglass and J. C. Dyre, Distance-as-time in physical aging, Phys. Rev. E106, 054615 (2022)
work page 2022
-
[30]
L. F. Cugliandolo and J. Kurchan, On the out-of-equilibrium relaxation of the Sherrington-Kirkpatrick model, J. Phys. A: Math. Gen.27, 5749 (1994)
work page 1994
- [31]
-
[32]
T. B. Schrøder and J. C. Dyre, Solid-like mean-square displacement in glass-forming liquids, J. Chem. Phys.152, 141101 (2020)
work page 2020
-
[33]
U. R. Pedersen, T. B. Schrøder, and J. C. Dyre, Phase diagram of Kob-Andersen-type binary Lennard-Jones mixtures, Phys. Rev. Lett.120, 165501 (2018)
work page 2018
-
[34]
N. P. Bailey, T. S. Ingebrigtsen, J. S. Hansen, A. A. Veldhorst, L. Bøhling, C. A. Lemarchand, A. E. Olsen, A. K. Bacher, L. Costigliola, U. R. Pedersen, H. Larsen, J. C. Dyre, and T. B. Schrøder, RUMD: A general purpose molecular dynamics package optimized to utilize GPU hardware down to a few thousand particles, Scipost Phys.3, 038 (2017)
work page 2017
-
[35]
L. F. Cugliandolo, J. Kurchan, and L. Peliti, Energy flow, partial equilibration, and effective temperatures in systems with slow dynamics, Phys. Rev. E55, 3898 (1997)
work page 1997
-
[36]
T. B¨ ohmer, J. P. Gabriel, L. Costigliola, J.-N. Kociok, T. Hecksher, J. C. Dyre, and T. Blochowicz, Time reversibility during the ageing of materials, Nat. Phys.20, 637 (2024)
work page 2024
-
[37]
C. Chamon and L. F. Cugliandolo, Fluctuations in glassy systems, J. Stat. Mech.7, P07022 (2007)
work page 2007
-
[38]
C. Toninelli, M. Wyart, L. Berthier, G. Biroli, and J.-P. Bouchaud, Dynamical susceptibility of glass formers: Contrasting the predictions of theoretical scenarios, Phys. Rev. E71, 041505 (2005)
work page 2005
-
[39]
G. Biroli, P. Charbonneau, G. Folena, Y. Hu, and F. Zamponi, Local dynamical heterogeneity in simple glass formers, Physical Review Letters128, 10.1103/PhysRevLett.128.175501 (2022)
-
[40]
H. E. Castillo and A. Parsaeian, Local fluctuations in the ageing of a simple structural glass, Nat. Phys.3, 26 (2007)
work page 2007
-
[41]
M. Ruiz-Ruiz, A. Vila-Costa , T. Bar , C. Rodriguez-Tinoco , M. Gonzalez-Silveira , J. A. Plaza , J. Alcala , J. Fraxedas , and J. Rodriguez-Viejo, Real-time microscopy of the relaxation of a glass, Nat. Phys.19, 1509 (2023)
work page 2023
-
[42]
C. Herrero, C. Scalliet, M. D. Ediger, and L. Berthier, Two-step devitrification of ultrastable glasses, PNAS120, e2220824120 (2023)
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.