Pith. sign in

REVIEW 3 major objections 4 minor 17 references

Adiabatic dynamics of quasiperiodic transverse Ising model

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read A quasiperiodic transverse Ising chain obeys Kibble-Zurek scaling across both of its critical lines, with defect densities decaying as $\tau^{-1/2}$ and $\tau^{-1/3}$.

desk verdict The Ising-line half of this Kibble–Zurek study is credible, but the new tau^{-1/3} claim for the z=2 critical line is under-supported by a factor-of-four tau window and no fitted exponent. read the letter →

arxiv 1908.01959 v1 pith:3SDG34J7 submitted 2019-08-06 cond-mat.stat-mech quant-ph

Revathy B. S. , Uma Divakaran This is my paper · ORCID
classification cond-mat.stat-mechquant-ph
keywords Kibble-ZurekscalingquasiperiodictransverseIsingmodelquantumphasetransitiondefectdensityuniversalityclassdualityincommensuratemodulationcriticalexponents
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper asks whether the Kibble-Zurek scaling law for defects generated by a slow sweep through a quantum critical point survives when the Hamiltonian is quasiperiodic rather than clean or disordered. It studies a transverse Ising chain with an incommensurate modulation, which has two distinct critical lines belonging to two different universality classes, and it claims that both lines obey the standard power law. For the Ising-type line the defect density falls as $n \sim \tau^{-1/2}$; for the second line the paper finds $n \sim \tau^{-1/3}$, matching the exponent predicted from the equilibrium critical exponents $\nu=1$, $z=2$. The authors verify both predictions by numerically evolving finite chains through two-step linear sweeps, and the agreement improves as the system size grows. If correct, this extends Kibble-Zurek scaling to quasiperiodic criticality, where excitations can be extended, localized, or multifractal.

What carries the argument

The load-bearing object is the Kibble-Zurek formula $n \sim \tau^{-\nu d/(\nu z + 1)}$, which ties the defect density to the equilibrium correlation-length and dynamical exponents of the critical point crossed. The numerical test is carried by a two-step quench protocol: the first step linearly sweeps one Hamiltonian parameter across the critical line at rate $1/\tau$, and the second step ramps the remaining quasiperiodic field to zero, leaving a clean Ising Hamiltonian in which defects are simply domain walls counted by $n = (1/L)\sum_i \langle \psi_f | \tfrac{1}{2}(1-\sigma_i^x \sigma_{i+1}^x) | \psi_f \rangle$. The mapping that connects the simulated plane to the known phase diagram is the Ising duality, which swaps the paramagnetic and ferromagnetic phases and moves the quasiperiodic modulation from the coupling $J$ to the field $h$ while, the paper asserts, leaving the dynamical nature of bulk single-particle excitations unchanged.

What would settle it

Run the same two-step linear sweeps directly in the original $J$-quasiperiodic Hamiltonian of Eq. (1), crossing both critical lines, and count domain-wall defects in the final clean Ising Hamiltonian; if the exponents are not $1/2$ and $1/3$, the duality mapping or the inherited exponents are wrong. Alternatively, compute the dynamical exponent of line B from the gap scaling in the $h$-plane and check whether $z=2$ holds exactly.

Watch

Extended reading notes

Core claim

The central discovery is that Kibble-Zurek scaling, $n \sim \tau^{-\nu d/(\nu z + 1)}$, holds separately for the two critical lines of the quasiperiodic transverse Ising model. Crossing the Ising-type critical line A, with $\nu=1$ and $z=1$, gives $n \sim \tau^{-1/2}$. Crossing the new critical line B, with $\nu=1$ and $z=2$, gives $n \sim \tau^{-1/3}$. The paper obtains these results in the $h$-$A_h$ plane, where the quasiperiodic modulation acts on the transverse field, and it uses an Ising duality to carry over the phase diagram and exponents established in the $J$-$A_J$ plane. Defects are counted as domain walls after a second sweep brings the system to a clean Ising Hamiltonian, and the numerics show power-law decay with exponents approaching the predicted values as $N$ increases.

Load-bearing premise

The prediction relies on the assumption that the Ising duality exactly maps the original $J$-quasiperiodic model to the numerically simulated $h$-quasiperiodic model and preserves the dynamical exponent of bulk single-particle excitations, so that the line-B value $z=2$ carries over unchanged.

Editorial extensions

If this is right

  • A slow sweep across the Ising-type critical line leaves $n \sim \tau^{-1/2}$ defects, confirming that quasiperiodicity does not alter the quantum Ising universality class.
  • A slow sweep across the new line leaves $n \sim \tau^{-1/3}$, which is a slower decay and therefore a larger residual defect density than at an Ising critical point for the same sweep rate.
  • The two-step protocol provides a practical way to extract defect densities from a final clean Ising Hamiltonian, making the quasiperiodic dynamics measurable through domain-wall counting.
  • Finite-size numerics approach the asymptotic exponents as the chain length increases, consistent with the equilibrium exponents $\nu=1$, $z=1$ and $\nu=1$, $z=2$ used in the prediction.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A direct simulation of the same sweeps in the original $J$-quasiperiodic Hamiltonian would test whether the Ising duality is exact for incommensurate modulations; if the exponents differ, the discrepancy would pinpoint where the mapping fails.
  • Because line B has $z=2$, the paper indirectly predicts how the defect exponent should change under nonlinear sweeps of the form $t^\alpha$, providing a sharper experimental test of the line's universality class than a single linear sweep.
  • The defect density alone may not distinguish localized, extended, or multifractal excitations; a natural extension would be to examine the full distribution of defects, which could carry signatures of the quasiperiodic eigenstate structure.
  • An optical-lattice realization of the quasiperiodic transverse Ising chain could measure the $\tau^{-1/3}$ exponent directly, separating line B from ordinary Ising criticality in an experiment.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper studies the linear quench dynamics of a one-dimensional quasiperiodic transverse Ising model, using two-step protocols in the h-quasiperiodic representation (Eq. 2) to cross two distinct critical lines. It claims that the defect density obeys Kibble-Zurek scaling with n ~ tau^{-1/2} for the Ising-type line A and n ~ tau^{-1/3} for the new line B, where the exponents are taken from the equilibrium analysis of Ref. 18. The paper concludes that this confirms both the KZ scaling and the equilibrium critical exponents for the quasiperiodic model.

Significance. If quantitatively established, the tau^{-1/3} result would be a valuable extension of KZ scaling to a quasiperiodic model with a new universality class. The paper's test is not circular by construction: the equilibrium exponents are independent inputs, and the defect densities are measured separately. The line-A data over roughly two decades are a useful confirmation. However, the central novel claim (line B) rests on a very narrow tau window, so the paper in its current form does not yet establish the exponent. The stress-test concern about circularity does not land, because the exponents are inputs from Ref. 18, not fitted here; the concern about the narrow line-B window does land. The duality mapping is standard and likely exact, but it should be made explicit to support the transfer of exponents from Eq. 1 to Eq. 2.

major comments (3)
  1. [§III, Fig. 3] The claimed tau^{-1/3} scaling for line B is not quantitatively established. The data cover only tau in [50, 200], a factor of 4, and no fitted exponent, error bars, or finite-size extrapolation are reported. Over such a short interval, a visually drawn tau^{-0.33} reference line cannot distinguish exponents differing by 0.1-0.2, and the largest-tau points are few. The statement that agreement with the exponent improves as the system size increases is asserted but not quantified. A local-slope analysis (binned log-log slopes as a function of tau and N) or a scaling collapse, together with confidence intervals, is required to support the central novel claim.
  2. [§III, first paragraph and Eq. (2)] The KZ predictions are evaluated for the h-quasiperiodic Hamiltonian in Eq. (2), but the exponents nu=1,z=1 and nu=1,z=2 are taken from the J-quasiperiodic Hamiltonian in Eq. (1) of Ref. 18. The paper asserts without derivation that the Ising duality swaps the phases and leaves the dynamics of bulk single-particle excitations unaltered. The standard Kramers-Wannier duality is expected to be exact for site-dependent couplings, so this is not a circularity problem; nevertheless, the mapping should be written explicitly for the incommensurate modulation, or the equilibrium gap scaling should be checked directly in Eq. (2), before the tau^{-1/3} prediction can be regarded as properly tested.
  3. [§III, Fig. 2] The line-A data are consistent with tau^{-1/2} over roughly two decades, which is encouraging, but the claim that deviations at large tau are due to finite-size effects is not supported by any quantitative comparison, such as a finite-size gap estimate or a scaling collapse. Since the paper's stated goal is to verify KZ scaling, a fit with confidence intervals for the line-A exponent should also be reported rather than relying on visual inspection.
minor comments (4)
  1. [§III, protocols] The two-step protocols are described verbally, but the time dependence is ambiguous (e.g., 'reduce Ah as t/tau from 3 to zero'), and the value of J used in the simulations is not stated. Explicit expressions for h(t) and Ah(t), the total sweep time, and the parameter values would make the calculations reproducible.
  2. [References] There are several typographical and bibliographic errors, including 'Chkarabarti' in Refs. 2 and 3, 'P J. D. Crowley' in Ref. 20, and missing volume/page information in Ref. 17; Ref. 16 lists Phys. Rev. Lett. 144, 083002 (2015), which should likely be volume 114.
  3. [Figs. 2-3] The reference lines in Figs. 2 and 3 appear to be guides rather than fits; the figure captions should state this explicitly and, ideally, report the fitted exponents and their uncertainties.
  4. [Abstract and Conclusion] The abstract and conclusion state that the power-law behavior is 'indeed obeyed' and 'clearly' shown; given the limited tau window for line B, this wording overstates the quantitative evidence and should be tempered.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the KZ exponents are independent equilibrium inputs, and defect densities are measured directly.

full rationale

The paper imports the equilibrium critical exponents nu=1, z=1 for line A and nu=1, z=2 for line B from Ref. 18, a separate equilibrium study, and then independently measures the defect density n as a function of tau by solving the Schrodinger equation for the quench protocols. Since the exponents are not fitted parameters and the defect densities are computed directly from the evolved wavefunction via Eq. 3, the resulting n versus tau curves provide an independent test of the Kibble-Zurek relation rather than reducing to the inputs by construction. The asserted Ising duality between Eq. 1 and Eq. 2 is an assumption about universality, not a circular step, and the self-citations (e.g., Ref. 19) are background references that do not carry the argument. The paper's central verification is therefore self-contained as a numerical test, with no fitted-input-called-prediction or self-citation chain forcing the result.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central claim (defect-density exponents 1/2 and 1/3) rests on four inputs the paper does not itself establish: the duality between the two quasiperiodic planes, the equilibrium exponents inherited from Ref. 18, the ground-state initialization, and an unstated free-fermion numerical scheme. No new entities are postulated and no parameter is fitted to data; all adjustable numbers are quench-protocol choices that do not enter the exponent.

free parameters (2)
  • Quench protocol endpoints (h0, Ah values) = Line A: start (h=5, Ah=0.1), then h to 0, then Ah to 0. Line B: h=0.1 fixed, Ah from 3 to 0, then h to 0.
    Hand-chosen to start deep in a phase and cross exactly one critical line during the sweep. They set the geometry of the protocol but do not enter the defect-density exponent.
  • Inverse quench rate range tau = Line A: tau in [1,100]. Line B: tau in [50,200].
    The narrow tau window for line B (a factor of 4) limits how strongly the exponent 1/3 can be distinguished from 1/2 in the data.
assumptions (4)
  • domain assumption Ising duality maps the quasiperiodic-J Hamiltonian (Eq. 1) to the quasiperiodic-h Hamiltonian (Eq. 2) and preserves the universality classes of critical lines A and B
    Invoked in Section III, first paragraph: the paper moves all numerics to the h-Ah plane and states that the phases swap while the dynamical nature of bulk single-particle excitations is unaltered. This is asserted, not demonstrated.
  • domain assumption Equilibrium critical exponents nu=1, z=1 for line A and nu=1, z=2 for line B from Ref. 18 are correct
    The KZ exponents 1/2 and 1/3 are computed from these inputs via nu d/(nu z+1); the paper does not re-derive them from first principles.
  • standard math The initial state at t=0 is the exact ground state of the starting Hamiltonian
    Standard KZ setup, stated at the start of Section III, with no justification for preparing this state numerically.
  • domain assumption The model is integrable via Jordan-Wigner fermionization, so the time evolution of up to 1600 spins is computed in single-particle space
    The Hamiltonian is quadratic, but the numerical method is never described in the paper. The feasibility of the 1600-site simulation presupposes free-fermion reducibility and an efficient single-particle evolution scheme.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Adiabatic dynamics of quasiperiodic transverse Ising model." pith.science (2026). https://pith.science/paper/3SDG34J7

@misc{pith2026190801959,
  author       = {Pith},
  title        = {Pith review of: Adiabatic dynamics of quasiperiodic transverse Ising model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3SDG34J7}},
  note         = {Machine review of arXiv:1908.01959}
}
read the original abstract

We study the non-equilibrium dynamics due to slowly taking a quasiperiodic Hamiltonian across its quantum critical point. The special quasiperiodic Hamiltonian that we study here has two different types of critical lines belonging to two different universality classes, one of them being the well known quantum Ising universality class. In this paper, we verify the Kibble Zurek scaling which predicts a power law scaling of the density of defects generated as a function of the rate of variation of the Hamiltonian. The exponent of this power law is related to the equilibrium critical exponents associated with the critical point crossed. We show that the power-law behavior is indeed obeyed when the two types of critical lines are crossed, with the exponents that are correctly predicted by Kibble Zurek scaling.

Figures

Figures reproduced from arXiv: 1908.01959 by the authors.

Figure 1
Figure 1. FIG. 1: Phase diagram of QPTIM consisting of FM, PM [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Variation of defect density with [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Variation of defect density with [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 16 canonical work pages

  1. [1]

    Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, UK, 1999)

    S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, UK, 1999). chakrabarti96 B. K. Chakrabarti, A. Dutta and P. Sen, Quantum Ising Phases and transitions in transverse Ising Models , m41 (Springer, heidelberg, 1996)

    work page 1999

  2. [2]

    Suzuki, J-i Inoue and Bikas K

    S. Suzuki, J-i Inoue and Bikas K. Chkarabarti, Quantum Ising Phases and Transitions in Transverse Ising Models (Springer, Lecture Notes in Physics, Vol. 862 (2013))

    work page 2013

  3. [3]

    Polkovnikov, K

    A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, Rev. Mod. Phys. 83 , 863 (2011)

    2011

  4. [4]

    Dziarmaga, Advances in Physics 59 , 1063 (2010)

    J. Dziarmaga, Advances in Physics 59 , 1063 (2010). dutta15 A. Dutta, G. Aeppli, B. K. Chakrabarti, U. Divakaran, T. Rosenbaum and D. Sen, Quantum Phase Transitions in Transverse Field Spin Models: From Statistical Physics to Quantum Information (Cambridge University Press, Cambridge, 2015). bloch08 I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 8...

    work page 2010

  5. [5]

    Mukherjee, U

    V. Mukherjee, U. Divakaran, A. Dutta, and D. Sen, Phys. Rev. B 76 , 174303 (2007); U. Divakaran, V. Mukherjee, A. Dutta, and D. Sen, J. Stat. Mech. P02007 (2009);U. Divakaran, A. Dutta, and D. Sen, Phys. Rev. B 78 , 144301 (2008); V. Mukherjee and A. Dutta, EPL 92 , 37004 (2010). mondal09 S. Mondal, K. Sengupta and D. Sen, Phys. Rev. B 79 , 045128 (2009)

    work page 2007

  6. [6]

    Aubry and G

    S. Aubry and G. Andr\'e, Ann. Israel Phys. Soc. 3 133 (1980)

    work page 1980

  7. [7]

    P. G. Harper, Proc. Phys. Soc. A 68, 874 (1955). anderson58 P. W. Anderson, Phys. Rev. 109 , 1492 (1958)

    work page 1955

  8. [8]

    Modugno, New J

    M. Modugno, New J. Phys. 11 033023 (2009)

    work page 2009

Show all 17 references
  1. [9]

    Deissler, E

    B. Deissler, E. Lucioni, M. Modugno, G. Roati, L. Tanzi, M. Zaccanti, M. Inguscio and G. Modugno, New J. Phys. 13, 023020 (2011). qp_lattice M. Schreiber, S. S. Hodgman, P. Bordia, H. P. L\"uschen, M. H. Fischer, R. Vosk, E. Altman, U. Schneider, and I. Bloch, Science 349 , 84...

    2011

  2. [10]

    Roati, C

    G. Roati, C. D Errico, L. Fallani, M. Fattori, C Fort, M

  3. [11]

    D. A. Huse, R. Nandkishore, V. Oganesyan, A. Pal and S. L

  4. [12]

    Igl\'oi, Zs

    F. Igl\'oi, Zs. Szatm\'ari, and Y.-C. Lin, Phys. Rev. B

  5. [13]

    Igl\'oi, G

    F. Igl\'oi, G. Ro\'osz and Y.-C. Lin, New J. Phys

  6. [14]

    P J. D. Crowley, A. Chandran and C. R. Laumann, Phys. Rev. Lett. 120 , 175702 (2018)

    2018

  7. [15]

    Dziarmaga, Phys

    J. Dziarmaga, Phys. Rev. Lett. 95 , 245701 (2005)

    2005

  8. [16]

    Caneva, R

    T. Caneva, R. Fazio, and G. E. Santoro Phys. Rev. B 76, 144427 (2007)

    2007

  9. [17]

    Sachdev, Quantum Phase Transitions (Cambridge

    S. Sachdev, Quantum Phase Transitions (Cambridge

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.