Polycyclic aromatic hydrocarbon (PAH) abundances in the disk around T Chamaeleontis (T Cha): PAH sizes, ionization fraction, and mass during JWST observations

技術要約：T Chamaeleontis 周辺円盤における多環芳香族炭化水素（PAH）の存在量

問題と背景
多環芳香族炭化水素（PAH）は、天体物理学的環境で観測される芳香族赤外線バンド（AIB）の担体として広く受け入れられている。原始惑星系円盤（PPD）において、PAHはガスや微細な円盤表面に動的に結合しており、ガスの加熱、熱化学、および光蒸発において重要な役割を果たしている。しかし、PPDにおけるPAHの存在量は十分に制約されていない。AIBはHerbig Ae/Be型円盤では頻繁に検出されるが、T Tauri型円盤では稀であり、これはT Tauri型星のUVフラックスの低さや、クラスター形成や大きな粒子への吸着によるPAHの枯渇に起因するとされている。T Chaは、ダストギャップを持つ既知の遷移円盤を有するT Tauri型星であり、近年のジェイムズ・ウェッブ宇宙望遠鏡（JWST）の観測において顕著なAIBが検出された特筆すべき事例である。課題は、これらの特徴を生み出しているPAH集団（サイズ、電離率、および質量）を特性化し、円盤構造内におけるその分布を理解することにある。

手法
著者らは、アーカイブされたJWST Mid-Infrared Instrument (MIRI) Medium Resolution Spectroscopy (MRS) データ、ならびに光度測定およびALMA、SEST、ATCAからのミリ波フラックスを利用した。本研究では、MCFOSTコードを用いたパラメトリック3Dダスト放射伝達モデルを採用した。主な手法の手順は以下の通りである：

1. スペクトル分解： 6.2、8.1、11.3 µmにおける観測されたAIBをDrudeプロファイルを用いてフィッティングし、基礎となるスペクトル成分を特定した。この解析により、標準的なモデルでよく用いられるクラスAのプロファイルとは異なる、クラスBの寄与を含むクラスCのプロファイルが明らかになった。
2. 不透明度の修正： 標準的な「アストロPAH」不透明度（Li & Draine 2001、主にクラスA/Bプロファイルを再現するもの）では不十分であった。著者らは、T Chaで観測された特定のクラスCの特徴を再現するために、Dr 最もDrudeプロファイル解析に基づいたPAH不透明度定式化を修正した。
3. 円盤モデリング： 内側円盤（膨らんだリムとシャドウ領域を持つ）、ダストギャップ、および外側円盤を含むマルチコンポーネント円盤モデルを構築した。外側円盤は、粒径の変化を考慮して「壁」（内縁）と「外側領域」にさらに細分化された。
4. 確率的加熱： PAHは確率的に加熱されるダスト粒子として扱われた。モデルでは、PAHのサイズ分布（炭素原子数、$N_C$）、電離率（$\phi$）、およびPAH対ダスト質量比（$M_{PAH}/M_{dust}$）を変化させ、中間赤外線（MIR）の連続成分、AIBのフラックス/比、および光学的/ミリ波フラックスを同時にフィッティングした。
5. 最適化： $\chi^2$最小化法および視覚的検査を用いてモデルのグリッドを探索し、Spitzer、VLT、およびALMAによる先行研究からの制約を取り入れることで、パラメータの縮退を減少させた。

主要な貢献と結果
本研究は、T ChaのJWSTスペクトル（MIR連続成分および顕著なAIBを含む）を再現することに成功し、以下の具体的な結果を得た：

• PAH集団の特性： 最良適合モデルは、約10から26個の炭素原子（$N_C \lesssim 26$）を含む小さなPAH集団を予測している。電離率は$\phi \approx 0.15$と推定され、これはPAH集団が主に中性（~85%）であることを示している。
• PAH存在量： 導出されたPAH対ダスト質量比は $\sim 7 \times 10^{-3}$ であり、これは中間層（ISM）の値の約17%に相当する。これは、モデル化された他の円盤と比較して高い値であるとされる。
• 円盤形態と粒径分布：
  • 外側円盤構造： 14–25 µmの間の急峻な連続成分勾配をフィッティングするために、小さなダスト粒子（最大4 µm）で構成される外側円盤の壁（31–36 au）をモデルが必要とした。外側領域（36–50 au）には、より大きな粒子（最大5 mm）が含まれており、緩やかなサイズ分布（$p = -3.3$）を示している。これは効率的な粒成長を示唆している。
  • ギャップの充填： Spitzerのデータで見られる深い落ち込みと比較して、JWSTスペクトルにおける10 µm付近の連続成分のプラトーが浅いことを説明するために、モデルはダストギャップ（1–31 au）内にサブミクロンサイズのダスト粒子（0.01–1 µm）が存在することを提案している。これらの粒子は、降着ガスと結合していると考えられる。
  • 内側円盤： 内側円盤のダスト質量は、減少した5–9 µmの連続成分に一致するように、従来のモデルよりも減少させた。これはT Chaで観測された「シーソー」型の変動性と一致している。
• スペクトル分類： T ChaのAIBプロファイルは、定性的にクラスC（一部にクラスBの兆候あり）として特徴付けられ、高度にUV処理された環境と比較して、脂肪族成分と芳香族成分の比率が高い、あるいは特定の脂肪族サブグループの寄与があることを示唆している。

意義と主張
著者らは、自らのモデリングが、詳細な特性評価が極めて稀であるT Tauri型円盤におけるPAH存在量の自己整合的な推定を提供すると主張している。

• 物理的含意： 高いPAH対ダスト質量比と、円盤表面における小さくほぼ中性のPAHの存在は、T Cha系においてFUV光蒸発が非常に効率的であることを示唆している。著者らは、これが$10^{-9}$から$10^{-8} M_{\odot} \text{yr}^{-1}$のオーダーの質量喪失率を駆動している可能性があると提案している。
• 円盤の進化： ギャップ内に小さな粒子が存在することは、破壊イベント後の内側円盤の補充理論を支持しており、外側円盤における粒径の層状構造は活発な粒成長を示している。
• 限界： 論文では、モデルが特にフレアリング指数や正確なPAHサイズ分布に関して縮退を示すことを控えめに認めている。著者らは、PAH存在量の推定が、修正された不透明度が特定のクラスCの物理を正確に捉えているという仮定に基づいていること、および9–14 µmの範囲における連続成分のフィッティングが完全ではなかったことを述べている。彼らは、結果は現在のパラメータ空間内では堅牢であるが、将来的なライン放射とダスト放射の同時モデリングが、これらの縮退をさらに減少させる可能性があると結論づけている。