顺着亮环与聚团阈值读原论文

先读真实亮度剖面

打开官方 DSHARP 的 Elias24 发布项。ASCII 产品有300个径向采样、八列，原始 README 解释半径、每合成波束强度、近似逐点误差，以及两种亮温换算。当前文件包含2019年的半径列修正，发布页于2020年补齐更新。

下图采用保存的全部强度，包括外盘含噪负值。预先指定45–65天文单位窗口，最低点在57.5；65–95窗口的最高点在76.5，两者亮度比为43.3。这是所选窗口中实测剖面的比例。

径向采样间隔为一天文单位，不代表每格都是独立观测，合成波束仍跨数个天文单位。图中的误差条保留来源逐点误差；若要拟合环宽并给置信区间，还需空间协方差与标定处理。

对应阅读：ALMA DSHARP official data release and revised radial profiles · DSHARP radial-profile original README · Huang et al.2018, DSHARP II Annular Substructures, v1 · Andrews et al.2018, DSHARP I Motivation, Sample, Calibration and Overview

换亮温前，先读波束

固定 FITS 主头给出频率 $231.943014843\,\mathrm{GHz}$，高斯波束宽度 $0.0367535\,\mathrm{arcsec}$ 与 $0.0342672\,\mathrm{arcsec}$。公开包保存精确主头块，没有下载大图像或其中像素。

把角宽转为弧度，计算 $\Omega_{\rm beam}=\pi\theta_{\rm maj}\theta_{\rm min}/(4\ln2)$。Jy/beam 除以波束立体角，再乘 $10^{-26}$，得到 SI 谱辐亮度。公开脚本重算 Rayleigh–Jeans 与完整 Planck 列，最大残差分别小于 $0.000854\,\mathrm K$ 与 $0.000738\,\mathrm K$，包含来源列的舍入。

Planck 反演只对正强度定义。文件里31个非正强度点保留在强度图中，亮温在输出里记为 null。把亮度解释成物质结构，还要温度与不透明度；仅凭发射环，不能确定它的形成机制。

对应阅读：DSHARP radial-profile original README · ALMA DSHARP official data release and revised radial profiles · Huang et al.2018, DSHARP II Annular Substructures, v1

用相应密度比读聚团图

再打开 Li 与 Youdin 2021第二版。第2节建立垂直分层模拟，图1以 Stokes 数和面密度金属丰度 $Z=\Sigma_{\rm d}/\Sigma_{\rm g}$ 标记强聚团区域。它与均匀线性计算里的局部体密度比 $\epsilon=\rho_{\rm d}/\rho_{\rm g}$ 是不同量。

先看原文怎样测尘层厚度，再看面丰度怎样换成中平面的密度比。把压强支撑和搅动条件留在阈值旁边；从一个参数切片读出的边界，移到另一种盘时，条件已经变了。

图中的判据指给定模拟里的强聚团。要形成有引力束缚的天体，还需引力、潮汐和粒子塌缩计算。集中起来与长成行星之间，仍有后续步骤。

对应阅读：Li & Youdin 2021, Thresholds for Particle Clumping, v2

比较线性图与非线性模拟

同一论文第3.5节、图7，给无分层线性模式加上波长与分辨率限制。检查波能否装进测得的尘层，以及网格能否分辨它；形式上任意短波长的最快根，未必在模拟中可用。

图8将相应线性率与分层聚团结果比较，直接检验“微扰长得快，最终就聚得强”的想法。原文发现相关性很差，并指出其中跨了三步：线性到非线性、无分层到有分层、模式增长到粒子集中。

这时可以把自己的教程放在旁边。周期无分层网格跟踪小波怎样随时间变化，聚团图则研究沉降与混合层的后续行为。两者提出的问题和观察的状态各有范围。

对应阅读：Li & Youdin 2021, Thresholds for Particle Clumping, v2 · Youdin & Goodman 2005, Streaming Instabilities, arXiv v1

追问扩散系数代表什么

打开 Umurhan、Estrada 与 Cuzzi2020第四版的第3.1–3.3节和附录A，分别找到气体粘性、粒子扩散与粒子压强。背景中还可能有粘性强迫，这些作用不只是在最后增长率上减去一个数。

接着看第5节的湍流稀释模型，怎样把搅动与局部体密度比连起来。改变湍流，既可能阻尼扰动，也可能改变背景尘层。只改一个扩散系数、其余都固定，是目的较窄的控制实验。

Li 与 Youdin 第3.5.3节还表明，简单扩散阻尼线性率不能预言其非线性分层聚团。比较不同阈值时，要带上扩散的假定方式、实测或拟合数值、几何和可用模式。

对应阅读：Umurhan, Estrada & Cuzzi 2020, Streaming Instability in Turbulent Disks, v4 · Li & Youdin 2021, Thresholds for Particle Clumping, v2

写一个能检验的下一步问题

对于 Elias24 实测亮环，可以继续问：气体运动是否支持同位置存在气压最高点？对于线性计算，问增长是否被分辨、在指定混合模型下是否仍在？对于稠密模拟团块，问自引力能否压过潮汐与速度弥散？这些问题各自指向不同观测或计算。

Huang 等人的 DSHARP II 第5节讨论环状结构的几种起因。把真实亮度剖面留在这些解释旁边：合成径向陷阱帮助理解输运，线性教程检验集体机制，真实环的成因则要继续用观测区分。

对应阅读：Huang et al.2018, DSHARP II Annular Substructures, v1 · Li & Youdin 2021, Thresholds for Particle Clumping, v2 · Youdin & Goodman 2005, Streaming Instabilities, arXiv v1