aiqbench analysis
This study investigates the transition to chaotic dynamics in the Yukawa-Sachdev-Ye-Kitaev model, a theoretical framework relevant to quantum simulation and non-Fermi liquid physics. The research quantifies how spectral correlations evolve from Poissonian to random-matrix statistics (GUE/GOE) as system parameters cross a critical threshold, offering insights into the onset of chaos in complex quantum systems.
Onset of random-matrix statistics in the Yukawa-Sachdev-Ye-Kitaev model: Spectral correlations and Krylov complexity
arXiv:2609.30333v1 Announce Type: cross Abstract: The Yukawa-Sachdev-Ye-Kitaev (YSYK) model, a solvable model of non-Fermi liquids with a proposed cavity-QED realization, couples $N$ complex fermions to $M$ bosons of frequency $\omega_0$ through random Yukawa couplings of strength $g$, and its dynamics depends only on $R\equiv\omega_0/g^{2/3}$. For $R\ll1$ the bosons are slow and act as a static random background in which the fermions are free, while for $R\gg1$ they are only virtually excited and generate a random four-fermion interaction of rank $M$, which is generically chaotic for $M\geq2$. Using exact diagonalization at infinite temperature, we follow the emergence of random-matrix statistics between these limits with unfolded spectral statistics and the Krylov complexity of the thermofield-double state. For complex couplings the level-spacing distribution evolves from Poisson to Gaussian-unitary-ensemble (GUE) statistics, the spectral form factor develops a correlation hole and a ramp, and the Krylov complexity evaluated on the unfolded spectrum develops a peak above its late-time plateau whose height reaches the GUE value. All diagnostics locate the onset around $R\approx0.3$, and adding boson modes moves it to smaller $R$. This delineates the parameter regime in which a cavity-QED realization would display random-matrix spectral correlations. On the other hand, real couplings give the same crossover towards Gaussian-orthogonal-ensemble (GOE) statistics, and along an interpolation between the two the peak height follows the mean gap ratio, so the peak identifies the symmetry class as well as the onset of chaos. Without unfolding the peak stays about 20$\%$ below its random-matrix value in both classes, which we trace to the non-semicircular density of states, and at larger $R$ the boson-number bands slow down the spreading of the state in Krylov space.
Abstract:The Yukawa-Sachdev-Ye-Kitaev (YSYK) model, a solvable model of non-Fermi liquids with a proposed cavity-QED realization, couples $N$ complex fermions to $M$ bosons of frequency $\omega_0$ through random Yukawa couplings of strength $g$, and its dynamics depends only on $R\equiv\omega_0/g^{2/3}$. For $R\ll1$ the bosons are slow and act as a static random background in which the fermions are free, while for $R\gg1$ they are only virtually excited and generate a random four-fermion interaction of rank $M$, which is generically chaotic for $M\geq2$. Using exact diagonalization at infinite temperature, we follow the emergence of random-matrix statistics between these limits with unfolded spectral statistics and the Krylov complexity of the thermofield-double state. For complex couplings the level-spacing distribution evolves from Poisson to Gaussian-unitary-ensemble (GUE) statistics, the spectral form factor develops a correlation hole and a ramp, and the Krylov complexity evaluated on the unfolded spectrum develops a peak above its late-time plateau whose height reaches the GUE value. All diagnostics locate the onset around $R\approx0.3$, and adding boson modes moves it to smaller $R$. This delineates the parameter regime in which a cavity-QED realization would display random-matrix spectral correlations. On the other hand, real couplings give the same crossover towards Gaussian-orthogonal-ensemble (GOE) statistics, and along an interpolation between the two the peak height follows the mean gap ratio, so the peak identifies the symmetry class as well as the onset of chaos. Without unfolding the peak stays about 20$\%$ below its random-matrix value in both classes, which we trace to the non-semicircular density of states, and at larger $R$ the boson-number bands slow down the spreading of the state in Krylov space.
Chaotic Dynamics (nlin.CD); Strongly Correlated Electrons (cond-mat.str-el); High Energy Physics - Theory (hep-th); Quantum Physics (quant-ph)
arXiv:2609.30333 [nlin.CD]
(or arXiv:2609.30333v1 [nlin.CD] for this version)
https://doi.org/10.48550/arXiv.2609.30333 Focus to learn more arXiv-issued DOI via DataCite (pending registration)
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.
Full text under the publisher's open licence · 376 words
Related coverage
Quantum Chaos Diagnostics for non-Hermitian Systems from Bi-Lanczos Krylov Dynamics
arXiv · Sep 17, 2026

From Randomness to Order: Physicists Map a Chaos to Integrability Transition
Quantum Zeitgeist · May 22, 2026
Entanglement entropy of fermions in a strange metal
arXiv · Aug 6, 2026
Unitary $k$-designs without Hamiltonian quenches
arXiv · Jul 24, 2026