We investigate the rate of thermalization of local operators in the one-dimensional anisotropic antiferromagnetic Heisenberg model with next-nearest neighbor interactions that break integrability. This is done by calculating the scaling of the difference of the diagonal and canonical thermal ensemble values as a function of system size, and by directly calculating the time evolution of the expectation values of the operators with the Chebyshev polynomial expansion. Spatial and spin symmetry is exploited and the Hamiltonian is divided into subsectors according to their symmetry. The rate of thermalization depends on the proximity to the integrable limit. When integrability is weakly broken thermalization is slow, and becomes faster the stronger the next-nearest neighbor interaction is. Three different regimes for the rate of thermalization with respect to the strength of the integrability breaking parameter are identified. These are shown to be directly connected with the relative strength of the low and higher energy difference off-diagonal operator matrix elements in the symmetry eigenbasis of the Hamiltonian. Close to the integrable limit the off-diagonal matrix elements peak at higher energies and high-frequency fluctuations are important and slow down thermalization. Away from the integrable limit a strong low-energy peak gradually develops that takes over the higher frequency fluctuations and leads to quicker thermalization.
We revisit the J1-J2 frustrated Heisenberg spin-1/2 chain with dimerization (δ) or modulation in the nearest-neighbor couplings to investigate its thermalization behavior. While the dimerization tends to induce localization, the next-nearest-neighbor interaction J2 generally favors thermalization, making the assessment of the model's compliance with the Eigenstate Thermalization Hypothesis (ETH) particularly subtle. The challenge is further compounded by the model's SU(2) symmetry; the study of ETH compliance is necessarily done for each symmetry sector but separating different sectors of this symmetry is known to be a computationally demanding task. The current study is driven by two main motivations: first, to explore whether the well-known ground-state phases of the model have any bearing on its thermalization properties; and second, to understand how the interplay between two competing factors, namely, the non-uniformity (via δ) and the beyond-nearest-neighbor interactions (via J2) governs the system's approach to thermal equilibrium. A systematic analysis shows that the ETH is most strongly satisfied for intermediate values of δ (~ 0.5) with J2 ranging from intermediate (~ 0.5) to large (~ 1)- a parameter regime falls within the spiral ground-state phase. It is also found that when the system is in the gapless ground-state phase (which falls within the N'eel phase), the ETH is more prone to violation. In the regime of large δ and small J2, the system is seen to enter a localized phase (characterized here by modulation in density-of-states; assessing ETH compliance is less meaningful for this phase.
We study thermalization slowing down of a quantum many-body spin system upon approach to two distinct integrability limits. Motivated by previous studies of classical systems, we identify two thermalization timescales: one quantum Lyapunov timescale is extracted by quantifying operator growth in time on an appropriately defined basis, while another ergodization timescale is related to the statistics of fluctuations of the time-evolved operator around its mean value based on the eigenstate thermalization hypothesis. Using a paradigmatic quantum Ising chain, we find that both timescales diverge upon approach to integrability. We investigate the relative strength of the divergence in the two limits and find that, despite significant qualitative differences in the mechanism of integrability breaking, the timescales diverge in a similar fashion. This allows us to establish a universality of integrability breaking in quantum spin dynamics.
We study how the proximity to an integrable point or to localization as one approaches the atomic limit, as well as the mixing of symmetries in the chaotic domain, may affect the onset of thermalization in finite one-dimensional systems. We consider systems of hard-core bosons at half-filling with nearest-neighbor hopping and interaction, and next-nearest-neighbor interaction. The latter breaks integrability and induces a ground-state superfluid to insulator transition. By full exact diagonalization, we study chaos indicators and few-body observables. We show that when different symmetry sectors are mixed, chaos indicators associated with the eigenvectors, contrary to those related to the eigenvalues, capture the onset of chaos. The results for the complexity of the eigenvectors and for the expectation values of few-body observables confirm the validity of the eigenstate thermalization hypothesis in the chaotic regime, and therefore the occurrence of thermalization. We also study the properties of the off-diagonal matrix elements of few-body observables in relation to the transition from integrability to chaos and from chaos to localization.
We study Heisenberg spin-1/2 and spin-1 chains with alternating ferromagnetic (
J1F) and antiferromagnetic (
J1A) nearest-neighbor interactions and ferromagnetic next-nearest-neighbor interaction (
J2F). In this model frustration is induced by the non-zero
J2F and the effect of dimerization is realized by changing the strength of
J1A for a fixed value of
J1F. The model with site spin
s behaves like a Haldane spin chain with site spin 2
s in the limit of vanishing
J2F and large
J1F/J1A. We show that the exact ground state of the model can be found along a line in the parameter space. The phase diagram in the space of
J1A−J2F is determined using numerical techniques complemented by some analytical calculations. A number of quantities, including the ground state spin, energy gap, structure factor, entanglement entropy, and zero temperature magnetization, are studied to understand the complete phase diagram. An interesting and application-wise important feature of this model is that it exhibits a macroscopic magnetization jump in the presence of a magnetic field; we study this using an effective Hamiltonian.
The thermal evolution of equal-time spin correlations in the antiferromagnetic isotropic spin-1/2 Heisenberg chain is examined by full exact diagonalization of a ten-site periodic cluster. The calculation retains the complete spectrum and evaluates canonical ensemble averages without decoupling approximations. Correlations through the maximum independent separation are reported over 0.10 ≤ kBT/J ≤ 5. The low-temperature state displays the expected alternating antiferromagnetic pattern and appreciable correlations across the finite ring. Heating suppresses the magnitude of every nonlocal correlator, with longer-range terms disappearing first. The nearest-neighbour correlation approaches the high-temperature limit continuously, while an effective finite-cluster correlation length falls rapidly around kBT/J of order unity. Low- and high-temperature consistency checks, SU (2) symmetry, and comparison with Bethe-ansatz expectations are used to validate the computation. The results provide a compact reproducible reference for interpreting thermal loss of short-range order, while clearly separating finite-size numerical evidence from thermodynamic-limit conclusions.