[论文解读] Supplementary Information for ``Rapid planetesimal formation in turbulent circumstellar discs''

Anders Johansen, Jeffrey S. Oishi+4 more|ArXiv.org|Aug 29, 2007

Astrophysics and Star Formation Studies参考文献 61被引用 3

一句话总结

本篇补充材料详细描述了在湍流、磁化原恒星盘中通过尘埃石块的引力坍缩实现快速行星esimal形成的数值方法与验证测试。采用高分辨率粒子-网格代码并辅以超扩散稳定化及自洽引力场，结果表明：即使在中等质量的盘中，流体不稳定性引发的团块化也会导致引力坍缩；随着分辨率提高，坍缩发生的面密度阈值降低，且碰撞速度低于理论预测——这表明存在一种稳健的、由自引力驱动的行星esimal形成路径，其不依赖于高效的聚集过程。

ABSTRACT

This document contains refereed supplementary information for the paper ``Rapid planetesimal formation in turbulent circumstellar discs''. It contains 15 sections (\S1.1 -- \S1.15) that address a number of subjects related to the main paper. We describe in detail the Poisson solver used to find the self-potential of the solid particles, including a linear and a non-linear test problem (\S1.3). Dissipative collisions remove energy from the motion of the particles by collisional cooling (\S1.4), an effect that allows gravitational collapse to occur in somewhat less massive discs (\S1.7). A resolution study of the gravitational collapse of the boulders is presented in \S1.6. We find that gravitational collapse can occur in progressively less massive discs as the grid resolution is increased, likely due to the decreased smoothing of the particle-mesh self-gravity solver with increasing resolution. In \S1.10 we show that it is in good agreement with the Goldreich & Ward (1973) stability analysis to form several-hundred-km-sized bodies, when the analysis is applied to 5 AU and to regions of increased boulder column density. \S11 is devoted to the measurement of random speeds and collision speeds between boulders. We find good agreement between our measurements and analytical theory for the random speeds, but the measured collision speeds are 3 times lower than expected from analytical theory. Higher resolution studies, and an improved analytical theory of collision speeds that takes into account epicyclic motion, will be needed to determine whether collision speeds have converged. In \S1.12 we present models with no magnetic fields. The boulder layer still exhibits strong clumping, due to the streaming instability, if the global solids-to-gas ratio is increased by a factor 3. Gravitational collapse occurs as readily as in magnetised discs.

研究动机与目标

• 验证模拟湍流、磁化盘中行星esimal形成所采用的数值方法。
• 测试在湍流与数值耗散存在下引力坍缩的鲁棒性。
• 评估碰撞冷却与粒子团簇化在降低固态面密度阈值下促成坍缩的作用。
• 将实测碰撞速度与理论预测进行比较，并评估分辨率收敛性。
• 通过在初始模型中排除凝聚与破碎过程，隔离自引力的影响。

提出的方法

• 采用六阶有限差分与三阶龙格-库塔时间积分的Pencil Code，以确保数值耗散最小化。
• 应用六阶超扩散（nabla^6）以稳定湍流并在小尺度耗散能量，同时不影响大尺度动力学。
• 采用粒子-网格方法结合泊松求解器计算自引力，并在线性与非线性测试问题上完成验证。
• 通过动量守恒算法计算气体与粒子间的曳力，包括插值、受力计算与反作用力分配。
• 通过统计分析追踪由流体不稳定性引发的石块密度增强区域，并测量随机速度与碰撞速度。
• 在64^3、128^3与256^3三种网格分辨率下开展分辨率研究，以评估引力坍缩与碰撞动力学的收敛性。

实验结果

研究问题

• RQ1 随着网格分辨率提高，石块引力坍缩是否在更低的固态面密度下发生，表明自引力求解器具有收敛性？
• RQ2 实测的石块间碰撞速度与基于湍流速度弥散的理论预测相比如何？
• RQ3 碰撞冷却在多大程度上使质量较小的盘中发生引力坍缩成为可能？
• RQ4 更强的径向压力支持如何影响引力坍缩与粒子团簇化的启动？
• RQ5 自引力是否足以独立驱动行星esimal形成，而无需依赖高效的凝聚或破碎过程？

主要发现

• 随着分辨率提高，引力坍缩发生的面密度阈值逐步降低，表明自引力求解器具有收敛性，且数值平滑效应减小。
• 实测碰撞速度约为理论预测值的三分之一，提示需改进模型以包含轨道摆动运动（epicyclic motion）。
• 碰撞冷却使引力坍缩在固态-气体质量比低于以往认为的阈值下成为可能，显著降低了坍缩所需的面密度。
• 即使无磁场存在，当固态-气体比提高三倍时，仍能产生强烈团簇化与坍缩，证实了流体不稳定性机制的鲁棒性。
• 在高密度区域（如坍缩团簇）中，凝聚时间尺度短于一个轨道周期，表明在高密度区可实现快速生长。
• 模型表明，行星esimal形成可通过自引力驱动的团簇化，随后快速凝聚实现，而无需在初始低密度盘阶段依赖高效凝聚。