© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2110.00017 · Journal of Statistical Mechanics Theory and Experiment

Thermalization of randomly coupled SYK models

Ramanjit Sohal, Laimei Nie, Xiao-Qi Sun, Eduardo Fradkin

Published Jan 1, 2022• 8 citations

Abstract

Abstract We investigate the thermalization of Sachdev–Ye–Kitaev (SYK) models coupled via random interactions following quenches from the perspective of entanglement. Previous studies have shown that when a system of two SYK models coupled by random two-body terms is quenched from the thermofield double state with sufficiently low effective temperature, the Rényi entropies do not saturate to the expected thermal values in the large- N limit. Using numerical large- N methods, we first show that the Rényi entropies in a pair SYK models coupled by two-body terms can thermalize, if quenched from a state with sufficiently high effective temperature, and hence exhibit state-dependent thermalization. In contrast, SYK models coupled by single-body terms appear to always thermalize. We provide evidence that the subthermal behavior in the former system is likely a large- N artifact by repeating the quench for finite N and finding that the saturation value of the Rényi entropy extrapolates to the expected thermal value in the N → ∞ limit. Finally, as a finer grained measure of thermalization, we compute the late-time spectral form factor of the reduced density matrix after the quench. While a single SYK dot exhibits perfect agreement with random matrix theory, both the quadratically and quartically coupled SYK models exhibit slight deviations.

DOI: 10.1088/1742-5468/ac416b

Related papers

Spread of entanglement in a Sachdev-Ye-Kitaev chain

Yingfei Gu, Andrew Lucas, Xiao-Liang Qi

We study the spread of Rényi entropy between two halves of a Sachdev-Ye-Kitaev (SYK) chain of Majorana fermions, prepared in a thermofield double (TFD) state. The SYK chain model is a model of chaotic many-body systems, which describes a one-dimensional lattice of Majorana fermions, with spatially local random quartic interaction. We find that for integer Rényi index n > 1, the Rényi entanglement entropy saturates at a parametrically smaller value than expected. This implies that the TFD state of the SYK chain does not rapidly thermalize, despite being maximally chaotic: instead, it rapidly approaches a prethermal state. We compare our results to the signatures of thermalization observed in other quenches in the SYK model, and to intuition from nearly-AdS2 gravity.

arXiv:1708.00871•Journal of High Energy Physics•Sep 1, 2017•109 citations

Subsystem Rényi entropy of thermal ensembles for SYK-like models

Pengfei Zhang, Chunxiao Liu, Xiao Chen

The Sachdev-Ye-Kitaev model is an N N -modes fermionic model with infinite range random interactions. In this work, we study the thermal Rényi entropy for a subsystem of the SYK model using the path-integral formalism in the large- N N limit. The results are consistent with exact diagonalization and can be well approximated by thermal entropy with an effective temperature when subsystem size M\leq N/2 M ≤ N / 2 . We also consider generalizations of the SYK model with quadratic random hopping term or U(1) U ( 1 ) charge conservation.

arXiv:2003.09766•SciPost Physics•Jun 26, 2020•43 citations

Eigenstate thermalization in the two-site SYK and SYK chain models

Seyyed M. H. Halataei

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.

arXiv:2104.05291•Physical review. D/Physical review. D.•Jun 1, 2025•2 citations

State-dependent thermalization in SYK model

Tousik Samui, Nilakash Sorokhaibam

We show that the integrable state of charged SYK model does not thermalize even though the state has non-zero Lyapunov exponent. The chaotic state always thermalizes which is expected. So whether the system thermalizes or not depends on the initial state. We also show that the introduction of random mass deformation (q=2 SYK term) slows down thermalization but the system thermalizes exponentially fast. This is observed despite the fact that large (q=2) SYK interaction forces spectral statistics to obey Poisson statistics. This means that a precise measure of chaos and its prediction of thermalization is still lacking. In the chaotic state, the effective temperature is non-monotonic. It has a bump at relatively long time before settling down to the final value. With non-zero chemical potential, the effective temperature oscillates noticeably before settling down to the final value. The spectral asymmetry frequency does not change during quantum quenches.

•Apr 29, 2020•1 citation

Thermalization of a closed Sachdev-Ye-Kitaev system in the thermodynamic limit

Santiago Salazar Jaramillo, Rishabh Jha, Stefan K. Kehrein

The question of thermalization of a closed quantum system is of central interest in nonequilibrium quantum many-body physics. Here we present one such study analyzing the dynamics of a closed coupled Majorana SYK system. We have a large- q SYK model prepared initially at equilibrium quenched by introducing a random hopping term, thus leading to nonequilibrium dynamics. We find that the final stationary state reaches thermal equilibrium with respect to the Green's functions and energy. Accordingly, the final state is characterized by calculating its final temperature and the thermalization rate. We provide a detailed review of analytical methods and derive the required Kadanoff-Baym equations, which are then solved using the algorithm developed in this work. Our results display rich thermalization dynamics in a closed quantum system in the thermodynamic limit.

arXiv:2411.12421•Physical review. B./Physical review. B•May 23, 2025•5 citations

Cooling the Sachdev-Ye-Kitaev model using thermofield double states

Thomas Schuster, Bryce Kobrin, Vincent P. Su +2 more

We analyze a simple and efficient experimental protocol to cool the Sachdev-Ye-Kitaev (SYK) model to low temperatures. The protocol utilizes local couplings between two copies of an SYK model to create a gapped adiabatic path, between a high temperature product state and a low temperature thermofield double state. By smoothly varying the coupling strength between these two limits, one efficiently cools the SYK model. We support these predictions-and demonstrate fast cooling to the low-temperature gravitational regime-via exact numerical solutions to the large-N equations of motion that govern the ground state and dynamical properties of the coupled system. Finally, we present a theoretical framework based upon eigenstate thermalization that provides a microscopic explanation for the efficacy of the cooling protocol; intriguingly, this suggests that the protocol may be applicable for cooling strongly-interacting quantum Hamiltonians more broadly.

quant-pharXiv:2511.09620•Nov 14, 2025•0 citations

Abstract We investigate the thermalization of Sachdev–Ye–Kitaev (SYK) models coupled via random interactions following quenches from the perspective of entanglement. Previous studies have shown that when a system of two SYK models coupled by random two-body terms is quenched from the thermofield double state with sufficiently low effective temperature, the Rényi entropies do not saturate to the expected thermal values in the large- N limit. Using numerical large- N methods, we first show that the Rényi entropies in a pair SYK models coupled by two-body terms can thermalize, if quenched from a state with sufficiently high effective temperature, and hence exhibit state-dependent thermalization. In contrast, SYK models coupled by single-body terms appear to always thermalize. We provide evidence that the subthermal behavior in the former system is likely a large- N artifact by repeating the quench for finite N and finding that the saturation value of the Rényi entropy extrapolates to the expected thermal value in the N → ∞ limit. Finally, as a finer grained measure of thermalization, we compute the late-time spectral form factor of the reduced density matrix after the quench. While a single SYK dot exhibits perfect agreement with random matrix theory, both the quadratically and quartically coupled SYK models exhibit slight deviations.

We study the spread of Rényi entropy between two halves of a Sachdev-Ye-Kitaev (SYK) chain of Majorana fermions, prepared in a thermofield double (TFD) state. The SYK chain model is a model of chaotic many-body systems, which describes a one-dimensional lattice of Majorana fermions, with spatially local random quartic interaction. We find that for integer Rényi index n > 1, the Rényi entanglement entropy saturates at a parametrically smaller value than expected. This implies that the TFD state of the SYK chain does not rapidly thermalize, despite being maximally chaotic: instead, it rapidly approaches a prethermal state. We compare our results to the signatures of thermalization observed in other quenches in the SYK model, and to intuition from nearly-AdS2 gravity.

The Sachdev-Ye-Kitaev model is an N N -modes fermionic model with infinite range random interactions. In this work, we study the thermal Rényi entropy for a subsystem of the SYK model using the path-integral formalism in the large- N N limit. The results are consistent with exact diagonalization and can be well approximated by thermal entropy with an effective temperature when subsystem size M\leq N/2 M ≤ N / 2 . We also consider generalizations of the SYK model with quadratic random hopping term or U(1) U ( 1 ) charge conservation.

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 show that the integrable state of charged SYK model does not thermalize even though the state has non-zero Lyapunov exponent. The chaotic state always thermalizes which is expected. So whether the system thermalizes or not depends on the initial state. We also show that the introduction of random mass deformation (q=2 SYK term) slows down thermalization but the system thermalizes exponentially fast. This is observed despite the fact that large (q=2) SYK interaction forces spectral statistics to obey Poisson statistics. This means that a precise measure of chaos and its prediction of thermalization is still lacking. In the chaotic state, the effective temperature is non-monotonic. It has a bump at relatively long time before settling down to the final value. With non-zero chemical potential, the effective temperature oscillates noticeably before settling down to the final value. The spectral asymmetry frequency does not change during quantum quenches.

The question of thermalization of a closed quantum system is of central interest in nonequilibrium quantum many-body physics. Here we present one such study analyzing the dynamics of a closed coupled Majorana SYK system. We have a large- q SYK model prepared initially at equilibrium quenched by introducing a random hopping term, thus leading to nonequilibrium dynamics. We find that the final stationary state reaches thermal equilibrium with respect to the Green's functions and energy. Accordingly, the final state is characterized by calculating its final temperature and the thermalization rate. We provide a detailed review of analytical methods and derive the required Kadanoff-Baym equations, which are then solved using the algorithm developed in this work. Our results display rich thermalization dynamics in a closed quantum system in the thermodynamic limit.

We analyze a simple and efficient experimental protocol to cool the Sachdev-Ye-Kitaev (SYK) model to low temperatures. The protocol utilizes local couplings between two copies of an SYK model to create a gapped adiabatic path, between a high temperature product state and a low temperature thermofield double state. By smoothly varying the coupling strength between these two limits, one efficiently cools the SYK model. We support these predictions-and demonstrate fast cooling to the low-temperature gravitational regime-via exact numerical solutions to the large-N equations of motion that govern the ground state and dynamical properties of the coupled system. Finally, we present a theoretical framework based upon eigenstate thermalization that provides a microscopic explanation for the efficacy of the cooling protocol; intriguingly, this suggests that the protocol may be applicable for cooling strongly-interacting quantum Hamiltonians more broadly.