Earth-density effects in long-baseline neutrino experiments in a four-flavor (3+1) sterile-neutrino framework

技術要約：4世代（3+1）ステライルニュートリノ・フレームワークにおける長基線ニュートリノ実験内の地球密度効果

問題提起
長基線（LBL）ニュートリノ実験（例：T2K、NO$\nu$A、DUNE）は、レプトンCP対称性の破れの位相 $\delta_{13}$ およびニュートリノ質量順序を決定するために、ミケエフ・スミルノフ・ヴォルフシュタイン（MSW）物質効果に依存している。標準的な3世代モデルにおける既知の系統的な制限は、地球の物質密度に対する扱いにある。先行研究（文献 [1]）では、地球の密度を一定の経路長平均として近似することは、$L \lesssim 5000$ km の基線に対しては有効な近似であるが、軌道が下部マントルや核を通過する場合（$L \gtrsim 9000$ km）、$\delta_{13}$ において重大な系統的バイアス（数十度）を導入することが示されている。

同時に、短基線実験（LSND、MiniBooNE、反応炉およびガリウム・アノマリー）における異常は、軽ステライルニュートリノ（$\Delta m^2 \sim \mathcal{O}(1)$ eV$^2$）の存在を動機付けている。ステライルニュートリノは中性カレント（NC）弱相互作用には関与しないが、能動的フレーバーとの混合により、追加の物質ポテンシャル項を導入する。本論文が扱う中心的な問題は、(3+1) フレームワークにおける軽ステライルニュートリノの導入が、標準的な3世代パラダイムで特定された地球密度の系統的バイアスを、どのように修正、増強、または抑制するかどうかである。

手法
著者らは、文献 [1] の解析を、フル $4 \times 4$ 行列指数関数法を用いて、4世代（3+1）フレームワークへと拡張している。主な手法上の更新点は以下の通りである：

1. 定式化： 本研究では、6つの回転によってパラメータ化されたユニタリな $4 \times 4$ 混合行列 $U$ を用いる。これには新たな混合角（$\theta_{14}, \theta_{24}, \theta_{34}$）とCP位相（$\delta_{14}$）が含まれる。ハミルトニアンには、荷電カレント（CC）および中性カレント（NC）の両方の物質ポテンシャルが含まれる。決定的なことに、NCポテンシャル $A_{NC}$ は能動的フレーバーには作用するがステライル状態には作用せず、中性子密度に比例するポテンシャルの差を生じさせる。
2. 解析的導出： 著者らは、$\delta_{13}$ 再構成バイアスの主要項となる解析的表現を導出している。$\Delta m^2_{41}$ に伴う高速振動を断熱的に積分除去することで、NCポテンシャルが標準的な3×3 能動的ハミルトニアンにおける実効的な非標準相互作用（NSI）様の項を生成することを証明している。この項は、標準的な $s_{13}$ 駆動のCP項と干渉する。
3. 地球密度モデリング： 先行研究で使用された離散的な4層地球モデルの代わりに、Dziewonski & Anderson (1981) による連続的な多項式 PREM（Preliminary Reference Earth Model）プロファイルを採用している。これにより、経路積分された密度のより精密な計算が可能となり、層の境界における人工的な不連続性が除去される。
4. パラメータ・スキャン： 標準的なベンチマーク・パラメータ（$\Delta m^2_{41} = 0.30$ eV$^2, \theta_{14} = 8.13^\circ, \theta_{24} = 5.40^\circ$）に加え、ステライルCP位相 $\delta_{14} \in \{0, \pi/2, \pi, 3\pi/2\}$ および混合角 $\theta_{34} \in \{0^\circ, 5^\circ\}$ をスキャンする。
5. 数値実装： 数値的な精度を確保し、スキャンのノイズを最小限に抑えるため、高解像度の $400 \times 400$ エネルギー-基線グリッド（$L = 1000–12000$ km）上で計算を行っている。

主要な貢献および解析的洞察
主要な理論的貢献は、$\delta_{13}$ バイアスを修飾するステライル誘起の干渉メカニズムの特定である。

• 実効的NSI項： 高速なステライル振動の積分により、$A_{NC}(x)$ に比例する実効的なハミルトニアン項が得られる。この項は、混合角の積 $s_{14}s_{24}$ と位相 $\delta_{14}$ を通じて $\nu_e$ と $\nu_\mu$ を結合させる。
• 干渉位相： バイアスは相対位相 $\Phi_{NC}(L) = \int_0^L A_{NC}(x) dx - \Delta_{31}(L)$ によって変調される。$A_{NC}(x)$ は地球の非線形な密度プロファイル（特に核・マントル境界および核の急激なジャンプ）を追跡する一方で、$\Delta_{13}(L)$ は $L$ に対して線形に増大するため、干渉項 $\cos[\delta_{13} - \delta_{14} + \Phi_{NC}(L)]$ は基線の増加とともに符号を変える。
• 再分配 vs 再スケーリング： 本論文は、ステライルセクターが単に3世代のバイアスを一様に増幅または抑制するのではないことを示している。むしろ、ステライルセクターはバイアスを非単調に再分配する。基線および $\delta_{14}$ と $\theta_{34}$ の具体的な値に応じて、バイアスは3世代の場合と比較して、抑制（破壊的干渉）または増強（建設的干渉）される。

結果

• 基線依存性： $L \lesssim 5000$ km では、4世代フレームワークにおけるバイアスは3世代のベンチマークと同等かそれ以下に留まる。軌道が下部マントルや核を通過する（$L > 6000$ km）と、バイアスは増大するが、その大きさと符号はステライル・パラメータに基づいて大きく変動する。
• 符号反転： 解析により、軌道が外核を横切る特定の基線（例：9000 km から 11000 km の間）において、ステライル誘起の修飾の符号が反転することが確認された。この反転は、数値スキャンのアーティファクトではなく、密度プロファイルの形状に直接起因するものである。
• パラメータのロバスト性： $\delta_{14}$ と $\theta_{34}$ のスキャンにより、再分配効果が(3+1)フレームワークの一般的な特徴であることが確認された。これらのパラメータを変化させると、符号反転が発生する基線はシフトするものの、非一様な再分配という定性的なパターン自体は消失しない。
• 定量的バイアス： 表IIに具体的なバイアス値を示す。例えば、$L=7000$ km において、4世代のバイアスは3世代の場合（後者は二次的な $\chi^2$ 極小値があることに注意、118.0$^\circ$）と比較して大幅に減少している（8.0$^\circ$）。一方で、$L=11000$ km では、バイアスが増強されている（40.0$^\circ$ に対し 32.0$^\circ$ ではない。※原文の数値関係に基づくと、増強の文脈での比較）。

意義および主張
本論文は、eVスケールのステライルニュートリノの存在が、LBL実験における地球密度の系統的な性質を根本的に変えることを主張している。ステライルセクターは単純なスケーリング因子として機能するのではなく、基線依存の干渉メカニズムを導入し、系統的バイアスを再分配する。

• 構造的普遍性： 著者らは、この再分配が、特定のパラメータ選択に紐付いた数値的な偶然ではなく、NCポテンシャルと地球の密度プロファイルの相互作用によって駆動される、(3+1)フレームワークのロバストで構造的に一般的な特徴であると主張している。
• 解析的根拠： 本研究は、バイアスがなぜ非単調に振る舞うのかについて、NC駆動の結合と標準的なCC駆動の振動との間の位相差に関連付けた解析的な説明（式6および7による）を提供している。
• 実験的含意： 結果は、ステライルニュートリノが存在する場合、一定密度を仮定したり、単純な3世代モデルを使用したりすることは、特に地球の核を通過する基線を持つ実験において、$\delta_{13}$ の再構成に複雑で直感に反する誤差をもたらす可能性があることを示唆している。

結論として、具体的な数値パターンは想定されるステライル・パラメータに依存するものの、再分配の根底にあるメカニズムは(3+1)モデルに固有のものであると述べている。今後の課題として、$\delta_{13}$ とステライル・パラメータの同時フィットを含めること、および大気ニュートリノやIceCube/KM3NeTのデータセットへの拡張が挙げられている。