声子谱看起来像一张曲线图，背后其实是一个很敏感的力常数问题。超胞太小，原子位移只感受到近邻相互作用；位移太大，会混入非线性响应；位移太小，又容易被数值噪声淹没。声子谱计算的核心不是把 phonopy 跑通，而是让位移产生的力能够代表真实的谐力常数。
这也是为什么同一个结构，有人算出无虚频，有人算出低频虚频。差别未必来自材料本身，可能来自超胞尺寸、位移幅度、收敛精度、二维真空层或是否充分优化。声子谱适合做稳定性检查，但前提是这些计算边界先站住。
声子谱为什么能判断稳定性？
声子谱可以先简单理解成“晶格振动模式的地图”。在一个已经优化好的结构里，原子不会完全静止，而是会围绕平衡位置做很小的振动。VASP 和 phonopy 做的事，不是模拟一段真实时间里的原子运动，而是从小位移产生的力响应里，推回这些振动模式对应的频率。
如果某个振动模式的频率是正常的正值，说明结构沿这个方向的小扰动会被拉回平衡位置；如果出现虚频，就意味着结构沿某个方向可能会滑向更低能构型，或者计算里的数值误差把低频模式推成了负值。声子谱判断的是结构在小位移扰动下的动力学稳定性，它和内聚能、形成能不是同一种稳定性证据。
图1 可以先帮我们把概念落到具体结构上：优化后的二维单层只是起点，接下来还要检查它在小振动扰动下会不会保持在原来的构型附近。
图1 1T′-HfCl₂ 单层结构与稳定性研究示意。声子谱可以理解为检查优化结构在小扰动下是否会保持原构型的一类工具。DOI: 10.1021/acsomega.4c10560。
为什么声子谱会牵涉到超胞？
声子谱不是只看某一个原子自己怎么动，而是看一个原子被轻轻推开后，周围原子受到怎样的力。这个力的响应会跨过多个原子，有时还会延伸到更远邻。若只用最小原胞，位移原子会和自己的周期镜像挤得太近，远程相互作用被截断，低频声学支和虚频判断就可能失真。
所以，超胞不是为了把图画得更大，而是为了让位移产生的力响应有足够空间展开。原胞决定你从哪个结构出发，超胞决定你保留多远的力常数相互作用。理解这一点后，再看位移幅度、力常数和虚频，就不会觉得这些参数是孤立的。
为什么不能直接用最小原胞？
有限位移法的思路，是轻轻推一下某个原子，然后看所有原子受到的力怎样响应。这个响应不是只发生在最近邻原子上，尤其是离子性材料、层状材料和低维材料，力常数可能延伸到更远邻。超胞的作用，就是给这些相互作用留出空间。
图2 给出了 1T′-HfCl₂ 单层的原胞结构。做声子谱前，要先看原胞大小、对称性和真空方向；如果直接用很小的超胞，周期镜像中的位移原子会彼此影响，低频声学支最容易被带偏。超胞选择本质上是在截断力常数相互作用范围，不是机械地把原胞复制几次。
图2 1T′-HfCl₂ 单层的俯视和侧视结构。声子超胞要从原胞、晶格对称性和相互作用范围出发选择。DOI: 10.1021/acsomega.4c10560。
把图2 这类原胞扩成声子超胞后，有限位移法真正读取的是“一个原子被推开后，其他原子受到多大力”的响应。这个响应可以写成力常数，后面所有声子频率都从它出发。
Φijαβ = – ∂Fiα / ∂ujβ 式(1)
式(1)里，Φijαβ 是原子 i 和原子 j 之间的力常数，Fiα 是 i 原子在 α 方向的力，ujβ 是 j 原子在 β 方向的位移。超胞越大，能保留的远程力常数越完整，但计算量也会明显增加。
多大的超胞才算够？
没有一个所有体系通用的超胞尺寸。三维小原胞材料常从 2×2×2 或 3×3×3 测起；二维材料常从 4×4×1、5×5×1 这类面内超胞测起；slab 和吸附体系则要根据表面周期、吸附物间距和固定层设置单独判断。
判断超胞够不够，不是看计算是否完成，而是看声子谱的关键特征是否稳定：低频声学支是否合理，虚频是否随超胞增大消失或收敛，主要光学支位置是否变化很小。声子超胞要用谱线收敛来定，而不是只按原胞大小猜。
位移太大会发生什么？
有限位移法默认原子只做很小的谐振动。如果位移幅度太大，计算到的力会包含明显的非谐贡献，力常数就不再代表平衡位置附近的二阶响应。对软材料、层状材料、吸附物和弱相互作用体系，这种影响更明显。
图3 展示了 1T′-HfCl₂ 的声子谱和热力学性质。无明显虚频说明该结构在对应计算设置下没有被小位移扰动推向不稳定方向；但这个结论仍然依赖超胞、位移和力收敛设置。无虚频是计算结果，不是脱离参数边界的永久标签。
图3 1T′-HfCl₂ 单层的声子谱和热力学性质。无明显虚频通常说明该结构在谐近似下具有动力学稳定性。DOI: 10.1021/acsomega.4c10560。
位移太小又会怎样？
位移太小时，原子受力差值很小，容易被电子步收敛误差、力精度和数值噪声干扰。最终表现为低频支抖动、Gamma 点附近声学支不贴零，或者出现很浅的虚频。常见位移幅度可以从 0.01 Å 左右开始测试，但弱键、重元素或软模体系要看力响应是否稳定。
F ≈ -KΔu 式(2)
式(2)是谐近似下力和位移的线性关系。Δu 太大，线性近似不再好；Δu 太小，F 的数值误差会占比变大。比较 0.005 Å、0.01 Å、0.02 Å 下的低频支，是判断位移幅度是否合适的简单办法。
为什么优化精度要比普通结构更严？
声子谱用的是力的二阶信息，所以对结构优化残余力非常敏感。一个普通吸附能计算里能接受的小残余力，放到声子谱里可能变成低频虚频。做声子前，最好先把结构优化到更严格的力标准，并在静态力计算中使用较高 ENCUT、较密 k 点和较严 EDIFF。
图4 对比了不同 ZrO₂ 晶相的结构、声子谱和声子态密度。不同相的软模特征不一样，说明虚频可能是数值问题，也可能是真实结构转变的信号。判断虚频前，先要确认力常数矩阵来自足够干净的力数据。
图4 四方和立方 ZrO₂ 的结构、谐声子谱与声子态密度。不同晶相的软模和虚频位置会改变稳定性判断。DOI: 10.1021/acs.cgd.2c01458。
Gamma 点附近小虚频一定要重算吗？
二维材料和 slab 模型常在 Gamma 点附近出现很浅的低频虚频，尤其是 ZA 支。这类虚频可能来自有限超胞、真空层、数值噪声或二维声学支处理方式。不能一看到负频就直接判定结构失稳，也不能完全忽略。
比较靠谱的做法，是看虚频大小、位置和对应振动模式。若虚频很浅、集中在 Gamma 附近，并且随超胞增大和收敛加严而减小，通常更像数值伪影；若虚频出现在特定 q 点，且对应清晰结构畸变，就要沿该模式重构并重新优化。虚频要结合模式位移图判断，只看曲线正负不够。
真空层和固定层会不会影响声子谱？
二维材料的 c 方向是真空，不应该像三维晶体那样复制超胞。声子超胞通常只在面内扩展，比如 4×4×1 或 5×5×1。若真空层太薄，周期镜像相互作用会影响低频模；若 slab 底层固定，声子计算的活动原子范围要和你的物理问题一致。
吸附体系更麻烦。若只关心吸附物的 ZPE 和振动熵，可以固定较深层基底，只让吸附物和表面几层参与频率计算；若要讨论整个 slab 的动力学稳定性，就不能只算一个局部频率。声子超胞的活动原子范围要和结论对象一致，否则局部结论会被误写成整体稳定性。
虚频可能指向真实重构吗？
图5 展示了二维 β-In₂Se₃ 中显著虚频、势能面和电荷密度波结构。这里的虚频不是简单的数值波动，而是和特定结构畸变、势能面和电荷密度波相关。遇到这种情况，应该沿虚频模式扰动结构，重新优化低对称相，再比较能量和声子谱。
图5 二维 β-In₂Se₃ 的虚频声子、势能面和电荷密度波结构。虚频模式需要结合位移方向和重构路径解释。DOI: 10.1038/s41467-024-54205-1。
这一步很关键：如果扰动后结构能量下降并且虚频消失，原结构很可能只是高对称鞍点或不稳定相；如果扰动后又回到原结构，或者虚频随数值设置减弱，就更可能是计算边界问题。
INCAR 参考示例可以怎样写？
SYSTEM = phonon force calculation
IBRION = -1
NSW = 0
ENCUT = 520
EDIFF = 1E-8
ISMEAR = 0
SIGMA = 0.05
LREAL = .FALSE.
这组设置用于位移结构的静态力计算。真正使用时，要保持所有位移结构的 INCAR、KPOINTS、POTCAR 完全一致，并确认每个 OUTCAR 都正常收敛。若其中一个位移结构电子步没收敛，整套力常数都会被污染。
检查顺序可以怎样安排？
实操上可以这样做：先把原胞或 slab 充分优化；再选一组超胞尺寸和位移幅度；生成位移结构后逐个跑静态力；最后用 phonopy 汇总 FORCE_SETS 或 force constants，画声子谱。第一轮只看是否有明显错误，第二轮再做超胞和位移幅度收敛。
D(q)eqν = ωqν2eqν 式(3)
式(3)里，D(q) 是动力学矩阵，ωqν 是 q 点、第 ν 个声子模式的频率，eqν 是对应本征矢。声子谱的每一条支，最终都来自力常数矩阵在不同 q 点的对角化。
图6 提醒我们，谐声子谱不是全部温度物理。若体系存在强非谐性、离子迁移或高温相稳定问题，只靠 0 K 谐声子谱会不够，需要进一步考虑非谐声子、分子动力学或准谐近似。声子谱是稳定性入口，不是所有温度效应的终点。
图6 非谐性与离子迁移示意。谐声子谱是稳定性入口，但温度效应、软模和非谐振动可能改变高温相判断。DOI: 10.1021/acs.cgd.2c01458。
要点：声子超胞决定力常数相互作用范围，不能只按最小原胞机械复制。
要点：位移幅度太大容易混入非谐响应，太小容易被数值噪声影响。
要点：虚频需要结合大小、位置、超胞收敛和振动模式一起判断。
要点：二维材料、slab 和吸附体系要特别说明真空层、活动原子范围和固定层边界。
