[论文解读] Towards extracting cosmic magnetic field structures from cosmic-ray arrival directions

M. Wirtz, Teresa Bister+1 more|arXiv (Cornell University)|Jan 8, 2021

Astrophysics and Cosmic Phenomena参考文献 35被引用 4

一句话总结

该论文提出COMPASS方法，一种新颖的似然比方法，通过将可旋转的椭圆密度分布拟合至邻近粒子，以检测超高能宇宙射线（UHECR）到达方向中的一致性磁场结构。该方法在不依赖银河系磁场（GMF）模型的前提下，识别出与天球相切的定向偏转模式，在天体物理模拟中，当源密度 ≥10⁻² Mpc⁻³ 时，拒绝各向同性的显著性水平超过4σ。

ABSTRACT

We present a novel method to search for structures of coherently aligned patterns in ultra-high energy cosmic-ray arrival directions simultaneously across the entire sky. This method can be used to obtain information on the Galactic magnetic field, in particular the integrated component perpendicular to the line of sight, from cosmic-ray data only. Using a likelihood-ratio approach, neighboring cosmic rays are related by rotatable, elliptically shaped density distributions and the significance of their alignment with respect to circular distributions is evaluated. In this way, a vector field tangential to the celestial sphere is fitted which approximates the local deflections in cosmic magnetic fields if significant deflection structures are detected. The sensitivity of the method is evaluated on the basis of astrophysical simulations of the ultra-high energy cosmic-ray sky, where a discriminative power between isotropic and signal-induced scenarios is found.

研究动机与目标

• 在不假设特定银河系磁场（GMF）模型的前提下，检测超高能宇宙射线（UHECR）到达方向中的一致性定向偏转模式。
• 开发一种模型无关的方法，仅从UHECR数据中推断GMF的垂直分量积分。
• 评估该方法在真实天体物理模拟中检测GMF诱导各向异性的灵敏度。
• 优化超参数（如GMF初始化和椭圆几何形状），以实现最大检测能力。

提出的方法

• 该方法在每个UHECR到达方向处使用一个可旋转的椭圆概率密度函数（PDF），以建模一致的偏转模式。
• 通过似然比比较椭圆PDF（信号假设）与圆形PDF（背景假设）的拟合效果，以评估对齐的显著性。
• 每个椭圆的方向由天球上的切向量场 û(ϑ, ϕ) 决定，该向量场通过TensorFlow使用梯度下降法进行优化。
• 该向量场通过球谐函数参数化，以确保空间变化的平滑性，并捕捉大尺度GMF结构。
• 通过超参数λF控制的惩罚项调节向量场的灵活性，平衡对真实信号的敏感性与对各向同性天空中虚假对齐的抑制。
• 该方法在能量 >40 EeV 的模拟UHECR数据上进行训练和验证，使用10⁴个各向同性天空实现以校准假阳性率。

实验结果

研究问题

• RQ1 是否可以在不假设特定银河系磁场模型的前提下，检测到UHECR到达方向中的一致性定向偏转模式？
• RQ2 COMPASS方法在真实天体物理模拟中检测GMF诱导各向异性的灵敏度如何？
• RQ3 超参数（如GMF初始化和椭圆几何形状）如何影响该方法区分信号与各向同性噪声的能力？
• RQ4 该方法能否可靠地从微弱的模拟UHECR源模式中重建局部GMF偏转方向？

主要发现

• 在天体物理模拟中，COMPASS方法对UHECR源密度ρS = 10⁻² Mpc⁻³ 的情况，拒绝各向同性的显著性水平超过4σ。
• 最低的各向同性偶然概率（pval）为2.8 × 10⁻³，对应椭圆几何形状(δmax, δmin) = (20°, 15°)，显著优于更小或更细长的椭圆。
• 在λF ≈ 0.5 处发现各向同性偶然概率的最小值，表明该值在模型灵活性与对各向同性涨落引起的虚假信号抑制之间达到最佳平衡。
• 该方法在模拟中成功重建了方向偏转模式，向量场方向与模拟源方向对齐。
• 即使假设的GMF模型（JF12）受到±45°的偶极不确定性扰动，该方法仍保持稳健，pval在最优λF附近保持稳定。
• 在似然比中使用圆形参考模型，能有效将一致的椭圆模式与随机过密区分开。