Modelling dust coagulation, dynamical drag and turbulent mixing during star and disc formation

技术摘要：模拟恒星与原行星盘形成过程中的尘埃凝聚、动力学阻力与湍流混合

问题陈述
尘埃在恒星和行星形成过程中起着至关重要的作用，影响着气体冷却、电离水平以及磁场与气体的耦合。虽然尘埃凝聚成更大颗粒（卵石）对于微行星形成和核心生长至关重要，但以往的三维流体动力学恒星形成模拟在很大程度上忽略了尘埃相对于气体的动力学演化。早期的研究方法（如 Bate 2022）将尘埃种群视为在每个流体单元（SPH 粒子）内独立演化的过程，并假设尘埃与气体完美耦合。这种近似在早期、低密度的坍缩阶段是成立的，但在更长的时间尺度上（特别是在密度更高、颗粒生长更快的原行星盘内），它无法解释尘埃的迁移、沉降和混合。此外，现有模型通常缺乏关于湍流如何同时驱动颗粒生长所需的相对速度以及驱动尘埃垂直混合的湍流扩散的自洽处理。

方法论
作者在平滑粒子流体动力学（SPH）代码 SPHNG 中实现了一种结合三种不同物理过程的新方法：

1. 尘埃凝聚： 利用离散化平均场 Smoluchowski 方程法（Hirashita & Yan 2009; Bate 2022），该代码追踪每个 SPH 气体粒子中 $N$ 个对数间隔尺寸间隔内的尘埃种群演化。该方法考虑了在特定凝聚阈值速度以下的双粒子碰撞，且碰撞粘附概率为 1。
2. 动力学阻力 (MULTIGRAIN)： 作者集成了 MULTIGRAIN 算法（Hutchison et al. 2018），该算法在终端速度近似下求解气体与多种尘埃组分的连续流体方程。这使得计算由于压力梯度和引力导致的尘埃不同尺寸之间与气体之间的整体差异速度成为可能。
3. 隐式积分与湍流扩散： 一个关键创新是使用最近开发的隐式积分方法（Elsender & Bate 2024）来求解尘埃阻力方程。这避免了处理具有长停止时间的颗粒时所面临的严格显式方法时间步长限制，并防止了不物理的负尘埃分数的出现。此外，作者引入了一种新的隐式方法来模拟湍流扩散。他们将尘埃的湍流搅拌视为一个扩散过程（$D_i = \alpha_{SS} c_s H / Sc$），其中施密特数（$Sc$）解释了不同尺寸颗粒的耦合效率。该扩散项使用高斯-赛德尔算子分裂法进行求解，允许尘埃在流体单元之间进行混合。

代码追踪了从分子云核坍缩、经过第一水文静态核心阶段，到原行星盘形成与演化的全过程，并利用汇粒子（sink particles）来追踪系统在恒星核心形成之后的演化。

主要贡献

• 首个三维耦合模拟： 本工作代表了首次在恒星及盘形成早期阶段同时模拟尘埃凝聚、动力学气体阻力和湍流混合的三维流体动力学模拟。
• 隐式扩散方案： 开发了一种鲁棒的湍流尘埃扩散隐式求解器，能够实现稳定的尘埃混合模拟，而不会受到与显式扩散方案相关的严重时间步长限制。
• 统一处理： 代码在一个统一的框架内整合了相对速度（驱动生长）和湍流扩散（驱动混合）的处理，并在测试用例中使用相同的湍流参数（$\alpha_{SS}$）。

结果
作者展示了来自测试计算（垂直柱中的尘埃沉降）和完整恒星形成模拟的结果：

• 测试计算： 代码成功恢复了尘埃沉降与湍流扩散之间的平衡解析解。测试确认，包含湍流扩散可以防止形成不切实际的极薄尘埃层，并保持有限的尘埃标高。
• 恒星形成模拟：
  • 第一水文静态核心： 与之前的工作（Bate 2022）一致，在第一水文静态核心内发生了快速的尘埃生长，在核心深处产生了局部单分散的种群，其尺寸超过 100 $\mu$m。由于演化时间尺度较短，加入阻力和扩散对这一早期阶段的影响可以忽略不计。
  • 原行星盘演化： 一旦盘形成，湍流扩散的引入显著改变了尘埃演化。
    • 增强生长： 湍流扩散显著增强了尘埃颗粒的生长速率。通过不断向高密度区域（即大颗粒生长的区域）供应小尺寸和中等尺寸的颗粒，扩散防止了凝聚过程中“食物”的局部耗尽。
    • 湍流依赖性： 湍流水平（$\alpha_{SS}$）至关重要。高湍流（$\alpha_{SS} = 10^{-3}$）导致整个盘内迅速生长为尺寸约为 30-100 $\mu$m 的单分散种群。较低的湍流（$\alpha_{SS} = 10^{-4}$）则导致生长较慢，并产生更具垂直分层的尘埃分布，较大的颗粒被限制在更靠近中平面处。
    • 尘埃沉降： 在弱湍流模拟中，尘埃沉降导致中平面尘埃/气体比显著增加（超过 50%），而小颗粒在靠近中平面处仍然稀少，并被限制在表面层。
  • 阻力 vs. 扩散： 对比包含或不包含特定过程的模拟发现，在这些时间尺度内，尘埃-气体动力学阻力对尘埃分布的影响很小，因为颗粒仍保持良好的耦合。相比之下，包含湍流扩散具有巨大的影响，它使尘埃分布更加平滑并加速了凝聚。

意义与主张
论文声称，这种新方法首次实现了在恒星及其盘形成过程中对尘埃生长、混合和迁移的自洽建模。作者强调，虽然加入阻力和扩散不会改变关于第一水文静态核心中尘埃生长的基本结论，但对于随后的原行星盘演化至关重要。具体而言，他们声称湍流尘埃扩散通过作为小颗粒向大颗粒形成的区域的来源，实质性地增强了尘埃颗粒的生长速率。

作者对研究范围保持谦逊，指出：

• 结果取决于所采用的湍流模型（例如 $\alpha$-盘或 MHD 湍流）。
• 目前的模拟只考虑了凝聚而未考虑破碎；不同的凝聚阈值可能会导致破碎并重新产生小颗粒。
• 辐射传输计算仍假设使用标准的星际消光曲线，而非使用演化中的局部粒径分布，这是作者确定的一个必要的未来扩展步骤。

最终，这项工作提供了一个工具，用于将三维流体动力学模拟与来自 ALMA、VLT-SPHERE 和 JWST 等设施的大量新观测数据（特别是关于边缘观测盘中的尘埃沉降和湍流搅拌）进行对比。