NuFast-Earth: Efficient Atmospheric, Solar, and Supernova Neutrino Propagation Through the Earth

DUNE や HyperK などの次世代実験の需要に応えるため、地球を通過する大気・太陽・超新星ニュートリノの振動確率を、任意の地球モデルに対して高速かつ正確に計算できる新しいアルゴリズム「NuFast-Earth」とその C++ 実装が開発された。

要約

🌍 物語の舞台：地球という「迷路」

ニュートリノは、太陽や大気、あるいは超新星爆発（星の死）から飛んでくる、正体不明の「幽霊のような粒子」です。これらは物質をすり抜けるのが得意で、地球の中心まで突き抜けてしまいます。

しかし、ニュートリノは地球の中を歩くとき、**「物質の密度」**という壁にぶつかります。

• 地表付近はスカスカ（空気）。
• 地下は岩盤。
• 中心は鉄の核。

このように、地球の内部は**「密度が激しく変化する迷路」**です。ニュートリノはこの迷路を通過する間に、正体（種類）を変化させます（これを「振動」と呼びます）。

🧮 従来の問題：計算が「重すぎる」

これまで、この迷路を通過するニュートリノの動きをシミュレーションするには、**「超高性能なスーパーコンピュータを何日も稼働させる」ほどの計算量が必要でした。
なぜなら、地球の密度は一定ではなく、層ごとに細かく変化するからです。まるで、「1000 段ある階段の、1 段 1 段の硬さを測りながら、1 歩 1 歩の動きを計算し直す」**ような作業でした。

これでは、将来の巨大実験（DUNE や HyperK など）で得られる膨大なデータを解析する際に、計算が追いついてしまい、実験結果が出る前に時代遅れになってしまいます。

🚀 解決策：「NuFast-Earth」という新アルゴリズム

この論文の著者たちは、**「NuFast-Earth（ニュウ・ファスト・アース）」**という新しい計算プログラムを開発しました。これは、計算を劇的に速くする「魔法の道具」です。

1. 「賢いメモ帳」を使う（キャッシュの活用）

このプログラムの最大の特徴は、**「同じ計算を二度しない」**ことです。

• 従来の方法： 角度を変えて計算するたびに、地球の全層を最初から計算し直す。
• NuFast-Earth： 「あ、この層の計算は前もやったな！」と記憶しておき、必要な部分だけを再利用します。
  • 比喩： 料理を作る際、毎回野菜を全部切り直すのではなく、「切り終わった野菜は冷蔵庫（メモ帳）に入れておき、必要な分だけ取り出して使う」ようなものです。これにより、計算速度が10 倍〜100 倍に向上します。

2. 「回転」を後回しにする（数学的な工夫）

ニュートリノの動きを計算する際、複雑な「回転（角度）」の要素が含まれます。

• 従来の方法： 計算の最初から最後まで、この複雑な回転を常に持ち運ぶ。
• NuFast-Earth： 地球の中を通過する間の計算は、回転を一旦外して「シンプルで実数だけの計算」で行い、最後に回転を戻すという手順をとります。
  • 比喩： 重い荷物を運ぶ際、荷物の重さ（回転）を常に背負って歩くのではなく、**「荷物を下ろして軽装で走り抜け、ゴールの手前でだけ荷物を背負う」**という戦略です。これにより、計算中の負荷が大幅に減ります。

3. 「近似」の賢い使い方

正確な答えを出すには、非常に複雑な数式（カルダノの公式など）を使いますが、これは計算に時間がかかります。
著者たちは、**「ほぼ正確で、かつ圧倒的に速い近似式」**を見つけ出し、それを採用しました。

• 比喩： 目的地までの距離を測る際、「1 ミリ単位まで正確に測る」のではなく、「10 秒で測れるが、誤差は髪の毛の太さ程度」という方法を使います。実験の精度には十分すぎるほど正確なのに、時間は 1/3 で済みます。

🌟 この技術がなぜ重要なのか？

この新しいプログラムを使えば、以下のことが可能になります。

1. 地球の内部を「透視」できる：
   ニュートリノの動きを精密に解析することで、地球の中心の核がどんな構造をしているか、地震学では見えない部分まで詳しく調べられるようになります（ニュートリノ・トモグラフィー）。
2. 太陽と超新星の謎を解く：
   太陽から来るニュートリノや、遠くの星の爆発から来るニュートリノが、地球を通過する際にどう変化するかを正確に予測できます。これにより、宇宙の秘密がより深く理解できます。
3. 実験の未来：
   次世代の実験装置が稼働すると、膨大なデータが生まれます。NuFast-Earth は、そのデータを瞬時に処理し、人類がまだ知らない「ニュートリノの性質」や「宇宙の法則」を見つけ出すための鍵となります。

📝 まとめ

この論文は、**「複雑怪奇な地球の迷路を、ニュートリノがどう通り抜けるかを計算する」という難問に対して、「賢いメモ帳」「軽装で走る」「速くて正確な近似」という 3 つのアイデアで、計算を「爆発的に速く」**したという報告です。

まるで、**「何日もかかる大工事だった計算を、数秒で終わらせる魔法のハンマー」**を手にしたようなものです。これにより、科学者たちはこれまで不可能だった精密な実験を行い、宇宙と地球の深淵な秘密に迫ることができるのです。

技術概要

以下は、提示された論文「NuFast-Earth: Efficient Atmospheric, Solar, and Supernova Neutrino Propagation Through the Earth」の技術的な要約です。

論文概要

タイトル: NuFast-Earth: Efficient Atmospheric, Solar, and Supernova Neutrino Propagation Through the Earth
著者: Peter B. Denton, Stephen J. Parke
日付: 2026 年 3 月 10 日（予稿）

この論文は、地球内部を通過する大気ニュートリノ、夜間の太陽ニュートリノ、および超新星ニュートリノの振動確率を、極めて高速かつ高精度に計算するための新しいアルゴリズムと C++ コード「NuFast-Earth」を提案するものです。次世代の実験（DUNE, HyperK, IceCube アップグレード, KM3NeT など）において、大量のイベント生成（モンテカルロ・シミュレーション）が必要となる中、計算コストの削減が急務となっています。

---

1. 背景と課題 (Problem)

ニュートリノ振動の確率計算において、物質効果（Wolfenstein 効果）は重要な要素ですが、地球のような急激に変化する密度プロファイル（マントル、外核、内核など）を通過する場合は計算が複雑になります。

• 計算コストの増大: 次世代実験では、Feldman-Cousins 法などの統計的手法を用いてイベント率を精密に予測するために、膨大な数の「投げ（throws）」が必要となります。各投げごとに振動確率を再計算する必要があり、これが計算のボトルネックとなっています。
• 既存手法の限界: 長基線加速器ニュートリノ用として開発された「NuFast-LBL」アルゴリズムは、一定密度の近似が効く場合や単一の距離に限定される場合は高速ですが、地球内部の急激な密度変化を扱う大気ニュートリノや、地球を通過する太陽・超新星ニュートリノには適していません。
• 精度と速度のトレードオフ: 従来の数値積分法や行列指数関数の直接計算は正確ですが、計算量が膨大です。一方、近似手法は速いですが、必要な精度を維持できるかが課題でした。

2. 手法とアルゴリズム (Methodology)

NuFast-Earth は、NuFast-LBL の要素を流用しつつ、地球の密度プロファイルに対応するためにゼロから構築された新しいアプローチを採用しています。

A. 計算の基盤となる「tilde 基底」

• 実数対称行列の活用: ハミルトニアンの計算において、$\theta_{23}$ と CP 位相 $\delta$ を含む回転行列を分離し、残りの部分（$\tilde{H}$）を実数対称行列として扱います。これにより、固有値と固有ベクトルがすべて実数となり、計算が大幅に簡略化されます。
• 固有値の高速計算: 特性方程式の解（固有値）を、カルダノの公式（厳密解）ではなく、優れた近似式からニュートン・ラフソン法（NR 法）を用いて反復計算することで求めます。
  • 結果: 2 回の NR 反復で倍精度浮動小数点数の精度に達し、カルダノ法よりも約 3 倍高速です。さらに、1 回の反復でも次世代実験に必要な精度を十分に満たします。
• 固有ベクトルの直接計算: 固有ベクトルを直接対角化するのではなく、固有値とハミルトニアンの随伴行列（Adjugate）の関係式を用いて効率的に計算します（Eigenvectors from Eigenvalues の手法）。

B. 地球通過経路の分割と対称性の利用

ニュートリノの軌道を 4 つのセグメントに分割し、計算を最適化します（図 1 参照）：

1. A (Air): 生成点から地表まで（真空中）。
2. S (Surface): 地表から検出器深度まで（反対側）。
3. I (In): 検出器深度から地球中心に最も近い点まで。
4. O (Out): 中心に最も近い点から検出器まで。

• 対称性の利用: 球対称な地球モデルにおいて、セグメント I と O の密度プロファイルは逆順で同一であるため、$\tilde{A}_O = \tilde{A}_I^T$（転置行列）の関係が成り立ちます。これにより、行列乗算の回数を約半分に削減できます。
• 大気層の扱い: 大気層は真空中とみなせるため、計算が容易で、生成高さの変更も容易です。

C. 確率への変換と「高速パラメータ」

• 地球通過の振幅計算の大部分は $\theta_{23}$ と $\delta$ に依存しない「tilde 基底」で行われます。
• 最終的に、$\theta_{23}$ と $\delta$ の回転を適用してフレーバー基底に戻し、確率を計算します。
• 高速パラメータ: 混合角 $\theta_{23}$、CP 位相 $\delta$、生成高さ、太陽密度などを変化させる場合、地球内部の複雑な計算（固有値・固有ベクトル・行列積）を再利用できるため、計算コストが極めて低くなります。

3. 主要な貢献 (Key Contributions)

1. NuFast-Earth コードの公開: GitHub で公開された C++ 実装。ユーザーマニュアル付き。
2. 計算速度の劇的な向上:
   • 既存のコード（nuSQuIDS, OscProb など）と比較して、特に「高速パラメータ」をループさせる場合、計算時間が桁違いに短縮されます（例：OscProb の約 40,000 ns に対し、NuFast-Earth は約 200 ns）。
   • 地球モデルの複雑さ（層の数）が増しても、$\theta_{23}$ や $\delta$ の変化に対する計算コストは増加しません。
3. 柔軟な地球モデルのサポート:
   • PREM（Preliminary Reference Earth Model）に基づく多層モデル、定数密度モデル、積分ベースモデルなど、任意の球対称密度プロファイルを扱えます。
   • 層の分割数をユーザーが制御でき、精度と速度のバランスを調整可能です。
4. 多様なニュートリノ源への対応:
   • 大気ニュートリノ（昼・夜）、太陽ニュートリノ（昼・夜）、銀河超新星ニュートリノのすべてを統一的な枠組みで計算可能。

4. 結果 (Results)

• 速度: コンパイラ最適化（-O3）を用いた場合、固有値・固有ベクトルの計算は厳密解（カルダノ法）よりも 3 倍以上高速です。1 つの軌道（コア通過など）の確率計算は、パラメータのループ条件により数 100 ns から数 $\mu$s のオーダーで完了します。
• 精度:
  • 地球モデルを 1000 層程度に分割し、固有値の NR 反復を 1 回行う設定では、平均誤差が $10^{-6}$程度、最大誤差が$10^{-3}$ 未満となり、次世代実験の要求精度を十分に満たします。
  • 検出器の深度（地表 vs 2km 深）による確率の差異は、低エネルギー領域や水平方向の軌道で最大 5% 程度現れることが確認され、高精度計算には深度の考慮が必須であることが示されました。
• 検証: 太陽ニュートリノの昼間・夜間の振動確率、DUNE のような長基線実験のシミュレーション結果など、既知の解析解や他のコードとの比較により、コードの正当性が検証されています。

5. 意義と結論 (Significance)

• 次世代実験への不可欠なツール: DUNE や HyperK などの次世代実験では、統計的な精度を高めるために膨大な数のシミュレーションが必要です。NuFast-Earth は、これらの計算を現実的な時間枠内で実行可能にするための重要なインフラとなります。
• 地球トモグラフィーへの貢献: 大気ニュートリノを用いた地球内部構造の探査（ニュートリノ・トモグラフィー）において、密度プロファイルの微細な変化を捉えるための高精度かつ高速な計算基盤を提供します。
• 計算手法の革新: 固有値・固有ベクトルの近似計算と、対称性を利用した行列演算の最適化、そして「再利用可能な計算」の概念を組み合わせることで、ニュートリノ振動計算のパラダイムシフトを促しています。

結論として、NuFast-Earth は、複雑な物質効果を持つ環境におけるニュートリノ振動確率計算において、**「速度」「精度」「柔軟性」**を同時に実現した画期的なツールであり、将来のニュートリノ物理学の進展に不可欠な基盤技術です。