How coherent confining strings reorganize into hadronic matter is a central open problem in nonequilibrium gauge dynamics. Using real-time tensor-network simulations of the massive Schwinger model, we identify a collective organization underlying string fragmentation. For intermediate masses, the dynamics simplifies and more than $99\%$ of the temporal evolution of the electric-field probability distribution is governed by a single collective mode whose decay tracks the survival probability of the parent string. The extracted survival law exhibits the same Gaussian dependence implied by the Lund world-sheet area law, providing microscopic support for its characteristic area law dependence without assuming the Lund fragmentation mechanism. The same mode controls the transition from the coherent string to the produced mesons, explains the universal non-thermal effective fluid equation of state, and essentially controls the growth of bipartite entanglement through the causal propagation of massive pairs. These results identify collective fragmentation as an organizing principle of coherent string breaking and provide a microscopic perspective on the dynamical origin of a central ingredient of Lund hadronization.
How coherent confining strings reorganize into hadronic matter is a central open problem in nonequilibrium gauge dynamics. Using real-time tensor-network simulations of the massive Schwinger model, we identify a collective organization underlying string fragmentation. For intermediate masses, the dynamics simplifies and more than 99% of the temporal evolution of the electric-field probability distribution is governed by a single collective mode whose decay tracks the survival probability of the parent string. The extracted survival law exhibits the same Gaussian dependence implied by the Lund world-sheet area law, providing microscopic support for its characteristic area law dependence without assuming the Lund fragmentation mechanism. The same mode controls the transition from the coherent string to the produced mesons, explains the universal non-thermal effective fluid equation of state, and essentially controls the growth of bipartite entanglement through the causal propagation of massive pairs. These results identify collective fragmentation as an organizing principle of coherent string breaking and provide a microscopic perspective on the dynamical origin of a central ingredient of Lund hadronization.
We simulate the real-time dynamics of a massive Schwinger model using the Time-Evolving Block Decimation tensor network algorithm. Starting from a non-equilibrium initial state with localized energy excitation on top of vacuum, we track the subsequent evolution to investigate two distinct physical phenomena. First, by analyzing the system's energy-momentum tensor, we show that this system exhibits hydrodynamic behavior analogous to Bjorken flow at large coupling-to-mass ratio, a signature that diminishes as the coupling weakens, or mass increases. Second, by examining the evolution of the electric field and charge density, we observe the signal of spontaneous parity symmetry breaking phase transition in a dynamical system. The parity-restored regime is marked by ''string breaking'' and efficient charge screening, while the parity-broken regime displays stable propagation of nearly free charges and persistent electric fields connecting them.
String breaking, the process by which flux tubes fragment into hadronic states, is a hallmark of confinement in strongly-interacting quantum field theories. A suite of quantum complexity measures is examined using Matrix Product States to characterize the string breaking process in the 1+1D Schwinger model. We demonstrate the presence of nonlocal quantum correlations along the string that may affect fragmentation dynamics, and show that entanglement and magic offer complementary perspectives on string formation and breaking beyond conventional observables.
Understanding the dynamics of bound state formation is one of the fundamental questions in confining quantum field theories such as Quantum Chromodynamics (QCD). One hadronization mechanism that has garnered significant attention is the breaking of a string initially connecting a fermion and an antifermion. Deepening our understanding of real-time string-breaking dynamics with simpler, lower dimensional models like the Schwinger model can improve our understanding of the hadronization process in QCD and other confining systems found in condensed matter and statistical systems. In this paper, we consider the string-breaking dynamics within the Schwinger model and investigate its modification inside a thermal medium, treating the Schwinger model as an open quantum system coupled to a thermal environment. Within the regime of weak coupling between the system and environment, the real-time evolution of the system can be described by a Lindblad evolution equation. We analyze the Liouvillian gaps of this Lindblad equation and the time dependence of the system's von Neumann entropy. We observe that the late-time relaxation rate decreases as the environment correlation length increases. Moreover, when the environment correlation length is infinite, the system exhibits two steady states, one in each of the sectors with definite charge-conjugation-parity (CP) quantum numbers. For parameter regimes where an initial string breaks in vacuum, we observe a delay of the string breaking in the medium, due to kinetic dissipation effects. Conversely, in regimes where an initial string remains intact in vacuum time evolution, we observe string breaking (melting) in the thermal medium. We further discuss how the Liouvillian dynamics of the open Schwinger model can be simulated on quantum computers and provide an estimate of the associated Trotter errors.
Abstract The spontaneous formation of particle pairs is one of the most intriguing phenomena in nature. A canonical example is when the potential energy between two elementary particles increases with their separation, as if the particles were confined to one another by a tense string. When the separation exceeds a critical value, new particle pairs can form, which causes the string to break. String-breaking dynamics in quantum chromodynamics, the theory of the strong force, play a fundamental role in high-energy particle collisions and the evolution of matter in the early universe. Simulating the evolution of strings and the formation of composite particles, or hadrons, is a grand challenge of modern physics. Quantum simulators are well suited to study dynamical processes. They are, therefore, expected to provide an edge over conventional methods based on classical computing. However, the experimental capabilities required to simulate the string-breaking phenomenon have not yet been demonstrated, even for simpler prototypical models of the strong force. In this work, we probe, for the first time on a quantum simulator, the spatiotemporally-resolved dynamics of string breaking in a (1+1)-dimensional Z2 lattice gauge theory, a theory that also exhibits confinement of charges. We employ a fully programmable trapped-ion quantum simulator, and emulate the effects of external static charges and strings via site-dependent control of magnetic fields, using a dual array of tightly focused laser beams targeting individual ions. First, we study the effect of confinement on the evolution of isolated charges. We find that these charges freely spread in the absence of string tension, but exhibit localized coherent oscillations as the string tension is increased. Then, we observe and characterize the breaking dynamics of a string initially stretched between two static charges, following an abrupt increase of the string tension. We find that charge pairs appear near the string edges and then spread out into the bulk, thereby identifying a route to dynamical string breaking that is distinct from the conventional Schwinger mechanism. This work, therefore, demonstrates that analog quantum simulators have achieved the control needed to uncover features of string-breaking dynamics, which may ultimately be relevant to nuclear and high-energy physics.
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.
Open-system lattice gauge theory (LGT) has so far been developed predominantly in Abelian settings, leaving open how genuinely non-Abelian gauge structure reshapes dissipative real-time dynamics. Here, we study a 1+1D SU(2) Yang--Mills LGT with dynamical matter coupled to a thermal scalar environment through a gauge-preserving Lindblad evolution, which we solve using tensor networks. Starting from a quark--antiquark pair connected by a chromoelectric flux string, we find that the thermal medium melts the string by delocalizing the color charges and screening the flux; on resonance, this dissipative melting competes with and delays coherent string breaking. The thermalization time is non-monotonic in the environment coupling, decreasing through environment-assisted transport at weak dissipation before increasing in a quantum-Zeno regime. In the strong-dephasing limit, a Schrieffer--Wolff expansion maps the dynamics to a classical exclusion process and yields the Liouvillian thermalization time analytically. Beyond these generic open-system effects, the non-Abelian matter structure produces a systematic mesonic bias in the steady state, while the thermalization time decreases with temperature, in contrast to the Abelian Schwinger model trend and in qualitative agreement with pNRQCD studies of the quark--gluon plasma. These results establish a gauge-preserving framework for thermalization and string dynamics in open non-Abelian lattice gauge theories.