Recognition: 2 theorem links
· Lean TheoremIdentifying the relevant parameters in design strategies for stable glasses
Pith reviewed 2026-05-13 03:34 UTC · model grok-4.3
The pith
Diameter dynamics, not optimized physical properties, produce ultrastable glasses.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors introduce computational methods to optimize physical quantities such as enhanced hyperuniformity at large scales and local ordering at small scales without any dynamical changes to particle diameters. The resulting glass configurations display values of these quantities far beyond those of bulk equilibrium states. However, the same configurations exhibit no increase in stability when assessed with standard metrics. This demonstrates that the targeted physical quantities are correlated with stability only in the presence of diameter dynamics and are not causally responsible for ultrastability.
What carries the argument
Computational optimization methods for hyperuniformity and local ordering that operate without changing particle diameters, used to separate the effect of these quantities from the dynamical process of diameter variation.
If this is right
- Design strategies for stable glasses must focus on the dynamical processes of preparation rather than on achieving specific final physical properties.
- Diameter-changing algorithms achieve enhanced stability through the act of varying diameters, not through the hyperuniformity or ordering they produce.
- Optimization of physical quantities alone, when decoupled from diameter dynamics, is insufficient to generate ultrastable glasses.
- Previous interpretations of design rules for stable glasses require revision to emphasize generating dynamics over target properties.
Where Pith is reading between the lines
- The same logic separating correlation from causation may apply to other material-design problems where algorithms optimize structural metrics through altered dynamics.
- Experimental vapor-deposition techniques could be re-examined to determine whether surface-level effective diameter variations contribute to observed stability.
- New algorithms that combine non-diameter optimizations with alternative dynamical variables could be tested to isolate the minimal requirements for stability enhancement.
Load-bearing premise
The stability of glass configurations produced by the new diameter-free optimization methods can be fairly compared to those from diameter-changing algorithms using identical stability metrics.
What would settle it
A simulation in which glasses optimized for hyperuniformity and local order without diameter changes display stability metrics equal to or higher than those achieved by diameter-changing methods when prepared and measured under the same conditions.
Figures
read the original abstract
A glass is conventionally obtained by cooling a bulk supercooled liquid through its glass transition temperature. The discovery of ultrastable glasses prepared using physical vapor deposition, together with the recent multiplication of numerical algorithms created to increase the stability of glasses, demonstrates the existence of a variety of strategies for designing glasses with different physical properties. This raises a broader question: which parameters most strongly govern the enhancement of glass stability? Existing computational strategies often produce highly stable glasses by optimizing certain physical properties through dynamical changes in particle diameters. We challenge the idea that these physical quantities are causally responsible for glass stability and suggest instead that diameter dynamics is the principal source of enhanced stability. To support our view, we introduce computational methods to optimize physical quantities without changing the particle diameters. Using the examples of enhanced hyperuniformity at large scale and local ordering at small scale, we design glass configurations with highly optimized values compared to bulk equilibrium states. However, these glasses do not show enhanced stability. The proposed physical quantities are correlated with glass stability, but are not causally responsible for ultrastability. These findings indicate that design rules for stable glasses should be reinterpreted in terms of the dynamical processes that generate stability, rather than the optimized physical quantities they target.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that common computational strategies for ultrastable glasses achieve enhanced stability primarily through dynamical changes in particle diameters rather than through the structural properties (such as large-scale hyperuniformity or small-scale local ordering) that those strategies optimize. To test this, the authors introduce new optimization algorithms that target the same structural quantities without allowing diameter changes. These yield glass configurations with substantially improved hyperuniformity and ordering relative to equilibrium bulk states, yet no corresponding gain in stability. The authors conclude that the targeted physical quantities are correlated with but not causally responsible for ultrastability, and that design rules should be reinterpreted in terms of the underlying dynamical processes.
Significance. If the central negative result holds under rigorous protocol controls, the work would meaningfully redirect research on amorphous materials by decoupling structural optimization from the dynamical mechanisms that actually confer stability. It supplies a concrete counter-example to property-targeted design heuristics and introduces diameter-fixed optimization techniques that could be reused. The manuscript is credited for its clear statement of a falsifiable distinction between correlation and causation and for the explicit construction of the new optimization routes.
major comments (2)
- [§3] §3 (new optimization methods): The manuscript does not demonstrate that the only controlled variable between the diameter-fixed optimizations and the reference diameter-changing algorithms is the presence or absence of diameter dynamics. In particular, it is not shown that both families of configurations are generated with equivalent effective cooling rates, sampled from the same ensemble of inherent structures, or evaluated with identical post-preparation stability protocols (e.g., the precise definition of fictive temperature or relaxation-time measurement). This comparison is load-bearing for the claim that the absence of stability enhancement is due to the lack of diameter changes rather than protocol mismatch.
- [Results] Results section (comparison of stability metrics): No quantitative data, error bars, or statistical tests are provided for the stability observables of the new diameter-fixed glasses versus the diameter-dynamic reference glasses. The central claim that the optimized configurations “do not show enhanced stability” therefore rests on an unverified computational outcome, undermining the strength of the negative result.
minor comments (2)
- [Abstract] The abstract states the negative result without any numerical values or error estimates; adding a single quantitative sentence would improve clarity for readers.
- Notation for the stability metric (e.g., whether it is potential energy per particle, fictive temperature, or relaxation time) should be defined once at first use and used consistently thereafter.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation of our work's significance and for the constructive major comments, which help clarify the robustness of our central negative result. We have revised the manuscript to address both points by adding explicit protocol comparisons and quantitative stability data with statistical support. Our responses to the major comments follow.
read point-by-point responses
-
Referee: [§3] §3 (new optimization methods): The manuscript does not demonstrate that the only controlled variable between the diameter-fixed optimizations and the reference diameter-changing algorithms is the presence or absence of diameter dynamics. In particular, it is not shown that both families of configurations are generated with equivalent effective cooling rates, sampled from the same ensemble of inherent structures, or evaluated with identical post-preparation stability protocols (e.g., the precise definition of fictive temperature or relaxation-time measurement). This comparison is load-bearing for the claim that the absence of stability enhancement is due to the lack of diameter changes rather than protocol mismatch.
Authors: We agree that explicit verification of protocol equivalence is essential to isolate the effect of diameter dynamics. Although the original implementations shared the same underlying Monte Carlo framework, cooling schedules, and inherent-structure sampling procedures, these equivalences were not stated with sufficient detail. In the revised §3 we now include a dedicated paragraph and accompanying table that document identical effective cooling rates (defined via the same temperature-step protocol and acceptance criteria), sampling from the same ensemble of inherent structures (via identical energy-minimization tolerances and quench rates), and identical post-preparation stability protocols (same definitions of fictive temperature via the intersection method and relaxation-time extraction from mean-squared displacement). A new supplementary figure directly compares these control quantities across the diameter-fixed and diameter-dynamic families, confirming no detectable mismatch. revision: yes
-
Referee: [Results] Results section (comparison of stability metrics): No quantitative data, error bars, or statistical tests are provided for the stability observables of the new diameter-fixed glasses versus the diameter-dynamic reference glasses. The central claim that the optimized configurations “do not show enhanced stability” therefore rests on an unverified computational outcome, undermining the strength of the negative result.
Authors: We acknowledge that the initial Results section presented stability comparisons only qualitatively. The revised manuscript now includes a new table and updated figures that report the quantitative stability observables (fictive temperatures and structural relaxation times) for the diameter-fixed optimized glasses, the diameter-dynamic reference glasses, and the equilibrium bulk states. All values are accompanied by standard errors obtained from at least ten independent realizations, and we have added two-sample t-tests (with p-values) demonstrating that the stability metrics of the diameter-fixed glasses are statistically indistinguishable from the bulk equilibrium values while those of the diameter-dynamic glasses are significantly lower. These additions make the negative result quantitatively rigorous. revision: yes
Circularity Check
No circularity: central claims rest on new independent simulations
full rationale
The paper's derivation proceeds by introducing new diameter-fixed optimization protocols, applying them to achieve extreme hyperuniformity and local ordering, and then directly comparing the resulting glass stabilities (via potential energy, fictive temperature, and relaxation times) against equilibrium reference states prepared under identical post-processing. These comparisons are performed via fresh simulations rather than any parameter fit, self-referential definition, or load-bearing prior result from the same authors. No step reduces a claimed prediction to an input by algebraic construction or by renaming a fitted quantity; the argument therefore remains externally falsifiable and self-contained.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Glass stability can be quantified and compared across preparation protocols using conventional metrics such as inherent-structure energy or relaxation time.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We challenge the idea that these physical quantities are causally responsible for glass stability and suggest instead that diameter dynamics is the principal source of enhanced stability.
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
optimizing hyperuniformity at large scale and local ordering at small scale... these glasses do not show enhanced stability
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]
16 FIG. 3. Equation of state in the 1/ϕ(free volume) versus 1/Z∼T /Prepresentation analogous to energy (or free vol- ume) versus temperature evolution for soft potentials. The equilibrium equation of state for the fluid ends near the computational glass transition (vertical dashed line). In the glass phase at larger pressure, the equation of state becomes...
- [2]
-
[3]
17 FIG. 5. Glass equations of state for glasses prepared using ei- ther conventional Monte Carlo methods (upwards triangles) or hyperuniform glasses prepared using random organization dynamics (downwards triangles), for three values ofZinit. Hy- peruniform and conventional glasses are equally stable. are many orders of magnitude larger than those we can s...
work page 2000
- [4]
-
[5]
17 FIG. 7. Glass equations of state for glasses prepared using ei- ther conventional Monte Carlo methods (upwards triangles) or optimized glasses prepared using biased Monte Carlo sim- ulations (downwards triangles), for three values ofZ init. Op- timized and conventional glasses are equally stable. ues used in Fig. 6(a), the reweighted equilibrium densit...
-
[6]
L. Berthier and M. D. Ediger, Facets of glass physics, Physics today69, 40 (2016)
work page 2016
-
[7]
C. A. Angell, Formation of glasses from liquids and biopolymers, Science267, 1924 (1995)
work page 1924
-
[8]
C. P. Royall, F. Turci, S. Tatsumi, J. Russo, and J. Robinson, The race to the bottom: approaching the ideal glass?, Journal of Physics: Condensed Matter30, 363001 (2018)
work page 2018
-
[9]
S. F. Swallen, K. L. Kearns, M. K. Mapes, Y. S. Kim, R. J. McMahon, T. Wu, L. Yu, and S. Satija, Organic glasses with exceptional thermodynamic and kinetic sta- bility, Science315, 353 (2007)
work page 2007
-
[10]
M. D. Ediger, Perspective: Highly stable vapor-deposited glasses, The Journal of Chemical Physics147, 210901 (2017)
work page 2017
-
[11]
C. Rodriguez-Tinoco, M. Gonzalez-Silveira, M. A. Ramos, and J. Rodriguez-Viejo, Ultrastable glasses: new perspectives for an old problem, La Rivista del Nuovo Cimento45, 325 (2022)
work page 2022
-
[12]
L. Berthier and D. R. Reichman, Modern computational studies of the glass transition, Nature Reviews Physics 5, 102 (2023)
work page 2023
-
[13]
J.-L. Barrat and L. Berthier, Computer simulations of the glass transition and glassy materials, Comptes Rendus. Physique24, 57 (2023)
work page 2023
-
[14]
Computational Methods towards Ultrastable Glasses
F. Leoni, M. Ozawa, J. Russo, T. Yanagishima, and A. Ninarello, Computational methods towards ultra- stable glasses, arXiv preprint arXiv:2605.02679 (2026)
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [15]
-
[16]
L. Berthier, P. Charbonneau, E. Flenner, and F. Zam- poni, Origin of ultrastability in vapor-deposited glasses, Phys. Rev. Lett.119, 188002 (2017)
work page 2017
- [17]
-
[18]
T. Yanagishima, J. Russo, R. P. A. Dullens, and H. Tanaka, Towards glasses with permanent stability, Phys. Rev. Lett.127, 215501 (2021)
work page 2021
-
[19]
V. F. Hagh, S. R. Nagel, A. J. Liu, M. L. Manning, and E. I. Corwin, Transient learning degrees of freedom for introducing function in materials, Proceedings of the Na- tional Academy of Sciences119, e2117622119 (2022)
work page 2022
-
[20]
J. R. Dale, J. D. Sartor, R. C. Dennis, and E. I. Corwin, Hyperuniform jammed sphere packings have anomalous material properties, Phys. Rev. E106, 024903 (2022)
work page 2022
-
[21]
X. Fan, D. Xu, J. Zhang, H. Hu, P. Tan, N. Xu, H. Tanaka, and H. Tong, Ideal non-crystals as a distinct form of ordered states without symmetry breaking, Na- ture Materials , 1 (2026)
work page 2026
-
[22]
Y. Wang, Z. Qian, H. Tong, and H. Tanaka, Hyperuni- form disordered solids with crystal-like stability, Nature 11 Communications16, 1398 (2025)
work page 2025
- [23]
-
[24]
V. M. Bolton-Lum, R. C. Dennis, P. K. Morse, and E. I. Corwin, Ideal glass and ideal disk packing in two dimen- sions, Phys. Rev. Lett.136, 058201 (2026)
work page 2026
-
[25]
G. Kapteijns, W. Ji, C. Brito, M. Wyart, and E. Lerner, Fast generation of ultrastable computer glasses by min- imization of an augmented potential energy, Phys. Rev. E99, 012106 (2019)
work page 2019
-
[26]
L. Berthier and M. Ediger, Designing disordered mate- rials beyond equilibrium: Glasses, Nature Materials , 1 (2026)
work page 2026
-
[27]
Torquato, Hyperuniform states of matter, Physics Re- ports745, 1 (2018)
S. Torquato, Hyperuniform states of matter, Physics Re- ports745, 1 (2018)
work page 2018
-
[28]
H. Tong and H. Tanaka, Revealing hidden structural or- der controlling both fast and slow glassy dynamics in supercooled liquids, Phys. Rev. X8, 011041 (2018)
work page 2018
-
[29]
D. Frenkel and B. Smit,Understanding Molecular Sim- ulation: From Algorithms to Applications, third edition ed. (Elsevier Academic Press, S.l., 2023)
work page 2023
-
[30]
T. S. Grigera and G. Parisi, Fast monte carlo algorithm for supercooled soft spheres, Phys. Rev. E63, 045102 (2001)
work page 2001
-
[31]
L. Berthier, D. Coslovich, A. Ninarello, and M. Ozawa, Equilibrium sampling of hard spheres up to the jamming density and beyond, Phys. Rev. Lett.116, 238002 (2016)
work page 2016
-
[32]
A. Ninarello, L. Berthier, and D. Coslovich, Models and algorithms for the next generation of glass transition studies, Phys. Rev. X7, 021039 (2017)
work page 2017
-
[33]
L. Berthier and W. Kob, The monte carlo dynamics of a binary lennard-jones glass-forming mixture, Journal of Physics: Condensed Matter19, 205130 (2007)
work page 2007
-
[34]
L. Berthier and T. A. Witten, Glass transition of dense fluids of hard and compressible spheres, Phys. Rev. E80, 021502 (2009)
work page 2009
-
[35]
Y. S. Elmatad, D. Chandler, and J. P. Garrahan, Cor- responding states of structural glass formers, The Jour- nal of Physical Chemistry B113, 5563 (2009), pMID: 19254014
work page 2009
-
[36]
L. Berthier and M. D. Ediger, How to “measure” a struc- tural relaxation time that is too long to be measured?, The Journal of Chemical Physics153, 044501 (2020)
work page 2020
-
[37]
L. Berthier, P. Chaudhuri, C. Coulais, O. Dauchot, and P. Sollich, Suppressed compressibility at large scale in jammed packings of size-disperse spheres, Phys. Rev. Lett.106, 120601 (2011)
work page 2011
-
[38]
A. B. Bhatia and D. E. Thornton, Structural aspects of the electrical resistivity of binary alloys, Phys. Rev. B2, 3004 (1970)
work page 1970
-
[39]
A. Ikeda and L. Berthier, Thermal fluctuations, mechan- ical response, and hyperuniformity in jammed solids, Phys. Rev. E92, 012309 (2015)
work page 2015
-
[40]
L. Berthier and G. Biroli, Theoretical perspective on the glass transition and amorphous materials, Rev. Mod. Phys.83, 587 (2011)
work page 2011
-
[41]
C. J. Fullerton and L. Berthier, Density controls the ki- netic stability of ultrastable glasses, Europhysics Letters 119, 36003 (2017)
work page 2017
- [42]
-
[43]
D. Hexner and D. Levine, Noise, diffusion, and hyperuni- formity, Phys. Rev. Lett.118, 020601 (2017)
work page 2017
- [44]
- [45]
-
[46]
L. Galliano, M. E. Cates, and L. Berthier, Two- dimensional crystals far from equilibrium, Phys. Rev. Lett.131, 047101 (2023)
work page 2023
-
[47]
L. Galliano and L. Berthier, Glass and jamming tran- sitions in a random organization model, arXiv preprint arXiv:2603.15519 (2026)
-
[48]
D. Landau and K. Binder,A guide to Monte Carlo simu- lations in statistical physics(Cambridge university press, 2021)
work page 2021
-
[49]
A. M. Ferrenberg and R. H. Swendsen, Optimized monte carlo data analysis, Phys. Rev. Lett.63, 1195 (1989)
work page 1989
-
[50]
L. Berthier, E. Flenner, C. J. Fullerton, C. Scalliet, and M. Singh, Efficient swap algorithms for molecular dynam- ics simulations of equilibrium supercooled liquids, Jour- nal of Statistical Mechanics: Theory and Experiment 2019, 064004 (2019)
work page 2019
-
[51]
Szamel, Theory for the dynamics of glassy mixtures with particle size swaps, Phys
G. Szamel, Theory for the dynamics of glassy mixtures with particle size swaps, Phys. Rev. E98, 050601(R) (2018)
work page 2018
- [52]
-
[53]
Y. Nishikawa, F. Ghimenti, L. Berthier, and F. van Wij- land, Irreversible swap algorithms for soft sphere glasses, Phys. Rev. E111, 045416 (2025)
work page 2025
-
[54]
F. Ghimenti, L. Berthier, and F. van Wijland, Irreversible monte carlo algorithms for hard disk glasses: From event- chain to collective swaps, Phys. Rev. Lett.133, 028202 (2024)
work page 2024
-
[55]
L. Berthier, F. Ghimenti, and F. van Wijland, Monte carlo simulations of glass-forming liquids beyond metropolis, The Journal of Chemical Physics161, 114105 (2024)
work page 2024
-
[56]
S. S. Dalal, D. M. Walters, I. Lyubimov, J. J. de Pablo, and M. Ediger, Tunable molecular orientation and ele- vated thermal stability of vapor-deposited organic semi- conductors, Proceedings of the National Academy of Sci- ences112, 4227 (2015)
work page 2015
-
[57]
L. Zhu, C. W. Brian, S. F. Swallen, P. T. Straus, M. D. Ediger, and L. Yu, Surface self-diffusion of an organic glass, Phys. Rev. Lett.106, 256103 (2011)
work page 2011
- [58]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.