Materials Studio计算声子谱：CASTEP模块的PHONON设置与分析

声子谱是评估材料动力学稳定性的直接判据，也是理解热传导、热容和晶格热膨胀的基础。在Materials Studio中，声子谱通过CASTEP的PHONON任务来完成，整个流程相比VASP的Phonopy路线更加集成化——从结构输入到色散曲线输出，全部在同一个图形界面内完成。但这并不意味着可以黑箱操作，几个关键参数的选择对结果精度有决定性影响。

声子计算的两条技术路线

在Materials Studio/CASTEP中，声子谱的计算有两种实现方式：

线性响应（DFPT）：又称密度泛函微扰理论，在倒空间中直接计算力常数矩阵，不需要构建超胞。这种方法对大体系更高效，但要求系统处于周期性基态，对含有缺陷或表面的体系适用性有限。

有限位移（超胞法）：构建超胞，对每个不等价原子施加小位移，计算各位移构型的Hellmann-Feynman力，再通过力常数提取声子。这种方法对超胞大小要求较高，计算量与超胞体积成正比，但对任意体系（含缺陷、合金等）均适用。

CASTEP同时支持这两种方法，在PHONON任务的设置面板中选择。对于简单体系（元素晶体、简单二元化合物），DFPT是首选；对于复杂体系，超胞法更稳健。

结构准备：声子计算对收敛的特殊要求

声子谱的计算精度对底层DFT计算的收敛要求比能量计算严格得多。原因是声子频率与力有关，而力是总能量对坐标的一阶导数——数值噪声在微分后会被放大。

计算声子谱前，结构优化必须达到比普通计算更严格的标准：

• 能量收敛精度（Convergence Threshold）建议设到 2×10⁻⁶ eV/atom（普通弛豫用 2×10⁻⁵ 即可）
• 最大力收敛到 0.005 eV/Å 以下（普通弛豫通常 0.03 eV/Å 就够）
• 最大应力分量小于 0.05 GPa

在CASTEP的Task选项中，先完成Geometry Optimization到高精度标准，再切换到PHONON任务——此时CASTEP会从优化好的结构出发，不会重新做结构弛豫（除非勾选了”Optimize structure before phonon calculation”）。

PHONON任务的关键参数

截断能（Cutoff Energy）：声子计算的截断能通常需要比普通计算高10-20%，以保证力的计算精度。如果原本用400 eV做能量计算，声子计算建议升至440-480 eV。

k点网格（k-point Sampling）：PHONON任务的k点密度直接影响力常数矩阵的精度。金属体系因为费米面的精细结构，需要比绝缘体更密的k点网格。通常建议k点间距不超过0.04 Å⁻¹。

声子超胞（Phonon Supercell）：对于超胞法，超胞的大小决定了能够描述的最长波长声子。超胞沿每个方向至少要扩展到使原子间的相互作用力衰减到可以忽略（通常需要扩展2-3倍原胞，即≥2×2×2超胞）。超胞不够大会导致色散曲线在边界k点附近出现不自然的折叠，特别是声学支在布里渊区边界处应趋近于零的行为会失真。

声子谱的物理解读

声子色散曲线是频率（或能量）对波矢q的函数，通常沿高对称方向显示。

动力学稳定性判据：如果声子谱中出现虚频（通常显示为负频率或在纵坐标零线以下的频率），说明该材料在对应波矢的扰动下是不稳定的——自发地向对应位移模式演化。布里渊区中心（Γ点）的虚频意味着铁电不稳定性；其他q点的虚频则对应相变驱动力或结构有序-无序转变。

声学支与光学支：对于n个原子的原胞，声子谱有3n支。其中3支是声学支（Γ点处频率为零），剩余3n-3支为光学支。两个声学支为横向（TA），一个为纵向（LA）；LA支的斜率决定了纵向声速，TA支斜率决定了横向声速。

声子态密度（PHDOS）：与电子态密度类似，PHDOS表示各频率下的振动模式数目，用于计算热容（Einstein模型或Debye模型的进阶版本）和声子贡献的熵。CASTEP的分析结果中，PHDOS和分原子PHDOS可以直接从Analysis面板提取。

科研学术网提供Materials Studio声子谱及热力学性质计算服务，欢迎访问 https://www.keyanxueshu.com 了解完整方案。