Empirical formula for total inelastic cross-section of proton-nucleus scattering

The authors propose a universal empirical formula for proton-nucleus total inelastic cross-sections, valid from 15 MeV to 1 TeV, which is parameterized using experimental data for Carbon and Aluminium and validated against a wide range of target nuclei and existing models like GEANT4.

The Gist

Imagine you are throwing a tennis ball (a proton) at a wall made of different materials (atomic nuclei). Sometimes, the ball bounces off cleanly (elastic scattering). But often, the ball hits the wall, shatters a few bricks, creates a cloud of dust, or even knocks a piece of the wall off entirely. This messy, destructive collision is called inelastic scattering.

Physicists need to know exactly how likely this "messy" collision is to happen. This likelihood is called the cross-section. It's like a target size: the bigger the cross-section, the more likely the ball is to hit something and cause chaos.

This paper is about creating a new, better "rulebook" (an empirical formula) to predict how big that target size is for a tennis ball hitting 33 different types of walls, from tiny hydrogen atoms to massive uranium atoms, across a huge range of speeds (energies).

Here is the breakdown of their work using simple analogies:

1. The Problem: The Old Maps Were Flawed

For decades, scientists used old maps (previous formulas) to predict these collisions.

• The "Power Law" Map: This was like a simple ruler. It said, "If the wall is twice as big, the target is twice as big." It worked okay for some walls but failed miserably for others.
• The "Letaw" Map: This was a more detailed map that worked well for fast balls but got confused when the ball was moving slowly. It didn't account for the fact that slow balls behave differently than fast ones.
• The "Complex" Maps: Some newer maps were so complicated (filled with thousands of tiny rules) that they were hard to use and still made mistakes with light or heavy walls.

The authors realized they needed a new map that was simple to use but accurate everywhere, from very slow balls to ultra-fast ones.

2. The Solution: A Three-Part Recipe

The authors built a new formula that acts like a three-layer cake to predict the collision size:

• Layer 1: The High-Speed Highway (High Energy)
  When the proton is moving super fast (like a bullet), it doesn't care about the tiny details of the nucleus. It just sees the overall size. The authors found that for these fast speeds, the collision size depends mostly on how heavy the nucleus is. They used a mathematical curve (a mix of a power law and a wavy line) to fit the data for heavy elements like Lead and Uranium. Think of this as measuring the "footprint" of the nucleus.

• Layer 2: The Speed Bump (Low Energy)
  When the proton is moving slower, things get weird. The collision size doesn't just stay the same; it dips down and then spikes up.

  • The Dip: Around a certain speed, the proton slips through the nucleus more easily.
  • The Spike: At very low speeds, the proton gets "stuck" or interacts more violently, causing a huge spike in the collision size.
    The authors added a special "wiggly" mathematical term to capture this dip and the spike, ensuring the formula doesn't just go flat when it should be jumping up.

• Layer 3: The Magnetic Repulsion (Coulomb Effect)
  Protons are positively charged. Nuclei are also positively charged. Like two magnets with the same pole facing each other, they repel.

  • If the proton is slow, this repulsion is strong, pushing the proton away before it can hit.
  • The authors added a "Coulomb term" to their formula. It acts like a force field that gets stronger for heavier nuclei (which have more charge) and slows down the proton. This explains why heavy elements behave differently than light ones at low speeds.

3. The Test Drive

To prove their new rulebook works, they took it for a spin against real-world data:

• The Track: They tested it on 33 different "walls" (nuclei), ranging from Deuterium (very light) to Uranium (very heavy).
• The Competitors: They compared their results against the old maps (Letaw, Shen, Tripathi, Nakano) and a high-tech computer simulation called GEANT4 (which is like a video game physics engine for particles).
• The Result: Their new formula was the most consistent. It didn't overestimate the size for heavy walls or underestimate it for light ones. It matched the real-world data better than the old formulas, especially in the tricky low-speed zone.

4. Why Does This Matter? (Real World Applications)

Why do we care about a tennis ball hitting a wall?

• Space Travel (Cosmic Rays): High-energy particles from space (cosmic rays) constantly bombard our atmosphere and spacecraft. To protect astronauts, we need to know how these particles interact with the hull of a spaceship. This new formula helps calculate how long a cosmic ray survives before hitting something, helping engineers design better shields.
• Medical Therapy: In proton therapy for cancer, doctors shoot protons at tumors. They need to know exactly how the protons will scatter inside the body to destroy the tumor without hurting healthy tissue.
• Particle Physics Experiments: In labs like Fermilab, scientists smash particles together to find new physics. They need to know the "background noise" (like antiprotons being created accidentally) to spot the rare signals they are looking for. This formula helps them predict that noise more accurately.

The Bottom Line

The authors created a universal, simple, and accurate calculator for proton collisions. It's like upgrading from a rough sketch of a map to a GPS that works perfectly whether you are driving a slow car through a city or a rocket through space. It bridges the gap between simple rules and complex physics, making it easier for scientists to design safer space missions, better cancer treatments, and more accurate experiments.

Technical Summary

1. Problem Statement

Proton-nucleus inelastic cross-sections are fundamental observables for understanding hadronic reactions in nuclear media, with critical applications in medical radiation therapy, accelerator shielding, space radiation protection, and cosmic ray physics. While several empirical models exist (e.g., Letaw et al., Shen, Tripathi, Nakano), they suffer from significant limitations:

• Lack of Universality: Existing models often fail to provide uniform precision across the entire energy spectrum (from low MeV to TeV) and for all target masses (from light Deuterium to heavy Uranium).
• Low-Energy Inaccuracy: Many models struggle to accurately predict cross-sections below the Coulomb barrier or for light nuclei, often missing characteristic peaks or oscillatory behaviors.
• Complexity vs. Simplicity: Some high-accuracy models (like Nakano et al.) are computationally complex with numerous parameters, while simpler models (like Power Law) lack precision for heavy nuclei.
• Data Gaps: Experimental data is sparse for certain nuclei and energy ranges, necessitating a reliable empirical formula to predict unmeasured cross-sections.

The authors aim to develop a generic, computationally simple empirical formula that offers high precision across a wide energy range (15 MeV to 1 TeV) for 33 different target nuclei.

2. Methodology

The proposed model is a factorized empirical formula that separates the cross-section into high-energy and low-energy regimes, incorporating mass number ($A$), atomic number ($Z$), and incident proton kinetic energy ($E$).

A. High-Energy Regime ($E > 5$ GeV)

In this region, the cross-section is dominated by nuclear geometry and becomes nearly energy-independent, approaching an asymptotic value based on the target mass.

• Formula:
  $$\sigma_h(A) = a A^b \left[ 1 - c \cos(d \sqrt{A}) + e A^f \right]$$
• Parameters: Determined via least chi-square fitting to experimental data (Bobchenko et al.) for 5–9 GeV.
  • $a, b$: Geometric approximation terms.
  • $c, d$: Oscillatory terms accounting for nuclear structure effects.
  • $e, f$: Correction terms specifically for heavy nuclei trends.
• Performance: Achieves a reduced chi-square per degree of freedom of 0.21, significantly outperforming the Power Law and Letaw models for heavy elements ($A > 130$).

B. Low-Energy Regime ($E < 2$ GeV)

At lower energies, the cross-section exhibits strong energy dependence, including a sharp peak below 200 MeV (often >60% higher than saturation) and a dip around 1 GeV.

• Factorization: The total inelastic cross-section is expressed as:
  $$\sigma_{inel}(Z, A, E) = \sigma_h(A) \cdot X(E) \cdot Y_{coul}(Z, E)$$
• Energy Correction Term ($X(E)$): Captures the oscillatory behavior and the dip.
  $$X(E) = 1 - \alpha e^{-E/\beta} \sin(\gamma E^\lambda) + \frac{1}{2} e^{-E/\epsilon}$$
  • Parameters ($\alpha, \beta, \gamma, \lambda, \epsilon$) were fitted using machine learning methods to extensive experimental data for Carbon ($^{12}$C) and Aluminium ($^{27}$Al).
• Coulomb Correction Term ($Y_{coul}(Z, E)$): Accounts for the repulsive Coulomb barrier at low energies, which varies with the target's atomic number ($Z$).
  $$Y_{coul}(Z, E) = 1 + k(Z) \cdot e^{-E/24}$$
  • The parameter $k(Z)$ is a step-function dependent on $Z$, derived to model the suppression or enhancement of the low-energy peak. For light nuclei ($Z \le 18$), $k(Z)$ is often positive (enhancing the peak), while for heavier nuclei, it is negative (suppressing the peak).

3. Key Contributions

1. Unified Formula: A single analytical expression covering 15 MeV to 1 TeV, eliminating the need for separate, disjointed models for different energy regimes.
2. Broad Applicability: Validated against experimental data for 33 target nuclei, ranging from light elements (Deuterium, Helium) to heavy elements (Lead, Uranium).
3. Simplicity: The model uses a minimal number of parameters compared to complex microscopic models (like Nakano et al.) while maintaining high accuracy.
4. Machine Learning Integration: Utilized machine learning techniques to optimize the low-energy fitting parameters based on Carbon and Aluminium data, ensuring robustness before applying the model to other nuclei.
5. Coulomb Parameterization: Introduced a systematic, $Z$-dependent function $k(Z)$ to accurately model the Coulomb barrier effects across the periodic table.

4. Results and Validation

The model was benchmarked against:

• Experimental Data: 33 nuclei (including $^2$H, $^4$He, $^{12}$C, $^{27}$Al, $^{56}$Fe, $^{208}$Pb, $^{238}$U).
• Existing Models: Letaw et al., Shen, Tripathi, and Nakano.
• Simulations: GEANT4 (using the FTFP_BERT physics list) and MCNP.

Key Findings:

• Accuracy: The proposed model consistently shows the smallest relative error (RE) across the majority of the 33 elements compared to other models. For example, for $^2$H, the RE is 0.59 (vs. 1.75 for Letaw); for $^{238}$U, it is 0.71 (vs. 0.90 for Letaw).
• Low-Energy Peaks: Unlike Letaw and Shen models, which often underestimate or miss the low-energy peak for heavy nuclei, the proposed model accurately reproduces the peak position and magnitude.
• High-Energy Saturation: The model correctly predicts the asymptotic saturation values for heavy nuclei, where previous models often showed large deviations.
• Specific Cases:
  • Light Nuclei ($^2$H, $^4$He): The model reproduces experimental data well in the low-energy region where other models fail (e.g., Nakano predicts vanishing cross-sections below 90 MeV for Deuterium).
  • Heavy Nuclei ($^{238}$U): The model avoids the "huge overestimated peak" seen in the Letaw model at low energies.

5. Significance and Applications

The paper demonstrates the utility of the new formula in two critical astrophysical and experimental physics scenarios:

1. Cosmic Ray Lifetime Estimation:

   • The formula was used to calculate the interaction rate and mean free path of cosmic rays traversing the Interstellar Medium (ISM).
   • Result: At 20 MeV, the model predicts a mean free path of $1.12 \times 10^7$ light years, showing a ~20% improvement in accuracy over Shen and Letaw models due to better low-energy peak modeling for light ISM targets (H and He).

2. Antiproton Background in Mu2e Experiment:

   • The formula was applied to estimate antiproton production rates in the Mu2e experiment (Fermilab), a critical background for muon-to-electron conversion searches.
   • Result: The model predicts an antiproton production rate of $6.620 \times 10^{-5}$ per proton on target (POT). This is consistent with the Letaw model (within 1%) but deviates by ~11% from the MCNP simulation. The authors argue their result is more reliable due to the refined parametrization of the inelastic cross-section.

Conclusion

The authors successfully propose a simple yet highly precise empirical formula for proton-nucleus inelastic cross-sections. By factorizing the problem into mass-dependent high-energy saturation and energy/Coulomb-dependent low-energy corrections, the model bridges the gap between computational simplicity and physical accuracy. It outperforms existing standard models across a vast range of energies and target masses, making it a valuable tool for nuclear physics phenomenology, radiation shielding design, and cosmic ray propagation studies.