Charging of rotating black holes: kinetic simulations of black hole magnetospheres

技術要約：回転ブラックホールの運動論的シミュレーションによる帯電

問題提起
本論文は、ブラックホール（BH）電磁力学における長年の未解決問題、すなわち、外部磁場に浸された回転ブラックホールが、ウォルドの平衡値 $Q_W = 2aMB$ （$M$ は質量、$a$ はスピンパラメータ、$B$ は磁場強度）まで帯電するかどうかという問題に取り組んでいる。ウォルドの解は、地平線と無限大の間の電位差が消失する平衡状態としての $Q_W$ を示唆しているが、これまでの理論的アプローチでは相反する予測が示されてきた。エネルギー最小化や選択的な降着に基づくモデルは、ウォルド値よりも低い平衡状態を示唆し、一方で完全な遮蔽（$Q \approx Q_W$）や振動挙動を主張するモデルもある。決定的なのは、既存の力学的なフォースフリー（force-free）モデルや理想MHDモデルでは、電荷分離に必要な平行電場を取り除いてしまう $\vec{E} \cdot \vec{B} = 0$ を強制するため、この問題を扱うことができない点である。中心となる問いは、自己整合的な対生成を伴うプラズマ充填された磁気圏内に置かれた天体物理学的ブラックホールが、$Q_W$ へと向かうのか、あるいは別の値へと進化するのか、そしてその最終状態が初期条件に依存するかどうかである。

手法
著者らは、GRZeltron コードを用いた、第一原理に基づく軸対称一般相対論的粒子・イン・セル（GRPIC）シミュレーションを用いている。このアプローチは、プラズマの導電率や供給メカニズムに関する仮定を排除している。

• セットアップ: シミュレーションでは、様々なスピンパラメータ（$a = 0.2$ から $0.95$）に対して、漸近的に一様なウォルド磁場をカー・ブラックホールの周囲に進化させる。
• 初期条件: 各スピンに対して、2つの異なる初期状態、すなわち中性ブラックホール（$Q_0 = 0$）と、ウォルド値まで事前帯電されたブラックホール（$Q_0 = Q_W$）をテストする。
• 含まれる物理: シミュレーションは、ソフトフォトンの逆コンプトン散乱による対生成、放射反作用、およびプラズマの完全な運動論的応答を自己整合的にモデル化している。
• 診断: 蓄積された電荷は、ガウスの法則によるイベントホライゾン（事象の地平線）を貫くフラックスを通じて測定され、径方向の電流フラックスの積分によって相互検証される。系は飽和状態に達するまで、約160重力時間（$t_g$）にわたって進化させる。

主要な貢献と解析的枠組み
本論文は、ウォルド問題における平衡ブラックホール電荷の、初の直接的なGRPICによる決定を提供している。さらに、著者らは磁束バランスの基準に基づき、平衡電荷（$\xi_{eq} \equiv Q_{eq}/Q_W$）の解析的な式を導出している。

• 解析的導出: 著者らは、正電荷が捕捉される極冠（polar cap）を貫く磁束が、負電荷が捕捉される赤道ベルト（equatorial belt）を貫く磁束と等しくなったときに平衡に達すると主張している。これには、$\vec{E} \cdot \vec{B} = 0$ によって定義されるセパラトリクス（分離線）が、地平線の磁束を二等分する必要がある。
• 導出された公式: これにより、スピン $a$ に依存する $\xi_{eq}$ の閉じた形式の式（式22）が得られる。このモデルは、低回転限界において $\xi_{eq} \approx 0.25$ となり、メイスナー効果（Meissner effect）による磁束の地平線からの排出により、$a \to 1$ に伴ってゼロに向かって減少することを予測している。

結果

1. 動力学的アトラクター: シミュレーションは、系が初期の電荷状態の記憶を失うことを示している。$Q_0 = 0$ と $Q_0 = Q_W$ の両方の分岐は、数十重力時間以内に同じ正の飽和電荷値へと収束する。これは、平衡電荷が運動論的磁気圏の動力学的アトラクターであることを裏付けている。
2. サブ・ウォルド平衡（Sub-Wald Equilibrium）: 平衡電荷は、一貫してウォルド限界を大きく下回っている。
   • $a \lesssim 0.7$ の場合、正規化された電荷は $\xi_{eq} \approx 0.30\text{--}0.32$ である。
   • 高スピン（$a = 0.9$）の場合、$\xi_{eq} \approx 0.20$ である。
   • 極限に近いスピン（$a = 0.95$）の場合、$\xi_{eq} \approx 0.14$ である。
   • いかなるケースにおいても、電荷は $\xi = 1$ や、他の文献で提案されている非常に小さな最小エネルギー値（$\xi \lesssim 0.05$）には到達しない。
3. 解析モデルとの比較: 測定値は、磁束バランスの予測と定性的に一致しており、低スピンでの平坦なプラトーと高スピンでの減少を再現している。しかし、シミュレートされた電荷は、中間スピンにおいて解析的予測よりも系統的に20〜30%高い。著者らは、この不一致の原因として、プラズマ電流が磁束のメイスナー排出を抑制し、真空のウォルド解が予測するよりも長く極の捕捉チャネルを開いたままにしていることを挙げている。

意義と結論
本論文は、天体物理学的ブラックホールの電荷状態が、真空静電学や初期データではなく、運動論的なプラズマプロセスによって駆動されることを結論づけている。

• 上限: ウォルド電荷 $Q_W$ は、一般的な平衡値ではなく、上限として機能する。
• エネルギー抽出への影響: $\xi_{eq} < 1$ であるため、地平線と無限大の間の電位差は完全に遮蔽されない。その結果、内側ライトサーフェス付近の加速ギャップが存続し、対放電が磁気圏に負荷を与え続け、ブラックホールからのエネルギー抽出のためのブランフォード・ズナジェック（Blandford-Znajek）機構が帯電によって停止することはない。
• 限界: 著者らは、自らの結果が軸対称構成に限定されており、電流シートの非軸対称不安定性や傾いた磁場幾何学が含まれていないことを指摘している。また、シミュレーションにおける重力半径とラーモア半径のスケール分離が、実際の天体源よりも小さいことも認めている。

本研究は、平衡電荷とは、対生成を行うプラズマにおける反対符号の電荷の地平線吸収の競合によって決定される運動論的な量であることを確立しており、ブラックホール磁気圏および超高エネルギー宇宙線の加速モデルに対する新たなベンチマークを提供するものである。