[論文レビュー] Statistical Methods for the Neutrino Mass Hierarchy

Emilio Ciuffoli|arXiv (Cornell University)|Apr 26, 2017

Neutrino Physics Research被引用数 4

ひとこと要約

本稿は、$Δ\chi^{2}$ を統計的検定量として用い、ニュートリノ質量階層を特定する統計的課題に取り組み、非ネストされた仮説のためウィルクスの定理が適用されないことを示している。代わりに、特定の条件下では、$Δ\chi^{2}$ は $σ = 2\sqrt{|¯Δ\chi^{2}|}$ を持つガウス分布に従うことが示され、頻度主義的およびベイズ的枠組みの両方で感受性の正確な評価が可能になる。また、周辺化と最小化の違いに関する重要な洞察と、非ネスト問題を解消するための新しいプルパラメータの提案がなされている。

ABSTRACT

In the next decade several experiments will attempt to determine the neutrino mass hierarchy, i.e. the sign of $Δm_{31}^2$. In the last years it was noticed that the two hierarchies are disjoint hypotheses and, for this reason, Wilks' theorem cannot be applied: this means that $Δχ^2=χ^2_{IH}-χ^2_{NH}$ does not follow a one-degree-of-freedom chi-square distribution. It was proven that, under certain assumptions, it follows instead a Gaussian distribution with $σ=2\sqrtμ$. I will present several possible definition of sensitivity and review the approaches proposed in the literature, both within the Bayesian and the frequentist framework, examining advantages and disadvantages and discussing how they should be modified if the conditions for Gaussianity are not fulfilled. I will also discuss the possibility of introducing a new pull parameter in order to avoid the issue related to the non-nested hypotheses and the differences between marginalization and minimization, showing under which conditions the two procedures yield the same $Δχ^2$.

研究の動機と目的

• 非ネストされた仮説（通常階層と逆転階層）のため、ウィルクスの定理が適用されないニュートリノ質量階層の決定における統計的課題に取り組む。
• 質量階層実験における感受性の評価・比較のための頻度主義的およびベイズ的アプローチを評価する。
• プルパラメータの周辺化と最小化が等価な$Δ\chi^{2}$値をもたらす条件を調査する。
• 非物理的プルパラメータ$η$を提案・分析し、非ネスト仮説問題をパラメータフィッティング問題に変換する。
• $Δ\chi^{2}$ がガウス分布に従う条件を明確にし、信頼性のある$n\sigma$レベルの解釈を可能にする。

提案手法

• $Δ\chi^{2} = \chi^{2}_{IH}(\hat{\theta}) - \chi^{2}_{NH}(\hat{\hat{\theta}})$ を検定量として用い、$\hat{\theta}$ と $\hat{\hat{\theta}}$ はそれぞれの階層下でのプルパラメータの最尤推定値である。
• アサイモフデータセットを用いて期待される$Δ\chi^{2}$を計算し、中央実験を$Δ\chi^{2} = \overline{Δ\chi^{2}}$ として定義する。
• 線形関数として$y_i(\theta_j)$ が$\theta_j$ に依存し、両階層の最小値近傍で超平面が平行である場合に限り、$Δ\chi^{2}$ のガウス近似を用い、$\mu = \overline{Δ\chi^{2}}$ および $\sigma = 2\sqrt{|\overline{Δ\chi^{2}}|}$ とする。
• p値と後確率を用いて頻度主義的およびベイズ的感受性を比較し、ベイズ的手法が$P(NH|D) + P(IH|D) = 1$ を満たす一貫した確率を提供することに注意する。
• 通常階層と逆転階層の間を補間する非物理的プルパラメータ$η$を導入し、モデル選択問題をパラメータフィッティング問題に変換する。
• ラプラス近似を用いて、周辺化と最小化が等価な$Δ\chi^{2}$ をもたらす条件を評価し、これは鋭い尤度関数と等しいヘッセ行列の行列式を要件とする。

実験結果

リサーチクエスチョン

• RQ1 ニュートリノ質量階層の非ネスト仮説検定において、$Δ\chi^{2}$ が$\chi^{2}$分布ではなくガウス分布に従う条件は何か？
• RQ2 頻度主義的およびベイズ的感受性評価アプローチにはどのような違いがあり、どのような点で補完的か？
• RQ3 ニュートリノ質量階層解析において、プルパラメータの周辺化と最小化が等価な$Δ\chi^{2}$値をもたらす条件は何か？
• RQ4 非物理的プルパラメータ$η$を用いて質量階層決定における非ネスト仮説問題を解消できるか？
• RQ5 反応炉ニュートリノ実験と加速器ニュートリノ実験の間で、$Δ\chi^{2}$ のガウス近似の性能はどのように異なるか？

主な発見

• 非ネスト仮説のため、検定統計量$Δ\chi^{2}$ は$\chi^{2}$分布に従わず、ウィルクスの定理の直接適用が不可能である。
• 適切な条件下では、$Δ\chi^{2}$ は$σ = 2\sqrt{|\overline{Δ\chi^{2}}|}$ を持つガウス分布でよく近似され、$n\sigma$レベルの解釈が可能になる。
• ベイズ的手法は一貫した確率$P(NH|D)$ を提供するが、頻度主義的手法では各仮説に対して別個のp値が必要であり、両者は補完的な情報を提供する。
• 周辺化と最小化が等価な$Δ\chi^{2}$ をもたらすのは、尤度関数が鋭くピークしている場合かつヘッセ行列の行列式が等しい場合に限られ、これは反応炉ニュートリノ実験では満たされるが、加速器ニュートリノモデルでは満たされない。
• 提案された$η$パラメータアプローチは、非ネスト問題をパラメータフィッティング問題に成功裏に変換するが、すべての状況に適用可能ではない（例：物質効果のため、加速器ニュートリノでは有効でない）。
• 中央感受性（頻度主義的）と交差感受性（ベイズ的）は異なる挙動を示し、ベイズ的手法は証拠の要約としてよりコン act で解釈可能な情報を提供する。