🔬 condensed matter

The origin of complex behavior of liquid carbon: an insight from computer simulation

本文利用机器学习势函数 GAP-20 对液态碳进行分子动力学模拟，揭示了该模型下碳具有相对较低的临界温度，这一发现可能影响石墨熔点实验测量结果的解读。

原作者： Yu. D. Fomin

发布于 2026-03-02

📖 1 分钟阅读☕ 轻松阅读

原作者： Yu. D. Fomin

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

这篇论文就像是在解决物理学界的一个“百年悬案”：石墨到底在多少度融化？

想象一下，石墨就是铅笔芯的主要成分。科学家们早就知道，如果你把石墨加热到几千度，它就会变成液态碳。但是，不同的实验室测出来的结果却大相径庭：有的说 4000 度就化了，有的说要到 6700 度才化。这就像大家测量一座山的高度，有人说是 1000 米，有人说是 3000 米，这中间的差距太大了，让人很困惑。

这篇论文的作者（来自俄罗斯科学院）决定用超级计算机模拟（就像在电脑里造一个微缩宇宙）来搞清楚到底发生了什么。他们发现，问题的根源可能在于液态碳的“脾气”太古怪了。

以下是用通俗语言和比喻对论文核心内容的解读：

1. 核心发现：石墨融化时，可能直接“跳”进了临界区

作者使用了一种名为 GAP-20 的高级人工智能模型（就像给原子装上了超级智能的导航系统）来模拟碳原子的行为。

传统观点：大家以为石墨融化后变成液体，然后随着温度继续升高，液体慢慢变成气体。
新发现：作者发现，石墨的熔点（开始融化的温度）和临界点（液体和气体分界消失的温度）靠得太近了！
- 比喻：想象你在爬一座山（升温过程）。通常，你会先到达山腰（融化成液体），然后继续往上爬，最后到达山顶（变成气体）。但作者发现，石墨这座山太特殊了，山腰和山顶几乎重叠在一起。当你把石墨加热到融化时，它瞬间就处于一种“既像液体又像气体”的混乱状态（物理学上叫超临界流体）。

2. 为什么实验结果会打架？（混乱的“临界点”效应）

既然熔点这么接近临界点，为什么测不准呢？因为在这个区域，物质变得极其敏感和情绪化。

比喻：想象一个极度敏感的气球（临界点附近的碳）。
- 在普通温度下，你轻轻捏一下（微小的实验误差），气球只会稍微变形。
- 但在临界点附近，你稍微吹一口气（哪怕压力或温度有一点点测量误差），气球就会剧烈膨胀或收缩，甚至爆炸。
后果：
- 热容量（吸热能力）乱变：在临界点附近，碳吸收热量的能力会突然飙升，像过山车一样。以前的实验假设碳吸热是恒定的，这就像假设气球无论怎么吹体积都不变，结果当然算错了。
- 密度乱跳：因为热胀冷缩效应在这里被无限放大，一点点温度变化就会导致密度剧烈波动。

结论：不同的实验团队，因为测量时的微小误差（比如温度计差了 1 度，或者压力计差了 1 个大气压），导致他们看到的“融化状态”完全不同。有的测到了剧烈波动的临界区，有的测到了稍微稳定一点的区域，所以得出了 4000K 到 6700K 这种天差地别的结果。

3. 计算机模拟的“局限性”与“启示”

作者很诚实，他们承认：

模型不是完美的：他们用的 GAP-20 模型主要是为了模拟“固态”和“液态”碳设计的，用来模拟“气态”或“临界点”可能不够精准。就像用造房子的砖块模型去模拟云朵，虽然能看出个大概，但细节可能有偏差。
计算出的温度偏低：他们的模拟显示临界点大概在 4235K，这比很多实验数据低。
但是，方向是对的：虽然具体数字可能有偏差，但定性结论（即“熔点和临界点靠得太近，导致物质行为极度不稳定”）极有可能是正确的。

4. 总结：我们学到了什么？

这篇论文并没有直接给出一个“标准答案”说石墨到底在多少度融化，但它给出了一个全新的视角：

以前的思路：试图通过更精确的测量来找到那个唯一的“融化温度”。
现在的思路：也许根本不存在一个单一的、稳定的“融化温度”。因为石墨一融化，就掉进了一个物理性质剧烈震荡的“雷区”。

一句话总结：
石墨融化不像水烧开那样平静，它更像是在悬崖边跳舞。因为熔点离“液体变气体”的临界点太近，导致物质在融化瞬间变得极度不稳定。以前的实验之所以结果打架，是因为大家没意识到自己正站在悬崖边，稍微动一下（测量误差），结果就天差地别。未来的研究需要考虑到这种“临界点的混乱”，才能解开石墨融化的谜题。

以下是基于 Yu. D. Fomin 所著论文《液体碳复杂行为的起源：计算机模拟的见解》（The origin of complex behavior of liquid carbon: an insight from computer simulation）的详细技术总结：

1. 研究背景与问题 (Problem)

石墨的熔化是物理学中长期存在的难题。尽管自 20 世纪初以来已有大量实验，但石墨熔点的测定结果存在巨大的离散性，范围从 4000 K 到 6700 K 不等。

核心矛盾：现有的实验数据差异巨大，无法达成共识。
现有理论的局限：
- 关于碳的气 - 液临界点（Critical Point）参数（如临界温度 $T_c$ 、临界密度 $\rho_c$ ）的估算非常不确定。
- 经典文献（如 Leider 等人，1973）基于一些假设（如熔化时体积跳跃为 20%、液态碳的热膨胀系数与石墨相同）估算出 $T_c \approx 6810$ K。
- 近期实验（如 Kondratyev 等人，2019）在 7000 K 下测量“液态碳”密度，但该温度可能已超过某些估算的临界温度，意味着物质处于超临界流体（SCF）状态而非纯液态。
- 实验分析中常忽略气相（升华碳）的影响，且常假设液态碳的热容和热膨胀系数为常数，这在临界点附近可能不成立。
研究目标：通过分子动力学模拟，重新评估碳的临界点位置，并探讨临界现象（如热力学响应函数的极大值）是否导致了石墨熔化实验数据的巨大分歧。

2. 方法论 (Methodology)

模拟系统：在周期性边界条件的立方盒子中模拟了 4096 个碳原子。
势函数：使用 GAP-20（基于机器学习的势函数），该势函数旨在重现密度泛函理论（DFT）计算结果，主要用于描述碳的凝聚相。
模拟设置：
- 系综：正则系综（NVT）。
- 时间步长：2.0 fs。
- 平衡时间：10 ps；生产运行时间：20 ps。
采样范围：
- 密度 ( $\rho$ )：0.05 g/ml 至 2.0 g/ml（步长 0.05 g/ml）。
- 温度 ( $T$ )：4000 K 至 5000 K（步长 200 K），并额外计算了 4100 K、6000 K 和 7000 K 的等温线。
数据处理：
- 采集压力 ( $P$ ) 和内能 ( $U$ ) 数据。
- 使用多项式拟合 ( $X = \sum a_{ij} T^i \rho^j$ ) 构建状态方程（EoS）。
- 基于拟合数据计算热力学响应函数：定压热容 ( $c_P$ )、热膨胀系数 ( $\alpha_P$ ) 和等温压缩率 ( $\beta_T$ )。
- 通过寻找等温线上的拐点确定临界点，并追踪响应函数极大值的位置（即 Widom 线）。

3. 主要结果 (Key Results)

临界点估算：
- 基于 GAP-20 模型，计算出的碳气 - 液临界点参数为：
  - $T_c \approx 4235$ K
  - $\rho_c \approx 0.1887$ g/ml
  - $P_c \approx 106.61$ bar
- 该临界温度远低于文献 [4] 估算的 6810 K，甚至低于部分实验测得的石墨熔点。
相态行为：
- 在 $T=4100$ K 和 $\rho=1.55$ g/ml 附近，观察到系统出现结晶迹象（向石墨相转变），表明临界点与熔化线在相图上非常接近。
- 在 $T=4200$ K 附近，系统处于超临界状态，表现出强烈的热力学涨落。
热力学响应函数的极大值（Widom 线）：
- 在临界点附近， $c_P$ 、 $\alpha_P$ 和 $\beta_T$ 均出现显著极大值。
- 定压热容 ( $c_P$ )：在临界点附近可达 $25 k_B$ /粒子，并随温度升高迅速下降。这直接反驳了“液态碳热容为常数”的假设。
- 热膨胀系数 ( $\alpha_P$ )：在临界点附近远大于石墨的热膨胀系数。
- 等温压缩率 ( $\beta_T$ )：极大值仅出现在略高于临界温度的区域，并迅速消失。
与实验数据的对比：
- 近期实验（Ivanov 等人，2026）测得石墨熔点约 6000 K，且液态碳的等容热容 $c_V$ 较高（约 4.19 $k_B$ ）。
- 模拟结果（ $T_c \approx 4235$ K）虽然数值偏低，但定性趋势表明：石墨熔化后系统可能立即进入强涨落的超临界流体区域，或者处于气 - 液两相区。

4. 关键贡献 (Key Contributions)

首次计算 GAP-20 模型下的碳 Widom 线：揭示了碳在超临界区域热力学响应函数的极大值分布。
提出实验数据分歧的新解释：
- 石墨熔点可能非常接近（或低于）真实的气 - 液临界温度。
- 当石墨熔化时，系统可能直接进入具有强热力学涨落的超临界流体（SCF）状态。
- 这种状态下的热容剧烈变化和热膨胀系数异常，导致实验测量（如量热法、密度测量、高温计读数）产生巨大误差。
- 实验分析中若忽略气相影响或假设热力学参数恒定，将导致熔点估算的严重偏差。
指出计算方法的局限性：承认 GAP-20 模型主要针对凝聚相优化，可能低估了真实碳的临界温度（真实 $T_c$ 可能更高，如 6000 K 以上），但模型揭示的“临界点靠近熔化线”这一物理图像可能是正确的。

5. 意义与结论 (Significance & Conclusion)

理论意义：该研究强调了在分析高温碳物性时，必须考虑临界现象和超临界流体的强涨落特性。传统的“液态”假设在石墨熔化附近可能失效。
实验指导：解释了为何不同实验组测得的石墨熔点差异巨大（4000 K - 6700 K）。这种差异可能源于实验条件（压力、温度）是否恰好落在临界点附近的强涨落区域，以及是否考虑了气相和热容的非恒定性。
未来展望：
- 需要构建能够同时精确描述碳的气相、液相和石墨相的新势函数模型。
- 需要更精确的量子模拟和实验来验证真实的临界点位置。
- 目前的模拟结果表明，即使临界点参数存在数值偏差，**“熔化点与临界点接近导致热力学性质剧烈波动”**这一机制是解释石墨熔化实验数据混乱的关键线索。

总结：本文通过分子动力学模拟指出，石墨熔化实验数据的巨大离散性可能源于熔化后的碳处于或接近气 - 液临界点，导致热力学响应函数（如热容、压缩率）出现极大值和强涨落。现有的实验分析往往忽略了这些临界效应，从而导致了测量结果的不可靠。