⚛️ lattice

An improved partial-wave projection of the one-particle exchange in relativistic three-body scattering

This paper presents an improved, finite-sum expression for the partial-wave projection of the one-particle exchange process in relativistic three-body scattering of massive, spinless particles, which accurately reproduces known physical properties while significantly enhancing computational efficiency for constructing robust amplitude analyses.

Original authors: Nicholas C. Chambers, Andrew W. Jackura

Published 2026-09-10

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Nicholas C. Chambers, Andrew W. Jackura

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the subatomic world, matter is not static; it is a constant, frenetic exchange of energy and particles. When physicists study how three particles scatter off one another, they are trying to understand the fundamental forces that hold the universe together, particularly the strong force that binds protons and neutrons. These interactions are incredibly complex because the particles do not just bounce off each other like billiard balls; they can temporarily swap partners, creating fleeting intermediate states before separating again. To make sense of the chaotic data from high-energy experiments, scientists break these interactions down into simpler components called partial waves. Think of these as sorting the chaotic noise of a storm into distinct, manageable frequencies, allowing researchers to isolate specific patterns and identify the short-lived particles, or resonances, that appear and disappear in the process. However, one specific mechanism—the exchange of a single particle between two pairs of colliding particles—has long been a mathematical stumbling block, creating a tangled web of calculations that becomes nearly impossible to untangle as the complexity of the collision increases.

Nicholas C. Chambers and Andrew W. Jackura have now cleared this path by developing a streamlined method to calculate this specific exchange process. In their work, they focused on a scenario where three massive, spinless particles interact, a simplified model that captures the essential kinematics of more complex real-world collisions. The core difficulty they addressed was how to describe the motion of these particles when they are viewed from different perspectives. In a three-particle collision, the particles are constantly changing their reference frames as they move and interact. To calculate the probability of a specific outcome, scientists must translate the motion of a particle from the perspective of one pair to the perspective of the entire system. Previously, this translation required a unique, custom-made calculation for every single combination of particle spins and orbital paths. As researchers tried to include more and more of these combinations to get a complete picture, the number of required calculations exploded, making it slow and prone to human error.

Chambers and Jackura solved this by finding a universal formula that replaces the need for endless custom calculations. They derived a compact expression that works for any total angular momentum and any arrangement of the particles' internal spins. Instead of building a new mathematical structure for every new scenario, their method uses a set of reusable building blocks, or coefficients, that can be generated recursively. This means that once the basic rules are established, the complex math of translating between different moving frames can be handled by a straightforward, finite sum. The researchers verified that their new formula correctly reproduces the known physical behaviors of these interactions, including how the particles behave at the very edge of their energy limits and how the mathematical singularities—points where the physics becomes extreme—appear and disappear. They confirmed that their approach captures the same critical features as previous, more cumbersome methods but does so with a clarity that allows for much larger and more robust sets of data to be analyzed.

The impact of this improvement is immediate for those studying the spectrum of hadrons, the family of particles that includes protons and neutrons. By removing the computational bottleneck, the new formula allows scientists to include a much wider range of particle interactions in their models without fear of the calculations becoming unmanageable. This is crucial because, as the energy of the collisions increases, the number of possible ways the particles can interact grows rapidly. In the past, researchers often had to limit their studies to a small, simplified set of interactions, potentially missing subtle signals of new particles. With this new tool, they can now construct comprehensive models that include higher-spin resonances and more complex orbital configurations. The authors demonstrated that their method works effectively for systems with total angular momentum up to six, a range that covers many of the particles of interest in modern hadron spectroscopy.

This work does not claim to have discovered a new particle or solved the mystery of the strong force, but rather provides a more reliable and efficient engine for the search. It is a foundational upgrade to the mathematical toolkit used by experimentalists and theorists alike. By ensuring that the description of the single-particle exchange is both accurate and computationally feasible, the researchers have removed a significant barrier to understanding the three-body dynamics that underpin much of nuclear physics. The result is a clearer view of the subatomic landscape, where the signals of rare and exotic resonances are no longer obscured by the limitations of the calculation methods used to find them. As the field moves toward analyzing even more complex systems, this streamlined approach ensures that the theoretical framework can keep pace with the growing precision of experimental data.

Technical Summary: An Improved Partial-Wave Projection of the One-Particle Exchange in Relativistic Three-Body Scattering

Problem Statement
Relativistic three-body scattering amplitudes are essential for interpreting nuclear and hadronic phenomena, particularly resonances decaying into three-particle final states. A critical component of these amplitudes is the one-particle exchange (OPE) process, where a particle is exchanged between incoming and outgoing two-body subsystems. The OPE introduces a kinematic singularity (a $u$ -channel pole) that must be accurately projected into partial waves to construct physical amplitudes consistent with unitarity and analyticity.

Previous methods for computing the partial-wave projection of the OPE (specifically the work in Ref. [63]) expressed the result in terms of coefficients that must be determined separately for every target total angular momentum and parity ( $J^P$ ) and every combination of internal angular momenta. This "case-by-case" approach becomes increasingly inefficient and error-prone as the partial-wave truncation is relaxed to include higher angular momenta, which is necessary for studying tensor and higher-spin resonances (e.g., $a_2(1320)$ ) and for reducing systematic uncertainties in lattice QCD amplitude analyses.

Methodology
The authors present an improved formulation that derives a compact, universal finite-sum representation for the partial-wave OPE, valid for three massive, spinless particles in arbitrary angular momentum configurations. The derivation proceeds in two main stages:

Frame Transformation and Boost Coefficients: The central challenge is relating angular functions defined in different rest frames (the total three-body center-of-momentum frame vs. the pair rest frames). The authors explicitly expand the Lorentz-boosted angular dependence of the two-body subsystems in terms of orthogonal functions defined in the three-body rest frame. This is achieved through "boost coefficients" ( $B^\ell_{S\lambda}$ ), which are defined as angular overlaps of spherical harmonics evaluated in different frames.
- The authors derive a recurrence relation (Eq. 11) and an explicit product representation (Eq. 41) for these coefficients, allowing them to be generated recursively from elementary rank-one coefficients.
- This approach avoids the need for separate angular integrations for each spin-orbit channel.
Partial-Wave Projection: Using the Neumann expansion for the pole term and the addition theorem for Wigner $D$ -matrices, the authors project the OPE function onto definite total angular momentum $J$ and parity $P$ . By utilizing standard properties of spherical harmonics and $SU(2)$ couplings (Clebsch-Gordan coefficients), the projection is reduced to a finite sum over an auxiliary angular momentum index $j$ .

Key Contributions
The primary result is a finite-sum expression for the spin-orbit partial-wave matrix elements of the OPE, $G^{JP}_{L'S', LS}(p, k)$ , given by Eq. (8):
$G^{JP}_{L'S', LS}(p, k) = \frac{1}{2pk} \sum_{j=0}^{j_{\text{max}}} (2j+1) Q_j(\zeta_{pk} + i\epsilon) H^{JP}_{j, L'S', LS}(p, k)$
where:

$Q_j$ are Legendre functions of the second kind, which capture the logarithmic singularities associated with the exchange pole.
$H^{JP}$ is a matrix containing the remaining angular-momentum dependence, constructed from the boost coefficients and Clebsch-Gordan coefficients (Eq. 9).
The sum terminates at $j_{\text{max}} = J + S + S'$ , ensuring the expression is a finite sum rather than an infinite series.

This formulation replaces the previous case-by-case construction with a single, implementation-ready formula valid for arbitrary $J^P$ and allowed spin-orbit transitions.

Results and Verification
The authors validate the new formulation through several checks:

Analytic Checks: The derived finite-sum expression is shown to reproduce the known threshold scaling behavior ( $G \sim p^{L'}k^L$ ) and the singularity structure of the OPE. The threshold behavior is verified analytically in Sec. III C by expanding the boost coefficients and solid harmonics in the limit of vanishing momenta.
Numerical Verification:
- The boost coefficients are shown to behave correctly in the limit where no boost is required ( $k \to 0$ ), reducing to identity mappings.
- The partial-wave OPE functions are plotted for various $J^P$ sectors (including $1^+, 1^-, 2^+, 3^-$ ) and higher angular momentum channels (up to $S=2$ and $L=3$ ). These plots confirm the expected branch point structures and unitarity constraints (non-zero imaginary parts in the physical region).
- The authors demonstrate that the finite-sum representation reproduces all formulae reported in Ref. [63] and the tabulated projections in Ref. [72] without performing separate integrations for each channel.
Truncation Analysis: Numerical studies of the relative magnitudes of subchannel contributions (Sec. IV C) reveal that there is no straightforward power-counting hierarchy for three-body subchannels across multiple kinematic variables. Higher-order partial waves (e.g., $1D_2 \to 3F_2$ ) can contribute significantly compared to lower-order diagonal terms, highlighting the necessity of robust, large wavesets in amplitude analyses.

Significance
The paper claims that this improved formulation significantly enhances the efficiency and robustness of three-body amplitude analyses. By providing a universal finite-sum representation, it eliminates the need for separate coefficient calculations for each spin-orbit channel, thereby reducing the risk of transcription errors and computational overhead. This capability is critical for:

Constructing robust wavesets for lattice QCD studies of three-hadron amplitudes, where systematic uncertainties must be quantified by varying truncation levels.
Extending analyses to higher total angular momenta ( $J$ ) required for the spectroscopy of tensor and higher-spin resonances.
Facilitating phenomenological studies of experimental data where large wavesets are necessary to disentangle resonance signals from kinematic singularities.

Although derived for spinless particles, the authors note that the underlying projection strategy provides a useful starting point for generalizations to particles with non-zero intrinsic spin.

More like this