⚛️ nuclear theory

A Practical Partial-Wave Method for Identifying Unstable Light Nucleus Resonances in Heavy-Ion Collisions

This paper proposes a practical partial-wave method based on the Lednický–Lyuboshitz framework, which incorporates experimental phase-shift data to extract resonance signals of unstable light nuclei from two-particle correlation functions in heavy-ion collisions, successfully demonstrating its application to the $^4$ Li and $^5$ Li systems using STAR experimental data.

Original authors: Junlin Wu, Hongchan Li, Ke Mi, Yaping Wang, Guannan Xie

Published 2026-09-04

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Junlin Wu, Hongchan Li, Ke Mi, Yaping Wang, Guannan Xie

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 extreme heat of a heavy-ion collision, where atomic nuclei smash together at nearly the speed of light, a fleeting soup of matter is born. This environment is so intense that it mimics the conditions of the universe just moments after the Big Bang. Within this searing chaos, scientists look for clues about how the fundamental building blocks of matter—protons and neutrons—stick together to form atomic nuclei. While stable nuclei like helium or carbon are well understood, the universe also contains unstable, short-lived versions of these particles. These unstable forms exist only as fleeting resonances, appearing for a fraction of a second before breaking apart. Detecting them is like trying to photograph a ghost; they leave no direct trace, only a subtle ripple in the way other particles scatter around them. Understanding these fleeting states is crucial because they reveal the forces that hold matter together and how the universe cooled down to form the atoms we see today.

For decades, researchers have studied stable light nuclei produced in these collisions, but the unstable ones have remained largely a mystery. The challenge is that these unstable particles decay so quickly they cannot be caught directly. Instead, scientists must infer their existence by looking at the patterns left behind by the particles they leave behind. When two particles, such as a proton and a helium nucleus, fly out of a collision, their paths are influenced by the forces between them. If an unstable resonance exists, it creates a specific, tell-tale bump in the way these particles correlate with one another. However, this signal is often buried under a sea of background noise from other types of interactions, making it incredibly difficult to isolate.

A team of researchers has now proposed a practical new method to pull these hidden signals out of the noise. They developed a technique that acts like a sophisticated filter, capable of separating the specific signature of an unstable resonance from the general background of particle interactions. By applying a mathematical framework originally designed for simple interactions and expanding it to handle more complex, spinning interactions, the team created a tool that can identify these short-lived states without needing to guess the details of the forces involved. Instead of building a theoretical model of the nuclear force, which can be uncertain, they used real-world data from low-energy scattering experiments as their guide. This approach allows them to map out the exact shape of the resonance signal, distinguishing it from the non-resonant background that usually obscures it.

The researchers tested their method on two specific cases: a system involving a proton and a helium-3 nucleus, and another with a proton and a helium-4 nucleus. These combinations correspond to the unstable ground states of lithium-4 and lithium-5, respectively. When they ran their calculations, the method worked exactly as hoped. For the lithium-4 system, a clear peak emerged in the correlation data at a relative momentum of approximately 72 MeV/c. For the lithium-5 system, a similar peak appeared at about 50 MeV/c. These positions match the known energy levels of these unstable nuclei, confirming that the method successfully isolates the resonance signal from the surrounding clutter. The results show that the resonance creates a distinct excess in the number of particle pairs detected at these specific momenta, a feature that would be impossible to see if one only looked at the total data without this specialized decomposition.

To show how this tool could be used in real experiments, the team applied their method to data from the STAR experiment at the Relativistic Heavy Ion Collider, which operates at a collision energy of 3 GeV. They used the measured spectra of protons and light nuclei from these collisions to predict what the yields of these unstable lithium isotopes would look like. The simulations produced detailed maps of how these particles would be distributed in momentum and speed across different collision scenarios. The results indicated that the production of these unstable nuclei is directly linked to the strength of the resonance, with the more unstable states appearing in proportion to their resonance characteristics. While these specific numbers are based on a demonstration with fixed parameters rather than a final measurement, the study proves that the underlying physics is sound and that the signal can be cleanly extracted.

This work offers a new, practical path for experimentalists to study the most elusive members of the nuclear family. By moving beyond simple models and using a method that respects the complex quantum nature of these interactions, scientists can now hope to identify a wider range of unstable nuclei, such as helium-6 or beryllium-8, in future collision experiments. The ability to separate the resonant signal from the background means that researchers can finally measure the properties of these short-lived states with greater precision. This could lead to a deeper understanding of how matter behaves under extreme conditions and how the forces of nature shape the building blocks of our world. The method provides a clear, reliable way to look for the ghosts in the machine, turning what was once a theoretical challenge into a manageable experimental task.

Technical Summary: A Practical Partial-Wave Method for Identifying Unstable Light Nucleus Resonances in Heavy-Ion Collisions

Problem Statement
The production of light nuclei in relativistic heavy-ion collisions (HICs) serves as a probe for the dynamics of hot, dense matter and freeze-out conditions. While stable light nuclei (e.g., deuterons, tritons, $^3$ He, $^4$ He) have been extensively studied using coalescence and statistical thermal models, unstable light nuclei remain largely unexplored. These unstable nuclei exist as short-lived resonance states (e.g., $^4$ Li and $^5$ Li) and decay via proton emission. Extracting their signals from invariant mass distributions is challenging because the resonance signals are often mixed with contributions from non-resonant partial waves and background correlations. Traditional methods, such as the standard Lednický–Lyuboshitz (LL) model, are primarily formulated for s-wave interactions and are insufficient for isolating resonant signals in higher-order partial waves (e.g., p-waves) where these unstable nuclei reside. Furthermore, previous attempts to model these correlations often relied on model-dependent potential parametrizations or failed to fully incorporate Coulomb interactions, which are critical for charged systems like p- $^3$ He and p- $^4$ He.

Methodology
The authors propose a partial-wave method based on an extension of the Lednický–Lyuboshitz (LL) framework to extract resonance signals from two-particle correlation functions. The theoretical approach involves the following key steps:

Extension of the LL Model: The standard LL model, which assumes s-wave dominance, is extended to arbitrary partial waves ( $l > 0$ ). The correlation function $C(k)$ is calculated using the Koonin-Pratt equation, integrating the squared asymptotic wave function $|\psi(\vec{r}, \vec{k})|^2$ over a Gaussian emission source $S(r)$ .
Asymptotic Wave Function with Coulomb Interaction: The method employs an asymptotic wave function that accounts for both short-range nuclear forces and long-range Coulomb interactions. To handle the divergence of the wave function at $r \to 0$ for higher partial waves, a cutoff parameter $r_c$ (set to 1 fm) is introduced, following a regularization scheme. The wave function is defined as the incoming wave for $r < r_c$ and the sum of incoming and scattered outgoing waves for $r > r_c$ .
Input from Experimental Phase Shifts: Instead of constructing parameterized nuclear potentials, the method directly incorporates experimental phase-shift data from low-energy nuclear scattering experiments (specifically p- $^3$ He and p- $^4$ He scattering). This allows for a model-independent determination of the scattering amplitudes.
Renormalized Scattering Amplitude: To address singularities in the scattering amplitude at low momentum ( $k \to 0$ ) in the presence of Coulomb interactions, a renormalized scattering amplitude is utilized. This involves an effective range expansion for the s-wave and direct interpolation of phase-shift data for p-waves.
Spin Weighting and Decomposition: The total correlation function is constructed as a spin-weighted sum of contributions from all relevant partial waves (s-wave and p-wave). The authors specifically isolate the resonant partial wave (e.g., $^3P_2$ for $^4$ Li and $P_{3/2}$ for $^5$ Li) from non-resonant background contributions (other s- and p-waves).
Numerical Implementation: A Monte Carlo sampling procedure is used to account for the momentum-space phase-space distribution and Lorentz contraction of the emission source. The method is applied to the p- $^3$ He and p- $^4$ He systems to simulate the production of $^4$ Li and $^5$ Li, respectively.

Key Contributions

Generalized LL Framework: The paper extends the LL model to higher-order partial waves while rigorously including Coulomb interactions, a necessary condition for analyzing charged light-nucleus systems.
Model-Independent Approach: By using experimental phase-shift data directly as input, the method avoids the need for detailed and often uncertain nuclear potential parametrizations.
Resonance Isolation: The approach enables the clean decomposition of resonant partial waves from non-resonant backgrounds, addressing the "snowball in the furnace" paradox and the difficulty of identifying short-lived resonances in invariant mass spectra.
Practical Tool for HICs: The authors provide a concrete framework for experimental collaborations to extract resonance signals from existing correlation function data.

Results
The method was applied to the p- $^3$ He and p- $^4$ He systems to demonstrate the extraction of $^4$ Li and $^5$ Li resonances:

Correlation Functions: Numerical results show a distinct peak in the correlation function corresponding to the resonance energy. For $^4$ Li (p- $^3$ He), a peak appears at $k \approx 72$ MeV/c. For $^5$ Li (p- $^4$ He), a peak appears at $k \approx 50$ MeV/c. These positions are consistent with known resonance parameters.
Background Subtraction: The study demonstrates that without subtracting the non-resonant background (other partial waves), the extracted signal in the low-momentum region could be negative due to repulsive Coulomb and s-wave interactions. The proposed method successfully isolates the resonant excess.
Spectra and Yields: Using measured proton and light-nuclei spectra from the STAR experiment at $\sqrt{s_{NN}} = 3$ GeV, the authors present estimated transverse momentum spectra ( $p_T$ ) and rapidity distributions ($dN/dy$) for $^4$ Li and $^5$ Li across different centrality bins. The results indicate that the resonance strength is positively correlated with the yield.
Parameter Sensitivity: The study analyzes the dependence of the correlation function on the source size ( $R_g$ ), the cutoff parameter ( $r_c$ ), and a scaling factor ( $g_p$ ) for p-wave contributions. It is noted that while $r_c$ significantly affects the magnitude of the p-wave contribution, the resonance peak position remains stable.

Significance and Claims
The paper claims to provide a "practical tool" for the experimental study of unstable light nuclei in relativistic heavy-ion collisions. The authors assert that their method allows for the effective isolation of resonance-induced correlation excesses, which is crucial for probing final-state interactions and freeze-out conditions. By extending the applicability of femtoscopic techniques to higher partial waves and charged systems, the method can be applied to a broader range of resonance states, including $^6$ He, $^6$ Li, $^7$ Li, and $^8$ Be. The authors emphasize that this approach helps resolve the challenges in measuring unstable nuclei, which have historically been difficult to detect due to their short lifetimes and the mixing of partial waves in invariant mass distributions. The work is presented as a demonstration of feasibility, using fixed parameters to show the potential of the method rather than making definitive quantitative predictions for specific experimental yields.

Technical Summary: A Practical Partial-Wave Method for Identifying Unstable Light Nucleus Resonances in Heavy-Ion Collisions

More like this