NeutrinoOsc3Flavor: CP Phase Dependence in Three-Flavor Neutrino Oscillations: A Numerical Study in Vacuum and Matter

技術要約：NeutrinoOsc3Flavor

問題提起
現代のニュートリノ振動研究、特に物質効果やCP対称性の破れの位相依存性を扱う研究では、実験シミュレーション用に設計された大規模かつ高度に最適化されたソフトウェアフレームワークへの依存が強まっている。これらのツールはデータ解析には不可欠であるが、基礎となる数値的手順、中間ステップ、および特定の物理的仮定が複雑なソフトウェアスタックの中に埋め込まれており、しばしば「ブラックボックス」として機能してしまう。この不透明さは、独立した検証、数式レベルの検査、および理論的なクロスチェックを困難にしている。研究者がニュートリノ振動の計算の全工程を明示的に追跡できるような、透明性が高く、最小限の構成であり、検証を目的としたリファレンス実装への明確なニーズが存在する。

手法
本論文では、純粋なPythonで記述された、軽量で完全に透明な計算フレームワークであるNeutrinoOsc3Flavorを提示する。これは、真空および一定密度の物質中における、厳密な3世代ニュートリノ振動に関するシュレディンガー発展方程式を解くように設計されている。

• 定式化: コードは、フレーバー基底におけるPMNS定式化の範囲内で動作する。真空および物質の場合の有効ハミルトニアン（$H_f$）を明示的に構築する。物質については、電子密度に起因する荷電流ポテンシャル（$V_{CC}$）を組み込み、中性流効果（振動確率において相殺されるため）およびニュートリノ吸収は無視する。
• 数値的手法: 本フレームワークは、NumPyライブラリを用いて有効ハミルトニアンを数値的に対角化し、固有値と固有ベクトルを取得することで、進化演算子 $S(L) = e^{-iH L}$ を構築する。これにより、振動確率および生存確率を直接計算することが可能となる。
• 解析的検証: この実装のユニークな特徴は、**カルダノの方法（Cardano method）**を用いた、物質による修正ハミルトニアンの固有値の独立した解析的決定である。コードは、$3\times3$ ハミルトニアン行列の特性多項式を解くことで固有値を解析的に導出し、これを数値的な対角化結果に対するベンチマークとして提供する。
• 依存関係と移植性: 本ソフトウェアは、線形代数用のNumPyと可視化用のMatplotlibのみに依存している。外部のC++/FORTRANコンパイラ、特殊なライブラリ、またはプラットフォーム固有のインストールを必要とせず、Linux、Windows、macOSでのネイティブ実行を保証する。
• CP位相への感度: フレームワークは、CP対称性の破れの位相（$\delta_{CP}$）を主要な変数として扱い、その影響（例：$0^\circ, 45^\circ, 90^\circ, 135^\circ, 180^\circ, 217^\circ$）を系統的に変化させることで、真空および物質中における振動確率への影響を研究する。

主な貢献

1. 透明なリファレンス実装: コードは、ハミルトニアンの構築から最終的な確率評価に至るまで、ニュートリノ振動計算の全工程に対して数式レベルのアクセシビリティを提供しており、不透明な大規模ツールキットとは一線を画している。
2. 内部一貫性チェック: カルダノの方法による解析的な固有値決定を数値的な対角化と併用することで、本フレームワークは堅牢な内部検証メカニズムを提供する。著者らは、解析的な固有値と数値的な固有値の間の優れた一致を示している。
3. CP位相診断ツール: 本フレームワークは、CP位相を固定されたパラメータとして扱うのではなく、数値的安定性や物理的正当性のための感度の高い診断指標として使用するように特別に構造化されている。
4. 教育的およびベンチマーク的有用性: 本コードは、完全な実験シミュレーションフレームワークに取って代わるものではなく、理論的なクロスチェック、教育的使用、およびより複雑なシミュレーションソフトウェアのベンチマークとして設計されている。

結果

• 固有値の一致: 解析的に導出された固有値（カルダノの方法による）が、数値的に計算された固有値（NumPyによる対角化）と高い精度（例：$10^{-16}$ から $10^{-13}$ のオーダーの差）で一致することを確認し、数値ソルバーの安定性を検証した。
• 真空 vs 物質中の挙動:
  • 真空中では、振動確率はCP位相によって大きく変化するが、生存確率は示された特定の構成において位相の変化による影響をほとんど受けない。
  • 物質中では、振動確率と生存確率の両方がCP位相に対して感度を示し、物質効果によって、真空の場合には見られない生存確率のわずかな変化が誘起される。
• 検証: 物質中におけるミューオンニュートリノの生存確率の結果をNuOscProbExactツールキットと比較し、一貫した挙動を示した。

意義
著者らは、NeutrinoOsc3Flavorを包括的な実験シミュレーションパッケージの代替ではなく、検証を目的とした特化したツールとして位置付けている。その意義は、その解析的な追跡可能性とミニマリストな設計にあり、これらが独立した検証への障壁を下げている。計算の「スループット」よりも計算の「プロセス（どのように計算するか）」を優先することで、本フレームワークは以下に信頼できる基盤を提供する：

• 複雑な振動コードの理論的なクロスチェック。
• シュレディンガー発展の明示的なステップの理解が極めて重要となる教育的応用。
• 特に3世代系におけるCP位相依存性に関する数値的安定性のベンチマーク。

特殊な依存関係なしに動作する本フレームワークの能力は、異なるコンピューティング環境間での再現性を保証し、ニュートリノ物理学における既存の大規模ツールに対する価値ある補完物となっている。