[论文解读] 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

一句话总结

该论文通过表明在无粘性原行星盘中，随时间变化的径向压力支撑可驱动角动量交换，从而解决了长期存在的行星-盘相互作用悖论，即使在无粘性或波阻尼效应下，也能实现共轨区域的质量排出。作者推导出解析的自相似解，显示间隙振幅随时间线性增长，且具有普遍的径向分布，2D模拟验证了浅间隙（≤20%深度）的情况。

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）的早期间隙演化解析框架。
• 通过一系列盘和行星参数的2D数值模拟验证该模型。
• 提供一种适用于低粘性、年轻原行星盘中尘埃动力学和观测解释的灵活线性解。

提出的方法

• 构建一个包含特定角动量时间导数（˙l）作为源项的线性化动量平衡方程，以反映随时间演化的径向压力支撑。
• 利用修改后的动量方程，推导出表面密度扰动δΣ的局部和全局解析解，以角动量沉积fdep(R)作为输入。
• 使用半解析技术（包括拉普拉斯变换和格林函数）求解得到的偏微分方程，适用于全局情形。
• 从行星驱动的密度波理论出发，构建角动量沉积函数fdep(R)，假设线性波传播并在lsh处形成激波。
• 使用Athena++在无粘性、等温条件下对2D流体动力学模拟进行验证，模拟对象为低质量行星。
• 通过推导粘性应力平衡的时间尺度tl−ν，将模型扩展至粘性盘，该时间尺度表示线性解仍有效的持续时间。

实验结果

研究问题

• RQ1 为何在无粘性盘中，尽管无角动量沉积，共轨区域仍发生质量排出？
• RQ2 在无粘性或波阻尼缺失的情况下，何种物理机制可实现早期间隙形成？
• RQ3 随时间变化的径向压力支撑如何影响间隙形成过程中的角动量平衡？
• RQ4 早期间隙演化是否可由具有时间线性增长的自相似解描述？
• RQ5 在早期间隙形成阶段，解析解在粘性盘中仍有多大的有效性？

主要发现

• 由于随时间演化的径向压力支撑所导致的特定角动量时间导数（˙l）是无粘性盘中驱动共轨区域质量排出的主导机制。
• 间隙开启始于自相似过程，其径向间隙轮廓保持普遍形状，且振幅随时间线性增长，特征间隙开启时间tgap ≈ 1.5 × (Hp/Rp)⁻¹ Ω⁻¹。
• 该解析模型预测的双谷型间隙结构与2D模拟结果高度一致，适用于相对间隙深度δΣ/Σ ≤ 20%的情况，且定量吻合极佳。
• 即使在粘性盘中，线性解在有限时间尺度tl−ν ≈ 0.5 × (ν/Hp²)⁻¹内仍有效，此后粘性应力开始占主导。
• 该模型成功解释了模拟中观测到的在≈0.65Rp处出现的表面密度小凹陷，而标准稳态模型未能捕捉到此现象。
• 通过将fdep(R)扩展为fdep(R,t)，该框架可推广至时间依赖的行星参数（如迁移、吸积），同时保持解析可解性。