弹性常数算出来以后，很多人会直接套 Born 条件：满足就说材料力学稳定，不满足就说结构不合理。这个方向大体没错，但很容易说过头。Born 判据检查的是平衡结构在小应变附近的二阶能量曲率是否为正。它回答的是“轻微拉伸、剪切或压缩时，结构会不会沿这个应变方向自发塌下去”。Born 条件是力学稳定性的必要检查，不是材料一定能合成、一定不会相变、一定适合催化环境的保证。
因此，力学稳定性判据要和声子谱、形成能、相竞争、AIMD、表面重构和化学环境分开读。若只把 Cij 满足条件写成材料合理，就会把小应变弹性问题扩展成了完整材料可行性问题。更准确的写法，是说明该结构在当前晶格、当前压力、当前交换关联泛函和当前弹性矩阵下通过了小应变力学检查。
它检查的是弹性能二阶项
晶体受到很小应变时，总能量可以在平衡结构附近展开。若二阶弹性矩阵在允许的应变空间中为正，说明任意小应变都会让能量升高，结构在力学意义上可以抵抗这些微小扰动。Born 判据看的不是结构会不会在真实环境中长期存在，而是小应变附近能量曲率是不是正。
ΔE = 1/2 V0 Σij Cij εi εj
εT C ε > 0
第二个式子表示弹性矩阵对任意非零小应变 ε 都应给出正的弹性能。不同晶系因为对称性不同，Cij 的独立分量不同，Born 条件也不同。立方晶系常见条件为 C11 – C12 > 0、C11 + 2C12 > 0、C44 > 0；六方、四方、正交晶系不能直接照搬立方晶系的三条条件。
图1 KMgX3 结构模型示意，弹性判据必须对应明确晶体结构和对称性。DOI：10.1039/d5ra04332h
判据要跟晶系和维度匹配
最常见的错误，是把立方晶系的 Born 条件套到非立方结构上，或者把三维体相的弹性常数逻辑直接套到二维材料与 slab 模型上。二维材料常需要使用二维弹性常数或对真空厚度做归一处理；表面 slab 由于含真空层，直接把 OUTCAR 中三维应力当作体相弹性常数也会产生口径问题。Born 条件不是一套通用口诀，必须绑定晶系、维度和弹性常数定义。
另外，弹性常数最好在接近平衡的结构上计算。若结构残余应力很大、原子力未降到合理范围，或者晶格参数还没有充分优化，得到的 Cij 会混入未释放的内应力。此时即使条件看起来满足，也只能说明这个受应力结构在当前设置下的曲率表现。
力学稳定不等于热力学相稳定
Born 条件不比较不同相之间谁能量更低。一个结构可以在小应变下不塌陷，但相对于另一种晶型、分解产物或元素相仍然能量偏高。对新材料预测，形成能、凸包能量、相图和竞争相才回答“这个相有没有热力学保留的可能”。对催化表面，表面能、缺陷形成能和吸附诱导重构也不由 Born 条件直接给出。
图2 CeCrO3 能量-体积曲线示意，说明相稳定与平衡体积需要独立考察。DOI：10.1021/acsomega.3c09292
例如一个高压相在某个压力下满足弹性条件，降到常压后却可能不再是最低能结构；一个二维材料在弹性矩阵上表现正常，也可能因边缘重构、层间反应或环境氧化而不适合实际使用。写结论时应把“力学上可承受小应变”和“材料在实验条件下可获得”分开。
力学稳定也不等于动力学没有软模
声子谱检查的是原子集体振动模式。Born 条件主要在长波极限和均匀应变附近发挥作用，而声子谱可以看到有限波矢处的软模。若声子谱存在明显虚频，说明结构可能沿某个原子位移模式降低能量，即使弹性条件满足，也需要重新检查结构、超胞、收敛和可能的相变路径。
图3 KMgX3 声子色散示意，可用于区分弹性小应变检查与动力学振动检查。DOI：10.1039/d5ra04332h
这对催化材料尤其重要。表面、缺陷和吸附态常常不是完美体相，局部原子位移可能比整体应变更敏感。一个体相弹性常数满足条件，并不说明带缺陷表面或吸附后表面不会发生重构。
应变幅度太大或太小都可能出问题
弹性常数来自能量-应变或应力-应变关系。应变太大时会进入非线性区，二阶弹性近似不再干净；应变太小时，数值噪声可能淹没能量差。常见做法是取几组正负小应变，检查能量或应力随应变是否接近二次或线性关系，而不是只跑一个应变点。
σi = Σj Cij εj
若用应力法，EDIFF、ENCUT、k 点、PREC、LREAL、ISIF 和 Pulay stress 都会影响结果。金属体系还要注意 ISMEAR 和 SIGMA 不宜让应力波动过大。若材料含磁性元素，弹性常数还应建立在合适磁性态上；把铁磁和反铁磁混在不同应变点中，会让拟合结果失真。
图4 KMgX3 的能带与力学性质图组示意，说明电子结构与弹性输出常需分层阅读。DOI：10.1039/d5ra04332h
压力条件要说清楚
Born 条件有时是在零压下写，有时要考虑有限压力下的修正弹性常数。若研究对象是高压相、外加应变材料或界面受限结构，就不能只套常压表达式。对 slab 和二维材料，也要说明真空方向是否参与弹性矩阵，避免把真空层当作材料的一部分来解释。
SYSTEM = elastic_static
ENCUT = 520
EDIFF = 1E-7
IBRION = -1
NSW = 0
ISMEAR = 0
SIGMA = 0.05
LREAL = .FALSE.
PREC = Accurate
# 对每个小应变结构分别做同口径静态计算
这段片段强调的是同口径静态计算。弹性常数不是靠某一行参数自动得出，而是靠一组小应变结构在相同计算精度下建立能量或应力关系。对于离子可弛豫弹性常数，还要明确应变后是否允许内部坐标弛豫。
把稳定性拆成几类回答
比较清楚的写法是：弹性常数满足对应晶系 Born 条件，说明结构在小应变力学扰动下没有直接失稳；声子谱无明显虚频，支持动力学层面的可保留性；AIMD 在目标温度下结构未发生明显破坏，说明有限温度短时尺度下未观察到剧烈重构；形成能或相图再用于讨论热力学竞争。
图5 KMgX3 AIMD 结果示意，可作为有限温度结构保持情况的补充检查。DOI：10.1039/d5ra04332h
若某一项不支持，不一定代表计算全错，而是说明材料稳定性只在某个层面成立。例如弹性通过但声子有软模，可写成该结构在均匀小应变下满足力学条件，但仍可能沿特定原子位移模式转变；AIMD 保持但相图能量偏高，可写成短时动力学可保留，但热力学竞争仍需谨慎。
用于催化研究时，体相力学稳定只能作为材料背景。真正的反应模型还要看表面能、吸附诱导重构、活性位迁移、溶剂和电位条件。Born 条件可以帮你排除明显力学不合理的体相结构，但不能替代催化工作态检查。
要点一：Born 条件检查的是小应变附近弹性能二阶项是否为正，不能扩展成所有稳定性都成立。
要点二：不同晶系、二维材料和 slab 模型应使用匹配的弹性常数定义与判据。
要点三：弹性通过不代表相稳定、动力学无软模、有限温度不重构或催化环境中不变化。
要点四：弹性常数对收敛精度、应变幅度、残余应力、磁性态和内部坐标处理都很敏感。
要点五：更严谨的结论应把弹性、声子、AIMD、相竞争和表面工作态分开说明。
