Quantum systems with strong long-range interactions are thought to robustly resist thermalization because their energy spectra organize into a discrete number of multiplets. We show that this nonthermal behavior is surprisingly fragile to an experimentally ubiquitous perturbation: a staggered magnetic field. Using self-consistent mean-field theory and exact diagonalization, we reveal that the energy spectrum, while composed of discrete subspaces or multiplets, collectively forms a dense spectrum. The equilibration time is independent of system size and depends only on the fluctuations in the initial state. For initial states at low to intermediate energies, the long-time average aligns with the S z sector-weighted microcanonical ensemble. However, for states in the middle of the spectrum, the long-time average depends on the initial state due to quantum scarlike eigenstates localized at unstable points in classical phase space. Our results can be readily tested on a range of experimental platforms, including Rydberg atoms or optical cavities.
Quantum systems with strong long-range interactions are thought to robustly resist thermalization because their energy spectra organize into a discrete number of multiplets. We show that this nonthermal behavior is surprisingly fragile to an experimentally ubiquitous perturbation: a staggered magnetic field. Using self-consistent mean-field theory and exact diagonalization, we reveal that the energy spectrum, while composed of discrete subspaces or multiplets, collectively forms a dense spectrum. The equilibration time is independent of system size and depends only on the fluctuations in the initial state. For initial states at low to intermediate energies, the long-time average aligns with the S z sector-weighted microcanonical ensemble. However, for states in the middle of the spectrum, the long-time average depends on the initial state due to quantum scarlike eigenstates localized at unstable points in classical phase space. Our results can be readily tested on a range of experimental platforms, including Rydberg atoms or optical cavities.
Using a Krylov-subspace time evolution algorithm, we simulate the real-time dynamics of translation invariant non-integrable finite spin rings to quite long times with high accuracy. We systematically study the finite-size deviation between the resulting equilibrium state and the thermal state, and we highlight the importance of the energy variance on the deviations. We find that the deviations are well described by the eigenstate thermalization hypothesis, and that the von Neumann entropy correction scaling is the square of the local operator scaling. We reveal also an area law contribution to the relaxed von Neumann entropy, which we connect to the mutual information between the considered subsystem and its immediate environment. We also find that local observables relax towards equilibrium exponentially with a relaxation time scale that grows linearly with system length and is somewhat independent of the local operator, but depends strongly on the energy of the initial state, with the fastest relaxation times found towards one end of the overall energy spectrum. To contrast this behaviour we also study domain wall initial states, which exhibit clear diffusive behaviour, with a Thouless time scaling quadratically with the systems size, leading to a rather precise estimate for the diffusion constant for states in the vicinity of the middle of the energy spectrum.
There is a dichotomy in the nonequilibrium dynamics of quantum many-body systems. In the presence of integrability, expectation values of local operators equilibrate to values described by a generalized Gibbs ensemble, which retains extensive memory about the initial state of the system. On the other hand, in generic systems such expectation values relax to stationary values described by the thermal ensemble, fixed solely by the energy of the state. At the heart of understanding, this dichotomy is the eigenstate thermalization hypothesis (ETH): individual eigenstates in nonintegrable systems are thermal, in the sense that expectation values agree with the thermal prediction at a temperature set by the energy of the eigenstate. In systems where ETH is violated, thermalization can be avoided. Thus, establishing the range of validity of ETH is crucial in understanding whether a given quantum system thermalizes. Here, we study a simple model with confinement, the quantum Ising chain with a longitudinal field, in which ETH is violated. Despite an absence of integrability, there exist rare (nonthermal) states that persist far into the spectrum. These arise as a direct consequence of confinement: pairs of particles are confined, forming new ``meson'' excitations whose energy can be extensive in the system size. We show that such states are nonthermal in both the continuum and in the low-energy spectrum of the corresponding lattice model. We highlight that the presence of such states within the spectrum has important consequences, with certain quenches leading to an absence of thermalization and local observables evolving anomalously.
Motivated by recent ion experiments on tunable long-range interacting quantum systems [Neyenhuis et al., Sci. Adv. 3, e1700672 (2017)SACDAF2375-254810.1126/sciadv.1700672], we test the strong eigenstate thermalization hypothesis for systems with power-law interactions ∼1/r^{α}. We numerically demonstrate that the strong eigenstate thermalization hypothesis typically holds, at least for systems with α≥0.6, which include Coulomb, monopole-dipole, and dipole-dipole interactions. Compared with short-range interacting systems, the eigenstate expectation value of a generic local observable is shown to deviate significantly from its microcanonical ensemble average for long-range interacting systems. We find that Srednicki's ansatz breaks down for α≲1.0, at least for relatively large system sizes.
We use large-scale exact diagonalization to study the quantum Ising chain in a transverse field with long-range power-law interactions decaying with exponent $\ensuremath{\alpha}$. We numerically study various probes for quantum chaos and eigenstate thermalization on the level of eigenvalues and eigenstates. The level-spacing statistics yields a clear sign towards a Wigner-Dyson distribution and therefore towards quantum chaos across all values of $\ensuremath{\alpha}>0$. Yet, for $\ensuremath{\alpha}<1$ we find that the microcanonical entropy is nonconvex. This is due to the fact that the spectrum is organized in energetically separated multiplets for $\ensuremath{\alpha}<1$. While quantum chaotic behavior develops within the individual multiplets, many multiplets do not overlap and do not mix with each other, as we analytically and numerically argue. Our findings suggest that a small fraction of the multiplets could persist at low energies for $\ensuremath{\alpha}\ensuremath{\ll}1$ even for large N, giving rise to ensemble inequivalence.
Understanding the stability of integrability in many-body quantum systems is key to controlling dynamics and predicting thermalization. While the breakdown of integrability in short-range interacting systems is well understood, the role of long-range couplings -- ubiquitous and experimentally realizable -- remains unclear. We show that in fully connected models, integrability is either robust or extremely fragile, depending on whether perturbations are non-extensive, extensive one-body, or extensive two-body. In contrast to finite short-range systems, where any of these perturbations can induce chaos at finite strength, in fully connected finite models, chaos is triggered by extensive two-body perturbations and even at infinitesimal strength. Chaos develops within energy bands defined by symmetries, leading to a fragmented realization of the eigenstate thermalization hypothesis and clarifying how microcanonical shells can be constructed in such systems. We also introduce a general symmetry-based framework that explains the stability of integrability.
Quantum many-body scars (QMBS) are exceptional energy eigenstates of quantum many-body systems associated with violations of thermalization for special nonequilibrium initial states. Their various systematic constructions require fine-tuning of local Hamiltonian parameters. In this work, we demonstrate that long-range interacting quantum spin systems generically host robust QMBS. We analyze spectral properties upon raising the power-law decay exponent α of spin-spin interactions from the solvable permutationally symmetric limit α = 0 . First, we numerically establish that, despite the fact that spectral signatures of chaos appear for infinitesimal α , the towers of α = 0 energy eigenstates with large collective spin are smoothly deformed as α is increased and exhibit characteristic QMBS features. To elucidate the nature and fate of these states in larger systems, we introduce an analytical approach based on mapping the spin Hamiltonian onto a relativistic quantum rotor nonlinearly coupled to an extensive set of bosonic modes. We analytically solve for the eigenstates of this interacting impurity model by means of a novel polaron-type canonical transformation and show their self-consistent localization in large-spin sectors of the original Hamiltonian for 0 < α < d (with d = spatial dimension of the lattice ). Our theory unveils the stability mechanism of such QMBS for an arbitrary system size and predicts instances of its breakdown, e.g., near dynamical critical points or in the presence of semiclassical chaos, which we verify numerically in long-range quantum Ising chains. As a by-product, we find a predictive criterion for the presence or absence of heating under periodic driving for 0 < α < d , beyond existing Floquet-prethermalization theorems.