Integrability breaking of the $E_8$ spectrum in the transverse-field Ising chain: heavy-particle lifetimes, lineshapes, and robustness of high-energy $E_8$ peaks in BaCo$_2$V$_2$O$_8$

Near the transverse-field-induced one-dimensional quantum critical point of the quasi-one-dimensional Ising antiferromagnet BaCo$_2$V$_2$O$_8$ ($H_c^{1D}\approx 4.7$ T), inelastic neutron scattering observes the particle spectrum of Zamolodchikov's integrable $E_8$ field theory. In that integrable theory only the three lightest particles lie below the lightest two-particle threshold $2m_1$; the heavier ones ($m_4\approx 2.405\,m_1$, $m_5\approx 2.956\,m_1$, ...) are stable only because of integrability. Yet in the material — which is not integrable — the high-energy single-particle peaks remain sharp instead of decaying into lower-mass multi-particle modes. The open problem, posed by the discovering paper and unresolved in the subsequent literature, is a quantitative theory of integrability-breaking effects on this spectrum: decay widths, lifetimes, lineshapes, and spectral-weight transfer of the heavy $E_8$ particles as field detuning, transverse anisotropy, and interchain coupling are increased, together with the experimental resolution of the higher $E_8$ modes and of the integrable-to-non-integrable crossover, and the response of the $E_8$ modes to time-dependent (quench, drive, transport) perturbations.

Background

A spin-1/2 Ising-like chain in a magnetic field transverse to its easy axis is the canonical model of a one-dimensional quantum phase transition: at a critical transverse field the excitation gap closes at a quantum critical point (QCP) whose long-wavelength theory is the Ising conformal field theory (CFT). Perturbing the critical Ising CFT by a magnetic (spin-density) perturbation — in the chain, an effective longitudinal field along the easy axis — yields Zamolodchikov's celebrated integrable quantum field theory, whose particle content is organized by the exceptional Lie algebra $E_8$: eight stable massive particles with exact mass ratios, beginning with $m_2/m_1=2\cos(\pi/5)\approx 1.618$ (the golden ratio) and $m_3/m_1=2\cos(\pi/30)\approx 1.989$. 'Integrability' means the theory has infinitely many conserved quantities, which make scattering exactly elastic and factorized and forbid particle decay. A key kinematic fact is that only the three lightest $E_8$ particles lie below the lightest two-particle threshold $2m_1$, the energy where the continuum of two-particle states begins; the heavier particles ($m_4\approx 2.405\,m_1$, $m_5\approx 2.956\,m_1$, up to $m_8$) are stable only because of integrability — in a generic non-integrable theory they would be resonances with finite decay widths. The standard tool for weak integrability breaking is form-factor perturbation theory (FFPT), perturbation theory built on the exactly known matrix elements ('form factors') of the integrable model; it was applied to the Ising field theory with both thermal (energy-density) and magnetic perturbations by Delfino, Mussardo, and Simonetti, and used to study the decay of above-threshold particles in that pure $1+1$-dimensional field theory.

The quasi-one-dimensional antiferromagnets CoNb$_2$O$_6$ and BaCo$_2$V$_2$O$_8$ (BCVO) realize this physics experimentally. In BCVO the Co$^{2+}$ screw chains are described at low energies by a four-site-periodic Heisenberg–Ising Hamiltonian in a transverse field. Nuclear magnetic resonance (NMR) measurements of the spin-lattice relaxation rate locate a 1D Ising QCP at $H_c^{1D}\approx 4.7$ T, hidden inside the field range where weak interchain coupling still produces three-dimensional magnetic order; in that ordered phase the interchain coupling acts on each chain as an effective staggered longitudinal field, which is precisely the magnetic perturbation generating the $E_8$ spectrum. Inelastic neutron scattering (INS) at $H_c^{1D}$ indeed finds the zone-center excitation energies matching the $E_8$ pattern: the transverse dynamical structure factor (DSF — the momentum- and frequency-resolved spin correlation function probed by INS) shows single-particle peaks through $m_6$ and multi-particle continua above the two-particle threshold, while terahertz (THz) spectroscopy had resolved up to $m_5$. Because the four Co$^{2+}$ ions per screw-chain repeat are shifted along the chain axis, the measured spectrum also contains 'zone-folded' shadow branches at shifted momenta, which complicate the identification of high-energy modes.

The puzzle is that the real material is not integrable — detuning the transverse field from $H_c^{1D}$, the transverse-exchange and staggered-transverse-field terms of the four-periodic chain, and finite interchain coupling all break integrability — yet the high-energy single-particle peaks above $2m_1$ remain sharp and do not decay into lower-mass multi-particle states, and the excitations stay non-diffusive over a finite field range away from the QCP. The discovering paper (Zou et al., Phys. Rev. Lett. 127, 077201 (2021)) states that the underlying physics 'needs to be further addressed theoretically' and identifies the study of such non-integrable effects, with BCVO as a test bed in which the effective model becomes a perturbed quantum $E_8$ model whose heavy particles begin to decay, as an important open direction, together with better experimental identification of the excitations and the non-equilibrium dynamics of the modes.

Problem Statement

The problem, as posed by Zou et al. (Phys. Rev. Lett. 127, 077201 (2021)) and still open, is to give a quantitative account of integrability-breaking effects on the E8 spectrum of the near-critical transverse-field Ising chain as realized in BaCo2V2O8, under one overarching objective: to understand, predict, and spectroscopically verify how the exactly solved E8 quasiparticle spectrum degrades as the integrability-breaking perturbations present in the material are turned on and increased. Its linked components are:

(1) Stability mechanism and decay theory. Explain quantitatively why the high-energy single-particle E8 peaks observed by INS at the 1D QCP ($H_c^{1D}\approx 4.7$ T) — in particular the above-threshold particles $m_4\approx 2.405\,m_1$, $m_5\approx 2.956\,m_1$, and $m_6$, noting that only $m_1$, $m_2$, $m_3$ lie below the lightest two-particle threshold $2m_1$ — remain sharp and do not rapidly decay into lower-mass multi-particle modes; and compute, as functions of the integrability-breaking couplings (transverse-field detuning from the QCP, the transverse-exchange and staggered-transverse-field terms of the four-periodic effective chain, and finite interchain coupling beyond its mean-field staggered-field action), the decay widths (lifetimes), lineshapes, and spectral-weight transfer of the heavy E8 particles, including where the single-particle peaks dissolve into the multi-particle continuum.

(2) Spectroscopic delineation. Provide the theoretical predictions needed to interpret and extend the identification of the E8 excitations in BCVO beyond the modes resolved to date — the single-particle peaks already reported through $m_6$ (and through $m_5$ by THz) — namely the mode energies and spectral weights of the still-unresolved higher modes $m_7$–$m_8$ together with the intrinsic linewidths of the above-threshold modes, at the level of the roughly 0.1 meV instrumental resolution and accounting for zone folding of the four-fold screw chain, together with the predicted crossover from the integrable E8 regime to non-integrable dynamics as the transverse field and the other perturbations are varied.

(3) Out-of-equilibrium response. Characterize the dynamics of the E8 modes under time-dependent perturbations (quench, periodic drive, or transport-relevant protocols) within the same effective model, clarifying how integrability-related features of the E8 spectrum manifest away from equilibrium.

The objective is completed when component (1) — the mandatory theoretical core — yields the predicted energies, widths, lineshapes, and spectral weights of the E8-related excitations as functions of the perturbation strengths, and these are quantitatively consistent with the spectroscopic data reported for BCVO (zone-center INS, Zou et al. 2021; THz, Zhang et al. 2020): the predicted peak positions and relative spectral weights must match those data within the uncertainties quoted there, and the predicted intrinsic linewidths of the above-threshold single-particle peaks must be quoted as definite numerical values together with their explicit dependence on the perturbation strengths, including in the regime where they lie below the roughly 0.1 meV instrumental resolution at which those peaks presently appear resolution-limited, thereby resolving the robustness puzzle. Agreement with the companion Ising-chain material CoNb$_2$O$_6$ (THz, Amelin et al. 2020) is an optional cross-check rather than part of the mandatory core. Completing component (1) suffices to complete the objective; components (2) and (3) are independent extensions that build on it: component (2) is completed to the extent that the theory supplies the predicted mode energies, intrinsic linewidths, and crossover behavior needed to interpret the existing data and to extend the spectroscopic identification towards the higher modes $m_7$–$m_8$, and component (3) is completed when the E8 response to at least one specified time-dependent protocol (a quench, a periodic drive, or a transport-relevant protocol) is computed within the same framework as an explicit mode-resolved spectral function — a time-dependent or frequency-domain response resolved by E8 particle mode. All three address the same physical question of how the integrable E8 structure survives or breaks in a real material.

References

1. H. Zou, Y. Cui, X. Wang, Z. Zhang, J. Yang, G. Xu, A. Okutani, M. Hagiwara, M. Matsuda, G. Wang, G. Mussardo, K. Hódsági, M. Kormos, Z. Z. He, S. Kimura, R. Yu, W. Yu, J. Ma, and J. Wu, 'E8 Spectra of Quasi-one-dimensional Antiferromagnet BaCo2V2O8 under Transverse Field', Phys. Rev. Lett. 127, 077201 (2021). DOI: 10.1103/PhysRevLett.127.077201; arXiv:2005.13302; https://arxiv.org/abs/2005.13302
2. A. B. Zamolodchikov, 'Integrals of motion and S-matrix of the (scaled) T = Tc Ising model with magnetic field', Int. J. Mod. Phys. A 4, 4235 (1989). DOI: 10.1142/S0217751X8900176X; https://doi.org/10.1142/S0217751X8900176X
3. R. Coldea, D. A. Tennant, E. M. Wheeler, E. Wawrzynska, D. Prabhakaran, M. Telling, K. Habicht, P. Smeibidl, and K. Kiefer, 'Quantum Criticality in an Ising Chain: Experimental Evidence for Emergent E8 Symmetry', Science 327, 177 (2010). DOI: 10.1126/science.1180085; https://doi.org/10.1126/science.1180085
4. G. Delfino, G. Mussardo, and P. Simonetti, 'Non-integrable quantum field theories as perturbations of certain integrable models', Nucl. Phys. B 473, 469 (1996). DOI: 10.1016/0550-3213(96)00265-9; arXiv:hep-th/9603011; https://arxiv.org/abs/hep-th/9603011
5. G. Delfino, P. Grinza, and G. Mussardo, 'Decay of particles above threshold in the Ising field theory with magnetic field', Nucl. Phys. B 737, 291 (2006). DOI: 10.1016/j.nuclphysb.2005.12.024; arXiv:hep-th/0507133; https://arxiv.org/abs/hep-th/0507133
6. Z. Zhang, K. Amelin, X. Wang, H. Zou, J. Yang, U. Nagel, T. Rõõm, T. Dey, A. A. Nugroho, T. Lorenz, J. Wu, and Z. Wang, 'Observation of E8 Particles in an Ising Chain Antiferromagnet', Phys. Rev. B 101, 220411 (2020). DOI: 10.1103/PhysRevB.101.220411; arXiv:2005.13772; https://arxiv.org/abs/2005.13772
7. K. Amelin, J. Engelmayer, J. Viirok, U. Nagel, T. Rõõm, T. Lorenz, and Z. Wang, 'Experimental Observation of Quantum Many-Body Excitations of E8 Symmetry in the Ising Chain Ferromagnet CoNb2O6', Phys. Rev. B 102, 104431 (2020). DOI: 10.1103/PhysRevB.102.104431; arXiv:2006.12956; https://arxiv.org/abs/2006.12956
8. X. Wang, H. Zou, K. Hódsági, M. Kormos, G. Takács, and J. Wu, 'Cascade of singularities in the spin dynamics of a perturbed quantum critical Ising chain', Phys. Rev. B 103, 235117 (2021). DOI: 10.1103/PhysRevB.103.235117; arXiv:2103.09128; https://arxiv.org/abs/2103.09128
9. X. Wang, K. Puzniak, K. Schmalzl, C. Balz, M. Matsuda, A. Okutani, M. Hagiwara, J. Ma, J. Wu, and B. Lake, 'Spin dynamics of the E8 particles', Science Bulletin 69, 2974 (2024). DOI: 10.1016/j.scib.2024.07.040; arXiv:2308.00249; https://arxiv.org/abs/2308.00249
10. M. Lencsés, G. Mussardo, and G. Takács, 'Quantum integrability vs experiments: correlation functions and dynamical structure factors', J. Phys. A: Math. Theor. 56, 383001 (2023). DOI: 10.1088/1751-8121/acf255; arXiv:2303.16556; https://arxiv.org/abs/2303.16556
11. R. G. Jha, A. Milsted, D. Neuenfeld, J. Preskill, and P. Vieira, 'Real-Time Scattering in Ising Field Theory using Matrix Product States', Phys. Rev. Research 7, 023266 (2025). DOI: 10.1103/9dxz-k5wb; arXiv:2411.13645; https://arxiv.org/abs/2411.13645
12. Y. Gao, Y. Jiang, and J. Wu, 'Mesons in a quantum Ising ladder', J. High Energy Phys. 2025, 72 (2025). DOI: 10.1007/JHEP07(2025)072; arXiv:2502.15463; https://arxiv.org/abs/2502.15463
13. X. Wang and J. Wu, 'Universal meson spectra near (1+1)-dimensional Ising criticality', arXiv:2608.20221 (2026); https://arxiv.org/abs/2608.20221
14. J. Yang, X. Wang, and J. Wu, 'Magnetic excitations in the one-dimensional Heisenberg–Ising model with external fields and their experimental realizations', J. Phys. A: Math. Theor. 56, 013001 (2023). DOI: 10.1088/1751-8121/acad48; arXiv:2301.05823; https://arxiv.org/abs/2301.05823
15. K. Amelin, J. Viirok, U. Nagel, T. Rõõm, J. Engelmayer, T. Dey, A. A. Nugroho, T. Lorenz, and Z. Wang, 'Quantum spin dynamics of quasi-one-dimensional Heisenberg-Ising magnets in a transverse field: confined spinons, E8 spectrum, and quantum phase transitions', J. Phys. A: Math. Theor. 55, 484005 (2022). DOI: 10.1088/1751-8121/aca6b8; https://doi.org/10.1088/1751-8121/aca6b8
16. B. Gabai and X. Yin, 'On The S-Matrix of Ising Field Theory in Two Dimensions', arXiv:1905.00710 (2019); https://arxiv.org/abs/1905.00710
17. D. Schuricht and F. H. L. Essler, 'Dynamics in the Ising field theory after a quantum quench', J. Stat. Mech. (2012) P04017. DOI: 10.1088/1742-5468/2012/04/p04017; arXiv:1203.5080; https://arxiv.org/abs/1203.5080
18. N. J. Robinson, I. Pérez Castillo, and E. Guzmán-González, 'Quantum quench in a driven Ising chain', Phys. Rev. B 103, L140407 (2021). DOI: 10.1103/PhysRevB.103.L140407; arXiv:2011.14345; https://arxiv.org/abs/2011.14345
19. N. Albert, Y. Zhang, and H.-H. Tu, 'Truncated Gaussian basis approach for simulating many-body dynamics', arXiv:2410.04204 (2024); https://arxiv.org/abs/2410.04204
20. H. A. Labib, G. C. Toga, J. K. Freericks, and A. F. Kemper, 'Probing bound state relaxation dynamics in systems out-of-equilibrium on quantum computers', APL Computational Physics 2, 016102 (2026). DOI: 10.1063/5.0301170; arXiv:2507.22988; https://arxiv.org/abs/2507.22988
21. X. He, X. Wang, and J. Wu, 'Numerical renormalization group integrated Hamiltonian truncation: Toward generic deformation of integrable lattice models', arXiv:2511.13218 (2025); https://arxiv.org/abs/2511.13218
22. J. Taylor, M. Yusuf, and Z. Papić, 'Fuzzy Spectroscopy of Bound States in Massive Quantum Field Theories', arXiv:2608.07655 (2026); https://arxiv.org/abs/2608.07655
23. J. Vovrosh, J. de Hond, S. Julià-Farré, J. Knolle, and A. Dauphin, 'Meson spectroscopy of exotic symmetries of Ising criticality in Rydberg atom arrays', arXiv:2506.21299 (2025); https://arxiv.org/abs/2506.21299