Theory of melting lines with a variable enthalpy of fusion

本文通过修正克劳修斯 - 克拉佩龙方程以计入由非谐固态效应驱动的熔化焓可变性，提出了一种新的熔化线解析模型，该模型导出了由基本热物理性质定义的抛物线解，并证实了近期关于液体熔化的通用模型。

要点

想象一下，你正在绘制一张固体转变为液体的地图（例如冰融化成水）。在物理学中，这张地图被称为熔化线。它展示了熔化某种物质所需的热量如何随施加的挤压程度（压力）而变化。

长期以来，科学家们使用一套简单的规则手册（克劳修斯 - 克拉佩龙方程）来绘制这些地图。但其中有个陷阱：他们假设熔化物质所需的“能量成本”（称为潜热）是恒定不变的，无论温度多高或压力多大。这套规则对于物质转变为气体（例如水沸腾）非常有效，但对于固体转变为液体来说，却是一个糟糕的猜测。因为固体和液体都具有高密度和粘性，所以其规律要复杂得多。

新观点：“可伸缩的”能量成本
本文提出了一种绘制该地图的新方法。作者安东尼·帕帕塔纳西乌（Anthony Papathanassiou）认为，熔化固体所需的能量成本并非一个固定数值；它更像是一条可伸缩的橡皮筋。随着固体被加热，其内部的原子开始剧烈抖动（非谐性），而将它们挣脱所需的能量量级，取决于原子所占的空间（体积）。

可以这样理解：

• 旧观点：想象试图将一只沉重的箱子推上斜坡。你假设箱子的重量在整个过程中始终保持完全不变。
• 新观点：这个箱子实际上是由一种特殊材料制成的，其重量会根据你移动它的速度以及你挤压它的程度而变轻或变重。要得到正确答案，你必须考虑这种变化的重量。

“体积”关联
本文运用了一个巧妙的技巧。它考察了固体受热时的膨胀程度（热膨胀）以及它所储存的热量。结果表明，在熔点附近，热容中“可伸缩”的部分，直接与固体和液体之间的尺寸差异相关联。

通过将这种“可伸缩”的能量概念代入旧有的规则手册，作者推导出了一个新的数学方程。

结果：完美的抛物线
当作者求解这个新方程时，熔化线的形状既不是直线，也不是奇怪的曲线。它实际上是一条抛物线（就像你向空中抛球时看到的 U 形轨迹）。

• 这为何令人兴奋？ 这意味着对于许多不同的材料（从氦到铁），压力与熔点之间的关系都遵循着同一条简单而弯曲的路径。
• “双重确认”：作者指出，另一位科学家（特拉琴科）最近也发现了完全相同的抛物线形状，但他们使用的是基于声波如何在液体中传播的截然不同的理论。这就像两个人从山的两侧攀登，最终在完全相同的顶峰相遇。这表明“抛物线熔化线”是自然界的基本真理，而不仅仅是一个幸运的猜测。

这张地图告诉了我们什么
本文声称，如果你了解关于某种材料的几个基本事实——它的可压缩性（体积模量）、受热时的膨胀程度以及它所储存的热量——你就可以预测其整个熔化曲线，而无需为每一个点运行昂贵的实验。

总结
本文指出：“停止假设熔化物质所需的能量是恒定的。它会随着原子的抖动和膨胀而变化。如果你考虑到这种变化，几乎所有材料的熔化线都是一条简单、可预测的曲线（抛物线），我们可以利用基本的物理属性来计算它。”

技术摘要

技术摘要：具有可变熔化焓的熔化线理论

问题陈述
与升华曲线和沸腾曲线相比，固 - 液相界（熔化线，MLs）的热力学描述在历史上一直具有挑战性。利用克劳修斯 - 克拉佩龙（CC）方程进行的常规推导通常依赖于恒定熔化潜热（$\delta H$）的近似，以获得解析解。虽然这一假设有助于推导气相转变的通用方程，但它无法捕捉固 - 液转变的复杂热力学，因为在该转变中，两相均为致密且强相互作用的凝聚态。现有模型往往仅关注固体的结构不稳定性（例如林德曼判据、西蒙 - 格拉策尔方程），或者需要关于液态理论的明确共识，而后者目前仍难以企及。因此，目前缺乏一种源自两相热力学框架、并考虑潜热可变性的熔化线通用解析方程。

方法论
本研究通过修改 CC 方程以纳入沿熔化线变化的熔化焓（$\delta H$），开发了一种两相解析模型。该方法根植于近期关于高温下晶体固体非谐性的理论进展。具体而言，该方法利用了一个发现：固体的等压热容（$C_P$）由一个与温度相关的振动分量和一个与体积相关（膨胀性）的非谐分量组成。

推导过程如下：

1. 可变潜热：熔化焓被表示为共存液相（$v_L$）和固相（$v_S$）比容的函数，该函数源自固体热容的体积依赖分量（$C_{TE} \propto \varepsilon \beta$）。由此得出关系式 $\delta H = \varepsilon \ln(v_L/v_S)$，其中 $\varepsilon$ 是一个与热容和热膨胀之比相关的固态参数。
2. 修正的 CC 方程：将此可变 $\delta H$ 代入 CC 方程（$dP/dT_m = \delta H / \delta v$），形成一个微分方程，其斜率取决于两相的比容。
3. 二阶微分方程：通过对修正后的 CC 方程关于温度求导，并应用热膨胀（$\beta$）和等温体积模量（$B$）的热力学关系，导出了支配熔化线的二阶线性非齐次常微分方程。
4. 边界条件与近似：解受到物理边界条件的约束：零压熔点（$T_0, 0$）和三相点（$T_c, P_c$）。分析假设参数 $\varepsilon$ 对压力的依赖性较弱，且项 $d$（与热膨胀和体积模量之差相关）显著大于项 $b$（与体积差的对数导数相关）。

主要结果
所导出的微分方程的解给出了熔化压力随温度变化的近似抛物线函数 $P(T_m)$。

• 抛物线标度：在 $d \gg b$ 的条件下，二阶微分方程简化为 $d^2P/dT_m^2 \approx d > 0$。通解为二次函数：$P(T_m) \approx \frac{1}{2}d T_m^2 + c_3 T_m + c_4$。
• 参数定义：该抛物线函数的系数完全由基本热物理性质定义，包括体积模量（$B_L, B_S$）、热膨胀系数（$\beta_L, \beta_S$）、比容（$v_L, v_S$）以及固体的等压热容（$C_{P,S}$）。未引入任何经验拟合参数。
• 斜率预测：熔化线的初始斜率推导为 $dP/dT_m \approx \beta_L B_L$。该结果表明，斜率由液体的热膨胀系数与其等温体积模量的乘积决定。
• 凸性与单调性：导出的函数在 $P \geq 0$ 和 $T \geq T_0$ 的所有范围内保持严格单调递增且凸的轮廓，这与熔化曲线的物理预期一致。

意义与主张
本文主张，该推导提供了一个通用的熔化线解析框架，克服了恒定潜热近似的局限性。其主要意义在于两个不同理论基础的汇聚：

1. 理论汇聚：此处从固体非谐性和热容分析推导出的抛物线标度律以及斜率的具体预测（$dP/dT_m \approx \beta_L B_L$），与 Trachenko 从液体声子理论（PTL）推导出的近期通用模型惊人地一致。
2. 两相本质：与早期仅关注固体不稳定性的一相模型不同，该方法是一种真正的两相理论，它通过比容和体积模量的差异明确纳入了液相的热力学。
3. 实验一致性：作者指出，预测的抛物线性质和斜率大小与各种材料的实验数据一致，包括 Ar、He、H$_2$、H$_2$O、In 和 Fe，正如在 Trachenko 工作的背景下先前所讨论的那样。

研究得出结论，熔化曲线的通用抛物线性质是一个稳健的特征，得到了独立理论途径的支持，为固 - 液转变的热力学提供了统一的视角，而无需依赖任意的经验拟合。