A Parameterized YSO Accretion Disk Model with Increasing Accretion Rate: Predicted Outburst Lightcurves

技術的サマリー： accretion rate が増加するパラメータ化された YSO 降着円盤モデル

問題提起
若い恒星天体（YSO）、特に FU オリオン型（FU Ori）亜種は、数ヶ月から数年にわたって光度が数桁増加し、その後数十年かけて減衰するという劇的な降着爆発を経験する。低状態の磁気圏降着から高状態の爆発的降着への遷移は既知の現象であるが、これらの事象における恒星光球、磁気圏衝撃波、粘性ガス円盤、受動的塵円盤といったシステム構成要素の具体的な進化を理解することは依然として困難である。既存のモデルは、しばしば複雑で計算コストの高い数値シミュレーション（MHD など）や、マルチバンドの光度曲線で観測される時間依存の形態変化を捉えられない定常状態の仮定に依存している。変化する降着率に応じたスペクトルエネルギー分布（SED）および光度曲線の時間的進化をシミュレートするために、これらの物理的構成要素を結合したパラメータ化された枠組みが必要とされている。

手法
著者らは、降着する若い星の時間系列測光とスペクトルをシミュレートするように設計された YSOpy コードに基づいたパラメータ化モデリングパイプラインを提示する。完全な放射輸送や自己整合的な数値シミュレーションとは異なり、このアプローチは 4 つの異なるシステム構成要素のスペクトルを個別に計算し、その後それらを合計するモジュール型インフラを使用する：

1. 恒星光球： 黒体または大気スラブとしてモデル化される。
2. 磁気圏降着衝撃波： 水素スラブ近似を用いてモデル化され、温度（$T_{hotspot}$）、電子密度（$n_e$）、光学的深さ（$\tau$）は、降着率（$\dot{M}$）と自由落下速度（$v_{ff}$）から自己整合的に導出される。
3. 粘性加熱ガス円盤： 内側円盤に対して修正された Shakura & Sunyaev (1973) の温度プロファイルを用いてモデル化され、1400 K 以上の温度には BT-Settl 大気グリッドを、より低温の領域には黒体をそれぞれ利用する。
4. 受動的加熱塵円盤： 照射された塵に対して Chiang & Goldreich (1997) の温度プロファイルを用いてモデル化され、内側境界は塵の昇華半径（$R_{sub}$）によって決定される。

このパイプラインは、時間依存の降着率プロファイル（$\dot{M}(t)$）を反復処理する。各時間ステップにおいて、モデルは磁気圏切断半径（$R_{trunc}$）、同期回転半径（$R_{co-rot}$）、内側塵円盤半径（$R_{dust,inner}$）を含む臨界半径を動的に計算する。内側ガス円盤半径（$R_{in}$）は磁気圧と降着圧のバランスによって決定され、低$\dot{M}$では$R_{co-rot}$から、高$\dot{M}$では$R_{trunc}$（または$R_*$）へと遷移する。モデルは、円盤が恒星表面に到達した際（$R_{in} \approx R_*$）、極性磁気圏降着が「オフ」になることを考慮する。最後に、時間分解された高分解能スペクトルを標準的なフィルター透過関数（例：Gaia、2MASS、WISE）と畳み込み、マルチバンドの光度曲線および色曲線を生成する。

主な貢献

• 統合されたパラメータ化枠組み： 本論文は、恒星、衝撃波、ガス、塵の構成要素を単一の時間進化モデルに結合する合理化されたパイプラインを導入し、3D MHD シミュレーションの計算コストなしに完全な爆発サイクルをシミュレート可能にする。
• 自己整合的な衝撃波物理学： 本モデルは、ホットスポットの温度と電子密度を降着率と自由落下速度から直接導出し、スケーリング因子$\chi$を介してエネルギー保存（$L_{slab} = L_{blackbody}$）を保証する。
• 動的な円盤境界： 本モデルは、内側ガス円盤の境界と塵の昇華前面の進化を$\dot{M}$の関数として明示的に追跡し、同期回転制限された円盤から磁気的に切断された円盤、そして最終的に恒星表面に接する円盤への遷移を捉える。
• マルチバンド診断ツール： このインフラは、可視光から赤外線バンドにわたる合成光度曲線と色曲線を生成し、観測された爆発の形態を解釈するための診断ツールを提供する。

結果
著者らは、このモデルを線形に増加する降着率と、観測された FU Ori 天体を模倣する 4 つの特定の降着プロファイル（HBC 722（3 区画線形）、V890 Aur（双曲線正接）、V960 Mon（両側指数関数）、Gaia 17bpi（ロジスティック×ガウス））に適用した。主要な知見は以下の通りである：

• 構成要素の支配的進化：

  • 低$\dot{M}$： システムは恒星光球と受動的塵円盤によって支配され、磁気圏衝撃波は無視できるか、光球よりも低温である。
  • 上昇する$\dot{M}$： $\dot{M}$が増加すると、磁気圏ホットスポットは急速に輝き、まず紫外線および青い可視光バンド（例：uバンド）を支配する。
  • ピーク/高$\dot{M}$： 内側ガス円盤の粘性加熱が、特に近赤外線（Gaia、2MASS バンド）において支配的なフラックス源となる。もし$\dot{M}$が$\sim 10^{-5} - 10^{-4} M_\odot$ yr$^{-1}$を超えると、円盤切断半径は恒星表面内側に移動し、極性降着とホットスポット放射を遮断する。
  • 中赤外挙動： 中赤外バンド（WISE）は、主に塵の昇華前面の加熱と拡大に応答し、内側塵円盤半径が外側へ移動するにつれて、可視光バンドよりも先に明るくなることが多い。

• 光度曲線の形態：

  • 赤い可視光および近赤外光度曲線は、一般的に入力された降着率プロファイルを追跡する。
  • 中赤外光度曲線は内側塵円盤の幾何学と加熱により敏感であり、可視光バンドに対して「前駆」的な明るさ増大や明確なタイミングの差を示すことがある。
  • 色の進化は変化する：システムは、ホットスポットと内側円盤の加熱により上昇期にはしばしば「明るくなると青くなる」傾向を示すが、粘性円盤が冷却され光球が再び現れる減衰期には赤くなる可能性がある。

• パラメータ感度：

  • 磁場強度： 強い磁場は円盤が恒星表面に接近するのを遅らせ、磁気圏降着相を延長し、紫外線/可視光において多峰性の光度曲線を生み出す。
  • 自転周期： 同期回転周期の変化は主に低状態の内側円盤半径に影響を与えるが、高降着爆発相における光度曲線への影響は無視できる。

意義と主張
著者らは、この作業を特定の個々の天体に対する予測モデルではなく、FU Ori 爆発の背後にある物理的メカニズムを解釈するための図示ツールとして位置づけている。彼らは、彼らのパラメータ化アプローチが、爆発サイクル全体を通じて異なるシステム構成要素が様々な波長でどのようにフラックスを支配するかを成功裏に実証したと主張する。

本論文は、波長を超えた明るさ増大のタイミングが診断的な力を持つことを強調している。紫外線/青い可視光の明るさ増大は磁気圏衝撃加熱の開始を示し、中赤外の明るさ増大は塵の昇華前面の拡大を信号する。著者らは、LSST や UVEX などの将来の施設が、これらの予測されたシグネチャを利用して初期段階の爆発を同定できることを提案している。さらに、このモデルは、観測された光度曲線と色曲線から基礎となる降着プロファイルと低状態の降着率を推定するための枠組みを提供し、単純な現象論的記述と複雑な物理シミュレーションの間のギャップを埋める。著者らは、彼らのモデルが磁場の再構成と赤道降着境界層を単純化しているものの、T タウリ型から FU Ori 型への降着状態への遷移の本質的な物理を捉えていることを認めている。