Kinetic approach of light-nuclei production in intermediate-energy heavy-ion collisions

This work presents a kinetic approach that dynamically incorporates transformations between nucleons and light nuclei as well as the Mott effect to successfully reproduce the experimental yields of light nuclei in medium-energy heavy-ion collisions, and attributes the enhancement of alpha-particle production at low energies to their high binding energy, which resists dissolution in the nuclear medium.

The Gist

Imagine a high-energy particle collision as a chaotic, super-fast game of cosmic billiard balls. Normally, physicists focus on the individual balls (protons and neutrons, or "nucleons") and the sparks they produce (pions). But in this work, the authors, led by Rui Wang and colleagues, decide to pay attention to something else: the temporary "clumps" or "teams" that form when these balls stick together. These clumps are light nuclei, such as deuterium (2 balls), tritium (3 balls), helium-3 (3 balls), and the alpha particle (4 balls sticking together).

Here is the story of their research, broken down into simple concepts:

1. The Problem: The "Teams" Are Ignored

In standard physics simulations of these collisions, scientists often treat each particle as a lone wolf. They calculate how individual balls bounce off one another. Yet, in the midst of a heavy crash (like the collision of two gold atoms), these balls often stay together to form small teams before flying apart again.

The authors argue that ignoring these teams is like watching a soccer match but only tracking individual players while disregarding the fact that they sometimes group together. To get the true picture, one must track the teams during the game, not just at the end.

2. The Solution: A New "Kinetic" Rulebook

The team developed a new set of rules (a "kinetic approach") to simulate these collisions. Imagine updating simulation software to recognize two new types of moves:

• Forming a Team: Two or more nucleons collide and stay together to become a light nucleus.
• Dissolving: A nucleon hits a light nucleus hard enough to tear it apart back into individual pieces.

They included all light nuclei up to the size of an alpha particle (4 nucleons). This allows the simulation to show how these teams are constantly created and destroyed during the crash.

3. The "Mott Effect": The Analogy of a Crowded Room

The most interesting part of their study is a phenomenon called the Mott effect.

Imagine a light nucleus (like an alpha particle) as a small group of friends holding hands in a crowded room.

• In an empty room (low density): The friends can easily hold hands and stay together.
• In a crowded room (high density): If the room is so full of other people (surrounding nucleons) that there is no room to move, the friends can no longer hold hands. They are forced to let go and drift apart as individuals.

In physical terms: When the density of the surrounding nuclear matter is too high, the "glue" holding the light nucleus together stops working, and the nucleus dissolves. The authors added a rule to their simulation: A light nucleus can only exist if the matter around it is not too dense.

4. The Alpha Particle Puzzle

The researchers compared their new simulation with real data collected by the FOPI collaboration, which collided gold atoms at various speeds.

They noticed something surprising: At lower collision speeds, there were far more alpha particles (4-nucleon teams) than expected. In fact, there were more alpha particles than helium-3 (3-nucleon teams).

Why?
The authors explain this again using the crowded room analogy.

• The alpha particle is like a very tightly bonded group of friends; they hold hands very firmly (high binding energy).
• The other light nuclei are like groups holding hands more loosely.
• When the "room" gets crowded, the loose groups let go immediately. But the tightly bonded alpha group is so strong that it can hold on even in a very crowded room.

Because the alpha particle is so robust, it survives the "Mott effect" (dissolution due to overcrowding) much better than the others. This explains why we see so many of them in the data.

5. The Result

By using their new simulation that tracks these teams and accounts for the "crowded room" rule (Mott effect), the authors successfully replicated the experimental results. They showed that the strange abundance of alpha particles is no puzzle; it simply stems from the fact that alpha particles are the "most robust" light nuclei and can survive in the dense, chaotic environment of a nuclear collision where others cannot.

In short: The work builds a better video game simulation of nuclear crashes. By allowing particles to form temporary teams and recognizing that some teams are too strong to be torn apart by the crowd, they finally solved the mystery of why so many alpha particles appear in these experiments.

Technical Summary

Technical Conclusion: Kinetic Approach to Light Nuclei Production in Intermediate-Energy Heavy-Ion Collisions

Problem Statement
Intermediate-energy heavy-ion collisions (from Fermi energy up to the GeV range) are crucial for investigating effective nucleon-nucleon interactions and the nuclear equation of state. While nucleon and pion observables (e.g., collective flow, spectral ratios) have been extensively analyzed, light nuclei ($A \le 4$), which are abundantly produced in these collisions, have historically received less theoretical attention. Most standard kinetic theory approaches, such as those based on the Boltzmann–Uehling–Uhlenbeck (BUU) equation, truncate interactions at two-body collisions and consider only nuclear degrees of freedom. Consequently, they fail to dynamically describe the production and dissociation of light nuclei (deuteron, triton, helium-3, and $\alpha$-particles) as well as the Mott effect—the phenomenon where light nuclei dissolve in a dense nuclear medium when the surrounding nucleon phase-space density becomes too high.

New experimental data from the INDRA and FOPI collaborations, particularly for central Au+Au collisions, show a significant increase in $\alpha$-particle yields at low beam energies. This observation suggests a weaker Mott effect for $\alpha$-particles compared to lighter clusters, likely due to their larger binding energy. However, existing theoretical models lack a unified dynamic framework that includes all light nuclei up to $A=4$ to explain these yields and the underlying medium effects.

Methodology
The authors develop a kinetic approach based on the real-time many-body Green's function formalism. This framework treats light nuclei as dynamic degrees of freedom alongside nucleons, pions, and $\Delta$-resonances.

1. Kinetic Equations: The time evolution of the Wigner functions (phase-space distributions) $f_\tau(\vec{r}, \vec{p}, t)$ for particle species $\tau$ (including $n, p, d, t, h, \alpha, \pi, \Delta$) is determined by coupled kinetic equations:
   $$(\partial_t + \vec{\nabla}_p \epsilon_\tau \cdot \vec{\nabla}_r - \vec{\nabla}_r \epsilon_\tau \cdot \vec{\nabla}_p)f_\tau = I^{coll}_\tau$$
   where $\epsilon_\tau$ is the single-particle energy derived from a Skyrme pseudopotential, and $I^{coll}_\tau$ represents the collision integral.

2. Collision Integrals and Reactions: The collision integral includes gain and loss terms derived from many-body Green's function expansions. The model explicitly integrates nucleon-induced catalytic reactions for the formation and decay of light nuclei:

   • $NNN \leftrightarrow Nd$
   • $NNNN \leftrightarrow Nt(h)$
   • $NNNNN \leftrightarrow N\alpha$
   • $NNt(h) \leftrightarrow N\alpha$
   • $N\alpha \leftrightarrow dt(h)$
     Transition amplitudes are calculated using the impulse approximation, where the squared modulus of the matrix element is decomposed into the product of the internal momentum-space wave function of the light nucleus (modeled as a Gaussian function) and the elastic nucleon-nucleon scattering amplitude. An energy-dependent factor $F(\sqrt{s})$ is introduced to correct for the neglect of elastic scattering and deficiencies in the approximation; it is adjusted to reproduce experimental nucleon-nucleus decay cross-sections.

3. Implementation of the Mott Effect: The Mott effect is implemented via a phase-space cutoff in the production term of the collision integral. A light nucleus of mass number $A$ can only form if the average nucleon phase-space density $\langle f_N \rangle_A$ in the surrounding medium is below a critical threshold $f^{cut}_A$:
   $$\langle f_N \rangle_A \equiv \int f_N\left(\frac{\vec{P}}{A} + \vec{p}\right)\rho_A(\vec{p})d\vec{p} \le f^{cut}_A$$
   Here, $\rho_A(\vec{p})$ is the internal momentum distribution of the nucleus. A smaller $f^{cut}_A$ corresponds to a stronger Mott effect (easier dissolution).

4. Numerical Solution: The equations are solved using the test-particle approach combined with the lattice Hamiltonian method for drift terms and a stochastic method for the collision integral.

Main Results
The model was applied to central Au+Au collisions at beam energies ranging from $0.25A$ to $1.0A$ GeV.

• Reproduction of Yields: With optimized threshold parameters ($f^{cut}_{A=2} = 0.11$, $f^{cut}_{A=3} = 0.16$, and $f^{cut}_{A=4} = 0.25$), the kinetic approach successfully reproduces the measured light nucleus yields from the FOPI collaboration across the entire energy range.
• Time Evolution: The simulation shows that light nuclei are abundantly produced during the early compression phase of the collision due to high nuclear density. This underscores the necessity of treating them dynamically throughout the entire collision, rather than only at kinetic freeze-out (as in coalescence models).
• Increase in $\alpha$-Particle Yields: The model successfully explains the observed increase in $\alpha$-particle yields at low energies, where $\alpha$ yields exceed those of helium-3. This is attributed to the weaker Mott effect on the $\alpha$-particle. Due to its significantly larger binding energy, the $\alpha$-particle requires a higher phase-space density to dissolve (a smaller Mott momentum) compared to deuterons, tritons, or helium-3. Consequently, a larger threshold parameter ($f^{cut}_{A=4} = 0.25$) is required to fit the data, whereas a smaller value ($f^{cut}_{A=4} = 0.15$) leads to a significant underestimation of the $\alpha$ yield.

Significance and Claims
The work claims to provide the first dynamic description of light nucleus production in intermediate-energy heavy-ion collisions that includes all species up to $A=4$ and explicitly accounts for the Mott effect. The study demonstrates that:

1. Light nuclei must be treated as dynamic degrees of freedom throughout the entire collision evolution to accurately describe yields.
2. The observed increase in $\alpha$-particle yields at low energies is a direct consequence of the Mott effect, specifically the resistance of the $\alpha$-particle to dissolution in the nuclear medium due to its high binding energy.
3. The kinetic approach enables the determination of threshold parameters $f^{cut}_A$, which serve as proxies for the strength of Mott effects in nuclear matter.

The authors note that while the current approach reproduces data for $A \le 4$, describing heavier fragments ($A \ge 5$) formed via spinodal decomposition requires an extension of the formalism to include fluctuations in nucleon phase-space distributions. They propose future applications of this framework to isoscaling, $\alpha$-cluster formation on heavy nuclei, and astrophysical scenarios involving core-collapse supernovae and compact stars.