A novel objective function minimizes resonant trapped energetic particle losses in stellarators

技術要約：共鳴による捕捉エネルギー粒子の損失を最小化する新しい目的関数を用いたステラレータの設計

問題提起
準オムニジェナス（near-omnigenous）なステラレータは、捕捉されたエネルギー粒子（EP）と磁場の非オムニジェナスな摂動との間の共鳴に起因するEP損失に対して脆弱である。既存のバウンス平均化された最適化目的関数（$\Gamma_c$、$\Gamma_\delta$、$\Gamma_\alpha$ など）は、ドリフト面と磁束面の非共鳴的な不一致を効果的に標的とするが、共鳴的な対流または拡散運動を捉えることはできない。これらの共鳴は、EPの歳差運動周波数とバウンス周波数の比が有理数となる場合に発生し、位相空間におけるアイランドの形成を招く。単一のアイランドは指向性のある「対流的」な径方向運動を生み出し、アイランドの重なりはカオス的な「拡散的」な運動を生み出す。両方のメカニズムが、特に共鳴の重なりが一般的な準軸対称（QA）構成や、ゼロ歳差運動共鳴が支配的な準ヘリカル（QH）および準アイソダイナミック（QI）構成において、プロンプトなEP損失に大きく寄与する。ガイドセンター・トレーシングに依存する既存の直接最適化（DO）手法は、計算コストが高く微分不可能であるため、勾配ベースの最適化ループにおける利用が制限されている。

手法
著者らは、準オムニジェナス磁場における捕捉されたEPのバウンス点に関する離散写像理論を開発した。作用量のような座標 $s$ と、角度のような磁力線ラベル $\eta$ を用いた座標系 $(s, \eta)$ を定義することで、捕捉粒子の運動をハミルトニアン離散写像として記述する。この枠組みにおいて、共鳴は歳差運動周波数 $\Omega_\eta$ が有理数 $p/q$ に等しいときに発生する。理論では、磁気アイランド幅の計算と同様に、これらの共鳴面を中心とした位相空間アイランドの幅（$\Delta s_W$）の式を導出している。

この理論に基づき、著者らは新しく微分可能なバウンス平均化された目的関数 $\Delta_{res}$ を導入した。$\Delta_{res}$ は、以下の手順で位相空間アイランドの幅を罰則化する：

1. $\Omega_\eta$ が有理数 ($p/q$) を横切る面を特定する。
2. 離散写像理論を用いて、アイランド幅のプロキシ（代理指標）を計算する。これは、ハミルトニアンの共鳴成分（$|f_r|$）と歳差運動周波数の勾配（$\Omega'_\eta$）に依存する。
3. 有理数横切りの近傍にある格子点を罰するために、径方向の重み関数を適用する。
4. 捕捉粒子に対する位相空間平均を実行し、磁束面あたりのスカラー指標を得る。

最適化中の数値的安定性を確保するため、実装には、準対称構成で頻繁に発生する $\Omega'_\eta(s)$ がゼロを通過する際の特異性を打ち消すための修正が含まれている。著者らは、DESCコードを利用して $\Delta_{res}$ を実装し、微分ベースの最適化を行った。

主な貢献

• 離散写像理論： 捕捉されたEPのバウンス平均化された離散写像理論を確立し、歳差運動周波数とバウンス周波数の間の共鳴から生じる位相空間ドリフトアイランドを定量化した。
• $\Delta_{res}$ 目的関数： 共鳴アイランドの幅を直接的に標的とし、位相空間の可積分性を促進する、計算効率が高く微分可能な指標 $\Delta_{res}$ の開発。
• バウンス平均化の検証： 有限軌道幅（FOW）効果がバウンス平均化と完全なガイドセンター・トレーシングの間の乖離を引き起こす可能性があるにもかかわらず、アイランド幅のバウンス平均化計算が共鳴位置とサイズの定性的に正確な予測を提供することを示し、その使用の妥当性を証明した。
• 最適化フレームワーク： 標準的な準対称性指標とともに、$\Delta_{res}$ をマルチオブジェクティブ最適化フレームワークに統合した。

結果
著者らは、準軸対称（QA）ステラレータ構成の最適化に $\Delta_{res}$ を適用した。この最適化は、$\Delta_{res}$ と2項の準対称性目的関数を組み合わせて行われた。

• 損失低減： 最適化された構成は、EP閉じ込めにおいて4倍の改善を達成した。3.5 MeVのアルファ粒子について、$\Delta_{res}$ なしの場合の損失率2.04%に対し、$\Delta_{res}$ ありの場合は0.57%（0.1秒間）へと減少した。
• 改善のメカニズム： 最適化は、主に2つのメカニズムを通じて閉じ込めを向上させた：
  1. 共鳴の移動： 低次の共鳴をプラズマ境界から遠ざけた。
  2. 輸送障壁の形成： 最適化された構成は輸送障壁（ドリフト面がよく閉じ込められた領域）を形成し、拡散領域が境界へと伝播するのを防いだ。
• 相関： 本研究は、特に共鳴的ドリフト対流または拡散が支配的な領域において、$\Delta_{res}$ 指標と実際のEP損失率との間に強い相関があることを示している。ポアンカレ・プロットにより、損失の減少が大きなドリフトアイランドの抑制およびコヒーレントなドリフト面の維持に対応していることが確認された。

意義
本論文は、$\Delta_{res}$ が、発電所レベルのステラレータにおける共鳴対流および拡散の両方の損失に対抗するための強力なツールであることを主張している。共鳴挙動を特異的に標的とする微分可能な指標を提供することで、$\Delta_{res}$ は、直接的なガイドセンター・トレーシングに伴う膨大な計算コストを回避しながら、勾配ベースの効率的な最適化を用いてEP損失を最小化することを可能にする。著者らは、$\Delta_{res}$ が（共鳴の重なりが主要な問題となる）QA構成において特に効果的であると述べているが、今後は、共鳴的対流アイランドが支配的なQIおよびQH構成における有効性についても検討する予定である。本手法は、EPの位相空間の可積分性を高めることにより、堅牢なEP閉じ込めを持つステラレータを設計するための道筋を提供するものである。