Generalized Gibbs Ensemble from Eigenstate Entanglement Hamiltonian

Relaxed quantum systems with conservation laws are believed to be approximated by the generalized Gibbs ensemble (GGE), which incorporates the constraints of certain conserved quantities serving as integrals of motion. By drawing an analogy between eigenstate reduced density matrix and GGE, we conjecture that a natural set of conserved quantities for GGE can emerge from the reduced density matrices of properly chosen eigenstates by the entanglement Hamiltonian superdensity matrix (EHSM) framework, and we demonstrate this explicitly for models mappable to free fermions. The framework proposes that such conserved quantities are linear superpositions of eigenstate entanglement Hamiltonians of a larger auxiliary system, where the eigenstates are Fock states occupying what we call the common eigenmodes, which remain eigenmodes when truncated within the physical subsystem. For 1D homogeneous free fermions with (anti)periodic boundary conditions, which maps to 1D hardcore bosons with nearest neighbor hoppings, these conserved quantities lead to a non-Abelian GGE, which predicts the relaxation of both fermion and boson bilinears more accurately than the conventional Abelian GGE. Generalization of this framework may provide novel numerical insights for quantum integrability.