CP2K做声子谱，在混合高斯-平面波方法的框架下走的是有限位移法的路线。这条路比DFPT（密度泛函微扰理论）在实现上直观得多——就是对每个原子沿三个方向给一个小位移，算出力常数矩阵，再对角化得到声子频率。但直观不代表简单，力常数矩阵的质量是整个计算的命门。
超胞大小是第一个需要慎重决策的变量。声子计算需要一个足够大的超胞来消除周期性镜像之间的虚假力常数耦合。经验是：晶体原胞沿每个方向的重复数乘以晶格常数后，超胞的最小边长不小于10Å。对层状材料来说，面外方向的真空层还要再加。做过一个MoS₂单层的声子谱计算，4×4的超胞在低频声学支上出现了明显的虚假软化——声子频率出现了虚频。扩展为6×6后虚频消失。这个翻车在CP2K的社区里是很标准的新手体验。
位移幅度是紧跟着的第二步决策。CP2K默认的位移值是0.01 Bohr（约0.0053Å），这个值对大多数体系是安全区间的下限。如果体系含氢原子这类轻元素，偏小的位移会让轻原子的力响应模模糊糊，导致高频光学支的信噪比下降。经验是把位移调到0.015 Bohr左右，但如果出现了非简谐效应导致的虚假硬化，就退回到默认值。这个参数需要在自己的体系上试。
CP2K在声子计算中的精度瓶颈常常不是位移法本身，而是SCF收敛。声子计算需要在每个位移组态下做高精度的SCF——EPS_SCF至少设到1e-8，比常规的结构优化紧两到三个数量级。因为这个步骤是每个位移结构独立算的，SCF不收敛产生的力误差会直接污染力常数矩阵，最终让低频声学支在Γ附近出现系统性偏差。
还有一个容易被忽略的效率优化点：对称性。如果晶体的空间群不是P1，利用对称性可以把需要独立计算的位移数量从3N（N为超胞原子数）大幅削减——立方体系可能只剩下几分之一。CP2K的声子计算脚本通常会自动检测和利用对称性，但前提是你输入的超胞结构确实保留了原始晶体的对称操作。结构优化阶段如果手工移动过某个原子导致对称性被破坏，到了声子计算阶段你会发现自己在一个不需要的P1对称性下多算了十倍的位移点。
声子谱算完之后的分析有三个层次：第一层看有无虚频（判断动力学稳定性）；第二层看Γ点附近的声学支是否线性（检验长波极限下的弹性性质）；第三层看光学支的分支是否解析（判断是否有Kohn异常或软模）。三个层次都过了，声子谱才算是”算准了”而不仅仅是”算出来了”。对于做热电材料的同行，声子群速度可以从声子色散的斜率中提取出来，用于计算晶格热导率——这才是声子谱计算的最终价值落脚点，而不是画完图就结束了。
声子谱算完之后还有一个重要的后处理步骤：从声子态密度（PhDOS）提取热力学性质。CP2K加phonopy的组合可以在声子频率算完之后直接积分得到等容热容Cv、Helmholtz自由能和熵随温度的变化。这个步骤的操作量很小——就一行phonopy的命令——但产出的信息量很大。尤其对于做相稳定性研究的项目，两个竞争相的自由能差在高温下由振动熵项主导，而振动熵只能从声子频率的温度积分中得到。做零点能修正已经算是”及格线”操作了，把声子谱延展到有限温度的振动自由能才算是把声子计算的潜力用足了。
https://www.keyanxueshu.com/category/dft/
