Abstract We analytically evaluate the moments of the spectral density of the q-body Sachdev-Ye-Kitaev (SYK) model, and obtain order 1/N 2 corrections for all moments, where N is the total number of Majorana fermions. To order 1/N, moments are given by those of the weight function of the Q-Hermite polynomials. Representing Wick contractions by rooted chord diagrams, we show that the 1/N 2 correction for each chord diagram is proportional to the number of triangular loops of the corresponding intersection graph, with an extra grading factor when q is odd. Therefore the problem of finding 1/N 2 corrections is mapped to a triangle counting problem. Since the total number of triangles is a purely graph-theoretic property, we can compute them for the q = 1 and q = 2 SYK models, where the exact moments can be obtained analytically using other methods, and therefore we have solved the moment problem for any q to 1/N 2 accuracy. The moments are then used to obtain the spectral density of the SYK model to order 1/N 2. We also obtain an exact analytical result for all contraction diagrams contributing to the moments, which can be evaluated up to eighth order. This shows that the Q-Hermite approximation is accurate even for small values of N.
Abstract We analytically evaluate the moments of the spectral density of the q-body Sachdev-Ye-Kitaev (SYK) model, and obtain order 1/N 2 corrections for all moments, where N is the total number of Majorana fermions. To order 1/N, moments are given by those of the weight function of the Q-Hermite polynomials. Representing Wick contractions by rooted chord diagrams, we show that the 1/N 2 correction for each chord diagram is proportional to the number of triangular loops of the corresponding intersection graph, with an extra grading factor when q is odd. Therefore the problem of finding 1/N 2 corrections is mapped to a triangle counting problem. Since the total number of triangles is a purely graph-theoretic property, we can compute them for the q = 1 and q = 2 SYK models, where the exact moments can be obtained analytically using other methods, and therefore we have solved the moment problem for any q to 1/N 2 accuracy. The moments are then used to obtain the spectral density of the SYK model to order 1/N 2. We also obtain an exact analytical result for all contraction diagrams contributing to the moments, which can be evaluated up to eighth order. This shows that the Q-Hermite approximation is accurate even for small values of N.
We analytically evaluate the moments of the spectral density of the q-body Sachdev-Ye-Kitaev (SYK) model, and obtain order 1/N2 corrections for all moments, where N is the total number of Majorana fermions. To order 1/N, moments are given by those of the weight function of the Q-Hermite polynomials. Representing Wick contractions by rooted chord diagrams, we show that the 1/N2 correction for each chord diagram is proportional to the number of triangular loops of the corresponding intersection graph, with an extra grading factor when q is odd. Therefore the problem of finding 1/N2 corrections is mapped to a triangle counting problem. Since the total number of triangles is a purely graph-theoretic property, we can compute them for the q = 1 and q = 2 SYK models, where the exact moments can be obtained analytically using other methods, and therefore we have solved the moment problem for any q to 1/N2 accuracy. The moments are then used to obtain the spectral density of the SYK model to order 1/N2. We also obtain an exact analytical result for all contraction diagrams contributing to the moments, which can be evaluated up to eighth order. This shows that the Q-Hermite approximation is accurate even for small values of N.
We derive an approximate analytical formula for the spectral density of the q-body Sachdev-Ye-Kitaev (SYK) model obtained by summing a class of diagrams representing leading intersecting contractions. This expression agrees with that of Q-Hermite polynomials, with Q a nontrivial function of $q\ensuremath{\ge}2$ and the number of Majorana fermions N. Numerical results, obtained by exact diagonalization, are in excellent agreement with this approximate analytical spectral density even for relatively small $N\ensuremath{\sim}8$. For $N\ensuremath{\gg}1$ and not close to the edge of the spectrum, we find that the approximate analytical spectral density simplifies to ${\ensuremath{\rho}}_{\mathrm{asym}}(E)=\mathrm{exp}[2{\mathrm{arcsin}}^{2}(E/{E}_{0})/\mathrm{log}\ensuremath{\eta}]$, where $\ensuremath{\eta}(N,q)$ is the suppression factor of the contribution of intersecting Wick contractions relative to nested contractions and E0 is the ground-state energy per particle. This spectral density reproduces the known result for the free energy in the large-q and large-N limit at arbitrary values of the temperature. In the infrared region, where the SYK model is believed to have a gravity dual, the analytical spectral density is given by $\ensuremath{\rho}(E)\ensuremath{\sim}\mathrm{sinh}[2\ensuremath{\pi}\sqrt{2}\sqrt{(1\ensuremath{-}E/{E}_{0})/(\ensuremath{-}\mathrm{log}\ensuremath{\eta})}]$. It therefore has a square-root edge, as in random matrix ensembles, followed by an exponential growth, a distinctive feature of black holes and also of low-energy nuclear excitations. Results for level statistics in this region confirm the agreement with random matrix theory. Physically this is a signature that, for sufficiently long times, the SYK model and its gravity dual evolve to a fully ergodic state whose dynamics only depends on the global symmetry of the system. Our results strongly suggest that random matrix correlations are a universal feature of quantum black holes and that the SYK model, combined with holography, may be relevant to modeling certain aspects of the nuclear dynamics.
Abstract We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations (of the order of the number of Majorana fermions N) from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being due to scale fluctuations for which we give a simple analytical estimate. Our numerical results for N = 32 show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion can be obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations in terms of a nonlinear σ-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlation functions to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the σ-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for q = 4 or higher even q-body interactions, but also for the q = 2 SYK model albeit with a much smaller Thouless energy while the correct result in this case should have been Poisson statistics. In our numerical studies we analyze the number variance and spectral form factor for N = 32 and q = 4. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings for the number variance or exponentially short times in the case of the spectral form factor.
Abstract In this paper we explain the relation between the free energy of the SYK model for N Majorana fermions with a random q-body interaction and the moments of its spectral density. The high temperature expansion of the free energy gives the cumulants of the spectral density. Using that the cumulants are extensive we find the p dependence of the 1/N 2 correction of the 2p-th moments obtained in [1]. Conversely, the 1/N 2 corrections to the moments give the correction (even q) to the β 6 coefficient of the high temperature expansion of the free energy for arbitrary q. Our result agrees with the 1/q 3 correction obtained by Tarnopolsky using a mean field expansion. These considerations also lead to a more powerful method for solving the moment problem and intersection-graph enumeration problems. We take advantage of this and push the moment calculation to 1/N 3 order and find surprisingly simple enumeration identities for intersection graphs. The 1/N 3 corrections to the moments, give corrections to the β 8 coefficient (for even q) of the high temperature expansion of the free energy which have not been calculated before. Results for odd q, where the SYK “Hamiltonian” is the supercharge of a supersymmetric theory are discussed as well.
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.
We present a diagrammatic expansion for thermal expectation values of low-weight Majorana strings in a fixed realization of the N-fermion Sachdev--Ye--Kitaev (SYK) model. At each order in the expansion, the observable is expressed as a finite sum of polynomials in the couplings multiplied by temperature-dependent kernels built from the large-N melonic propagator. At fixed βJ, the expansion is organized in powers of N−1/2, and each order can be efficiently evaluated with classical resources polynomial in N. The expansion also extends to correlations at arbitrary imaginary times and sufficiently short real times. For the four-body SYK model, we show that the diagrammatic predictions give an excellent account of exact diagonalization results at sizes up to N=24. We also apply the method to the Maldacena--Qi (MQ) model formed from two identical copies of a fixed SYK realization. This enables the identification of a nearby SYK-like Hamiltonian whose thermofield double state has higher fidelity with the Maldacena--Qi ground state than the thermofield double of the original SYK Hamiltonian.