3. 晶体对称性与独立弹性常数数目

弹性常数张量 \(C_{ijkl}\) 在三维空间中名义上有 81 个分量，但由于应力张量与应变张量的对称性，即 \(C_{ijkl}=C_{jikl}=C_{ijlk}=C_{klij}\)，分量数先降到 36 个，可用 Voigt 记法写成 6×6 的 \(C_{mn}\) 矩阵。真正决定“独立”常数数目的，是晶体点群的对称操作：当晶体绕某轴旋转、镜面反射或反演后与自身重合时，弹性常数张量必须保持不变。每施加一个对称操作，就相当于给 36 个分量加上一组线性约束方程，能自由取值的分量随之减少。因此，独立弹性常数数目不是人为规定的，而是“点群对称性”直接推导出来的结果。

一、从 36 到独立数目：对称操作的约束机制

设点群中的对称操作用正交矩阵 \(R\) 表示，则弹性常数张量满足变换关系：

C'_{ijkl} = R_{ip} R_{jq} R_{kr} R_{ls} C_{pqrs} = C_{ijkl}

把 6×6 的 \(C_{mn}\) 矩阵按 Voigt 指标展开，逐个代入 \(R\)，就能得到分量之间的等式与零约束。例如立方晶系含四重旋转轴，绕 z 轴旋转 90° 会把 x 与 y 互换、把 xy 与 −xy 互换，于是 \(C_{11}=C_{22}\)、\(C_{44}=C_{55}\)，同时 \(C_{16}\)、\(C_{26}\) 等剪切耦合项被强制为零。把所有生成元操作依次作用，最终只剩 3 个独立常数：\(C_{11}\)、\(C_{12}\)、\(C_{44}\)。这就是“对称性越高，独立常数越少”的定量来源。

二、七大晶系的独立弹性常数数目

按对称性由高到低，七大晶系对应的独立弹性常数数目如下表。表中“独立数目”指在完整点群约束下，6×6 矩阵中真正可独立取值的分量个数；各晶系内部不同点群可能略有差异，表中给出的是该晶系典型点群的结果。

需要特别说明：四方与三方晶系中，部分点群因缺少水平镜面或反演，会多出一个独立分量（如四方 4、−4、422、4mm 为 6 个，4/m、4/mmm 为 6 个，而 −42m 为 6 个；三方 3、32 为 7 个，3m、−3m 为 6 个）。因此不能只凭“晶系”二字就断定数目，必须落到具体点群。第一性原理计算中，通常先由晶体结构确定空间群，再取空间群的平移子群商掉后得到点群，最后查表或程序自动生成独立常数集合。

三、立方、六方、四方、正交的矩阵形式与推导要点

下面给出四个常用晶系的 6×6 弹性矩阵形式，空白处为零。这些矩阵是后续设置应变、提取 \(C_{ij}\) 的直接依据。

立方（3 个独立）：
[ C11 C12 C12  0   0   0 ]
[ C12 C11 C12  0   0   0 ]
[ C12 C12 C11  0   0   0 ]
[  0   0   0  C44  0   0 ]
[  0   0   0   0  C44  0 ]
[  0   0   0   0   0  C44]

六方（5 个独立）：
[ C11 C12 C13  0   0   0 ]
[ C12 C11 C13  0   0   0 ]
[ C13 C13 C33  0   0   0 ]
[  0   0   0  C44  0   0 ]
[  0   0   0   0  C44  0 ]
[  0   0   0   0   0 (C11-C12)/2 ]

四方 4/mmm（6 个独立）：
[ C11 C12 C13  0   0   0 ]
[ C12 C11 C13  0   0   0 ]
[ C13 C13 C33  0   0   0 ]
[  0   0   0  C44  0   0 ]
[  0   0   0   0  C44  0 ]
[  0   0   0   0   0  C66]

正交（9 个独立）：
[ C11 C12 C13  0   0   0 ]
[ C12 C22 C23  0   0   0 ]
[ C13 C23 C33  0   0   0 ]
[  0   0   0  C44  0   0 ]
[  0   0   0   0  C55  0 ]
[  0   0   0   0   0  C66]

以六方晶系为例说明推导要点：六重轴沿 z 方向，绕 z 轴旋转 60° 使 x、y 方向等价，因此 \(C_{11}=C_{22}\)、\(C_{13}=C_{23}\)、\(C_{44}=C_{55}\)；同时 \(C_{66}\) 并非独立，而由 \(C_{66}=(C_{11}-C_{12})/2\) 给出，这是六方对称性带来的额外关系。四方 4/mmm 与六方的差别在于四重轴只要求 90° 等价，不再强制 \(C_{66}\) 与 \(C_{11}\)、\(C_{12}\) 绑定，因此多出 1 个独立常数。正交晶系三个方向互不等价，故 \(C_{11}\)、\(C_{22}\)、\(C_{33}\) 与 \(C_{44}\)、\(C_{55}\)、\(C_{66}\) 都各自独立，共 9 个。

四、在第一性原理计算中如何利用这些约束

实际计算时，并不需要把 21 个或 13 个分量全部独立算一遍。以立方晶系为例，只需施加三组独立应变：

1. 施加体积应变或单轴应变 ε₁，由应力 σ₁ 得到 C₁₁，由 σ₂ 得到 C₁₂；
2. 施加剪切应变 ε₄，由 σ₄ 得到 C₄₄；
3. 用 C₁₁、C₁₂、C₄₄ 即可写出完整 6×6 矩阵，其余分量由对称性自动确定。

对六方晶系，需要独立确定 C₁₁、C₁₂、C₁₃、C₃₃、C₄₄ 五个量，C₆₆ 由 (C₁₁−C₁₂)/2 得到。对四方 4/mmm，则需额外确定 C₆₆，共 6 个。对正交晶系，9 个分量各自独立，通常需要沿三个主轴分别施加正应变和剪切应变。计算流程中，先由晶体结构文件确定空间群，再查点群表得到独立常数集合，最后只对这些分量设计应变-应力或应变-能量拟合方案，可以显著减少计算量，也避免因误设应变而引入非物理的耦合项。

总结起来，晶体对称性对弹性常数的约束可以概括为一条主线：点群操作 → 张量不变条件 → 分量等式与零约束 → 独立常数数目。七大晶系中，立方 3 个、六方 5 个、四方 6 个、三方 6 或 7 个、正交 9 个、单斜 13 个、三斜 21 个。掌握这条主线，就能在拿到任意晶体结构时，快速判断该算哪些分量、不该算哪些分量，并为后续的弹性常数计算与验证提供明确的对称性判据。