REVIEW 3 major objections 4 minor 1 cited by
Quenched Dynamics of Artificial Spin Ice: Coarsening versus Kibble-Zurek
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read For deep quenches of square ice, defect density falls as the inverse square root of quench time, and the mechanism is coarsening, not Kibble-Zurek.
desk verdict Clean simulation result, but the case for coarsening over Kibble-Zurek rests on an unverified universality class for the colloidal model. 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 vertex types of square ice, labeled by the number of colloids pointing toward a vertex: $N_0,N_4$ double monopoles, $N_1,N_3$ monopoles, biased $N_{2,\mathrm{biased}}$, and ground-state $N_{2,\mathrm{gs}}$; non-ground-state vertices are the defects whose density is measured. The argument runs on scaling comparisons: the measured defect-density exponent $\beta \simeq 1/2$ and time-rescaling exponent $\alpha=3/4$ are checked against the Kibble-Zurek formulas $\beta = D\nu/(1+z\nu)$ and $\alpha = z\nu/(1+z\nu)$ for the 2D Ising class ($\nu=1$, $z=2$, $D=2$), each predicting $2/3$. The coarsening alternative uses the standard growth law $R(t)\sim t^{1/z}$ for ordered regions to say that defect density should fall as $1/R(t)$, giving $\rho_d\sim \tau_Q^{-1/2}$.
What would settle it
Measure the time-resolved domain radius $R(t)$ directly in the same colloidal square-ice simulations or experiments during the ordered-phase part of the quench. If $R(t)$ does not grow as $t^{1/2}$, or if the defect density is not proportional to $1/R(t)$, the coarsening attribution for $\rho_d\sim \tau_Q^{-1/2}$ is unsupported. Alternatively, run quenches much faster than 10 s: if the defect density then crosses to the $2/3$ exponent, the coarsening-only interpretation would need revision.
Extended reading notes
Core claim
In square colloidal ice, quenches from $B=0$ to $B=40$ mT leave a fraction $N_{ngs}/N$ of non-ground-state vertices, monopoles and other ice-rule violations, that decreases with quench duration. Fitting $N_{ngs}/N \sim \tau_Q^{-\beta}$ at fixed field gives $\beta \simeq 0.45$ for fields above about 12 mT, close to $1/2$, while the Kibble-Zurek mechanism for the 2D Ising universality class ($\nu=1$, $z=2$, $D=2$) predicts $\beta = D\nu/(1+z\nu)=2/3$. A rescaling collapse of the time traces uses exponent $\alpha=3/4$, also distinct from the KZ value $z\nu/(1+z\nu)=2/3$. The paper concludes that the defect population is governed by coarsening: after crossing the transition, ordered domains grow as $R(t)\sim t^{1/z}$ with $z=2$, the defect density tracks $1/R(t)$, and defect-defect Coulomb interactions are strong enough to keep annihilating defects during the ordered-phase part of the quench, which violates the KZ adiabatic-impulse assumption.
Load-bearing premise
The central assumption is that after the transition the ordered patches grow in size as the square root of time and the defect count falls as one over that patch size, so the measured exponent can be read as coarsening; the paper does not measure the patch growth directly.
Editorial extensions
If this is right
- Deep quenches in colloidal square ice produce a defect density that follows $\rho_d \sim \tau_Q^{-1/2}$ over the experimentally accessible range of quench times from 10 to 6000 s.
- Because the 2D Ising Kibble-Zurek prediction is $2/3$, square-ice systems with long-ranged monopole interactions should show the coarsening exponent rather than the KZ exponent.
- Quenches ending close to the critical field, between about 10 and 12 mT, should show a crossover regime with smaller exponents between 0.2 and 0.3, where coarsening and Kibble-Zurek compete.
- Hexagonal colloidal ice, which lacks an ordering transition, should show no power-law decay of defects with quench rate, since its ice-rule violations are not topologically protected and can disappear individually.
- Magnetic nanoisland or superconducting artificial ices with weaker defect-defect interactions, or much faster quenches, may access a genuine Kibble-Zurek regime.
Reading between the lines
- If the coarsening attribution holds, the $\tau_Q^{-1/2}$ defect scaling should be largely independent of the microscopic realization of square ice, provided the monopole interactions remain long-ranged; a test would be to compare colloid, nanomagnetic, and superconducting square ices at matched reduced quench rates.
- The unexplained time-rescaling exponent $\alpha=3/4$ may carry information beyond simple domain growth; measuring the domain-size distribution $P(R,t)$ directly during quenches could reveal whether the collapse reflects growing correlations, active annihilation, or the approach to the critical point.
- A natural extension is to inject controlled quenched disorder: if disorder slows monopole motion more than it slows domain growth, the coarsening exponent should cross over toward logarithmic decay and the Kibble-Zurek regime might become visible at slower quench rates.
- The square-versus-hexagonal contrast suggests that quench-protocol experiments in other frustrated geometries could use the presence or absence of a power-law defect density as a probe of whether the system has a true ordering critical point.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports molecular-dynamics simulations of quenches in colloidal square and hexagonal artificial spin ice, sweeping the magnetic field from B=0 to B=40 mT at ramp durations tau_Q between 10 and 6000 s. For square ice, the fraction of non-ground-state vertices Nngs/N at fixed final fields decays as a power law in tau_Q for fields above about 12 mT, with fitted exponent beta approximately 0.45, and the full Nngs/N-versus-time curves collapse when time is rescaled by tau_Q^0.75. The authors interpret these observations as evidence that deep quenches are governed by critical coarsening rather than by the Kibble-Zurek mechanism, because the 2D Ising KZ predictions for both exponents are 2/3. In the hexagonal ice, no power-law defect decay is found, which is attributed to the absence of a true ordering transition. The paper concludes that the quenched square ice provides a test bed for coarsening-dominated defect dynamics.
Significance. If the conclusion is correct, this is a valuable contribution to nonequilibrium dynamics in artificial spin ice and to the KZ-versus-coarsening debate: it is a particle-resolved MD study with 100 realizations per ramp over nearly three decades in tau_Q, directly mimicking an experimentally realizable colloidal system, and it includes the hexagonal ice as a control. The empirical scaling results and the defect-vertex characterization are strengths. The interpretation is, however, currently not secure because the KZ null hypothesis is computed with 2D Ising exponents that are assumed rather than measured for the colloidal model, the final defect density is measured after substantial post-transition dynamics, and the coarsening attribution is not verified by direct domain-growth measurements. These gaps are central to the paper's headline conclusion.
major comments (3)
- [Kibble-Zurek Mechanism (p. 3)] The rejection of the KZ mechanism rests on comparing the measured exponents alpha=0.75 and beta approximately 0.45 with the 2D Ising predictions z*nu/(1+z*nu)=2/3 and D*nu/(1+z*nu)=2/3. The manuscript assumes nu=1 and z=2 for the colloidal square ice, but on p. 4 it states that the colloidal ice can only be mapped exactly into a magnetic square ice at equilibrium and that its out-of-equilibrium kinetics may differ. Static universality does not fix the dynamic exponent, and Ref. [39] is an equilibrium result that cannot establish z=2 for this particle-based model. If the equilibrium critical point of this model had nu approximately 0.9 and z approximately 3.3, the KZ predictions would be z*nu/(1+z*nu)=0.75 and 2*nu/(1+z*nu)=0.45, matching the reported values almost exactly. The manuscript must therefore measure, or convincingly bound, nu and z for the colloidal model, or else the data cannot be said to rule out KZ.
- [Results, Kibble-Zurek paragraph (p. 4)] The positive claim that coarsening governs the defect density is not directly tested. The paper assumes ordered regions grow as R(t) approximately t^(1/z) with z=2 and that rho_d is approximately 1/R(t), but no measurement of R(t) is presented, and the collapse exponent alpha=0.75 is an empirical fit rather than a coarsening prediction. In addition, the measured Nngs/N is the final value at B=40 mT, after the system has spent considerable time on the ordered side of the transition; a KZ test based on freeze-out scaling would require either measuring the defect density at the freeze-out time or explicitly accounting for the subsequent coarsening. A direct measurement of domain growth and a demonstration that Nngs*R(t) stays constant, or a quantitative coarsening prediction for rho_d(tau_Q), is needed before the claim that coarsening governs can be asserted.
- [Abstract, Results (Fig. 4d), Conclusion] The central scaling statement is given as rho_d proportional to tau_Q^(-1/2) in the abstract and conclusion, while the fitted exponent reported in the text and in Fig. 4(d) is beta=0.45. This discrepancy matters because the KZ-versus-coarsening discrimination is exponent-based. Please report the fitted exponent with its uncertainty and justify any rounding to 1/2; if 0.45 is the best estimate, the paper should state rho_d proportional to tau_Q^(-0.45).
minor comments (4)
- [Conclusion (p. 5)] The phrase 'the university class of the square ice' should read 'universality class.'
- [Results, Fig. 4(d)] Figure 4(d) shows exponents for B<9 mT even though the text says the system does not order for B<9 mT; please clarify whether these small-field exponents are meaningful or are fitting artifacts.
- [Hexagonal system (p. 5)] The inset of Fig. 6 reports a collapse with alpha=0.88 for a system with no critical point; the text does not explain why a clean time-rescaling collapse is expected in a crossover regime, and a sentence on this point would help the reader.
- [Introduction (p. 2)] Minor grammar: 'there is a underlying second-order phase transition' should be 'there is an underlying second-order phase transition.'
Circularity Check
No significant circularity: the measured quench exponents are compared against external KZ and coarsening predictions, not derived from them.
full rationale
The central result, rho_d ∝ tau_Q^{-1/2}, is obtained by fitting the simulated defect density versus quench time in Fig. 4(c); it is a measurement, not an output of the coarsening or KZ formulas. The KZ null hypothesis is constructed from stated external inputs: the paper takes the square ice to be in the 2D Ising universality class with nu=1 and z=2, citing Fan and Wu [39], and computes z*nu/(1+z*nu)=2/3 and D*nu/(1+z*nu)=2/3. The measured alpha=3/4 and beta≈0.45 are then compared with those fixed numbers. Similarly, the coarsening expectation is quoted from Hohenberg and Halperin [40] (R(t)∝t^{1/z}, z=2 for Ising) and the qualitative relation rho_d∝1/R(t); neither is fitted from the same data and then renamed a prediction. The paper's self-citations (Refs. [14,35,43]) support the colloidal model, the earlier observation of monopole Coulomb dynamics, and the equilibrium-mapping caveat; none of them supplies a uniqueness theorem or a fitted parameter that forces the quench scaling, and the Methods section gives the full simulation parameters independently. The genuine caveat is a correctness risk, not circularity: the rejection of KZ depends on the assumed Ising values nu=1, z=2 for this particle-based model, and the paper itself notes that the colloidal square ice can only be mapped exactly at equilibrium. If the relevant equilibrium exponents were instead nu≈0.9 and z≈3.3, KZ would predict alpha=z*nu/(1+z*nu)=0.75 and beta=2*nu/(1+z*nu)=0.45, matching the data; the paper does not measure nu and z. This makes the KZ-versus-coarsening discrimination less secure, but the measured exponents are not defined in terms of the conclusion, so the derivation is not circular.
Assumptions & free parameters
free parameters (1)
- B_c (critical magnetic field) =
9 mT
assumptions (3)
- domain assumption The overdamped Langevin equation with harmonic double-well traps and r^{-4} repulsive colloid interactions faithfully represents the experimental colloidal artificial spin ice dynamics.
- domain assumption The equilibrium order-disorder transition of the colloidal square ice belongs to the 2D Ising universality class with nu = 1 and z = 2.
- ad hoc to paper In the coarsening regime, the defect density is proportional to 1/R(t) where the ordered region radius R(t) grows as t^{1/z} with z = 2.
Cite this review
Pith. "Pith review of Quenched Dynamics of Artificial Spin Ice: Coarsening versus Kibble-Zurek." pith.science (2026). https://pith.science/paper/5342TP44
@misc{pith2026190805328,
author = {Pith},
title = {Pith review of: Quenched Dynamics of Artificial Spin Ice: Coarsening versus Kibble-Zurek},
year = {2026},
howpublished = {\url{https://pith.science/paper/5342TP44}},
note = {Machine review of arXiv:1908.05328}
}
read the original abstract
Artificial spin ices are ideal frustrated model systems in which to explore or design emergent phenomena with unprecedented characterization of the constituent degrees of freedom. In square spin ice, violations of the ice rule are topological excitations essential to the kinetics of the system, providing an ideal testbed for studying the dynamics of such defects under varied quench rates. In this work we describe the first test of the Kibble-Zurek mechanism and critical coarsening in colloidal square and colloidal hexagonal ice under quenches from a weakly interacting liquid state into a strongly interacting regime. As expected, for infinitely slow quenches, the system is defect free, while for increasing quench rate, an increasing number of defects remain in the sample. For square ice, we find regimes in which the defect population decreases as a power law with decreasing quench rate. A detailed scaling analysis shows that for a wide range of parameters, including quench rates that are accessible by experiments, the behavior is described by critical coarsening rather than by the Kibble-Zurek mechanism, since the defect-defect interactions are long ranged. For quenches closer to the critical point, however, there can be a competition between the two mechanisms.
Figures
Figures from the paper (3 more)
Forward citations
Cited by 1 Pith paper
-
Transition between ground states in square anisotropic artificial colloidal ice
In square artificial colloidal ice, fast field rotation reaches the anisotropic 4-in/4-out ground state by diffusionless synchronized flips, while slow rotation traps the system in a partially ordered metastable state.
Reference graph
Works this paper leans on
- [39]
-
[1]
R. F. Wang, C. Nisoli, R. S. Freitas, J. Li, W. McConville, B. J. Cooley, M. S. Lund, N. Samarth, C. Leighton, V. H. Crespi, et al., Nature (London) 439, 303 (2006)
work page 2006
- [2]
- [3]
-
[4]
Y. Qi, T. Brintlinger, and J. Cumings, Phys. Rev. B 77, 094418 (2008)
work page 2008
- [5]
-
[6]
E. Mengotti, L. J. Heyderman, A. F. Rodr´ ıguez, F. Nolt- ing, R. V. H¨ ugli, and H.-B. Braun, Nature Phys. 7, 68 (2011)
work page 2011
- [7]
Show all 51 references
-
[8]
J. P. Morgan, A.Stein, S. Langridge, and C. H. Marrows, Nature Phys. 7, 75 (2011), URL http://dx.doi.org/ 10.1038/nphys1853
2011 doi
-
[9]
Wang, Z.-L
Y.-L. Wang, Z.-L. Xiao, A. Snezhko, J. Xu, L. Ocola, R. Divan, J. Pearson, G. Crabtree, and W.-K. Kwok, Science 352, 962 (2016)
2016
-
[10]
Lib´ al, C
A. Lib´ al, C. Reichhardt, and C. J. O. Reichhardt, Phys. Rev. Lett. 97, 228302 (2006). 7
2006
-
[11]
Ortiz-Ambriz and P
A. Ortiz-Ambriz and P. Tierno, Nature Commun. 7, 10575 (2016), URL http://dx.doi.org/10.1038/ ncomms10575
2016
-
[12]
Loehr, A
J. Loehr, A. Ortiz-Ambriz, and P. Tierno, Phys. Rev. Lett. 117, 168001 (2016)
2016
-
[13]
D. Y. Lee and P. Tierno, Phys. Rev. Mater. 2, 112601 (2018)
2018
-
[14]
Lib´ al, D
A. Lib´ al, D. Y. Lee, A. Ortiz-Ambriz, C. Reichhardt, C. J. O. Reichhardt, P. Tierno, and C. Nisoli, Nature Commun. 9, 4146 (2018)
2018
-
[15]
Lib´ al, C
A. Lib´ al, C. J. O. Reichhardt, and C. Reichhardt, Phys. Rev. Lett. 102, 237004 (2009)
2009
-
[16]
M. L. Latimer, G. R. Berdiyorov, Z. L. Xiao, F. M. Peeters, and W. K. Kwok, Phys. Rev. Lett. 111, 067001 (2013)
2013
-
[17]
Trastoy, M
J. Trastoy, M. Malnou, C. Ulysse, R. Bernard, N. Bergeal, G. Faini, J. Lesueur, J. Briatico, and J. E. Villegas, Nature Nanotechnol. 9, 710 (2014)
2014
-
[18]
Xue, J.-Y
C. Xue, J.-Y. Ge, A. He, V. S. Zharinov, V. V. Moshchalkov, Y. H. Zhou, A. V. Silhanek, and J. V. de Vondel, Phys. Rev. B 96, 024510 (2017)
2017
-
[19]
Y.-L. Wang, X. Ma, J. Xu, Z.-L. Xiao, A. Snezhko, R. Di- van, L. E. Ocola, J. E. Pearson, B. Janko, and W.-K. Kwok, Nature Nanotechnol. 13, 560 (2018)
2018
-
[20]
M¨ oller and R
G. M¨ oller and R. Moessner, Phys. Rev. Lett.96, 237202 (2006)
2006
-
[21]
Lib´ al, C
A. Lib´ al, C. Nisoli, C. Reichhardt, and C. Reichhardt, Phys. Rev. Lett. 120, 027204 (2018)
2018
-
[22]
Sendetskyi, V
O. Sendetskyi, V. Scagnoli, N. Leo, L. Anghinolfi, A. Al- berca, J. L¨ uning, U. Staub, P. M. Derlet, and L. J. Hey- derman, Phys. Rev. B 99, 214430 (2019)
2019
-
[23]
T. W. B. Kibble, J. Phys. A: Math. Gen. 9, 1387 (1976)
1976
-
[24]
W. H. Zurek, Nature (London) 317, 505 (1985)
1985
-
[25]
Monaco, J
R. Monaco, J. Mygind, R. J. Rivers, and V. P. Koshelets, Phys. Rev. B 80, 180501(R) (2009)
2009
-
[26]
M. J. Bowick, L. Chandar, E. A. Schiff, and A. M. Sri- vastava, Science 263, 943 (1994)
1994
-
[27]
C. N. Weiler, T. W. Neely, D. R. Scherer, A. S. Bradley, M. J. Davis, and B. P. Anderson, Nature (London) 455, 948 (2008)
2008
-
[28]
Lamporesi, S
G. Lamporesi, S. Donadello, S. Serafini, F. Dalfovo, and G. Ferrari, Nature Physics 9, 656 (2013), URL http: //dx.doi.org/10.1038/nphys2734
2013 doi
-
[29]
Navon, A
N. Navon, A. L. Gaunt, R. P. Smith, and Z. Hadzibabic, Science 347, 167 (2015)
2015
-
[30]
S. Ulm, J. Roßnagel, G. Jacob, C. Deg¨ unther, S. T. Dawkins, U. G. Poschinger, R. Nigmatullin, A. Retzker, M. B. Plenio, F. Schmidt-Kaler, et al., Nature Commun. 4, 2290 (2013)
2013
-
[31]
K. Pyka, J. Keller, H. L. Partner, R. Nigmatullin, T. Burgermeister, D. M. Meier, K. Kuhlmann, A. Ret- zker, M. B. Plenio, W. H. Zurek, et al., Nature Commun. 4, 2291 (2013)
2013
-
[32]
S. M. Griffin, M. Lilienblum, K. T. Delaney, Y. Kumagai, M. Fiebig, and N. A. Spaldin, Phys. Rev. X 2, 041022 (2012)
2012
-
[33]
del Campo and W
A. del Campo and W. H. Zurek, Int. J. Mod. Phys. A 29, 1430018 (2014)
2014
-
[34]
Biroli, L
G. Biroli, L. F. Cugliandolo, and A. Sicilia, Phys. Rev. E 81, 050101(R) (2010)
2010
-
[35]
Lib´ al, C
A. Lib´ al, C. Nisoli, C. Reichhardt, and C. Reichhardt, Sci. Rep. 7, 651 (2017)
2017
-
[36]
C. W. Liu, A. Polkovinkov, and A. W. Sandvik, Phys.Rev. B 89, 054307 (2014)
2014
-
[37]
Budrikis, K
Z. Budrikis, K. L. Livesey, J. P. Morgan, J. Akerman, A. Stein, S. Langridge, C. H. Marrows, and R. L. Stamps, New J. Phys. 14, 035014 (2012)
2012
-
[38]
Drisko, T
J. Drisko, T. Marsh, and J. Cumings, Nature Commun. 8, 14009 (2017)
2017
-
[40]
P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49, 435 (1977)
1977
-
[41]
Wu, Physical Review Letters 22, 1174 (1969)
F. Wu, Physical Review Letters 22, 1174 (1969)
1969
-
[42]
Nisoli, New Journal of Physics 16, 113049 (2014)
C. Nisoli, New Journal of Physics 16, 113049 (2014)
2014
-
[43]
Nisoli, Phys
C. Nisoli, Phys. Rev. Lett. 120, 167205 (2018)
2018
-
[44]
Levis, L
D. Levis, L. F. Cugliandolo, L. Foini, and M. Tarzia, Phys. Rev. Lett. 110, 207206 (2013)
2013
-
[45]
Kapaklis, U
V. Kapaklis, U. B. Arnalds, A. Farhan, R. V. Chopdekar, A. Balan, A. Scholl, L. J. Heyderman, and B. Hj¨ orvarsson, Nature Nanotechnol.9, 514 (2014)
2014
-
[46]
Anghinolfi, H
L. Anghinolfi, H. Luetkens, J. Perron, M. G. Flokstra, O. Sendetskyi, A. Suter, T. Prokscha, P. M. Derlet, S. L. Lee, and L. J. Heyderman, Nature Commun. 6, 8278 (2015)
2015
-
[47]
Farhan, C
A. Farhan, C. F. Petersen, S. Dhuey, L. Anghinolfi, Q. H. Qin, M. Saccone, S. Velten, C. Wuth, S. Gliga, P. Mel- lado, et al., Nature Commun. 8, 995 (2017)
2017
-
[48]
Y. Lao, F. Caravelli, M. Sheikh, J. Sklenar, D. Gardeaza- bal, J. D. Watts, A. M. Albrecht, A. Scholl, K. Dahmen, C. Nisoli, et al., Nature Phys. 14, 723 (2018)
2018
-
[49]
Sklenar, Y
J. Sklenar, Y. Lao, A. Albrecht, J. D. Watts, C. Nisoli, G.-W. Chern, and P. Schiffer, Nature Phys. 15, 191 (2019)
2019
-
[50]
Perrin, B
Y. Perrin, B. Canals, and N. Rougemaille, Nature (Lon- don) 540, 410 (2016)
2016
-
[51]
Farhan, M
A. Farhan, M. Saccone, C. F. Petersen, S. Dhuey, R. V. Chopdekar, Y.-L. Huang, N. Kent, Z. Chen, M. J. Alava, T. Lippert, et al., Sci. Adv. 5, 6380 (2019)
2019
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.