pith. machine review for the scientific record. sign in

arxiv: 2605.12127 · v1 · submitted 2026-05-12 · ❄️ cond-mat.stat-mech · cond-mat.dis-nn· cond-mat.soft· physics.comp-ph

Recognition: 2 theorem links

· Lean Theorem

Identifying the relevant parameters in design strategies for stable glasses

Authors on Pith no claims yet

Pith reviewed 2026-05-13 03:34 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech cond-mat.dis-nncond-mat.softphysics.comp-ph
keywords glass stabilityultrastable glasseshyperuniformityparticle diameter dynamicsglass design strategiesoptimization methodssupercooled liquids
0
0 comments X

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.

The paper examines what controls the creation of highly stable glasses beyond ordinary cooling. Existing algorithms improve stability by dynamically varying particle diameters while targeting properties such as large-scale hyperuniformity and small-scale local order. The authors develop alternative optimization procedures that achieve the same property improvements but keep all particle diameters fixed. These diameter-free glasses reach extreme property values yet show no gain in stability relative to ordinary equilibrium glasses. The finding indicates that the dynamical process of diameter variation itself, rather than the final optimized properties, is what generates ultrastability.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2605.12127 by Leonardo Galliano, Ludovic Berthier.

Figure 1
Figure 1. Figure 1: FIG. 1. Equilibrium slow dynamics of the bulk model. (a) [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Static bulk properties in thermal equilibrium. (a) [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Optimizing hyperuniformity at large length scales. [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Glass equations of state for glasses prepared using ei [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [§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.
  2. [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)
  1. [Abstract] The abstract states the negative result without any numerical values or error estimates; adding a single quantitative sentence would improve clarity for readers.
  2. 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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

The work relies on standard domain assumptions in statistical mechanics of glasses and on the validity of computational optimization techniques; no free parameters or invented entities are introduced in the abstract.

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.
    Invoked when claiming that the optimized configurations show no enhanced stability.

pith-pipeline@v0.9.0 · 5524 in / 1123 out tokens · 90841 ms · 2026-05-13T03:34:30.407207+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

58 extracted references · 58 canonical work pages · 1 internal anchor

  1. [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. [2]

    8 1/Z 1/φ

    8 Zinit = 16. 8 1/Z 1/φ

  3. [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...

  4. [4]

    3 1/Z 1/φ

    8 Zinit = 14. 3 1/Z 1/φ

  5. [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. [6]

    Berthier and M

    L. Berthier and M. D. Ediger, Facets of glass physics, Physics today69, 40 (2016)

  7. [7]

    C. A. Angell, Formation of glasses from liquids and biopolymers, Science267, 1924 (1995)

  8. [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)

  9. [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)

  10. [10]

    M. D. Ediger, Perspective: Highly stable vapor-deposited glasses, The Journal of Chemical Physics147, 210901 (2017)

  11. [11]

    Rodriguez-Tinoco, M

    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)

  12. [12]

    Berthier and D

    L. Berthier and D. R. Reichman, Modern computational studies of the glass transition, Nature Reviews Physics 5, 102 (2023)

  13. [13]

    Barrat and L

    J.-L. Barrat and L. Berthier, Computer simulations of the glass transition and glassy materials, Comptes Rendus. Physique24, 57 (2023)

  14. [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)

  15. [15]

    Zhang, F

    G. Zhang, F. H. Stillinger, and S. Torquato, The perfect glass paradigm: Disordered hyperuniform glasses down to absolute zero, Scientific reports6, 36963 (2016)

  16. [16]

    Berthier, P

    L. Berthier, P. Charbonneau, E. Flenner, and F. Zam- poni, Origin of ultrastability in vapor-deposited glasses, Phys. Rev. Lett.119, 188002 (2017)

  17. [17]

    Brito, E

    C. Brito, E. Lerner, and M. Wyart, Theory for swap ac- celeration near the glass and jamming transitions for con- tinuously polydisperse particles, Phys. Rev. X8, 031050 (2018)

  18. [18]

    Yanagishima, J

    T. Yanagishima, J. Russo, R. P. A. Dullens, and H. Tanaka, Towards glasses with permanent stability, Phys. Rev. Lett.127, 215501 (2021)

  19. [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)

  20. [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)

  21. [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)

  22. [22]

    Y. Wang, Z. Qian, H. Tong, and H. Tanaka, Hyperuni- form disordered solids with crystal-like stability, Nature 11 Communications16, 1398 (2025)

  23. [23]

    Leoni, J

    F. Leoni, J. Russo, F. Sciortino, and T. Yanagishima, Generating ultrastable glasses by homogenizing the local virial stress, Phys. Rev. Lett.134, 128201 (2025)

  24. [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)

  25. [25]

    Kapteijns, W

    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)

  26. [26]

    Berthier and M

    L. Berthier and M. Ediger, Designing disordered mate- rials beyond equilibrium: Glasses, Nature Materials , 1 (2026)

  27. [27]

    Torquato, Hyperuniform states of matter, Physics Re- ports745, 1 (2018)

    S. Torquato, Hyperuniform states of matter, Physics Re- ports745, 1 (2018)

  28. [28]

    Tong and H

    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)

  29. [29]

    Frenkel and B

    D. Frenkel and B. Smit,Understanding Molecular Sim- ulation: From Algorithms to Applications, third edition ed. (Elsevier Academic Press, S.l., 2023)

  30. [30]

    T. S. Grigera and G. Parisi, Fast monte carlo algorithm for supercooled soft spheres, Phys. Rev. E63, 045102 (2001)

  31. [31]

    Berthier, D

    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)

  32. [32]

    Ninarello, L

    A. Ninarello, L. Berthier, and D. Coslovich, Models and algorithms for the next generation of glass transition studies, Phys. Rev. X7, 021039 (2017)

  33. [33]

    Berthier and W

    L. Berthier and W. Kob, The monte carlo dynamics of a binary lennard-jones glass-forming mixture, Journal of Physics: Condensed Matter19, 205130 (2007)

  34. [34]

    Berthier and T

    L. Berthier and T. A. Witten, Glass transition of dense fluids of hard and compressible spheres, Phys. Rev. E80, 021502 (2009)

  35. [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

  36. [36]

    Berthier and M

    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)

  37. [37]

    Berthier, P

    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)

  38. [38]

    A. B. Bhatia and D. E. Thornton, Structural aspects of the electrical resistivity of binary alloys, Phys. Rev. B2, 3004 (1970)

  39. [39]

    Ikeda and L

    A. Ikeda and L. Berthier, Thermal fluctuations, mechan- ical response, and hyperuniformity in jammed solids, Phys. Rev. E92, 012309 (2015)

  40. [40]

    Berthier and G

    L. Berthier and G. Biroli, Theoretical perspective on the glass transition and amorphous materials, Rev. Mod. Phys.83, 587 (2011)

  41. [41]

    C. J. Fullerton and L. Berthier, Density controls the ki- netic stability of ultrastable glasses, Europhysics Letters 119, 36003 (2017)

  42. [42]

    Corte, P

    L. Corte, P. M. Chaikin, J. P. Gollub, and D. J. Pine, Random organization in periodically driven systems, Na- ture Physics4, 420 (2008)

  43. [43]

    Hexner and D

    D. Hexner and D. Levine, Noise, diffusion, and hyperuni- formity, Phys. Rev. Lett.118, 020601 (2017)

  44. [44]

    Hexner, P

    D. Hexner, P. M. Chaikin, and D. Levine, Enhanced hy- peruniformity from random reorganization, Proceedings of the National Academy of Sciences114, 4294 (2017)

  45. [45]

    Wilken, R

    S. Wilken, R. E. Guerra, D. Levine, and P. M. Chaikin, Random close packing as a dynamical phase transition, Phys. Rev. Lett.127, 038002 (2021)

  46. [46]

    Galliano, M

    L. Galliano, M. E. Cates, and L. Berthier, Two- dimensional crystals far from equilibrium, Phys. Rev. Lett.131, 047101 (2023)

  47. [47]

    Galliano and L

    L. Galliano and L. Berthier, Glass and jamming tran- sitions in a random organization model, arXiv preprint arXiv:2603.15519 (2026)

  48. [48]

    Landau and K

    D. Landau and K. Binder,A guide to Monte Carlo simu- lations in statistical physics(Cambridge university press, 2021)

  49. [49]

    A. M. Ferrenberg and R. H. Swendsen, Optimized monte carlo data analysis, Phys. Rev. Lett.63, 1195 (1989)

  50. [50]

    Berthier, E

    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)

  51. [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)

  52. [52]

    Ikeda, F

    H. Ikeda, F. Zamponi, and A. Ikeda, Mean field theory of the swap monte carlo algorithm, The Journal of Chemical Physics147, 234506 (2017)

  53. [53]

    Nishikawa, F

    Y. Nishikawa, F. Ghimenti, L. Berthier, and F. van Wij- land, Irreversible swap algorithms for soft sphere glasses, Phys. Rev. E111, 045416 (2025)

  54. [54]

    Ghimenti, L

    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)

  55. [55]

    Berthier, F

    L. Berthier, F. Ghimenti, and F. van Wijland, Monte carlo simulations of glass-forming liquids beyond metropolis, The Journal of Chemical Physics161, 114105 (2024)

  56. [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)

  57. [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)

  58. [58]

    Daley, Z

    C. Daley, Z. Fakhraai, M. Ediger, and J. Forrest, Com- paring surface and bulk flow of a molecular glass former, Soft Matter8, 2206 (2012)