如何从声子算热力学性质？

说明：本文华算科技主要介绍声子频率和声子态密度怎样生成零点能、振动内能、Helmholtz 自由能、振动熵和定容热容，以及准谐近似中的体积依赖怎样产生热膨胀、定压热容和有限温相稳定结果。

声子热力学量来源与声子态密度定义

结构优化得到的是势能面上的一个参考构型。原子在平衡位置附近偏移时，能量对位移的二阶导数构成力常数矩阵；对质量加权后的动力学矩阵求本征值，便得到每个波矢 q 和支序 ν 的频率 ω_qν。含 N 个原子的原胞共有 3N 条声子支，其中三条在 Γ 点附近属于声学支，其余为光学支。

声子色散保留了频率随波矢变化的信息，热力学积分则常把全部 q 点和支序压缩成声子态密度 g(ω)。若按一个原胞计数，g(ω) 对全频区的积分应为 3N；若结果按化学式单元或摩尔给出，态密度、能量和热容必须采用同一计量基准。归一化差一个原胞倍数，所得热容量级也会按同一倍数偏离。

实际积分使用有限 q 网格，经插值或展宽形成连续的 g(ω)。低温热容集中采样 Γ 点附近的低频声学模，粗网格会漏掉这部分相空间，即使高温热容已经收敛，低温 T³ 区仍可能偏差较大。检查对象应同时包含总态数积分、低频态密度形状和温度曲线收敛。

谐近似把每个正常模视为频率固定的量子谐振子。该模式的能级间隔为 ℏω，温度通过无量纲量 x=ℏω/(k_BT) 控制热占据数。低频模式具有较小的 x，在较低温度便出现可观占据；高频模式直到 k_BT 接近 ℏω 时才明显参与熵和热容。由此，同一份 g(ω) 会被不同温度核函数重新加权，并生成一组随温度变化的振动热力学量。

LiGa₂Ir 的色散和分波态密度给出了这种转换：色散中的每个本征频率都进入右侧频率统计，高频区主要含 Li 振动，较低频区含较多 Ga、Ir 分量。热力学积分并不读取曲线颜色本身，而是读取各频段内的状态数和频率。元素分波态密度适合追踪某一温区由哪些原子振动主导；总态密度才对应整套振动自由度的配分函数。

振动态密度怎给自由能熵定容热容？

在给定体积下，所有模式的配分函数相乘，取对数后便转化为对态密度的积分。零温仍保留每个模式的基态能量，振动零点能写成：

E_ZP = 1/2 ∫ ℏωg(ω) dω

有限温度振动内能在零点项上加入玻色占据，Helmholtz 振动自由能则含有热激发带来的熵项：

U_vib(T) = ∫ g(ω)ℏω[1/2 + 1/(e^x−1)] dω

F_vib(T) = ∫ g(ω){1/2ℏω + k_BT ln[1−e^−x]} dω

对 F_vib 作温度导数得到 S_vib=−(∂F_vib/∂T)_V，对 U_vib 作温度导数得到定容热容。写成频率积分后，单个模式对 C_V 的权重为 k_Bx²e^x/(e^x−1)²。温度导数改变了各频段的权重：零点能偏重高频模式，低温热容偏重声学低频模式，振动熵则随可热激发状态数增加而上升。

LiGa₂Ir 的红色曲线直接由声子态密度积分得到。低温区曲线迅速离开零值，反映低频声学模开始热激发；温度继续升高后，高频 Li 模逐步贡献热容。Debye–Einstein 表示分别近似低频连续谱和一组较集中的高频模，完整态密度积分则保留了各频段的真实分布。

同一套频率还能生成振动熵、Helmholtz 自由能和定容热容曲线。总 Helmholtz 自由能常写为静态电子能与 F_vib 之和，不同结构必须采用一致的原子数和静态能参照。自由能随温度下降的幅度由 TS_vib 控制，频谱较软的结构往往拥有较大的振动熵，有限温相序因而可能不同于静态总能排序。

比较两个相的稳定温区时，应对每个相分别计算 E₀+F_vib，并保持化学计量和能量零点一致。相间自由能差常小于各自的振动自由能，低频软模、零点能或少量未收敛的虚频都可能移动交点温度。相变温度的数值应伴随 q 网格、超胞和体积采样的收敛结果。

低温极限还带有明确的频谱特征。三维普通晶体的声学支在 Γ 点附近近似线性色散，对应低频态密度随 ω² 增长，积分后得到 C_V∝T³。温度远高于主要声子能标后，每个模式的热容趋近 k_B，一个含 N 个原子的原胞趋近 3N k_B。这一高温平台是经典 Dulong–Petit 极限，并未包含热膨胀产生的 C_p−C_V 差值。

苯到并五苯的结果呈现了频谱位置对升温过程的影响。低频晶格模较多的体系在低温获得更大的归一化热容；分子内高频振动的 x 仍较大时，对热容贡献接近冻结。随温度升高，更多模式进入热激发区，C_V/(3Nk_B) 逐步接近 1。实验若测得定压热容，高温区还会含体积膨胀和显式非谐效应，数值可高于谐近似的 C_V。

热膨胀和定压热容为何需算多体积？

单一体积的谐声子计算给出 F_vib(T)、S_vib(T) 和 C_V(T)，体积在整段温区保持固定。热膨胀涉及平衡体积随温度移动，必须获得频率对体积的响应。准谐近似在若干体积上分别优化内部坐标并计算谐声子，形成 ω_qν(V) 和 F_vib(T,V)。

静态能和振动自由能相加得到 F(T,V)=E₀(V)+F_vib(T,V)。给定外压 p 后，对 F+pV 关于体积取最小值，得到 V(T,p)；体膨胀系数来自 α_V=(1/V)(∂V/∂T)_p。频率随体积软化或硬化的方向决定各模式对热压和膨胀的贡献，静态状态方程则抵抗体积偏离零温平衡值。

Si 的三条态密度曲线随体积发生频移，g(ω,V) 已成为体积的函数。把每个体积的态密度分别代入振动自由能积分，便得到一组温度相关的 F–V 曲线。单相准谐模型中，各模式的 Grüneisen 参数和定容热容权重共同决定膨胀系数；这里采用的 Si 模型还含 Si-I 与亚稳 Si-IV 微观态的构型熵，后面的线膨胀曲线同时包含该项。

定压热容可沿平衡体积路径对焓求导，也可由 C_p−C_V=α_V²B_TVT 计算，其中 B_T 为等温体模量。该式要求 α_V、B_T、V 与热容量采用一致的每原胞、每摩尔或每质量基准。单体积声子输出中的 C_V 缺少膨胀功，与常压量热得到的 C_p 比较时应补上这一项。

立方晶体常用一组等比例缩放体积。层状、正交或低对称晶体可能沿不同晶轴呈现相反的膨胀，单一体积变量只能给出体积平均响应。此时可在独立应变分量上构建自由能面，并由平衡晶格参数求线膨胀张量。体膨胀系数与各晶轴膨胀系数的实验对象不同，对照数值时应保持定义一致。

何情况谐声子热力学偏离真实材料？

谐配分函数要求参考构型附近的二阶势能为正定。真实虚频对应负本征值，势能沿该正常模向下弯曲，围绕该构型建立的谐振子配分函数失去物理定义。直接删除虚频会删去该自由度的熵和热容，得到的自由能缺少确定的统计系综。Γ 点附近很小的声学负频若随超胞、q 点和声学和规则处理而消失，可归入数值误差；持续存在的软模则指向结构畸变、相变或温度稳定化。

准谐近似允许频率通过体积随温度变化，却仍把每个固定体积上的振动写成独立谐振子。显式声子—声子相互作用未进入 ω(V)。接近熔点、强软模、超离子迁移、转动无序或动态稳定高对称相中，频率本身会随温度重整化，热膨胀也可能反馈到非抛物势能区。

MgO 在较低温区的准谐结果与实验接近，约 1000 K 以上体积偏差逐渐增大，定压热容和膨胀系数也出现系统偏离。采用温度依赖的声子准粒子频率后，高温曲线更接近实验数据。此类差异来自固定体积频率随温度发生的非谐重整化，增加体积采样点并不会补回这一物理项。

声子积分给出的是振动对总自由能和热容的贡献。金属的电子热激发、磁性材料的磁熵、多组分固溶体的构型熵、缺陷浓度随温度的变化均属于另外的自由度。相稳定计算应在统一的原子数、组分、压力和能量零点下，把适用的电子、磁性与构型项加入 F(T,V)。若软模随温度重整化，则应采用自洽声子、温度依赖有效势、从头算分子动力学或热力学积分得到相应的有限温自由能。