Interaction of accelerator neutrinos with energies up to 55 MeV with ${}^{127}$I nuclei

This study investigates the interaction of accelerator neutrinos up to 55 MeV with ${}^{127}$I nuclei at the SNS, calculating cross sections that reveal the dominant contribution of the Gamow-Teller resonance GTR-1 (60–80%) alongside significant effects from higher resonances GTR-2 and AR-2, with theoretical results showing agreement with experimental data below the neutron separation threshold.

The Gist

Imagine the universe is filled with invisible messengers called neutrinos. They are like cosmic ghosts: they zip through everything—stars, planets, and even your body—without ever saying "hello." To catch them, scientists build massive detectors, like giant fishing nets made of specific atoms.

This paper is about a new, very sensitive fishing net made of Iodine-127 (a specific type of iodine atom). The scientists wanted to see how well this net catches "fast" neutrinos coming from a powerful machine called the Spallation Neutron Source (SNS) in the US, which shoots neutrinos at speeds up to 55 MeV (much faster than the slow neutrinos coming from our Sun).

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

1. The "Bouncy Castle" Analogy (The Nucleus)

Think of the Iodine nucleus as a bouncy castle.

• Solar Neutrinos (Slow): These are like toddlers gently bouncing on the castle. They don't have enough energy to break anything. They just make the castle wiggle a little. Scientists have known exactly how the castle wiggles for a long time.
• Accelerator Neutrinos (Fast): These are like professional parkour athletes jumping onto the castle. They hit it so hard that the castle doesn't just wiggle; it starts to shake violently, and pieces might even fly off (neutrons popping out).

The problem? No one had ever mapped out exactly how the castle behaves when hit by these "parkour athletes" (energies above 20 MeV). The scientists had to predict the physics of these violent collisions.

2. The "Musical Resonance" (The Strength Function)

When you hit a bell, it rings at a specific pitch. If you hit it harder, it might ring at a different, higher pitch or even shatter. In physics, this "ringing" is called a resonance.

The authors calculated a "sound map" (called the Strength Function) for the Iodine nucleus. They discovered that the nucleus has different "notes" it can play:

• The Main Chord (GTR-1): The loudest, most common note. This accounts for about 60–80% of the reaction. It's the main way the nucleus absorbs the hit.
• The High Notes (GTR-2 & AR-2): These are higher-pitched, harder-to-reach notes that only get played when the neutrino hits really hard (high energy). The scientists found these new "high notes" for the first time. They contribute about 12% and 10% respectively.

The Analogy: Imagine trying to predict how a car crash will damage a vehicle. You know how it handles a fender bender (low energy), but you needed a new model to predict what happens in a high-speed crash (55 MeV). They found that the car has specific "weak spots" (resonances) that only break under extreme speed.

3. The "Missing Puzzle Piece" (The Discrepancy)

Here is where the story gets tricky. The scientists ran their calculations (their "crash simulation") and compared them to real-world data from an experiment called COHERENT.

• The Good News: When the neutrinos hit the Iodine gently (below the energy needed to knock a neutron out), the math matched the experiment perfectly. The "fender bender" predictions were spot on.
• The Bad News: When the neutrinos hit hard enough to knock a neutron out (the "high-speed crash"), the math and the experiment disagreed wildly.
  • The experiment saw fewer "broken pieces" (neutrons) than the math predicted.
  • The math predicted a lot more "broken pieces" than the experiment found.

Why the mismatch?
The authors suggest a few reasons:

1. The "Volume Knob" (Quenching): Maybe the nucleus is "quieter" than we thought. In physics, there's a rule that says how much energy a nucleus can absorb. It seems the nucleus might be "turning down the volume" (a phenomenon called quenching) when hit hard, which our current models don't fully account for.
2. The "Blind Spot": We don't have a complete map of the "high notes" (resonances) above 20 MeV. We are guessing how the castle behaves at the highest energies, and our guess might be wrong.

4. The Conclusion: "We Need a Better Map"

The paper ends with a call to action. The scientists say:

"Our math works great for slow neutrinos, but for the fast ones, we are missing data. We need to build a better 'sound map' of the Iodine nucleus at high energies."

They suggest doing a new, more precise experiment (using a different type of particle beam, like a 3He beam) to measure exactly how the Iodine nucleus reacts when hit hard.

Summary in One Sentence

Scientists modeled how a fast-moving neutrino hits an Iodine atom, finding that while their math works perfectly for gentle hits, it fails to predict the results of violent hits, suggesting we need new experiments to understand the "high-energy music" of the atomic nucleus.

Technical Summary

1. Problem Statement

The paper addresses a critical gap in understanding neutrino-nucleus interactions at high energies (up to 55 MeV), specifically relevant for the Spallation Neutron Source (SNS) accelerator neutrino experiments (e.g., the COHERENT collaboration).

• The Gap: While neutrino interactions with $^{127}\text{I}$ are well-understood for solar neutrinos (energies < 20–30 MeV), experimental data and theoretical calculations for the charge-exchange strength function $S(E)$ are lacking for excitation energies above 20 MeV.
• The Discrepancy: Recent experimental measurements by the COHERENT collaboration showed a significant discrepancy between measured and theoretically calculated cross-sections for neutrino capture events involving neutron emission ($\langle\sigma(\ge 1n)\rangle$) at energies above the neutron separation threshold ($E_x > S_{1n}$).
• The Need: To resolve this, a comprehensive microscopic calculation of the strength function $S(E)$ up to 60 MeV is required to account for high-lying resonances that dominate the interaction at these energies.

2. Methodology

The authors employed a microscopic theory of finite Fermi systems (TFFS) to model the interaction.

• Theoretical Framework:
  • Calculations were based on the Saxon–Woods model for one-particle structure and the Landau–Migdal local nucleon-nucleon interaction ($F_\omega$).
  • The system of secular equations for the effective field was solved for allowed transitions.
  • Resonance Identification: The study identified not only the well-known low-lying resonances (Gamow–Teller GTR-1, Analog AR-1, and Pygmy PR) but also predicted and calculated high-lying secondary resonances:
    • GTR-2: Associated with single-particle transitions where the main quantum number changes by $\Delta n = 1$ (e.g., $2p_{1/2} \to 3p_{3/2}$).
    • AR-2: Associated with transitions like $(n-p)$ with $\Delta n = 1$.
• Strength Function $S(E)$:
  • Constructed using Breit–Wigner formulas to envelope the calculated matrix elements.
  • Quenching Effect: The authors applied a quenching factor ($q \approx 0.81$, corresponding to an effective charge $e_q = 0.9$) to normalize the theoretical sum rule to experimental data from the $^{127}\text{I}(p,n)^{127}\text{Xe}$ reaction up to 20 MeV.
• Cross-Section Calculation:
  • The neutrino capture cross-section $\sigma(E_\nu)$ was calculated using the standard weak interaction formula, integrating the strength function $S(E)$ over the neutrino energy spectrum.
  • Fermi Function: Due to the lack of tabulated Fermi function values for high electron energies ($E_e > 15$ MeV), the authors used a linear extrapolation of existing tables.
  • Spectrum Averaging: Cross-sections were averaged over the specific energy spectrum of electron neutrinos from $\pi^+$ decay at the SNS.

3. Key Contributions

• Extension of $S(E)$ to High Energies: The paper provides the first theoretical calculation of the charge-exchange strength function for $^{127}\text{I}$ extending up to 60 MeV, covering the full energy range of SNS accelerator neutrinos.
• Identification of High-Lying Resonances: The study explicitly calculates the contributions of GTR-2 and AR-2 resonances, which arise from transitions with $\Delta n = 1$. These were previously unaccounted for in standard models for this energy range.
• Quantitative Breakdown of Contributions: The authors provide a detailed decomposition of how different resonances contribute to the cross-section in three distinct energy regimes:
  1. $E_x < S_{1n}$ (No neutron emission)
  2. $S_{1n} < E_x < S_{2n}$ (One neutron emission)
  3. $E_x > S_{2n}$ (Two neutron emission)
• Resolution of Quenching: The paper rigorously addresses the normalization of the strength function, determining the optimal quenching parameter ($q=0.81$) that aligns theory with low-energy experimental data.

4. Results

• Resonance Contributions:
  • GTR-1 remains the dominant contributor across all energy ranges, accounting for 60–80% of the total cross-section.
  • GTR-2 contributes approximately 12% to the total cross-section, primarily at higher energies.
  • AR-2 contributes $\le 10$%, with its effect becoming noticeable only above the two-neutron separation threshold ($E_x > S_{2n}$).
  • AR-1 contributes significantly (34%) in the single-neutron emission region ($S_{1n} < E_x < S_{2n}$).
• Cross-Section Comparison:
  • Low Energy ($E_x < S_{1n}$): The calculated spectrum-averaged cross-section without neutron emission ($\langle\sigma(0n)\rangle = 1.94 \times 10^{-40} \text{ cm}^2$) shows excellent agreement with both previous experimental data (Los Alamos) and COHERENT measurements.
  • High Energy ($E_x > S_{1n}$): The calculated cross-sections for events with neutron emission ($\langle\sigma(\ge 1n)\rangle$) align well with MARLEY Monte Carlo simulations but show a strong discrepancy with the experimental measurements reported by the COHERENT collaboration.
• Specific Values (in units of $10^{-40} \text{ cm}^2$):
  • $\langle\sigma(0n)\rangle_{\text{TFFS}}$: 1.94
  • $\langle\sigma(1n)\rangle_{\text{TFFS}}$: 16.07
  • $\langle\sigma(2n)\rangle_{\text{TFFS}}$: 2.88

5. Significance and Conclusion

• Validation of Low-Energy Models: The agreement at energies below the neutron separation threshold confirms that the microscopic TFFS approach, when properly normalized, accurately describes neutrino capture on $^{127}\text{I}$ for bound-state transitions.
• Highlighting Theoretical/Experimental Gaps: The persistent discrepancy at higher energies ($E_x > S_{1n}$) suggests that current theoretical models (even with the inclusion of GTR-2/AR-2) or experimental analyses may be missing critical physics. Potential causes include:
  • Uncertainties in the extrapolation of the Fermi function at high energies.
  • Missing nuclear structure effects (e.g., multi-particle configurations) not captured by the current TFFS model.
  • Experimental systematic errors in the COHERENT measurement of neutron-emission events.
• Future Directions: The authors conclude that a new high-resolution measurement of the strength function $S(E)$ is urgently needed, specifically using reactions like $^{127}\text{I}(^3\text{He}, t)^{127}\text{Xe}$ at energies above 20 MeV, to resolve the discrepancy and improve the accuracy of neutrino detectors for accelerator-based experiments.

In summary, this paper provides a robust theoretical baseline for high-energy neutrino interactions with Iodine, confirming the dominance of GTR-1 while highlighting a critical unresolved conflict between theory and experiment in the neutron-emission regime.