HL 星の原始惑星系円盤外縁部における非熱的速度分散に関する論文の技術的サマリー

1. 研究の背景と問題提起

原始惑星系円盤内の乱流は、円盤構造の進化や惑星形成プロセスに決定的な役割を果たす。しかし、ガス降着-envelope（包囲ガス）に囲まれた若い埋没円盤（Class 0/I 段階）における乱流の強度は依然として不明確である。

• 既存の知見: Class II 段階（より進化が進んだ円盤）では、分子線幅の測定やダストの垂直スケール高さからの推定により、乱流パラメータ $\alpha$ が $10^{-4} \sim 10^{-3}$ と弱いことが示されている。
• 未解決の課題: より若い円盤、特に降着流（infall）や包囲ガスと関連している円盤における乱流の直接測定は困難であり、その実態が解明されていない。
• 対象天体: HL 星（HL Tau）は、タウラス分子雲に位置する Class I/II 段階の原始星であり、複数のリングとギャップ構造を持つ塵円盤を持ちつつも、まだ包囲ガスからの降着流と関連している。これまでに円盤内縁部（$r \lesssim 100$ au）では非常に弱い乱流（$\alpha \sim 3 \times 10^{-4}$）が報告されているが、外縁部（$r \gtrsim 80$ au）の乱流特性は未解明であった。

2. 研究方法と手法

本研究では、ALMA のアーカイブデータを用い、HL 星の円盤外縁部における非熱的速度分散を初めて直接測定した。

2.1. 観測データ

• 分子線: 甲醛（H$_2$CO）の $3_{1,2}-2_{1,1}$ 遷移（225.7 GHz）を使用。
• データソース: ALMA アーカイブ（2018 年観測、プロジェクト番号 2018.1.01037.S）。
• 選択理由: H$_2$CO は円盤面（midplane）付近をトレースし、降着包囲ガスやアウトフローからの汚染を受けにくいことが既知である。これに対し、CO や HCO$^+$ は包囲ガスやアウトフローを強くトレースする。
• 処理: 自己較正（self-calibration）を適用し、Briggs 重み付け（robust=0.5）を用いて画像化。空間分解能は約 0.3 角秒（約 44 au）。

2.2. 解析手法：パラメトリックモデルフィッティング

観測されたスペクトル線幅は、熱的広がりだけでなく、有限の空間分解能によるケプラー運動のせん断（Keplerian shear）の影響を強く受ける。これを分離するために、以下の手法を採用了：

1. パラメトリック円盤モデルの構築:
   • 粘性円盤モデル（Shakura-Sunyaev）を基礎とし、表面密度、温度、幾何学的パラメータ（傾斜角、位置角、中心星質量など）を定義。
   • 局所的な線幅（$\Delta V$）を半径のべき乗則関数としてパラメータ化（$\Delta V(R) = \Delta V_0 (R/100\text{ au})^l$）。
2. 放射伝達計算の近似:
   • 円盤の厚みによる速度せん断を考慮するため、「チャンネルベースの 3 層近似（channel-based three-layer approximation）」を導入。
   • 各速度チャンネルにおいて、塵層の前後にあるガス層を独立して扱い、放射伝達方程式を解析的に解くことで計算コストを削減しつつ、円盤厚みの影響を正確に評価。
3. MCMC によるフィッティング:
   • 観測された H$_2$CO の速度チャンネルマップをモデルにフィッティングし、パラメータの事後分布をサンプリング。
   • 吸収や降着ストリーマーによる汚染が疑われる特定の速度チャンネルは除外。
4. 非熱的速度分散の導出:
   • 得られた局所線幅 $\Delta V$ から、既知の温度分布（ダスト連続スペクトルや HCO$^+$ 線に基づく）を用いて推定した熱的速度成分（$v_{\text{th}}$）を差し引く。
   • 式：$\Delta V = \sqrt{v_{\text{th}}^2 + v_{\text{nth}}^2}$ より $v_{\text{nth}}$ を算出。

3. 主要な結果

3.1. 非熱的速度分散の測定値

• 半径依存性: 半径 80–180 au の範囲で、非熱的速度分散は半径とともにわずかに増加する傾向を示す。
• 平均値: 半径 80–180 au における平均非熱的速度分散は $\sim 0.15 \text{ km s}^{-1}$ である。
• 音速との比較: 仮定した温度（20 K）における等温音速（$c_s$）の約 0.6 倍（$v_{\text{nth}} \sim 0.6 c_s$）に相当する。
• 乱流パラメータ: この値は、乱流が完全に乱流起因であり $\alpha \sim \mathcal{M}^2$（$\mathcal{M}$ はマッハ数）と仮定した場合、$\alpha \sim 0.16$、または乱流マッハ数 $\mathcal{M} \sim 0.4$ に相当する。これは内縁部で報告された値（$\alpha \sim 10^{-4}$）に比べて極めて大きい。

3.2. 放射面の高度

• 速度チャンネルマップの形状（馬蹄型構造）とモデルの比較から、H$_2$CO 放射は円盤面から $z/R \lesssim 0.1$ の高度に主に由来することが確認された。これは、H$_2$CO が円盤の内部領域をトレースしていることを示唆し、外縁部の降着流や風からの影響が小さいことを裏付ける。

3.3. 誤差評価

• 統計的誤差と温度分布の不確実性を考慮しても、測定された非熱的速度（$\sim 0.15 \text{ km s}^{-1}$）は有意である（誤差範囲 $\sim 30 \text{ m s}^{-1}$）。
• モックデータを用いた検証により、この解析手法が乱流速度を $\sim 20 \text{ m s}^{-1}$ の精度で回復できることが確認された。

4. 物理的解釈と議論

4.1. 乱流駆動メカニズム

観測された大きな非熱的運動（$\alpha \sim 0.16$）を説明する可能性のあるメカニズムとして以下の二つが最も自然であると考えられる：

1. 重力不安定（GI）:
   • HL 星の外縁部円盤は質量が大きく、トーマー Q パラメータが 1.7 以下と推定されており、重力不安定を起こしやすい状態にある。
   • 数値シミュレーションでは、円盤面付近で GI による乱流が $v_{\text{turb}} \sim 0.2-0.6 c_s$ を引き起こすことが示されており、観測値と一致する。
   • 円盤質量の増加（降着）と GI による角運動量輸送のバランスが、この高い乱流を維持している可能性がある。
2. 包囲ガスからの降着（Infall）:
   • HL 星は包囲ガスからの降着流（ストリーマー）と強く関連している。
   • 円盤と降着ガスの衝突による衝撃波が螺旋密度波を生成し、レイノルズ応力を通じて円盤全体に乱流を駆動するシナリオが支持される。
   • 降着率と円盤質量の比（$\dot{M}_{\text{inf}}/M_{\text{disk}}$）は、シミュレーションで $\alpha \sim 0.1$ の乱流を生成する条件と整合的である。

4.2. 他の可能性の排除

• MRI（磁気回転不安定）: 円盤面付近では電離度が低く、MRI がこれほど強い乱流を駆動するのは困難。
• 円盤風: H$_2$CO は円盤面付近をトレースしており、風基底（$z/H \gtrsim 5$）の影響は受けない。
• 圧力勾配による非ケプラー運動: ギャップ付近での圧力勾配による速度変化は数%程度（$\sim 0.06-0.1 \text{ km s}^{-1}$）であり、観測された広がり全体を説明するには不十分。

5. 結論と意義

• 主要な結論: HL 星の円盤外縁部（80–180 au）では、内縁部とは対照的に非常に強い乱流（$\alpha \sim 0.16$）が存在する。これは、重力不安定または包囲ガスからの降着流が主要な駆動源である可能性が高い。
• 技術的貢献: 埋没円盤において、ケプラーせん断を精密にモデル化し、H$_2$CO 線を用いて非熱的速度分散を直接測定する手法を確立した。
• 科学的意義:
  • 円盤内の乱流強度が半径方向に大きく変化すること（内縁部は静穏、外縁部は乱流）を初めて実証した。
  • 惑星形成の初期段階において、重力不安定や降着流が円盤ダイナミクスに与える影響が重要であることを示唆。
  • 円盤の進化モデルや惑星形成シナリオ（特にストリーミング不安定性による微惑星形成など）への制約を提供する。

本研究は、若い原始惑星系円盤における乱流の性質と駆動メカニズムに関する理解を深める重要なステップであり、将来の高分解能分子線観測による螺旋構造の直接検出への道を開くものである。