Matter-Induced CPT Violation and Earth-Density Stratification Effects in Long-Baseline Neutrino Oscillation Experiments

技術要約：長基線ニュートリノ振動実験における物質誘起CPT破れと地球密度層状化効果

問題提起
長基線（LBL）ニュートリノ振動実験（T2K、NO$\nu$A、DUNE、Hyper-Kamiokande）は、レプトンCP対称性の破れの位相 $\delta_{CP}$ を精密に測定し、ニュートリノ質量階層を決定することを目的としている。しかし、これらの測定は、地球の物質中を伝播する際に生じる、物理的に結合した2つの系統誤差によって複雑化している：

1. 外因的CPT破れ： 地球の物質中に電子が存在すること（陽電子は存在しない）により、物質ポテンシャル（$V_f$）が生じ、これがニュートリノと反ニュートリノの間で符号を変化させる。これにより、真空中でCP対称性が正確に成立していたとしても、非ゼロのアシンメトリー（$A^{CPT}_{\mu e} \neq 0$）が誘起される。この「偽の」アシンメトリーは、真のCP対称性の破れを孤立させるために差し引かれる必要がある。
2. 地球密度層状化： 基線長 $L \gtrsim 5000$ km の場合、ニュートリノの軌道は下部マントルおよび外核を通過する。標準的な一定の平均密度近似では、これらの変動を捉えることができないため、$\delta_{CP}$ の再構成に系統的なバイアスをもたらす可能性がある。

核心となる問題は、両方の効果が同じ物質ポテンシャル・ハミルトニアン（$H_f \propto \rho(x)$）に依存していることである。したがって、密度のモデリングにおける誤差は、$A^{CPT}_{\mu e}$ の推定値と $\delta_{CP}$ の再構成の両方に同時に伝播する。先行研究では、これらの系統誤差を個別に扱ってきたが、本論文はこれらを結合してモデル化すべきであると主張している。

手法
著者らは、統一された理論的および統計的フレームワークを用いている：

• 理論的フレームワーク： 解析には、厳密な3世代行列指数関数プロパゲータを用いて、進化方程式 $i d|\nu(x)\rangle/dx = H_f(x)|\nu(x)\rangle$ を解いている。ハミルトニアンには、真空項（NuFit 5.3 パラメータを使用）と、参照地球モデル（PREM）から導出されたMSW物質ポテンシャルが含まれる。
• シミュレーション： ニュートリノ振動確率を T2K、NO$\nu$A、DUNE、Hyper-Kamiokande に対して計算している。研究では、外因的CPTアシンメトリー $A^{CPT}_{\mu e} = (P_{\mu e} - P_{\bar{e}\bar{\mu}}) / (P_{\mu e} + P_{\bar{e}\bar{\mu}})$ を評価している。
• 密度不確かさ： 地球物理学的な不確かさに対する $A^{CPT}_{\mu e}$ の感度は、分数密度シフト $\rho \to \rho(1 + \epsilon)$（ここで $\epsilon = \pm 1\%, \pm 2\%, \pm 5\%$）を適用することで定量化されている。
• 統計解析： PREM プロファイルを用いた「真値」に対する定数密度モデルを用いた際の $\delta_{CP}$ 再構成のバイアスを分析するために、ポアソン対数尤度 $\chi^2$ 統計量を構築している。解析には、フラックス、断面積、および検出器効率に関する剰余パラメータのプル項が含まれる。極めて重要な点として、著者らは密度不確かさを $\delta_{CP}$ と $A^{CPT}_{\mu e}$ の両方に影響を与える相関した系統誤差として扱うための、密度プル項（$\xi_\rho$）を含む「拡張 $\chi^2$」へと拡張している。

主な貢献
本研究は、物質誘起の外因的CPT破れと地球密度層状化を、単一の統計的フレームワーク内で相関した系統誤差として扱う最初の体系的な研究の一つである。具体的な貢献は以下の通りである：

1. 結合的な取り扱い： $\delta_{CP}$ 再構成と $A^{CPT}_{\mu e}$ 推定の両方に影響を与える密度不確かさを同時に考慮する、統一された $\chi^2$ フレームワーク。
2. 階層依存性： 正の階層（NH）と負の階層（IH）が $A^{CPT}_{\mu e}$ にどのように影響するかを定量的に特徴付け、DUNE の基線において $\sim 9\%$ の差があることを明らかにしている。
3. 3D 表面解析： $A^{CPT}_{\mu e}(E, \delta_{CP})$ 表面を計算し、固有のCP効果と外因的CPT効果が絡み合う領域（低エネルギー）と、物質が支配的な領域（高エネルギー）を特定した。
4. バイアス定量化： 現在の実験から仮説的な超長基線（VLB）提案に至るまで、幅広い基線長にわたって、定数密度近似によって導入される $\delta_{CP}$ 再構成のバイアスを厳密に評価した。

主な結果

• 外因的CPTアシンメトリー： ピークエネルギーにおいて、$A^{CPT}_{\mu e}$ は正の階層の下で、0.022 (T2K/Hyper-K) から 0.180 (DUNE) の範囲にある。このアシンメトリーは $\theta_{23}, \theta_{13}, \Delta m^2_{31}$ に敏感であり、DUNE において NH と IH の間で $\sim 9\%$ の差がある。
• 密度感度： $\pm 5\%$ の密度不確かさは、DUNE において $|\Delta A^{CPT}| \lesssim 0.008$（約 3–4% の相対値）のシフトを誘起する。$A^{CPT}_{\mu e}$ の比の定義は、密度の誤差に対して部分的な自己キャンセル機能を提供するが、その影響は無視できない。
• 3D 表面領域： DUNE ($L=1285$ km) では、$A^{CPT}_{\mu e}$ 表面は $E \gtrsim 1.5$ GeV で $\delta_{CP}$ に対して平坦であり、単純な差し引きが可能である。しかし、$E \lesssim 0.5$ GeV では、太陽振動項により表面が $\delta_{CP}$ と共に大きく傾斜するため、バイアスを避けるために結合された CP–CPT 尤度解析が必要となる。
• 再構成バイアス：
  • 現在の実験 ($L \le 1285$ km) では、定数密度モデルを使用することによる $\delta_{CP}$ 再構成のバイアスは無視可能（$< 0.3^\circ$）であり、現在の実験精度目標を十分に下回っている。
  • 提案されている VLB ファシリティ ($L > 5000$ km) では、バイアスは壊滅的となる。$L=7000$ km ではバイアスは $17.8^\circ$（DUNE Phase II の目標の 3.6 倍）となり、$L=12000$ km ではバイアスは $172.2^\circ$ に達し、再構成された位相の符号を事実上反転させる。
• 統計的フレームワーク： 密度プル項を含む拡張 $\chi^2$ は、現在の基線においては定数密度近似が妥当であることを裏付けている。しかし、VLB 提案においては、密度不確かさを相関した系統誤差として扱うことが不可欠である。

意義と主張
本論文は、現世代の LBL 実験（T2K, NO$\nu$A, DUNE, Hyper-K）にとって、定数密度近似は CP および CPT 解析の両方において十分であり、物質密度の不確かさは $A^{CPT}_{\mu e}$ に対して管理可能な $\sim 1.3\%$ の系統誤差しか寄与しないと主張している。

しかし、著者らは、将来の超長基線（VLB）提案（例：大陸間基線 $L > 5000$ km）については、地球密度層状化と物質誘起CPT破れは独立した改良点ではなく、「同一の支配的な系統誤差の二つの側面」であると強調している。論文は、そのような施設に関する将来の解析においては、有効な CP および CPT 測定を保証するために、空間的に解像された PREM プロファイルを用い、密度不確かさを結合して扱う、提示された厳密かつ相関した統計的フレームワークを使用することが必須であると結論付けている。