pymatgen2026.9.24：从实际硅应变数据到弹性张量

解压进入pymatgen-elastic，在已分配CPU作业内运行bash run.sh；或sbatch job.slurm。CPU分区按56核整节点分配和计费，分析仅使用一个线程，不启动新的DFT。不要在登录节点运行生产任务。

inputs/actual-stress-strain.json来自20个已完成的集群VASP硅计算：固定/弛豫离子各包含xx和xy五个正负小应变。应力已经由ASE转换为拉伸为正的GPa，顺序xx,yy,zz,yz,xz,xy；不要再反号。inputs/Si-POSCAR是同一实际未变形8原子晶胞。

本例按立方近似分别拟合C11、C12、C44，并重建张量。实际参考晶胞有微小非立方残差，这不是一般晶系21个独立常数的完整拟合；不得把小残余分量当成已独立计算的其他弹性常数。工程剪切gamma对应epsilon_xy=gamma/2。

ElasticTensor给出Voigt/Reuss/Hill、柔度、方向量、声速和弹性Debye温度；Strain/Deformation/Stress演示工程剪切、Green-Lagrange应变和应力映射。directional_elastic_mod返回C_nnnn方向纵向刚度，不是自由横向收缩的方向杨氏模量；后者由1/S_nnnn得到。y_mod的单位是Pa，Cij和B/G为GPa。连续介质声速、弹性Debye温度不是声子色散。

数据与完整张量在results/，实际新图在plots/，原始参考图在reference-plots/。换体系时用新的数据单独建目录，匹配完整应变集合、对称性、单位和零应力参考，不沿用固定硅参考。
