String-breaking statics and dynamics in a (1+1)D SU(2) lattice gauge theory

This paper develops and applies a gauge-invariant tensor-network toolkit based on the loop-string-hadron formulation to systematically study both static string tension and the complex dynamical processes of string breaking, including particle production and entanglement spreading, in a (1+1)D SU(2) lattice gauge theory with dynamical fermions.

The Gist

Imagine you are watching a game of tug-of-war, but instead of people pulling on a rope, you have invisible particles called quarks pulling on a strand of energy called a string.

In the world of particle physics (specifically inside the protons and neutrons that make up our universe), these strings are the glue that holds matter together. But here's the weird part: if you try to pull the two quarks apart, the string doesn't just snap like a rubber band. Instead, the energy in the string gets so intense that it spontaneously creates new particles out of thin air, effectively "breaking" the string into two shorter strings. This process is called string breaking, and it's the reason why we never see a single quark floating alone in nature.

For decades, scientists have struggled to simulate this process on computers because the math is incredibly messy and involves "sign problems" that make standard calculations crash.

This paper is like a team of engineers building a brand new, specialized toolkit to finally watch this tug-of-war in slow motion. Here is how they did it, explained simply:

1. The New Toolkit: "Loop-String-Hadron"

Think of the old way of doing this as trying to describe a complex dance by listing every single muscle movement of every dancer. It's overwhelming and prone to errors.

The authors used a new method called the Loop-String-Hadron (LSH) formulation. Imagine instead that you describe the dance by looking at the patterns the dancers make:

• Loops: Dancers holding hands in a circle.
• Strings: Dancers holding hands in a line.
• Hadrons: Dancers hugging in a tight group.

By focusing on these simple, local patterns, they created a "dictionary" that translates the complex, confusing rules of particle physics into a language that computers (specifically Tensor Networks) can actually understand without getting a headache. It's like translating a foreign language into English so you can read the story.

2. The Static Experiment: Measuring the Rope's Tension

First, the team wanted to know: "How strong is this rope?"
They set up a simulation where they held two quarks still and pulled them apart. They measured the energy required to stretch the string.

• The Result: They calculated the exact "tension" of the string. This is like knowing exactly how many pounds of force it takes to stretch a rubber band before it starts to snap. They found that as the string gets longer, the energy rises linearly until it hits a breaking point.

3. The Dynamic Experiment: The "Quench"

Next, they wanted to see the string break in real-time. They didn't just pull the quarks apart slowly; they gave the system a sudden "kick" (called a quench). Imagine snapping a guitar string.

They watched what happened in two different scenarios, using two different "weights" for the particles at the ends of the string:

• Scenario A: The Light Quarks (The Energetic Kids)
  When the particles at the ends were light, they were very active. As the string stretched, it didn't just break cleanly. It got chaotic!

  • The Analogy: Imagine a stretched rubber band snapping. The energy flies everywhere, creating a shower of new, smaller rubber bands and little knots.
  • What happened: The string stretched, then suddenly "dissociated." The energy rushed outward, creating a cascade of new particles (a "particle shower"). The original string dissolved completely, and the energy spread out like a ripple in a pond. The system became very "entangled" (the particles became deeply connected in a quantum way).

• Scenario B: The Heavy Quarks (The Heavyweights)
  When the particles at the ends were heavy, they moved sluggishly.

  • The Analogy: Imagine trying to snap a thick, heavy steel cable. It's hard to get it moving.
  • What happened: The string stretched, but it didn't break as violently. The energy stayed more localized. The "shower" of new particles was much smaller. The string remained somewhat intact for longer, and the energy didn't spread out as fast. It was a much quieter, more controlled event.

4. Why This Matters

This paper is a breakthrough for a few reasons:

• No More Guessing: Previous models of how particles break apart (hadronization) were based on guesswork and fitting parameters to experiments. This study calculates it from first principles—pure math and physics.
• A New Lens: The "Loop-String-Hadron" view gives scientists a clear, intuitive way to see exactly what is happening inside the quantum soup. They can now see the "microscopic" steps: the string stretching, the particles splitting, and the energy flowing.
• Future Tech: This toolkit is designed to be scalable. Just as they used it for a simple 1D line of particles, they can now use the same logic to study more complex, 3D versions of these theories. This paves the way for understanding how the universe formed right after the Big Bang and how particle colliders (like the Large Hadron Collider) create new matter.

In summary: The authors built a new, clearer pair of glasses (the LSH formulation) to watch a chaotic quantum dance (string breaking). They found that light particles make a messy, energetic explosion when the string breaks, while heavy particles make a quieter, slower event. This helps us understand the fundamental rules of how matter is created and destroyed in our universe.

Technical Summary

1. Problem Statement

The hadronization process in high-energy particle collisions, where color-charged quarks and gluons transform into color-neutral hadrons, is fundamentally driven by string breaking. In Quantum Chromodynamics (QCD), this involves the formation of a flux tube (string) between color charges, which eventually breaks via the Schwinger mechanism (pair production) when the energy density is sufficient.

However, simulating these real-time (Minkowski-time) dynamics from first principles is notoriously difficult:

• Sign Problem: Conventional Monte Carlo methods for Lattice Gauge Theories (LGTs) fail in real-time evolution due to the sign problem.
• Non-Abelian Complexity: While Tensor Networks (TNs) have successfully simulated Abelian theories (like the Schwinger model), extending them to non-Abelian theories (like SU(2) or SU(3)) is challenging. The primary obstacle is the non-commutativity of Gauss's laws, which makes it difficult to construct basis states that are simultaneously local and gauge-invariant.
• Existing Limitations: Previous studies often relied on truncated Hilbert spaces, Gaussian variational ansatzes, or lacked fully dynamical fermions, limiting their ability to capture the full microscopic complexity of string breaking and particle production.

2. Methodology

The authors employ a Tensor Network (TN) approach based on the Loop-String-Hadron (LSH) formulation of the (1+1)D SU(2) Lattice Gauge Theory coupled to dynamical fermions.

A. The Loop-String-Hadron (LSH) Formalism

Instead of using the standard Kogut-Susskind formulation, the authors utilize the LSH reformulation, which:

• Localizes Degrees of Freedom: The Hilbert space at each lattice site is spanned by local, gauge-invariant "snapshots" of non-local observables:
  • Loops ($n_l$): Local snapshots of Wilson loops or background flux.
  • Strings ($n_i, n_o$): Incoming and outgoing string excitations representing fermions sourcing flux.
  • Hadrons: Local meson/baryon-like excitations.
• Simplifies Constraints: The non-Abelian Gauss's law constraints are intrinsically built into the local basis states. The only constraints remaining are Abelian Gauss's Laws (AGLs) ensuring electric flux continuity between sites. This allows the use of standard TN packages that handle Abelian symmetries efficiently.
• Hamiltonian: The Hamiltonian is expressed entirely in terms of LSH operators (number operators and ladder/string operators), eliminating the need for non-local gauge transformations during time evolution.

B. Tensor Network Ansatzes

• Static Studies: The authors use Matrix Product States (MPS) and Matrix Product Operators (MPO) with the Density Matrix Renormalization Group (DMRG) algorithm to compute ground states and static potentials.
• Dynamical Studies: Real-time evolution is performed using the Time-Dependent Variational Principle (TDVP) algorithm (specifically 1-site and 2-site variants) coupled with a Krylov-basis expansion to dynamically adjust the bond dimension ($D_{max}$).
• Symmetry Handling: Global symmetries (total fermion number and flux difference) are explicitly incorporated into the MPS/MPO tensors via block-sparse structures, ensuring gauge invariance without penalty terms.

3. Key Contributions

1. First Comprehensive LSH-TN Study of SU(2) Dynamics: This work extends previous LSH applications to a full dynamical study of string breaking with fully dynamical fermions, moving beyond static charge approximations or Gaussian ansatzes.
2. Microscopic Diagnosis of Processes: The LSH framework provides a unique, intuitive language to diagnose specific physical processes during string breaking, such as:
   • String expansion/contraction.
   • Endpoint splitting and particle showers.
   • Inelastic scattering (string dissociation and recombination).
   • Particle production (mesons and baryons).
3. Systematic Control of Truncation: The authors demonstrate the critical importance of the electric flux cutoff ($j_{max}$). They show that low cutoffs yield qualitatively different (and incorrect) dynamics compared to higher cutoffs, establishing a benchmark for convergence in non-Abelian TN simulations.
4. Continuum and Thermodynamic Limits: The study successfully extrapolates static string tension to the continuum and thermodynamic limits, providing a precise benchmark value.

4. Key Results

A. Static String Breaking

• Static Potential: The authors computed the static potential between two charges. It exhibits a linear rise (confinement) followed by a plateau (string breaking) at a critical separation.
• String Tension: In the continuum limit for a fermion mass $m/g = 0.5$, the string tension was determined to be:
  $$\frac{\kappa}{g^2} = 0.330(3)$$
• Mass Dependence: The string tension increases with fermion mass, approaching the pure Yang-Mills value ($\approx 0.375$) for heavy masses, as fermions become effectively static.

B. Dynamical String Breaking (Quench Dynamics)

The authors initiated a "meson-string" state (two dynamical fermions connected by a flux tube) from the interacting vacuum and evolved it in time. They compared two regimes: light fermion mass ($m/g = 0.2$) and heavy fermion mass ($m/g = 2.0$).

• Light Mass Regime (Inelastic/Dissociative):

  • Mechanism: The string endpoints move rapidly, leading to string dissociation (breaking) and recombination.
  • Particle Production: There is abundant production of bare mesons and baryons. The flux in the center of the string is rapidly depleted (screened).
  • Transport: Energy transport is ballistic.
  • Entanglement: Entanglement entropy grows rapidly, and the system exhibits complex scattering events where wavefronts collide, creating "showers" of excitations.
  • Correlations: Long-range correlations are suppressed as the flux tube is screened; the system transitions to a state of separated mesonic pairs.

• Heavy Mass Regime (Elastic/Suppressed):

  • Mechanism: The heavy mass penalizes particle creation. The string endpoints move slower, and string dissociation is suppressed.
  • Particle Production: Very little particle production occurs. The string remains largely intact, stretching and contracting but not fully breaking into new hadrons.
  • Transport: Energy transport is non-ballistic and non-diffusive.
  • Flux: The electric flux in the interior is only partially screened; the flux tube remains connected over long distances.
  • Correlations: Strong long-range correlations persist between the endpoints, indicating the survival of the stretched flux tube.

C. Convergence and Truncation Effects

• Flux Cutoff ($j_{max}$): Simulations with minimal cutoffs ($j_{max}=1/2$) showed oscillatory, confined dynamics that differed drastically from the converged results ($j_{max}=5/2$). This highlights that capturing the correct real-time physics in non-Abelian theories requires a sufficiently large Hilbert space.
• Bond Dimension: The results were shown to be converged for $D_{max} = 600$ for the heavy mass case, while the light mass case required careful monitoring of relative errors (limited to $gt=8$ due to rapid entanglement growth).

5. Significance and Outlook

• Benchmarking: This work provides a rigorous first-principles benchmark for string-breaking dynamics in non-Abelian gauge theories, essential for validating phenomenological models (like the Lund model) used in collider physics.
• Methodological Advancement: It demonstrates that the LSH formulation combined with TNs is a powerful tool for overcoming the sign problem and handling non-Abelian constraints without sacrificing locality.
• Future Directions:
  • The framework is directly applicable to (1+1)D SU(3) and can be extended to (2+1)D using PEPS or Tree TNs.
  • It paves the way for studying scattering processes and hadronization rates in more complex theories.
  • The results serve as a benchmark for quantum simulation experiments on quantum computers and simulators, which are beginning to tackle similar lattice gauge theory problems.

In summary, this paper successfully bridges the gap between static confinement studies and real-time hadronization dynamics in non-Abelian gauge theories, offering a microscopic, gauge-invariant, and computationally robust description of how strings break and particles are produced.