VASP 结构性质结果中，C₁₁、C₁₂、C₄₄、B、G、E、ν 经常并列列出，因此容易被误解为同一层级的力学量。二者对应的物理层级不同：Cᵢⱼ 描述材料对不同方向、不同剪切方式的微小应变如何响应；B、G、E、ν 则由弹性张量经过平均或方向投影得到。弹性常数是张量元素，弹性模量是由张量换算得到的力学响应指标。
这个区别在各向同性材料里不明显，因为少数几个常数就能描述整体响应；在晶体、二维材料、层状材料和各向异性催化载体里，它会直接影响结果解释。仅报告杨氏模量较高，不能说明哪个晶向抗拉更强；仅列出 C₁₁、C₁₂、C₄₄，也不能完整回答材料整体抗压、抗剪或延展趋势。
图1 弹性常数高通量计算与筛选流程。图中从结构优化、施加应变到过滤弹性张量，说明 Cᵢⱼ 是后续模量分析的基础输入。DOI: 10.1038/sdata.2015.9。
在小形变范围内，晶体通常采用线弹性近似。VASP 或其他第一性原理程序常通过施加一组很小的应变，计算体系应力或应变能变化，再反推出弹性张量。对三维晶体，常用 Voigt 记号把六个应变分量写成向量，弹性张量就表示为 6 × 6 矩阵。Cᵢⱼ 的核心含义是某个应变分量引起某个应力分量的响应强度。
σᵢ = Σⱼ Cᵢⱼ εⱼ 式(1)
式(1)中 σᵢ 是应力分量，εⱼ 是应变分量，Cᵢⱼ 是弹性常数。C₁₁ 常和 x 方向拉伸或压缩响应有关，C₁₂ 体现一个方向拉伸对另一个方向应力的耦合，C₄₄ 则常对应剪切响应。不同晶系独立常数数量不同，立方晶系只需 C₁₁、C₁₂、C₄₄ 三个独立常数，低对称结构会有更多独立项。
C₁₁ 较大只能说明某个坐标方向上的纵向响应较强，不等于材料在所有方向都硬，也不等于整体抗剪能力强。层状材料可能面内 C₁₁ 很高、层间相关常数较低；二维材料可能面内拉伸很强，但面外方向并不适合直接按三维体模量理解。
因此，弹性常数提供方向分辨的响应信息。它保留了晶体对称性、晶向差异和拉伸-剪切耦合关系。如果研究对象是单晶、取向薄膜、二维片层或界面应变，Cᵢⱼ 比单个平均模量更适合描述方向相关形变。
图2 不同维度材料的弹性常数分布。图中展示 3D、2D、1D 和 0D 材料的 Cᵢⱼ 分布差异，说明弹性常数保留了维度和方向响应信息。DOI: 10.1103/physrevb.98.014107。
弹性模量不是独立于张量之外的计算量，而是从弹性常数或柔度矩阵换算出来的力学指标。体模量 B 描述材料抵抗体积压缩的能力；剪切模量 G 描述抵抗形状改变的能力；杨氏模量 E 描述单轴拉伸或压缩时的刚度；泊松比 ν 描述横向收缩与纵向拉伸的比例。B、G、E、ν 分别对应抗压、抗剪、抗拉和横向变形耦合。
S = C⁻¹ 式(2)
式(2)里的 C 是弹性常数矩阵，S 是柔度矩阵。很多方向相关的杨氏模量和泊松比要从 S 中得到；多晶平均模量则常用 Voigt、Reuss 和 Hill 平均。Cᵢⱼ 对应完整张量数据，B/G/E/ν 对应面向抗压、抗剪、抗拉和横向变形耦合的派生量。
因为同一弹性张量可以在不同平均假设下得到不同模量。Voigt 平均假设应变均匀，Reuss 平均假设应力均匀，Hill 平均取二者平均，常用于多晶近似。对各向异性单晶，还可以沿不同方向计算方向杨氏模量，因此一个材料不只有一个 E 值。
E = 9BG / (3B + G), ν = (3B − 2G) / [2(3B + G)] 式(3)
式(3)是各向同性或多晶平均口径下常用的关系。它说明 E 和 ν 可以由 B 与 G 推出，但前提是已经选择了平均方式。如果没有说明 Voigt、Reuss、Hill 或方向投影口径，只给出一个模量数值，信息是不完整的。
图3 体模量、剪切模量、泊松比和原子体积的分布图。图中把弹性常数进一步转化为模量和泊松比，用于比较大量晶体的整体力学响应。DOI: 10.1038/sdata.2015.9。
在立方晶系中，C₁₁、C₁₂、C₄₄ 是最常见的三项。C₁₁ 反映沿晶轴方向拉伸或压缩时的纵向刚度；C₁₂ 反映不同正应变方向之间的耦合；C₄₄ 反映剪切形变响应。三者都合格，才构成立方晶体弹性张量的基本描述。
B = (C₁₁ + 2C₁₂) / 3 式(4)
式(4)只适用于立方晶系的体模量表达。它也提醒我们：体模量不是 C₁₁ 的同义词，而是 C₁₁ 与 C₁₂ 共同决定的平均压缩响应。C₁₂ 代表方向耦合，忽略它会把单轴刚度误读成整体抗压能力。
C₄₄ 对剪切变形敏感，和滑移、层间剪切、结构软化以及某些力学失稳模式有关。一个材料可以有较高的 C₁₁，却因为剪切响应较弱而表现出较低的 G。对层状材料、二维材料和弱相互作用体系，剪切相关常数与面内拉伸常数可能存在显著差异。
在催化和表面模型里，C₄₄ 不一定直接决定吸附能，但它会影响载体抗形变能力、界面应变释放和结构重构趋势。若研究应变调控催化活性，单看面内拉伸刚度不够，还要看剪切和方向耦合是否会改变局域配位。
图4 不同数据库中 Voigt 体模量和剪切模量的对比。体模量与剪切模量来自同一弹性张量的不同组合，二者不应混为一个力学指标。DOI: 10.1103/physrevb.98.014107。
VASP 相关后处理通常先得到弹性常数矩阵，再由 VASPKIT、pymatgen、ELATE 或自写脚本换算模量。合理的结果组织顺序是：先说明结构已经充分优化，再给出 Cᵢⱼ 矩阵和晶系独立项，随后说明采用哪种平均方式得到 B、G、E 和 ν。没有弹性张量来源的模量表，不能单独支撑力学性质判断。
如果讨论单晶或二维材料，建议保留方向信息。例如二维材料常报告 C₁₁、C₂₂、C₁₂、C₆₆，必要时给出 x/y 方向杨氏模量。三维多晶近似则可以给出 Bᵥ、Bᵣ、Bₕ、Gᵥ、Gᵣ、Gₕ，并说明最终讨论采用 Hill 平均。
二维材料的弹性常数常有 N/m 和 GPa 两种写法。N/m 是按二维面内刚度理解；GPa 往往需要假定一个有效厚度。若不同论文采用不同厚度定义，直接比较 GPa 数值会引入额外误差。对单层材料，报告 N/m 通常更符合二维力学量的定义。
这也是弹性常数和模量容易混乱的地方：同一个符号 C₁₁，在三维晶体里通常是 GPa，在二维材料里可能是 N/m。写结果时要把维度、单位和平均口径说清楚，否则看似相同的数值可能并不在同一比较条件下。
弹性模量描述小形变线性响应，不等于断裂强度，也不等于材料一定能合成。B 大说明抗体积压缩能力强，G 大说明抗剪切能力强，E 大说明单轴刚度高，ν 和 B/G 常被用来辅助判断延展或脆性倾向，但这些都是线弹性范围内的指标。
B/G > 1.75 常被用作延展倾向的经验参考 式(5)
式(5)属于 Pugh 比经验口径，不是绝对判据。它适合做材料筛选和趋势讨论，不能单独替代断裂、缺陷、温度、相稳定和实验制备条件。弹性模量能回答小形变刚度问题，不能直接等同于强度、韧性或可合成性。
吸附与催化体系里，弹性信息常用于三个地方。第一，判断载体或二维片层在应变调控下是否容易产生较大结构响应。第二，解释异质结晶格失配后的应变分配，哪一层更容易释放形变。第三，辅助讨论活性位附近重构趋势，尤其是缺陷、边缘和单原子锚定位点。
因此，弹性结果应与具体形变类型对应：如果研究界面应变，重点分析方向相关 Cᵢⱼ 和面内模量；如果研究高压或体相稳定性，体模量更重要；如果研究层间滑移或二维材料形变，剪切相关常数不能省略。应先明确材料承受的形变类型，再选择对应的弹性量进行解释。
图5 Pugh-Pettifor 经验区间与泊松比分布。图中展示 B/G 和 ν 可用于延展或脆性趋势分析，但这类模量派生指标应结合材料类别与结构证据使用。DOI: 10.1103/physrevb.98.014107。
概括而言，C₁₁、C₁₂、C₄₄ 是弹性张量中的局部坐标响应；B、G、E、ν 是用于描述抗压、抗剪、抗拉和横向收缩的换算指标。分析 Cᵢⱼ 可以保留方向信息，分析模量可以获得便于横向比较的材料性质。
要点：弹性常数 Cᵢⱼ 是应力-应变关系中的张量元素，保留晶向、剪切和耦合信息。
要点：弹性模量 B、G、E 和 ν 是从弹性张量或柔度矩阵换算出的整体或方向力学指标。
要点：C₁₁ 不等于杨氏模量，C₄₄ 不等于剪切模量；它们需要结合晶系和平均方式理解。
要点：二维材料、层状材料和各向异性晶体必须说明单位、方向和 Voigt/Reuss/Hill 平均口径。
