この論文「Bridging Scales in Black Hole Accretion and Feedback: Subgrid Prescription from First Principles（ブラックホール降着とフィードバックにおけるスケール間の架け橋：第一原理に基づくサブグリッド処方）」の技術的な要約を以下に記します。

1. 研究の背景と課題

超大質量ブラックホール（BH）とその宿主銀河の共進化を理解する上で、イベントホライズン（事象の地平面）から銀河スケール（ボンディ半径）にわたる広範な空間・時間スケールを跨ぐ「双方向の流入・流出」を記述することは中心的な課題です。

• 既存の課題:
  • GRMHD シミュレーション: 事象の地平面付近の降着やジェット噴射を詳細に解像できるが、計算コストの制約から通常は数百年のシュワルツシルト半径程度までしか進化させられず、銀河スケールでの動的平衡に達していない。
  • 宇宙論シミュレーション・半解析モデル: 銀河スケールを扱えるが、未解像の BH 降着とフィードバックを記述するための「サブグリッド処方（モデル）」が必要であり、現在は理想化された解析式や小領域 GRMHD 結果の単純な外挿に依存している。
  • 問題点: 小領域の GRMHD 結果を銀河スケールに直接適用することの信頼性が検証されておらず、特に初期条件（ガス回転や磁場強度）に依存した結果が長期的な平衡状態を正しく反映していない恐れがある。

2. 手法：マルチゾーン法を用いた第一原理シミュレーション

本研究では、グループが開発した「マルチゾーン法（multizone method）」を用いて、事象の地平面からボンディ半径（$R_B$）を超えた領域までを接続する長期 GRMHD（一般相対性磁気流体力学）シミュレーションを実施しました。

• シミュレーション設定:
  • コード: KHARMA（適応メッシュ細分化付き一般相対性 MHD コード）。
  • パラメータ空間:
    • ブラックホールスピン $|a_*| = 0, 0.1, 0.3, 0.5, 0.7, 0.9$（最大 0.97 の追加実験あり）。
    • ボンディ半径 $R_B/r_g = 4\times10^2 \sim 2\times10^6$（$r_g$ は重力半径）。
  • 初期条件: 球対称ボンディ流を基本とし、ガス角運動量や磁場強度（プラズマ$\beta$）を変化させた実験も実施。
  • 時間スケール: 各シミュレーションは $700 t_B$（ボンディ時間）まで進化させ、動的平衡状態に達するまで追跡。
• マルチゾーン法の利点:
  • 異なる時間ステップを順次反復させることで、ホライズン付近の極めて短い時間スケールと銀河スケールの長い時間スケールの両方を効率的に計算可能。
  • 小スケールと大スケールの間の十分な情報伝達を可能にし、システムが全体的な動的平衡に達することを保証する。

3. 主要な成果とサブグリッド処方

シミュレーション結果から、低エディントン比（$f_{\rm Edd} \lesssim 10^{-3}$、高温・非放射効率降着流：ADAF）における BH 降着とフィードバックのための、スピン依存のサブグリッド処方を導出しました。

A. 時間平均された処方（Deterministic Prescription）

降着率 $\dot{M}$ とフィードバック電力 $\dot{E}_{\rm fb}$ に対する解析的な近似式を導出しました。これらはボンディ降着率 $\dot{M}_B$ に対してスピン $a_*$ とボンディ半径 $R_B$ の関数として表されます。

• 降着率の抑制:
  $$\frac{\dot{M}}{\dot{M}_B} \approx (-3|a_*| + 4.7) \left(\frac{R_B}{r_g}\right)^{-0.5}$$
  スピンが高いほど降着率は抑制され、$R_B$ が増大するほど降着率は減少します。
• フィードバック電力:
  $$\frac{\dot{E}_{\rm fb}}{\dot{M}_B c^2} \approx (0.98|a_*| + 0.13) \left(\frac{R_B}{r_g}\right)^{-0.6}$$
  フィードバック電力はスピンにほぼ比例し、$R_B$ に対してべき乗則で減少します。
• フィードバック効率 ($\eta$):
  • 高速回転 BH（$a_* \gtrsim 0.7$）では、$\eta \approx 0.3$ となり、$R_B$ に依存しません。
  • 低速回転 BH（$a_* \lesssim 0.5$）では、$R_B$ に対して弱い依存性（$\eta \propto R_B^{-0.2}$）を示します。

B. 確率的処方（Stochastic Prescription）

降着とフィードバックは本質的に時間変動が激しいため、その確率分布も提供しました。

• 分布特性: 対数降着率 $\ln(\dot{M}/\dot{M}_B)$ と対数効率 $\ln \eta$ は、いずれも対数正規分布に従います。
• 変動幅: 分布の幅（標準偏差）は、ボンディ半径 $R_B$ が大きくなるにつれて増加します。
• 意義: 宇宙論シミュレーションにおいて、時間平均値だけでなく、バースト的な挙動を再現する際にこの分布情報を利用できます。

C. ブラックホールスピンの進化

低エディントン比環境における BH スピンの進化を評価しました。

• スピンダウン: すべてのケースで BH はスピンの減衰（スピンダウン）を経験します。スピンアップパラメータ $s(a_*)$ は以下の線形関係で近似されます。
  $$s(a_*) \simeq -3.7 a_*$$
• スピン凍結: スピンダウンの時間スケール $t_s$ は非常に長く、
  $$t_s \simeq 12 \left(\frac{10^{-3}}{f_{\rm Edd}}\right) \text{Gyr}$$
  となります。$f_{\rm Edd} \lesssim 10^{-3}$ の静穏な降着期では、宇宙論的な時間スケールにおいて BH スピンは実質的に「凍結」されていると結論付けられます。スピンの変化は、主に高降着率期や BH 合体時に起こると考えられます。

4. 既存研究との比較と重要な知見

• 初期条件への独立性: 従来の小領域シミュレーションでは、初期のガス回転方向（順行/逆行）や磁場強度によって結果が大きく異なることがありました。しかし、マルチゾーン法による長期シミュレーションでは、磁場が蓄積し、降着流が非軸対称化して回転方向の記憶を失うため、最終的な平衡状態は初期条件に依存しないことが確認されました。
• 小領域シミュレーションの限界: 従来の小領域 GRMHD 結果（数百年 $r_g$ 程度）をそのままサブグリッドモデルに適用すると、長期的な平衡状態（強い変動性やジェット間欠性）を捉えきれず、降着率やフィードバック効率を過大または過小評価する恐れがあることが示されました。

5. 意義と応用

• 理論的基盤の確立: 観測データに合わせたチューニングを行わず、第一原理（GRMHD）から導かれた、スピン依存のサブグリッド処方を初めて提供しました。
• 次世代シミュレーションへの適用: この処方は、コズモロジーシミュレーションや半解析モデル（SAM）に直接組み込むことが可能であり、特に BH スピン進化を追跡するモデル（例：NEWHORIZON, AREPO などの最新モデル）にとって重要な基盤となります。
• 物理的洞察: 低エディントン比の降着において、BH スピンが長期的に固定されること、およびフィードバック効率がスピンに強く依存することなど、銀河進化モデルにおける重要な物理的制約を提供しました。

この研究は、ブラックホールと銀河の共進化を記述するモデルにおいて、小スケールの物理を大スケールに正しく橋渡しするための、堅牢で物理的に裏付けられた枠組みを提供するものです。