Higher order quantization conditions for two-body scattering with spin

This paper derives and validates high-order Lüscher quantization conditions for spinless and spin-1/2 particle scattering in both cubic and elongated periodic boxes, providing a robust framework for precision finite-volume studies of systems with half-integer spin, such as meson-baryon scattering.

The Gist

Imagine you are trying to understand how two billiard balls bounce off each other. In the real, infinite world, you could roll them across a giant table, watch them collide, and measure exactly how they scatter. This is what physicists call "scattering," and it tells us about the forces holding the universe together.

However, in the world of Quantum Chromodynamics (QCD)—the theory that explains how quarks and gluons stick together to form protons and neutrons (baryons)—we can't just roll infinite billiard balls. We have to simulate the universe on a computer.

The Problem: The Tiny Room

To simulate this, physicists put the particles inside a tiny, invisible, cubic "room" (a lattice) with walls that wrap around. If a particle hits the wall, it instantly reappears on the opposite side, like in the video game Pac-Man.

Because the room is so small, the particles can't move freely. They are forced to dance in specific, quantized patterns, like a guitar string that can only vibrate at specific notes. This creates a list of allowed energy levels.

The Challenge: We want to know how these particles scatter in the real, infinite world, but we only have their energy levels from the tiny, fake room. How do we translate the "notes" from the tiny room to the "scattering behavior" of the real world?

The Solution: The Lüscher Map

In 1991, a physicist named Martin Lüscher invented a mathematical "dictionary" (called a Quantization Condition) that translates the energy levels in the box to the scattering data in the infinite world.

Think of it like this: If you know the exact pitch of a note a bell makes when it's inside a small box, Lüscher's formula tells you exactly how that bell would ring if it were in a vast open field.

What This Paper Does: Adding Spin and Complexity

Previous versions of this "dictionary" worked great for simple particles (like ping-pong balls with no spin). But nature is more complicated. Many particles, like protons and neutrons, have spin (they are like spinning tops).

This paper, by Lucas Chandler, Frank X. Lee, and Andrei Alexandru, does three major things:

1. Updates the Dictionary for Spinning Tops: They rewrote the mathematical rules to handle particles that spin (specifically, a spin-0 particle hitting a spin-1/2 particle). This is crucial for studying meson-baryon scattering (like a pion hitting a proton), which is essential for understanding nuclear physics.
2. Expands the Room Shapes: They didn't just look at cubic rooms. They looked at elongated rooms (like a shoebox) and moving rooms (where the whole box is zooming through space). This gives physicists more flexibility to find the specific "notes" they need.
3. Checks the Math (The "Stress Test"): This is the most important part. They didn't just write the formulas; they built a simulated universe to test them.
   • They created a fake potential (a made-up force) between two particles.
   • They calculated the energy levels in the box using one method (solving the Schrödinger equation).
   • They calculated the infinite-world scattering using a different method.
   • They fed the infinite-world data into their new formulas to see if they could predict the box energy levels.
   • The Result: The formulas worked perfectly, matching the simulation to six decimal places.

The "Spin-Orbit" Twist

The paper also deals with spin-orbit coupling. Imagine two dancers: one is spinning, and the other is moving in a circle. Their spins interact, changing how they dance together. This makes the math much harder because the "notes" in the box get mixed up. The authors showed how to untangle this mess, ensuring that even with these complex interactions, the translation from the box to the real world remains accurate.

Why This Matters

This work is like building a high-precision bridge.

• On one side: We have the raw data from supercomputers (the energy levels in the box).
• On the other side: We have the real physics of the universe (how protons and pions interact).
• The Bridge: These new, high-order quantization conditions.

Without this bridge, the data from supercomputers is just a list of numbers with no connection to reality. With this bridge, physicists can finally extract precise details about how matter is built, potentially leading to breakthroughs in understanding nuclear forces and the structure of the universe.

In short: The authors took a complex mathematical tool, upgraded it to handle spinning particles and weird room shapes, and then proved it works by running a massive simulation. They have handed the physics community a sharper, more reliable tool to decode the secrets of the atomic world.

Technical Summary

1. Problem Statement

The paper addresses the challenge of extracting infinite-volume scattering phase shifts from finite-volume energy spectra in Lattice Quantum Chromodynamics (LQCD), specifically for systems involving half-integer spin (e.g., meson-baryon scattering).

• Context: The standard Lüscher method connects finite-volume energy levels to scattering phase shifts. While well-established for spinless particles, extending this to systems with spin (specifically spin-0 and spin-1/2) is complex due to the mixing of partial waves and the necessity of using double-cover group theory.
• Limitations of Previous Work: Existing literature often covers only low-order partial waves (up to $J=7/2$) or specific geometries. There was a lack of a comprehensive, high-order framework covering:
  • Total angular momentum up to $J = 11/2$.
  • Both cubic and elongated (z-elongated) box geometries.
  • Both rest frames and moving frames (non-zero total momentum).
  • Systems with unequal masses.
• Goal: To derive, validate, and provide a robust framework for high-order quantization conditions (QCs) for spin-0/spin-1/2 scattering to enable precision studies in the baryon sector.

2. Methodology

The authors employ a multi-step approach combining non-relativistic quantum mechanics, group theory, and numerical validation.

A. Theoretical Derivation

1. Framework: The derivation is performed in a non-relativistic framework (Schrödinger equation) but is shown to be easily transformable to relativistic kinematics.
2. Formalism:
   • The scattering problem is formulated in a periodic box of size $L$ (with elongation factor $\eta$ in the $z$-direction).
   • The wavefunction is matched between the interior (interaction region) and exterior (free region) using Green's functions.
   • The resulting condition is a determinant equation involving a matrix $M(k, L)$ constructed from generalized zeta functions ($Z_{lm}$) and phase shifts ($\delta_{Jl}$).
3. Spin and Symmetry:
   • Spin-1/2 Inclusion: The formalism couples orbital angular momentum $\ell$ with spin $s=1/2$ to form total angular momentum $J = \ell \pm 1/2$. This requires double-cover group theory (e.g., $2O_h$, $2D_{4h}$) to handle the spinor nature of the wavefunctions.
   • Irreducible Representations (Irreps): The continuous rotational symmetry is broken by the box geometry. The authors project the infinite tower of partial waves onto the irreps of the specific symmetry group (cubic or elongated, rest or moving frame).
   • Moving Frames: For moving frames (total momentum $\vec{P} \neq 0$), the symmetry is reduced to "little groups" (e.g., $2C_{4v}$, $2C_{3v}$, $2C_{2v}$). The formalism accounts for the complex phase factors in the boundary conditions and the mixing of even and odd partial waves due to the loss of parity symmetry in moving frames.

B. Numerical Validation Strategy

To verify the derived QCs, the authors developed a rigorous cross-checking method:

1. Independent Energy Spectrum: They solved the Schrödinger equation numerically for a specific test potential (including spin-orbit coupling) in a finite periodic box.
   • They utilized a reduced basis ($|P, r, m_s\rangle$) to reduce the computational cost from $N^6$ to $N^3$, allowing for high-precision calculations on large lattices.
   • They employed high-order finite-difference stencils (7-point for Laplacian) to minimize discretization errors ($O(a^6)$) and performed continuum extrapolations.
2. Independent Phase Shifts: They calculated infinite-volume phase shifts for the same potential using the variable phase method.
3. Convergence Test:
   • They fed the infinite-volume phase shifts into the derived QCs to predict the finite-volume energy levels.
   • They compared these predictions against the numerically computed box energy levels.
   • Metric: A numerical $\chi^2$ measure was defined to quantify the convergence order by order (adding higher $J$ values). A QC is considered converged when the predicted levels match the box levels within machine precision.

3. Key Contributions

1. High-Order Derivations: The paper provides explicit quantization conditions up to $J = 11/2$ for 19 distinct scenarios, including:
   • Rest frames in cubic and elongated boxes.
   • Moving frames with boost vectors $\vec{d} = (0,0,1)$, $(1,1,1)$, and $(1,1,0)$.
   • Both cubic ($2O_h$) and elongated ($2D_{4h}$) geometries.
2. Spin-Orbit Coupling Formalism: A detailed derivation of how spin-orbit interactions are incorporated into the matrix elements $M_{Jl, J'l'}$ within the double-cover group framework.
3. Validation of Convergence: The study demonstrates that:
   • Truncating QCs at low orders (standard Lüscher formula) can lead to significant errors, especially in moving frames where partial wave mixing is severe.
   • Convergence requires including higher partial waves (up to $J=11/2$) to accurately reproduce energy levels, particularly for "pinched" levels near non-interacting poles.
   • Kramers Degeneracy: The authors explicitly verify that 1D irreps in moving frames (K1 and K2) form degenerate pairs due to time-reversal symmetry, a feature naturally handled by their formalism.
4. Public Resources: The authors provide a supplemental file containing Python code to reconstruct the full matrix elements for any QC and order, as well as a PDF with all convergence data, ensuring reproducibility.

4. Results

• Convergence Analysis:
  • Rest Frames: Convergence is generally faster. For example, in the $G_{1g}$ irrep, the lowest 5 levels converged by Order 3 (including $J=1/2, 7/2, 9/2$).
  • Moving Frames: Convergence is slower due to the mixing of even and odd $\ell$ (loss of parity). For the $G_2$ irrep in a moving frame, convergence required up to Order 5 ($J=11/2$) to achieve precision.
  • Pinched Levels: The study highlights that levels extremely close to non-interacting poles ("pinched" levels) require higher-order partial waves to be "unpinched" and correctly predicted.
• Elongated Boxes: The results show that elongated boxes perform comparably to cubic boxes in terms of accuracy but offer a cost-effective way to access a wider kinematic range with fewer lattice points.
• Matrix Elements: The paper presents the explicit matrix elements for the $M$ matrix in various irreps (Tables 9–15 in the Appendix), which are essential for implementing these conditions in LQCD codes.

5. Significance

• Enabling Precision Baryon Physics: This work removes a major bottleneck for LQCD studies of meson-baryon scattering (e.g., $\pi N$, $K N$). By providing validated, high-order QCs, researchers can now reliably extract phase shifts and resonance parameters (like the $\Delta(1232)$ or $\Lambda(1405)$) from finite-volume data.
• Handling Partial Wave Mixing: The paper clarifies the impact of partial wave mixing in moving frames, demonstrating that neglecting higher orders leads to unreliable phase shift extraction. This guides future LQCD simulations on which irreps and energy levels to prioritize.
• Generalizability: Although derived non-relativistically, the framework is explicitly shown to be adaptable to relativistic kinematics, making it immediately applicable to modern LQCD calculations.
• Resource for the Community: The provision of automated tools (Python code) and comprehensive tables lowers the barrier to entry for other groups wishing to apply Lüscher's method to spin-1/2 systems.

In summary, this paper establishes a rigorous, validated, and accessible foundation for the next generation of finite-volume scattering studies involving baryons, ensuring that the extraction of hadronic interaction parameters is performed with the necessary precision and theoretical completeness.