Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2210.05288 · International Symposium on Lattice Field Theory (LATTICE)

Quantum chaos in supersymmetric Yang-Mills-like model: equation of state, entanglement, and spectral form-factors

P. V. Buividovich

Published Nov 7, 2022• 4 citations
arXivHTMLPDF

Abstract

We analyze in detail a sharp transition between the low-energy, low-dimensional eigenstates and the high-energy chaotic bulk of the spectrum for a simple supersymmetric quantum-mechanical model with Hamiltonian $\hat{H}_S = \left(\hat{p}_1^2 + \hat{p}_2^2 + \hat{x}_1^2 \, \hat{x}_2^2\right) \otimes I + \hat{x}_1 \otimes \sigma_1 + \hat{x}_2 \otimes \sigma_3$, which mimics the structure of the Banks-Fischler-Susskind-Stanford (BFSS) matrix model, the spatially compactified $\mathcal{N} = 1$ super-Yang-Mills theory. We conjecture that this transition might be similar to the transition between the $D0$-brane and $M$-theory regimes in the BFSS model, and find that it does not lead to irregularities in the thermodynamic equation of state. We demonstrate that real-time spectral form-factor for our supersymmetric model exhibits the ``ramp'' behavior typical for quantum chaos. We also analyze the entanglement entropy and the spectrum of the reduced density matrix of the eigenstates of $\hat{H}_S$, considering one of the bosonic degrees of freedom as a subsystem. The entanglement entropy of low-energy eigenstates appears to be practically energy-independent. Exactly at the onset of random-matrix-type level spacing fluctuations, this behavior rapidly changes into a steady growth of entanglement with energy. We demonstrate that the spectrum of the reduced density matrix also exhibits universal level-spacing fluctuations towards its higher end, even for the ground state of the supersymmetric model. Thus even the regularly spaced, non-chaotic eigenstates contain some information about semi-classical chaotic dynamics at high energies.

DOI
10.22323/1.430.0246

Related papers

Quantum chaos in supersymmetric quantum mechanics: An exact diagonalization study

P. V. Buividovich

We use exact diagonalization to study energy level statistics and out-of-time-order correlators (OTOCs) for the simplest supersymmetric extension ${\stackrel{^}{H}}_{S}={\stackrel{^}{H}}_{B}\ensuremath{\bigotimes}I+{\stackrel{^}{x}}_{1}\ensuremath{\bigotimes}{\ensuremath{\sigma}}_{1}+{\stackrel{^}{x}}_{2}\ensuremath{\bigotimes}{\ensuremath{\sigma}}_{3}$ of the bosonic Hamiltonian ${\stackrel{^}{H}}_{B}={\stackrel{^}{p}}_{1}^{2}+{\stackrel{^}{p}}_{2}^{2}+{\stackrel{^}{x}}_{1}^{2}{\stackrel{^}{x}}_{2}^{2}$. For a long time, this bosonic Hamiltonian was considered one of the simplest systems which exhibit dynamical chaos both classically and quantum-mechanically. Its structure closely resembles that of spatially compactified pure Yang-Mills theory. Correspondingly, the structure of our supersymmetric Hamiltonian is similar to that of spatially compactified supersymmetric Yang-Mills theory, also known as the Banks-Fischler-Shenker-Susskind (BFSS) model. We present numerical evidence that a continuous energy spectrum of the supersymmetric model leads to monotonous growth of OTOCs down to the lowest temperatures, a property that is also expected for the BFSS model from holographic duality. We find that this growth is saturated by low-energy eigenstates with effectively one-dimensional wave functions and a completely nonchaotic energy level distribution. We observe a sharp boundary separating these low-energy states from the bulk of chaotic high-energy states. Our data suggest, although with limited confidence, that at low temperatures the OTOC growth might be exponential over a finite range of time, with the corresponding Lyapunov exponent scaling linearly with temperature. In contrast, the gapped low-energy spectrum of the bosonic Hamiltonian leads to oscillating OTOCs at low temperatures without any signatures of exponential growth. We also find that the OTOCs for the bosonic Hamiltonian are never sufficiently close to the classical Lyapunov distance. On the other hand, the OTOCs for the supersymmetric system agree with the classical limit reasonably well over a finite range of temperatures and evolution times.

arXiv:2205.09704•Physical review. D/Physical review. D.•Aug 1, 2022•14 citations

Universal spectral correlations in the chaotic wave function and the development of quantum chaos

Xiao Chen, Andreas W. W. Ludwig

We investigate the appearance of quantum chaos in a single many-body wave function by analyzing the statistical properties of the eigenvalues of its reduced density matrix ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{A}$ of a spatial subsystem $A$. We find that (i): the spectrum of the density matrix is described by so-called Wishart random matrix theory, which (ii): exhibits besides level repulsion, spectral rigidity, and universal spectral correlations between eigenvalues separated by distances ranging from one up to many mean level spacings, which we investigate. We use these universal spectral characteristics of the reduced density matrix as a definition of chaos in the wave function. A simple and precise characterization of such universal correlations in a spectrum is a segment of strictly linear growth at sufficiently long times, recently called the ``ramp,'' of the spectral form factor which is the Fourier transform of the correlation function between a pair of eigenvalues. It turns out that Wishart and standard random matrix theory have the same universal ``ramp.'' Specifically, here numerical results for the spectral form factor of the density matrix of generic nonintegrable many-body systems, such as one-dimensional quantum Ising and Floquet spin models, are found to exhibit a universal ``ramp'' identical to that appearing for a ``random pure state'' (``Page state,'' or ``Haar state''). The density matrix of the latter is precisely the Wishart random matrix, the reduced density matrix of a completely random wave function. In addition, we study the development of chaos in the wave function by letting an initial direct product state evolve under the unitary time evolution. We find that the universal spectral correlations as manifested by the ``ramp'' set in as soon as the entanglement entropy begins to grow, and first develop for the eigenvalues at the top of the spectrum of the density matrix ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{A}$, subsequently spreading over the entire spectrum at later times. Finally, we study a prethermalized regime described by a generalized Gibbs ensemble, which develops in a rapidly driven Floquet model at intermediate times. We find that the prethermalized regime exhibits no chaos, as evidenced by the absence of a ``ramp'' in the spectral form factor of the density matrix, while the universal spectral correlations start to develop when the prethermalized regime finally relaxes at late times to the fully thermalized (infinite temperature) chaotic regime.

arXiv:1710.02686•Physical review. B./Physical review. B•Aug 30, 2018•58 citations

Quantum chaos, thermalization, and entanglement generation in real-time simulations of the Banks-Fischler-Shenker-Susskind matrix model

P. V. Buividovich, Masanori Hanada, Andreas Schäfer A

We study numerically the onset of chaos and thermalization in the Banks-Fischler-Shenker-Susskind (BFSS) matrix model with and without fermions, considering Lyapunov exponents, entanglement generation, and quasinormal ringing.We approximate the real-time dynamics in terms of the most general Gaussian density matrices with parameters which obey self-consistent equations of motion, thus extending the applicability of real-time simulations beyond the classical limit.Initial values of these Gaussian density matrices are optimized to be as close as possible to the thermal equilibrium state of the system.Thus attempting to bridge between the low-energy regime with a calculable holographic description and the classical regime at high energies, we find that quantum corrections to classical dynamics tend to decrease the Lyapunov exponents, which is essential for consistency with the Maldacena-Shenker-Stanford bound at low temperatures.The entanglement entropy is found to exhibit an expected "scrambling" behavior-rapid initial growth followed by saturation.At least at high temperatures the entanglement saturation time appears to be governed by classical Lyapunov exponents.Decay of quasinormal modes is found to be characterized by the shortest timescale of all.We also find that while the bosonic matrix model becomes nonchaotic in the low-temperature regime, for the full BFSS model with fermions the leading Lyapunov exponent, entanglement saturation time, and decay rate of quasinormal modes all remain finite and nonzero down to the lowest temperatures.

arXiv:1810.03378•Physical review. D/Physical review. D.•Feb 19, 2019•29 citations

Chaotic Many-Body Quantum Dynamics, Spectral Correlations, and Energy Diffusion

J. T. Chalker, Dominik Hahn

We study chaotic many-body quantum dynamics in a minimal model with spatial structure and local interactions. It has a time-independent Hamiltonian, in contrast to quantum circuits and Brownian models, and is simple at the single-site level, in contrast to Sachdev-Ye-Kitaev chains. It is analytically tractable for large local Hilbert space dimension and weak intersite coupling. In this limit, we show that energy dynamics is described by a classical master equation and is diffusive. We also show that the spectral form factor can be expressed exactly in terms of the solution to this master equation. For a two-site system, we obtain closed-form expressions for both the two-point correlator of energy density and the spectral form factor, in essentially perfect agreement with numerical simulations. For an L -site system we show at late times how a linear ramp emerges in the spectral form factor, as universally expected from level repulsion in chaotic quantum systems. Conversely, at earlier times we identify two distinct mechanisms for an increase of the spectral form factor above its ramp value. One of these is associated with energy diffusion and is effective until the Thouless time, which varies as L 2 . The other involves contributions like those that would appear if the system were composed of many uncoupled subsystems: They generate a large enhancement of the spectral form factor and are suppressed on a timescale varying as ( ln L ) 2 . Besides being exact for the limit considered, we believe our approach provides the natural approximation even for small local Hilbert space dimension and strong intersite coupling. We present a numerical study of a spin-half chain, finding an early-time enhancement of the spectral form factor which is qualitatively similar to that in our solvable model.

arXiv:2510.02198•Physical Review X•Aug 25, 2026•1 citation

Quantum Chaos and Eigenstate Thermalization in Unified Matrix--QCA Universe\\ \large Spectral Rigidity, Random Matrix Universality and Entanglement Growth under Unified Time Scale

Haobo Ma, Wenlin Zhang

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.

•Nov 19, 2025•0 citations

Quantum chaos at finite temperature in local spin Hamiltonians

Christopher M. Langlett, Cheryne Jonay, Vedika Khemani +1 more

Understanding the emergence of chaos in many-body quantum systems away from semi-classical limits, particularly in spatially local interacting spin Hamiltonians, has been a long-standing problem. In these intrinsically quantum regimes, quantum chaos has been primarily understood through the correspondence between the eigensystem statistics of midspectrum eigenstates and the universal statistics described by random matrix theory (RMT). However, this correspondence no longer holds for finite-temperature eigenstates. Here we show that the statistical properties of finite-temperature eigenstates of quantum chaotic Hamiltonians can be accurately described by pure random states constrained by a local charge, with the average charge density of the constrained random state ensemble playing the same role as the average energy density of the eigenstates. By properly normalizing the energy density using a single Hamiltonian-dependent parameter that quantifies the typical energy per degree of freedom, we find excellent agreement between the entanglement entropy statistics of eigenstates and that of constrained random states. Interestingly, in small pockets of Hamiltonian parameter phase space which we previously identified as `maximally chaotic' [PRX 14, 031014 (2024)], we find excellent agreement not only at the level of the first moment, including O(1) corrections, but also at the level of statistical fluctuations. These results show that notions of maximal chaos -- in terms of how much randomness eigenstates contain -- can still be defined at finite temperature in physical Hamiltonian models away from semi-classical and large-$N$ limits.

arXiv:2501.13164•Jan 22, 2025•0 citations
We analyze in detail a sharp transition between the low-energy, low-dimensional eigenstates and the high-energy chaotic bulk of the spectrum for a simple supersymmetric quantum-mechanical model with Hamiltonian H^S=(p^12+p^22+x^12 x^22)⊗I+x^1⊗σ1+x^2⊗σ3\hat{H}_S = \left(\hat{p}_1^2 + \hat{p}_2^2 + \hat{x}_1^2 \, \hat{x}_2^2\right) \otimes I + \hat{x}_1 \otimes \sigma_1 + \hat{x}_2 \otimes \sigma_3H^S​=(p^​12​+p^​22​+x^12​x^22​)⊗I+x^1​⊗σ1​+x^2​⊗σ3​, which mimics the structure of the Banks-Fischler-Susskind-Stanford (BFSS) matrix model, the spatially compactified N=1\mathcal{N} = 1N=1 super-Yang-Mills theory. We conjecture that this transition might be similar to the transition between the D0D0D0-brane and MMM-theory regimes in the BFSS model, and find that it does not lead to irregularities in the thermodynamic equation of state. We demonstrate that real-time spectral form-factor for our supersymmetric model exhibits the ``ramp'' behavior typical for quantum chaos. We also analyze the entanglement entropy and the spectrum of the reduced density matrix of the eigenstates of H^S\hat{H}_SH^S​, considering one of the bosonic degrees of freedom as a subsystem. The entanglement entropy of low-energy eigenstates appears to be practically energy-independent. Exactly at the onset of random-matrix-type level spacing fluctuations, this behavior rapidly changes into a steady growth of entanglement with energy. We demonstrate that the spectrum of the reduced density matrix also exhibits universal level-spacing fluctuations towards its higher end, even for the ground state of the supersymmetric model. Thus even the regularly spaced, non-chaotic eigenstates contain some information about semi-classical chaotic dynamics at high energies.
We use exact diagonalization to study energy level statistics and out-of-time-order correlators (OTOCs) for the simplest supersymmetric extension ${\stackrel{^}{H}}_{S}={\stackrel{^}{H}}_{B}\ensuremath{\bigotimes}I+{\stackrel{^}{x}}_{1}\ensuremath{\bigotimes}{\ensuremath{\sigma}}_{1}+{\stackrel{^}{x}}_{2}\ensuremath{\bigotimes}{\ensuremath{\sigma}}_{3}$ of the bosonic Hamiltonian ${\stackrel{^}{H}}_{B}={\stackrel{^}{p}}_{1}^{2}+{\stackrel{^}{p}}_{2}^{2}+{\stackrel{^}{x}}_{1}^{2}{\stackrel{^}{x}}_{2}^{2}$. For a long time, this bosonic Hamiltonian was considered one of the simplest systems which exhibit dynamical chaos both classically and quantum-mechanically. Its structure closely resembles that of spatially compactified pure Yang-Mills theory. Correspondingly, the structure of our supersymmetric Hamiltonian is similar to that of spatially compactified supersymmetric Yang-Mills theory, also known as the Banks-Fischler-Shenker-Susskind (BFSS) model. We present numerical evidence that a continuous energy spectrum of the supersymmetric model leads to monotonous growth of OTOCs down to the lowest temperatures, a property that is also expected for the BFSS model from holographic duality. We find that this growth is saturated by low-energy eigenstates with effectively one-dimensional wave functions and a completely nonchaotic energy level distribution. We observe a sharp boundary separating these low-energy states from the bulk of chaotic high-energy states. Our data suggest, although with limited confidence, that at low temperatures the OTOC growth might be exponential over a finite range of time, with the corresponding Lyapunov exponent scaling linearly with temperature. In contrast, the gapped low-energy spectrum of the bosonic Hamiltonian leads to oscillating OTOCs at low temperatures without any signatures of exponential growth. We also find that the OTOCs for the bosonic Hamiltonian are never sufficiently close to the classical Lyapunov distance. On the other hand, the OTOCs for the supersymmetric system agree with the classical limit reasonably well over a finite range of temperatures and evolution times.
We investigate the appearance of quantum chaos in a single many-body wave function by analyzing the statistical properties of the eigenvalues of its reduced density matrix ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{A}$ of a spatial subsystem AAA. We find that (i): the spectrum of the density matrix is described by so-called Wishart random matrix theory, which (ii): exhibits besides level repulsion, spectral rigidity, and universal spectral correlations between eigenvalues separated by distances ranging from one up to many mean level spacings, which we investigate. We use these universal spectral characteristics of the reduced density matrix as a definition of chaos in the wave function. A simple and precise characterization of such universal correlations in a spectrum is a segment of strictly linear growth at sufficiently long times, recently called the ``ramp,'' of the spectral form factor which is the Fourier transform of the correlation function between a pair of eigenvalues. It turns out that Wishart and standard random matrix theory have the same universal ``ramp.'' Specifically, here numerical results for the spectral form factor of the density matrix of generic nonintegrable many-body systems, such as one-dimensional quantum Ising and Floquet spin models, are found to exhibit a universal ``ramp'' identical to that appearing for a ``random pure state'' (``Page state,'' or ``Haar state''). The density matrix of the latter is precisely the Wishart random matrix, the reduced density matrix of a completely random wave function. In addition, we study the development of chaos in the wave function by letting an initial direct product state evolve under the unitary time evolution. We find that the universal spectral correlations as manifested by the ``ramp'' set in as soon as the entanglement entropy begins to grow, and first develop for the eigenvalues at the top of the spectrum of the density matrix ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\rho}}}_{A}$, subsequently spreading over the entire spectrum at later times. Finally, we study a prethermalized regime described by a generalized Gibbs ensemble, which develops in a rapidly driven Floquet model at intermediate times. We find that the prethermalized regime exhibits no chaos, as evidenced by the absence of a ``ramp'' in the spectral form factor of the density matrix, while the universal spectral correlations start to develop when the prethermalized regime finally relaxes at late times to the fully thermalized (infinite temperature) chaotic regime.
We study numerically the onset of chaos and thermalization in the Banks-Fischler-Shenker-Susskind (BFSS) matrix model with and without fermions, considering Lyapunov exponents, entanglement generation, and quasinormal ringing.We approximate the real-time dynamics in terms of the most general Gaussian density matrices with parameters which obey self-consistent equations of motion, thus extending the applicability of real-time simulations beyond the classical limit.Initial values of these Gaussian density matrices are optimized to be as close as possible to the thermal equilibrium state of the system.Thus attempting to bridge between the low-energy regime with a calculable holographic description and the classical regime at high energies, we find that quantum corrections to classical dynamics tend to decrease the Lyapunov exponents, which is essential for consistency with the Maldacena-Shenker-Stanford bound at low temperatures.The entanglement entropy is found to exhibit an expected "scrambling" behavior-rapid initial growth followed by saturation.At least at high temperatures the entanglement saturation time appears to be governed by classical Lyapunov exponents.Decay of quasinormal modes is found to be characterized by the shortest timescale of all.We also find that while the bosonic matrix model becomes nonchaotic in the low-temperature regime, for the full BFSS model with fermions the leading Lyapunov exponent, entanglement saturation time, and decay rate of quasinormal modes all remain finite and nonzero down to the lowest temperatures.
We study chaotic many-body quantum dynamics in a minimal model with spatial structure and local interactions. It has a time-independent Hamiltonian, in contrast to quantum circuits and Brownian models, and is simple at the single-site level, in contrast to Sachdev-Ye-Kitaev chains. It is analytically tractable for large local Hilbert space dimension and weak intersite coupling. In this limit, we show that energy dynamics is described by a classical master equation and is diffusive. We also show that the spectral form factor can be expressed exactly in terms of the solution to this master equation. For a two-site system, we obtain closed-form expressions for both the two-point correlator of energy density and the spectral form factor, in essentially perfect agreement with numerical simulations. For an L -site system we show at late times how a linear ramp emerges in the spectral form factor, as universally expected from level repulsion in chaotic quantum systems. Conversely, at earlier times we identify two distinct mechanisms for an increase of the spectral form factor above its ramp value. One of these is associated with energy diffusion and is effective until the Thouless time, which varies as L 2 . The other involves contributions like those that would appear if the system were composed of many uncoupled subsystems: They generate a large enhancement of the spectral form factor and are suppressed on a timescale varying as ( ln L ) 2 . Besides being exact for the limit considered, we believe our approach provides the natural approximation even for small local Hilbert space dimension and strong intersite coupling. We present a numerical study of a spin-half chain, finding an early-time enhancement of the spectral form factor which is qualitatively similar to that in our solvable model.
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.
Understanding the emergence of chaos in many-body quantum systems away from semi-classical limits, particularly in spatially local interacting spin Hamiltonians, has been a long-standing problem. In these intrinsically quantum regimes, quantum chaos has been primarily understood through the correspondence between the eigensystem statistics of midspectrum eigenstates and the universal statistics described by random matrix theory (RMT). However, this correspondence no longer holds for finite-temperature eigenstates. Here we show that the statistical properties of finite-temperature eigenstates of quantum chaotic Hamiltonians can be accurately described by pure random states constrained by a local charge, with the average charge density of the constrained random state ensemble playing the same role as the average energy density of the eigenstates. By properly normalizing the energy density using a single Hamiltonian-dependent parameter that quantifies the typical energy per degree of freedom, we find excellent agreement between the entanglement entropy statistics of eigenstates and that of constrained random states. Interestingly, in small pockets of Hamiltonian parameter phase space which we previously identified as `maximally chaotic' [PRX 14, 031014 (2024)], we find excellent agreement not only at the level of the first moment, including O(1) corrections, but also at the level of statistical fluctuations. These results show that notions of maximal chaos -- in terms of how much randomness eigenstates contain -- can still be defined at finite temperature in physical Hamiltonian models away from semi-classical and large-NNN limits.