[论文解读] Measurement of the UHECR energy spectrum using data from the Surface Detector of the Pierre Auger Observatory

M. Roth|arXiv (Cornell University)|Jun 14, 2007

Dark Matter and Cosmic Phenomena被引用 38

一句话总结

本文利用皮埃尔·让天文台表面探测器（SD）对超高能宇宙射线（UHECR）能谱进行了高统计精度测量。通过混合事件（即荧光探测器与表面探测器同时记录同一空气簇射）校准SD信号，建立了一套无需模拟的稳健能量标度，得到的能谱系统不确定度为22%，统计不确定度为6%，在约10^20 eV处显示出能量谱的陡降，与GZK截止理论一致。

ABSTRACT

At the southern site of the Pierre Auger Observatory, which is close to completion, an exposure that significantly exceeds the largest forerunner experiments has already been accumulated. We report a measurement of the cosmic ray energy spectrum based on the high statistics collected by the surface detector. The methods developed to determine the spectrum from reconstructed observables are described. The energy calibration of the observables, which exploits the correlation of surface detector data with fluorescence measurements in hybrid events, is presented in detail. The methods are simple and robust, exploiting the combination of fluorescence detector (FD) and surface detector (SD) and do not rely on detailed numerical simulation or any assumption about the chemical composition. Besides presenting statistical uncertainties, we address the impact of systematic uncertainties.

研究动机与目标

• 利用皮埃尔·让天文台表面探测器对超高能宇宙射线能谱进行高统计精度测量。
• 开发一种稳健的、无需模拟的SD信号绝对能量校准方法，利用混合事件实现。
• 量化并最小化SD能量标度中的系统不确定度，特别是荧光产额、校准和重建相关的不确定度。
• 提供一种对成分或强子模型假设最少的能谱。
• 为未来基于估计能量分辨率的能谱去卷积提供支持。

提出的方法

• 表面探测器通过水箱中的信号幅度测量簇射粒子的横向分布，参数化为S(1000)，即距簇射核心1000米处的信号幅度。
• 采用似然法，利用多个SD站点的数据，重构簇射轴向、S(1000)和前缘曲率。
• 恒定强度法通过拟合观测事件率随cos²θ的变化，提取衰减曲线CIC(θ)，以校正天顶角依赖性。
• S₃₈°定义为S(1000)/CIC(θ)，并归一化至参考天顶角38°用于能量校准。
• 利用同时由FD和SD记录的混合事件，通过线性拟合将S₃₈°校准为FD测热能量E_FD：log E_FD = A + B·log S₃₈°。
• 系统不确定度从FD能量、VAOD测量和重建中传播而来，通过平方和开方组合，估算σ_E_FD和σ_E_SD。

实验结果

研究问题

• RQ1 皮埃尔·让天文台表面探测器在高统计量下测得的UHECR能谱形状如何？
• RQ2 如何在不依赖详细模拟或成分假设的前提下，可靠地校准SD能量标度？
• RQ3 SD能量测量中的主要系统不确定度来源是什么，如何对其进行量化？
• RQ4 观测到的能谱在多大程度上支持GZK截止假说？
• RQ5 基于SD数据推导出的能谱与仅使用混合事件分析的结果相比如何？

主要发现

• 最佳拟合校准关系为log E_FD = 17.08 ± 0.03 + (1.13 ± 0.02)·log S₃₈°，自由度调整卡方值为1.3，表明存在强线性相关性。
• 在10²⁰ eV能量处，SD能量估计的统计不确定度约为5%，而总系统不确定度为22%。
• 最大的系统不确定度来源于绝对荧光产额（14%）、FD校准（9.5%）和重建方法（10%）。
• 能谱在约10²⁰ eV处出现陡降，与GZK截止一致，如相对于E⁻².⁶幂律的分数偏差所示。
• 探测体积为5165 km² sr yr，超过AGASA的三倍以上，支持高统计量分析。
• 该方法稳健且无需模拟，仅依赖混合事件中FD与SD之间的相关性，对成分或强子模型的假设极少。