REVIEW 3 major objections 5 minor 46 references
Dynamic heterogeneity in the self-induced spin glass state of elemental neodymium
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Slow and fast magnetic regions coexist in neodymium's spin glass.
desk verdict A potentially important first real-space look at dynamic heterogeneity in a spin glass, but the key experimental classification is done by eye and needs a quantitative check before the claim is solid. 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 central objects are the nearly degenerate magnetic wave-vector pockets $Q_A$, $Q_B$, and $Q_C$ read from the static magnetic structure factor $S(Q)$, obtained by Fourier transforming spin-polarized scanning tunneling microscopy images of the frozen magnetization. The experimental protocol is a three-step dynamical cycle: zero-field cool through the Néel and glass transitions to 1.3 K, apply magnetic-field cycles to induce dynamics and re-image at zero field, and warm through the ordered phase to reinitialize. Slow and fast dynamics are mapped by dividing images into $12 \times 12$ nm boxes and sorting each box by eye, using criteria in Supplementary Text S5, into stable, changing, or undetermined categories; the Jensen-Shannon divergence $\mathcal{D}_{JS}$ quantifies the similarity between successive structure factors. The simulations use atomistic spin dynamics based on a Heisenberg Hamiltonian with density-functional-theory-derived exchange interactions, tracking the magnetization over hundreds of picoseconds, and use a subtraction-filtering method to identify regions of change.
What would settle it
Reanalyze the same before-and-after magnetization images with an automated, pre-registered classifier that extracts the dominant Q-states in each box without human judgment and compare the resulting slow/fast maps to the by-eye maps of Fig. 4 and Fig. S6; if the automated maps do not reproduce the coexistence of stable and changing regions across multiple samples and field cycles, the claim of dynamic heterogeneity would fail. A complementary check is to compute the four-point dynamic susceptibility from the image series and look for a growing dynamical correlation length, whose absence would contradict the dynamic-heterogeneity interpretation.
Extended reading notes
Core claim
The central claim is that the self-induced spin glass state of elemental neodymium exhibits dynamic heterogeneity, both in spin-polarized scanning tunneling microscopy experiments and in atomistic spin dynamics simulations. In the frozen glass at 1.3 K, the magnetization is a patchwork of locally ordered multi-Q patterns; after dynamics are triggered by cycling the magnetic field to 2-7 T and the system is re-imaged at zero field, comparison of before-and-after images shows spatial regions whose patterns are unchanged (slow dynamics) beside regions whose patterns have changed (fast dynamics). This coexistence persists across repeated cycles, different field values, and several samples, and the simulations reproduce it, showing regions of change that nucleate at particular locations and then grow. In parallel, the magnetic structure factor $S(Q)$ evolves systematically under field cycling: the $Q_A$ pocket moves radially outward while $Q_B$ and $Q_C$ move inward, and the evolution saturates after a few cycles. Zero-field cooling reproducibly imprints a similar initial $S(Q)$, showing that the freshly cooled glass is an intermediate metastable state that ages under field cycling and can be thermally reinitialized.
Load-bearing premise
The classification of real-space regions into slow and fast dynamics is done by eye using criteria in Supplementary Text S5, and the box size of 12 nm is chosen in S7 because it maximizes the fraction of boxes that can be sorted into one of the two categories; if that visual sorting is biased or coincidental, the central observation of dynamic heterogeneity could be an artifact of the analysis.
Editorial extensions
If this is right
- A correct account of aging in this spin glass must include local length scales, not only the global, length-invariant picture of mean-field spin-glass theory.
- The systematic, sample-independent evolution of $S(Q)$ implies the zero-field-cooled state is an intermediate metastable state that ages toward a different set of periodicities under repeated field cycling.
- Thermal cycling above the glass transition resets the periodicities but not the specific real-space arrangement of patterns, so reinitialization is statistical rather than a recovery of the exact configuration.
- The coexistence of slow and fast magnetic regions provides a concrete experimental bridge between spin-glass aging and dynamic heterogeneity in structural glasses.
- The simulations reproduce spatially heterogeneous dynamics without an accompanying change of $S(Q)$, suggesting that structural-factor aging and heterogeneous real-space dynamics are distinct aspects of the glassy behavior.
Reading between the lines
- A testable extension beyond this paper would be to apply the same imaging protocol to other frustrated but disorder-free magnets: if dynamic heterogeneity is generic to self-induced spin glasses, similar slow/fast coexistence should appear in those systems too.
- The by-eye box classification could be replaced by an automated local-FFT classifier that assigns dominant Q-states per box; if such a classifier reproduces the slow/fast maps across samples, it would put the heterogeneity claim on a more quantitative footing and could yield a four-point-like dynamic correlation length.
- The separation between reproducible periodicities and non-reproducible spatial patterns suggests a distinction between thermodynamic memory and configurational memory; partial warming excursions could reveal which Q-pockets reset first, something the present study does not test.
- If field-history can reliably set and reset specific Q-pockets, repeated cycling might be used as a form of magnetic state writing, analogous to phase-change memory; the paper mentions this pathway but does not demonstrate it.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports spin-polarized scanning tunneling microscopy (SP-STM) measurements and atomistic spin dynamics (ASD) simulations of the self-induced spin glass state of elemental Nd(0001). The central claim is the observation of dynamic heterogeneity (DH): a spatial coexistence of slow and fast evolving magnetic patterns in the frozen spin glass state, reminiscent of DH in structural glasses. The paper also reports that zero-field cooling imprints a preferred set of metastable magnetic periodicities (Q-pockets) that evolve systematically under repeated magnetic-field cycling and can be thermally reinitialized; this evolution is quantified with the Jensen-Shannon divergence between successive magnetic structure factors. The DH claim is supported by qualitative slow/fast classification of real-space boxes by eye (Supplementary Text S5), by a box size chosen to maximize classifiable boxes (S7), and by threshold-based subtraction in the simulations (S10).
Significance. If the DH claim is established quantitatively, this would be a significant real-space observation of dynamic heterogeneity in a spin glass, bridging glassy dynamics in magnetic systems and structural glasses, and providing evidence beyond the mean-field description. The S(Q)-evolution and reinitialization part is more quantitative and rests on multiple samples, 27-hour stability checks, and a well-defined similarity metric; this portion is convincing and provides a solid basis for the metastable-periodicity claim. The simulations reproduce spatially heterogeneous dynamics and glassy autocorrelation functions, and the paper includes a data-availability statement and reproducible methodology via UppASD. The principal weakness is that the experimental DH classification is by eye and lacks a quantitative mobility metric or null-model test, so the core claim is not yet rigorously supported.
major comments (3)
- [Main text, 'Real-space imaging of DH'; Supplementary Text S5 and S7] The central claim of dynamic heterogeneity rests on classifying 12-nm boxes as slow or fast by eye, as explicitly stated in S5 ('Each box was examined by eye'). S7 states that L = 12 nm is selected because it maximizes the fraction of boxes that can be sorted into slow/fast categories. No quantitative local mobility metric (e.g., an overlap function or a local autocorrelation), no null-hypothesis test (e.g., random permutation of box labels or a comparison to a Poisson process), and no repeatability/inter-rater check is provided. The spatially clustered appearance in Fig. 4C,F and Fig. S6G-I could therefore be produced by the sorting procedure itself rather than by a physical property of the spin glass. The remark that a four-point correlation function is difficult to apply does not remove the need for a quantitative demonstration that the slow/fast labels are spatially correlated beyond what a random arrangement would produce.
- [Fig. 5 and Supplementary Text S10] The simulated evidence for DH uses a subtraction method that is defined by multiple hand-set parameters (top 30% of extremal magnetization values, a Gaussian filter, threshold tau_M = 0.25, hole filling, and neglect of the smallest regions). The main text claims that the evolution shown in Fig. 5D is 'neither consistent with stochastic behavior (i.e. random switching) nor with the clear movement of a favorable domain,' but no statistical comparison to a stochastic null model or to a coarsening model is presented. Because the same subtraction defines both the changed regions and their spatial clustering, the simulation analysis is circular unless the observed clustering is compared with, for example, randomly shuffled spin flips or randomly placed changed regions. The ASD snapshots are valuable and show qualitative heterogeneous dynamics, but they do not yet quantitatively establish DH or the nucleation claim.
- [Supplementary Text S7 and Fig. 4] The choice of L = 12 nm is justified only by maximizing the fraction of boxes that can be classified as slow or fast, which conflates the analysis resolution with the physical length scale of the heterogeneous dynamics. A quantitative measure of the spatial correlation of the slow/fast labels, such as a spatial correlation function of the binary map, a cluster-size distribution compared to a random permutation, or a four-point susceptibility computed on the box labels, is needed to establish that the heterogeneity has a characteristic length scale. Without such a measure, the statement that the DH 'persisted' for variable box sizes (as long as L stays smaller than the pattern size) does not distinguish a real dynamical length scale from an artifact of the grid.
minor comments (5)
- [Abstract] The phrase 'ubiquitous behavior of dynamic heterogeneity' is unclear and overbroad: the paper studies a single material, and 'ubiquitous' is not defined. Suggest rewording to 'we demonstrate dynamic heterogeneity in the self-induced spin glass state of elemental neodymium.'
- [Supplementary Text S7 and Fig. S7 caption] The text in S7 says 'for diminishing box sizes L,' while the Fig. S7 caption says 'increasing box sizes L.' Please reconcile the direction of the series.
- [Main text, Fig. 5D description] The main text refers to the 'subtraction method in S9' when describing Fig. 5D, but the subtraction method is described in S10. Please correct the cross-reference.
- [Conclusion and References] The conclusion ends with an orphan citation '(43-46)' that is not attached to a sentence or in-text discussion. These references are used in the Methods and Supplementary, so they should not appear as a dangling citation in the main text.
- [Supplementary Text S10] The threshold tau_M = 0.25 is introduced without a definition of how it is applied to the Gaussian-filtered image. Please specify the exact operation (e.g., pixels above 0.25 times the maximum intensity) so that the method is reproducible.
Circularity Check
No significant circularity: the DH claim is a direct observation from new imaging and simulation data.
full rationale
The paper's central claim, dynamic heterogeneity in the self-induced spin glass state of elemental neodymium, is an experimental and simulational observation rather than a quantity derived from fitted parameters. The slow/fast classification in Fig. 4C,F and Supplementary Text S5 is a direct visual comparison of local periodicities before and after field cycling; the box size L = 12 nm is selected in S7 because it minimizes the fraction of 'undetermined' boxes, but both slow and fast boxes are identified at all tested box sizes, so the coexistence of the two categories is not manufactured by that choice. The S(Q) evolution is quantified with the Jensen-Shannon divergence computed directly from filtered FFT maps (S2), and the filtering threshold n is an analysis choice, not a parameter fitted to force the claimed radial shifts, which are also visible in the raw comparisons of Fig. 2G-H. The ASD simulations use the same DFT-derived Heisenberg Hamiltonian as prior work (Ref. 24) and reproduce glassy autocorrelation behavior (S12); the subtraction-based change maps of S10-S11 are again direct comparisons of simulated snapshots, not outputs of a fitted model. Self-citations to Refs. 20-25 establish the background existence and characterization of the self-induced spin glass state, but the new dynamic-heterogeneity observation does not reduce to those citations or to any equation in this paper. No fitted parameter is renamed as a prediction, and no quantity is defined in terms of the conclusion it supports. The by-eye classification and threshold choices are legitimate robustness concerns, but they are not circularity.
Assumptions & free parameters
free parameters (4)
- background threshold coefficient n =
1-3.5 (adjusted per sample)
- box size L =
12 nm
- extremal magnetization top fraction =
30%
- gaussian filter threshold tau_M =
0.25
assumptions (4)
- domain assumption The Heisenberg Hamiltonian with DFT-derived exchange interactions from Ref 24 (no anisotropy) describes the magnetic behavior of dhcp Nd.
- domain assumption SP-STM images at B=0 and T=1.3 K capture the frozen spin texture without perturbing it, and magnetic field cycles induce the dynamics.
- domain assumption The Nd(0001) surface magnetization represents bulk-like spin glass behavior of the islands.
- domain assumption Visual classification of local patterns into Q-states is a valid way to identify slow and fast dynamics.
Cite this review
Pith. "Pith review of Dynamic heterogeneity in the self-induced spin glass state of elemental neodymium." pith.science (2026). https://pith.science/paper/4WJVHFE6
@misc{pith2026241215916,
author = {Pith},
title = {Pith review of: Dynamic heterogeneity in the self-induced spin glass state of elemental neodymium},
year = {2026},
howpublished = {\url{https://pith.science/paper/4WJVHFE6}},
note = {Machine review of arXiv:2412.15916}
}
read the original abstract
Spin glasses are magnetic materials exhibiting numerous magnetization patterns, that randomly vary both in real space and in time. To date, it is still not well understood what the nature of these spatiotemporal dynamics is, namely if they are completely random or if there are links between given time and length scales. Here we show the ubiquitous behavior of dynamic heterogeneity in the self-induced spin glass state of elemental neodymium. We used spin-polarized scanning tunneling microscopy in combination with atomistic spin dynamics simulations to image the locally ordered magnetic patterns in the glass state, and tracked the induced spatiotemporal dynamics in response to external perturbations. We observed that the real space magnetization exhibited a coexistence of slow and fast dynamics reminiscent of dynamic heterogeneity in structural glasses. Furthermore, we found that zero-field cooling imprints a specific set of metastable periodicities into the spin glass, which evolved during aging and could be thermally reinitialized. These results demonstrate the importance of local length scales for the understanding of aging dynamics in spin glasses and provide a link to the more general picture of true glasses.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
- [1]
-
[2]
K. H. Fischer, J. A. Hertz, Spin glasses. (Cambridge University Press, Cambridge, 1993)
work page 1993
- [3]
- [4]
-
[5]
S. F. Edwards, P. W. Anderson, Theory of Spin Glasses. J. Phys. F Met. Phys. 5, 965- 974 (1975)
work page 1975
-
[6]
S. F. Edwards, P. W. Anderson, Theory of Spin-Glasses. II. J. Phys. F Met. Phys. 6, 1927-1937 (1976)
work page 1976
-
[7]
D. Sherrington, S. Kirkpatrick, Solvable Model of a Spin-Glass. Phys. Rev. Lett. 35, 1792-1796 (1975)
work page 1975
-
[8]
Parisi, Infinite Number of Order Parameters for Spin-Glasses
G. Parisi, Infinite Number of Order Parameters for Spin-Glasses. Phys. Rev. Lett. 43, 1754-1756 (1979)
work page 1979
Show all 46 references
-
[9]
Parisi, Order Parameter for Spin-Glasses
G. Parisi, Order Parameter for Spin-Glasses. Phys. Rev. Lett. 50, 1946-1948 (1983). 11
1983
-
[10]
Mezard, G
M. Mezard, G. Parisi, N. Sourlas, G. Toulous e, M. Virasoro, Nature of the Spin-Glass Phase. Phys. Rev. Lett. 52, 1156-1159 (1984)
1984
-
[11]
Lefloch, J
F. Lefloch, J. Hammann, M. Ocio, E. Vincent, Can Aging Phenomena Discriminate between the Droplet Model and a Hierarchical Description in Spin-Glasses. Europhys. Lett. 18, 647-652 (1992)
1992
-
[12]
Jonason, E
K. Jonason, E. Vincent, J. Hammann, J. P. Bouchaud, P. Nordblad, Memory and chaos effects in spin glasses. Phys. Rev. Lett. 81, 3243-3246 (1998)
1998
-
[13]
Baity-Jesi, E
M. Baity-Jesi, E. Calore, A. Crvz, L. A. Fernandez, J. M. Gil-Narvion, I. G. A. Pemartin, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V . Martin- Mayor, J. Moreno-Gordo, A. M. Sudupe, D. Navarro, I. Paga, G. Parisi, S. Perez- Gaviro, F. Ricci-Tersenghi, J. J. R...
2023
-
[14]
Y . G. Joh, R. Orbach, G. G. Wood, J. Hammann, E. Vincent, Extraction of the spin glass correlation length. Phys. Rev. Lett. 82, 438-441 (1999)
1999
-
[15]
Jonason, P
K. Jonason, P. Nordblad, E. Vincent, J. Hammann, J. P. Bouchard, Memory interference effects in spin glasses. Eur. Phys. J. B 13, 99-105 (2000)
2000
-
[16]
Dupuis, F
V . Dupuis, F. Bert, J. P. Bouchaud, J. Ha mmann, F. Ladieu, D. Parker, E. Vincent, Aging, rejuvenation and memory phenomena in spin glasses. Pramana-J. Phys. 64, 1109-1119 (2005)
2005
-
[17]
Refregier, E
P. Refregier, E. Vincent, J. Hammann, M. Ocio, Aging Phenomena in a Spin-Glass - Effect of Temperature-Changes Below Tg. J. Phys. (Paris) 48, 1533-1539 (1987)
1987
-
[18]
V . S. Dotsenko, Fractal Dynamics of Spin-Glasses. J. Phys. C Solid State 18, 6023- 6031 (1985)
1985
-
[19]
Vincent, in Ageing and the Glass Transition, M
E. Vincent, in Ageing and the Glass Transition, M. Henkel, M. Pleimling, R. Sanctuary, Eds. (Springer Berlin Heidelberg, Berlin, Heidelberg, 2007), pp. 7-60
2007
-
[20]
Principi, M
A. Principi, M. I. Katsnelson, Stripe glasses in ferromagnetic thin films. Phys. Rev. B 93, 054410 (2016)
2016
-
[21]
Principi, M
A. Principi, M. I. Katsnelson, Self-Induced Glassiness and Pattern Formation in Spin Systems Subject to Long-Range Interactions. Phys. Rev. Lett. 117, 137201 (2016)
2016
-
[22]
Mauri, M
A. Mauri, M. I. Katsnelson, Frustrated ma gnets in the limit of infinite dimensions: Dynamics and disorder-free glass transition. Phys. Rev. B 109, 144414 (2024)
2024
-
[23]
Kolmus, M
A. Kolmus, M. I. Kasnelson, A. A. Kh ajetoorians, H. J. Kappen, Atom-by-atom construction of attractors in a tunable finite size spin array. New. J. Phys. 22 , 023038 (2020)
2020
-
[24]
Kamber, A
U. Kamber, A. Bergman, A. Eich, D. Iusan, M. Steinbrecher, N. Hauptmann, L. Nordström, M. I. Katsnelson, D. Wegner, O. Eriksson, A. A. Khajetoorians, Self- induced spin glass state in elemental and crystalline neodymium. Science 368, eaay6757 (2020)
2020
-
[25]
Verlhac, L
B. Verlhac, L. Niggli, A. Bergman, U. Ka mber, A. Bagrov, D. Iusan, L. Nordstrm, M. I. Katsnelson, D. Wegner, O. Eriksson, A. A. Khajetoorians, Thermally induced magnetic order from glassiness in elemental neodymium. Nat. Phys. 18, 905-911 (2022)
2022
-
[26]
Dynamical Heterogeneities in Glasses, Colloids, and Granular Media. L. Berthier, G. Biroli, J.-P. Bouchaud, L. Cipelletti, W. van Saarloos, Eds., (Oxford University Press, 2011)
2011
-
[27]
E. R. Weeks, J. C. Crocker, A. C. Levitt, A. Schofield, D. A. Weitz, Three-dimensional direct imaging of structural relaxation near the colloidal glass transition. Science 287, 627-631 (2000). 12
2000
-
[28]
Dauchot, G
O. Dauchot, G. Marty, G. Biroli, Dyna mical heterogeneity close to the jamming transition in a sheared granular material. Phys. Rev. Lett. 95, (2005)
2005
-
[29]
S. C. Glotzer, Spatially heterogeneous dynamics in liquids: insights from simulation. J. Non-Cryst. Solids 274, 342-355 (2000)
2000
-
[30]
V ollmayr-Lee, W
K. V ollmayr-Lee, W. Kob, K. Binder, A. Zippelius, Dynamical heterogeneities below the glass transition. J. Chem. Phys. 116, 5158-5166 (2002)
2002
-
[31]
W. K. Kegel, A. van Blaaderen, Direct observation of dynamical heterogeneities in colloidal hard-sphere suspensions. Science 287, 290-293 (2000)
2000
-
[32]
S. C. Glotzer, N. Jan, T. Lookman, A. B. MacIsaac, P. H. Poole, Dynamical heterogeneity in the Ising spin glass. Phys. Rev. E 57, 7350-7353 (1998)
1998
-
[33]
H. E. Castillo, C. Chamon, L. F. Cugliandol o, J. L. Iguain, M. P. Kennett, Spatially heterogeneous ages in glassy systems. Phys. Rev. B 68, 134442 (2003)
2003
-
[34]
Belletti, A
F. Belletti, A. Cruz, L. A. Fernandez, A. Gordillo-Guerrero, M. Guidetti, A. Maiorano, F. Mantovani, E. Marinari, V . Martin-Mayor, J. Monforte, A. M. Sudupe, D. Navarro, G. Parisi, S. Perez-Gaviro, J. J. Ruiz-Lorenzo, S. F. Schifano, D. Sciretti, A. Tarancon, R. Tripiccione, ...
2009
-
[35]
S. C. Glotzer, V . N. Novikov, T. B. Sc hroder, Time-dependent, four-point density correlation function description of dynamical heterogeneity and decoupling in supercooled liquids. J. Chem. Phys. 112, 509-512 (2000)
2000
-
[36]
Donati, S
C. Donati, S. Franz, S. C. Glotzer, G. Parisi, Theory of non-linear susceptibility and correlation length in glasses and liquids. J. Non-Cryst. Solids 307, 215-224 (2002)
2002
-
[37]
Berthier, G
L. Berthier, G. Biroli, J. P. Bouchaud, L. Ci pelletti, D. El Masri, D. L'Hôte, F. Ladieu, M. Pierno, Direct experimental evidence of a growing length scale accompanying the glass transition. Science 310, 1797-1800 (2005)
2005
-
[38]
A. S. Keys, A. R. Abate, S. C. Glot zer, D. J. Durian, Measurement of growing dynamical length scales and prediction of the jamming transition in a granular material. Nat. Phys. 3, 260-264 (2007)
2007
-
[39]
Lacevic, F
N. Lacevic, F. W. Starr, T. B. Schroder, S. C. Glotzer, Spatially heterogeneous dynamics investigated via a time-dependent four-point density correlation function. J. Chem. Phys. 119, 7372-7387 (2003)
2003
-
[40]
M. T. Cicerone, M. D. Ediger, Relaxa tion of Spatially Heterogeneous Dynamic Domains in Supercooled Ortho-Terphenyl. J. Chem. Phys. 103, 5684-5692 (1995)
1995
-
[41]
H. E. Castillo, C. Chamon, L. F. Cugliandol o, M. P. Kennett, Heterogeneous aging in spin glasses. Phys. Rev. Lett. 88, 237201 (2002)
2002
-
[42]
Baity-Jesi, E
M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, I. Gonzalez- Adalid Pemartin, A. Gordillo-Guerrero, D. Iniguez, A. Maiorano, E. Marinari, V . Martin-Mayor, J. Moreno-Gordo, A. Munoz Sudupe, D. Navarro, I. Paga, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersen...
2024
-
[43]
J. H. Lin, Divergence Measures Based on the Shannon Entropy. Ieee. T. Inform. Theory. 37, 145-151 (1991)
1991
-
[44]
D. M. Endres, J. E. Schindelin, A ne w metric for probability distributions. Ieee. T. Inform. Theory. 49, 1858-1860 (2003)
2003
-
[45]
https://github.com/UppASD/UppASD, (2020)
Uppsala Atomic Spin Dy namics (UppASD) software. https://github.com/UppASD/UppASD, (2020)
2020
-
[46]
What can we ‘learn’ with atoms?
O. Eriksson, A. Bergman, L. Bergqvist, J. Hellsvik, Atomistic Spin Dynamics: Foundations and Applications. (Oxford University Press, 2017). Acknowledgments: We thank O. Eriksson and A. Mauri for fruitful discussions. 13 Funding: This project was supported by the European Resea...
2017
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.