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

Thermalization from quantum entanglement: Jet simulations in the massive Schwinger model

Adrien Florio, David Frenklakh, Sebastian Grieninger, Dmitri E. Kharzeev, Andrea Palermo, Shuzhe Shi

Published Nov 7, 2025• 10 citations
arXivHTMLPDF

Abstract

We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.

DOI
10.1103/sgrx-jpp9

Related papers

Quantum real-time evolution of entanglement and hadronization in jet production: Lessons from the massive Schwinger model

Adrien Florio, David Frenklakh, Kazuki Ikeda +4 more

The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024

arXiv:2404.00087•Physical review. D/Physical review. D.•Nov 15, 2024•23 citations

Quantum simulation of entanglement and hadronization in jet production: lessons from the massive Schwinger model

Shuzhe Shi, Adrien Florio, David Frenklakh +4 more

In this work, we investigate the potential connection between entanglement and thermalization, as well as the dynamics of hadronization in QCD physics. Through simulations of quantum dynamics on classical hardware, we analyze the real-time response of the massive Schwinger model coupled to external sources, which serves as an analog for the production and fragmentation of quark jets. This analogy is justified by the shared properties of confinement and chiral symmetry breaking between the Schwinger model and QCD. We observe the growth of time-dependent entanglement between the produced jets, which is driven by the increasing number of contributing eigenstates of the reduced density matrix with sufficiently large and closely spaced eigenvalues. This growth is intrinsically linked to the process of thermalization. The emergent eigenstates are identified as meson-like bound states, and by tracking their temporal evolution, we are able to observe real-time hadronization. At mid-rapidity, the long-time values of local observables approach approximately constant values, indicating the onset of equilibrium and the approach to thermalization.

Conference on Quark Confinement and the Hadron Spectrum (Confinement)•Oct 24, 2025•0 citations

Onset of Bjorken Flow in Quantum Evolution of the Massive Schwinger Model

Shao, Haiyang, Shile Chen, Shuzhe Shi

The onset of hydrodynamics in the hot medium created in relativistic heavy-ion collisions is a crucial theoretical question. A first-principle simulation requires a real-time, non-perturbative calculation of the quantum system. In this Letter, we perform such simulations using the tensor network method, which enables large-scale quantum many-body simulations by retaining only the most essential quantum states for collective behaviors. We focus on the massive Schwinger model, a low-dimensional analog of quantum chromodynamics (QCD), as they share important properties such as confinement and chiral symmetry breaking. Starting from an initial state that puts a localized excitation atop the vacuum and mimics the energy deposition from colliding nuclei, we observe hydrodynamic behavior consistent with Bjorken flow in all relevant degrees of freedom: energy density, fluid velocity, and bulk pressure. The time scale for hydrodynamic onset aligns with the thermalization time of the quantum distribution function.

arXiv:2509.10855•Sep 13, 2025•0 citations

Quantum thermalization of Quark-Gluon Plasma

Shile Chen, Li Yan, Shuzhe Shi

The thermalization of quark gluon plasma created in relativistic heavy-ion collisions is a crucial theoretical question in understanding the onset of hydrodynamics, and in a broad sense, a key step to the exploration of thermalization in isolated quantum systems. Addressing this problem theoretically, in a first principle manner, requires a real-time, non-perturbative method. To this end, we carry out a fully quantum simulation on a classical hardware, of a massive Schwinger model, which well mimics QCD as it shares the important properties such as confinement and chiral symmetry breaking. We focus on the real-time evolution of the Wigner function, namely, the two-point correlation function, which approximates quark momentum distribution. In the context of the eigenstate thermalization hypothesis and the evolution of entropy, our solution reveals the emergence of quantum thermalization in quark-gluon plasma with a strong coupling constant, while thermalization fails progressively as a consequence of the gradually increased significance of quantum many-body scar states in a more weakly coupled system. More importantly, we observe the non-trivial role of the topological vacuum in thermalization, as the thermalization properties differ dramatically in the parity-even and parity-odd components of the Wigner function.

arXiv:2412.00662•Dec 1, 2024•0 citations

Dynamics of entanglement in expanding quantum fields

Jürgen Berges

Here, we develop a novel real-time approach to computing the entanglement between spatial regions for Gaussian states in quantum field theory. The entanglement entropy is characterized in terms of local correlation functions on space-like Cauchy hypersurfaces. The framework is applied to explore an expanding light cone geometry in the particular case of the Schwinger model for quantum electrodynamics in 1+1 space-time dimensions. We observe that the entanglement entropy becomes extensive in rapidity at early times and that the corresponding local reduced density matrix is a thermal density matrix for excitations around a coherent field with a time dependent temperature. Since the Schwinger model successfully describes many features of multiparticle production in e+e- collisions, our results provide an attractive explanation in this framework for the apparent thermal nature of multiparticle production even in the absence of significant final state scattering.

•May 24, 2022•65 citations

Dynamics of entanglement in expanding quantum fields

Jürgen Berges, Stefan Floerchinger, Raju Venugopalan

Abstract We develop a functional real-time approach to computing the entanglement between spatial regions for Gaussian states in quantum field theory. The entanglement entropy is characterized in terms of local correlation functions on space-like Cauchy hypersurfaces. The framework is applied to explore an expanding light cone geometry in the particular case of the Schwinger model for quantum electrodynamics in 1+1 space-time dimensions. We observe that the entanglement entropy becomes extensive in rapidity at early times and that the corresponding local reduced density matrix is a thermal density matrix for excitations around a coherent field with a time dependent temperature. Since the Schwinger model successfully describes many features of multiparticle production in e + e − collisions, our results provide an attractive explanation in this framework for the apparent thermal nature of multiparticle production even in the absence of significant final state scattering.

arXiv:1712.09362•Journal of High Energy Physics•Apr 1, 2018•79 citations
We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.
The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024
In this work, we investigate the potential connection between entanglement and thermalization, as well as the dynamics of hadronization in QCD physics. Through simulations of quantum dynamics on classical hardware, we analyze the real-time response of the massive Schwinger model coupled to external sources, which serves as an analog for the production and fragmentation of quark jets. This analogy is justified by the shared properties of confinement and chiral symmetry breaking between the Schwinger model and QCD. We observe the growth of time-dependent entanglement between the produced jets, which is driven by the increasing number of contributing eigenstates of the reduced density matrix with sufficiently large and closely spaced eigenvalues. This growth is intrinsically linked to the process of thermalization. The emergent eigenstates are identified as meson-like bound states, and by tracking their temporal evolution, we are able to observe real-time hadronization. At mid-rapidity, the long-time values of local observables approach approximately constant values, indicating the onset of equilibrium and the approach to thermalization.
The onset of hydrodynamics in the hot medium created in relativistic heavy-ion collisions is a crucial theoretical question. A first-principle simulation requires a real-time, non-perturbative calculation of the quantum system. In this Letter, we perform such simulations using the tensor network method, which enables large-scale quantum many-body simulations by retaining only the most essential quantum states for collective behaviors. We focus on the massive Schwinger model, a low-dimensional analog of quantum chromodynamics (QCD), as they share important properties such as confinement and chiral symmetry breaking. Starting from an initial state that puts a localized excitation atop the vacuum and mimics the energy deposition from colliding nuclei, we observe hydrodynamic behavior consistent with Bjorken flow in all relevant degrees of freedom: energy density, fluid velocity, and bulk pressure. The time scale for hydrodynamic onset aligns with the thermalization time of the quantum distribution function.
The thermalization of quark gluon plasma created in relativistic heavy-ion collisions is a crucial theoretical question in understanding the onset of hydrodynamics, and in a broad sense, a key step to the exploration of thermalization in isolated quantum systems. Addressing this problem theoretically, in a first principle manner, requires a real-time, non-perturbative method. To this end, we carry out a fully quantum simulation on a classical hardware, of a massive Schwinger model, which well mimics QCD as it shares the important properties such as confinement and chiral symmetry breaking. We focus on the real-time evolution of the Wigner function, namely, the two-point correlation function, which approximates quark momentum distribution. In the context of the eigenstate thermalization hypothesis and the evolution of entropy, our solution reveals the emergence of quantum thermalization in quark-gluon plasma with a strong coupling constant, while thermalization fails progressively as a consequence of the gradually increased significance of quantum many-body scar states in a more weakly coupled system. More importantly, we observe the non-trivial role of the topological vacuum in thermalization, as the thermalization properties differ dramatically in the parity-even and parity-odd components of the Wigner function.
Here, we develop a novel real-time approach to computing the entanglement between spatial regions for Gaussian states in quantum field theory. The entanglement entropy is characterized in terms of local correlation functions on space-like Cauchy hypersurfaces. The framework is applied to explore an expanding light cone geometry in the particular case of the Schwinger model for quantum electrodynamics in 1+1 space-time dimensions. We observe that the entanglement entropy becomes extensive in rapidity at early times and that the corresponding local reduced density matrix is a thermal density matrix for excitations around a coherent field with a time dependent temperature. Since the Schwinger model successfully describes many features of multiparticle production in e+e- collisions, our results provide an attractive explanation in this framework for the apparent thermal nature of multiparticle production even in the absence of significant final state scattering.
Abstract We develop a functional real-time approach to computing the entanglement between spatial regions for Gaussian states in quantum field theory. The entanglement entropy is characterized in terms of local correlation functions on space-like Cauchy hypersurfaces. The framework is applied to explore an expanding light cone geometry in the particular case of the Schwinger model for quantum electrodynamics in 1+1 space-time dimensions. We observe that the entanglement entropy becomes extensive in rapidity at early times and that the corresponding local reduced density matrix is a thermal density matrix for excitations around a coherent field with a time dependent temperature. Since the Schwinger model successfully describes many features of multiparticle production in e + e − collisions, our results provide an attractive explanation in this framework for the apparent thermal nature of multiparticle production even in the absence of significant final state scattering.