Improving divergence cleaning in cosmological SPMHD simulations

技術要約：宇宙論的SPMHDシミュレーションにおけるダイバージェンス・クリーニングの改善

問題提起
宇宙論的な平滑化粒子磁気流体力学（SPMHD）シミュレーションにおいて、ソレノイド条件（$\nabla \cdot \mathbf{B} = 0$）を維持することは、持続的な数値的課題である。Gadgetなどのコードで採用されている標準的なスキームは、通常、Powellら（1999）による「8波（8-wave）」アプローチに依存している。この手法はダイバージェンス誤差を流体の流れとともに移流させるが、それらを能動的に除去するものではない。SPHベースのコードでは、残留する $\nabla \cdot \mathbf{B}$ 誤差がオーダー1に達することがあり、これが非物理的な磁場増幅や、磁気圧がガス圧を上回る際に運動方程式に偽のソース項が支配的となる「テンシル不安定性（tensile instability）」を引き起こす。Dednerら（2002）によって提案された双曲型／放物線型ダイバージェンス・クリーニング・スキームは解決策を提供しているが、そのオリジナルの定式化は、動的で個別のタイムステップを用いる環境、特に強い密度ジャンプが存在する環境下のSPMHDにおいては不安定であることが判明しており、非物理的な磁気エネルギーの増大源として作用することがある。

手法
著者らは、制約付きかつエネルギー保存型のDedner型クリーニング・アルゴリズム（Tricco & Price 2012; Tricco et al. 2016a）を、宇宙論的SPMHDコードであるOpenGadget3に実装した。主な手法的な適応は以下の通りである：

• エネルギー保存： 全エネルギー保存を保証するために、支配方程式を再定式化した。これには、補助的なスカラー $\psi$ 単体ではなく $\psi/c_h$（ここで $c_h$ はクリーニング速度）としてクリーニング変数場をエボリューションさせ、進化方程式に特定のソース項 $-\frac{1}{2}\psi(\nabla \cdot \mathbf{v})$ を含めることが含まれる。これにより、クリーニング場の放物線的な減衰が、無制限な成長を招くことなく磁気エネルギーを除去することを保証している。
• 宇宙論的適応： 方程式を共動座標および個別の粒子タイムステップに適応させた。クリーニング速度 $c_h$ は局所的な高速磁気音速（$\sqrt{v_A^2 + c_s^2}$）に設定されており、これによりクリーニングが標準的なSPHタイムステップの制約内で本質的にCourant安定であることを保証している。
• シミュレーション設定： 著者らは、一様な磁場（共動座標で$10^{-14}$ G）が埋め込まれた巨大銀河団（$M_{200c} \approx 10^{15} M_\odot$）を用いた一連のシミュレーションを実施した。5つの異なる実行（ラン）を比較している：
  1. NIFTY-noclean: Powellの移流スキームのみを使用したリファレンス・ラン。
  2. NIFTY-clean: エントロピー・バックリアクションを含む、フル機能の制約付きクリーニング・スキームを用いた標準的なラン。
  3. NIFTY-clean2: 上記と同様だが、オーバークリーニング係数 $\zeta=2$（局所的な波速の2倍でクリーニング）を用いたもの。
  4. NIFTY-advect: 速度ダイバージェンスのソース項をオフにしたクリーニング・スキーム。
  5. NIFTY-heat: 散逸したクリーニングエネルギーをエントロピー方程式にフィードバックせず、破棄するクリーニング・スキーム。

主な結果

• ダイバージェンス誤差の低減： 主要な結果として、Powellのみのスキームと比較して、クラスター全域における相対的なダイバージェンス誤差（$|\nabla \cdot \mathbf{B}|/B$）が2〜3桁減少した。この減少は、クラスターのコア部と外縁部の両方で一様に観察された。
• 磁場増幅： 磁場が弱く解像度が低いクラスターの外縁部（$r \approx 1\text{--}3 \, h^{-1}$ Mpc）において、能動的なクリーニング・スキームは、Powellのみの実行と比較して5〜10倍の磁場増幅を可能にした。著者らは、これを、以前はダイバージェンス誤差による非物理的なローレンツ力が、これらの領域における小規模ダイナモを抑制していたためであると説明している。
• コアの安定性： クラスターコアにおける全域的な密度構造と磁場幾何学は、クリーニングを行う実行と行わない実行の間でほとんど変化しておらず、これはクリーニングがコアのダイナモの飽和状態を人工的に変えていないことを示している。
• 熱力学的整合性： エントロピー・バックリアクション（散逸したクリーニングエネルギーを熱として戻すこと）の包含は、熱力学的整合性を保つために極めて重要であることが示された。このバックリアクションがない実行（NIFTY-heat）では、コアの断熱圧縮により、中心部のエントロピーが低くなり、中心部の飽和磁場が高くなる結果となった。
• 計算コスト： この実装による実行時間のオーバーヘッドは計測可能な範囲では存在しない。改良されたクリーニングは、既存のCourantタイムステップの制約内で実質的に「無料」で達成されている。

意義と主張
本論文は、ダイバージェンス誤差を正確に制御することが、銀河団の低密度で解像度の低い外縁部における磁場増幅を捉えるために不可避であると主張している。著者らは、標準的なPowellスキームにおける制御不能なダイバージェンス誤差が、サチュレーション（飽和）前の領域におけるダイナモを抑制する数値的ノイズとして作用していると述べている。制約付きのエネルギー保存型クリーニング・スキームを実装することで、シミュレーションは計算上のペナルティなしに、より忠実なMHD物理を実現している。

著者らは、制約付きクリーニング・スキームについて、以下の結論を導いている：

1. 数値的なダイバージェンス誤差を劇的に減少させる（2〜3桁）。
2. 解像度の低い領域における磁場増幅のより正確な表現を可能にする。
3. より滑らかでコヒーレントな磁場構造を生み出す。
4. 異なる解像度間での磁場進化の安定性を高め、中程度の解像度であっても物理的な忠実度を向上させる。

本研究は、宇宙論的なSPMHDアプリケーションに対し、自然なクリーニング速度（$\zeta=1$）を使用し、エネルギー収支を閉じるために放物線的なエネルギー・バックリアクションを有効にした状態で、このクリーニング・スキームをデフォルトとして採用することを推奨している。この実装は現在、OpenGadget3において公開されている。