[论文解读] High energy extension of the FLUKA atmospheric neutrino flux

G. Battistoni, Anna Ferrari+2 more|ArXiv.org|May 19, 2003

Astrophysics and Cosmic Phenomena被引用 11

一句话总结

该论文将FLUKA大气中微子通量模拟扩展至10^4 GeV，使高能中微子实验（如向上穿过μ子的探测）能够实现精确建模。通过在能量区间内并行模拟并采用加权求和方法，实现了与能量无关的低统计不确定性，揭示出在1 TeV以上能量区域，初级通量不确定性主导系统误差，超过20%。

ABSTRACT

The atmospheric neutrino flux calculated with FLUKA was originally limited to 100-200 GeV for statistical reasons. In order to make it available for the analysis of high energy events, like upward through-going muons detected by neutrino telescopes, we have extended the calculation so to provide a reliable neutrino yield per primary nucleon up to about 10**6 GeV/nucleon, as far as the interaction model is concerned. We point out that the primary flux model above 100 GeV/nucleon still contributes with an important systematic error to the neutrino flux.

研究动机与目标

• 将FLUKA大气中微子通量计算扩展至100 GeV以上，以支持高能中微子实验（如向上穿过μ子的探测）。
• 克服以往FLUKA模拟中因计算成本过高而导致结果受限于约100 GeV的统计局限性。
• 量化初级宇宙射线通量不确定性对高能中微子通量预测的影响，特别是在1 TeV以上能量区域。
• 为不同初级能谱（包括2001年Bartol拟合与旧版Bartol拟合）提供可靠且高精度的中微子通量表格。

提出的方法

• 在从~30 GeV到10^6 GeV/nucleon的窄能量区间内进行并行模拟运行，以确保统计不确定性与能量无关。
• 采用加权求和方案，重建任意所需初级全核素能谱的中微子通量，避免仅依赖单一能谱采样。
• 保持与先前FLUKA研究相同的物理模型和相互作用算法，利用其高精度的强子和电磁级联簇射模型。
• 在单个半球内记录10个天顶角区间的中微子角分布，以支持方向性通量分析。
• 将模拟能量上限扩展至10^6 GeV/nucleon，以确保在10^4 GeV能量区域的中微子产额可靠。
• 对比使用两种不同初级通量模型的结果：2001年Bartol拟合与1996年Bartol拟合。

实验结果

研究问题

• RQ1 如何将FLUKA大气中微子通量扩展至适用于向上穿过μ子实验的能量范围，该实验要求通量达到至少10^4 GeV？
• RQ2 不同初级宇宙射线能谱对高能大气中微子通量预测有何影响？
• RQ3 在1 TeV以上能量区域，初级宇宙射线通量不确定性在高能中微子通量系统误差中占多大比重？
• RQ4 能否通过并行模拟与加权求和方法，将高能中微子通量计算的统计不确定性最小化？

主要发现

• FLUKA大气中微子通量已成功扩展至10^4 GeV，使其可用于高能向上穿过μ子事例的分析。
• 由于采用并行模拟与加权求和方法，中微子通量的统计不确定性现几乎与能量无关。
• 初级全核素能谱的选择对高能中微子通量有显著影响，在10^4 GeV能量区间内差异超过两倍。
• 在1 TeV以上能量区域，初级通量不确定性估计对中微子通量预测的总系统误差贡献超过20%。
• 2001年Bartol拟合的初级能谱比旧版拟合更偏软，导致高能中微子通量显著降低，尤其对νμ和ν̄μ更为明显。
• 在假设初级通量为幂律谱（E^-γ）且γ值在2.6至3.0之间的外推结果表明，初级通量的高能行为对最终中微子通量有强烈影响。