Eigenstate Entanglement Entropy in Random Quadratic Hamiltonians

Single-quasiparticle eigenstate thermalization

Quadratic Hamiltonians that exhibit single-particle quantum chaos are called quantum-chaotic quadratic Hamiltonians. One of their hallmarks is single-particle eigenstate thermalization introduced in Łydżba et al. [Phys. Rev. B 104, 214203 (2021)2469-995010.1103/PhysRevB.104.214203], which describes statistical properties of matrix elements of observables in single-particle eigenstates. However, the latter has been studied only in quantum-chaotic quadratic Hamiltonians that obey the U(1) symmetry. Here, we focus on quantum-chaotic quadratic Hamiltonians that break the U(1) symmetry and, hence, their "single-particle" eigenstates are actually single-quasiparticle excitations introduced on the top of a many-body state. We study their wave functions and matrix elements of one-body observables, for which we introduce the notion of single-quasiparticle eigenstate thermalization. Focusing on spinless fermion Hamiltonians in three dimensions with local hopping, pairing, and on-site disorder, we also study the properties of disorder-induced near zero modes, which give rise to a sharp peak in the density of states at zero energy. Finally, we numerically show equilibration of observables in many-body eigenstates after a quantum quench. We argue that the latter is a consequence of single-quasiparticle eigenstate thermalization, in analogy to the U(1) symmetric case from Łydżba et al. [Phys. Rev. Lett. 131, 060401 (2023)0031-900710.1103/PhysRevLett.131.060401].

Eigenstate entanglement entropy in the integrable spin-1/2 XYZ model

We study the average and the standard deviation of the entanglement entropy of highly excited eigenstates of the integrable interacting spin-1/2 XYZ chain away from and at special lines with U(1) symmetry and supersymmetry. We universally find that the average eigenstate entanglement entropy exhibits a volume-law coefficient that is smaller than that of quantum-chaotic interacting models. At the supersymmetric point, we resolve the effect that degeneracies have on the computed averages. We further find that the normalized standard deviation of the eigenstate entanglement entropy decays polynomially with increasing system size, which we contrast with the exponential decay in quantum-chaotic interacting models. Our results provide state-of-the art numerical evidence that integrability in spin-1/2 chains reduces the average and increases the standard deviation of the entanglement entropy of highly excited energy eigenstates when compared with those in quantum-chaotic interacting models.

Generalized Thermalization in Quantum-Chaotic Quadratic Hamiltonians

Thermalization (generalized thermalization) in nonintegrable (integrable) quantum systems requires two ingredients: equilibration and agreement with the predictions of the Gibbs (generalized Gibbs) ensemble. We prove that observables that exhibit eigenstate thermalization in single-particle sector equilibrate in many-body sectors of quantum-chaotic quadratic models. Remarkably, the same observables do not exhibit eigenstate thermalization in many-body sectors (we establish that there are exponentially many outliers). Hence, the generalized Gibbs ensemble is generally needed to describe their expectation values after equilibration, and it is characterized by Lagrange multipliers that are smooth functions of single-particle energies.

Average entanglement entropy of midspectrum eigenstates of quantum-chaotic interacting Hamiltonians

To which degree the average entanglement entropy of midspectrum eigenstates of quantum-chaotic interacting Hamiltonians agrees with that of random pure states is a question that has attracted considerable attention in the recent years. While there is substantial evidence that the leading (volume-law) terms are identical, which and how subleading terms differ between them is less clear. Here we carry out state-of-the-art full exact diagonalization calculations of clean spin-1/2 XYZ and XXZ chains with integrability breaking terms to address this question in the absence and presence of U(1) symmetry, respectively. We first introduce the notion of maximally chaotic regime, for the chain sizes amenable to full exact diagonalization calculations, as the regime in Hamiltonian parameters in which the level spacing ratio, the distribution of eigenstate coefficients, and the entanglement entropy are closest to the random matrix theory predictions. In this regime, we carry out a finite-size scaling analysis of the subleading terms of the average entanglement entropy of midspectrum eigenstates when different fractions ν of the spectrum are included in the average. We find indications that, for ν→0, the magnitude of the negative O(1) correction is only slightly greater than the one predicted for random pure states. For finite ν, following a phenomenological approach, we derive a simple expression that describes the numerically observed ν dependence of the O(1) deviation from the prediction for random pure states.

Chaos and bipartite entanglement between Bose-Josephson junctions

The entanglement between two weakly coupled bosonic Josephson junctions is studied in relation to the classical mixed phasespace structure of the system, containing symmetry-related regular islands separated by chaos. The symmetry-resolved entanglement spectrum and bipartite entanglement entropy of the system's energy eigenstates are calculated and compared to their expected structure for random states that exhibit complete or partial ergodicity. The entanglement spectra of chaos-supported eigenstates match the microcanonical structure of a Generalized Gibbs Ensemble due to the existence of an adiabatic invariant that restricts ergodization on the energy shell. The symmetry-resolved entanglement entropy of these quasistochastic states consists of a mean-field maximum entanglement term and a fluctuation correction due to the finite size of the constituent subsystems. The total bipartite entanglement entropy of the eigenstates correlates with their chaoticity. Island-supported eigenstates are macroscopic Schrödinger cat states for particles and excitations with substantially lower entanglement.

Tight-binding billiards

Recent works have established universal entanglement properties and demonstrated validity of single-particle eigenstate thermalization in quantum-chaotic quadratic Hamiltonians. However, a common property of all quantum-chaotic quadratic Hamiltonians studied in this context so far is the presence of random terms that act as a source of disorder. Here we introduce tight-binding billiards in two dimensions, which are described by noninteracting spinless fermions on a disorder-free square lattice subject to curved open (hard-wall) boundaries. We show that many properties of tight-binding billiards match those of quantum-chaotic quadratic Hamiltonians: The average entanglement entropy of many-body eigenstates approaches the random matrix theory predictions and one-body observables in single-particle eigenstates obey the single-particle eigenstate thermalization hypothesis. On the other hand, a degenerate subset of single-particle eigenstates at zero energy (i.e., the zero modes) can be described as chiral particles whose wave functions are confined to one of the sublattices.

Entanglement of midspectrum eigenstates of chaotic many-body systems: Reasons for deviation from random ensembles

Eigenstates of local many-body interacting systems that are far from spectral edges are thought to be ergodic and close to being random states. This is consistent with the eigenstate thermalization hypothesis and volume-law scaling of entanglement. We point out that systematic departures from complete randomness are generically present in midspectrum eigenstates, and focus on the departure of the entanglement entropy from the random-state prediction. We show that the departure is (partly) due to spatial correlations and due to orthogonality to the eigenstates at the spectral edge, which imposes structure on the midspectrum eigenstates.

Entanglement distribution in the quantum symmetric simple exclusion process

We study the probability distribution of entanglement in the quantum symmetric simple exclusion process, a model of fermions hopping with random Brownian amplitudes between neighboring sites. We consider a protocol where the system is initialized in a pure product state of M particles, and we focus on the late-time distribution of Rényi-q entropies for a subsystem of size ℓ. By means of a Coulomb gas approach from random matrix theory, we compute analytically the large-deviation function of the entropy in the thermodynamic limit. For q>1, we show that, depending on the value of the ratio ℓ/M, the entropy distribution displays either two or three distinct regimes, ranging from low to high entanglement. These are connected by points where the probability density features singularities in its third derivative, which can be understood in terms of a transition in the corresponding charge density of the Coulomb gas. Our analytic results are supported by numerical Monte Carlo simulations.