晶格振动说的是晶体里原子围绕平衡位置做集体运动；声子则是这种集体振动在量子化描述中的能量单位和模式标签。可以先粗略理解为：晶格振动是现象，声子是描述这种振动的语言。声子不是额外多出来的粒子，而是晶格集体振动的量子化表达。
在 VASP 和 Phonopy 里，我们常看到声子谱、声子态密度、虚频、Γ 点频率这些词。它们背后都在回答同一个问题：原子如果轻微离开平衡位置，体系会以哪些方式振动，这些振动频率是多少，结构会不会沿某个方向自发变形。
图1 声子性质计算流程示意，展示从力常数到声子谱和热力学量的连接。DOI: 10.1038/sdata.2018.65。
晶格振动先描述什么？
原子是在随机乱动吗？
晶体中的原子并不是固定不动的点。即使在低温下，它们也会围绕平衡位置振动；温度升高时，振动幅度和可激发模式会改变。由于原子之间有化学键和力常数，一个原子动起来会牵动邻近原子，因此振动往往是集体性的。
如果只看一个原子上下摆动，很难描述整个晶体。更自然的做法是把许多原子的耦合运动分解成一组独立的正则模。每个正则模都有特定频率、相位关系和原子位移方向。晶格振动的核心，是原子之间通过力常数耦合形成的集体运动。
为什么要用谐振近似？
声子计算通常从平衡结构附近的小位移出发，把势能面近似成二次曲面。这样就能用力常数矩阵描述原子位移和回复力的关系，再求出振动频率和模式。这个近似叫谐振近似，是多数基础声子谱计算的起点。
Fi = −Σj Φij uj 式(1)
式(1)里 Fi 是第 i 个自由度上的力，uj 是第 j 个自由度的位移，Φij 是力常数。力常数越大，说明原子被拉回平衡位置的趋势越强，对应频率通常越高。
声子怎样对应一组振动模式？
一个 q 点代表什么？
在周期晶体中，振动可以用波矢 q 来标记。q 描述振动在不同晶胞之间的相位变化：Γ 点表示所有晶胞同相振动；非 Γ 点则表示相邻晶胞之间存在相位差。声子谱沿高对称路径画出来，就是不同 q 点下振动频率的变化。
每个 q 点上会有多条声子分支。若原胞里有 n 个原子，就有 3n 条分支，其中 3 条是声学支，其余是光学支。声学支对应长波极限下整个晶格的平移或低频集体运动，光学支常对应原胞内部不同原子之间的相对振动。声子谱上的每一条线，本质上是一类集体振动模式随 q 的频率变化。
图2 不同 SiO₂ 相的声子能带和声子态密度对比，说明声子谱能反映结构振动特征。DOI: 10.1038/sdata.2018.65。
声子态密度又在看什么？
声子谱告诉我们频率随 q 点怎样变化，声子态密度则把所有 q 点和所有分支的频率分布统计在一起。它不再强调某条路径上的色散细节，而是回答哪些频率范围内可用振动模式更多。
这对热力学量很有用。热容、振动熵、零点能都和频率分布有关。低频模式多，熵贡献往往更明显；高频键伸缩模式则会显著影响零点能。声子态密度把分散在整个布里渊区的振动模式汇总成频率分布。
声子和普通分子频率有什么不同？
为什么晶体要看整条路径？
分子是有限体系，振动频率是一组离散正则模；晶体是周期体系，振动还要考虑不同晶胞之间的相位关系，所以需要 q 点。只算 Γ 点频率，相当于只看所有晶胞同相振动的一小部分模式，不能代表完整晶体的动力学稳定性。
这也是为什么材料声子谱常沿 Γ-X-M、Γ-K-M 等高对称路径展示。路径上的频率变化能告诉我们某些长波或短波扰动是否会让结构不稳定。对二维材料、表面 slab 和软晶格材料，低频区域尤其需要谨慎看。
图3 单层石墨烯和 MoS₂ 的声子色散关系示例，展示二维材料中声学支和光学支的差异。DOI: 10.3762/bjnano.8.136。
声子是不是只能用于固体？
声子这个概念主要用于周期固体、二维材料、表面和有序结构。吸附物频率计算常更接近分子振动分析，但当吸附物和表面强耦合时，也会出现吸附物振动与表面低频模式混合的情况。
因此，催化计算中要区分两类任务：一类是材料声子谱，用来判断结构动力学稳定性；另一类是吸附态频率，用来估算 ZPE 和振动熵。两者都和原子振动有关，但计算对象和解释口径不同。
虚频说明结构一定不稳定吗？
虚频来自什么？
在声子谱中，如果某个模式频率写成负值或虚数，意味着该方向上的二阶曲率为负。直观地说，结构沿这个振动方向轻微位移后，能量不是被拉回去，而是继续降低。对体相晶体来说，明显虚频通常提示当前结构不是该条件下的动力学稳定相。
ω² ∝ λ(D) 式(2)
式(2)中 D 是动力学矩阵，λ 是本征值。若 λ 为负，就会表现为虚频。需要注意，声子软件常把虚频画到负频率区域，方便可视化；这不是频率真的小于零，而是负曲率的标记。虚频本质上提示某个振动模式方向存在负曲率。
为什么小虚频要结合体系判断？
二维材料的 ZA 声学支、表面 slab 的低频摆动、超胞不够大或数值精度不足，都可能带来靠近 Γ 点的小虚频。此时不能只看一眼负值就下结论，而要检查位移模式、超胞尺寸、力收敛、真空层和截断能。
如果虚频模式对应整体波纹、表面层轻微滑动或有限尺寸效应，可能需要更大超胞和更高精度；如果虚频对应清晰的键断裂、相变位移或原子成排偏移，就更可能是真实不稳定。虚频判断必须回到振动模式本身。
图4 二维异质界面中的声子色散和振动模式示例，说明局域结构畸变会影响声子模式。DOI: 10.1126/sciadv.abi6699。
VASP和Phonopy怎样连接？
有限位移在做什么？
有限位移法会在超胞里轻轻移动某些原子，然后用 VASP 计算受力。Phonopy 根据这些力反推出力常数矩阵，再构建动力学矩阵，最后得到声子频率和模式。位移幅度太小会受数值噪声影响，太大则可能超出谐振近似。
超胞大小也很重要。力常数有空间范围，超胞太小会把远程相互作用截断得不合理，尤其对极性材料、二维材料和软晶格体系影响更明显。声子谱质量首先取决于结构、力常数和超胞是否可靠。
SYSTEM = phonon_force_displacement
ENCUT = 520
EDIFF = 1E-8
IBRION = -1
NSW = 0
ISMEAR = 0
SIGMA = 0.05
LREAL = .FALSE.
ADDGRID = .TRUE.
# 参考示例：Phonopy 生成位移结构，VASP 只负责高精度单点力计算
怎样把声子和晶格振动说清？
更自然的表达是：晶格振动是原子集体运动，声子是这些集体运动模式的量子化描述；声子谱展示不同 q 点的振动频率，声子态密度统计所有模式的频率分布；虚频提示某个振动方向可能引发结构变形。
如果研究催化表面，还要说明计算的是体相、二维材料、干净 slab 还是吸附态。不同对象对应不同解释：体相声子看相稳定，slab 声子看表面低频和约束边界，吸附态频率更多服务于 ZPE 和熵修正。
声子计算前的结构优化精度要比普通几何优化更严格。若残余力偏大，力常数会把未收敛结构的斜率带进去，低频区域尤其容易出问题。常见做法是把 EDIFF 收紧，把最终力降到较小范围，并关闭会引入不稳定噪声的近似设置。
图5 石英低能声子色散与模式性质示例，展示声子分支可携带更丰富的振动信息。DOI: 10.1038/s41586-023-06016-5。
对极性材料，还要留意 LO-TO 分裂和非解析项修正；对二维材料，要处理真空层、偶极和面外声学支的数值行为；对表面 slab，要说明哪些层参与振动、哪些层被固定。声子谱不是一个独立按钮，而是和模型边界强相关的振动分析。
在催化自由能中使用频率时，也要避免把材料声子谱和吸附物频率混用。*H、*OH、*CO 这类吸附态通常只需要关注反应区振动对 ZPE 和熵的贡献；判断材料动力学稳定性时，则要看完整周期结构的声子谱是否存在真实虚频。
