声子谱（phonon dispersion）描述晶格振动频率随波矢的色散关系，是连接晶格动力学与宏观热物性的桥梁。从声子谱可获得声速、德拜温度、热容、热导率（结合弛豫时间）以及结构相变的软模信息，是材料热学性质与稳定性分析的基础。实验上声子谱可由非弹性中子/X 射线散射测量，而计算模拟提供了互补的原子级预测手段。
LAMMPS 的 fix phonon 命令可在经典 MD 框架下通过分析速度自相关函数（VACF）或格林函数方法获得声子谱。本项目以面心立方镍（Ni）为对象，用 LAMMPS（EAM 势）计算其声子色散关系与声子态密度，并与实验中子散射数据对比。
对晶体中传播的平面波振动 u = u₀exp[i(k·r − ωt)]，声子频率 ω 与波矢 k 的色散关系 ω(k) 由动力矩阵本征方程确定。FCC 晶体每个原胞含 1 个原子、3 个自由度，因此有 3 条声学支（1 纵 + 2 横），无光学支。在 Brillouin 区边界（X、W、L 点），横声学支可发生简并。Ni 的声子谱最高频率约 9 THz。
LAMMPS fix phonon 基于线性响应/格林函数方法：施加小幅外力扰动（如单原子位移），监测体系对扰动的响应，从位移-力关联中提取动力矩阵，进而对角化获得各波矢的声子频率。另一种常用方法是速度自相关函数的傅里叶变换获得声子态密度 g(ω)。fix phonon 需要较长的平衡轨迹以保证统计精度。
声子态密度 g(ω) 描述频率分布，其低频段呈德拜二次方行为（g(ω) ∝ ω²）。由 g(ω) 可计算晶格热容（Debye 模型修正）、德拜温度与零点能。声子谱中的虚频（ω² < 0）是结构失稳的标志，可用于判断相稳定性。
图 11 的左面板给出 FCC Ni 沿 Γ–X–W–L–Γ 高对称路径的声子色散关系。图中 3 条声学支从 Γ 点（k=0）零频率出发：纵声学支（LA，红色）沿 Γ→X 快速上升至约 7.5 THz，横声学支（TA，蓝色/绿色）上升较缓至约 4.2 THz；在 W、L 点 TA 支发生部分简并。由于 Ni 原胞含 1 个原子，无光学支，这与金刚石（2 原子原胞，有光学支）形成对比。所有频率均为正（无虚频），证明 EAM 势描述的 FCC Ni 结构动力学稳定。最大频率约 9 THz，与实验非弹性中子散射结果（约 8.7–9.2 THz）一致，验证了势函数的可靠性。声学支在 Γ 点附近的线性斜率给出声速（纵波约 4900 m/s、横波约 2900 m/s），与 Ni 的力学实验值同量级。
图 11 的右面板给出声子态密度 g(ω)（棕色曲线）。g(ω) 在低频段近似抛物线上升（德拜二次方行为，红色阴影示意），对应声学支的连续分布；在 4–6 THz 出现第一个宽峰（横声学支的态密度极大），7–9 THz 出现第二个峰（纵声学支的范霍夫奇点特征），9 THz 处态密度截止。态密度的峰结构直接对应色散曲线的平缓区域（高密度态）。由 g(ω) 可积分获得 Ni 的德拜温度（约 375–450 K）与 300 K 晶格热容（接近经典值 3R），为热物性预测提供了基础。
两面板结合，左面板给出声子频率随波矢的色散（用于声速、稳定性与中子散射对比），右面板给出频率分布密度（用于热物性与热容计算），共同完成 Ni 晶格动力学的系统表征，验证 EAM 势并支撑热学性质预测。
本项目用 LAMMPS 计算了 FCC Ni 的声子谱：3 条声学支色散正常、最大频率约 9 THz 与实验一致，声子态密度呈德拜二次方加范霍夫峰结构。案例图以色散曲线与态密度双面板呈现。
后续可拓展：(1) 结合玻尔兹曼输运方程与三阶力常数计算晶格热导率；(2) 扫描温度研究非谐效应对声子谱的软化；(3) 对 Ni 基合金（Ni-Al、Ni-Ti）计算声子谱，研究合金化对热导率与相稳定性的影响。
