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:2007.13837 · Physical review. D/Physical review. D.

Sparse Sachdev-Ye-Kitaev model, quantum chaos, and gravity duals

Antonio M. Garcı́a-Garcı́a, Yiyang Jia, Dario Rosa, J. J. M. Verbaarschot

Published May 3, 2021• 69 citations
arXivHTMLPDF

Abstract

We study a sparse Sachdev-Ye-Kitaev (SYK) model with $N$ Majoranas where only $\ensuremath{\sim}kN$ independent matrix elements are nonzero. We identify a minimum $k\ensuremath{\gtrsim}1$ for quantum chaos to occur by a level statistics analysis. The spectral density in this region, and for a larger $k$, is still given by the Schwarzian prediction of the dense SYK model, though with renormalized parameters. Similar results are obtained for a beyond linear scaling with $N$ of the number of nonzero matrix elements. This is a strong indication that this is the minimum connectivity for the sparse SYK model to still have a quantum gravity dual. We also find an intriguing exact relation between the leading correction to moments of the spectral density due to sparsity and the leading $1/d$ correction of Parisi's U(1) lattice gauge theory in a $d$-dimensional hypercube. In the $k\ensuremath{\rightarrow}1$ limit, different disorder realizations of the sparse SYK model show emergent random matrix statistics that for fixed $N$ can be in any universality class of the tenfold way. The agreement with random matrix statistics is restricted to short-range correlations, no more than a few level spacings, in particular in the tail of the spectrum. In addition, emergent discrete global symmetries in most of the disorder realizations for $k$ slightly below one give rise to ${2}^{m}$-fold degenerate spectra, with $m$ being a positive integer. For $k=3/4$, we observe a large number of such emergent global symmetries with a maximum ${2}^{8}$-fold degenerate spectra for $N=26$.

DOI
10.1103/physrevd.103.106002

Related papers

Q -Laguerre spectral density and quantum chaos in the Wishart-Sachdev-Ye-Kitaev model

Lucas Sá, Antonio M. Garcı́a-Garcı́a

We study the Wishart-Sachdev-Ye-Kitaev (WSYK) model consisting of two $\stackrel{^}{q}$-body Sachdev-Ye-Kitaev (SYK) models with general complex couplings, one the Hermitian conjugate of the other, living in off-diagonal blocks of a larger WSYK Hamiltonian. The spectrum is positive with a hard edge at zero energy. We employ diagrammatic and combinatorial techniques to compute analytically the low-order moments of the Hamiltonian. In the limit of large number $N$ of Majoranas, we have found striking similarities with the moments of the weight function of the Al-Salam-Chihara $Q$-Laguerre polynomials. For $\stackrel{^}{q}=3$, 4, the $Q$-Laguerre prediction, with $Q=Q(\stackrel{^}{q},N)$ also computed analytically, agrees well with exact diagonalization results for $30<N\ensuremath{\le}34$ while we observe some deviations for $\stackrel{^}{q}=2$. The most salient feature of the spectral density is that, for odd $\stackrel{^}{q}$, low-energy excitations grow as a stretched exponential, with a functional form different from that of the supersymmetric SYK model. For $\stackrel{^}{q}=4$, a detailed analysis of level statistics reveals quantum chaotic dynamics even for time scales substantially shorter than the Heisenberg time. More specifically, the spacing ratios in the bulk of the spectrum and the microscopic spectral density and the number variance close to the hard edge are very well approximated by that of an ensemble of random matrices that, depending on $N$, belong to the chiral or superconducting universality classes. In particular, we report the first realization of level statistics belonging to the chiral Gaussian Unitary Ensemble universality class, which completes the tenfold-way classification in the SYK model.

arXiv:2104.07647•Physical review. D/Physical review. D.•Jan 5, 2022•12 citations

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

Antonio M. Garcı́a-Garcı́a, Chang Liu, J. J. M. Verbaarschot

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body (q>2) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N≤64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N=64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature.

arXiv:2311.00639•Physical Review Letters•Aug 26, 2024•11 citations

The Wishart-Sachdev-Ye-Kitaev model: Q-Laguerre spectral density and quantum chaos.

Lucas Sá, Antonio M. Garcı́a-Garcı́a

We study the Wishart-Sachdev-Ye-Kitaev (WSYK) model consisting of two $\hat{q}$-body Sachdev-Ye-Kitaev (SYK) models with general complex couplings, one the complex conjugate of the other, living in off-diagonal blocks of a larger WSYK Hamiltonian. In the limit of large number $N$ of Majoranas and $\hat{q} \propto \sqrt N$, we employ diagrammatic and combinatorial techniques to show that the spectral density of the model can be approximated by the weight function of the Al-Salam-Chihara Q-Laguerre polynomials. For $\hat{q} = 2, 3, 4$, the analytical spectral density, with $Q=Q(\hat{q},N)$ computed analytically, agrees well with numerical results obtained by exact diagonalization. The spectrum is positive with a hard edge at zero energy and, for odd $\hat{q}$, low-energy excitations grow as a stretched exponential, though the functional form is different from that of the supersymmetric SYK model. The microscopic spectral density and level statistics close to the hard edge are very well approximated by that of an ensemble of random matrices belonging to non-standard universality classes, which indicates quantum chaotic dynamics.

•Apr 15, 2021•1 citation

Binary-coupling sparse Sachdev-Ye-Kitaev model: An improved model of quantum chaos and holography

Masaki Tezuka, Onur Oktay, Enrico Rinaldi +2 more

The sparse version of the Sachdev-Ye-Kitaev (SYK) model reproduces essential features of the original SYK model while reducing the number of disorder parameters. In this Research Letter, we propose a further simplification of the model which we call the binary-coupling sparse SYK model. We set the nonzero couplings to be $\ifmmode\pm\else\textpm\fi{}1$, rather than being sampled from a continuous distribution such as Gaussian. Remarkably, this simplification turns out to be an improvement: The binary-coupling model exhibits strong correlations in the spectrum, which is the important feature of the original SYK model that leads to the quick onset of the random-matrix universality, more efficiently in terms of the number of nonzero terms. This model is better suited for analog or digital quantum simulations of quantum chaotic behavior and holographic metals due to its simplicity and scaling properties.

arXiv:2208.12098•Physical review. B./Physical review. B•Feb 6, 2023•27 citations

Universality and Thouless energy in the supersymmetric Sachdev-Ye-Kitaev model

Antonio M. Garcı́a-Garcı́a, Yiyang Jia, J. J. M. Verbaarschot

We investigate the supersymmetric Sachdev-Ye-Kitaev (SYK) model, $N$ Majorana fermions with infinite range interactions in $0+1$ dimensions. We have found that, close to the ground state $E\ensuremath{\approx}0$, discrete symmetries alter qualitatively the spectral properties with respect to the non-supersymmetric SYK model. The average spectral density at finite $N$, which we compute analytically and numerically, grows exponentially with $N$ for $E\ensuremath{\approx}0$. However the chiral condensate, which is normalized with respect the total number of eigenvalues, vanishes in the thermodynamic limit. Slightly above $E\ensuremath{\approx}0$, the spectral density grows exponentially with the energy. Deep in the quantum regime, corresponding to the first $O(N)$ eigenvalues, the average spectral density is universal and well described by random matrix ensembles with chiral and superconducting discrete symmetries. The dynamics for $E\ensuremath{\approx}0$ is investigated by level fluctuations. Also in this case we find excellent agreement with the prediction of chiral and superconducting random matrix ensembles for eigenvalue separations smaller than the Thouless energy, which seems to scale linearly with $N$. Deviations beyond the Thouless energy, which describes how ergodicity is approached, are universally characterized by a quadratic growth of the number variance. In the time domain, we have found analytically that the spectral form factor $g(t)$, obtained from the connected two-level correlation function of the unfolded spectrum, decays as $1/{t}^{2}$ for times shorter but comparable to the Thouless time with $g(0)$ related to the coefficient of the quadratic growth of the number variance. Our results provide further support that quantum black holes are ergodic and therefore can be classified by random matrix theory.

arXiv:1801.01071•Physical review. D/Physical review. D.•May 4, 2018•62 citations

Quantum chaos on edge

Alexander Altland, Kun Woo Kim, Tobias Micklitz +3 more

Recently, the physics of many-body quantum chaotic systems close to their ground states has come under intensified scrutiny. Such studies are motivated by the emergence of model systems exhibiting chaotic fluctuations throughout the entire spectrum [the Sachdev-Ye-Kitaev (SYK) model being a renowned representative] as well as by the physics of holographic principles, which likewise unfold close to ground states. Interpreting the edge of the spectrum as a quantum critical point, here we combine a wide range of analytical and numerical methods to the identification and comprehensive description of two different universality classes: the near edge physics of “sparse” and the near edge of “dense” chaotic systems. The distinction lies in the ratio between the number of a system's random parameters and its Hilbert space dimension, which is exponentially small or algebraically small in the sparse and dense case, respectively. Notable representatives of the two classes are generic chaotic many-body models (sparse) and invariant random matrix ensembles or chaotic gravitational systems (dense). While the two families share identical spectral correlations at energy scales comparable to the level spacing, the density of states and its fluctuations near the edge are different. Considering the SYK model as a representative of the sparse class, we apply a combination of field theory and exact diagonalization to a detailed discussion of its edge spectrum. Conversely, Jackiw-Teitelboim gravity is our reference model for the dense class, where an analysis of the gravitational path integral and random matrix theory reveal universal differences to the sparse class, whose implications for the construction of holographic principles we discuss. Published by the American Physical Society 2024

arXiv:2403.13516•Physical Review Research•Sep 12, 2024•14 citations
We study a sparse Sachdev-Ye-Kitaev (SYK) model with NNN Majoranas where only \ensuremath∼kN\ensuremath{\sim}kN\ensuremath∼kN independent matrix elements are nonzero. We identify a minimum k\ensuremath≳1k\ensuremath{\gtrsim}1k\ensuremath≳1 for quantum chaos to occur by a level statistics analysis. The spectral density in this region, and for a larger kkk, is still given by the Schwarzian prediction of the dense SYK model, though with renormalized parameters. Similar results are obtained for a beyond linear scaling with NNN of the number of nonzero matrix elements. This is a strong indication that this is the minimum connectivity for the sparse SYK model to still have a quantum gravity dual. We also find an intriguing exact relation between the leading correction to moments of the spectral density due to sparsity and the leading 1/d1/d1/d correction of Parisi's U(1) lattice gauge theory in a ddd-dimensional hypercube. In the k\ensuremath→1k\ensuremath{\rightarrow}1k\ensuremath→1 limit, different disorder realizations of the sparse SYK model show emergent random matrix statistics that for fixed NNN can be in any universality class of the tenfold way. The agreement with random matrix statistics is restricted to short-range correlations, no more than a few level spacings, in particular in the tail of the spectrum. In addition, emergent discrete global symmetries in most of the disorder realizations for kkk slightly below one give rise to 2m{2}^{m}2m-fold degenerate spectra, with mmm being a positive integer. For k=3/4k=3/4k=3/4, we observe a large number of such emergent global symmetries with a maximum 28{2}^{8}28-fold degenerate spectra for N=26N=26N=26.
We study the Wishart-Sachdev-Ye-Kitaev (WSYK) model consisting of two $\stackrel{^}{q}$-body Sachdev-Ye-Kitaev (SYK) models with general complex couplings, one the Hermitian conjugate of the other, living in off-diagonal blocks of a larger WSYK Hamiltonian. The spectrum is positive with a hard edge at zero energy. We employ diagrammatic and combinatorial techniques to compute analytically the low-order moments of the Hamiltonian. In the limit of large number NNN of Majoranas, we have found striking similarities with the moments of the weight function of the Al-Salam-Chihara QQQ-Laguerre polynomials. For $\stackrel{^}{q}=3$, 4, the QQQ-Laguerre prediction, with $Q=Q(\stackrel{^}{q},N)$ also computed analytically, agrees well with exact diagonalization results for $30<N\ensuremath{\le}34$ while we observe some deviations for $\stackrel{^}{q}=2$. The most salient feature of the spectral density is that, for odd $\stackrel{^}{q}$, low-energy excitations grow as a stretched exponential, with a functional form different from that of the supersymmetric SYK model. For $\stackrel{^}{q}=4$, a detailed analysis of level statistics reveals quantum chaotic dynamics even for time scales substantially shorter than the Heisenberg time. More specifically, the spacing ratios in the bulk of the spectrum and the microscopic spectral density and the number variance close to the hard edge are very well approximated by that of an ensemble of random matrices that, depending on NNN, belong to the chiral or superconducting universality classes. In particular, we report the first realization of level statistics belonging to the chiral Gaussian Unitary Ensemble universality class, which completes the tenfold-way classification in the SYK model.
The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body (q>2) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N≤64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N=64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature.
We study the Wishart-Sachdev-Ye-Kitaev (WSYK) model consisting of two q^\hat{q}q^​-body Sachdev-Ye-Kitaev (SYK) models with general complex couplings, one the complex conjugate of the other, living in off-diagonal blocks of a larger WSYK Hamiltonian. In the limit of large number NNN of Majoranas and q^∝N\hat{q} \propto \sqrt Nq^​∝N​, we employ diagrammatic and combinatorial techniques to show that the spectral density of the model can be approximated by the weight function of the Al-Salam-Chihara Q-Laguerre polynomials. For q^=2,3,4\hat{q} = 2, 3, 4q^​=2,3,4, the analytical spectral density, with Q=Q(q^,N)Q=Q(\hat{q},N)Q=Q(q^​,N) computed analytically, agrees well with numerical results obtained by exact diagonalization. The spectrum is positive with a hard edge at zero energy and, for odd q^\hat{q}q^​, low-energy excitations grow as a stretched exponential, though the functional form is different from that of the supersymmetric SYK model. The microscopic spectral density and level statistics close to the hard edge are very well approximated by that of an ensemble of random matrices belonging to non-standard universality classes, which indicates quantum chaotic dynamics.
The sparse version of the Sachdev-Ye-Kitaev (SYK) model reproduces essential features of the original SYK model while reducing the number of disorder parameters. In this Research Letter, we propose a further simplification of the model which we call the binary-coupling sparse SYK model. We set the nonzero couplings to be $\ifmmode\pm\else\textpm\fi{}1$, rather than being sampled from a continuous distribution such as Gaussian. Remarkably, this simplification turns out to be an improvement: The binary-coupling model exhibits strong correlations in the spectrum, which is the important feature of the original SYK model that leads to the quick onset of the random-matrix universality, more efficiently in terms of the number of nonzero terms. This model is better suited for analog or digital quantum simulations of quantum chaotic behavior and holographic metals due to its simplicity and scaling properties.
We investigate the supersymmetric Sachdev-Ye-Kitaev (SYK) model, NNN Majorana fermions with infinite range interactions in 0+10+10+1 dimensions. We have found that, close to the ground state $E\ensuremath{\approx}0$, discrete symmetries alter qualitatively the spectral properties with respect to the non-supersymmetric SYK model. The average spectral density at finite NNN, which we compute analytically and numerically, grows exponentially with NNN for $E\ensuremath{\approx}0$. However the chiral condensate, which is normalized with respect the total number of eigenvalues, vanishes in the thermodynamic limit. Slightly above $E\ensuremath{\approx}0$, the spectral density grows exponentially with the energy. Deep in the quantum regime, corresponding to the first O(N)O(N)O(N) eigenvalues, the average spectral density is universal and well described by random matrix ensembles with chiral and superconducting discrete symmetries. The dynamics for $E\ensuremath{\approx}0$ is investigated by level fluctuations. Also in this case we find excellent agreement with the prediction of chiral and superconducting random matrix ensembles for eigenvalue separations smaller than the Thouless energy, which seems to scale linearly with NNN. Deviations beyond the Thouless energy, which describes how ergodicity is approached, are universally characterized by a quadratic growth of the number variance. In the time domain, we have found analytically that the spectral form factor g(t)g(t)g(t), obtained from the connected two-level correlation function of the unfolded spectrum, decays as 1/t21/{t}^{2}1/t2 for times shorter but comparable to the Thouless time with g(0)g(0)g(0) related to the coefficient of the quadratic growth of the number variance. Our results provide further support that quantum black holes are ergodic and therefore can be classified by random matrix theory.
Recently, the physics of many-body quantum chaotic systems close to their ground states has come under intensified scrutiny. Such studies are motivated by the emergence of model systems exhibiting chaotic fluctuations throughout the entire spectrum [the Sachdev-Ye-Kitaev (SYK) model being a renowned representative] as well as by the physics of holographic principles, which likewise unfold close to ground states. Interpreting the edge of the spectrum as a quantum critical point, here we combine a wide range of analytical and numerical methods to the identification and comprehensive description of two different universality classes: the near edge physics of “sparse” and the near edge of “dense” chaotic systems. The distinction lies in the ratio between the number of a system's random parameters and its Hilbert space dimension, which is exponentially small or algebraically small in the sparse and dense case, respectively. Notable representatives of the two classes are generic chaotic many-body models (sparse) and invariant random matrix ensembles or chaotic gravitational systems (dense). While the two families share identical spectral correlations at energy scales comparable to the level spacing, the density of states and its fluctuations near the edge are different. Considering the SYK model as a representative of the sparse class, we apply a combination of field theory and exact diagonalization to a detailed discussion of its edge spectrum. Conversely, Jackiw-Teitelboim gravity is our reference model for the dense class, where an analysis of the gravitational path integral and random matrix theory reveal universal differences to the sparse class, whose implications for the construction of holographic principles we discuss. Published by the American Physical Society 2024