2026-05-08 | | Total: 7
Quantum thermalization describes how interacting quantum systems relax toward thermal equilibrium, a central problem in modern physics. Yet most experimental information on many-body systems comes from short-time transition spectroscopy, typically interpreted within Kubo's linear-response framework. These perspectives - long-time equilibration versus short-time response - seem fundamentally disconnected. Here we establish an exact link between them: the Kubo-Thermalization correspondence, which connects long-time thermalized magnetization under weak driving to short-time linear-response spectra for a spin coupled to a thermal bath. The correspondence holds even when the steady state differs substantially from the initial state and when each regime is individually difficult to describe theoretically. We experimentally confirm the correspondence using effective spin-1/2 impurities realized with ultracold fermions in two internal states coupled to a Fermi sea. Our results provide a rare exact statement about quantum thermalization and offer a novel route to infer thermalization dynamics from equilibrium response measurements in strongly interacting quantum systems, independent of microscopic details of the system-bath coupling.
Hexagonal optical lattices, emulating graphene and hexagonal boron nitride (h-BN) structures, provide a versatile platform for exploring strongly correlated quantum matter. Using continuous-space exact diagonalization and quantum Monte Carlo simulations, we investigate the phase diagrams of ultracold bosons in honeycomb and h-BN lattices. For the honeycomb lattice, we find significant deviations from the standard Bose-Hubbard model even for strong lattice amplitudes. We observe suppressed Mott insulator lobes and the absence of higher-order insulating phases, attributed to strong density-assisted tunneling effects. In the h-BN case, a rich phase diagram emerges, featuring multiple Mott lobes with various sublattice occupations, driven by the interplay of lattice asymmetry, interactions, and particle filling. Our results highlight the necessity of continuous-space treatments for capturing the full complexity of bosonic quantum phases in hexagonal geometries, paving the way for experimental realizations with ultracold atoms and further theoretical work.
Quantum many-body scars enable persistent non-ergodic dynamics in otherwise thermalizing systems, yet their stabilization typically relies on fine-tuned initial states or engineered Hamiltonian perturbations. Here we show that lattice geometry alone can serve as a powerful and experimentally accessible control knob for inducing and enhancing scarring. By transforming a one-dimensional chain into a quasi-one-dimensional triangle-decorated lattice, we find that the fully polarized state -- normally thermalizing in the PXP model -- exhibits pronounced fidelity revivals, slow entanglement growth, and strong overlap with a tower of weakly entangled eigenstates. We trace this behavior to a geometry-induced restructuring of the constrained Hilbert space, whereby the adjacency graph decomposes into hypercube subgraphs that enforce coherent population transfer and stabilize an emergent approximate $\mathrm{su}(2)$ algebra. We propose a direct implementation in programmable arrays of tweezer-trapped Rydberg atoms, where the triangle-decorated geometry can be realized using spatial light modulators and the resulting scarring dynamics probed via time-resolved measurements of excitation density. Our results establish lattice connectivity as a design principle for engineering non-ergodic dynamics in constrained quantum systems.
Many-body localized (MBL) systems are known to thermalize in periodically driven systems. In this work, we demonstrate that under proper driving protocol, this thermalization this thermalization can be resisted such that the MBL phase turns into a non-ergodic extended phase, known as the many-body critical (MBC) phase. Considering a kicked quasiperiodic Ising chain, we show that while at high-frequency driving the ergodic, MBL, and the MBC phases coexist, at moderate driving frequencies the MBL phase is completely suppressed and the MBC phase proliferates in the parameter space. Using quasienergy statistics, Floquet eigenstates, autocorrelation dynamics, and entanglement growth, we characterize the emergent phases and identify non-monotonic signatures revealing richness of the nonergodic phases. Our results establish Floquet driving as a powerful route to stabilizing nonergodic extended many-body phases beyond the conventional Floquet-MBL paradigm.
Periodically driven quantum systems, known as Floquet systems, provide a versatile platform for engineering novel topological phases absent in static settings. However, dynamically characterizing these non-equilibrium topological invariants remains a challenge. Here, we develop a Floquet perturbation theory in the extended Hilbert space to analytically describe the center-of-mass (CoM) dynamics of a wave packet. When applied to the driven Su-Schrieffer-Heeger model, our theory reveals that the CoM exhibits multi-frequency Zitterbewegung oscillations, whose spectral composition and phase are directly tied to the system's Floquet band structure. Crucially, we find that band inversions at topological phase transitions imprint distinct signatures in the CoM dynamics, including the emergence of low-frequency modes and phase shifts of the oscillatory trajectory. These dynamical signatures offer a practical protocol for detecting Floquet topological invariants, which we demonstrate for both high-frequency and strongly driven regimes. Our work establishes CoM dynamics as a simple and experimentally accessible probe for exploring topological phase transitions in Floquet systems.
Models of many-body localization (MBL) exhibit slow numerical drifts towards delocalization with increasing system size, for which no satisfactory theory exists. Numerics indicates that these drifts are driven by the proliferation of many-body resonances at intermediate disorder strengths. We develop a statistical method to predict the distribution of resonance oscillation frequencies which captures how the formation of resonances at larger frequency scales subsequently affects the formation of resonances at lower frequencies. Working within the statistical Jacobi approximation (SJA), we derive a flow equation for a power-law exponent $θ(w)$ characterizing the density of resonances at frequency scale $w$. A localized phase is described by a line of fixed points with $θ(w)>0$, while an instability of the flow signals resonance proliferation and the onset of thermalization. The predicted $θ(w)$ matches numerics on the Anderson model on random regular graphs and the Lévy-Rosenzweig-Porter random matrix ensemble, both of which host resonance-driven delocalization transitions. We further connect the flow to eigenstate properties such as the participation ratio and to dynamical observables such as the return probability. The predicted $θ(w)$ also matches what is numerically measured in real-space models of MBL at intermediate disorder strengths, representing a significant step towards explaining the finite-size drifts observed in MBL.
In a free Fermi gas at temperature $T$ much higher than the Fermi temperature one expects that the fluctuations of the number of particles in a given region has Poissonian/classical statistics. On the other hand at low temperature the Pauli exclusion principle leads to non trivial counting statistics. It is of great interest from a theoretical and experimental point of view to characterize the crossover between these two limits. Here we focus on the hole probability $P(R,T)$, i.e. the probability that a region of size $R$ is devoid of particles, in dimension $d$, and on the case of a spherical region of large radius $R$. We show that at low temperature it takes the scaling form $P(R,T)\sim \exp\big[-(k_F R)^{d+1}Φ_d(u=2R\,T/k_F)\big],$ where $k_F$ is the Fermi momentum. By mapping the problem to an effective Coulomb gas, we compute exactly the scaling function $Φ_d(u)$ in any dimension. Remarkably, it exhibits a transition of order $\tfrac{3}{2}(d+1)$ at the universal critical value $u_c=2/π$, signaling a sharp change in the mechanism of rare fluctuations, associated with the emergence of a macroscopic gap in the optimal density of the associated Coulomb gas. Our analytical predictions are supported by precise numerical evaluations of the corresponding Fredholm determinants.