REVIEW 4 major objections 4 minor 47 references
Jamming Energy Landscape is Hierarchical and Ultrametric
T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Jammed packings have a hierarchical energy landscape that becomes precisely ultrametric as system size grows.
desk verdict A clean numerical study that likely shows real hierarchical structure in the jamming landscape, but the headline claim of 'precisely ultrametric' in the thermodynamic limit is stronger than the data support. 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 key object is the subdominant ultrametric $d_<$, the closest ultrametric to a given metric, which the paper obtains by computing the minimum spanning tree of the pairwise-distance graph and setting $d_<(a,b)$ to the largest edge weight on the tree path from $a$ to $b$. The underlying metric $d$ counts differences in the stable contact vector network, so that $d(a,b)\approx\sqrt{\text{(number of changed contacts)}}$ for distances below $\sqrt{N}$. Comparing $d$ with $d_<$ through $D$ measures how close the landscape is to being ultrametric, and the scaling of $D$ with $N$ and pressure is what carries the argument.
What would settle it
Run the same protocol with a perturbation cutoff an order of magnitude larger, $\varepsilon_{\max} = 4/\sqrt{N}$, and compare the collapse of $D(N,p)$. If the plateau value of $D$ changes or no longer collapses onto the same master curve, the ultrametricity is an artifact of the chosen sampling radius rather than a structural property of the landscape. A complementary check: compute $D$ for a crystalline or very high-pressure packing, where the Gardner phase should be absent; if $D$ also vanishes, the measure does not distinguish marginal from non-marginal systems.
Extended reading notes
Core claim
The central claim is that the energy landscape of jammed soft-sphere packings is asymptotically ultrametric: for any three nearby minima $a$, $b$, $c$, the contact-based distance satisfies $d(a,c) \le \max\{d(a,b), d(b,c)\}$ in the limit of infinite system size. The paper demonstrates this by constructing the distance metric from stable contact vector networks for sets of 500–5000 minima found by perturbing a single packing and re-minimizing, and comparing this metric to its subdominant ultrametric $d_<$, built from a minimum spanning tree. The average deviation $D = \sqrt{\langle (d-d_<)^2\rangle}$ collapses onto a master curve as a function of $N^2p$ and approaches a plateau near 2.7, while typical distances grow as $\sqrt{N}$; the fractional deviation therefore vanishes as $N\to\infty$. The authors take this as direct evidence for a marginal Gardner phase along the zero-temperature jamming line, arising from geometry alone rather than thermal fluctuations.
Load-bearing premise
The load-bearing premise is that the set of nearby minima obtained by perturbing a single arbitrary packing with lengths up to $0.4/\sqrt{N}$ fairly represents the local energy landscape, so that the observed ultrametricity reflects the landscape rather than the sampling method.
Editorial extensions
If this is right
- The Gardner phase predicted by mean-field replica theory is not restricted to high dimensions: it appears in ordinary three-dimensional athermal packings.
- Near jamming, the number of nearby minima proliferates and the region that can be densely sampled shrinks, consistent with a fractal basin structure.
- The contact-network metric gives experimentalists a way to search for ultrametricity in granular or colloidal packings without any thermal sampling.
- The finite-size scaling of $D$ means that larger simulated or experimental packings should show cleaner ultrametric signatures, making the prediction testable at accessible system sizes.
Reading between the lines
- An extension not pursued here: polydisperse packings or Hertzian (non-harmonic) contacts may show the same ultrametricity, but the paper tests only monodisperse harmonic spheres, so the universality claim is still open.
- A stricter validation would vary the perturbation cutoff $\varepsilon_{\max}$ well beyond $0.4/\sqrt{N}$; if $D$ were to change appreciably, the hierarchy would be an artifact of the sampling radius rather than an intrinsic landscape property.
- If the landscape is truly ultrametric, relaxation dynamics at low temperature should show hierarchical barrier crossing: a minimum in one sub-basin must pass through the parent basin before reaching a different sub-basin, a signature that could be sought in particle-tracking experiments.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that the energy landscape of three-dimensional athermal jammed packings is hierarchical and ultrametric, providing direct evidence for a marginal Gardner phase along the zero-temperature jamming line. The authors generate sets of nearby minima by repeatedly perturbing a single initial minimum with random displacements and re-minimizing, then define a metric based on differences of stable contact vectors. They construct the subdominant ultrametric from this metric via a minimum spanning tree and quantify the discrepancy with a quantity D. They report that D collapses as a function of N^2p and reaches a plateau of about 2.7 as the scaled pressure goes to zero. Since pairwise distances between minima scale as sqrt(N), they conclude that the fractional excess vanishes in the thermodynamic limit and that the metric becomes 'precisely ultrametric' for all pressures explored.
Significance. If established, this result would be a significant step: it would connect mean-field full replica symmetry breaking (Gardner) predictions to finite-dimensional athermal packings, where no thermal exploration is available. The paper's strengths include the high-precision minimization (force threshold 10^-20), the use of stable contact vectors to avoid rattler ambiguities, and the reported collapse of D across system sizes. However, the current evidence is not yet conclusive because the sampling protocol uses a single parent minimum per state point and lacks any null model or control comparison, so the measured D has no baseline meaning. The central observation is interesting and potentially important, but additional controls are needed before the strong 'precisely ultrametric' claim can be accepted.
major comments (4)
- [Figure 5 and Conclusions] The claim that the metric 'becomes precisely ultrametric' in the thermodynamic limit is stronger than the measurement. Figure 5 shows D reaching a plateau of about 2.7 as N^2p goes to zero, and this plateau is independent of N, while typical pairwise distances grow as sqrt(N). Since D itself does not decrease with N, the absolute violation of the ultrametric inequality does not vanish; only the relative violation D/d tends to zero. To support 'precisely ultrametric', the authors would need to show that D itself tends to zero, or reformulate the conclusion as asymptotic approximate ultrametricity in a relative sense.
- [Sampling procedure, Figures 4 and 5] All nearby minima are generated by perturbing a single arbitrary initial minimum, and no averaging over independent parent packings is reported. The subdominant ultrametric and D are computed from that one local sample, so the results could reflect the structure of only a single basin rather than the landscape as a whole. The conclusion that the jamming energy landscape is hierarchical and ultrametric requires parent-to-parent averaging and a demonstration that D is independent of the chosen parent minimum.
- [Figure 2 and definition of epsilon_max] The perturbation cutoff epsilon_max = 0.4/sqrt(N) is selected from the observed distance-versus-perturbation curves in Figure 2, which are obtained from the same systems that are later analyzed. This couples the sampling protocol to the measured outcome and leaves open the possibility that a different cutoff would change the value of D or the collapse. The authors should test the robustness of the plateau and the collapse to the choice of epsilon_max, or choose the cutoff from independent systems.
- [Equations (2)-(3) and null model] No null model is provided for the quantity D. Because D compares the measured metric to the subdominant ultrametric constructed from the same data, a value of D about 2.7 has no baseline meaning without comparison to, for example, random distance matrices or samples of uncorrelated minima. Such a control is necessary to establish that the observed hierarchy is specific to jammed systems and not a generic property of any distance matrix processed through the same algorithm.
minor comments (4)
- [Equation (1) text] In the sentence containing Equation (1), 'the distanced between' contains a typo; it should read 'the distance d between'.
- [Figure 3 caption] The phrase 'labelled by square rooted numbers' is awkward; consider rewriting as 'the color scale labels are square roots of the metric distance' or similar.
- [Notation for epsilon_max] The perturbation epsilon is first defined with normalization by N^{-1/3}, but the cutoff is later written as epsilon_max = 0.4/sqrt(N). Please clarify the relationship between these normalizations, since the figures appear to use the sqrt(N) convention.
- [References [31] and [32]] The author names in references [31] and [32] contain garbled accented characters in the preprint text; the final version should use properly encoded names.
Circularity Check
No significant circularity: the ultrametricity measurement is an internal diagnostic, not a fitted prediction, and self-citations are not load-bearing.
full rationale
The paper's central quantitative step is the computation of D in Eqs. (2)-(3), which compares the measured distance metric d to the subdominant ultrametric d< constructed from that same metric via a minimum spanning tree. This is a standard, self-contained diagnostic of ultrametricity: d< is by definition the closest ultrametric to d, so D measures how far the empirical metric is from the best ultrametric approximation. This is not circular in the sense of defining the conclusion in terms of the input; the conclusion that the metric becomes 'precisely ultrametric' in the thermodynamic limit is a separate inference from the observed plateau of D and the growth of pairwise distances as sqrt(N). The parameter eps_max = 0.4/sqrt(N) is selected empirically from Figure 2, but it is a sampling cutoff, not a fitted parameter that is later relabeled as a prediction. The measured D value is not forced by that choice. Self-citations (e.g., refs. [24], [27], [28], [33], [34], [36]) are used for simulation methods, scaling relations, and contextual support, not as the load-bearing justification for the ultrametricity claim. The absence of a null model for D is a statistical robustness concern, not a circularity of the kind defined in the reviewing rules. Overall, no step in the derivation reduces to its own inputs by construction, and no load-bearing self-citation chain is present.
Assumptions & free parameters
free parameters (1)
- epsilon_max (perturbation cutoff) =
0.4/sqrt(N)
assumptions (4)
- domain assumption Mean-field replica theory of structural glasses, in which the marginal Gardner phase is characterized by an ultrametric free energy landscape.
- domain assumption The contact vector metric d(a,b) in Equation 2 faithfully represents the relevant configurational distance between energy minima.
- domain assumption The scaling variable N^2 p measures distance to jamming.
- domain assumption Minima found by small random perturbations of one initial minimum and re-minimization are representative of the local energy landscape.
Cite this review
Pith. "Pith review of Jamming Energy Landscape is Hierarchical and Ultrametric." pith.science (2026). https://pith.science/paper/4FY7LB7D
@misc{pith2026190808159,
author = {Pith},
title = {Pith review of: Jamming Energy Landscape is Hierarchical and Ultrametric},
year = {2026},
howpublished = {\url{https://pith.science/paper/4FY7LB7D}},
note = {Machine review of arXiv:1908.08159}
}
read the original abstract
The free energy landscape of mean field marginal glasses is ultrametric. We demonstrate that this feature remains in finite three dimensional systems by finding sets of minima which are nearby in configuration space. By calculating the distance between these nearby minima, we produce a small region of the distance metric. This metric exhibits a clear hierarchical structure and shows the signature of an ultrametric space. That such a hierarchy exists for the jamming energy landscape provides direct evidence for the existence of a marginal phase along the zero temperature jamming line.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
Spin glasses with p-spin interactions,
E. Gardner, “Spin glasses with p-spin interactions,” Nu- clear Physics B 257, 747–765 (1985)
work page 1985
- [2]
-
[3]
Spin-glass theory for pedestrians,
Tommaso Castellani and Andrea Cavagna, “Spin-glass theory for pedestrians,” J. Stat. Mech. 2005, P05012 (2005)
work page 2005
-
[4]
Jorge Kurchan, Giorgio Parisi, Pierfrancesco Urbani, and Francesco Zamponi, “Exact Theory of Dense Amor- phous Hard Spheres in High Dimension. II. The High Density Regime and the Gardner Transition,” J. Phys. Chem. B 117, 12979–12994 (2013)
work page 2013
-
[5]
Frac- tal free energy landscapes in structural glasses,
Patrick Charbonneau, Jorge Kurchan, Giorgio Parisi, Pierfrancesco Urbani, and Francesco Zamponi, “Frac- tal free energy landscapes in structural glasses,” Nature Communications 5, 3725 (2014)
work page 2014
-
[6]
Patrick Charbonneau, Jorge Kurchan, Giorgio Parisi, Pierfrancesco Urbani, and Francesco Zamponi, “Exact theory of dense amorphous hard spheres in high dimen- sion. III. The full replica symmetry breaking solution,” J. Stat. Mech. 2014, P10009 (2014)
work page 2014
-
[7]
Glass and Jamming Transitions: From Exact Results to Finite- Dimensional Descriptions,
Patrick Charbonneau, Jorge Kurchan, Giorgio Parisi, Pierfrancesco Urbani, and Francesco Zamponi, “Glass and Jamming Transitions: From Exact Results to Finite- Dimensional Descriptions,” Annual Review of Condensed Matter Physics 8, 265–288 (2017)
work page 2017
-
[8]
Ergodicity breaking transition in a glassy soft sphere system at small but non-zero temperatures,
Moumita Maiti and Michael Schmiedeberg, “Ergodicity breaking transition in a glassy soft sphere system at small but non-zero temperatures,” Scientific Reports 8, 1837 (2018)
work page 2018
Show all 47 references
-
[9]
Mean-field the- ory of hard sphere glasses and jamming,
Giorgio Parisi and Francesco Zamponi, “Mean-field the- ory of hard sphere glasses and jamming,” Rev. Mod. Phys. 82, 789–845 (2010)
2010
-
[10]
Following the Evolution of Hard Sphere Glasses in Infinite Dimensions under Ex- ternal Perturbations: Compression and Shear Strain,
Corrado Rainone, Pierfrancesco Urbani, Hajime Yoshino, and Francesco Zamponi, “Following the Evolution of Hard Sphere Glasses in Infinite Dimensions under Ex- ternal Perturbations: Compression and Shear Strain,” Phys. Rev. Lett. 114, 015701 (2015)
2015
-
[11]
Following the evolution of glassy states under external perturba- tions: the full replica symmetry breaking solution,
Corrado Rainone and Pierfrancesco Urbani, “Following the evolution of glassy states under external perturba- tions: the full replica symmetry breaking solution,” J. Stat. Mech. 2016, 053302 (2016)
2016
-
[12]
Shear Yielding and Shear Jamming of Dense Hard Sphere Glasses,
Pierfrancesco Urbani and Francesco Zamponi, “Shear Yielding and Shear Jamming of Dense Hard Sphere Glasses,” Phys. Rev. Lett. 118, 038001 (2017)
2017
-
[13]
Liu-Nagel phase diagrams in infinite dimension,
Giulio Biroli and Pierfrancesco Urbani, “Liu-Nagel phase diagrams in infinite dimension,” SciPost Physics 4, 020 (2018)
2018
-
[14]
Marginally stable phases in mean-field structural glasses,
Camille Scalliet, Ludovic Berthier, and Francesco Zam- poni, “Marginally stable phases in mean-field structural glasses,” Phys. Rev. E 99, 012107 (2019)
2019
-
[15]
Breakdown of elasticity in amorphous solids,
Giulio Biroli and Pierfrancesco Urbani, “Breakdown of elasticity in amorphous solids,” Nature Physics12, 1130– 1133 (2016)
2016
-
[16]
Growing timescales and lengthscales characterizing vi- brations of amorphous solids,
Ludovic Berthier, Patrick Charbonneau, Yuliang Jin, Giorgio Parisi, Beatriz Seoane, and Francesco Zamponi, “Growing timescales and lengthscales characterizing vi- brations of amorphous solids,” PNAS 113, 8397–8401 (2016)
2016
-
[17]
Absence of Marginal Stability in a Structural Glass,
Camille Scalliet, Ludovic Berthier, and Francesco Zam- poni, “Absence of Marginal Stability in a Structural Glass,” Phys. Rev. Lett. 119, 205501 (2017)
2017
-
[18]
Spin-glass-like aging in colloidal and granular glasses,
Beatriz Seoane and Francesco Zamponi, “Spin-glass-like aging in colloidal and granular glasses,” Soft Matter 14, 5222–5234 (2018)
2018
-
[19]
Gardner Transition in Physical Dimensions,
C.L. Hicks, M.J. Wheatley, M.J. Godfrey, and M.A. Moore, “Gardner Transition in Physical Dimensions,” Phys. Rev. Lett. 120, 225501 (2018)
2018
-
[20]
A stability-reversibility map uni- fies elasticity, plasticity, yielding, and jamming in hard sphere glasses,
Yuliang Jin, Pierfrancesco Urbani, Francesco Zamponi, and Hajime Yoshino, “A stability-reversibility map uni- fies elasticity, plasticity, yielding, and jamming in hard sphere glasses,” Science Advances 4, eaat6387 (2018)
2018
-
[21]
Hierarchical Land- scape of Hard Disk Glasses,
Qinyi Liao and Ludovic Berthier, “Hierarchical Land- scape of Hard Disk Glasses,” Phys. Rev. X 9, 011049 (2019)
2019
-
[22]
An exploratory study of the glassy landscape near jamming,
Claudia Artiaco, Paolo Baldan, and Giorgio Parisi, “An exploratory study of the glassy landscape near jamming,” arXiv:1908.06127 [cond-mat] (2019), arXiv: 1908.06127
2019
-
[23]
Experimental Evidence of the Gardner Phase in a Granular Glass,
A. Seguin and O. Dauchot, “Experimental Evidence of the Gardner Phase in a Granular Glass,” Phys. Rev. Lett. 117, 228001 (2016)
2016
-
[24]
Seeing through a glass clearly: Experimental observation of the marginal glass phase,
Andrew P. Hammond and Eric I. Corwin, “Seeing through a glass clearly: Experimental observation of the marginal glass phase,” arXiv:1908.08152 [cond-mat] (2019), arXiv: 1908.08152
2019 arXiv
-
[25]
Jamming at zero temperature and zero applied stress: The epitome of disorder,
Corey S. OHern, Leonardo E. Silbert, Andrea J. Liu, and Sidney R. Nagel, “Jamming at zero temperature and zero applied stress: The epitome of disorder,” Phys. Rev. E 6 68, 011306 (2003)
2003
-
[26]
Jamming Transition in Granular Systems,
T. S. Majmudar, M. Sperl, S. Luding, and R. P. Behringer, “Jamming Transition in Granular Systems,” Phys. Rev. Lett. 98, 058001 (2007)
2007
-
[27]
Universal Microstructure and Mechanical Stability of Jammed Packings,
Patrick Charbonneau, Eric I. Corwin, Giorgio Parisi, and Francesco Zamponi, “Universal Microstructure and Mechanical Stability of Jammed Packings,” Phys. Rev. Lett. 109, 205501 (2012)
2012
-
[28]
Geometric signa- tures of jamming in the mechanical vacuum,
Peter K. Morse and Eric I. Corwin, “Geometric signa- tures of jamming in the mechanical vacuum,” Phys. Rev. Lett. 112, 115701 (2014)
2014
-
[29]
Replica Symmetry Breaking in Short-Range Spin Glasses: The- oretical Foundations and Numerical Evidences,
Enzo Marinari, Giorgio Parisi, Federico Ricci-Tersenghi, Juan J. Ruiz-Lorenzo, and Francesco Zuliani, “Replica Symmetry Breaking in Short-Range Spin Glasses: The- oretical Foundations and Numerical Evidences,” Journal of Statistical Physics 98, 973–1074 (2000)
2000
-
[30]
On Ultrametricity, Data Coding, and Computation,
Fionn Murtagh, “On Ultrametricity, Data Coding, and Computation,” Journal of Classification 21, 167–184 (2004)
2004
-
[31]
Ultrametricity property of energy landscapes of multi- disperse packing problems,
Johannes J. Schneider, Andr Mller, and Elmar Schmer, “Ultrametricity property of energy landscapes of multi- disperse packing problems,” Phys Rev E Stat Nonlin Soft Matter Phys 79, 031122 (2009)
2009
-
[32]
Structural Relaxation Made Simple,
Erik Bitzek, Pekka Koskinen, Franz Ghler, Michael Moseler, and Peter Gumbsch, “Structural Relaxation Made Simple,” Physical Review Letters 97 (2006), 10.1103/PhysRevLett.97.170201
2006 doi
-
[33]
Universal Non- Debye Scaling in the Density of States of Amorphous Solids,
Patrick Charbonneau, Eric I. Corwin, Giorgio Parisi, Alexis Poncet, and Francesco Zamponi, “Universal Non- Debye Scaling in the Density of States of Amorphous Solids,” Phys. Rev. Lett. 117, 045503 (2016)
2016
-
[34]
Echoes of the Glass Transition in Athermal Soft Spheres,
Peter K. Morse and Eric I. Corwin, “Echoes of the Glass Transition in Athermal Soft Spheres,” Phys. Rev. Lett. 119, 118003 (2017)
2017
-
[35]
P Allen and author.) Tildesley, D
M. P Allen and author.) Tildesley, D. J, Computer sim- ulation of liquids , second edition ed. (Oxford : Oxford University Press, 2017)
2017
-
[36]
Jamming Criticality Revealed by Removing Localized Buckling Excitations,
Patrick Charbonneau, Eric I. Corwin, Giorgio Parisi, and Francesco Zamponi, “Jamming Criticality Revealed by Removing Localized Buckling Excitations,” Phys. Rev. Lett. 114, 125504 (2015)
2015
-
[37]
Enu- meration of distinct mechanically stable disk packings in small systems,
G.-J. Gao, J. Blawzdziewicz, and C. S. O’Hern, “Enu- meration of distinct mechanically stable disk packings in small systems,” Philosophical Magazine 87, 425–431 (2007)
2007
-
[38]
Calculations of the structure of basin volumes for mechanically stable packings,
S. S. Ashwin, Jerzy Blawzdziewicz, Corey S. O’Hern, and Mark D. Shattuck, “Calculations of the structure of basin volumes for mechanically stable packings,” Phys. Rev. E 85, 061307 (2012)
2012
-
[39]
Direct De- termination of the Size of Basins of Attraction of Jammed Solids,
Ning Xu, Daan Frenkel, and Andrea J. Liu, “Direct De- termination of the Size of Basins of Attraction of Jammed Solids,” Phys. Rev. Lett. 106, 245502 (2011)
2011
-
[40]
Microscopic theory of the jamming transition of harmonic spheres,
Ludovic Berthier, Hugo Jacquin, and Francesco Zam- poni, “Microscopic theory of the jamming transition of harmonic spheres,” Phys. Rev. E 84, 051103 (2011)
2011
-
[41]
On the shortest spanning subtree of a graph and the traveling salesman problem,
Joseph B. Kruskal, “On the shortest spanning subtree of a graph and the traveling salesman problem,” Proc. Amer. Math. Soc. 7, 48–50 (1956)
1956
-
[42]
On the degree of ultrametricity,
R. Rammal, J.C. Angles d’Auriac, and B. Doucot, “On the degree of ultrametricity,” Journal de Physique Lettres 46, 945–952 (1985)
1985
-
[43]
A Survey of Recent Advances in Hierar- chical Clustering Algorithms,
F. Murtagh, “A Survey of Recent Advances in Hierar- chical Clustering Algorithms,” Comput J 26, 354–359 (1983)
1983
-
[44]
Finite-Size Scaling at the Jamming Transition,
Carl P. Goodrich, Andrea J. Liu, and Sidney R. Nagel, “Finite-Size Scaling at the Jamming Transition,” Phys. Rev. Lett. 109, 095704 (2012). 7 Supplementary Materials THE SUBDOMINANT UL TRAMETRIC We present here a simple outline of the algorithm for creating the subdominant ult...
2012
-
[45]
This can be reinterpreted as an edge-weighted graph where the nodes are the minima and the edge weights are the distances between minima
The metric, d, is a symmetric matrix of pairwise distances between minima. This can be reinterpreted as an edge-weighted graph where the nodes are the minima and the edge weights are the distances between minima
-
[46]
The minimum spanning tree is unique [41] and has the property that every pair of nodes has only one path connecting them
We compute the minimum spanning tree of this graph, which is simply the network with the minimum possible total edge weight (sum of distance values) which connects every node into a single tree. The minimum spanning tree is unique [41] and has the property that every pair of n...
-
[47]
The hierarchical nature of d< derives from that of the minimum spanning tree
The subdominant ultrametric, d<, is created as a symmetric matrix with entriesd< ij determined by the maximum edge weight in the path from node i to node j in the minimum spanning tree. The hierarchical nature of d< derives from that of the minimum spanning tree. The maximum c...
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.