REVIEW 3 major objections 5 minor 64 references
Discrete Time Crystal in quantum Sherrington-Kirkpatrick model
T0 review · 3 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read This paper argues that a periodically driven quantum Sherrington-Kirkpatrick spin glass can host a stable discrete time crystal phase despite its all-to-all random Ising interactions.
desk verdict Careful ED study of DTC in the SK model, but the thermodynamic-phase claim outruns L<=14 and 600T; worthwhile refereeing. 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 load-bearing object is the two-step Floquet unitary $U_F = e^{-iH_2} e^{-iH_1}$, with $H_1 = \frac{\pi}{2}\sum_i \sigma_i^x$ acting as a perfect spin flip and $H_2$ the quantum SK Hamiltonian, namely $\sum_{i<j} J_{ij}\sigma_i^z\sigma_j^z$ plus either a transverse field $\sum_i h_i^x \sigma_i^x$ or nearest-neighbour XY coupling $\sum_i g_i(\sigma_i^x\sigma_{i+1}^x+\sigma_i^y\sigma_{i+1}^y)$. The diagnostic is the stroboscopic spin autocorrelation $C(t) = \langle \psi_0|\sigma_i^z(t)\sigma_i^z(0)|\psi_0\rangle$, whose persistent period-2 oscillations with a sharp Fourier peak at half the drive frequency define the DTC. The mechanism invoked is eigenstate localization of $H_2$: the normalized Shannon entropy $S(E_n)/\ln(2^L)$, with $S(E_n)=-\sum_i |\psi_n(i)|^2 \ln|\psi_n(i)|^2$, and the level-spacing ratio $r_n$ are used to show that non-ergodic eigenstates keep the spin-flip information from heating, and the DTC order parameter $O = \frac{1}{N}\sum_n [(-1)^n C(nT)-C(nT)]$ tracks the low-entropy regime.
What would settle it
At a claimed DTC point, for example J=2.0 with random $h^x$ of strength 0.2, measure the decay time of the period-doubling autocorrelation for L=10, 12, 14, and 16: if the decay time does not grow with system size, or if the eigenstates of the full Floquet unitary $U_F$ show level repulsion instead of two separated localized sectors whose eigenphases differ by $\pi$, then the order is prethermal, not a stable DTC.
Extended reading notes
Core claim
The central claim is that the periodically driven quantum SK model exhibits a genuine DTC phase for a finite range of transverse-field or XY-coupling strengths, even though every spin pair is randomly coupled. The stability condition is traced to the eigenstates of H2, the undriven quantum SK Hamiltonian in the second half of the drive: when those eigenstates are many-body localized, with low normalized Shannon entropy and Poisson level-spacing ratio, an initial state overlaps only a small fraction of them, so the pure spin-flip pulse in the first half of the drive produces period-2 autocorrelation oscillations that persist without visible decay; when the eigenstates are ergodic, the oscillations decay. The DTC occupies only part of the non-ergodic region, and the DTC order parameter closely follows the eigenstate Shannon entropy of H2. In a comparison with power-law interacting models, uniform long-range couplings fail to stabilize the DTC under most protocols, while adding random XY couplings can induce it even there.
Load-bearing premise
The load-bearing premise is that the undriven spin-glass Hamiltonian H2 remains non-thermalizing for system sizes larger than the L=14 simulated here, so that it never heats up under the periodic drive and the period-doubling seen in finite-time simulations is a true phase rather than an extremely long transient.
Editorial extensions
If this is right
- If the central claim is correct, the class of DTC-supporting systems extends from short-range disordered chains to all-to-all random spin glasses, so short-range interactions are not a necessary ingredient.
- The DTC order parameter closely following the Shannon entropy of H2 means the DTC parameter window can be predicted from eigenstate localization of the undriven Hamiltonian alone.
- In the third driving protocol the DTC phase width in the spin-flip error $\epsilon$ is roughly independent of the SK coupling strength J, unlike protocols I and II, so the robustness window is set by the drive itself.
- Random XY couplings stabilize period doubling even when the long-range Ising part is uniform, so a short-range random spin-flip term can rescue DTC order in otherwise uniform long-range models.
Reading between the lines
- An extension the authors do not pursue: the phase boundary $h_c(J)$ could be estimated from the mobility edge of H2; if the Shannon-entropy contour of 0.5 at fixed J coincides with the DTC boundary at larger L, the correspondence would be a predictive tool.
- If the finite-size non-ergodicity does not survive the thermodynamic limit, the observed DTC would instead be a very long prethermal plateau typical of all-to-all models; a trapped-ion or Rydberg simulator with programmable all-to-all couplings could distinguish the two by measuring autocorrelation decay over tens of thousands of cycles.
- The random-XY result suggests a broader design rule: adding short-range random hopping to a non-thermalizing long-range spin Hamiltonian may generically stabilize period doubling, connecting DTC physics to random-hopping localization.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies the periodically driven quantum Sherrington-Kirkpatrick (SK) model under three two-step driving protocols: a spin-flip pulse followed by evolution under H2, which contains random all-to-all Ising zz couplings plus either a transverse field (protocol I), nearest-neighbor XY coupling (protocol II), or only Ising zz interactions (protocol III). Starting from random product states in the z basis, the authors compute stroboscopic spin autocorrelations and a DTC order parameter, observing period-2T oscillations with a peak at half the drive frequency. They identify DTC regions in the J-hx and J-g planes and correlate them with the level-spacing ratio and Shannon entropy of the eigenstates of H2. The central claim is that the quantum SK model hosts a robust DTC phase despite all-to-all random long-range interactions, that randomness in the transverse field or XY coupling broadens the DTC region, and that uniform power-law spin-preserving interactions do not stabilize a DTC except in the presence of random XY couplings.
Significance. If correct, the result would substantially extend the DTC paradigm from short-range disordered systems to all-to-all random spin glasses and would identify random spin-flip terms as an important stabilizing ingredient. The paper's strengths are the use of three driving protocols, a direct comparison with uniform power-law models, and the use of independent static diagnostics (r-ratio and Shannon entropy) that are not fitted to the dynamical order parameter. However, the evidence for a thermodynamic DTC phase is currently indirect: the system sizes are limited to L=14 (L=20 for protocol III), the evolution times are finite (150-600 periods, with a mention of 10^4 periods but no data shown), no Floquet spectral analysis is presented, and no error bars are given for the phase boundaries. The paper is therefore a useful numerical study whose central claim needs additional support before it can be accepted at the level of a robust phase.
major comments (3)
- [Section III, protocols I-III (Figs. 2, 3, 6, 8)] The inference from finite-time period doubling to a 'robust DTC phase' in the thermodynamic limit is not yet supported. The largest systems are L=14 (L=20 only for protocol III), the displayed evolution times are 150-600 periods, and the statement in the text that oscillations persist to a time scale of 10^4 is not backed by a figure or a quantitative analysis. Since prethermal time-crystalline plateaus can persist for times exponential in a control parameter, the observed C(t) and order parameter O do not by themselves distinguish a genuine DTC from a transient. I request a finite-size scaling of the period-doubling lifetime at fixed parameters (for example, the time at which O decays below a threshold as a function of L) or an analysis of the Floquet spectrum, showing symmetry-related Floquet eigenstate pairs and long-range order in the eigenstates. Without one of these, the thermodynamic 'phase' claim should be softened to a finite-time, finite-size observation.
- [Section III, phase diagrams (Figs. 3, 5, 6, 8)] The claimed phase boundaries hc, gc, and the epsilon-window in protocol III are extracted by eye from color plots of the order parameter, with no quantitative criterion and no statistical uncertainty. The order parameter is computed for a single initial state in the phase diagrams, averaged over 25-50 disorder realizations, and the threshold used to call a region 'DTC' (the text mentions O>1/2) is not defined precisely. Since a central conclusion is that the DTC region is broader for random transverse fields and random XY couplings, the authors should define the criterion (for example, time-averaged O>0.5 over a fixed number of periods) and report the disorder and initial-state statistics, including the variance of O across realizations.
- [Section III and Appendix A (Figs. 4, 7, 10)] The non-ergodicity of H2 is established only at L=10-14 via the level-spacing ratio and the normalized Shannon entropy, and the paper then uses this to argue that localized eigenstates of H2 stabilize the DTC. For all-to-all couplings, the number of perturbative resonances grows with system size, and the manuscript itself cites works (refs. [45,51]) showing a delocalization tendency of random long-range spin-preserving interactions. Poisson level statistics at L<=14 and S/ln(2^L)<0.5 are therefore not sufficient to conclude MBL-like non-ergodicity in the thermodynamic limit. The sentence in Section III, 'If H2 represents a non-ergodic system whose eigenstates are MBL', is an assumption rather than a demonstrated fact. Please add finite-size scaling of the r-ratio and Shannon entropy at fixed parameters, or explicitly restrict the claims to finite-size non-ergodicity. This is load-bearing because the stability of the DTC is attributed to this non-ergodicity.
minor comments (5)
- [Section III, Fig. 3] The statement that the DTC order parameter 'follows Shannon entropy very closely' is qualitative; a scatter plot of O versus S/ln(2^L) across the parameter plane, or a correlation measure, would make this quantitative.
- [Section II, Eq. (4)] The choice N=300 in the time-averaged order parameter is not discussed; please state the time window in periods and justify that the average has converged.
- [Section IV, Fig. 9] The negative results for power-law interacting systems are presented for alpha=1,2 and L=14; a brief comment on finite-size effects for these conclusions would be helpful, since the absence of period doubling is itself a finite-size statement.
- [Appendix B, Fig. 11] The uniform-XY data in Appendix B are averaged over only 10 disorder realizations; it would be useful to report error bars or at least state that the qualitative conclusions are stable across the realizations.
- [General] The text uses 'MBL' for an all-to-all model without spatial locality; conventional MBL refers to local Hamiltonians, so 'non-ergodic' or 'localized' would be safer terminology unless the authors define precisely what they mean by MBL in this context.
Circularity Check
No significant circularity: the DTC order parameter and the eigenstate diagnostics are independently computed, and no fitted quantity is renamed as a prediction.
full rationale
None of the paper's load-bearing steps reduces to its own inputs. The DTC order parameter O (Eq. 4) is computed directly from the stroboscopic autocorrelation C(nT), which is obtained by exact-diagonalization time evolution under the Floquet unitary. The Shannon entropy and level-spacing ratio are independent static diagnostics of the eigenstates of H2, computed from a separate spectral analysis. The observed overlap between the O>0.5 region and the S/ln(2^L)<0.5 region is presented as an empirical correlation, not as a derivation of one from the other; neither quantity is fitted to reproduce the other. The paper's self-citations (refs. 48, 50, 51, 57) supply background on SK-model non-ergodicity and long-range MBL, but the central DTC evidence is this paper's own numerical simulation, so those citations are not load-bearing. Concerns about finite-size scaling, prethermal transients, and the absence of a Floquet spectral analysis are legitimate scientific risks, but they are not circularity.
Assumptions & free parameters
assumptions (5)
- standard math Exact diagonalization gives exact finite-size dynamics for the driven SK model.
- domain assumption H2 has a non-ergodic eigenstate phase in the thermodynamic limit.
- domain assumption Non-ergodic eigenstates of H2 prevent heating under the periodic drive and stabilize the DTC.
- domain assumption Watanabe-Oshikawa no-go theorem does not forbid discrete time crystals in long-range systems.
- domain assumption The order parameter O with N=300 and the O>0.5 threshold define the DTC phase.
Cite this review
Pith. "Pith review of Discrete Time Crystal in quantum Sherrington-Kirkpatrick model." pith.science (2026). https://pith.science/paper/WTNODECN
@misc{pith2026250419378,
author = {Pith},
title = {Pith review of: Discrete Time Crystal in quantum Sherrington-Kirkpatrick model},
year = {2026},
howpublished = {\url{https://pith.science/paper/WTNODECN}},
note = {Machine review of arXiv:2504.19378}
}
abstract
Discrete time crystals (DTC) have emerged as a significant phase of matter for out-of-equilibrium many-body systems. We study how long-range interactions and disorder contribute to the stability of the DTC phase. Generally, a stable DTC phase is believed to be realized in disordered systems with short-range interactions. In this work, we study periodically driven quantum Sherrington-Kirkpatrick (SK) model of Ising spin-glass in which all spins are randomly coupled. We investigate the possibilities of DTC phase in the SK model within three different driving protocols and found that the quantum SK model exhibits a robust DTC phase despite the long-range nature of interactions. The DTC phase in quantum SK model persists for a larger range of parameters if the $XY$ coupling or a transverse field is also random. This suggests that disorder in the $XY$ coupling or transverse field is also crucial for stabilising a broad DTC phase, despite the SK model's random couplings. Our analysis shows that the stability of the DTC phase is determined by the non-ergodic nature of the eigenstates of the quantum SK model, with the DTC order-parameter closely following the Shannon entropy of eigenstates. We compare the periodically driven SK model to alternative models of long-range interactions with uniform coefficients and found that the DTC phase is absent in these models for most of the driving protocols.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
Wilczek, Phys
F. Wilczek, Phys. Rev. Lett. 109, 160401 (2012)
2012
- [2]
-
[3]
Nozi` eres, Europhysics Letters103, 57008 (2013)
P. Nozi` eres, Europhysics Letters103, 57008 (2013)
2013
-
[4]
Watanabe and M
H. Watanabe and M. Oshikawa, Phys. Rev. Lett. 114, 251603 (2015)
2015
-
[5]
D’Alessio and M
L. D’Alessio and M. Rigol, Phys. Rev. X 4, 041048 (2014)
2014
-
[6]
D. M. Basko, I. L. Aleiner, and B. L. Altshuler, Annals of Physics 321, 1126 (2006)
2006
-
[7]
I. V. Gornyi, A. D. Mirlin, and D. G. Polyakov, Phys. Rev. Lett. 95, 206603 (2005)
2005
-
[8]
E. Altman and R. Vosk, Annual Review of Condensed Matter Physics 6, 383 (2015)
work page 2015
Show all 64 references
-
[9]
Alet and N
F. Alet and N. Laflorencie, Comptes Rendus Physique 19, 498 (2018)
2018
-
[10]
Nandkishore and D
R. Nandkishore and D. A. Huse, Annual Review of Con- densed Matter Physics 6, 15 (2015)
2015
-
[11]
D. A. Abanin, E. Altman, I. Bloch, and M. Serbyn, Reviews of Modern Physics 91, 021001 (2019)
2019
-
[12]
Sierant, M
P. Sierant, M. Lewenstein, A. Scardicchio, L. Vidmar, and J. Zakrzewski, Reports on Progress in Physics 88, 026502 (2025)
2025
-
[13]
D. A. Abanin, W. De Roeck, and F. Huveneers, Annals of Physics 372, 1 (2016)
2016
-
[14]
Ponte, A
P. Ponte, A. Chandran, Z. Papi´ c, and D. A. Abanin, Annals of Physics 353, 196 (2015)
2015
-
[15]
Ponte, Z
P. Ponte, Z. Papi´ c, F. m. c. Huveneers, and D. A. Abanin, Phys. Rev. Lett. 114, 140401 (2015)
2015
-
[16]
Khemani, A
V. Khemani, A. Lazarides, R. Moessner, and S. L. Sondhi, Phys. Rev. Lett. 116, 250401 (2016)
2016
-
[17]
D. V. Else, B. Bauer, and C. Nayak, Physical Review Letters 117, 090402 (2016)
2016
-
[18]
N. Y. Yao, A. C. Potter, I.-D. Potirniche, and A. Vish- wanath, Phys. Rev. Lett. 118, 030401 (2017)
2017
-
[19]
S. Choi, J. Choi, R. Landig, G. Kucsko, H. Zhou, J. Isoya, F. Jelezko, S. Onoda, H. Sumiya, V. Khemani, C. von Keyserlingk, N. Y. Yao, E. Demler, and M. D. Lukin, Nature 543, 221–225 (2017)
2017
-
[20]
A com- ment on
V. Khemani, R. Moessner, and S. L. Sondhi, “A com- ment on ”discrete time crystals: rigidity, criticality, and realizations”,” (2021), arXiv:2109.00551 [cond-mat.dis- nn]
2021 arXiv
-
[21]
Ippoliti, K
M. Ippoliti, K. Kechedzhi, R. Moessner, S. Sondhi, and V. Khemani, PRX Quantum 2, 030346 (2021)
2021
-
[22]
B. Li, J. S. Van Dyke, A. Warren, S. E. Economou, and E. Barnes, Phys. Rev. B 101, 115303 (2020)
2020
-
[23]
Kyprianidis, F
A. Kyprianidis, F. Machado, W. Morong, P. Becker, K. S. Collins, D. V. Else, L. Feng, P. W. Hess, C. Nayak, G. Pagano, N. Y. Yao, and C. Monroe, Science 372, 1192 (2021)
2021
-
[24]
Beatrez, C
W. Beatrez, C. Fleckenstein, A. Pillai, J. Ye, B. Song, W. S. Bakr, V. Khemani, and B. L. Lev, Nature Physics 19, 407 (2023)
2023
-
[25]
Stasiuk and P
A. Stasiuk and P. Cappellaro, Phys. Rev. X 13, 041016 (2023)
2023
-
[26]
Prethermal time-crystalline corner modes,
S. Jiang, D. Yuan, W. Jiang, D.-L. Deng, and F. Machado, “Prethermal time-crystalline corner modes,” (2024), arXiv:2406.01686 [quant-ph]
2024 arXiv
-
[27]
Unpolarized prethermal discrete time crystal,
T. Yokota and T. N. Ikeda, “Unpolarized prethermal discrete time crystal,” (2025), arXiv:2501.09461 [cond- mat.quant-gas]
2025 arXiv
-
[28]
Zhang, P
J. Zhang, P. W. Hess, A. Kyprianidis, P. Becker, A. Lee, J. Smith, G. Pagano, I. D. Potirniche, A. C. Potter, A. Vishwanath, N. Y. Yao, and C. Monroe, Nature 543, 217–220 (2017)
2017
-
[29]
Rovny, R
J. Rovny, R. L. Blum, and S. E. Barrett, Phys. Rev. Lett. 120, 180603 (2018)
2018
-
[30]
Rovny, R
J. Rovny, R. L. Blum, and S. E. Barrett, Phys. Rev. B 97, 184301 (2018)
2018
-
[31]
S. Pal, N. Nishad, T. S. Mahesh, and G. J. Sreejith, Phys. Rev. Lett. 120, 180602 (2018)
2018
-
[32]
B. Liu, L. Zhang, Q. Wang, et al. , Nature Communica- tions 15, 9730 (2024)
2024
-
[33]
Frey and S
P. Frey and S. Rachel, Science Advances 8, eabm7652 (2022)
2022
-
[34]
Randall, C
J. Randall, C. E. Bradley, F. V. van der Gronden, A. Galicia, M. H. Abobeih, M. Markham, D. J. Twitchen, F. Machado, N. Y. Yao, and T. H. Taminiau, Science 374, 1474 (2021)
2021
-
[35]
Khemani, R
V. Khemani, R. Moessner, and S. Sondhi, arXiv preprint arXiv:1910.10745 (2019)
2019 arXiv
-
[36]
M. P. Zaletel, M. Lukin, C. Monroe, C. Nayak, F. Wilczek, and N. Y. Yao, Rev. Mod. Phys. 95, 031001 (2023)
2023
-
[37]
M. P. Estarellas, T. Osada, V. M. Bastidas, B. Renoust, K. Sanaka, W. J. Munro, and K. Nemoto, Science Ad- vances 6, eaay8892 (2020)
2020
-
[38]
Discrete time crystal sensing,
L. J. I. Moon, P. M. Schindler, R. J. Smith, E. Druga, Z.-R. Zhang, M. Bukov, and A. Ajoy, “Discrete time crystal sensing,” (2024), arXiv:2410.05625 [quant-ph]
2024 arXiv
-
[39]
Yousefjani, K
R. Yousefjani, K. Sacha, and A. Bayat, Phys. Rev. B 111, 125159 (2025)
2025
-
[40]
Montenegro, M
V. Montenegro, M. G. Genoni, A. Bayat, et al., Commu- nications Physics 6, 304 (2023)
2023
-
[41]
Machado, D
F. Machado, D. V. Else, G. D. Kahanamoku-Meyer, C. Nayak, and N. Y. Yao, Phys. Rev. X 10, 011043 (2020)
2020
-
[42]
C. Lyu, S. Choudhury, C. Lv, Y. Yan, and Q. Zhou, Phys. Rev. Res. 2, 033070 (2020)
2020
-
[43]
Pizzi, J
A. Pizzi, J. Knolle, and A. Nunnenkamp, Nature Com- munications 12, 2341 (2021)
2021
-
[44]
A. L. Burin, Phys. Rev. B 91, 094202 (2015)
2015
-
[45]
N. Y. Yao, C. R. Laumann, S. Gopalakrishnan, M. Knap, M. M¨ uller, E. A. Demler, and M. D. Lukin, Phys. Rev. Lett. 113, 243002 (2014)
2014
-
[46]
D. B. Gutman, I. V. Protopopov, A. L. Burin, I. V. Gornyi, R. A. Santos, and A. D. Mirlin, Phys. Rev. B 93, 245427 (2016)
2016
-
[47]
K. S. Tikhonov and A. D. Mirlin, Phys. Rev. B 97, 214205 (2018)
2018
-
[48]
Nag and A
S. Nag and A. Garg, Phys. Rev. B 99, 224203 (2019)
2019
-
[49]
Roy and D
S. Roy and D. E. Logan, SciPost Physics 7, 042 (2019)
2019
-
[50]
A. Jana, V. R. Chandra, and A. Garg, Phys. Rev. B 109, 214209 (2024). 11
2024
-
[51]
Prasad and A
Y. Prasad and A. Garg, Phys. Rev. B103, 064203 (2021)
2021
-
[52]
Sherrington and S
D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 1792 (1975)
1975
-
[53]
Y. Y. Goldschmidt and P.-Y. Lai, Phys. Rev. Lett. 64, 2467 (1990)
1990
-
[54]
Lancaster and F
D. Lancaster and F. Ritort, Journal of Physics A: Math- ematical and General 30, L41 (1997)
1997
-
[55]
Takahashi, Phys
K. Takahashi, Phys. Rev. B 76, 184422 (2007)
2007
-
[56]
A. P. Young, Phys. Rev. E 96, 032112 (2017)
2017
-
[57]
Mukherjee, S
S. Mukherjee, S. Nag, and A. Garg, Phys. Rev. B 97, 144202 (2018)
2018
-
[58]
C. L. Baldwin, C. R. Laumann, A. Pal, and A. Scardic- chio, Phys. Rev. Lett. 118, 127201 (2017)
2017
-
[59]
E. Lieb, T. Schultz, and D. Mattis, Annals of Physics 16, 407 (1961)
1961
-
[60]
Theodorou and M
G. Theodorou and M. H. Cohen, Phys. Rev. B 13, 4597 (1976)
1976
-
[61]
T. P. Eggarter and R. Riedinger, Phys. Rev. B 18, 569 (1978)
1978
-
[62]
Dhar, Physica A: Statistical Mechanics and its Appli- cations 102, 370 (1980)
D. Dhar, Physica A: Statistical Mechanics and its Appli- cations 102, 370 (1980)
1980
-
[63]
Fleishman and D
L. Fleishman and D. C. Licciardello, Journal of Physics C: Solid State Physics 10, L125 (1977)
1977
-
[104]
This shows that the long-range nature of spin-preserving in- teractions enhance the non-ergodicity which is consistent with earlier works on the quantum SK model [57, 58]
and survives for a broad range of parameters. This shows that the long-range nature of spin-preserving in- teractions enhance the non-ergodicity which is consistent with earlier works on the quantum SK model [57, 58]. System with Uniform T ransverse Field: So far we have prese...
2000
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.