从一粒尘埃到一个独立天体
年轻恒星周围常有气体和尘埃组成的盘。ALMA拍到的明暗环、间隙和不对称结构，说明固体物质并非始终均匀铺开。但这些图像还没有直接展示：微小颗粒怎样变成依靠自身重力聚在一起的天体。
图像署名：ALMA (ESO/NAOJ/NRAO)；CC BY 4.0。这张2014年观测图提供研究对象的直观样貌，后文讨论的2025年气体运动样本是另一组观测。8
星子是行星形成过程中的早期自引力天体。这里关心的是尘埃与较大颗粒怎样聚成这种天体，不是已经形成的行星胚胎如何继续吞并物质。1
这条路径包含不同层次的变化：颗粒可以先靠接触粘附长大；许多颗粒可以在某片区域集中；如果集中足够强，整团颗粒才可能在自引力作用下收缩。颗粒变大、空间密度增高、形成自束缚天体，是需要分别解释的三步。
为什么不断碰撞并不能保证一路长大
最初的小颗粒通过碰撞发生接触，接触点可以粘住。随着集合体变大，结果却越来越取决于撞击速度、材料、疏松程度和大小比例。有些碰撞继续粘结，有些弹开，有些使颗粒重排压密，还有些造成侵蚀或碎裂。1
因此，增长会遇到几种不同的障碍：
| 障碍 | 发生的事情 | 需要改变什么才可能继续 |
|---|---|---|
| 反弹 | 颗粒相遇，却很少净增加质量 | 接触性质、结构或碰撞速度分布 |
| 碎裂与侵蚀 | 新增加的质量被破坏过程抵消，甚至净损失 | 撞击能量、材料承受能力或大小配比 |
| 径向漂移 | 颗粒还没来得及在该区域继续长大，就已向内迁移 | 增长与流失的时间比较、压力结构与气固反馈 |
第三种障碍不是粘不牢。即使每次碰撞都很理想，也可能赶不上物质离开这片区域的速度。反过来，把尘埃留在原地，也不保证碰撞就会变成增长。1
这正是形成机制需要相互配合的原因。能够克服一个障碍的过程，可能仍被另一个障碍卡住。
气体怎样让尘埃向内漂移
原行星盘通常向外压力降低。气体受到向外的压力支撑，绕恒星运动时可以略慢于纯粹由引力决定的圆轨道速度，也就是开普勒速度$v_K$。
固体颗粒不承受气体那样的整体压力支撑，却会与气体交换动量。这种气动阻力改变颗粒的角动量，造成径向漂移。这里不能把气体当成静止的背景：颗粒受到拖曳，同时也把力传回气体。12
描述颗粒受气体影响程度，常用斯托克斯数
$t_{\rm stop}$是阻力显著改变相对速度的时间，$\Omega$是当地轨道角频率。St很小时，颗粒与气体耦合紧；St变大以后，相对运动更明显。同一毫米大小的颗粒，在不同气体密度、位置和颗粒疏松程度下，可以有不同St。因此不能把“St=0.01”永久翻译成某个固定毫米数。
再定义两个容易混淆的含量指标：
$\epsilon$比较某处尘埃与气体的体积密度；$Z$比较把盘厚度积分后的柱密度。颗粒向中平面沉降，可以提高局部$\epsilon$，而不增加该柱中的总固体质量。即使都叫$\epsilon$，整层平均值、细丝内部的典型值和单个网格的峰值也不是同一量。
在单一颗粒大小、局部稳态、无额外黏性吸积流的气固模型中，尘埃的径向速度为
$\eta v_K$表示背景压力支撑的速度尺度，通常的向外递减压力对应$\eta>0$。于是尘埃向内、气体受到反作用而向外，局部模型的总径向质量通量相抵。这是该模型的平衡，不代表真实盘不存在向恒星的整体吸积。2
例如，固定St=0.1，$\epsilon$从0.01变到3时，模型的尘埃漂移速度从约$0.194\eta v_K$降到约$0.0125\eta v_K$。颗粒密集到足以影响气体以后，它们所面对的流动环境也会改变。
漂移变慢只是开头，怎样才能让小差异自行放大
“密集区域走得慢，所以后面的颗粒堵上来”是一个有用的直觉，但它没有完整解释流式不稳定性。
Youdin与Johansen在2007年的原始研究中专门指出：如果只把上面的稳态漂移速度代入尘埃连续性方程，小扰动表现为传播的密度波；非线性时波形可以变陡，但这一步本身并不会产生所需的线性指数增长。2
流式不稳定性，即streaming instability，涉及完整的反馈：颗粒密度扰动改变气体受到的拖曳，气体速度和压力扰动又改变颗粒运动；旋转盘中的径向与方位运动还通过科里奥利作用耦合。各变量并不是同时到达局部稳态，它们之间的相位关系决定颗粒会不会不断向密度增强处汇聚。
本库保留气体与尘埃各自的密度、三方向速度，建立八变量线性模型。气体等温、颗粒停止时间固定，不加入颗粒自引力、碰撞增长、垂直分层或外部湍流。
| 基准情况 | St | 局部密度比$\epsilon$ | 波数$(K_x,K_z)$ | 本库增长率$/\Omega$ |
|---|---|---|---|---|
| linA | 0.1 | 3 | (30,30) | 约0.4190 |
| linB | 0.1 | 0.2 | (6,6) | 约0.0155 |
波数按$\eta r$归一化。这两组增长率与2007年基准接近；保留或省略拖曳系数中的气体密度扰动，会造成约$10^{-5}\Omega$的差别，复算记录保留了具体处理与原表符号问题。2
在$\Omega t=20$时，linA的扰动增幅约为4360倍，但从百万分之一开始，最终也只是约0.44%的相对密度扰动。这个例子说明了“差异怎样被放大”，还远没有计算到浓密团块形成。
两个对照进一步限定了机制。在重新计算相应基态后，让颗粒只受气体拖曳而不反作用于气体，两组基准不再有正增长；取消背景压力驱动，结果也是如此。前者是用于比较的单向耦合近似，不再是一套完整双向动量交换模型。
阻力在这里依然耗散相对运动的能量。对整个气固系统，拖曳耗散率为
不稳定性利用的是背景压力梯度持续支持的相对运动，而不是阻力凭空制造能量。本库也核对了稳态中的压力做功、拖曳耗散和角动量输运相关的能量项相抵。2
从颗粒集中到引力塌缩，还缺一个条件判断
颗粒集中以后，自引力才有机会成为主角。但它不仅要对抗恒星造成的潮汐拉扯，还要对抗把颗粒重新分散的随机运动和扩散。
常用的密度尺度为
通常称Hill密度，也常被称为Roche密度。它提供的是自引力与潮汐之间的比较。若内部扩散或随机运动很强，仅仅某个网格超过这条线，还不足以说明一团有限大小的颗粒能持续收缩。3
团块的大小也有关系：在扩散系数相近时，扩散跨过一个小团块所需的时间，比跨过大团块短。即使密度相同，小团块也可能在重力完成聚集之前被摊开；实际收缩速度又会受到气动阻力影响。因此需要同时比较密度、尺寸与聚集和扩散的时间，而不是只看密度峰值。
Lim等人在2024年使用三维、垂直分层的局部盘模拟，同时加入外部强迫湍流与颗粒自引力。他们先让颗粒运动形成相应结构，再开启颗粒自引力，结合密度随时间的快速上升和紧凑团块的空间形态，把计算分类为出现或未出现塌缩。气体自引力没有加入，颗粒大小在每次运行中固定。3
本库从其原始表格中选出St=0.01时，每个湍流水平下相邻的未塌缩与塌缩计算点：
| 湍流扩散参数$\alpha_D$ | 已试验且未塌缩的$Z$ | 相邻已试验且塌缩的$Z$ |
|---|---|---|
| $10^{-4}$ | 4% | 6.5% |
| $10^{-3.5}$ | 10% | 12.5% |
| $10^{-3}$ | 20% | 25% |
较强湍流会把颗粒层搅厚、把局部集中扩散开，在这些模拟条件下，出现塌缩需要更高的固体含量。原文摘要用约6%、10%、20%概括拟合的临界水平；表中实际相邻计算点则如上所列。这里优先保留具体运行结果，而不把约数误写成精确门槛。
这个研究比只观察密度越线更进一步，因为它确实演化了颗粒自引力。但它依旧是具有固定颗粒大小、有限计算域、分辨率与湍流驱动方式的数值实验，不能把该表直接当作每一处原行星盘的标准配方。
颗粒增长与流式不稳定，能否互相创造有利条件
如果只研究固定大小的颗粒，会漏掉一个反馈：集中提高碰撞机会，运动环境变化也可能降低碰撞破坏；颗粒继续长大以后，又会改变与气体的耦合方式。
Carrera等人2025年的研究把这条想法具体化。他们使用2024年三维模拟的数据，拟合细丝内部典型的尘埃密度比，再把该拟合与颗粒增长、碎裂壁垒及漂移壁垒连接起来。这是一套半解析反馈模型，没有在同一项新模拟中同时求解全部碰撞、气固动力学和引力塌缩。4
先看其中的碎裂限制。设没有质量加载影响时的尺寸上限为
这里$v_{\rm frag}$是模型采用的碎裂碰撞速度，$c_s$是气体声速，$\alpha$用于描述参考湍动速度。研究者假设有限的湍动能在气体与尘埃之间分配，于是局部尘埃越多，气体与颗粒的随机运动可以减弱，碎裂上限提高。其公式为
当$\epsilon=0$时，它回到原来的$\mathrm{St}_x$。然后通过经验拟合，把更大的颗粒对应到更强的细丝富集；再用更高的局部密度更新碰撞与尺寸上限，反馈便闭合起来。4
例如，固定模型$\alpha=10^{-4}$，裸碎裂上限为St=0.01时：$Z=2\%$的轨迹停在St约0.0144、局部$\epsilon$约0.45；$Z=2.65\%$时，轨迹则持续增强，达到本库保留的St=0.1校准边界。对后一条轨迹，如果固定碰撞壁垒，St仍为0.01；如果不增加流式不稳定造成的富集，St只到约0.013。这个例子说明，两个环节相互反馈，可以产生单独改变某个环节时没有的增幅。
但要特别注意三点。
第一，固定总湍动能怎样在两相间分配，是这个模型的假设。实际湍流的能量注入、耗散和颗粒自身激发的运动可能同时变化，不能由“颗粒多了”自动推出碰撞一定更温和。
第二，模型的初始$\epsilon$来自颗粒层厚度估计，后续$\epsilon$是细丝内部的典型值；而2024年形成判据使用中平面平均密度。原论文用两倍判据作余量比较，这并不构成三种量之间的严格换算。本库因此不在图上画一条声称“越过就形成星子”的线。
第三，经验富集拟合来自限定的St、$Z$和湍流范围。图中的St=0.1是本库停止外推的位置，不是物理上禁止颗粒更大。原文的漂移受限分支也没有显示同样明显的密度增强，压力俘获等过程仍可能是必要补充。4
不同文章中的$\alpha$还可能分别表示速度方差、扩散系数或黏性应力。2025模型的$\alpha$与2024表中的$\alpha_D$不能只因符号相近就逐个等同；数值对应需要相关时间、方向和校准方式。
新的凝聚模拟说明，协同效果也会改变阶段
2026年的另一项研究，用无垂直分层、轴对称的局部模型，把代表性颗粒质量的演化加入气固方程。它跟踪单一代表尺寸，没有完整演化所有大小颗粒的分布；为构造可比较的平衡基态，质量演化方程中还加入了抵消初始增长的损失项。5
在一些中等固体含量的计算中，较密区域先凝聚得更快，颗粒大小和漂移进一步产生不均匀；随后流式不稳定在这些结构内发展得更快。这支持了“颗粒增长改变流动不稳定性”的机制联系。
不过，在已经很富尘、流式不稳定本来就很快的条件下，额外凝聚对增长和最终饱和密度的影响很小。颗粒长到较大St后，最后的状态也可能接近固定大小的流式不稳定计算。该研究还发现，协同过程能够激发更强随机运动，这又可能限制后续碰撞增长。5
所以，协同作用不能简单等同于“从头到尾都更容易形成”。某个过程可能主要帮助进入下一个阶段，而到达该阶段后，另一套动力学开始主导。
这些计算没有垂直重力，也没有直接形成自引力星子，且没有采用2025模型的全部质量加载降速假设。它们提供的是互补的条件检验，而不是对上一套半解析轨迹的逐项独立验证。
压力陷阱怎样留住颗粒，观测支持到哪一步
如果盘中出现局部压力峰，峰两侧的压力梯度方向相反，原本向内漂移的颗粒可能被收集到峰附近。在适当的近圆轨道、局部平衡条件下，这就是压力俘获。湍流扩散则倾向于把集中区域重新摊宽。1
单看一个亮环，还不能确定它一定由压力峰造成：辐射强度同时依赖温度、颗粒光学性质和光学厚度。更有区分力的做法，是同时读取气体的运动。
2025年exoALMA研究从$^{12}\mathrm{CO}$和$^{13}\mathrm{CO}$谱线的速度信息测量旋转曲线，寻找与压力变化相符的偏离。在已经能分辨和评价的子结构中，$^{12}\mathrm{CO}$的21个尘埃环有16个符合相应旋转梯度特征，12个间隙有10个符合；$^{13}\mathrm{CO}$分别为14/17和8/10。7
这里的分母不是全部天体，更不是整个宇宙中的所有尘埃环。无法充分分辨的结构未纳入该对应关系的统计；两种CO也观察同一批盘，不能当作完全独立的两份样本。
这项支持较强，因为压力特征来自谱线的速度质心，而不仅是某个位置亮不亮。但它依赖盘倾角、发射高度、恒星质量和运动分量的提取。CO常来自盘表层，毫米尘埃更接近中平面，把表层速度推到中平面压力，还需要温度和垂直结构模型。7
因此，观测支持“许多可分辨尘埃环与气体压力结构相联系，俘获机制有现实依据”，尚不能据此确认环内已经发生了哪一种小尺度不稳定，更不能数出形成了多少星子。
压力峰自身的来源也是下一层问题。它可以受到磁流体或流体结构、盘的演化和已有天体影响。若某个模型靠已经存在的行星制造压力峰，它能解释后续一批固体的聚集，却还没有解释第一批行星种子从哪里来。
三维结构与气体消散，还会重新安排形成机会
压力峰也不是永远固定的容器。尘埃集中后对气体的反作用，可以改造压力曲线，产生新的小环，改变峰的位置或让原本集中区域变宽。
2026年Tatarelli等人的研究比较了二维与三维盘中的尘埃陷阱。在所试的有黏性条件下，二维中的含尘罗斯贝波不稳定性能够造成明显结构，到了三维、尘埃沉降到中平面以后，上方仍较少受尘埃影响的气层会抑制这种发展。无黏性情况下则可以形成多个细小子环；增加尘埃输入有时主要拓宽环的分布，而不是让同一个团块密度无限增高。6
这一研究也尝试让气体逐渐减少：在尘埃得以保留的设定中，气体消散改变固体相对含量、扩散和沉降，使局部尘埃进一步集中。这提供了盘演化较晚阶段仍能打开聚集机会的一条思路。
但其压力隆起由外加力矩维持，气体移除采用规定的指数衰减，模型不含自引力；它说明的是这些设定下密度怎样变化。文中把后续可能达到Hill密度与塌缩联系起来，仍属于推断。部分后期沉降还受到垂直分辨率限制。6
这也说明，减少气体不能仅被理解为“把同一团尘埃除以一个更小分母”。如果要帮助形成星子，还要看尘埃是否被保留下来、实际空间密度是否增加、随机运动是否降低，以及团块能否在有限时间内自束缚。
目前最清楚的机制图景
从尘埃到星子，可以有一条多环节协同路径：碰撞先产生更适合气动集中与沉降的颗粒；压力结构延长它们的停留时间；气固共同运动放大局部差异；足够浓密而不易被扩散拆散的团块，才可能进一步在自引力下收缩。
这些环节各有依据：气固线性不稳定有可核的动力学；部分三维模拟已经演化出自引力塌缩；盘的气体旋转和尘埃环为压力俘获提供观测约束。尚未完全闭合的是：实际颗粒的材料与尺寸分布、随位置变化的湍流和盘结构，能否在真实演化历史中同时满足这些条件。
理解这道题时，可以追踪每个过程改变了什么：颗粒大小、气固相对运动、停留时间、局部密度，还是团块抵抗潮汐与扩散的能力。不同环境可能在不同环节受限，也可能需要不同的形成路线。
复算与回查
本库代码为research/calculations/planetesimal_formation.py。完整输入与结果、气固基准比较、稳态漂移表和原模拟相邻点保留了参数与检查输出。
气固算子采用固定停止时间，直接核基态受力、动量交换和能量平衡，并用矩阵指数检查线性演化。颗粒增长部分仅重建2025作者稿的碎裂受限分支，使用规定的富集时间尺度；减半更新步长不改变所选轨迹的终止类型，校准边界处密度变化约在0.32%以内。没有独立重跑大规模三维模拟、完整碰撞分布或引力塌缩，图中的模型轨迹也没有承担这些结论。