Measurement of the forward angle $^{12}$C+$^{12}$C fragmentation differential cross sections at 62 MeV/nucleon

Technical Summary: Measurement of Forward-Angle $^{12}\text{C} + ^{12}\text{C}$ Fragmentation Differential Cross Sections at 62 MeV/nucleon

Problem Statement
The Fermi energy regime (approximately 10–100 MeV/nucleon) represents a transitional zone in nuclear reaction dynamics, bridging direct fragmentation reactions and compound-nucleus decay mechanisms. While substantial data exists for incident energies above 100 MeV/nucleon, the experimental database for the $^{12}\text{C} + ^{12}\text{C}$ system within the Fermi-energy regime remains poorly constrained, particularly at forward angles ($\theta_{lab} < 10^\circ$). This region is critical as it accounts for the predominant portion of the total reaction yield. The lack of high-precision differential cross-section data in this range hinders the accurate modeling of nuclear fragmentation, which is essential for fundamental reaction dynamics and practical applications such as carbon-ion radiotherapy (CIRT). In CIRT, precise knowledge of fragmentation cross-sections is required to model secondary particle damage and verify beam range using positron-emitting fragments (e.g., $^{11}\text{C}$, $^{10}\text{C}$) via PET imaging.

Methodology
The study utilized the FAZIA (Forward mass-number and charge-number Identification Array) detector system at the INFN-Laboratori Nazionali del Sud (LNS).

• Experimental Setup: A $^{12}\text{C}$ beam at 62 MeV/nucleon was directed onto a 240 $\mu\text{g/cm}^2$ $^{12}\text{C}$ target. The detection system, positioned 1005 mm downstream, consisted of four block modules covering polar angles from $2^\circ$ to $8^\circ$. Each block contained 16 telescopes, with each telescope comprising a three-layer sequence: a 300 $\mu\text{m}$ Si detector (Si1), a second Si detector (Si2, either 300 or 500 $\mu\text{m}$), and a 10 cm CsI(Tl) scintillator.
• Particle Identification (PID): Identification of fragments with $Z=1$ to $6$ was achieved by combining the $\Delta E-E$ technique with Pulse Shape Analysis (PSA) in the CsI(Tl) crystals. The $\Delta E-E$ method was used for low-energy regions and heavier fragments ($Z \ge 5$), while PSA (utilizing the ratio of fast to total light output) was effective for higher energies and light isotopes.
• Data Processing:
  • Efficiency Correction: Incomplete Energy Deposition (IED) events, where particles scatter out of the crystal or produce secondary particles, were identified using the PSA–$\Delta E$ correlation and corrected via Geant4 simulations.
  • Energy Calibration: A self-consistent calibration was performed by treating detector active thicknesses as free parameters. Calibration points were derived from $\alpha$ sources, elastically scattered beam particles, and punch-through points (PTPs) to minimize extrapolation errors.
  • Normalization: Absolute differential cross-sections were determined by normalizing fragment yields to the elastic scattering of the beam, calculated using the FRESCO code with optical model potentials.

Key Contributions and Results

• High-Precision Data: The authors extracted angular distributions for 16 isotopes ($^1\text{H}$ through $^{13}\text{C}$) in the range $2^\circ \le \theta_{lab} \le 8^\circ$. The results show clear symmetry around $\theta_{lab} = 0^\circ$, validating the consistency of the calibration and corrections.
• Model Comparison:
  • For heavier fragments (e.g., $^{10}\text{B}$, $^{11}\text{C}$), the angular distributions are better described by the Van Bibber formulation, which includes a transverse momentum width component ($W_2$) accounting for Coulomb and nuclear deflection, rather than the standard Goldhaber model.
  • However, the data suggests that the observed distributions are dominated by a wide component from dissipative processes (de-excitation of pre-fragments) rather than orbital deflection, as deflection effects for light targets like carbon are negligible.
  • $\alpha$-Particle Anomaly: $\alpha$ particles exhibit an anomalously narrow angular distribution compared to both models. The authors attribute this to the intrinsic cluster structure of $^{12}\text{C}$, where pre-formed $\alpha$ clusters avoid the random momentum summation of independent nucleons.
• Kinematic Scaling and Energy Dependence:
  • Comparisons with existing data at 62 MeV/nucleon show excellent agreement, with the current study providing a more complete set of species and higher precision.
  • A kinematic scaling analysis was proposed to compare data at 50, 62, and 95 MeV/nucleon.
  • At 95 MeV/nucleon, the scaled theoretical distributions successfully reproduce light isotope yields.
  • At 50 MeV/nucleon, the scaling fails: experimental yields for light fragments are lower than predicted, while heavier fragments are higher. This suggests that at lower energies, the de-excitation (decay) of primary fragments is reduced, leading to a higher survival probability of heavy fragments.

Significance
The paper provides a critical dataset of high-precision forward-angle fragmentation cross-sections for the $^{12}\text{C} + ^{12}\text{C}$ system at 62 MeV/nucleon, filling a gap in the Fermi-energy regime. The results offer new insights into the reaction mechanisms, specifically highlighting the dominance of direct fragmentation and pre-fragment de-excitation over simple orbital deflection for light targets. The observed narrow distribution of $\alpha$ particles supports the role of intrinsic cluster structures in the fragmentation process. These findings contribute to the development of more accurate theoretical models for light nuclear fragmentation and provide necessary data for improving the precision of range verification in carbon-ion radiotherapy. The authors note that further systematic comparisons with transport-model calculations are required to fully elucidate the underlying dynamical mechanisms.