《Computational Condensed Matter》：Critical importance of k-point convergence in supercell calculations of mechanical instabilities: Implications for shape-memory alloys
编辑推荐：
在超胞计算中，k点采样收敛性对机械稳定性预测至关重要。研究人员针对形状记忆合金的力学不稳定性研究，通过高NbRu合金案例，揭示了不充分k点采样会导致定性错误的预测，并提出并验证了保持倒空间采样密度恒定的普适性标度关系。这项工作为涉及大规模模拟的计算材料设计，特别是机械不稳定性决定功能特性的系统，提供了关键的方法学指导。
在计算材料科学中，通过第一性原理计算模拟材料的性质与相变已成为常规工具。然而，当研究涉及机械不稳定性——例如驱动形状记忆合金发生马氏体相变的关键软模时，计算参数的精确收敛尤为重要。尽管对原胞进行k点采样收敛性测试是标准操作，但在使用超胞（研究缺陷、化学无序和长波声子模所必需）时，情况变得复杂。为了节省计算时间，在基于超胞的ab initio分子动力学(AIMD)模拟中，仅使用Γ点（即倒空间原点）采样已成为常见做法。这种便利做法是否存在风险？当研究的核心是那些由微小能量差决定稳定性的力学不稳定性时，不充分的k点采样是否会得出完全错误的物理结论？这直接关系到高通量计算筛选和材料设计的可靠性。
为了回答这些问题，研究人员以高温形状记忆合金(high-temperature shape-memory alloy, HTSMA) NbRu为具体案例，系统研究了k点采样密度在超胞计算中对机械不稳定性预测的定性及定量影响。他们发现，对于大型超胞（如128原子甚至432原子），仅使用Γ点采样会完全颠倒对驱动马氏体转变的关键不稳定性波的稳定性预测，错误地得出材料稳定的结论。相反，遵循其提出的标度关系进行充分收敛的计算则能正确揭示不稳定性。这项研究强调了在涉及力学不稳定性的计算中，必须将k点采样收敛性测试的重点放在相关的不稳定模能量上，而非仅关注总能。相关研究成果发表在《Computational Condensed Matter》上。
研究人员开展此项研究主要应用了第一性原理计算、超胞构建与标度分析、以及针对特定声子模的“冻结声子”能量计算等关键技术方法。所有密度泛函理论(density functional theory, DFT)计算均使用Vienna Ab initio Simulation Package (VASP)软件完成。研究以B2结构的NbRu为模型体系，重点关注其<110>[1?0]不稳定波（M点声子模）。通过构建不同尺寸的超胞（16原子、128原子、432原子），并系统改变k点网格密度，计算了在特定位移波扰动下的能量变化ΔE，以此评估体系的稳定性。研究中所用的计算参数（除k点网格外）与参考文献保持一致，确保了对比研究的基础一致性。
研究结果
- •
2.1. Brillouin zone folding and k-point density
研究从基本原理出发，阐述了超胞计算中k点收敛的标度关系。当超胞包含N个原胞时，其布里渊区(Brillouin zone, BZ)体积缩小为原胞的1/N。为了在倒空间保持等效的采样密度，超胞计算所需的总k点数应与原胞k点数除以N成正比（Nksupercell= Nkprimitive/ N）。对于各向同性的超胞，沿每个倒易方向的k点网格数应按N1/3的倒数进行标度。
- •
3.3. Convergence in small supercells
在16原子超胞中，对不同k点网格（从2×2×2到16×16×16）下<110>[1?0]波的ΔE进行计算。所有网格均显示系统在位移下失稳（ΔE < 0），表明定性行为一致。然而，能量差的幅度变化显著，定量收敛在8×8×8 k点时达到（与更密网格的差异远低于1 meV/原子的标准DFT精度极限）。
- •
3.4. Qualitative errors in large supercells
在128原子超胞中，仅使用Γ点采样会错误地预测系统在位移下稳定（ΔE > 0），完全颠倒了物理事实。即使对于432原子超胞，Γ点采样也仅显示微弱稳定。而按照标度关系进行正确收敛的计算（如128原子用4×4×4网格，432原子用2×2×2网格）则一致显示失稳。对原子受力的进一步分析表明，Γ点采样给出近乎线性的恢复力响应，暗示模式稳定；而收敛的k点网格则揭示出非线性力响应，这是由电子-声子耦合导致的声子软化的特征，正确地反映了不稳定性。
- •
4. Practical guidelines
基于研究发现，文章为研究接近机械或动力学不稳定性的超胞计算提出了实用的k点收敛指导原则：1. 确定最重要的不稳定声子模，并构建可容纳该波的最小超胞；2. 在该最小超胞中，收敛该不稳定模的能量所需的k点密度；3. 应用上述普适标度关系来确定更大超胞所需的最小k点网格；4. 对关键计算进行明确的收敛性测试。
研究结论与意义
本研究表明，在超胞计算中，正确的k点收敛对于研究机械不稳定性至关重要，它不仅影响定量精度，更关系到定性预测的正确性。文章提出并验证了一个简单的普适标度关系：为保持等效的倒空间采样密度，超胞计算所需的k点数应与原胞k点数除以超胞包含的原胞数N成正比。这一关系必须应用于晶体中最重要的不稳定位移波的能量景观，而非仅看总能。
以高温形状记忆合金NbRu为具体案例的验证表明，不充分的k点采样（如仅用Γ点）可以在大型超胞计算中导致完全错误的稳定性预测，从而可能在高通量计算筛选中错误地排除有前景的材料候选物（如NbRu），或错误地选中不稳定的材料。这些发现对快速发展的计算材料设计领域，特别是涉及大规模、省时模拟的研究具有直接意义。研究强调，对于金属体系，由于费米面效应和电子结构细节主导结构稳定性，必须确保足够的k点采样密度以捕捉电子-声子耦合。在基于机器学习(ML)的材料设计和高通量数据库中，若训练数据因k点采样不足而包含系统性误差，所得的模型无论架构如何先进，都将继承这些缺陷。因此，在涉及力学不稳定性的材料计算中，遵循本文的标度关系进行充分的k点收敛测试，是避免在材料信息学流程中传播定性错误的重要投资。
