[論文レビュー] Early stages of gap opening by planets in protoplanetary discs

Amelia J. Cordwell, Roman R. Rafikov|arXiv (Cornell University)|Jul 1, 2024

Molecular Spectroscopy and Structure被引用数 1

ひとこと要約

本論文は、粘性または波動減衰効果が存在しない状況でも、共軌道領域からの質量の排出が生じるという長年のパラドックスを解消した。非粘性原始惑星系円盤における時間的に変化する径方向圧力支持が、角運動量の交換を駆動し、その結果、共軌道領域からの質量の排出が可能になることを示した。著者らは、自己相似的な解析解を導出し、ギャップの振幅が時間に比例して線形に増大し、普遍的な径方向プロファイルを示すことを明らかにした。このモデルは、浅いギャップ（≤20%の深さ）に対して2次元シミュレーションと一致する。

ABSTRACT

ABSTRACT Annular substructures in protoplanetary discs, ubiquitous in sub-mm observations, can be caused by gravitational coupling between a disc and its embedded planets. Planetary density waves inject angular momentum into the disc leading to gap opening only after travelling some distance and steepening into shocks (in the absence of linear damping); no angular momentum is deposited in the planetary co-orbital region, where the wave has not shocked yet. Despite that, simulations show mass evacuation from the co-orbital region even in inviscid discs, leading to smooth double-trough gap profiles. Here, we consider the early time-dependent stages of planetary gap opening in inviscid discs. We find that an often-overlooked contribution to the angular momentum balance caused by the time-variability of the specific angular momentum of the disc fluid (caused, in turn, by the time-variability of the radial pressure support) plays a key role in gap opening. Focusing on the regime of shallow gaps with depths of $\lesssim 20~{{\ m per\ cent}}$, we demonstrate analytically that early gap opening is a self-similar process, with the amplitude of the planet-driven perturbation growing linearly in time and the radial gap profile that can be computed semi-analytically. We show that mass indeed gets evacuated from the co-orbital region even in inviscid discs. This evolution pattern holds even in viscous discs over a limited period of time. These results are found to be in excellent agreement with 2D numerical simulations. Our simple gap evolution solutions can be used in studies of dust dynamics near planets and for interpreting protoplanetary disc observations.

研究の動機と目的

• 非粘性円盤において角運動量の注入がないにもかかわらず、流体ダイナミクスシミュレーションで共軌道領域からの質量排出が観測されるというパラドックスを解消すること。
• 初期ギャップ形成の駆動要因として、見過ごされていた特定角運動量（l）および径方向圧力支持の時間変動の役割を同定すること。
• 角運動量の時間微分（˙l）を運動量バランスに組み込む解析的枠組みを構築し、初期ギャップ進化を記述すること。
• さまざまな円盤および惑星パラメータに対して、2次元数値シミュレーションと照合してモデルを検証すること。
• 低粘性・若い原始惑星系円盤におけるダストダイナミクスおよび観測解釈に応用可能な柔軟な線形解を提供すること。

提案手法

• 時間微分された特定角運動量（˙l）をソース項として含む線形化された運動量バランス方程式を定式化し、変化する径方向圧力支持を反映する。
• 運動量方程式を修正した上で、表面密度の摂動δΣについて局所的およびグローバルな解析解を導出。角運動量注入関数fdep(R)を入力として用いる。
• ラプラス変換およびグリーン関数を用いた半解析的手法により、偏微分方程式を解く。グローバルな場合に適用。
• 惑星駆動波理論に基づき、線形波伝播およびlshにおける衝撃形成を仮定して、角運動量注入関数fdep(R)を構築する。
• 非粘性・等温条件の下で、Athena++を用いた2次元流体ダイナミクスシミュレーションと照合して予測を検証。低質量惑星を対象とする。
• 粘性円盤へのモデルの拡張を図り、粘性応力バランスに基づいて、線形解が有効である時間スケールtl−νを導出する。

実験結果

リサーチクエスチョン

• RQ1 非粘性円盤において角運動量の注入がないにもかかわらず、なぜ共軌道領域からの質量が排出されるのか？
• RQ2 粘性または波動減衰が存在しない状況で、初期ギャップ形成を可能にする物理的メカニズムは何か？
• RQ3 径方向圧力支持の時間変動が、ギャップ形成期における角運動量バランスにどのように影響するか？
• RQ4 初期ギャップ進化は、時間に比例して線形に増大する自己相似解で記述可能か？
• RQ5 初期ギャップ形成段階において、粘性円盤においても解析解はどの程度有効に保たれるか？

主な発見

• 時間的に変化する径方向圧力支持に起因する特定角運動量（˙l）の時間微分が、非粘性円盤における共軌道領域からの質量排出を駆動する主要因である。
• ギャップ形成は自己相似的プロセスとして開始され、径方向ギャッププロファイルは普遍的な形状を保ち、振幅は時間に比例して線形に増大する。特徴的なギャップ形成時間はtgap ≈ 1.5 × (Hp/Rp)⁻¹ Ω⁻¹である。
• 解析モデルは、相対的ギャップ深さδΣ/Σ ≤ 20%の範囲で、2次元シミュレーションと良好な定量的整合性を示す二つの谷を持つギャップ構造を予測する。
• 粘性円盤においても、線形解は有限の時間スケールtl−ν ≈ 0.5 × (ν/Hp²)⁻¹の間有効であり、その後粘性応力が支配的になる。
• 本モデルは、シミュレーションで観測されるが、標準的ステディステートモデルでは捉えきれない、約0.65Rpにおける表面密度の小さな谷をうまく説明する。
• fdep(R) → fdep(R,t)とすることで、時間依存の惑星パラメータ（例えば、移動、降着）にもモデルを拡張可能であり、解析的取り扱いの柔軟性を保つ。