We show that the emergence of a ``magic barrier'', characterized by a transient maximum of nonstabilizerness during rapid entanglement growth, is a universal feature of thermalization within quantum ergodic theories, which is intimately connected to random matrix spectral properties. Starting from a pure product state, we analytically determine the time evolution of entanglement entropy, antiflatness, capacity of entanglement, and nonlocal trace distance magic under unitary dynamics generated by random matrix Hamiltonians drawn from all three Gaussian ensembles. We find that the transient maxima of antiflatness, capacity of entanglement, and nonlocal magic occur on the same characteristic timescale as the maximum rate of entanglement growth. These results establish that the magic barrier is a universal feature of RMT thermalization, suggesting universality across quantum chaotic systems. We conclude that quantitative simulations of thermalization in ergodic quantum theories inherently require quantum computation.
We show that the emergence of a ``magic barrier'', characterized by a transient maximum of nonstabilizerness during rapid entanglement growth, is a universal feature of thermalization within quantum ergodic theories, which is intimately connected to random matrix spectral properties. Starting from a pure product state, we analytically determine the time evolution of entanglement entropy, antiflatness, capacity of entanglement, and nonlocal trace distance magic under unitary dynamics generated by random matrix Hamiltonians drawn from all three Gaussian ensembles. We find that the transient maxima of antiflatness, capacity of entanglement, and nonlocal magic occur on the same characteristic timescale as the maximum rate of entanglement growth. These results establish that the magic barrier is a universal feature of RMT thermalization, suggesting universality across quantum chaotic systems. We conclude that quantitative simulations of thermalization in ergodic quantum theories inherently require quantum computation.
We investigate the time dependence of antiflatness in the entanglement spectrum, a measure for nonstabilizerness and lower bound for nonlocal quantum magic resource, on a subsystem of a linear SU(2) plaquette chain during thermalization. Tracing the time evolution of a large number of initial states, we find that the antiflatness exhibits a barrierlike maximum during the time period when the entanglement entropy of the subsystem grows rapidly from the initial value to the microcanonical entropy. The location of the peak is strongly correlated with the time when the entanglement exhibits the strongest growth. This behavior is found for generic highly excited initial computational basis states and persists for coupling constants across the ergodic regime, revealing a universal structure of the entanglement spectrum during thermalization. We conclude that quantitative simulations of thermalization for non-Abelian gauge theories require quantum computing. We speculate that this property generalizes to other quantum chaotic systems, a conjecture supported by analogous behavior observed in real-time simulations of the mixed-field Ising model.
Quantum entanglement and magic are complementary resources underlying quantum computational advantage, yet their dynamical relation in many-body systems remains poorly understood. In this Letter, we show that the mechanism of bipartite entanglement growth is encoded in the relative timescale between the entropy-growth-rate peak and the magic barrier, defined as the transient peak of the anti-flatness of the entanglement spectrum. When entanglement is locally built, the same microscopic process increases the entropy and reshapes the Schmidt spectrum, so the magic-barrier peak occurs in the time window of maximal entropy growth. When entanglement is mainly transported or redistributed, entropy can grow before appreciable spectral non-flatness is generated, naturally separating the two peak times. We demonstrate this distinction in the random-field XXZ chain: the two peaks remain strongly correlated in the thermal regime, while their separation grows systematically across the thermal--MBL crossover. We further validate this theoretical framework by employing Bell-pair initial states alongside a tunable SWAP--Haar random circuit. Our results reveal an intrinsic dynamical connection between entanglement and magic, establishing the magic barrier as a powerful spectral diagnostic of how quantum information is generated, transported, and reshaped.
Entanglement and magic are resources that reveal complementary aspects of quantum many-body systems. Their interplay is captured by nonlocal magic, the magic that survives arbitrary local changes of basis. Yet their markedly different dynamical behavior leaves open how this irreducible component of magic spreads. Here we connect nonlocal magic to the capacity of entanglement, a tractable quantity measuring fluctuations of the entanglement Hamiltonian. This connection enables analytically controlled predictions across a wide range of many-body dynamics, from chaotic to integrable systems, which we investigate also using large-scale numerical simulations. In chaotic systems, nonlocal magic exhibits a transient buildup, with logarithmic growth in time followed by decay to a size-independent value. For integrable systems, we develop a quasiparticle picture in quantitative agreement with numerics, showing that the same initial growth instead leads to saturation at a value logarithmic in subsystem size. These contrasting behaviors have an operational consequence for entanglement embezzlement: the strongest scramblers embezzle only transiently, whereas free-fermionic dynamics can sustain universal embezzlement in the steady state.
We show that, starting from a product state, local unitary quantum evolutions generate intermediate states which exhibit a multiscale entanglement-spectrum structure required for universal entanglement embezzlement. This facilitates entanglement extraction from a catalyst many-body state while leaving it asymptotically unchanged. Remarkably, these atypical structures emerge generically at intermediate stages, well before reaching maximum entropy where thermalization has flattened out the spectral hierarchy. The resulting state is accompanied by nonlocal nonstabilizerness that diverges with the system size, consistent with a recently established equivalence between universal embezzlement and divergent nonlocal magic. Thus, without any fine tuning, a chaotic quantum evolution generates intermediate states which form a universal family of catalytic reservoirs.
Thermalization conventionally describes local properties of isolated many-body systems at equilibrium. Magic, or nonstabilizerness $\unicode{x2013}$ the resource enabling universal quantum computation $\unicode{x2013}$ is by contrast encoded in the global structure of the many-body wavefunction. We show that, despite its global nature, the magic of equilibrium pure states of chaotic many-body systems, including late-time evolved states and energy eigenstates, is universal: it is captured by the thermal Scrooge ensemble, the minimally informative ensemble of pure states consistent with the Gibbs state at the same effective temperature. Therefore, for systems with no conserved quantities other than the total energy, equilibrium magic is a function of temperature alone, independent of the initial state and other microscopic features of the equilibrium state. This yields concrete universal predictions for the stabilizer Rényi entropies (SREs). At infinite temperature, the SRE is set by Haar-like fluctuations of the Pauli spectrum, while at finite temperature energy conservation induces a volume-law thermodynamic correction controlled by the thermal Pauli spectrum. We support these predictions with analytical arguments and extensive numerical simulations. We further show that chaotic many-body systems at high temperatures possess long-range magic and entanglement that cannot be removed by finite-depth local quantum circuits. Our results establish magic as a thermodynamic property of chaotic many-body systems and suggest that Scrooge ensembles may provide a unified framework for quantum many-body resources.
Under the framework of unified time scale, boundary time geometry, Matrix Universe THE--MATRIX, and Quantum Cellular Automaton (QCA) Universe, we construct a structural unified theory for ``Quantum Chaos and Eigenstate Thermalization''. The Bohigas--Giannoni--Schmit (BGS) conjecture suggests that spectral statistics of quantum chaotic systems follow Random Matrix Theory (RMT), while the Eigenstate Thermalization Hypothesis (ETH) explains thermalization of isolated quantum systems from the perspective of matrix elements of local observables. Based on the unified time scale mother formula equation \kappa(\omega) =\varphi'(\omega){\pi} =\rho_{rel}(\omega) =1{2\pi}tr Q(\omega), equation this paper interprets quantum chaos and thermalization as: ergodicity of scattering phase flow in the Matrix Universe and operator mixing in the QCA Universe; and further regards spectral rigidity and RMT universality as statistical inevitable results of unified time scale density under maximum entropy principle. Specifically: 1. In the Matrix Universe, defining the ``scattering phase flow'' dynamical system on the energy shell, we prove that when the scattering matrix S(\omega) exhibits chaotic behavior (positive Lyapunov exponent or mixing), the spectral correlation function of the unified time scale \kappa(\omega) asymptotically approaches the RMT (GOE/GUE) kernel function equation R_2(\omega_1,\omega_2) \simeq 1-\Bigl(\sin[\pi\bar\kappa(\omega_1-\omega_2)]{\pi\bar\kappa(\omega_1-\omega_2)}\Bigr)^2, equation thereby providing a scattering--spectral unified proof for the BGS conjecture. 2. In the QCA Universe, considering the operator algebra \mathcal A_{qloc} of local cell observables and the Heisenberg evolution generated by the unitary update U, we prove that for generic chaotic QCA (satisfying algebraic mixing condition), local operators spread to global operators under time evolution, and their matrix elements in energy eigenstates satisfy the ETH ansatz equation \langle E_a|\mathcal O|E_b\rangle \simeq \mathcal O(E)\delta_{ab} +e^{-S(E)/2}f_{\mathcal O}(E,\Delta E)R_{ab}, equation where R_{ab} is a random variable, establishing the microscopic mechanism of thermalization in the discrete universe. 3. Investigating the growth of entanglement entropy S_{ent}(t) in QCA time evolution, we find that under unified time scale control, the entanglement growth rate is bounded by the Lieb--Robinson velocity and related to the ``boundary time area law''. For chaotic systems, linear growth S_{ent}(t)\propto t is consistent with the entropy production rate of the scattering matrix. 4. Under unified time scale constraints, writing the spectral form factor (SFF) K(t) as the Fourier transform of the two-point correlation of \kappa(\omega), we demonstrate the ``slope--dip--plateau'' structure of SFF in the Matrix--QCA universe, interpreting the linear ramp region as a direct manifestation of the spectral rigidity of the unified time scale. 5. The appendix systematically organizes: standard definitions of quantum chaos, RMT spectral statistics and SFF; operator spreading and out-of-time-order correlator (OTOC) description in QCA; rigorous formulation of ETH; and the maximum entropy derivation of unified time scale spectral correlations. Results indicate a purely structural unified picture: quantum chaos and eigenstate thermalization can be viewed as the statistical behavior of scattering phase flow in the Matrix Universe and operator mixing in the QCA Universe under the unified time scale, thereby answering the structural version of ``why quantum chaos and thermalization occur'' within the unified universe mother structure.