弹性张量输出里出现负值时，最容易发生两种极端判断：一种是直接认为材料力学不稳定，另一种是把负值全部当成计算错误。VASP 给出的 Cᵢⱼ 来自微小应变下的应力或能量响应，负值可能对应真实的软化方向，也可能来自结构、应变、对称性或拟合口径。因此它不是一个单独结论，而是需要定位来源的计算现象。
独立常数和派生模量一样吗？
弹性常数矩阵 C 描述应力 σ 与应变 ε 的线性关系。对不同晶系，矩阵中允许出现的独立项数量和稳定性条件不同。某些非对角项为负并不必然代表力学失稳，因为它们描述方向耦合；而某些剪切项或稳定性组合为负，则可能对应小形变下能量曲率为负。
σᵢ = Σⱼ Cᵢⱼεⱼ 式(1)
图1 展示高通量弹性常数计算会经历结构优化、施加应变、张量过滤等步骤。负值出现在弹性张量输出中时，应该沿这个流程回看：是原始结构已经不在局部极小值，还是应变后的应力不在线性区间，或者是张量过滤和对称化环节改变了独立常数。随后再确认负值属于 Cᵢⱼ 元素、稳定性组合，还是由 Cᵢⱼ 换算出的模量。
图1 弹性常数高通量计算与筛选流程。图中从结构优化、施加应变到过滤弹性张量，说明负值排查应沿计算流程逐步定位。DOI: 10.1038/sdata.2015.9。
方向响应为什么会有正负？
图2 把杨氏模量、线性压缩率、剪切模量和泊松比分别画成方向相关性质。弹性张量不是一个标量，它描述不同方向和剪切方式下的响应。某些材料存在负泊松比或负线性压缩率等特殊方向响应，但这类派生量和 Cᵢⱼ 负值不是同一个概念。
因此，看到负值后要先查单位和物理对象。例如 C₁₂ 为负，在某些晶系中可能表示横向耦合方向相反；C₄₄ 为负，通常更接近剪切失稳信号；Hill 剪切模量为负，则说明张量整体平均已经不合理。不同位置的负值有不同含义，不能把所有负号解释成同一种失稳。
图2 晶体方向相关弹性性质示意，包括杨氏模量、线性压缩率、剪切模量和泊松比。该图说明弹性结果需要区分张量元素和方向派生量。DOI: 10.1039/c9sc01682a。
残余应力从哪里来？
弹性常数计算要求基态结构处在足够接近局部能量极小的位置。如果晶格常数、原子坐标或磁矩还没有稳定，施加很小应变后得到的应力就会叠加残余应力，拟合出的斜率可能偏离真实线性响应。对低对称结构、层状材料、含空位或掺杂超胞，这个问题更常见。
例如，结构优化阶段只放松了离子位置，没有充分优化晶胞形状，随后直接计算弹性张量。此时某些剪切项为负，可能不是材料真实剪切失稳，而是初始晶胞带有未释放的形变。再如，磁性氧化物在不同 MAGMOM 初值下收敛到不同局域磁矩，弹性计算继承了较高能磁态。弹性负值排查要先确认结构、应力和磁态已经进入同一基态假设。
优化精度怎样设置？
结构优化通常需要更严格的力和应力收敛。可从 EDIFF = 10⁻⁶ eV、力收敛阈值约 -0.01 到 -0.02 eV/Å 量级开始测试，是否采用 ISIF = 3 或分步优化要按体系决定。对于层状或二维体系，还要处理真空方向与面内应变的边界，不能把非物理真空响应混入三维弹性常数。
若结构含缺陷、吸附物或低对称界面，建议先检查未应变结构的压力张量和关键键长，再施加正负应变。若正应变和负应变给出的能量曲线明显不对称，说明当前构型可能不适合直接做线性弹性拟合。
应变太小或太大会怎样？
应变幅度太小，应力差接近数值噪声；应变幅度太大，材料响应进入非线性区间。两种情况都可能让拟合斜率异常，尤其是弱相互作用层状体系、柔性框架材料或软模接近零的结构。合理做法是测试多个应变幅度，例如 ±0.5%、±1.0%、±1.5%，检查 Cᵢⱼ 是否在一个区间内稳定。
图3 展示大规模材料中体模量、剪切模量、泊松比和体积的分布。真实材料的弹性响应横跨很宽范围，软材料的斜率本来就小，因此更容易受应力精度影响。越软的体系越需要检查应变幅度、电子收敛和应力张量噪声。
图3 体模量、剪切模量、泊松比和原子体积分布。图中说明材料弹性响应跨度很大，软材料的应力拟合更容易受数值噪声影响。DOI: 10.1038/sdata.2015.9。
对称性会不会改写张量？
VASP 计算中若保留对称性，结构的微小畸变可能被强制映射回较高对称；若关闭对称性，低对称噪声又可能让等价张量项不再相等。对弹性常数而言，晶系判断、独立常数识别和张量对称化都很重要。错误地把低对称结构按高对称晶系套公式，会直接导致 Born 判据和模量换算出错。
以掺杂超胞为例，母相可能是立方晶系，但单个掺杂原子已经降低局域对称性。若仍强行使用立方晶系的 C₁₁、C₁₂、C₄₄ 三项判据，负值或不满足条件的来源就会被误判。弹性张量应按实际计算结构的对称性解释，而不是按理想母相套用。
负值和稳定性条件等价吗？
Born 稳定性要求弹性能量二次型在允许应变下为正。对立方晶系，常见条件包括 C₁₁ – C₁₂ > 0、C₁₁ + 2C₁₂ > 0、C₄₄ > 0；对六方、四方、正交等晶系，条件不同。可见，判断稳定性看的并不总是单个元素是否为正，而是张量组合是否让二次型保持正定。
ΔE/V₀ = 1/2 ΣᵢΣⱼ Cᵢⱼεᵢεⱼ 式(2)
式(2)说明，弹性稳定性本质上来自应变能的曲率。若某个稳定性组合为负，小形变会沿某个模式降低能量，材料在该结构附近不具备线弹性稳定性。若只是某个非对角耦合项为负，则要结合完整矩阵本征值和对应晶系条件判断。
各向异性会放大问题吗？
图4 展示弹性各向异性的分布。高度各向异性的材料中，某些方向非常硬，另一些方向接近软化，少量应力噪声就可能让方向模量或剪切响应出现异常。二维材料、层状氧化物、分子晶体和含空隙框架尤其要把方向响应分开看。
一个有用的检查是计算弹性矩阵本征值或用 ELATE、pymatgen 等工具查看方向模量。如果某个方向的剪切模量接近零，再出现轻微负值，就要进一步结合声子谱和结构重优化确认。Born 判据给出小形变稳定性边界，不能单独替代声子和热稳定性检查。
图4 无机化合物弹性各向异性分布及其与泊松比的关系。高各向异性体系中，方向相关软化更容易影响负值判断。DOI: 10.1039/c9sc01682a。
什么时候是真实失稳？
若结构已经充分弛豫，应变幅度测试稳定，对称性和晶系条件正确，且弹性矩阵存在负本征值或关键 Born 组合为负，那么负值可以解释为该结构在小形变范围内力学不稳定。此时继续报告杨氏模量或泊松比意义不大，应先说明该相或该构型不满足力学稳定条件。
图5 把声子谱和态密度作为稳定性分析的一部分。真实失稳往往不只出现在弹性张量中，也可能伴随 Γ 点或有限波矢的软模。若弹性判据失败而声子谱也出现明显虚频，结构重优化后转向低对称相的可能性更高。真实失稳需要由弹性曲率、声子响应和结构弛豫共同支撑。
图5 满足热力学、力学和动力学稳定条件结构的声子谱与声子态密度。该图说明弹性稳定性应与声子检查共同使用。DOI: 10.1038/s41598-025-02046-3。
什么时候更像数值问题？
若负值只在某个应变幅度出现，换更严格 EDIFF 或更密 k 点后消失，或者正负应变拟合曲线明显散乱，通常应先定位为数值或拟合问题。若含磁性元素的体系在不同初始磁矩下总能相差明显，也要先确定磁基态再重复弹性计算。
最终结果可以分成三类：一是数值设置导致的负值，修正收敛和拟合后消失；二是结构假设导致的负值，重新优化或降低对称性后进入新结构；三是真实力学失稳，需明确该构型不适合作为后续吸附、催化或输运分析的基底。负值排查的目标是给出来源分类，而不是只在后处理中把负号抹掉。
要点：弹性常数负值可能来自真实失稳，也可能来自结构弛豫、应变幅度、对称性或拟合口径。
要点：非对角项负值不一定等于失稳，关键要看晶系 Born 条件、矩阵本征值和派生模量。
要点：软材料、层状体系和低对称缺陷超胞更需要测试应变幅度、k 点和应力收敛。
要点：若弹性稳定性失败，应结合声子谱和结构重优化判断是否存在真实软模。
