Algebraic Structure of Three-Flavor Neutrino Oscillations in Constant-Density Matter: Cayley--Hamilton Evolution, DMP Resummation, and Closed-Form Uncertainty Propagation

技術要約：定数密度物質中における3世代ニュートリノ振動の代数的構造

問題提起
ニュートリノ振動物理学は精密測定の時代へと移行しているが、物質中での振動確率の計算は、DUNE、Hyper-Kamiokande (HK)、JUNOといった次世代の実験にとって計算上のボトルネックとなっている。これらの実験では、質量順序、$\theta_{23}$ オクタント、およびレプトンCP対称性の破れ（$\delta_{CP}$）を解明するために、$O(10^6)$–$O(10^7)$ 回の評価をフィットごとに必要とする。既存のツール、例えば Denton–Parke NuFast-LBL アルゴリズムは高速かつ高精度であるが、反復解法に依存しており、解析的な透明性に欠けている。具体的には、既存の文献では、進化演算子の代数的構造、MSW共鳴（$\hat{A} \to 1$）付近で Denton–Minakata–Parke (DMP) 近似が摂動的発散を回避する具体的なメカニズム、および NuFIT 6.0 パラメータ基底における厳密な不確かさ伝播に必要な閉形式のヤコビアンが、しばに明示されていない。本論文は、これらの要素を統一された代数的枠組みへと統合することを目的とする。

手法
著者らは、定数密度物質中における3世代ニュートリノ振動に関する包括的な解析的手法を、以下の柱に基づいて展開する：

1. Cayley–Hamilton による代数的進化： 本論文は、フレーバー基底のハミルトニアン（$H_f$）に対して Cayley–Hamilton の定理を明示的に適用している。これにより、進化演算子 $S(L) = \exp(-iH_f L)$ が $H_f$ の厳密な二次多項式であることを示している。この多項式の係数は、Vandermonde 系を通じて $H_f$ の3つの実固有値によって決定される。著者らは、この定式化が従来のニュートリノ文献では明示的に記述されていないことを指摘している。
2. スペクトル解と DMP 再総和化（Resummation）： 固有値は Cardano の三角公式を用いて導出される（Zaglauer–Schwarzer の式を回収する）。本論文は、Denton–Minakata–Parke (DMP) 近似を単なる摂動展開としてではなく、**再総和化（resummation）**として再定義し、分析する。1–3 セクターにおいて物質混合角 $\theta^M_{13}$ による厳密な回転を行うことで、DMP 法は、素朴な展開が $\hat{A} \to 1$ で発散する原因となる近接退縮を取り除いている。これにより、非有界な項 $(1-\hat{A})^{-1}$ を、あらゆるエネルギー領域で一様に有界な有効展開パラメータ $\epsilon_0 \lesssim 0.015$ に置き換えている。
3. 密度行列とデコヒーレンス： この定式化は、Lindblad マスター方程式と $su(3)$ における Bloch ベクトル表現を用いた開いた系の動力学へと拡張される。これにより、波束のデコヒーレンス補正、すなわち地上ベースラインにおける物質修飾されたコヒーレンス長（$L_{coh}$）の計算が可能となる。
4. 不確かさの伝播： 著者らは、6つの NuFIT 6.0 パラメータ（$\theta_{12}, \theta_{13}, \theta_{23}, \delta_{CP}, \Delta m^2_{21}, \Delta m^2_{31}$）に対する振動確率の閉形式のヤコビアンを導出している。これらは、物理的境界（例：$\theta_{13} \to 0$）や循環変数（$\delta_{CP}$）の近傍での妥当性を保証するため、Feldman–Cousins プロファイル尤度マッピングによって検証されている。

主な結果

• 進化演算子の構造： 本論文は、3つの固有値を二次進化演算子の係数に結びつける明示的な Vandermonde 系（式 15）を提供しており、反復的な行列指数関数計算を経ることなく、固有値から確率への直接的な代数的経路を提供している。
• DMP の精度とメカニズム： DMP 近似は、地上エネルギーにおいて $10^{-4}$ 未満の分数精度を達成することが示されている。本論文は、この精度が再総和化によるものであることを明らかにしており、その結果、展開パラメータが有界な $\epsilon_0$ となることを示している。表 6 では、展開パラメータが $O((1-\hat{A})^{-1})$ でスケールする他のスキーム（Madrid, AJLOS, MP）と対比させている。
• デコヒーレンスの大きさ： T2K/HK や DUNE といった地上長基線実験において、ベースラインとコヒーレンス長の比（$L/L_{coh}$）は $10^{-3}$ から $10^{-2}$ の範囲にあると計算されている。デコヒーレンスは現在の精度においては無視できる程度であるが、本論文は、それが文献で想定されているほど極端に微小なものではないことを指摘している。
• CP 非対称性の分解： 全体の CP 非対称性は、「真の」項（Jarlskog 不変量 $J_{CP}$ に依存）と「偽の」項（物質誘起、$\hat{A}$ に関して奇関数）に解析的に分割される。
• ヤコビアンのインフラストラクチャ： $\partial P_{\mu e}/\partial p_k$ に関する閉形式の表現（式 51–56）が提供されている。これにより、有限差分近似を用いることなく、単一の行列積による線形化された不確かさ伝播が可能となる。

意義と主張
本論文は、個々の結果（例：固有値解や DMP 近似そのもの）自体は新しいものではないと明言している。主要な貢献は、統合と代数的透明性にある：

• 統一された言語： Cayley–Hamilton 構造、DMP の再総和化解釈、および密度行列デコヒーレンスを、単一の整合的な枠組みへと結びつけている。
• 解析的な連続性： スピードを最優先する反復解法（NuFast-LBL など）とは異なり、本解析的枠組みは、パラメータの連続性、極限ケース（真空、高エネルギー、2世代）、および閉形式のヤコビアンを提供する。
• 解析における役割： 著者らは、本研究を生産的なフィット（スピードが最も重要となる場面）における NuFast-LBL アルゴリズムの代替品としてではなく、解釈、検証、および不確かさ伝播のためのツールとして位置づけている。これは、反復的な解法では構造上不可能な、限界の理解や相関のある不確かさの閉形式での伝播に必要な解析的インフラを提供するものである。
• 拡張性： 代数的な論理は、非標準相互作用（NSI）やステライルニュートリノのシナリオ（$4\times4$ 行列に対する三次 Cayley–Hamilton 多項式を必要とする）へ自然に拡張できることが示されている。

要約すれば、本論文は、進化演算子とその導関数を明示的にし、ニュートリノ振動物理学の精密測定時代における「帳簿管理」と代数的基礎を提供するものである。