Real-time Scattering in \phi^4 Theory using Matrix Product States

This paper demonstrates that uniform matrix product states combined with the time-dependent variational principle can effectively simulate real-time two-particle scattering in (1+1)-dimensional $\phi^4$ theory, revealing distinct elastic and inelastic behaviors across different phases and identifying a dynamical divergence near the critical point as a signature of the closing mass gap.

The Gist

The Big Picture: Simulating the Universe on a Computer

Imagine you are a physicist trying to understand how tiny particles interact. In the real world, we smash particles together in giant machines like the Large Hadron Collider (LHC) to see what happens. But sometimes, the math is so incredibly complex (involving infinite possibilities and "strong" forces) that our best supercomputers can't solve the equations directly.

This paper is about a new way to simulate these particle collisions using a clever mathematical trick called Matrix Product States (MPS). Think of this not as a standard computer simulation, but as building a "digital model" of the universe that is smart enough to ignore the noise and focus only on the most important connections between particles.

The authors are studying a specific type of particle interaction (called $\phi^4$ theory) in a simplified 1-dimensional world (like a single string of beads). They wanted to answer two big questions:

1. Where is the tipping point? (At what point does the system change from one state to another?)
2. What happens when particles crash into each other? (Do they bounce off cleanly, or do they shatter into new pieces?)

---

Analogy 1: The "Sandwich" Protocol

To study collisions, the researchers used a method they call a "Sandwich Geometry."

Imagine a long, infinite loaf of bread (representing the empty vacuum of space). In the middle of this loaf, you place a small, distinct filling (the two particles you want to collide).

• The Bread: This is the "ground state" or the calm, empty background.
• The Filling: These are two little "wave packets" (like ripples in a pond) moving toward each other from the left and right.

The researchers set up this digital sandwich and then hit "play" on the simulation. They watched what happened when the two ripples met in the middle. Did they bounce off? Did they merge? Did they explode?

Analogy 2: The "Traffic Jam" vs. The "Ghost"

The paper looks at three different "neighborhoods" or phases of this particle world. Here is how the collisions behaved in each:

1. The Symmetric Phase (The Chaotic City)

• The Setting: Imagine a busy, chaotic city where the rules of physics are "soft" and fluctuating.
• The Crash: When the two particles collided here, it was a messy crash. Instead of bouncing off cleanly, they smashed together and created a lot of debris (new particles).
• The Result: The simulation showed that only about 71% of the time did the particles bounce off as they started. The rest of the time, they turned into something else. It was highly "inelastic."

2. The Broken Phase (The Calm Countryside)

• The Setting: Imagine a quiet, frozen lake. The rules here are rigid and stable.
• The Crash: When the particles collided here, they acted like ghosts or perfect billiard balls. They passed through each other or bounced off with almost zero energy loss.
• The Result: The collision was 99.9% elastic. They didn't break apart; they just kept going. This tells us that in this phase, the particles are very stable and don't like to change into other things.

3. The Critical Point (The Edge of the Cliff)

• The Setting: This is the exact moment of transition between the chaotic city and the frozen lake. It's like standing on the edge of a cliff where the ground is about to give way.
• The Crash: This is where the magic happened. When they tried to simulate a collision here, the simulation broke down.
• Why? In this critical state, the "correlation length" (how far a ripple can travel before fading) becomes infinite. It's like trying to measure the ripple of a pebble in an ocean that is infinitely deep and wide. The "sandwich" (the finite window of the simulation) was too small to contain the chaos.
• The Discovery: The fact that the simulation failed to produce a clean collision was actually a success. It proved they had found the exact critical point. The "failure" was a signature that the system had become infinitely connected and sensitive.

The Tool: "Finite-Entanglement Scaling"

How did they find the exact spot of the critical point without knowing the answer beforehand?

They used a technique called Finite-Entanglement Scaling.

• The Analogy: Imagine trying to draw a perfect circle on a piece of graph paper. If the grid is too coarse (low resolution), your circle looks like a jagged square. If you zoom in and use a finer grid (higher resolution), it looks more like a circle.
• The Science: The researchers ran their simulation with different "grid sizes" (called bond dimensions). They noticed that as they made the grid finer, the results started to converge on a specific number. By looking at how the results changed as they zoomed in, they could mathematically pinpoint the exact location of the critical point (the "tipping point") with extreme precision.

Why Does This Matter?

1. It's a New Way to Do Physics: Usually, to study particle collisions, you need massive particle accelerators or very rough approximations. This paper shows you can use "tensor networks" (a type of AI/math structure) to simulate these collisions on a regular computer with high precision.
2. It Maps the Unknown: They created a detailed map of the $\phi^4$ theory, showing exactly where the "chaotic" phase ends and the "stable" phase begins.
3. It Detects Critical Points: They proved that you can find the "edge of the cliff" (the critical point) not just by looking at static pictures, but by watching how the system behaves when you try to crash things into each other. If the crash looks weird and the simulation struggles, you've found the critical point.

Summary

The authors built a digital "sandwich" of particles and watched them crash.

• In the chaotic phase, the crash was messy and created new particles.
• In the stable phase, the crash was clean and perfect.
• At the critical point, the crash was so complex that the simulation couldn't handle it, which served as a perfect signal that they had found the exact tipping point of the universe they were modeling.

They successfully used this method to map out the behavior of a fundamental quantum field theory, proving that these mathematical tools are powerful enough to explore the deepest secrets of particle physics without needing a billion-dollar collider.

Technical Summary

1. Problem Statement

The paper addresses the challenge of simulating real-time scattering dynamics and critical behavior in interacting relativistic quantum field theories (QFTs) directly from first principles.

• Limitations of Existing Methods: Perturbative techniques (Feynman diagrams) fail in strongly coupled regimes and cannot easily access real-time dynamics. While analytic S-matrix methods exist for integrable models in 1+1 dimensions, they do not apply to non-integrable theories like the interacting $\phi^4$ model.
• The Gap: There is a need for non-perturbative methods that can handle the thermodynamic limit, locate quantum critical points (QCPs), and simulate particle collisions (scattering) to extract observables like elastic probabilities and time delays.
• Target System: The study focuses on the $(1+1)$-dimensional $\phi^4$ theory, a canonical model for studying phase transitions (Ising universality) and spontaneous symmetry breaking.

2. Methodology

The authors employ Uniform Matrix Product States (uMPS) combined with the Time-Dependent Variational Principle (TDVP).

A. Discretization and Hamiltonian

The continuous Lagrangian is discretized onto a 1D lattice with spacing $a=1$. The resulting Hamiltonian is:
$$H_{\phi^4} = \sum_j \left( \frac{1}{2}\pi_j^2 + \frac{1}{2}(\phi_{j+1} - \phi_j)^2 + \frac{1}{2}\mu^2\phi_j^2 + \frac{\lambda}{4!}\phi_j^4 \right)$$
where $\mu^2$ is the bare mass-squared and $\lambda$ is the quartic coupling. The local Hilbert space is truncated to dimension $d$, and the system is renormalized by subtracting the local ground-state energy density.

B. Finite-Entanglement Scaling (FES) for Criticality

To locate the quantum critical point ($\mu_c^2$) without finite-size scaling:

• They utilize the relationship between the bond dimension $D$ and the effective correlation length $\xi_D$ near a conformal fixed point.
• The von Neumann entropy $S$ scales as $S \simeq \frac{c}{6} \log \xi_D + \text{const}$, where $c$ is the central charge.
• By sweeping $\mu^2$ and varying $D$, they extract the effective central charge $c_{\text{eff}}$. The critical point is identified where $c_{\text{eff}} \to 0.5$ (Ising universality).

C. Real-Time Scattering Protocol ("Sandwich Geometry")

To simulate scattering in the thermodynamic limit:

1. Vacuum Preparation: A uniform uMPS ground state is computed for a specific phase (symmetric, broken, or near-critical).
2. Wave Packet Imprinting: Two localized wave packets with opposite momenta ($\pm \kappa$) are imprinted onto the uniform background using a "sandwich" geometry. This involves creating a non-uniform MPS segment embedded in the uniform vacuum.
3. TDVP Evolution: The state is evolved in real-time using TDVP, which projects the Schrödinger equation onto the tangent space of the uMPS manifold, ensuring unitary evolution within the variational ansatz.
4. Observables Extraction:
   • Elastic Scattering Probability ($P_{11\to11}$): Calculated by comparing the spectral weight of the diagonal momentum components before and after collision.
   • Wigner Time Delay ($\Delta t$): Extracted from the phase shift of the elastic S-matrix element $S_{11\to11}(E)$ using finite differences in momentum.

3. Key Contributions

• Extension of TDVP-uMPS: Successfully extends scattering protocols (previously applied to Ising field theory) to the interacting, non-integrable $\phi^4$ theory.
• Precision Critical Point Mapping: Provides a quantitatively controlled estimate of the critical mass-squared $\mu_c^2$ using FES, avoiding ad-hoc parameter choices.
• Dynamical Signature of Criticality: Demonstrates that the breakdown of the "sandwich" scattering protocol at the critical point is not a numerical failure but a physical signature of the diverging correlation length.
• Phase-Dependent Scattering Dynamics: Systematically maps how scattering behavior transitions from highly inelastic (symmetric phase) to nearly perfectly elastic (broken phase).

4. Key Results

A. Critical Point Determination

Using FES at $\lambda = 0.8$ with bond dimensions $D \in \{36, 48, 64\}$:

• The effective central charge $c_{\text{eff}}$ converges to $0.5$ near $\mu^2 \approx -0.25945$.
• The critical mass-squared is bounded as: $\mu_c^2 \in ]-0.2595, -0.2594[$.
• This confirms the system belongs to the 2D Ising universality class.

B. Excitation Spectra

• Symmetric Phase ($\mu^2 = +0.2$): Gapped spectrum with a clear separation between single-particle and multi-particle bands.
• Critical Point ($\mu^2 \approx -0.2594$): The mass gap closes ($m=0$), and the spectrum becomes gapless with linear dispersion ($E \propto |p|$), consistent with Conformal Field Theory (CFT).
• Broken Phase ($\mu^2 < 0$): The gap reopens, and a massive single-particle band emerges.

C. Scattering Dynamics

The study simulated two-particle collisions across different phases:

1. Symmetric Phase ($\mu^2 = +0.2$):
   • Result: Strongly inelastic scattering.
   • Metrics: $P_{11\to11} \approx 0.712$, Time delay $\Delta t \approx -158$.
   • Interpretation: Significant particle production and radiation due to the non-integrable nature of the theory in this regime.
2. Spontaneously Broken Phase ($\mu^2 = -0.1$ and $-0.5$):
   • Result: Almost perfectly elastic scattering.
   • Metrics: $P_{11\to11} \approx 1$.
   • Interpretation: Excitations behave as stable, weakly radiating quasiparticles. The time delay becomes more negative (e.g., $\Delta t \approx -177.78$ for $\mu^2=-0.5$) as the gap increases.
3. Near Critical Point ($\mu^2 \approx \mu_c^2$):
   • Result: Breakdown of the scattering protocol.
   • Observation: Instead of a clean "X" collision pattern, the field expectation value $\langle \phi_n(t) \rangle$ shows a long-wavelength drift across the entire window. No distinct incoming/outgoing packets separate.
   • Cause: The assumption of a finite correlation length (required to isolate the scattering region from the asymptotic vacuum) fails because $\xi \to \infty$ at criticality. This serves as a dynamical signature of the quantum critical point.

5. Significance and Conclusion

• Non-Perturbative Benchmarking: The work demonstrates that uMPS-TDVP can effectively probe non-perturbative scattering and critical dynamics in lattice field theories with controlled entanglement truncation.
• Unified Framework: It unifies static criticality analysis (via FES) and dynamic scattering analysis within a single numerical framework.
• Physical Insight: The results highlight a fundamental difference in quasiparticle stability: the symmetric phase supports rich inelastic channels, while the broken phase supports stable soliton-like or particle-like excitations.
• Future Directions: The authors suggest that increasing bond dimensions and refining wave-packet constructions will allow for more precise phase shift calculations and the study of more complex QFTs. The "failure" of the sandwich geometry at criticality is identified not as a bug, but as a powerful diagnostic tool for detecting quantum phase transitions dynamically.