NEPHONON：基于神经进化势的高效声子计算程序

《Computer Physics Communications》：NEPHONON: An Efficient Phonon Calculator Based on Neuroevolution Potentials

　　

声子计算对于理解材料的热力学与动力学性质至关重要，然而密度泛函理论(DFT)方法在摩尔超晶格等大尺度体系上常面临难以承受的计算代价。研究人员提出NEPHONON，一款开源高性能程序包，旨在以近DFT精度快速计算扩展体系的声子性质。通过利用机器学习神经进化势(Neuroevolution Potentials, NEPs)的计算效率，NEPHONON即便在含数千原子的体系中也能快速生成二阶力常数。除标准声子能带结构、态密度与群速度外，NEPHONON实现了等频声子面与非弹性中子散射谱S(Q,E)的模拟。该程序还提供完整声子本征矢访问，支持拓扑手性声子分析与振动模式可视化。通过在扭曲双层黑磷、石墨烯与铜上的基准测试，研究人员展示了该程序的性能，证实NEPHONON是快速探索材料中复杂声子现象的有力工具。

研究背景方面，晶格振动量子化为声子，主导比热、热膨胀、晶格热导、相稳定性与超导等材料性质。实验手段如非弹性中子散射(INS)、高分辨非弹性X射线散射、电子能量损失谱与拉曼光谱可直接获取声子色散，但解释常依赖理论计算。传统密度泛函微扰理论(DFPT)与有限位移法虽为标准第一性原理工具，但DFT随体系尺寸扩展代价剧增，常规精确计算仅限数十原子单元，严重制约摩尔超晶格、无序合金、复杂化合物等千原子级体系研究。经典力场虽廉价但精度与可迁移性不足。机器学习原子间势(MLIPs)特别是神经进化势(NEPs)以DFT数据集训练神经网络，兼顾量子精度与经典效率，为大体系声子计算提供新路径。为此研究人员发布NEPHONON，采用Modern Fortran写晶格动力学核心、C++后端算力的混合架构，对接BLAS/LAPACK做动力学矩阵装配与对角化，以预训练NEP算力，填补大尺度近DFT精度声子计算工具空白。该文发表于《Computer Physics Communications》。

关键技术方法上，研究人员采用有限位移法，对超胞内每个原子沿三笛卡尔轴正负向位移，共6n次位移，由CPU版NEP计算器基于预训练nep.txt输出原子力，中心差分近似求二阶力常数Φ_αβ(iκ,jκ')；力常数可选迭代施加平移对称（声学求和规则）与置换对称。动力学矩阵D_κκ'^αβ(q)由力常数对晶格矢R傅里叶求和并除质量平方根构造，用LAPACK对角化得频率ω_qν与本征矢e_qν；群速度由Hellmann-Feynman定理解析求导动力学矩阵获得；长程偶极作用经非解析修正(NAC)用Born有效电荷与高频介电张量补正；INS谱S(Q,ω)按相干单声子动态结构因子公式算单晶高对称路径与粉末球平均（黄金比序列采样），δ函数高斯展宽拟实验分辨率。基准体系含90°扭曲双层黑磷3×3×1超胞2520原子、单层石墨烯、面心立方铜，NEP模型取自官网预训练权重。

研究结果部分，第二节力常数中，研究人员阐明谐波近似下Φ_αβ定义与中心差分公式，指出NEP算力廉价故不做对称性约简即可跑满6n位移，并施迭代对称约束满足声学求和与置换对称。第三节声子色散与本征分析中，给出D(q)构造与对角化方程，强调本征矢可用于手性/拓扑声子分析；群速度表达式避免能带交叉处有限差分失效；NAC项在q→0极限公式明确含单位胞体积Ω、Born张量Zⁱ、ε^∞；INS模块输出单晶S(Q,ω)与粉末S(|Q|,ω)，含Bose因子n(ω)、相干散射长度b_κ、Debye-Waller因子W_κ(Q)。第四节安装与工作流中，程序依赖gfortran/g++与BLAS/LAPACK，Makefile编译产单可执行文件；输入inp.control定结构与参数，放nep.txt算力常数，或读PHONOPY格式力常数，inp.lotosplitting给NAC参数，标志位控流程。第五节应用实例中，5.1扭曲双层黑磷用2520原子超胞单核少于1.5小时跑完，能带显示光学支扁平化、30–42 meV大声子带隙、Γ与M点低能模，态密度印证上述特征；5.2石墨烯等频面110 meV与50 meV以FermiSurfer格式输出，闭环/开放曲线反映群速度各向异性；5.3铜的INS谱单核少于2秒算出，与J-PARC实验纵向/横向声学支及强度变化吻合。第六节总结中研究人员指出NEPHONON以Fortran核心+C++轻接口调NEP经有限差分产生二阶力常数（或读Phonopy格式），LAPACK/BLAS对角化，覆盖能带、态密度高斯展宽、解析群速、等频面、单晶S(Q,E)与粉末S(|Q|,E)、可选Γ附近LO-TO非解析修正，namelist输入与模块化设计利于扩展，适千原子级复杂材料高通量声子性质计算。

讨论与结论翻译：研究人员总结NEPHONON为一款开源可移植声子计算程序，含Fortran核心与轻量C++接口调用NEP通过有限差分生成二阶力常数（或读取PHONOPY格式FORCE_CONSTANTS），使用LAPACK/BLACK进行稳健对角化。工作流涵盖用户定义路径声子能带、高斯展宽态密度、解析群速度、声子等频面、单晶S(Q,E)与粉末S(|Q|,E)中子散射强度，及利用介电张量与Born有效电荷在Γ附近做非解析LO-TO修正。强调可复现与可移植，采用简单Fortran namelist输入、清晰文件I/O、兼容常见后处理工具输出，模块化便于新力常数源、对称约束与分析输出接入。该方法适用于复杂材料大超胞谐波声子研究（乃至数千原子），可作为跨材料数据库计算声子相关性质的高通量工具。