白矮星里的晶体为什么会浮起来

一块晶体形成时，元素分别去了哪里

先看Bauer的图三。固态内区的氧比例升高，外侧液体则更富碳。随着前沿外移，整条成分剖面都在变；这不是在温度曲线上额外贴一段延迟。

富碳液层可能因浮力而不稳定。程序混合上方液体，直到组成梯度符合采用的稳定条件；搬运前后各元素的总量相等。先读公式三至六怎样处理质量不相等的网格区，再回来看彩色剖面是怎样算出来的。

碳和氧都大致每两个核子配一个电子，单靠理想电子压，难以算出分离的能量。较小的、依赖组成的等离子体项，以及随后的结构调整，都有作用。论文调用状态方程并求解恒星结构；本站固定密度的纯碳模型既没有这套成分搬运，也没有它产生的加热。

来源：Bauer 2023, Carbon-Oxygen Phase Separation in MESA White Dwarf Models, arXiv2303.10110v2, Eq1–8 and §§2–5 · Bauer 2023, fixed MESA inputs and starting models, Zenodo 7846751, CC BY 4.0

物理相界与数值开关，要分开读

Bauer把无量纲自由能差定义为 $\Delta f=f_{i,\rm liquid}-f_{i,\rm solid}$ 。它为正时，固态的自由能较低。公式一、二采用下式平滑变量，并取 $w=10^{-2}$ 。两相自由能相等仍发生在 $\phi=0.5$ ；平滑让熵的变化能连续进入恒星演化的时间步。

程序却等到 $\phi_{\rm sep}=0.9$ 才开始改变成分。这样做，是为了避免成分跳变打断平滑区内的潜热释放。论文报告，经过这段区间时光度约变动 $0.04\,\mathrm{dex}$ ，冷却延迟最多可能因此高估约一成。这是一项需要评估后果的数值选择，没有把物理结晶条件改到零点九。

平滑自由能中的权重也不能写反：液态权重减小，固态权重增大。由这套完整热力学求导，已经含有潜热；再独立加一份潜热源会重复计算。这里的简化模型采用另一种明确拆法：连续的显热支路，加上单独的潜热储库。

\phi=\frac{\exp(\Delta f/w)}{1+\exp(\Delta f/w)},\qquad f_{i,\rm smooth}=(1-\phi)f_{i,\rm liquid}+\phi f_{i,\rm solid}

这是Bauer 2023第二版的公式一、二，包含原文自由能差的正负约定。这是状态方程的相变量；径向模型的固相质量占比则来自球壳积分。

来源：Bauer 2023, Carbon-Oxygen Phase Separation in MESA White Dwarf Models, arXiv2303.10110v2, Eq1–8 and §§2–5

分离加热要怎样进入能量方程

Bauer的公式八先固定时间步起点的温度与压力，再改变组成。密度会略微调整。改变前后的比内能差，给出这一操作的加热率，随后作为源项进入恒星结构求解。

这份重排释放的能量，叠加在结晶潜热之外。把它概括成引力能释放，有助于建立图景；但在不变的密度剖面上猜一项引力加热，不能代替这套计算。那样的近似也必须有自己一致的结构、组成与能量收支。

读图四和后面的程序比较时，先找初始碳氧剖面、氢氦层、导热不透明度和沉降处理。Bauer报告的MESA、BaSTI与LPCODE氖沉降延迟并不完全相同，原文没有给出一种通用的差值修正。保留这些模型差异，比把它们塞进任意高斯误差更有信息。

\epsilon_k=\frac{e(\rho_{k,\rm start},T_k,\{X_i\}_{k,\rm start})-e(\rho_{k,\rm end},T_k,\{X_i\}_{k,\rm end})}{\delta t}

这里的

e

是比内能，因此

\epsilon_k

的单位是每单位质量、每单位时间的能量。前后密度不必相等。这份源项进入随后的恒星演化求解。要检验全局能量，还需计入其他能量变化与损失。

来源：Bauer 2023, Carbon-Oxygen Phase Separation in MESA White Dwarf Models, arXiv2303.10110v2, Eq1–8 and §§2–5 · MIST white dwarf models, arXiv2509.21717v1, cooling clocks, surface boundary and model comparisons

氖把浮沉问题改了

再翻到Bédard、Blouin与Cheng的物理图景和Methods。新生晶体可能比周围液体更贫氖，在提出的条件下会上浮、融化，较重液体则下沉。反复发生的分离把氖向内集中，并释放能量。这就是这里说的蒸馏。这一阶段，产晶区整体仍是液态，没有像普通结晶那样不断积起固态内核。

作者修改STELUM演化计算来描述这个过程。初始氧剖面和氦层质量会影响结果。论文图四的默认一族使用薄氦层和指定氧剖面；本课程选到的三个文件则属于30PR30标准非均匀氧剖面、较厚氦层的一族。它们出自同一研究，却不能悄悄互换延迟。

还要看版本。早期STELUM网格没有包含碳氧相分离，不等于后来每一版STELUM计算都没有。应按2024论文Methods确认这次实际采用的物理。读者可以重算公开表面序列之间的差值；要重新演化恒星，还缺作者未公开的程序。

来源：Bédard, Blouin and Cheng 2024, Buoyant crystals halt the cooling of white dwarf stars, arXiv2409.04419v1, author manuscript of the Nature study, main text and Methods · Bédard, Blouin and Cheng 2024, fixed cooling sequences, Zenodo 10201676, CC BY 4.0 · Bauer 2023, Carbon-Oxygen Phase Separation in MESA White Dwarf Models, arXiv2303.10110v2, Eq1–8 and §§2–5

从起点处方，读到最后的组成

Methods的蒸馏起始拟合使用离子数目占比，声明的低氖范围是 $x_{^{22}\mathrm{Ne}}<0.04$ 。模型随后规定，中间组成随碳的库仑耦合参数线性变化，走向约 $\Gamma_{\rm C}=208$ 、 $x_{\rm C}=0.8$ 、 $x_{\rm O}=0$ 、 $x_{^{22}\mathrm{Ne}}=0.2$ 的终点。低氖起始拟合不能一路外推到这个终点。

数目占比和质量占比是两件事。终点只剩碳与氖时，用下面的原子质量换算，氖质量占比约为 $0.31$ 。固定PDF中紧跟其后的上标八是文献编号，不是第三位小数。只读抽取文字、不看印刷页面，可能把引用变成一场并不存在的数值矛盾。

从起始到中央最终组成的路径、蒸馏何时停止，都仍含模型假设。Methods还说明，出于数值原因，后续内侧约四分之一质量的结晶被人为放慢。解读晚期曲线时，应把这项处理放在旁边；仅凭表格，分不清一个小特征是否已经解析了内部物理。

最后，公开输出经过约十亿年的平滑。把相邻行插得再密，也不会获得更细的物理时间分辨率。比较表面序列时，应同时标出原生时钟和平滑尺度；要画内部结晶前沿，还需要径向数据。

X_{^{22}\mathrm{Ne}}=\frac{22x_{^{22}\mathrm{Ne}}}{12x_{\rm C}+22x_{^{22}\mathrm{Ne}}}=\frac{22(0.2)}{12(0.8)+22(0.2)}\simeq0.314

这个换算检查终点占比的含义，并没有推导出终点，也没有验证中间的蒸馏路径。

来源：Bédard, Blouin and Cheng 2024, Buoyant crystals halt the cooling of white dwarf stars, arXiv2409.04419v1, author manuscript of the Nature study, main text and Methods · Bédard, Blouin and Cheng 2024, fixed cooling sequences, Zenodo 10201676, CC BY 4.0

做一个现有数据真能回答的比较

把教程的原生比较命令用于三条固定的一倍太阳质量30PR30序列。在 $\log_{10}(L/L_\odot)=-3$ 处，选定的部分与完整蒸馏序列比标准模型晚约 $1.764\,\mathrm{Gyr}$ 和 $5.417\,\mathrm{Gyr}$ 到达。程序从真实时间、光度行计算交点，保留局部回头。这是在所选平滑模型一族内可以复算的差值。

这一步没有重现论文图四的默认模型，没有求出径向固相比例，也没有量出某颗星“缺失”的年龄。若要独立检验能量来源，还需要状态方程、相图、搬运与混合规则、结构响应，以及完整加热和损失项。几列表面量不能补出这些材料。

在更大的演化计算之前，还有一个值得问的问题：处理大气分类、双星和选样函数后，观测中的堆积还剩多少？条件Gaia练习提供原始字段比较，也保留拟合失败，却没有解决这个总体问题。冷却变慢、变慢的机制，以及有多少星停在那里，是相互关联、但各需证据的三件事。

python -W error scripts/astronomy/ac31/run_native_comparison.py

来源：Bédard, Blouin and Cheng 2024, Buoyant crystals halt the cooling of white dwarf stars, arXiv2409.04419v1, author manuscript of the Nature study, main text and Methods · Bédard, Blouin and Cheng 2024, fixed cooling sequences, Zenodo 10201676, CC BY 4.0 · Sandbox Physics, standalone white-dwarf scientific reference, original model inputs and independent tests