---
title: "白矮星里的晶体为什么会浮起来"
description: "沿固定版Bauer与Bédard原文，从碳氧分离读到氖蒸馏，保留成分处方、数值平滑和待解问题。"
canonical_url: "https://sandboxphysics.com/zh/stories/white-dwarf-cooling/white-dwarf-phase-separation-distillation/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/white-dwarf-cooling/white-dwarf-phase-separation-distillation/"
last_modified: "2026-11-03"
---

[白矮星的热与时钟](https://sandboxphysics.com/zh/stories/white-dwarf-cooling/)

L4 · 白矮星冷却与结晶

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

Sandbox Physics · 2026-11-03

结晶会释放潜热，也可能重新分配整颗星里的元素。碳、氧和氖进入晶体的比例，未必等于周围液体中的比例。沿两篇固定版论文追下去：先看MESA里的碳氧分离，再看一种可能让大质量白矮星多亮数十亿年的氖蒸馏。

## 开始前需要什么

先读《结晶怎样拖慢白矮星冷却》和原生模型教程。阅读版固定为Bauer的arXiv2303.10110v2，以及Bédard、Blouin、Cheng的arXiv2409.04419v1作者稿；后者对应2024年三月发表于Nature的研究。把论文中的完整演化计算、可取得的标量表和这里的受控热学模型分开。

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

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

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

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

来源：[Bauer 2023, Carbon-Oxygen Phase Separation in MESA White Dwarf Models, arXiv2303.10110v2, Eq1–8 and §§2–5](https://arxiv.org/abs/2303.10110v2) · [Bauer 2023, fixed MESA inputs and starting models, Zenodo 7846751, CC BY 4.0](https://doi.org/10.5281/zenodo.7846751)

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

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](https://arxiv.org/abs/2303.10110v2)

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

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](https://arxiv.org/abs/2303.10110v2) · [MIST white dwarf models, arXiv2509.21717v1, cooling clocks, surface boundary and model comparisons](https://arxiv.org/abs/2509.21717v1)

## 氖把浮沉问题改了

再翻到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](https://arxiv.org/abs/2409.04419v1) · [Bédard, Blouin and Cheng 2024, fixed cooling sequences, Zenodo 10201676, CC BY 4.0](https://doi.org/10.5281/zenodo.10201676) · [Bauer 2023, Carbon-Oxygen Phase Separation in MESA White Dwarf Models, arXiv2303.10110v2, Eq1–8 and §§2–5](https://arxiv.org/abs/2303.10110v2)

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

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](https://arxiv.org/abs/2409.04419v1) · [Bédard, Blouin and Cheng 2024, fixed cooling sequences, Zenodo 10201676, CC BY 4.0](https://doi.org/10.5281/zenodo.10201676)

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

把教程的原生比较命令用于三条固定的一倍太阳质量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](https://arxiv.org/abs/2409.04419v1) · [Bédard, Blouin and Cheng 2024, fixed cooling sequences, Zenodo 10201676, CC BY 4.0](https://doi.org/10.5281/zenodo.10201676) · [Sandbox Physics, standalone white-dwarf scientific reference, original model inputs and independent tests](https://sandboxphysics.com/stories/white-dwarf-cooling/materials/white-dwarf-cooling-science.zip)

## 自己推一步

平滑相变量到了零点九，是否表示两相自由能相等的物理边界也移到了那里？

没有。自由能相等仍在二分之一处。较晚的开关让成分分离避开潜热平滑区，它的数值后果需要另行评估。

新生晶体会上浮、融化，是否意味着蒸馏阶段中央固相质量一定增加？

不一定。提出的阶段可以在产晶区整体仍为液态时搬运氖。用普通向外生长的固态前沿去画，会给错物理图景。

三个选定的平滑表面序列，能证明蒸馏能量收支闭合吗？

它们能复算表面时钟差与停留时间。完整能量检验还需要原生内部、引力、组成和损失项，这些标量文件没有提供。

### [比较原生标准序列与延迟序列](https://sandboxphysics.com/zh/labs/white-dwarf-crystallization/)

1. 选择已发表平台，把30PR30模型族、原生时间和平滑说明一起读。
2. 用独立参考，在同一光度处比较STD、部分和完整蒸馏；不要把受控核心的固相比例抄到这个视图。
3. 回到受控球壳结晶，说明还要补哪些组成与结构方程，才可以称为氖蒸馏计算。

## 原始来源

1. [Sandbox Physics, standalone white-dwarf scientific reference, original model inputs and independent tests](https://sandboxphysics.com/stories/white-dwarf-cooling/materials/white-dwarf-cooling-science.zip)
2. [Bauer 2023, Carbon-Oxygen Phase Separation in MESA White Dwarf Models, arXiv2303.10110v2, Eq1–8 and §§2–5](https://arxiv.org/abs/2303.10110v2)
3. [Bauer 2023, fixed MESA inputs and starting models, Zenodo 7846751, CC BY 4.0](https://doi.org/10.5281/zenodo.7846751)
4. [MIST white dwarf models, arXiv2509.21717v1, cooling clocks, surface boundary and model comparisons](https://arxiv.org/abs/2509.21717v1)
5. [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](https://arxiv.org/abs/2409.04419v1)
6. [Bédard, Blouin and Cheng 2024, fixed cooling sequences, Zenodo 10201676, CC BY 4.0](https://doi.org/10.5281/zenodo.10201676)

[上一篇：用原始Gaia数据比较白矮星冷却模型](https://sandboxphysics.com/zh/stories/white-dwarf-cooling/mesa-gaia-cooling-model-comparison/)