Kink-kink correlations in nonlinear quenches across a quantum critical point
Pith reviewed 2026-05-08 16:08 UTC · model grok-4.3
The pith
For algebraic quenches in the transverse-field Ising model, kink-kink correlations depend only on the Kibble-Zurek length when the quench is superlinear; otherwise they require an additional dephasing length and decay as a compressed-expon-
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using adiabatic perturbation theory, analytical arguments and numerical integration, the authors establish that the kink-kink correlation function after an algebraic quench depends only on the KZ length for superlinear protocols. For linear and sublinear protocols an additional dephasing length is required. The dephased correlator decays as exp[−(r/ξ)^α] where the exponent α varies continuously with the quench exponent.
What carries the argument
Adiabatic perturbation theory applied to the Bogoliubov-de Gennes equations of the TFIM, which yields the mode occupation numbers and the relative phases between low-energy modes that determine the kink-kink correlator.
If this is right
- The kink-kink correlator becomes a universal function of the scaled distance for all superlinear algebraic quenches.
- An extra dephasing length appears and controls the correlator when the quench exponent is less than or equal to one.
- The decay of the dephased correlator is a compressed exponential whose stretching exponent interpolates continuously between different values as the quench exponent is varied.
- The KZ mechanism must be supplemented by phase information to capture two-point functions in general nonlinear quenches.
Where Pith is reading between the lines
- Analogous dephasing scales and compressed-exponential forms may appear in other integrable critical systems with linear dispersion.
- Ultracold-atom or ion-trap experiments could measure the continuous change of the decay exponent by varying the quench protocol.
- The results suggest that coherence between distant modes survives the quench and shapes spatial correlations beyond simple defect counting.
Load-bearing premise
Adiabatic perturbation theory accurately captures the non-adiabatic transitions and the phases of the low-lying modes produced by the algebraic quench across the critical point.
What would settle it
A numerical computation of the exact time evolution for the TFIM under an algebraic quench that reveals either dependence on a length scale other than KZ and dephasing, or a decay form that is not compressed exponential with quench-dependent exponent.
Figures
read the original abstract
When a quantum system exhibiting a second order phase transition is quenched across the critical point in large but finite time, the dynamics are not adiabatic in the critical region and the Kibble-Zurek (KZ) mechanism provides a framework to determine local observables such as the mean defect density. However, to find higher-point functions, one has to go beyond the KZ paradigm asshown in recent works on one-dimensional transverse field Ising model (TFIM) following a linear quench. It has been found that (i) besides the KZ scale, the quench dynamics depend on another length scale that arises due to the finite phase difference between the low energy modes, and (ii) contrary to the expectations based on the KZ mechanism, in general, the correlation functions do not decay exponentially with distance. Motivated by these results for the linear quench, we are interested in understanding if these properties are universal, and consider the 1D TFIM when the transverse field varies algebraically in the vicinity of the critical field. We focus on the equal-time,longitudinal kink-kink correlation function at the end of the quench from the paramagnetic to the ferromagnetic phase, and find that (i) the correlator depends only on the KZ length for superlinear quenches, otherwise an additional dephasing length is required to describe it, and (ii) the dephased correlator decays as a compressed exponential with an exponent that changes continuously with the quench exponent. Our results are obtained using an adiabatic perturbation theory, analytical arguments and exact numerical integration of the relevant equations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the equal-time longitudinal kink-kink correlation function in the 1D transverse-field Ising model after algebraic (nonlinear) quenches of the transverse field across the quantum critical point from paramagnetic to ferromagnetic phase. Employing adiabatic perturbation theory to obtain mode occupations and phases, together with scaling arguments and direct numerical integration of the time-dependent Bogoliubov-de Gennes equations, the authors conclude that the correlator depends only on the Kibble-Zurek length for superlinear quenches; for sublinear quenches an additional dephasing length is required, and the dephased correlator decays as a compressed exponential whose exponent varies continuously with the quench exponent.
Significance. If the central claims hold, the work usefully generalizes prior linear-quench results to the algebraic case, demonstrating that the necessity of a dephasing scale and the compressed-exponential form are robust features of the post-quench correlations. The combination of analytic perturbation theory with exact numerical solution of the mode equations supplies an internal consistency check that is stronger than typical KZ studies relying solely on scaling or mean-field approximations.
minor comments (2)
- The abstract contains the typographical error 'asshown' (should be 'as shown').
- Section 3 (or wherever the dephasing length is introduced) would benefit from an explicit equation defining the dephasing length in terms of the mode phases before it is used in the scaling arguments.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our manuscript and the recommendation to accept. The referee's summary correctly captures our central results on the longitudinal kink-kink correlator after algebraic quenches in the 1D TFIM, including the distinction between superlinear and sublinear cases and the form of the dephased decay.
Circularity Check
No significant circularity; derivation self-contained via standard TFIM methods
full rationale
The paper's central results on kink-kink correlations for algebraic quenches follow from adiabatic perturbation theory applied to the TFIM mode equations, supplemented by scaling arguments and direct numerical solution of the time-dependent Bogoliubov-de Gennes equations. These steps are independent of the claimed functional forms and scaling behaviors; the numerics provide an external consistency check rather than a fit or redefinition. The reference to prior linear-quench results is motivational and not load-bearing for the nonlinear case. No self-definitional reductions, fitted inputs renamed as predictions, or uniqueness theorems imported from self-citations appear in the derivation chain.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math The 1D TFIM is exactly solvable by Jordan-Wigner transformation allowing numerical integration of time-dependent modes
- domain assumption Adiabatic perturbation theory captures leading non-adiabatic corrections near the critical point
Reference graph
Works this paper leans on
-
[1]
and emerging quantum technologies [ 2–4]. As an example, consider adiabatic quantum computing [ 3] in which starting from an easy-to-prepare ground state of a quantum system, the Hamiltonian is varied slowly so as to arrive at the ground state of a Hamiltonian that encodes the solution of a complex problem. However, if a phase transition occurs during suc...
work page 2000
-
[2]
Similarly, the second term in Eq
= eiΘ1cτ ω + 1 cτ ω (ω + 1)2 sin( π 2ω )F ω (B.9) Assuming that the contour can be deformed such that it passes through the saddle points where the integral C1 does not diverge with R, the function h1 is evaluated at those saddle points where cos Θ1 < 0. Similarly, the second term in Eq. ( B.1) gives b2 = 2 sin( π 2ω )G which yields h2(x∗
-
[3]
= eiΘ1cτ ω + 1 cτ ω (ω + 1)2 sin( π 2ω )G ω (B.10) For the third term in Eq. ( B.1), b3 = −ie iπ 2ω (F + G) so that the saddles occur at x∗ 3 ∝ ei(2m+1)π(ω+1)e− iπ(ω+1) 2ω and Eq. ( B.7) then gives h3(x∗
-
[4]
Finally, for the fourth term in Eq
= eiΘ3cτ ω + 1 cτ ω (ω + 1)(F + G) ω (B.11) where Θ3 = (2m + 1)πω with non-negative integer m chosen such that cos Θ3 < 0. Finally, for the fourth term in Eq. ( B.1), b4 = ie− iπ 2ω (F + G), the saddles occur at x∗ 4 ∝ e iπ(ω+1) 2ω , and Eq. ( B.7) yields h4(x∗
-
[5]
( B.1) by its modulus, and the rest of the discussion above remains unchanged
= − cτ ω + 1 cτ ω (ω + 1)(F + G) ω (B.12) When sin( π 2ω ) < 0, the excitation probability is obtained on replacing sin ( π 2ω ) in Eq. ( B.1) by its modulus, and the rest of the discussion above remains unchanged. For sin ( π 2ω ) = 0 , pk ≈ π2 36 [2 + eiX(−1)n(F +G) + h.c.] for ω = (2 n)−1, n = 1 , 2, 3.... The first term in pk gives R π 0 dk 2 cos(kr) ...
work page 2000
-
[6]
Sachdev, Cambridge University Press (2011)
S. Sachdev, Cambridge University Press (2011)
work page 2011
-
[7]
I. M. Georgescu, S. Ashhab, and F. Nori, Rev. Mod. Phys. 86, 153 (2014)
work page 2014
- [8]
-
[9]
M. O. Scully and W. G. Unruh, Science 390, 998 (2025)
work page 2025
-
[10]
T. W. B. Kibble, J. Phys. A 9, 1387 (1976)
work page 1976
-
[11]
W. H. Zurek, Nature 317, 505 (1985)
work page 1985
-
[12]
Damski, Physical review letters 95, 035701 (2005)
B. Damski, Physical review letters 95, 035701 (2005)
work page 2005
-
[13]
B. Damski and W. H. Zurek, Physical Review A—Atomic, Molecular, and Optical Physics 73, 063405 (2006)
work page 2006
-
[14]
D. Sadhukhan, A. Sinha, A. Francuz, J. Stefaniak, M. M. Rams, J. Dziarmaga, and W. H. Zurek, Phys. Rev. B 101, 144429 (2020)
work page 2020
-
[15]
F. Suzuki and W. H. Zurek, Proceedings of the National Academy of Sciences 122, e2523903122 (2025)
work page 2025
- [16]
- [17]
-
[18]
H. Ricateau, L. F. Cugliandolo, and M. Picco, J. Stat. Mech. , P013201 (2018)
work page 2018
- [19]
- [20]
-
[21]
Priyanka, S. Chatterjee, and K. Jain, Journal of Statistical Mechanics: Theory and Experi- ment 2021, 033208 (2021)
work page 2021
- [22]
- [23]
-
[24]
C. Tang, K. Sfairopoulos, W. Zhang, and C. Ding, Physical Review E 112, 024111 (2025)
work page 2025
- [25]
-
[26]
D. Sen, K. Sengupta, and S. Mondal, Physical Review Letters 101, 016806 (2008)
work page 2008
-
[27]
M. Kolodrubetz, B. K. Clark, and D. A. Huse, Phys. Rev. Lett. 109, 015701 (2012)
work page 2012
- [28]
- [29]
- [30]
-
[31]
M. Schmitt, M. M. Rams, J. Dziarmaga, M. Heyl, and W. H. Zurek, Science Advances 8, eabl6850 (2022)
work page 2022
- [32]
- [33]
- [34]
-
[35]
J. Soto-Garcia and N. Chepiga, Physical Review B 113, 085430 (2026)
work page 2026
-
[36]
S. Maegochi, K. Ienaga, and S. Okuma, Physical Review Letters 129, 227001 (2022)
work page 2022
-
[37]
K. Du, X. Fang, C. Won, C. De, F.-T. Huang, W. Xu, H. You, F. J. Gómez-Ruiz, A. del Campo, and S.-W. Cheong, Nature Physics , 1 (2023)
work page 2023
-
[38]
H. Wang, X. Li, and C. Li, Nature Communications 16, 10584 (2025)
work page 2025
-
[39]
B. Ko, J. W. Park, and Y.-i. Shin, Nature physics 15, 1227 (2019)
work page 2019
-
[40]
K. Lee, S. Kim, T. Kim, and Y.-i. Shin, Nature Physics 20, 1570 (2024)
work page 2024
- [41]
-
[42]
Z.-J. Cheng, Y.-K. Wu, S.-J. Li, and et al., Science Advances 10, eadr9527 (2024) , https://science.org
work page 2024
-
[43]
S. Deutschländer, P. Dillmann, G. Maret, and P. Keim, Proc. Natl. Acad. Sci. 112, 20 , 6925 (2015)
work page 2015
- [44]
-
[45]
L. Wang, C. Zhou, T. Tu, H.-W. Jiang, G.-P. Guo, and G.-C. Guo, Phys. Rev. A 89, 022337 (2014)
work page 2014
- [46]
-
[47]
A. D. King, S. Suzuki, J. Raymond, and et al., Nature Physics 18, 1324 (2022)
work page 2022
-
[48]
A. Miessen, D. J. Egger, I. Tavernelli, and G. Mazzola, PRX Quantum 5, 040320 (2024)
work page 2024
- [49]
-
[50]
A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, Rev. Mod. Phys. 83, 863 (2011)
work page 2011
- [51]
-
[52]
A. del Campo and W. H. Zurek, Int. J. Mod. Phys. A 29, 1430018 (2014)
work page 2014
- [53]
- [54]
-
[55]
R. Barankov and A. Polkovnikov, Physical review letters 101, 076801 (2008)
work page 2008
-
[56]
R. W. Cherng and L. S. Levitov, Phys. Rev. A 73, 043614 (2006)
work page 2006
- [57]
-
[58]
K. Roychowdhury, R. Moessner, and A. Das, Physical Review B 104, 014406 (2021)
work page 2021
-
[59]
R. J. Nowak and J. Dziarmaga, Physical Review B 104, 075448 (2021)
work page 2021
- [60]
- [61]
-
[62]
C. De Grandi and A. Polkovnikov, Quantum Quenching, Annealing and Computation , 75 (2010)
work page 2010
- [63]
-
[64]
G. B. Mbeng, A. Russomanno, and G. E. Santoro, SciPost Physics Lecture Notes , 082 (2024)
work page 2024
- [65]
-
[66]
P. Calabrese, F. H. L. Essler, and M. Fagotti, Physical Review Letters 106, 227203 (2011)
work page 2011
-
[67]
P. Calabrese, F. H. L. Essler, and M. Fagotti, Journal of Statistical Mechanics: Theory and Experiment 2012, P07016 (2012)
work page 2012
-
[68]
O. V. Ivakhnenko, S. N. Shevchenko, and F. Nori, Physics Reports 995, 1 (2023)
work page 2023
-
[69]
C. M. Bender and S. A. Orszag, McGraw-Hill Book Company,New York, NY (1978)
work page 1978
- [70]
-
[71]
L. Vidmar and M. Rigol, Journal of Statistical Mechanics: Theory and Experiment 2016, 064007 (2016)
work page 2016
-
[72]
M. Heyl, A. Polkovnikov, and S. Kehrein, Phys. Rev. Lett. 110, 135704 (2013)
work page 2013
-
[73]
Heyl, Reports on Progress in Physics 81, 054001 (2018)
M. Heyl, Reports on Progress in Physics 81, 054001 (2018) . 31
work page 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.