In the Ising spin chain with transverse and longitudinal fields, the longitudinal field confines domain walls into mesonic bound states, and quench dynamics proceeds in stages: rapid relaxation to a prethermal Gibbs ensemble with approximately conserved meson number, followed by rare-event avalanches, and finally true thermalization driven by meson-number violation through a mechanism akin to the Schwinger effect of pair creation from the vacuum by a strong field, at times exponential in the inverse longitudinal field. This problem formulates the open theory question left by Birnkammer, Bastianello, and Knap (Nature Communications 13, 7663 (2022)): constructing a quantitative kinetic theory that incorporates microscopic meson creation–annihilation rates and describes the drift from the prethermal meson-conserving ensemble to the true thermal state, together with the extension of the prethermal description to direct meson–meson interactions beyond the dilute hard-rod approximation, in particular under power-law long-range couplings.

Background

Far-from-equilibrium quantum many-body systems often approach thermal equilibrium in stages rather than monotonically. Prethermalization is the phenomenon in which observables first plateau at values captured by a statistical ensemble reflecting an approximate conservation law, with full thermalization deferred to much later times. A clean theoretical setting exhibiting this is the ferromagnetic Ising chain in both a transverse and a longitudinal field, H^=−∑j[σ^j+1zσ^jz+h⊥σ^jx+h∥σ^jz],\hat H=-\sum_j\left[\hat\sigma^z_{j+1}\hat\sigma^z_j+h_\perp\hat\sigma^x_j+h_\parallel\hat\sigma^z_j\right], where the transverse field h⊥h_\perp makes domain walls mobile (at h∥=0h_\parallel=0 the model maps to free fermions by a Jordan–Wigner transformation, i.e. a rewriting of spins in terms of fermionic variables, and is ordered for ∣h⊥∣<1|h_\perp|<1), while the longitudinal field h∥h_\parallel explicitly breaks the Z2\mathbb{Z}_2 spin-flip symmetry and confines domain walls: separating two domain walls by a distance ℓ\ell costs the string energy 2h∥σˉ ℓ2h_\parallel\bar\sigma\,\ell, with σˉ=(1−h⊥2)1/8\bar\sigma=(1-h_\perp^2)^{1/8} the renormalized ground-state magnetization. The spectrum therefore reorganizes into mesons — bound states of domain walls held together by the linear confining potential, in close analogy to confinement in particle physics. Starting from one of the two symmetry-broken ground states at h∥=0h_\parallel=0 and suddenly switching on both field components (e.g. from (h⊥,h∥)=(0.4,0)(h_\perp,h_\parallel)=(0.4,0) to h⊥∈{0.65,0.75,0.85,0.95}h_\perp\in\{0.65,0.75,0.85,0.95\}, h∥=0.05h_\parallel=0.05 in the reference protocol), Birnkammer, Bastianello, and Knap established the following multi-stage picture, combining tensor-network simulations (variational representations of many-body quantum states as networks of small tensors, enabling simulations of chains far beyond exact diagonalization) with semiclassical analytics. First, local observables rapidly relax to a prethermal state described by a Gibbs ensemble of a gas of mesons with an exactly conserved meson number; its thermodynamics is that of a multi-species hard-rod gas (mesons treated as impenetrable rods of length ℓ(J,k)\ell(J,k), i.e. an excluded-volume approximation), with meson energies quantized semiclassically by the Bohr–Sommerfeld rule J=2π(n−1/2)J=2\pi(n-1/2), which assigns energy levels to classical periodic orbits of the confining motion. The hard-rod description is controlled only at low meson density; deviations are already visible at ρ≈0.4\rho\approx0.4. Second, meson number is only approximately conserved: meson pairs are slowly created out of the ordered background by a mechanism akin to the Schwinger effect — pair creation from the vacuum by a strong field, here the longitudinal field acting on the domain-wall 'charges'. Rare nucleation events trigger an avalanche of meson collisions that sets the mesons in motion and drives relaxation to the prethermal Gibbs ensemble on a timescale tPreTh∼(ρ2dmax⁡v)−1t_{\rm PreTh}\sim(\rho^2 d_{\max}v)^{-1} that is algebraic in the model parameters (with dmax⁡∼4h⊥/(h∥σˉ)d_{\max}\sim 4h_\perp/(h_\parallel\bar\sigma) the avalanche length and vv the meson velocity). Third, true thermalization, in which meson number is violated by Schwinger-type processes, requires a time tTht_{\rm Th} that is exponential in 1/h∥1/h_\parallel — a parametric separation exploited through the truncated Wigner approximation (TWA, a semiclassical phase-space method that samples quantum states as classical distributions and evolves them classically, accurate when quantum interference is negligible) applied to a projected dynamics that enforces meson-number conservation. Related earlier results set the stage: entanglement dynamics in confining spin chains was computed by Scopa, Calabrese, and Bastianello; numerical studies of the tilted (confined) Ising chain by Verdel, Zhu, and Heyl found a dynamical localization transition in string breaking, i.e. a dichotomy between fast and dynamically-suppressed string breaking depending on quench parameters; and Surace et al. mapped the confinement physics of Z2\mathbb{Z}_2 lattice gauge theories onto Rydberg-atom quantum simulators. What the reference work did not provide is a kinetic theory — a Boltzmann-type description with microscopic meson-number-changing rates — nor any treatment of meson–meson interactions beyond the dilute hard-rod approximation, which is the problem formulated here.

Problem Statement

Construct a quantitative theory of thermalization for one-dimensional quantum systems with confined excitations, beyond the dilute meson-number-conserving approximation, as left open by Birnkammer, Bastianello, and Knap. The single overarching objective is a controlled, quantitative description of how the prethermal meson-conserving Gibbs ensemble drifts to the true thermal state, comprising three jointly supporting components: (1) derive from the microscopic confined-chain Hamiltonian (the tilted Ising chain H^=−∑j[σ^j+1zσ^jz+h⊥σ^jx+h∥σ^jz]\hat H=-\sum_j[\hat\sigma^z_{j+1}\hat\sigma^z_j+h_\perp\hat\sigma^x_j+h_\parallel\hat\sigma^z_j] being the reference realization, and generic one-dimensional systems with linearly confined domain-wall excitations) the rates of meson-number-changing processes — Schwinger-type creation of meson pairs out of the ordered background and inelastic fusion or fission of mesons — going beyond the order-of-magnitude estimates of the source work; (2) build from these rates a kinetic theory (a Boltzmann-type description, or an alternative controlled framework) that quantitatively describes the drift from the prethermal plateau to thermal equilibrium: the time dependence of local observables such as the longitudinal magnetization and energy-density observables through the crossover, and the thermalization time tTht_{\rm Th} with its parametric dependence on h∥h_\parallel and meson density, including how the avalanche mechanism terminates in the number-violating regime; and (3) extend the prethermal description to direct meson–meson interactions beyond the hard-rod excluded-volume approximation, in particular for power-law long-range couplings, determining how such interactions modify the prethermal regime, the avalanche, and the timescale separations. The three components form one coherent objective — a quantitative theory of confined-chain thermalization beyond the number-conserving dilute approximation — and component (3) is the same extension needed wherever (2) is applied at non-dilute meson densities. The method is otherwise unconstrained: semiclassical treatments must come with controlled quantum corrections (the regime where Bohr–Sommerfeld quantization of the meson spectrum and the truncated Wigner approximation fail, i.e. when the fermion bandwidth is comparable to the longitudinal field, must be handled), and purely numerical demonstrations do not by themselves constitute the requested theory. A complete answer demonstrates quantitative validity against specified benchmarks such as the source quench protocol (from (h⊥,h∥)=(0.4,0)(h_\perp,h_\parallel)=(0.4,0) to h⊥=0.65h_\perp=0.65, h∥=0.05h_\parallel=0.05, and neighboring values) where exact-diagonalization or tensor-network data exist. The tilted Ising chain is the reference realization and a complete answer must handle it explicitly; the generic confined class and the power-law long-range couplings of component (3) must be covered by the same theory but may be argued analytically or qualitatively, since quantitative benchmarks exist only for the reference chain.

References

  1. Stefan Birnkammer, Alvise Bastianello, Michael Knap, 'Prethermalization in one-dimensional quantum many-body systems with confinement', Nature Communications 13, 7663 (2022). DOI: 10.1038/s41467-022-35301-6; arXiv:2202.12908. https://arxiv.org/abs/2202.12908
  2. Stefano Scopa, Pasquale Calabrese, Alvise Bastianello, 'Entanglement dynamics in confining spin chains', Physical Review B 105, 125413 (2022). DOI: 10.1103/PhysRevB.105.125413. https://doi.org/10.1103/PhysRevB.105.125413
  3. Federica M. Surace, Paolo P. Mazza, Giuliano Giudici, Alessio Lerose, Andrea Gambassi, Marcello Dalmonte, 'Lattice Gauge Theories and String Dynamics in Rydberg Atom Quantum Simulators', Physical Review X 10, 021041 (2020). DOI: 10.1103/PhysRevX.10.021041. https://doi.org/10.1103/PhysRevX.10.021041
  4. Jacob H. Robertson, Riccardo Senese, Fabian H. L. Essler, 'Decay of long-lived oscillations after quantum quenches in gapped interacting quantum systems', Physical Review A 109, 032208 (2024). DOI: 10.1103/PhysRevA.109.032208; arXiv:2307.04466. https://arxiv.org/abs/2307.04466
  5. Roberto Verdel, Guo-Yi Zhu, Markus Heyl, 'Dynamical Localization Transition of String Breaking in Quantum Spin Chains', Physical Review Letters 131, 230402 (2023). DOI: 10.1103/PhysRevLett.131.230402; arXiv:2304.12957. https://arxiv.org/abs/2304.12957
  6. Matteo Gallone, 'On Exponentially Long Prethermalization Timescales in Isolated Quantum Systems', arXiv:2604.13781 (2026). https://arxiv.org/abs/2604.13781
  7. Nishan Ranabhat, Mario Collura, 'Thermalization of long range Ising model in different dynamical regimes: A full counting statistics approach', SciPost Physics Core 7, 017 (2024). DOI: 10.21468/SciPostPhysCore.7.2.017; arXiv:2212.00533. https://arxiv.org/abs/2212.00533
  8. Nishan Ranabhat, Alessandro Santini, Emanuele Tirrito, Mario Collura, 'Dynamical deconfinement transition driven by density of excitations', Physical Review B 111, 024304 (2025). DOI: 10.1103/PhysRevB.111.024304; arXiv:2310.02320. https://arxiv.org/abs/2310.02320
  9. Leonardo Biagetti, Maciej Lebek, Milosz Panfil, Jacopo De Nardis, 'Generalised BBGKY hierarchy for near-integrable dynamics', arXiv:2408.00593 (2024). https://arxiv.org/abs/2408.00593
  10. Florian Kranzl, Stefan Birnkammer, Manoj K. Joshi, Alvise Bastianello, Rainer Blatt, Michael Knap, Christian F. Roos, 'Observation of Magnon Bound States in the Long-Range, Anisotropic Heisenberg Model', Physical Review X 13, 031017 (2023). DOI: 10.1103/PhysRevX.13.031017; arXiv:2212.03899. https://arxiv.org/abs/2212.03899
  11. Elizabeth R. Bennewitz, Brayden Ware, Alexander Schuckert, Alessio Lerose, Federica M. Surace, Ron Belyansky, William Morong, De Luo, Arinjoy De, Kate S. Collins, Or Katz, Christopher Monroe, Zohreh Davoudi, Alexey V. Gorshkov, 'Simulating Meson Scattering on Spin Quantum Simulators', Quantum 9, 1773 (2025). DOI: 10.22331/q-2025-06-17-1773; arXiv:2403.07061. https://arxiv.org/abs/2403.07061
  12. Joseph Vovrosh, Julius de Hond, Sergi Julià-Farré, Johannes Knolle, Alexandre Dauphin, 'Meson spectroscopy of exotic symmetries of Ising criticality in Rydberg atom arrays', arXiv:2506.21299 (2025). https://arxiv.org/abs/2506.21299
  13. Yaohua Li, Devendra Singh Bhakuni, Yong-Chun Liu, Marcello Dalmonte, 'Metastable confinement in Rydberg lattice gauge theories', arXiv:2602.22890 (2026). https://arxiv.org/abs/2602.22890
  14. Jaka Vodeb, Jean-Yves Desaules, Andrew Hallam, Andrea Rava, Gregor Humar, Dennis Willsch, Fengping Jin, Madita Willsch, Kristel Michielsen, Zlatko Papić, 'Stirring the false vacuum via interacting quantized bubbles on a 5,564-qubit quantum annealer', Nature Physics (2024). DOI: 10.1038/s41567-024-02765-w; arXiv:2406.14718. https://arxiv.org/abs/2406.14718
  15. Francis A. Bayocboc Jr., Jacek Dziarmaga, Marek M. Rams, Jaka Vodeb, 'Post-Critical Meson Dynamics of Kibble-Zurek Excitations in a 5,564-Qubit Quantum Annealer', arXiv:2607.13842 (2026). https://arxiv.org/abs/2607.13842
  16. Jean-Yves Desaules, Guo-Xian Su, Ian P. McCulloch, Bing Yang, Zlatko Papić, Jad C. Halimeh, 'Ergodicity Breaking Under Confinement in Cold-Atom Quantum Simulators', Quantum 8, 1274 (2024). DOI: 10.22331/q-2024-02-29-1274; arXiv:2301.07717. https://arxiv.org/abs/2301.07717
  17. Vaibhav Sharma, Kaden R. A. Hazzard, 'Meson dynamics from locally exciting a particle-conserving Z2 lattice gauge theory', Quantum 9, 1872 (2025). DOI: 10.22331/q-2025-10-01-1872; arXiv:2503.02791. https://arxiv.org/abs/2503.02791
  18. Ana Hudomal, Ryan Smith, Andrew Hallam, Zlatko Papić, 'Integrability Breaking and Bound States in Google's Decorated XXZ Circuits', PRX Quantum 5, 010316 (2024). DOI: 10.1103/PRXQuantum.5.010316; arXiv:2307.13042. https://arxiv.org/abs/2307.13042
  19. Milena Horvath, Alvise Bastianello, Sudipta Dhar, Rebekka Koch, Yanliang Guo, Jean-Sébastien Caux, Manuele Landini, Hanns-Christoph Nägerl, 'Observing Bethe strings in an attractive Bose gas far from equilibrium', arXiv:2505.10550 (2025). https://arxiv.org/abs/2505.10550
  20. Gesualdo Delfino, Marianna Sorba, 'On unitary time evolution out of equilibrium', Nuclear Physics B 1014, 116587 (2024). DOI: 10.1016/j.nuclphysb.2024.116587; arXiv:2403.13477. https://arxiv.org/abs/2403.13477