Pith. sign in

REVIEW 4 major objections 4 minor 61 references

The staggered-field XXZ spin chain, deep in its antiferromagnetic phase, stops being chaotic: spinon pairs bind into stable mesons whose low-lying energies follow a parameter-free Airy ladder that matches exact diagonalization closely.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review

2026-08-03 21:33 UTC pith:Y55256RR

load-bearing objection Solid ED study of confinement-induced ergodicity breaking in the staggered-field XXZ chain; the meson-ladder test is a real, externally grounded result, but the 'parameter-free' label oversells it because a global alignment removes the absolute offset. the 4 major comments →

arxiv 2511.14847 v2 pith:Y55256RR submitted 2025-11-18 cond-mat.str-el

Breakdown of Quantum Chaos in the Staggered-Field XXZ Chain: Confinement and Meson Formation

Julia Wildeboer , Marton Lajer , Robert M. Konik This is my paper
classification cond-mat.str-el
keywords confinementmeson formationXXZ spin chainstaggered fieldquantum chaosHilbert space fragmentationmany-body scarsAiry ladder
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper argues that confinement, not disorder or fine-tuning, destroys quantum chaos in the spin-1/2 XXZ chain with a staggered field. Deep in the antiferromagnetic Ising regime, spinons bind into domain-wall mesons, and the domain-wall number W emerges as a quasi-conserved quantum number. This shows up as a crossover from GOE to Poisson-like level statistics, as eigenstates banding by W in spin correlations and entanglement, and as a low-lying one-meson spectrum that matches an analytic Airy ladder with no fitted parameters. The low-lying W=2 meson energies at zero momentum agree closely with exact diagonalization, and the number of stable mesons below the two-meson threshold is predicted to within one or two levels. If right, the result makes confinement a generic, materials-relevant route to non-ergodicity, bridging Hilbert-space fragmentation and quantum many-body scars.

Core claim

The central claim is that the staggered-field XXZ chain realizes confinement-induced non-ergodicity: in the gapped antiferromagnetic phase the staggered field produces a linear potential between spinons (domain walls), binding them into mesons. The authors show that the low-lying W=2 (one-meson) spectrum at total momentum P=0 obeys, near the two-spinon continuum threshold, the parameter-free Airy ladder En^(+)(0) = Econt(0) + A(0) z_n + f a_+(0) + O(f^(4/3)), where Econt is the deconfined two-spinon edge, A(0) is set by the curvature of the spinon dispersion and scales as f^(2/3), z_n are the Airy zeros, f is the string tension of the staggered field, and a_+(0) = -gamma_E/pi is a universal

What carries the argument

The load-bearing object is the Airy quantization of the relative motion of two spinons in a linear confining potential. The staggered field creates a string tension f = 2h sigma-bar(eta), and near the two-spinon threshold the relative-coordinate bound states are controlled by Airy zeros z_n, yielding the meson ladder En = Econt + A(P) z_n + f a_+(P) + O(f^(4/3)) with A(P) = f^(2/3) [epsilon''(0|P)/2]^(1/3). The single-spinon dispersion omega(q) = I sqrt(1 - k^2 cos^2 q) is parameterized elliptically, with I fixed by the elliptic integral K(k) and sinh(eta), and Delta = -cosh(eta). The channel shift a_+(0) = -gamma_E/pi is universal, independent of J, eta, and k. The companion diagnostic is t

Load-bearing premise

The quantitative match assumes that N=22 periodic-boundary eigenstates selected as W=2 one-meson states are accurately described by the thermodynamic-limit two-spinon bound states of the analytic ladder, so that 1/N momentum quantization and O(f^(4/3)) corrections are negligible for the low-lying levels compared.

What would settle it

Plot (En - Econt(0) - f a_+(0))/A(0) against the Airy magnitudes z_n for N=22, Delta=-4.5, h=0.1 in the (S_z, P, I, C_flip) = (0,0,+1,+1) sector: if the first three W=2 levels do not collapse onto the line y=z within the reported fit residuals, or if the lowest W=4 level does not sit within one level spacing of the analytic two-meson threshold, the quantitative claim fails. The paper's own counting gives n*_ED = 8 versus n*_th = 10, so checking whether the n=9 and n=10 levels hybridize with the W=4 sector is a direct test.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Level statistics in the fixed symmetry sector cross from GOE to Poisson-like behavior as Delta is made more negative; a finite-size extrapolation locates the thermodynamic onset of non-ergodicity at Delta_tra = -4.19 (h=0.3).
  • The domain-wall number W acts as an emergent quasi-conserved quantum number deep in the AF regime, producing flat bands in the nearest-neighbor correlator C_zz and in the staggered magnetization.
  • The half-chain entanglement entropy reorganizes from a random-state dome into suppressed, band-resolved plateaus, with the W=2 one-meson band at markedly sub-random-state values.
  • The low-lying W=2 spectrum follows the parameter-free Airy ladder with an offset-removed collapse onto y=z and spacings proportional to A(0)(z_{n+1} - z_n); the number of stable one-meson levels below the two-meson threshold is n* ~ 8-10, and n* grows as f^(-2/3) as the staggered field weakens.
  • The correlation banding, sub-random-state entanglement plateaus, and near-threshold Airy spectra provide a template for meson spectroscopy in quasi-one-dimensional magnets with Ising anisotropy and effective staggered fields.

Where Pith is reading between the lines

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

  • If the Airy ladder is as reliable as the N=22 data suggest, the same offset-removed collapse should survive at larger N with residuals scaling as 1/N and f^(4/3); a direct finite-size scaling of the fitted slope A_fit against A(0) would sharpen the comparison.
  • The n*_ED = 8 versus n*_th = 10 discrepancy is a natural testbed: it likely reflects the two-meson threshold's own finite-size correction, so a systematic study of the lowest W=4 level versus N could reconcile the counting and pin down where the n=9 and n=10 resonances fully hybridize with the W=4 sector.
  • Because the mechanism is energetic rather than algebraic, similar meson banding should appear in other one-dimensional Ising-like models with a staggered field, and the C_zz banding could serve as a model-agnostic experimental signature.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper studies the spin-1/2 XXZ chain in a staggered magnetic field, focusing on the gapped antiferromagnetic regime (Δ<−1). It argues that spinon confinement produces meson-like bound states that break ergodicity. Using exact diagonalization in symmetry-resolved sectors, the authors report (i) a crossover in the mean adjacent-gap ratio from GOE to 'Poisson-like' statistics with increasing |Δ|, (ii) banding of eigenstates in nearest-neighbor correlations and bipartite entanglement, labeled by an emergent domain-wall number W, and (iii) quantitative agreement between the low-lying W=2 (one-meson) spectrum at P=0 and the analytic Airy ladder of Rutkevich, including continuum-relative bindings, an offset-removed Airy collapse, and two-meson thresholds. The paper connects these results to weak Hilbert-space fragmentation and quantum-scar-like nonergodicity.

Significance. The central asset is the quantitative meson spectroscopy in Sec. VI: the benchmark is an external, non-self-cited analytic theory (Rutkevich, Phys. Rev. B 106, 134405 (2022)), the ladder is used with no fitted parameters in its construction, and the comparison honestly reports residuals, the n*=8 vs 10 count discrepancy, and the limited low-n agreement window. If the level-statistics and banding evidence is robust, the paper would provide a valuable confinement-based route to weak fragmentation and scar phenomenology in a materials-relevant model, with concrete diagnostics (Airy scaling, sub-Page entanglement bands) for neutron-scattering and quench experiments. The main limitations are that the 'parameter-free' test is weakened by an n=1 energy alignment, the W=2 state selection is not fully specified, and the level-statistics crossover is supported by only three system sizes and no error bars.

major comments (4)
  1. [Sec. VI B, Eq. (6.11), Fig. 6(a)] The comparison is presented as parameter-free, but in Fig. 6(a) the ED levels are shifted by a single constant so that n=1 matches theory, which removes the absolute offset E_cont(0)+f a_+(0) from Eq. (6.11). What remains is a test of relative Airy spacings and of the y=z functional form in Fig. 6(c), not of the full absolute prediction. The manuscript should state this limitation explicitly, or better, also report the unshifted ED n=1 energy against E_cont(0)+f a_+(0), and/or use the finite-N PBC version of Rutkevich's theory that is acknowledged in Sec. VI A.
  2. [Sec. VI A/VI B, Fig. 6(b)] The comparison uses thermodynamic-limit formulas (Eqs. (6.8)-(6.11)) while a finite-N PBC version is acknowledged but never employed. The need to account for finite-size effects is already visible: residuals grow with n, and the stable-meson count is n*_ED=8 versus n*_th=10. These discrepancies are attributed to finite-size and O(f^{4/3}) corrections, but without quantifying them the 'close quantitative agreement' is confined to the lowest few levels. Please either use the finite-N expressions or provide a quantitative estimate of the 1/N and subleading corrections in the comparison window.
  3. [Sec. VI B, Sec. V, Eq. (4.3)] The selection of W=2 one-meson states via the correlation measure C_zz ≈ +4 is not specified with a tolerance or a counting procedure. Since this same observable is used to demonstrate banding in Sec. V, the ordered ED sequence compared with theory in Figs. 6, 9-11 is not independently reproducible from the text. The authors should state the selection window, the number of states per band, and how the W=2 states are separated from higher-W branches.
  4. [Sec. III, Figs. 2-3, Table III] The level-statistics evidence for a crossover to nonergodic behavior is weaker than the text suggests. At the representative deep-AF point (Δ=-4.5, h=0.3, N=22), ⟨r⟩≈0.47, which is intermediate between GOE (≈0.536) and Poisson (≈0.386); there are no error bars on ⟨r⟩. The thermodynamic crossover Δ_tra(∞)=−4.19±0.115 is obtained from a linear fit to only three system sizes (N=18,20,22) using an ad hoc crossover criterion. Please provide bootstrap or binning errors, and a more robust finite-size scaling analysis, before claiming a thermodynamic crossover to Poisson-like statistics.
minor comments (4)
  1. [Sec. V] Typo: 'degrees fo freedom' should be 'degrees of freedom'.
  2. [Fig. 3 caption] The caption gives Δ_tra(∞)=−4.187±0.115 while the text states −4.19±0.115; use consistent significant figures.
  3. [Sec. VII] The interpretation in terms of weak fragmentation and scar-like dynamics is inferential; no dynamical probe (quench, revival, autocorrelation) is presented. Consider softening the wording or adding a short dynamical diagnostic to support the 'long-lived subspaces' claim.
  4. [Appendix C, Figs. 9-11] The residual and deviation panels would benefit from error bars or finite-size comparisons; currently it is difficult to judge whether the deviations at larger n are statistically meaningful.

Circularity Check

0 steps flagged

No significant circularity: the analytic meson ladder is an external, parameter-free theory and the ED comparisons are transparently acknowledged as offset-aligned and low-order fits.

full rationale

Section VI's analytic ladder is taken from Rutkevich's external theory (Ref. [50]) and evaluated with no free parameters; E_cont, A, and a_+ are computed from Δ, h, and J via Eqs. (6.8)-(6.14). The ED data are independent numerical spectra in a fixed symmetry sector. Although Fig. 6(a) aligns the n=1 ED level to theory and Fig. 6(b) performs a two-parameter linear fit over the lowest three levels, both procedures are disclosed explicitly ('the ED series is shifted so that the first ED level (n=1) coincides with the analytic value'; 'fit parameters determined by the lowest three z_n points'), and the claimed parameter-free line is the unadjusted theory, not the fit. The offset-removed Airy collapse in Eq. (6.27) is for the theory an identity by construction, as the caption states ('theory lies on y=z by construction'), but the ED collapse is a genuine test of Airy zero spacings; the same is true of the spacing comparison, which cancels E_cont and f a_+. The two-meson threshold is either computed from the same external analytic ladder or from the lowest W=4 ED level, and the near-coincidence is an independent cross-check, not an input. The C_zz ≈ +4 criterion is a state-selection convention; the spectroscopy comparison itself uses ED energies, not the same C_zz values as predicted outputs. The self-citations (Refs. [25]-[29]) motivate the fragmentation/scar interpretation but are not used to derive Eq. (6.11). Finite-size limitations and the n* = 8 vs 10 difference are acknowledged as deviations. No step in the derivation chain reduces to its own inputs, so the paper is not circular.

Axiom & Free-Parameter Ledger

3 free parameters · 6 axioms · 1 invented entities

The analytic meson-ladder comparison is genuinely parameter-free in its construction (A(0), E_cont, f a_+ all computed from Δ, h, J), which is why the circularity burden is low. The fitted quantities are the comparison protocol's alignment constant, the (A_fit, c_fit) decomposition, and the hand-chosen crossover criterion for Δ_tra(N). The main invented quantity is the quasi-conserved W, evidenced internally. The finite-size-to-thermodynamic-limit transfer and the integrability confound are the sharpest unproven premises.

free parameters (3)
  • global energy alignment ε_align (per parameter set) = not tabulated; applied so that the n=1 ED meson level coincides with the analytic value
    Fig. 6(a) and App. C: ED energies are shifted by a single constant before level-by-level comparison; this removes the absolute offset from the test of the theory's constant term.
  • (A_fit, c_fit) — linear fit of ED bindings vs Airy magnitudes = not tabulated in text; reported as 'differ modestly' from (A(0), f a_+(0))
    Fig. 6(b) and App. C: two-parameter linear fit over the lowest three z_n points, used to decompose residuals into slope and intercept parts; disclosed as a fit, not a prediction.
  • crossover criterion defining Δ_tra(N) = Δ_tra(18), Δ_tra(20), Δ_tra(22) used in Fig. 3; values not tabulated
    Sec. III: 'values at which ⟨r⟩ begins to deviate from the GOE value' is a hand-chosen, not algorithmically specified, criterion; the resulting 2-parameter linear fit through three points yields Δ_tra(∞) = −4.19 ± 0.115 with no uncertainties on the input points.
axioms (6)
  • domain assumption The h=0 spinon dispersion ω(q)=I√(1−k²cos²q) with the elliptic-modulus relations (B2)–(B4), and the string tension f=2hσ̄(η) with σ̄ the spontaneous-magnetization product (B13), correctly describe the h≠0 staggered-field chain near threshold.
    Taken from Rutkevich [50] (Bethe ansatz); adopted as inputs in Sec. VI A. The assumption that the deconfined two-spinon kinematics survives the staggered field as a linear confining potential is the physical content being tested.
  • standard math Airy quantization of a linear-potential bound state: E_n = E_cont + A z_n + f a_± + O(f^{4/3}) (Eq. 6.11), with A ∝ f^{2/3}.
    Standard WKB/Airy result for a linear potential; the O(f^{4/3}) corrections are neglected in the comparison and invoked to explain residuals.
  • domain assumption The '+' (symmetric/C_flip=+1) branch is the one realized in the (S_z, P, I, C_flip) = (0,0,+1,+1) sector.
    Assumed in Sec. VI A; if the wrong meson parity were identified, the comparison would be off by the channel shift difference f(a_+ − a_−) ≈ 2f ln2/π.
  • domain assumption Finite-size (N=22, PBC) ED data can be compared with thermodynamic-limit analytic formulas, with 1/N and momentum-quantization corrections subleading for the lowest levels.
    Stated in Sec. VI and App. B; the n* = 8 vs 10 discrepancy and the growing residuals with n show the assumption is imperfect.
  • standard math Standard elliptic-function, theta-function, digamma, and product identities used in App. B.4–B.5 to evaluate a_±(0) = −γ_E/π and −(γ_E + 2 ln2)/π.
    Unproved background results, cited to NIST DLMF and Whittaker–Watson; standard within the frame of the paper.
  • domain assumption The level-statistics crossover at large |Δ| is caused by confinement (quasi-conserved W) rather than by proximity to the integrable h=0 XXZ point.
    Interpretive premise in Sec. III; the paper does not directly test the confound that deep in the AF phase the model is a small perturbation of the integrable XXZ chain.
invented entities (1)
  • emergent quasi-conserved domain-wall number W no independent evidence
    purpose: Labels eigenstate bands and explains suppressed entanglement, band-resolved correlations, and weakened level repulsion in the AF regime
    Domain walls are standard in spin-chain physics, but the quasi-conservation of W in this model is established only through the paper's own diagnostics (C_zz banding, ⟨r⟩, S_vN); no falsifiable handle outside the paper is supplied beyond the model itself.

pith-pipeline@v1.3.0-alltime-deepseek · 28483 in / 25118 out tokens · 252390 ms · 2026-08-03T21:33:25.894664+00:00 · methodology

0 comments
read the original abstract

Confinement of fractionalized excitations can strongly restructure many-body spectra. We investigate this phenomenon in the gapped spin-$\frac{1}{2}$ XXZ chain subject to a staggered field, where spinons bind into domain-wall ``mesons'' deep in the antiferromagnetic phase. We present evidence that this non-integrable model exhibits both Hilbert space fractionalization and quantum scar formation as controlled by the anisotropy parameter $\Delta$. Exact diagonalization across symmetry-resolved sectors reveals a crossover from Gaussian-orthogonal (chaotic) level statistics at weak anisotropy $\Delta \sim 1$ to non-ergodic behavior deep in the antiferromagnetic regime $\Delta\gg 1$ through scrutinizing the adjacent gap ratios, accompanied by a striking banding of eigenstates by domain-wall number in correlation and entanglement measures. The Page-like entanglement dome characteristic of chaotic spectra gives way to suppressed, band-resolved entanglement consistent with emergent quasi-conservation of domain walls. To investigate further the formation mechanism of mesonic scar states, we carry out meson spectroscopy near the two-spinon threshold and compare with the analytic ladder predicted by Rutkevich [Phys. Rev. B 106, 134405 (2022)]. We test the theory through continuum-relative bindings, an offset-removed Airy scaling collapse, and explicit two-meson thresholds that determine the number of stable meson levels. The low-lying spectrum shows close quantitative agreement, while deviations at higher energies are consistent with finite-size and subleading corrections. These results establish a unified account of confinement-induced nonergodicity and provide a template for quantitative meson spectroscopy in quantum spin chains.

Figures

Figures reproduced from arXiv: 2511.14847 by Julia Wildeboer, Marton Lajer, Robert M. Konik.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) shows the finite-size dependence of ⟨r⟩ as a function of the anisotropy factor ∆ for system sizes N = h ⟨r⟩ (N=20) ⟨r⟩ (N=22) 0.1 0.437191 0.443036 0.2 0.428236 0.472977 0.3 0.452605 0.496836 0.4 0.455710 0.500026 0.5 0.456062 0.505295 TABLE III. Mean adjacent–gap ratio ⟨r⟩ in the fixed domain-wall sector W = 10 for ∆ = −4.5 within the symmetry block (S z total, P, I, Cflip) = (0, 0, +1, +1). Values ar… view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: (b), defined by subtracting the edge: bn ≡ En − Econt(0). The parameter–free prediction b th n = A(0) zn + f a+(0) is shown together with a linear regression of the ED bindings b ED n ≈ Afit zn + cfit (fit performed versus the Airy magnitudes zn over the indicated low–lying window n ≤ 3). The lower panel displays the binding residual δ (b) n ≡ b ED n − b th n = (Afit − A(0)) zn + [PITH_FULL_IMAGE:figures/… view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p020_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11 [PITH_FULL_IMAGE:figures/full_fig_p020_11.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

61 extracted references · 1 linked inside Pith

  1. [1]

    J. M. Deutsch, Quantum statistical mechanics in a closed sys- tem, Phys. Rev. A43, 2046 (1991)

    2046

  2. [2]

    Srednicki, Chaos and quantum thermalization, Phys

    M. Srednicki, Chaos and quantum thermalization, Phys. Rev. E 50, 888 (1994)

    1994

  3. [3]

    D. M. Basko, I. L. Aleiner, and B. L. Altshuler, Metal–insulator transition in a weakly interacting many-electron system with localized single-particle states, Annals of Physics321, 1126 (2006)

    2006

  4. [4]

    Nandkishore and D

    R. Nandkishore and D. A. Huse, Many-body localization and thermalization in quantum statistical mechanics, Annual Re- view of Condensed Matter Physics6, 15 (2015)

    2015

  5. [5]

    Khemani, M

    V . Khemani, M. Hermele, and R. Nandkishore, Localization from hilbert space shattering: From theory to physical realiza- tions, Phys. Rev. B101, 174204 (2020)

    2020

  6. [6]

    L. Zhao, P. R. Datla, W. Tian, M. M. Aliyu, and H. Loh, Ob- servation of quantum thermalization restricted to hilbert space fragments and𭟋 2k scars, Phys. Rev. X15, 011035 (2025)

    2025

  7. [7]

    P. Sala, T. Rakovszky, R. Verresen, M. Knap, and F. Pollmann, Ergodicity breaking arising from hilbert space fragmentation in dipole-conserving hamiltonians, Phys. Rev. X10, 011047 (2020)

    2020

  8. [8]

    Moudgalya, B

    S. Moudgalya, B. A. Bernevig, and N. Regnault, Quantum many-body scars and hilbert space fragmentation: a review of exact results, Rep. Prog. Phys.85, 086501 (2022)

    2022

  9. [9]

    Z.-C. Yang, F. Liu, A. V . Gorshkov, and T. Iadecola, Hilbert- space fragmentation from strict confinement, Phys. Rev. Lett. 124, 207602 (2020)

    2020

  10. [10]

    Kohlert, S

    T. Kohlert, S. Scherg, P. Sala, F. Pollmann, B. H. Madhusud- hana, I. Bloch, and M. Aidelsburger, Exploring the regime of fragmentation in strongly tilted fermi-hubbard chains, Phys. Rev. Lett.130, 010201 (2023)

    2023

  11. [11]

    L. Zhao, P. R. Datla, W. Tian, M. M. Aliyu, and H. Loh, Ob- servation of quantum thermalization restricted to hilbert space fragments and2Kscars, Phys. Rev. X15, 011035 (2025)

    2025

  12. [12]

    Khemani, M

    V . Khemani, M. Hermele, R. Nandkishore, and R. Moessner, Localization from hilbert space fragmentation, Nature Physics 15, 1231 (2019)

    2019

  13. [13]

    Moudgalya, E

    S. Moudgalya, E. J. Bergholtz, N. Regnault, and S. L. Sondhi, Quantum many-body scars and hilbert space fragmentation: A review of exact results, Reviews of Modern Physics94, 045005 (2022)

    2022

  14. [14]

    Chandran and V

    A. Chandran and V . Khemani, Quantum many-body scars and hilbert space fragmentation: Past, present, and future, Annual Review of Condensed Matter Physics14, 443 (2023)

    2023

  15. [15]

    Bernien, S

    H. Bernien, S. Schwartz, A. Keesling, H. Levine, A. Omran, H. Pichler, S. Choi, A. S. Zibrov, M. Endres, M. Greiner, V . Vuleti´c, and M. D. Lukin, Probing many-body dynamics on a 51-atom quantum simulator, Nature551, 579 (2017)

    2017

  16. [16]

    Bluvstein, A

    D. Bluvstein, A. Omran, H. Levine, A. Keesling, G. Semegh- ini, T. Wang, S. Ebadi, H. Bernien, A. S. Zibrov, H. Pichler, S. Choi, J. Cui, M. Rossignolo, M. Kalinowski, J. Kolder, L.-P. Henry, A. Vishwanath, and M. D. Lukin, Controlling quantum many-body dynamics in driven rydberg atom arrays, Science 371, 1355 (2021)

    2021

  17. [17]

    C. J. Turner, A. A. Michailidis, D. A. Abanin, M. Serbyn, and Z. Papi´c, Weak ergodicity breaking from quantum many-body scars, Nature Physics14, 745 (2018)

    2018

  18. [18]

    C. J. Turner, A. A. Michailidis, D. A. Abanin, M. Serbyn, and Z. Papi´c, Quantum scarred eigenstates in a rydberg atom chain: Entanglement, breakdown of thermalization, and stability to perturbations, Phys. Rev. B98, 155134 (2018)

    2018

  19. [19]

    S. Choi, C. J. Turner, H. Pichler, W. W. Ho, A. A. Michailidis, Z. Papi´c, M. Serbyn, M. D. Lukin, and D. A. Abanin, Emergent su(2) dynamics and perfect quantum many-body scars, Phys. Rev. Lett.122, 220603 (2019)

    2019

  20. [20]

    W. W. Ho, S. Choi, H. Pichler, and M. D. Lukin, Periodic orbits, entanglement, and quantum many-body scars in constrained models: Matrix product state approach, Phys. Rev. Lett.122, 040603 (2019)

    2019

  21. [21]

    Shiraishi and T

    N. Shiraishi and T. Mori, Systematic construction of counterex- amples to the eigenstate thermalization hypothesis, Physical Review Letters119, 030601 (2017)

    2017

  22. [22]

    At more negative∆, all sizes continue to trend downward toward Poisson-like statistics, consistent with confinement-induced nonergodicity

    the GOE plateau persists up to roughly∆≈ −2.5, beyond which a clear crossover toward non-GOE values sets in. At more negative∆, all sizes continue to trend downward toward Poisson-like statistics, consistent with confinement-induced nonergodicity. Across all panels (a)-(c), these field-resolved trends dovetail with the correlation/entanglement “banding” d...

  23. [23]

    Schecter and T

    M. Schecter and T. Iadecola, Weak ergodicity breaking and quantum many-body scars in spin-1 xy magnets, Physical Re- view Letters123, 147201 (2019)

    2019

  24. [24]

    Moudgalya, S

    S. Moudgalya, S. Rachel, B. A. Bernevig, and N. Regnault, Ex- act excited states of nonintegrable models, Physical Review B 98, 235155 (2018)

    2018

  25. [25]

    S. B. Rutkevich, Energy spectrum of bound-spinons in the quantum ising spin chain with a longitudinal magnetic field, Journal of Statistical Mechanics: Theory and Experiment , P07015 (2010)

    2010

  26. [26]

    A. J. A. James and R. M. Konik, Quantum dynamics of excita- tions in strongly confined ising spin chains, Physical Review B 87, 241103 (2013)

    2013

  27. [27]

    N. J. Robinson, A. J. A. James, and R. M. Konik, Signatures of rare states and thermalization in a theory with confinement, Phys. Rev. B99, 195108 (2019)

    2019

  28. [28]

    A. J. A. James, R. M. Konik, and N. J. Robinson, Nonther- mal states arising from confinement in one and two dimensions, Phys. Rev. Lett.122, 130603 (2019)

    2019

  29. [29]

    H. Kim, T. N. Ikeda, and D. A. Huse, Testing whether all eigen- states obey the eigenstate thermalization hypothesis, Phys. Rev. E90, 052105 (2014)

    2014

  30. [30]

    A. J. A. James, R. M. Konik, P. Lecheminant, N. J. Robin- son, and A. M. Tsvelik, Non-perturbative methodologies for low-dimensional strongly-correlated systems: From non- abelian bosonization to truncated spectrum methods, Reports on Progress in Physics81, 046002 (2018)

    2018

  31. [31]

    Moudgalya, N

    S. Moudgalya, N. Regnault, and B. A. Bernevig, Exact excited states of nonintegrable models and quantum many-body scars, Physical Review B98, 235156 (2018)

    2018

  32. [32]

    Affleck, T

    I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Rigorous re- sults on valence-bond ground states in antiferromagnets, Phys- ical Review Letters59, 799 (1987)

    1987

  33. [33]

    O’Dea, F

    N. O’Dea, F. J. Burnell, A. Chandran, and V . Khemani, From quantum chaos to localization transition in the spin-1 aklt chain, Physical Review Letters124, 040605 (2020). 19

    2020

  34. [34]

    L. L. Kish, A. Weichselbaum, D. M. Pajerowski, A. T. Savici, A. Podlesnyak, L. Vasylechko, A. Tsvelik, R. Konik, and I. A. Zaliznyak, High-temperature quantum coherence of spinons in a rare-earth spin chain, Nature Communications16, 10.1038/s41467-025-61715-z (2025)

    doi:10.1038/s41467-025-61715-z 2025

  35. [35]

    L. Wu, S. E. Nikitin, M. Brando, L. Vasylechko, G. Ehlers, M. Frontzek, A. T. Savici, G. Sala, A. D. Christianson, M. D. Lumsden,et al., Antiferromagnetic ordering and dipolar inter- actions of ybalo 3, Physical Review B99, 195117 (2019)

    2019

  36. [36]

    L. Wu, S. Nikitin, Z. Wang, W. Zhu, C. D. Batista, A. Tsvelik, A. M. Samarakoon, D. A. Tennant, M. Brando, L. Vasylechko, et al., Tomonaga–luttinger liquid behavior and spinon confine- ment in ybalo3, Nature communications10, 698 (2019)

    2019

  37. [37]

    Breunig, M

    O. Breunig, M. Garst, E. Sela, B. Buldmann, and P. e. a. Becker, The spin- 1 2 xxz chain system cs2cocl4 in a transverse magnetic field, Phys. Rev. Lett.111, 187202 (2013)

    2013

  38. [38]

    Coldea, D

    R. Coldea, D. A. Tennant, E. M. Wheeler, E. Wawrzynska, D. Prabhakaran, M. T. F. Telling, K. Habicht, P. Smeibidl, and K. Kiefer, Quantum criticality in an ising chain: Experimental evidence for emergente 8 symmetry, Science327, 177 (2010)

    2010

  39. [39]

    H. Zou, Y . Cui, X. Wang, Z. Zhang, J. Yang, G. Xu, A. Okutani, M. Hagiwara, M. Matsuda, G. Wang, G. Mussardo, K. H´ods´agi, M. Kormos, Z. He, S. Kimura, R. Yu, J. Ma, and J. Wu, E8 spec- tra of quasi-one-dimensional antiferromagnet baco2v2o8 under transverse field, Phys. Rev. Lett.127, 077201 (2021)

    2021

  40. [40]

    Z. Wang, T. Lorenz, D. I. Gorbunov, P. T. Cong, Y . Kohama, S. Niesen, O. Breunig, J. Engelmayer, A. Hermann, J. Wu, K. Kindo, J. Wosnitza, S. Zherlitsyn, and A. Loidl, Quantum criticality of an ising-like spin- 1 2 antiferromagnetic chain in a transverse magnetic field, Phys. Rev. Lett.120, 207205 (2018)

    2018

  41. [41]

    Faure, S

    Q. Faure, S. Takayoshi, S. Petit, V . Simonet, S. Raymond, L.- P. Regnault, M. Boehm, J. S. White, M. M ˚ansson, C. R ¨uegg, P. Lejay, B. Canals, T. Lorenz, S. C. Furuya, T. Giamarchi, and B. Grenier, Topological quantum phase transition in the ising- like antiferromagnetic spin chain baco 2v2o8, Nature Physics 14, 716 (2018)

    2018

  42. [42]

    Matsuda, H

    M. Matsuda, H. Onishi, A. Okutani, J. Ma, P. Agrawal, T. Hong, D. M. Pajerowski, J. R. D. Copley, K. Okunishi, M. Mori, S. Kimura, and M. Hagiwara, Magnetic structure and disper- sion relation of the quasi-one-dimensional ising-like antiferro- magnet baco2v2o8 in a transverse magnetic field, Phys. Rev. B 96, 024439 (2017)

    2017

  43. [43]

    Kimura, H

    S. Kimura, H. Yashiro, K. Okunishi, M. Hagiwara, Z. He, K. Kindo, T. Taniyama, and M. Itoh, Field-induced order- disorder transition in antiferromagnetic baco 2v2o8 driven by a softening of spinon excitations, Phys. Rev. Lett.99, 087602 (2007)

    2007

  44. [44]

    Grenier, S

    B. Grenier, S. Petit, V . Simonet, E. Can ´evet, L.-P. Regnault, S. Raymond, B. Canals, C. Berthier, and P. Lejay, Longitudinal and transverse zeeman ladders in the ising-like chain antiferro- magnet baco2v2o8, Phys. Rev. Lett.114, 017201 (2015)

    2015

  45. [45]

    A. K. Bera, B. Lake, F. H. L. Essler, L. Vanderstraeten, C. Hu- big, U. Schollw ¨ock, A. T. M. N. Islam, A. Schneidewind, and D. L. Quintero-Castro, Spinon confinement in a quasi-one- dimensional anisotropic heisenberg magnet, Phys. Rev. B96, 054423 (2017)

    2017

  46. [46]

    Z. Wang, M. Schmidt, A. K. Bera, A. T. M. N. Islam, B. Lake, A. Loidl, and J. Deisenhofer, Spinon confinement in the one- dimensional ising-like antiferromagnet srco 2v2o8, Phys. Rev. B91, 140404(R) (2015)

    2015

  47. [47]

    D’Alessio, Y

    L. D’Alessio, Y . Kafri, A. Polkovnikov, and M. Rigol, From quantum chaos and eigenstate thermalization to statistical me- chanics and thermodynamics, Advances in Physics65, 239 (2016), https://doi.org/10.1080/00018732.2016.1198134

    arXiv 2016

  48. [48]

    Oganesyan and D

    V . Oganesyan and D. A. Huse, Localization of interacting fermions at high temperature, Phys. Rev. B75, 155111 (2007)

    2007

  49. [49]

    (4.1) without the alternating factor(−1) i is zero since we work in theS z total = 0sector

    The average longitudinal magnetization defined as Eq. (4.1) without the alternating factor(−1) i is zero since we work in theS z total = 0sector

  50. [50]

    The energy shiftE shift is introduced to avoid the data belonging to the different values of anisotropy∆(partially) lying on top of each other

  51. [51]

    S. B. Rutkevich, Spinon confinement in the gapped antiferro- magnetic xxz spin- 1 2 chain, Phys. Rev. B106, 134405 (2022)

    2022

  52. [52]

    For zero staggered fieldh= 0the infinite-chain Hamiltonian is written with aJ/2prefactor and sums overj∈Z; a finite-N version with PBC is also given and reduces to the same results in theN→ ∞limit

  53. [53]

    Fonseca and A

    P. Fonseca and A. Zamolodchikov, Ising field theory in a mag- netic field: Analytic properties of the free energy, (2001), arXiv:hep-th/0112167

    Pith/arXiv arXiv 2001

  54. [54]

    6(b) use the same reference Econt(0)

    All thresholds and bindings in Fig. 6(b) use the same reference Econt(0). This avoids ambiguities when comparing ED and the- ory lines

  55. [55]

    Rutkevich’s conven- tion employs a coupling constant ofJ/2instead ofJand spe- cific expressions need to be adjusted to accommodate for that choice

    We stress that the expression are valid for our choice of cou- pling constantJin the Hamiltonian (2.1). Rutkevich’s conven- tion employs a coupling constant ofJ/2instead ofJand spe- cific expressions need to be adjusted to accommodate for that choice

  56. [56]

    gov/(2010), online companion to the NIST Handbook of Mathematical Functions

    NIST Digital Library of Mathematical Functions, NIST Digital Library of Mathematical Functions,https://dlmf.nist. gov/(2010), online companion to the NIST Handbook of Mathematical Functions

    2010

  57. [57]

    E. T. Whittaker and G. N. Watson,A Course of Modern Anal- ysis, 4th ed. (Cambridge University Press, Cambridge, 1996) reprint of the 4th edition (1927)

    1996

  58. [58]

    Abramowitz and I

    M. Abramowitz and I. A. Stegun, eds.,Handbook of Mathe- matical Functions with Formulas, Graphs, and Mathematical Tables(Dover, New York, 1965)

    1965

  59. [59]

    F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, eds.,NIST Handbook of Mathematical Functions(Cambridge University Press, New York, 2010)

    2010

  60. [60]

    F. W. J. Olver,Asymptotics and Special Functions(AK Peters, Wellesley, MA, 1997) reprint of the 1974 Academic Press edi- tion

    1997

  61. [61]

    Gasper and M

    G. Gasper and M. Rahman,Basic Hypergeometric Series, 2nd ed., Encyclopedia of Mathematics and its Applications, V ol. 96 (Cambridge University Press, Cambridge, 2004). 20 FIG. 9.ED vs analyticW= 2meson spectrum atP= 0forN= 22,∆ =−4.5,h= 0.05(additional data).(a) Aligned absolute energies: ED levels are shifted by a single constant so that then=1meson coin...

    2004

This paper was first reviewed by deepseek-v4-flash on August 3, 2026.