说明：本文华算科技主要介绍虚频的数学来源、q 点与结构周期的关系、沿虚频模式寻找低能构型的方法，以及温度、压力和非谐自由能对虚频归属的影响。
虚频的定义和数学来源
给定一个原子构型 R0，能量在小位移 u 附近的二阶项写成 E≈E0+uTHu/2。H 是 Hessian；晶体计算把实空间力常数做质量加权和傅里叶变换，得到每个 q 点的动力学矩阵 D(q)。对角化所得本征值 λqν 满足 ωqν2=λqν，本征矢 eqν 记录对应的集体位移。
λqνqν=i√|λqν|。声子图通常把 i√|λ| 的模长放到零线下方，形成“负频率”分支。晶格不存在沿负频率往复的周期振动；零线下的数值是负曲率的绘图记号，对应参考构型沿 eqν 位移后能量下降。
虚频模长反映参考点附近负曲率与参与原子质量的组合。|ω| 较大，对应谐展开中的下坡更陡；双阱深度、远处极小值和相变能垒由三阶、四阶及更高阶项决定，单个 i100 cm−1 数值并不包含势阱深度。THz、cm−1 和 meV 只是同一模能量的不同单位。
非线性分子在极小值处有 3N−6 个内部振动，平移和转动对应六个零模；线性分子有 3N−5 个内部振动。晶体在 Γ 点有三个声学零模。有限力误差、残余应力或平移不变性未完全满足时，零模附近会留下很小的负值。孤立的 −0.1 cm−1 与延伸到有限 q 区间的负分支属于不同的数值形态。
虚频出现何q点且为何改结构周期？
晶体本征模在相邻原胞间带有相位因子 exp(iq·R)。Γ 点 q=0 的位移在每个原胞同相，畸变可以保留原胞平移周期；q=(1/2,0,0) 的相邻原胞符号交替，沿第一晶格方向至少需要两倍周期。虚频所在 q 点直接给出畸变波长，M、R、N 等高对称 q 点的软模常对应倍周期倾转、剪切或层错图样。
有限位移超胞只容纳与其周期相容的一组 q 点，力常数截断又会改变长波支。光激发 Si 的收敛结果采用 12×12×12 动量网格与 8×8×8 有限位移超胞；4×4×4 网格和 3×3×3 超胞给出的低频 TA 曲线、负频范围和态密度均发生变化。负频分支应随超胞与 q 分辨率复算，零线附近的单点符号才有稳定含义。
负曲率还可能局限在很窄的倒空间区间。光激发载流子在 Si 中产生近 Γ 点 Kohn 异常；电子温度 Te=789 K 时，软化集中在有限 q，Γ 点声速仍为正。Te=1579 K 的占据展宽改变异常形状，曲线可在更宽区间进入零线以下。虚频宽度会随电子占据模型变化，它与离子温度属于两项独立条件。
N1×N2×N3 超胞允许的波矢分量由 mj/Nj 离散取值决定。若最深负频位于 1/4 附近，2×2×2 盒子没有对应模，沿该波长的重排会被周期条件排除。原胞声子谱与实际畸变超胞必须波矢相容，否则沿本征矢生成的结构会在超胞接缝处失配。
怎样沿虚频模式找到新的构型？
本征矢给出参考构型处的下坡方向，长度由归一化约定决定。质量加权模可写成 Rκ(Q)=Rκ,0+Qeκ/√Mκ；周期模还要把 q 与 −q 组合成实位移。e 与 −e 描述同一条坐标的两个方向，简并模则张成二维或三维子空间。本征矢是局部切向量，低能终态的位置需要有限振幅能量曲线和结构弛豫确定。
沿 Q 方向的局部势能可写成 U(Q)=−aQ2+bQ4。负二次项产生虚频，正四次项把能量限制在有限范围，两个极小值位于 ±√(a/2b)。光激发 Si 的模式模型把激发前谐势与激发后双阱并列，ΔQ 和 ΔU 分别记录畸变幅度与势阱深度，二者都不是虚频模长本身。
固定其他自由度扫描 E(Q) 可检验二阶负曲率是否持续到有限振幅；从 +Q 和 −Q 构型分别解除对称性并弛豫，可得到相同极小值、对称等价极小值或不同亚稳态。多个负模相互耦合时，最低路径会偏离任意单一 eqν，二维 E(Q1,Q2) 面可区分竞争畸变和协同畸变。
虚频模与真实动力学路径可用轨迹投影比较。光激发 Si 的 512 与 1728 原子轨迹中，虚频模只占全部本征模的 2.3%，约 700 fs 时却分别贡献 97% 和 87% 的归一化位移。轨迹位移在虚频子空间中的权重把谐近似的负曲率与后续原子重排联系起来。
分子反应中的一阶鞍点恰有一个负 Hessian 本征值，该虚频本征矢应沿成键、断键或原子迁移坐标。正负位移后的内禀反应坐标分别通向反应物与产物；两个以上内部虚频对应高阶鞍点或尚未收敛的约束构型。晶体软模的 q 周期信息在孤立分子 Hessian 中没有对应项。
温度自由能为何会改变虚频的解释？
常规谐声子是固定晶格常数、电子占据和 0 K 参考构型附近的二阶响应。热膨胀会改变体积，三声子与四声子相互作用会移动模频，量子零点涨落也会重塑有效势能面。有限温自洽声子得到的是温度相关的重整化频率 Ωqν(T)，它与 0 K 谐频率属于不同近似。
Fm3̅m-LuH3 的谐声子在多个 q 区域出现负值。SSCHA 纳入量子非谐晶格效应后，0 K 时大部分分支转为实频，X 与 K 附近仍保留负模；温度升到 150 K 以上后这些失稳受到抑制。同一结构的虚频会随 T 和 p 改变，图中 4.4 GPa、150 K 条件下的低频失稳已消失。
在压力、应变或电子占据改变势能面后，固定实验晶格常数与完全弛豫晶格常数可能给出不同的频率符号；金属费米面附近的电子展宽会移动 Kohn 异常；极性晶体的 Γ 点纵横光学模还受非解析项控制。结构、体积、磁序、电荷态与电子温度应和虚频数值写在同一结果记录中。
所有模转为实频只确定该构型在给定 T、p 下的局部动力学稳定状态。候选相是否出现由 Gibbs 自由能比较决定，振动熵、非谐修正、电子自由能和 pV 项都会进入相对自由能。Zr、ZrO2 与 GeTe 的计算把虚频区、频率重整化过渡区和全实频区同 ΔF(T) 放在一张图内。
BCC Zr 的重整化负模约在 700 K 后转为实频，实验 HCP→BCC 转变位于约 1135 K。只含准谐项的 ΔFQHA 在 1500 K 内没有交叉；加入非谐自由能后，ΔFQHA+anh 在约 1100 K 穿过零线。局部曲率转正与竞争相自由能换序相差约 400 K，两项温度回答的是不同物理量。
