Bivariate incomplete-Bessel kernels for the first nonlinear Vlasov-Maxwell response

技術要約：第1次非線形ブラーソ・マックスウェル応答における二変量不完全ベッセル核

問題の定義
磁化プラズマにおける非線形波結合の運動論的理論は、伝統的に、均質な平衡状態の周りでの弱非線形展開に依存している。標準的な手法では、プラズマ応答を二重サイクロトロン調和展開として表現する。この表現は明示的ではあるものの、電磁波が背景磁場に対して任意の角度で伝播する場合、代数的に非常に煩雑になる。第1次の非線形補正において、応答は二つの結合した調和指数を伴う。一つは出力モードの外的な特性積分から生じるものであり、もう一つはソースとして作用する内的な線形応答から生じるものである。この構造は、LiuおよびTripathi (1986) やその後の研究に見られるような、ベッセル関数の複雑な二重和と入れ子状の共鳴分母をもたらす。

手法
本論文では、ラーモア位相を調和級数へと即座に展開することを避け、代わりに位相を摂動のない軌道に保持し、特性積分を直接評価する代替的な定式化を提案する。

1. 線形ベースライン： 著者らは、特性積分が不完全ベッセル関数 $G_\mu(z, \psi)$ によって評価される最近の線形定式化（Ricci 2026）を回顧する。この関数は、感受率テンソルが投影される前に、無限のサイクロトロン調和級数の代わりを務める。
2. 非線形拡張： 第1次の非線形ブラーソ方程式は、同じ特性法を用いて解かれる。非線形ソース項は、あるモードのローレンツ力が別のモードの線形応答に作用するものとして構成される。内的な線形応答を不完全ベッセル形式で表現することで、外的な特性積分はシフトされた $G_{\mu_q}$ 応答に対して作用することになる。
3. 核（カーネル）の導出： このプロセスは自然に、二変量不完全ベッセル関数（$G^{(r)}_{\mu,\nu}(z, \psi; w, \chi)$ と表記）を生成する。この関数は、横方向の速度因子やジャイロ位相微分によって生じる有限の調和シフトを整数 $r$ で考慮した、軌道分解能核として機能する。
4. 古典的な結果の回収： 著者らは、これらの二変量関数を用いて非線形分布関数を明示的に導出する。その後、これらの関数をヤコビ・アンガー公式およびTurkinの関数を用いて展開することで、伝統的な二重調和のLiu–Tripathi公式が回収されることを示し、この新しい定式化の妥当性を検証する。
5. 感受率の投影： 最後に、二変量分布を速度モーメントへと投影することで、非線形感受率テンソルを導出する。これは、二変量関数の角度フーリエ投影である「投影された二変量核」を通じて達成され、完全な二重調和展開に戻る必要を回避する。

主な貢献と結果

• 二変量核の定義： 本論文は、$G^{(r)}_{\mu,\nu}$ を第1次非線形特性問題から生成される自然な軌道積分として導入する。この関数は単なる表記上の便宜ではなく、特性法の反復によって生じる基礎的な対象であることを確立している。
• 構造的再編成： 非線形分布関数は、これらの二変量核に作用する微分演算子を用いてコンパクトに表現される。この定式化は、角度の縮退を、再利用可能な一連の特殊関数ブロックへと孤立させている。
• 数学的恒等式： 二変量関数に関する包括的な恒等式のセット（重ね合わせの性質、調和展開、漸化式、および微分方程式（外側ベッセル方程式および内側Turkin漸化式を含む）を含む）が導出されている。
• 整合性の検証： 著者らは、二変量定式化が古典的な二重調和アプローチと数学的に等価であることを示す厳密な導出を提供している。中間調和指数を縮約するために、二変量関数を展開し、Grafの加法定理を適用することで、伝統的な公式が回収される。
• 明示的な感受率の公式： 非線形感受率テンソルは、9つの基本二変量ブロック（$a, b \in \{-1, 0, 1\}$ でインデックス付けされる）に対する有限和と、それに続く二次元速度積分として提示される。角度依存性は、投影された核 $K^{(r,a)}_N$ およびその微分によってカプセル化されている。

意義と主張
本論文は、この定式化の主な利点は構造的、かつ潜在的に計算的であると主張している：

• 構造的： これは、第1次の非線形サイクロトロン応答の基礎となる特定の特殊関数オブジェクト（二変量不完全ベッセル核）を特定し、調和和の代数的な複雑さではなく、特性流の幾何学に基づいて理論を整理するものである。
• 計算的： 著者らは、数値実装において、二変量定式化がより構造化された対象を提供することを示唆している。シフトされた指数のベッセル関数の長い振動的な二重和を扱う代わりに、投影された核とその微分をテーブル化して再利用することができる。この柔軟性は、直接的な調和の総和が収束しにくい、あるいは打ち消し合いのために数値的にデリケートな領域において有利となる可能性がある。
• スケーラビリティ： 著者らは、このアプローチが再帰的に高次へと拡張可能であると述べている。$n$ 次の分布は、$n$ 変量軌道核を用いて表現できると予想され、これは非線形相互作用ごとに一つの特性時間積分を加えることで生成される。

本論文は、数値的な利点が普遍的な計算コストの削減ではないことを認め、プラズマパラメータ（ラーモア半径、デチューニング、減衰）に依存することを明記しており、控えめな立場を維持している。また、この構成は固定特性摂動階層に適用されるものであり、特性演算自体を修正するような非線形軌道補正（トラッピングやポンデロモーティブ・ドリフトなど）を考慮していないことを明示している。