We compute exactly the von Neumann entanglement entropy of the eta-pairing states - a large set of exact excited eigenstates of the Hubbard Hamiltonian. For the singlet eta-pairing states the entropy scales with the logarithm of the spatial dimension of the (smaller) partition. For the eta-pairing states with finite spin magnetization density, the leading term can scale as the volume or as the area-times-log, depending on the momentum space occupation of the Fermions with flipped spins. We also compute the corrections to the leading scaling. In order to study the eigenstate thermalization hypothesis (ETH), we also compute the entanglement Rényi entropies of such states and compare them with the corresponding entropies of thermal density matrix in various ensembles. Such states, which we find violate strong ETH, may provide a useful platform for a detailed study of the time-dependence of the onset of thermalization due to perturbations which violate the total pseudospin conservation.
We compute exactly the von Neumann entanglement entropy of the eta-pairing states - a large set of exact excited eigenstates of the Hubbard Hamiltonian. For the singlet eta-pairing states the entropy scales with the logarithm of the spatial dimension of the (smaller) partition. For the eta-pairing states with finite spin magnetization density, the leading term can scale as the volume or as the area-times-log, depending on the momentum space occupation of the Fermions with flipped spins. We also compute the corrections to the leading scaling. In order to study the eigenstate thermalization hypothesis (ETH), we also compute the entanglement Rényi entropies of such states and compare them with the corresponding entropies of thermal density matrix in various ensembles. Such states, which we find violate strong ETH, may provide a useful platform for a detailed study of the time-dependence of the onset of thermalization due to perturbations which violate the total pseudospin conservation.
A prime characterization of many-body localized (MBL) systems is the entanglement of their eigenstates; in contrast to the typical ergodic phase whose eigenstates are volume law, MBL eigenstates obey an area law. In this paper, we show that a spin-disordered Hubbard model has both a large number of area-law eigenstates as well as a large number of eigenstates whose entanglement scales logarithmically with system size (log-law). This model provides a microscopic Hamiltonian which is neither ergodic nor MBL. We establish these results through a combination of analytic arguments based on the eta-pairing operators combined with a numerical analysis of eigenstates. In addition, we describe and simulate a dynamic time evolution approach starting from product states through which one can separately probe the area-law and log-law eigenstates in this system.
We construct exact eigenstates of multicomponent Hubbard models in arbitrary dimensions by generalizing the η -pairing mechanism. Our models include the SU ( N ) Hubbard model as a special case. Unlike the conventional two-component case, the generalized η -pairing mechanism permits the construction of eigenstates that feature off-diagonal long-range order and magnetic long-range order. These states form fragmented fermionic condensates due to a simultaneous condensation of multicomponent η pairs. While the η -pairing states in the SU(2) Hubbard model are based on the η -pairing symmetry, the exact eigenstates in the N -component system with N ≥ 3 arise not from symmetry of the Hamiltonian but from a spectrum generating algebra defined in a restricted Hilbert space. We exploit this fact to show that the generalized η -pairing eigenstates do not satisfy the eigenstate thermalization hypothesis and serve as quantum many-body scar states. This result indicates a weak breakdown of ergodicity in the N -component Hubbard models for N ≥ 3 . Furthermore, we show that these exact eigenstates constitute integrable subsectors in which the Hubbard Hamiltonian effectively reduces to a noninteracting model. This partial integrability causes various multicomponent Hubbard models to weakly break ergodicity. We propose a method of harnessing dissipation to distill the integrable part of the dynamics and elucidate a mechanism of nonthermalization caused by dissipation. This work establishes the coexistence of off-diagonal long-range order and SU ( N ) magnetism in excited eigenstates of the multicomponent Hubbard models, which presents a possibility of novel out-of-equilibrium pairing states of multicomponent fermions. These models unveil a unique feature of quantum thermalization of multicomponent fermions, which can experimentally be tested with cold-atom quantum simulators. Published by the American Physical Society 2024
law is enhanced by a logarithmic factor. In Ref. [1], we use multi-dimensional bosonization to provide a simple derivation of this result, and show that the logarithimic factor has a 1D origin. More importantly the bosonization technique allows us to take into account the Fermi liquid interactions, and obtain the leading scaling behavior of the entanglement entropy of Fermi liquids. The central result of our work is that Fermi liquid interactions do not alter the leading scaling behavior of the entanglement entropy, and the logarithmic enhancement of area law is a robust property of the Fermi liquid phase. In sharp contrast to the fermioic systems with Fermi surfaces, quantum critical (or gapless) bosonic systems do not violate the area law above 1D (except for the case discussed below). The fundamental difference lies in the fact that gapless excitations live near a single point (usually origin of momentum space) in such bosonic systems, while they live around an (extended) Fermi surface in Fermi liquids. In Ref. [2], we studied entanglement properties of some specific examples of the so called Bose metal states, in which bosons neither condense (and become a superfluid) nor localize (and insulate) at T=0. The system supports gapless excitations around “Bose surfaces,” instead of isolated points in momentum space. We showed that similar to free Fermi gas and Fermi liquids, these states violate the entanglement area law in a logarithmic fashion. Compared to ground states, much less is known concretely about entanglement in (highly) excited states. Going back to free fermion systems, in Ref. [3] we show that there exists a duality relation between ground and excited states, and the area law obeyed by ground state turns into a volume law for excited states, something widely expected but hard to prove. Most importantly, we find in appropriate limits the reduced density matrix of a subsystem takes the form of thermal density matrix, providing an explicit example of the eigenstate thermalization hypothesis. Our work [3] explicitly demonstrates how statistical physics emerges from entanglement in a single eigenstate.
With increasing subsystem size and energy, bipartite entanglement entropies of energy eigenstates cross over from the groundstate scaling to a volume law. In previous work, we pointed out that, when strong or weak eigenstate thermalization (ETH) applies, the entanglement entropies of all or, respectively, almost all eigenstates follow a single crossover function. The crossover functions are determined by the subsystem entropy of thermal states and assume universal scaling forms in quantum-critical regimes. This was demonstrated by field-theoretical arguments and the analysis of large systems of non-interacting fermions and bosons. Here, we substantiate such scaling properties for integrable and non-integrable interacting spin-1/2 chains at criticality using exact diagonalization. In particular, we analyze XXZ and transverse-field Ising models with and without next-nearest-neighbor interactions. Indeed, the crossover of thermal subsystem entropies can be described by a universal scaling function following from conformal field theory. Furthermore, we analyze the validity of ETH for entanglement in these models. Even for the relatively small system sizes that can be simulated, the distributions of eigenstate entanglement entropies are sharply peaked around the subsystem entropies of the corresponding thermal ensembles.
For the typical quantum many-body systems that obey the eigenstate thermalization hypothesis (ETH), we argue that the entanglement entropy of (almost) all energy eigenstates is described by a single crossover function. The ETH implies that the crossover functions can be deduced from subsystem entropies of thermal ensembles and have universal properties. These functions capture the full crossover from the ground-state entanglement regime at low energies and small subsystem size (area or log-area law) to the extensive volume-law regime at high energies or large subsystem size. For critical one-dimensional systems, a universal scaling function follows from conformal field theory and can be adapted for nonlinear dispersions. We use it to also deduce the crossover scaling function for Fermi liquids in d>1 dimensions. The analytical results are complemented by numerics for large noninteracting systems of fermions in d≤3 dimensions and have also been confirmed for bosonic systems and nonintegrable spin chains.
The eta-pairing state is the eigenstate of the hypercubic Hubbard model, which exhibits anomalous logarithmic scaling of entanglement entropy. In multi-band systems, eta-pairing can be exact eigenstate when the band is flat without interband coupling. However, typical flatband systems such as Lieb and Kagome lattices often feature band touchings, where interband coupling effects are non-negligible. Using the Creutz ladder, we investigate the deformation of eta-pairing states under the interband coupling effect. Our results show corrections to entanglement entropy scaling, with modified eta-pairing states displaying broadened doublons, nonuniform energy spacing, and deviations from exact behavior for configurations with more than one eta-pair, even in the large band gap limit, except at t = 0. Through a Schrieffer-Wolff transformation, we quantify corrections to the spectrum generating algebra, offering insights into the interplay between interaction-driven phenomena and band structure effects. These findings illuminate the robustness and limitations of eta-pairing in realistic flatband systems.