[论文解读] Catching drifting pebbles II. A stochastic equation of motions for pebbles

Chris W. Ormel, Beibei Liu|arXiv (Cornell University)|Mar 16, 2018

Astrophysics and Star Formation Studies参考文献 88被引用 31

一句话总结

本文为原行星盘中 Pebble 尺寸的颗粒推导出一个随机运动方程（SEOM），结合湍流气体阻力与扩散效应，以模拟其垂直输运与吸积效率。SEOM 表明，湍流主要通过垂直扩散机制抑制 Pebble 吸积，而在高相对速度下则通过沉降机制失效实现抑制，且在外部盘区与低质量恒星周围效应更显著。

ABSTRACT

Turbulence plays a key role in the transport of pebble-sized particles. It also affects the ability of pebbles to be accreted by protoplanets, because it stirs pebbles out of the disk midplane. In addition, turbulence can suppress pebble accretion once the relative velocities become too large for the settling mechanism to be viable. Following Paper I, we aim to quantify these effects by calculating the pebble accretion efficiency $\varepsilon$ using three-body simulations. To model the effect of turbulence on the pebbles, we derive a stochastic equation of motion (SEOM) applicable to stratified disk configurations. In the strong coupling limit (ignoring particle inertia) the limiting form of this equation agrees with previous works. We conduct a parameter study and calculate $\varepsilon$ in 3D, varying pebble and gas (turbulence) properties and accounting for the planet inclination. We find that strong turbulence suppresses pebble accretion through turbulent diffusion, agreeing within factors of order unity with previous works. Another reduction of $\varepsilon$ occurs when the turbulent rms motions are large and the settling mechanism fails. Efficiency-wise, the outer disk regions are more affected by turbulence than the inner regions. At the location of the H$_2$O iceline, planets around low-mass stars achieve much higher efficiencies. Including the results from Paper I, we present a framework to obtain $\varepsilon$ under general circumstances.

研究动机与目标

• 建模 Pebble 在湍流原行星盘中的随机运动，考虑其与气体的气动耦合及湍流扩散效应。
• 量化湍流如何通过破坏中平面聚集与增加相对速度，降低 Pebble 吸积效率（ε）。
• 推导一个通用的 ε 框架，整合粒子停止时间、气体扩散率、关联时间与引力作用。
• 研究盘层结、湍流特性（如 αz）与行星倾角在塑造吸积结果中的作用。
• 通过提供适用于多种盘与行星条件的 3D、湍流解析的 ε 预测公式，拓展论文 I 的结果。

提出的方法

• 推导 Pebble 的随机运动方程（SEOM），在牛顿动力学基础上引入由湍流气体耦合产生的随机力。
• 引入两个关键湍流参数：气体扩散率（D_gas）与关联时间（t_corr），以控制粒子运动的统计特性。
• 将 SEOM 应用于具有垂直变化湍流扩散率（αz(z)）的层结盘配置，允许中平面聚集与表面扩散。
• 通过三体模拟数值求解 SEOM，计算在不同湍流强度、粒子尺寸（τs）与行星质量下的 Pebble 吸积效率（ε）。
• 在强耦合极限（t_stop → 0）下验证 SEOM，结果与假设白噪声（不相关）湍流的先前研究一致。
• 将结果与论文 I 的二维表达式整合，构建统一的、通用的 ε 预测公式，作为 Pebble、行星与盘参数的函数。

实验结果

研究问题

• RQ1 湍流扩散如何影响原行星盘中 Pebble 的垂直分布与中平面聚集？
• RQ2 当沉降机制失效时，湍流相对速度在多大程度上抑制 Pebble 吸积？
• RQ3 湍流特性（如 αz、t_corr）的变化如何影响不同盘区与行星质量下的吸积效率 ε？
• RQ4 湍流扩散与湍流速度弥散在降低 ε 中的相对重要性如何？
• RQ5 在存在湍流的情况下，吸积效率 ε 如何依赖于行星质量、倾角与恒星质量，特别是在不同区域？

主要发现

• 湍流扩散是抑制 Pebble 吸积效率（ε）的主要机制，因其将 Pebble 从中平面抬升，降低其局部密度。
• 对于 τs ≈ 0.1–1 的颗粒，理想 MRI 模拟中 ε 低于无碰撞扩散模拟，原因在于更强的湍流速度效应，其中 σ ≳ v_*，导致沉降机制边缘失效。
• SEOM 框架在 t_stop → 0 极限下恢复强耦合近似（SCA），验证了其与假设白噪声湍流的先前模型的一致性。
• 由于更强的引力聚焦与更高的 Pebble 密度，内盘中 Pebble 吸积更高效，而外盘区域则因湍流抑制更严重。
• 围绕低质量恒星的行星在 H2O 冰线处可实现更高的 ε，原因在于更大的 Pebble 捕获半径与更低的湍流速度弥散。
• 所提出的框架可实现通用的、三维的 ε 预测公式，综合考虑湍流效应（扩散与速度）与引力聚焦，适用于多种盘与行星条件。