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 $\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 $\alpha >0$. Yet, for $\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 $\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 $\alpha \ll 1$ even for large $N$, giving rise to ensemble inequivalence.

• Received 22 December 2020
• Revised 3 August 2021
• Accepted 16 September 2021

DOI:https://doi.org/10.1103/PhysRevB.104.094309

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Open access publication funded by the Max Planck Society.

Published by the American Physical Society