计算声子谱：从晶格动力学到热力学性质的量化路径

计算声子谱是评估晶体动力学稳定性与预测热力学性质的关键手段。当项目需要判断一种预测的新型晶体结构能否在室温下稳定存在时，声子谱中的虚频是否出现就成了决定性的判据。

计算声子谱有两条主流路线：密度泛函微扰理论（DFPT）和有限位移法。DFPT通过线性响应直接计算力常数矩阵，物理上更严格，但计算量随超胞尺寸增长显著。有限位移法通过对每个原子施加微小位移并计算受力变化来构建力常数，实现简单且与任意DFT代码兼容。项目最终认定有限位移法更适合当前体系——并非因为它更精确，而是它配合Phonopy工具链的自动化程度更高，能够系统性管理位移模式与力常数提取。

项目使用VASP进行力计算，配合Phonopy后处理。非解析项校正（NAC）对极性晶体的声子谱至关重要——缺乏NAC校正时，Γ点附近光学声子支的LO-TO劈裂会完全消失。项目对含B-N键的体系启用了NAC，Born有效电荷通过DFPT单独计算获得。

超胞尺寸在声子谱计算中影响最大。力常数矩阵的作用范围随原子间距衰减，超胞必须足够大以容纳非零力常数的全部作用对。项目从2×2×2超胞起步逐步扩大，发现2×2×2超胞下声学声子在Γ点附近出现2.3 THz的虚频——这是一个危险的信号，意味着结构可能动力学不稳定。然而扩大至3×3×3超胞后，虚频消失，声子谱全频段为正。这个对比说明，超胞尺寸不足会引入虚假虚频，导致对结构稳定性的误判。

项目最终采用4×4×4超胞，对应包含256个原子的体系。位移幅度设为0.01 Å，这一取值参考了Togo在Phonopy开发文档中的建议值——过大的位移引入非简谐效应污染力常数，过小则力信号被数值噪声淹没。每个位移模式的力计算采用与结构优化一致的精度参数：ENCUT=520 eV，EDIFF=10⁻⁸ eV，k点密度在超胞中相应缩减以保持总计算量可控。

计算声子谱最具诊断价值的场景是虚频分析。项目在处理一种高温相变材料时，发现声子谱在M点出现明显的虚频支，频率为-1.5 THz。虚频意味着晶体在该波矢方向存在能量更低的畸变构型——系统倾向于沿这个模式发生结构相变。项目沿虚频本征向量对原子坐标施加扰动并重新优化，果然得到一个对称性更低的稳定构型，能量比原结构低8 meV/atom。

这一发现让项目团队既兴奋又警觉：声子谱不仅仅是稳定性判据，它还能指向未知的结构相变路径。正如Physical Review Letters多项研究所揭示的，声子虚频分析已成为预测结构相变的重要工具。但项目也注意到，虚频的物理来源需要谨慎甄别——超胞尺寸不足、k点采样稀疏或收敛精度不够都可能产生数值假虚频，只有排除这些因素后才能将虚频归因于真实的结构不稳定性。

声子态密度不仅是结构信息，还是热力学性质的输入。项目通过Phonopy计算了晶格振动对自由能、熵和热容的贡献。德拜温度从声子谱低频声学支的声速拟合得到，结果为385 K，与实验测量值392 K偏差不到2%。恒容热容Cv在300 K时趋近于经典Dulong-Petit极限，低温区遵循T³规律，声子谱计算给出的德拜T³系数与实验数据定量吻合。

值得警醒的是，声子谱计算基于简谐近似，无法描述高温下的热膨胀效应。项目在500 K以上的热容预测中观察到与实验值的系统性偏离，偏差随温度升高而扩大。准谐近似（QHA）通过在不同体积下重复声子谱计算部分修正了这一问题，但计算成本翻了五倍。对于涉及显著非简谐效应的体系，如软模相变材料，QHA本身也可能失效。

回过头看，计算声子谱的价值远超出”有没有虚频”这一简单判断。这个项目在超胞收敛、NAC校正与虚频溯源中建立的流程，使声子谱从单一的结构验证工具拓展为连接晶格动力学与宏观热力学性质的量化桥梁，每个参数设定都有明确的物理动机和收敛验证支撑。