Quantum simulation of baryon scattering in SU(2) lattice gauge theory

This paper presents the first real-time tensor-network simulation of hadronic scattering in a (1+1)-dimensional SU(2) lattice gauge theory, revealing that while meson-meson and baryon-baryon collisions exhibit predominantly elastic dynamics similar to the Schwinger model, meson-baryon interactions display unique entanglement-driven behavior characterized by the spatial delocalization of slower wavepackets.

The Gist

Imagine the universe is built from tiny, invisible Lego bricks. Some of these bricks are "quarks," and when they snap together, they form larger structures called hadrons (like protons and neutrons, which make up your body).

Physicists have always wanted to film a "movie" of these Lego bricks crashing into each other in real-time to see how they bounce, stick, or break apart. But there's a problem: the rules of the game (Quantum Mechanics) are so complex and the forces so strong that our current supercomputers get stuck. It's like trying to solve a Rubik's cube that changes its rules every time you touch it.

This paper is about a new way to simulate these crashes using a clever mathematical trick called Tensor Networks. Think of this not as a brute-force computer, but as a highly efficient "compression algorithm" that can predict the outcome of the crash without needing to calculate every single impossible detail.

Here is the breakdown of their experiment, explained simply:

1. The Setup: A One-Dimensional Race Track

The researchers built a simplified model of the universe.

• The Track: Instead of a 3D room, they used a 1D line (a single file of people holding hands).
• The Players: They used two types of particles:
  • Mesons: Like a pair of dancers holding hands (a quark and an antiquark).
  • Baryons: Like a trio of dancers (three quarks).
• The Force: They used a specific set of rules (SU(2) gauge theory) that mimics the "strong force" holding real protons together.

2. The Experiment: Three Different Collisions

They sent these particle "wave packets" (think of them as fuzzy clouds of probability) crashing into each other from opposite ends of the line. They tested three scenarios:

Scenario A: The Polite Bounce (Meson vs. Meson)

• What happened: Two pairs of dancers ran into each other, bounced off, and kept running.
• The Result: It was a perfect, elastic collision. They didn't change shape, they didn't get tangled, and they didn't leave a mess. It was exactly like two billiard balls hitting each other.
• The Analogy: Imagine two people walking through a crowded hallway. They bump, say "excuse me," and keep walking exactly as they were before.

Scenario B: The Bouncy Castle Effect (Baryon vs. Baryon)

• What happened: Two groups of three dancers collided.
• The Result: Just like the first scenario, they bounced off each other cleanly.
• Why? The energy wasn't high enough to break them apart or create new particles. They were too heavy and moving too slowly to do anything fancy.

Scenario C: The Weird Mix (Meson vs. Baryon)

• What happened: A pair of dancers (Meson) crashed into a trio of dancers (Baryon).
• The Result: This is where it got weird.
  • The trio (Baryon) stayed mostly in one spot, like a heavy anchor.
  • The pair (Meson) didn't just bounce off. It got entangled with the trio.
  • Instead of separating cleanly, the Meson's "cloud" spread out across the entire line, while the Trio stayed put. They became a single, messy, collective state.
• The Analogy: Imagine a fast-moving skateboarder (the Meson) crashing into a heavy, stationary boulder (the Baryon). Instead of bouncing off, the skateboarder's energy spreads out like a ripple in a pond, wrapping around the boulder. They become "stuck" together in a quantum sense, even though they aren't physically glued.

3. The New Tool: The "Information Lattice"

The researchers didn't just watch the crash; they used a new tool called the Information Lattice to see how the particles were connected.

• Standard View: Usually, we measure "Entanglement Entropy," which tells us how much two things are connected (like a score of "how tangled are they?").
• The New View: The Information Lattice is like a microscope for relationships. It asks: "Is this connection just between two people, or is it a complex group hug involving three or four people?"
• The Discovery: In the weird "Meson vs. Baryon" crash, the microscope showed that the particles didn't form a neat new structure (like a new molecule). Instead, they just became a messy, spread-out cloud of shared information.

Why Does This Matter?

1. It's a Test Run: This is the first time anyone has successfully simulated this specific type of non-Abelian (complex) collision in real-time using these new tools.
2. It's a Bridge: It proves that we can use "Quantum Information Science" (the math behind quantum computers) to solve problems that traditional supercomputers can't touch.
3. The Future: If we can simulate these simple 1D crashes, we are one step closer to simulating the messy, 3D collisions that happen inside the Large Hadron Collider or inside the hearts of neutron stars.

In a nutshell: The researchers used a smart mathematical shortcut to film a quantum crash. They found that while some particles bounce off politely, others get so tangled that they lose their individual identity, creating a new, messy state of matter that only exists in the quantum world.

Technical Summary

1. Problem Statement

The paper addresses a central challenge in quantum field theory (QFT): understanding the real-time, non-perturbative dynamics of hadronic scattering.

• Limitations of Traditional Methods: Perturbation theory fails at strong coupling, and Euclidean lattice QCD suffers from the "sign problem," which prevents the direct simulation of real-time evolution.
• The Gap: While tensor network methods have successfully simulated scattering in Abelian ($U(1)$) gauge theories (like the Schwinger model), extending these simulations to non-Abelian models (specifically $SU(2)$) is crucial for capturing genuinely QCD-like phenomena such as confinement, baryon formation, and flavor dynamics.
• Specific Goal: To perform the first real-time study of hadronic scattering in a $(1+1)$-dimensional $SU(2)$ lattice gauge theory with fundamental fermions, focusing on sectors with different global baryon numbers ($B$).

2. Methodology

The authors employ Tensor Network (TN) techniques, specifically Matrix Product States (MPS), within a gaugeless Hamiltonian formulation.

A. Theoretical Framework

• Model: $(1+1)$D $SU(2)$ gauge theory coupled to fundamental fermions (the simplest non-Abelian generalization of the Schwinger model).
• Hamiltonian Formulation:
  • Starts with the Kogut–Susskind staggered-fermion lattice Hamiltonian in temporal gauge ($A^a_0 = 0$).
  • Gauge Elimination: Unlike $U(1)$ models, the $SU(2)$ gauge field acts as a non-dynamical potential. By exploiting residual gauge freedom and performing site-dependent color rotations, the gauge links are integrated out. This yields a purely fermionic Hamiltonian with long-range interactions mediated by electric flux, requiring no truncation of the electric flux Hilbert space.
  • Gauss's Law: Physical states are constrained by $G^a_n |\text{phys}\rangle = 0$, where $G^a_n$ involves the electric field and color charge operators.
• Spin Mapping:
  • A Jordan-Wigner transformation maps the two-component staggered fermions at each site to a spin-1/2 chain.
  • The resulting spin Hamiltonian ($H_s$) includes kinetic energy, mass terms, and long-range color-exchange interactions ($H_E$).
  • Qubit Encoding: The physical content is determined by site parity ($n \mod 4$): sites $n \equiv 1, 2$ encode antiquarks ($\bar{r}, \bar{g}$), while $n \equiv 3, 0$ encode quarks ($r, g$).

B. Simulation Setup

• Software: Numerical simulations were performed using the iTensor library.
• Parameters:
  • System size: $N = 60$ qubits.
  • Coupling: Strong coupling regime ($ga = 5$, $ma = 0.2$).
  • Algorithms:
    • Ground State Preparation: Density Matrix Renormalization Group (DMRG) with 2000 sweeps and maximum bond dimension $\chi_{\max} = 80$ (checked up to 200).
    • Real-Time Evolution: Time-Dependent Variational Principle (TDVP) with $\chi_{\max} = 80$ and time step $a\Delta t = 0.1$.
• Initial State: Hadronic wavepackets were prepared as Gaussian-modulated superpositions of local creation operators acting on the vacuum $|\Omega\rangle$. Two wavepackets with opposite momenta were initialized on opposite sides of the chain.

3. Key Contributions

1. First Non-Abelian Scattering Simulation: This work presents the first real-time tensor network simulation of scattering in a non-Abelian ($SU(2)$) gauge theory, moving beyond the previously studied Abelian ($U(1)$) models.
2. Gaugeless Formulation Implementation: Successfully demonstrated the integration of gauge fields in $(1+1)$D $SU(2)$ to create a purely fermionic spin Hamiltonian, avoiding the need for electric flux truncation.
3. Novel Observables: Introduced the Information Lattice as a diagnostic tool alongside standard entanglement entropy to characterize correlations and distinguish between bound states and collective entangled states.

4. Results

The study investigated three scattering sectors based on the total global baryon number $B$:

Sector $B=0$ (Meson–Meson)

• Dynamics: The scattering is predominantly elastic.
• Behavior: Two mesons pass through each other with their internal structures unchanged.
• Entanglement: The entanglement entropy returns to its initial profile after the collision.
• Information Lattice: Shows a characteristic transient peak at correlation length $\ell \approx 3-4$ during overlap, which restores afterward. This resembles the dynamics of the $U(1)$ Schwinger model.

Sector $B=1$ (Meson–Baryon)

• Dynamics: Reveals qualitatively new, non-Abelian physics.
• Behavior:
  • The baryon remains approximately localized.
  • The meson undergoes partial reflection and partial tunneling.
  • Crucial Finding: Post-collision, the meson probability density spreads across the lattice. The meson and baryon wavepackets do not separate cleanly; instead, they become entangled and form a single collective state.
• Information Lattice: No new distinct length scales ($\ell$) are significantly populated, indicating the collective state does not form a new bound state with a distinct internal correlation structure, but rather a delocalized entangled system.

Sector $B=2$ (Baryon–Baryon)

• Dynamics: The scattering is elastic, closely resembling the $U(1)$ Schwinger model.
• Reason: Inelastic channels (e.g., producing two mesons from a baryon-antibaryon pair) are kinematically blocked because the maximum baryon momentum is far below the threshold required for such production.

5. Significance

• Validation of QIS Tools: The study demonstrates that Quantum Information Science (QIS) inspired tools (tensor networks, information lattices) are effective for probing complex hadronic dynamics in non-Abelian gauge theories.
• Insight into Confinement: The results highlight the unique role of baryon number conservation in non-Abelian theories, showing how it leads to entangled collective states in mixed sectors ($B=1$) that are absent in Abelian models.
• Pathway to QCD: This work serves as a critical stepping stone toward simulating full QCD-like phenomena on quantum computers and classical tensor network simulators, bridging the gap between simplified models and realistic high-energy physics.

In conclusion, the paper successfully simulates real-time hadronic scattering in $SU(2)$, identifying a distinct non-Abelian regime where meson-baryon collisions result in entangled, delocalized collective states, a phenomenon not observed in Abelian analogues.