[论文解读] Planetesimal Formation without Thresholds. I: Dissipative Gravitational Instabilities and Particle Stirring by Turbulence

Andrew N. Youdin|arXiv (Cornell University)|Aug 30, 2005

Astro and Planetary Science参考文献 1被引用 4

一句话总结

本文表明，由于气体阻力引起的耗散性引力不稳定性，行星esimal的形成无需依赖传统的临界条件（如Toomre的Q < 1或Roche密度）。研究显示，通过气体阻力耗散角动量，可在比以往认为更低的密度和更高的速度弥散下实现引力坍缩，且当粒子驻留时间与轨道周期相当时，增长速度最快。

ABSTRACT

We analyze the gravitational collapse of solids subject to gas drag in a protoplanetary disk. We also study the stirring of solids by turbulent fluctuations to determine the velocity dispersion and thickness of the midplane particle layer. The usual thresholds for determining gravitational instability in disks, Toomre's criterion and/or the Roche density, do not apply. Dissipation of angular momentum allows instability at longer wavelengths, lower densities, and higher velocity dispersions than without drag. Small solids will slowly leak into axisymmetric rings since initial collapse occurs over many orbits. Growth is fastest when particle stopping times are comparable to orbital times. Our analysis of particle stirring by turbulence is consistent with previous results for tightly coupled particles, but is generalized to loose coupling where epicyclic motions contribute to random velocities. A companion paper applies these results to turbulent protoplanetary disks.

研究动机与目标

• 挑战传统观点，即原行星盘中的引力不稳定性必须满足严格的阈值条件（如Toomre的Q < 1或超过Roche密度）。
• 分析气体阻力如何通过耗散角动量，使固体在低密度和高速度弥散下仍能发生引力坍缩。
• 量化不同耦合状态下由湍流涨落激发的中平面粒子层的速度弥散和标高。
• 确立行星esimal形成比以往认为的更具鲁棒性，尤其在真实湍流和多尘盘条件下。
• 为后续关于湍流原行星盘模型中此类不稳定性可动态运作的理论研究奠定基础。

提出的方法

• 基于黏性阻力定律，推导受气体阻力作用的类流体固体盘的线性稳定性，将Toomre判据推广至包含耗散效应的情形。
• 提出引力不稳定性增长速率的改进分析，表明耗散效应使不稳定性可在更长波长和更低密度下发生，优于无耗散系统。
• 采用随机强迫方法模拟固体的湍流激发，推导出平衡速度弥散和标高作为湍流扩散率与粒子驻留时间的函数。
• 将粒子激发处理推广至弱耦合粒子情形，其中椭圆运动对随机速度的贡献显著，而不仅限于湍流涡旋。
• 推导长波和短波模式下最快增长模式的近似增长率，显示其对$ Q_T $、$ Q_R $和$ \tau_s $的依赖关系。
• 考虑气动漂移和非弹性碰撞的作用，表明在大多数情况下湍流速度主导于径向漂移，且碰撞耗散仅在极端条件下才显著。

实验结果

研究问题

• RQ1 原行星盘中固体的引力不稳定性是否可在不依赖标准阈值（如$ Q_T < 1 $）或超过Roche密度的条件下发生？
• RQ2 气体阻力如何改变固体引力坍缩的稳定性判据，特别是在角动量耗散和增长率标度方面？
• RQ3 由湍流涨落激发的中平面粒子层的平衡速度弥散和标高是多少，尤其在弱耦合区域？
• RQ4 引力不稳定性增长率如何依赖于粒子驻留时间$ \tau_s $，在何种条件下增长最快？
• RQ5 与湍流激发相比，气动漂移和非弹性碰撞在多大程度上影响粒子层的速度弥散和稳定性？

主要发现

• 耗散性引力不稳定性可在比标准判据预测更低的密度和更高的速度弥散下增长，从而无需严格依赖$ Q_T < 1 $等阈值。
• 当粒子驻留时间$ \tau_s $与轨道 timescale 相当时，增长速度最快，表明此时固气耦合最有利于不稳定性增长。
• 对于长波模式，增长率近似为$ \gamma_{\rm f} \approx 1/(\tau_s Q_T^2) $，表明中等阻力和低Toomre参数可增强不稳定性。
• 对于短波模式，增长率$ \gamma_{\rm f} \approx 2/(\tau_s (Q_R - 2)) $，在$ Q_R > 2 $时有效，表明不稳定性即使在Roche极限以下仍可维持。
• 湍流激发产生的速度弥散与紧密耦合粒子的先前结果一致，但该模型已推广至弱耦合粒子情形，其中椭圆运动对速度弥散有显著贡献。
• 在大部分参数空间中，气动漂移速度远低于湍流速度，即使在弱湍流下，由于粒子层变薄和致密化，漂移可忽略不计，对不稳定性触发无显著影响。