[论文解读] Embryos grown in the dead zone: Assembling the first protoplanetary cores in low mass self-gravitating circumstellar disks of gas and solids

Wladimir Lyra, Anders Johansen+2 more|ArXiv.org|Jul 16, 2008

Astrophysics and Star Formation Studies参考文献 20被引用 47

一句话总结

该论文提出，在低质量、自引力的原行星盘中，死区边缘的罗斯比涡旋中被捕获的固态物质的引力坍缩可迅速形成火星大小的原行星胚胎。通过包含自引力和多种粒子种类的全局模拟，研究显示固体在5个轨道内达到临界密度，坍缩形成38个引力束缚的胚胎——其中一半质量超过火星，这是由于高效的粒子捕获和阻力冷却所致，速度弥散度随时间从约10 m s⁻¹降至小于1 m s⁻¹。

ABSTRACT

In the borders of the dead zones of protoplanetary disks, the inflow of gas produces a local density maximum that triggers the Rossby wave instability. The vortices that form are efficient in trapping solids. We aim to assess the possibility of gravitational collapse of the solids within the Rossby vortices. We perform global simulations of the dynamics of gas and solids in a low mass non-magnetized self-gravitating thin protoplanetary disk with the Pencil code. We use multiple particle species of radius 1, 10, 30, and 100 cm. The dead zone is modeled as a region of low viscosity. The Rossby vortices excited in the edges of the dead zone are very efficient particle traps. Within 5 orbits after their appearance, the solids achieve critical density and undergo gravitational collapse into Mars sized objects. The velocity dispersions are of the order of 10 m/s for newly formed embryos, later lowering to less than 1 m/s by drag force cooling. After 200 orbits, 38 gravitationally bound embryos were formed inside the vortices, half of them being more massive than Mars. The embryos are composed primarily of same-sized particles. We conclude that the presence of a dead zone naturally gives rise to a population of protoplanetary cores in the mass range of 0.1-0.6 Earth masses, on very short timescales.

研究动机与目标

• 研究位于原行星盘死区边界处的罗斯比涡旋中被捕获的固体是否能发生引力坍缩，形成原行星胚胎。
• 评估自引力在低质量、自引力盘中，由厘米至米级固体形成束缚团块过程中的作用。
• 考察气体阻力和粒子尺寸分布对形成中胚胎的集中度、速度弥散度和稳定性的影响。
• 确定死区边界处的罗斯比波不稳定性是否为形成寡头式胚胎提供了一条可行路径，而无需预先存在的小行星。
• 在包含多种粒子种类和自引力的全局盘模拟框架下，评估在真实条件下胚胎形成的鲁棒性。

提出的方法

• 使用Pencil代码对薄型、自引力的原行星盘中的气体和固体进行全局模拟，其中死区被建模为低粘性区域。
• 采用混合的埃普斯坦-斯托克斯阻力模型计算气体阻力，并包含粒子阻力对气体的反作用。
• 基于多重傅里叶变换的粒子-网格泊松求解器，用于计算固体的自引力，从而实现对引力坍缩的追踪。
• 模拟四种不同粒子尺寸（1、10、30和100 cm半径）以研究尺寸依赖的动力学行为和分异。
• 实施垂直平均的流体动力学模型，采用具有陡峭过渡的粘性分布以模拟死区边界。
• 采用高分辨率模拟（最高达512²个网格单元和4×10⁵个粒子），以解析涡旋结构和胚胎形成动力学。

实验结果

研究问题

• RQ1 在自引力盘中，死区边缘的罗斯比涡旋是否能高效捕获固体并触发其引力坍缩，形成原行星胚胎？
• RQ2 在该类涡旋中，自引力在由厘米至米级固体形成束缚团块的过程中起什么作用？
• RQ3 粒子尺寸分布和气动分选如何影响形成中胚胎的速度弥散度和结构一致性？
• RQ4 气体阻力和气态涡旋的潮汐力在多大程度上影响固体团块的稳定性和生长？
• RQ5 所描述的机制能否产生质量与火星相当的胚胎群体，与寡头式形成情景一致？

主要发现

• 在涡旋形成后的5个轨道内，罗斯比涡旋中的固体达到临界密度并发生引力坍缩，形成束缚胚胎。
• 200个轨道后共形成38个引力束缚胚胎，其中一半质量超过火星（0.1–0.6 M⊕）。
• 胚胎的初始速度弥散度约为10 m s⁻¹，随时间因气体阻力冷却而降低至1 m s⁻¹以下。
• 最重的胚胎主要由单一粒子尺寸组成（例如，98%为30 cm粒子），表明存在强烈的尺寸分异。
• 粒子尺寸分布更均匀的胚胎表现出更低的速度弥散度，表明其稳定性与结构一致性更高。
• 自引力和粒子反作用的存在改变了涡旋的演化，使其在自引力模型中比非自引力模型更稳定，不易被破坏。