DIPLODOCUS I: Framework for the evaluation of relativistic transport equations with continuous forcing and discrete particle interactions

技術要約：DIPLODOCUS I：連続的な強制力および離散的な粒子相互作用を伴う相対論的輸送方程式の評価のためのフレームワーク

問題提起
ブラザール・ジェットのような高エネルギー現象の現在の天体物理学的モデリングは、マクロな流体力学モデル（GRMHD）またはミクロな運動論モデルのいずれかに大きく依存している。GRMHDは局所熱平衡と等方性を仮定しており、詳細な運動量空間分布を無視するため、放射スペクトルを直接シミュレートすることには失敗する。対照的に、粒子・イン・セル（PIC）シミュレーションのような微視的なアプローチは、計算上の制約により有限の粒子数と小さな空間スケールに限定されており、マクロな適用が困難である。さらに、従来の運動論モデルは、次元を削減するために、運動量空間の等方性や空間の均質性（単一／マルチゾーンモデル）といった厳格な対称性を課すことが多い。これらの仮定は、決定的な物理的洞察を覆い隠す可能性がある。特に、近年の研究によって、等方的なモデルでは説明できないことが示唆されている、放射における異方的な粒子進化や空間的な変位（例：Terzan 5のガンマ線異常）に関する問題である。したがって、制限的な対称性の仮定に頼ることなく、連続的な力と離散的な相互作用を扱うことができる、メゾスコピックなフレームワークが必要とされている。

手法
本論文では、天体物理学的システムのメゾスコピックなモデリングのために設計されたフレームワークであるDIPLODOCUS（Distribution-In-PLateaux methODOlogy for the CompUtation of transport equationS）を紹介する。本手法は、主に以下の3つの段階で構成される：

1. 座標不変な定式化： 著者らは、位相空間を表す微分可能多様体 $M$ 上における、相対論的輸送方程式（ボルツマン方程式）の厳密かつ座標不変な共変的定式化を確立した。このアプローチは、粒子輸送を、質量による制約を含む7次元位相空間（4次元時空 ＋ 3次元運動量）を通る世界線の流れとして扱う。本フレームワークは、リー微分を利用して位相空間の体積形式の進化を記述し、保存的な力（ハミルトニアン的）と非保存的な力（非ハミルトニアン的、例えば放射反作用など）の両方を明示的に考慮する。
2. 座標依存の実装： 天体物理学的モデリングへの適用を可能にするため、本フレームワークは特定の座標系と基底ベクトル（テトラド）を導入する。曲がった時空の接続係数と任意の四元力を組み込み、制限された7次元多様体 $M_m$ 上の位相流ベクトル $v$ を定義する。また、静止（エウラー的）観測者のための体積形式および超曲面要素を導出し、同時事象と非同時事象を区別する。
3. 衝突積分と離散化： 本フレームワークは、離散的な相互作用（衝突、散乱、崩壊）を、世界線の終端と開始をカウントする衝突積分 $C(x, p)$ として組み込んでいる。この積分は、遷移確率と散乱行列から導出される。高次元の積分を計算可能にするため、本論文ではDistribution-In-Plateaux (DIP) 法を提案する。これは、位相空間の体積を共通の境界を持つ連続的な部分領域（プラトー）へと離散化するものである。この離散化により、積分輸送方程式を、境界を横切るフラックスと内部の衝突が均衡する保存的な数値スキームへと変換し、等方性を仮定することなく異方的な相互作用の評価を可能にする。

主な貢献

• 一般化された輸送フレームワーク： 連続的な強制力（ローレンツ力、放射反作用など）と離散的な粒子相互作用を、単一の共変的な枠組みの中で統一する、包括的な相対論的輸送方程式の導出を提供する。
• 非保存的な力の取り扱い： 保存的な系のみに焦点を当てた従来の多くの導出とは異なり、本研究は $d\omega \neq 0$ となる非保存的な力を含むように、リウヴィルの定理を明示的に一般化しており、これは天体プラズマにおける放射損失のモデリングにおいて極めて重要である。
• 対称性の制約のない異方性： DIP手法は、運動量空間の等方性や空間の均質性を仮定せずに輸送を評価するように設計されている。これにより、異方的な粒子分布とその進化のモデリングが可能となり、現在の単一／マルチゾーンモデルの限界を解消する。
• 保存的な数値スキーム： 離散化の手順は、位相空間ドメイン全体における粒子数（世界線）の保存を保証するものであり、これは堅牢な長期シミュレーションに不可なく欠な条件である。

結果と範囲
本論文はシリーズの第一報であり、数値シミュレーションや特定の天体物理学的ケーススタディは提示しない。代わりに、「結果」は、解析的および数値的なフレームワーク自体の導出と定義である。本論文は、7次元位相空間を通じて粒子分布関数を進化させるために必要な数学的バックボウンを確立する。また、異方性と非保存的な力の導入が、数学的に一貫しており、現代のコンピューティング能力の範囲内で計算可能であることを示しており、その証拠として、対生成プラズマ風や宇宙線伝播における軸対称分布のモデリングにおける最近の成功を挙げている。

意義
著者らは、DIPLODOCUSが、現在のゾーンモデルに内在する等方性や均質性の仮定を緩和する必要性に対処するものであると主張している。粒子輸送の全容を扱うことができるフレームワークを提供することで、粒子加速の起源やジェット物質の組成（レプトン的かハドロン的か）といった、天体データの解釈における曖昧さを解決することを目指している。本フレームワークは、ジェット源を超えた幅広い天体現象に適用可能な一般的なツールとして提示されている。本論文は、数値実装、微視的能力（相互作用、放射、強制力）、および巨視的能力（空間構造と移流）を詳述する後続の著作（論文IIおよびIII）への舞台を整えるものである。