Investigations on the SYK Model and its Dual Gravity Theory

The Sachdev–Ye–Kitaev (SYK) model is a prominent model of strongly interacting fermions that serves as a toy model of quantum gravity and black hole physics. In this work, we study the Trotter error and gate complexity of the quantum simulation of the SYK model using Lie–Trotter–Suzuki formulas. Building on recent results by Chen and Brandão \cite{Anthony} — in particular their uniform smoothing technique for random matrix polynomials — we derive bounds on the first- and higher-order Trotter error of the SYK model, and subsequently find near-optimal gate complexities for simulating these models using Lie–Trotter–Suzuki formulas. For the k -local SYK model on n Majorana fermions, at time t , our gate complexity estimates for the first-order Lie–Trotter–Suzuki formula scales with O ~ ( n k + 5 2 t 2 ) for even k and O ~ ( n k + 3 t 2 ) for odd k , and the gate complexity of simulations using higher-order formulas scales with O ~ ( n k + 1 2 t ) for even k and O ~ ( n k + 1 t ) for odd k . Given that the SYK model has Θ ( n k ) terms, these estimates are close to optimal. These gate complexities can be further improved upon in the context of simulating the time evolution of an arbitrary fixed input state | ψ ⟩ , leading to a O ( n 2 ) -reduction in gate complexity for first-order formulas and O ( n ) -reduction for higher-order formulas.We also apply our techniques to the sparse SYK model, which is a simplified variant of the SYK model obtained by deleting all but a Θ ( n ) fraction of the terms in a uniformly i.i.d. manner. We find the average (over the random term removal) gate complexity for simulating this model using higher-order formulas scales with O ~ ( n 1 + 1 2 t ) for even k and O ~ ( n 2 t ) for odd k . Similar to the full SYK model, we obtain a O ( n ) -reduction simulating the time evolution of an arbitrary fixed input state | ψ ⟩ .Our results highlight the potential of Lie–Trotter–Suzuki formulas for efficiently simulating the SYK and sparse SYK models, and our analytical methods can be naturally extended to other Gaussian random Hamiltonians.

This dissertation presents a comprehensive investigation of fermionic quantum matter into the intricate relationship between quantum chaos, thermalization, and their connection to gravitational duals, with particular emphasis on the Sachdev-Ye-Kitaev (SYK) model and its extensions. Through an integrated approach combining analytical solutions in the large-N limit with advanced numerical techniques, we explore fundamental questions in quantum many-body physics, including the interplay of chaos and ergodicity, universal transport in non-Fermi liquids, and the role of gravitational analogs in strongly correlated systems. We first demonstrate that coupled SYK lattices evade instantaneous thermalization despite weak charge fluctuations, governed by a discrete wave equation dependent solely on hopping strength, independent of on-site interactions. Thermodynamically, coupled SYK models exhibit critical exponents universal to van der Waals fluids and AdS black holes, yet diverge from the Hawking-Page transition at ultralow temperatures, revealing a continuum of first-order phase transitions. They simultaneously retain the quantum chaotic properties of both a single-dot SYK model and black holes. In extended SYK chains, we identify a universal DC conductivity bound across interaction regimes, identifying crossovers between insulating, Fermi liquid, strange metal, and bad metal phases with signatures of holographic insulator behavior. Quench dynamics in mixed SYK systems reveal rapid thermalization without prethermal plateaus, enabling direct observation of thermalization in closed quantum systems in the thermodynamic limit. To diagnose ergodicity-breaking, we introduce the Krylov variance, a probe of operator localization in Krylov space, validated through mass-deformed SYK and quantum East models. Finally, we uncover Page curve entanglement dynamics and temporal quantum phase transitions in entanglement Hamiltonians of interacting fermionic chains, where symmetry-sector crossings highlight limitations of perturbative calculations performed around initial times, such as in Hawking’s analysis of black hole evaporation. We provide a condensed matter analog of the holographic “island formula” — a mechanism invoking quantum extremal surfaces to resolve the black hole information paradox. Methodologically, this work integrates large-N-type expansions, Keldysh contour techniques, and numerical methods to bridge condensed matter physics, quantum information, and gravitational duality. The results establish SYK-type systems as a universal framework for studying non-Fermi liquids, quantum criticality, and black hole analogs, while paving the way for exploring bosonic extensions and generalizations to open quantum systems.

The Sachdev-Ye-Kitaev (SYK) model describes Majorana fermions with random interaction, which displays many interesting properties such as non-Fermi liquid behavior, quantum chaos, emergent conformal symmetry and holographic duality. Here we consider a SYK model or a chain of SYK models with N Majorana fermion modes coupled to another SYK model with N2 Majorana fermion modes, in which the latter has many more degrees of freedom and plays the role as a thermal bath. For a single SYK model coupled to the thermal bath, we show that although the Lyapunov exponent is still proportional to temperature, it monotonically decreases from 2π/β (β = 1/(kBT), T is temperature) to zero as the coupling strength to the thermal bath increases. For a chain of SYK models, when they are uniformly coupled to the thermal bath, we show that the butterfly velocity displays a crossover from a sqrt{T} -dependence at relatively high temperature to a linear T-dependence at low temperature, with the crossover temperature also controlled by the coupling strength to the thermal bath. If only the end of the SYK chain is coupled to the thermal bath, the model can introduce a spatial dependence of both the Lyapunov exponent and the butterfly velocity. Our models provide canonical examples for the study of thermalization within chaotic models.

The Sachdev–Ye–Kitaev (SYK) model is a strongly coupled, quantum many-body system that is chaotic, nearly conformally invariant, and exactly solvable. This remarkable and, to date, unique combination of properties have driven the intense activity surrounding the SYK model and its applications within both high energy and condensed matter physics. In this review we give an introduction to the SYK model and recent developments connected to it. We discuss: SYK and tensor models as a new class of large N quantum field theories, the near-conformal invariance in the infrared, the computation of correlation functions, generalizations of the SYK model, and applications to AdS/CFT and strange metals.

We compute the thermodynamic properties of the Sachdev-Ye-Kitaev (SYK) models of fermions with a conserved fermion number, $\mathcal{Q}$. We extend a previously proposed Schwarzian effective action to include a phase field, and this describes the low temperature energy and $\mathcal{Q}$ fluctuations. We obtain higher-dimensional generalizations of the SYK models which display disordered metallic states without quasiparticle excitations, and we deduce their thermoelectric transport coefficients. We also examine the corresponding properties of Einstein-Maxwell-scalar theories on black brane geometries which interpolate from either AdS$_4$ or AdS$_5$ to an AdS$_2\times \mathbb{R}^2$ or AdS$_2\times \mathbb{R}^3$ near-horizon geometry. These provide holographic descriptions of non-quasiparticle metallic states without momentum conservation. We find a precise match between low temperature transport and thermodynamics of the SYK and holographic models. In both models the Seebeck transport coefficient is exactly equal to the $\mathcal{Q}$-derivative of the entropy. For the SYK models, quantum chaos, as characterized by the butterfly velocity and the Lyapunov rate, universally determines the thermal diffusivity, but not the charge diffusivity.

Nonperturbative formulation of the Sachdev-Ye-Kitaev (SYK) model is discussed. The partition function of the model can be represented as a functional integral over the Grassmann variables in Euclidean time which is well defined but it diverges after the transformation to the fermion bilocal fields. We point out that the generating functional of the SYK model in real time is well defined even after the transformation to the bilocal fields and it can be used for nonperturbative investigations of its properties. The SYK model in zero dimensions is studied, its large N expansion is evaluated and phase transitions are investigated.

Finding the ground state of strongly-interacting fermionic systems is often the prerequisite for fully understanding both quantum chemistry and condensed matter systems. The Sachdev–Ye–Kitaev (SYK) model is a representative example of such a system; it is particularly interesting not only due to the existence of efficient quantum algorithms preparing approximations to the ground state such as Hastings–O'Donnell (STOC 2022), but also known no-go results for many classical ansatzes in preparing low-energy states. However, this quantum-classical separation is known to not persist when the SYK model is sufficiently sparsified, i.e., when terms in the model are discarded with probability 1 − p , where p = Θ ( 1 / n 3 ) and n is the system size. This raises the question of how robust the quantum and classical complexities of the SYK model are to sparsification.In this work we initiate the study of the sparse SYK model where p ∈ [ Θ ( 1 / n 3 ) , 1 ] and show there indeed exists a certain robustness of sparsification. We prove that with high probability, Gaussian states achieve only a Θ ( 1 / n ) -factor approximation to the true ground state energy of sparse SYK for all p ≥ Ω ( log ⁡ n / n 2 ) , and that Gaussian states cannot achieve constant-factor approximations unless p ≤ O ( log 2 ⁡ n / n 3 ) . Additionally, we prove that the quantum algorithm of Hastings–O'Donnell still achieves a constant-factor approximation to the ground state energy when p ≥ Ω ( log ⁡ n / n ) . Combined, these show a provable separation between classical algorithms outputting Gaussian states and efficient quantum algorithms for the goal of finding approximate sparse SYK ground states whenever p ≥ Ω ( log ⁡ n / n ) , extending the analogous p = 1 result of Hastings–O'Donnell.

We study the chaos exponent of some variants of the Sachdev-Ye-Kitaev (SYK) model, namely, the mathcal{N} = 1 supersymmetry (SUSY)-SYK model and its sibling, the (N|M)-SYK model which is not supersymmetric, for arbitrary interaction strength. We find that for large q the chaos exponent of these variants, as well as the SYK and the mathcal{N} = 2 SUSY-SYK model, all follow a single-parameter scaling law. By quantitative arguments we further make a conjecture, i.e. that the found scaling law might hold for general one-dimensional (1D) SYK-like models with large q. This points out a universal route from maximal chaos towards completely regular or integrable motion in the SYK model and its 1D variants.

The Sachdev-Ye-Kitaev (SYK) model is a system of N Majorana fermions with random interactions and strongly chaotic dynamics, which at low energy admits a holographically dual description as two-dimensional Jackiw-Teitelboim gravity. Hence the SYK model provides a toy model of quantum gravity that might be feasible to simulate with near-term quantum hardware. Motivated by the goal of reducing the resources needed for such a simulation, we study a sparsified version of the SYK model, in which interaction terms are deleted with probability 1−p. Specifically, we compute numerically the spectral form factor (SFF, the Fourier transform of the Hamiltonian’s eigenvalue pair correlation function) and the nearest-neighbor eigenvalue gap ratio r (characterizing the distribution of gaps between consecutive eigenvalues). We find that when p is greater than a transition value p1, which scales as 1/N3, both the SFF and r match the values attained by the full unsparsified model and with expectations from random matrix theory (RMT). But for p < p1, deviations from unsparsified SYK and RMT occur, indicating a breakdown of holography in the highly sparsified regime. Below an even smaller value p2, which also scales as 1/N3, even the spacing of consecutive eigenvalues differs from RMT values, signaling a complete breakdown of spectral rigidity. Our results cast doubt on the holographic interpretation of very highly sparsified SYK models obtained via machine learning using teleportation infidelity as a loss function.

A study of Rényi entanglement entropy in the Sachdev-Ye-Kitaev (SYK) chain of Majorana fermions suggested that the model does not rapidly thermalize, despite being maximally chaotic. In this work, I examine the eigenstate thermalization hypothesis (ETH) for both the SYK chain and the two-site SYK models using exact diagonalization. I focus on two physically motivated few-body observables: the one-body particle number operator and the two-body hopping operator. These operators provide insight into local occupation and, in certain contexts, may reflect aspects of transport. Both operators have been employed in previous studies of the ETH in related SYK models. I show that single realizations of both models approximately satisfy ETH conditions for these operators, while ensemble averages strictly satisfy the ETH. Therefore, I conclude that the finite-size SYK chain and two-site SYK models can rapidly thermalize with respect to these observables through the ETH mechanism. This suggests that the subthermal behavior observed in previous studies of Rényi entanglement entropy does not manifest in finite-size systems and that these systems can thermalize rapidly via their light modes. It also indicates that the proposed gravitational dual may undergo rapid thermalization.

We study the Lindbladian dynamics of the Sachdev-Ye-Kitaev (SYK) model, where the SYK model is coupled to Markovian reservoirs with jump operators that are either linear or quadratic in the Majorana fermion operators. Here, the linear jump operators are non-random while the quadratic jump operators are sampled from a Gaussian distribution. In the limit of large $N$, where $N$ is the number of Majorana fermion operators, and also in the limit of large $N$ and $M$, where $M$ is the number of jump operators, the SYK Lindbladians are analytically tractable, and we obtain their stationary Green's functions, from which we can read off the decay rate. For finite $N$, we also study the distribution of the eigenvalues of the SYK Lindbladians.

We develop a unified minimal scheme to classify quantum chaos in the Sachdev-Ye-Kitaev (SYK) and supersymmetric (SUSY) SYK models and also work out the structure of the energy levels in one periodic table. The SYK with even q-body or SUSY SYK with odd q-body interaction, with N even or odd number of sites, are put on an equal footing in the minimal Hilbert space; N (mod 8), q (mod 4) double Bott periodicity, and a reflection condition are identified. Exact diagonalizations (EDs) are performed to study both the bulk energy level statistics and hard-edge behaviors. Excellent agreements between the ED results and the symmetry classifications are demonstrated. Our compact and systematic methods can be transformed to map out more complicated periodic tables of SYK models with more degrees of freedom, tensor models, or symmetry protected topological phases.

The random matrix theory (RMT) can be used to classify both topological phases of matter and quantum chaos. We develop a systematic and transformative RMT to classify the quantum chaos in the colored Sachdev-Ye-Kitaev (SYK) model first introduced by Gross and Rosenhaus. Here we focus on the 2-colored case and 4-colored case with balanced number of Majorana fermion $N$. By identifying the maximal symmetries, the independent parity conservation sectors, the minimum (irreducible) Hilbert space, and especially the relevant anti-unitary and unitary operators, we show that the color degree of freedoms lead to novel quantum chaotic behaviours. When $N$ is odd, different symmetry operators need to be constructed to make the classifications complete. The 2-colored case only show 3-fold Wigner-Dyson way, and the 4-colored case show 10-fold generalized Wigner-Dyson way which may also have non-trivial edge exponents. We also study 2- and 4-colored hybrid SYK models which display many salient quantum chaotic features hidden in the corresponding pure SYK models. These features motivate us to develop a systematic RMT to study the energy level statistics of 2 or 4 un-correlated random matrix ensembles. The exact diagonalizations are performed to study both the bulk energy level statistics and the edge exponents and find excellent agreements with our exact maximal symmetry classifications. Our complete and systematic methods can be easily extended to study the generic imbalanced cases. They may be transferred to the classifications of colored tensor models, quantum chromodynamics with pairings across different colors, quantum black holes and interacting symmetry protected (or enriched) topological phases.

We present a complete symmetry classification of the Sachdev-Ye-Kitaev (SYK) model with mathcal{N} = 0, 1 and 2 supersymmetry (SUSY) on the basis of the Altland-Zirnbauer scheme in random matrix theory (RMT). For mathcal{N} = 0 and 1 we consider generic q-body interactions in the Hamiltonian and find RMT classes that were not present in earlier classifications of the same model with q = 4. We numerically establish quantitative agreement between the distributions of the smallest energy levels in the mathcal{N} = 1 SYK model and RMT. Furthermore, we delineate the distinctive structure of the mathcal{N} = 2 SYK model and provide its complete symmetry classification based on RMT for all eigenspaces of the fermion number operator. We corroborate our classification by detailed numerical comparisons with RMT and thus establish the presence of quantum chaotic dynamics in the mathcal{N} =2 SYK model. We also introduce a new SYK-like model without SUSY that exhibits hybrid properties of the mathcal{N} = 1 and mathcal{N} = 2 SYK models and uncover its rich structure both analytically and numerically.

The Sachdev-Ye-Kitaev (SYK) model attracts attention in the context of information scrambling, which represents delocalization of quantum information and is quantified by the out-of-time-ordered correlators (OTOC). The SYK model contains $N$ fermions with disordered and four-body interactions. Here, we introduce a variant of the SYK model, which we refer to as the Wishart SYK model. We investigate the Wishart SYK model for complex fermions and that for hard-core bosons. We show that the ground state of the Wishart SYK model is massively degenerate and the residual entropy is extensive, and that the Wishart SYK model for complex fermions is integrable. In addition, we numerically investigate the OTOC and level statistics of the SYK models. At late times, the OTOC of the fermionic Wishart SYK model exhibits large temporal fluctuations, in contrast with smooth scrambling in the original SYK model. We argue that the large temporal fluctuations of the OTOC are a consequence of a small effective dimension of the initial state. We also show that the level statistics of the fermionic Wishart SYK model is in agreement with the Poisson distribution, while the bosonic Wishart SYK model obeys the GUE or the GOE distribution.