VASP弹性常数计算： ISYM = 0 IBRION = 6 # CG Geom Opt ISIF = 3 # Relax atoms and cell NFREE = 4 要在不同的系统上运行它，只需要修改文件 init.mod （温度 结构数据） 和 potential.mod （势文件）。 为了正确计算弹性常数，必须注意在 init.mod 中指定正确的单位（units、cfac 和 cunits）。 验证应力分量的 MD 采样是否生成准确的统计平均值也很重要。 一个迹象是弹性常数对 init.mod 中变量 ${up} 的选择不敏感。 另一个是检查有限大小效应。 relaxBox 0压0温度/data/gpfs/assoc/matersimul/kunluo/Mg/dpgen/MD/relaxBox（in.elastic里有relaxbox的命令，所以不是必须的） 计算elastic 0温度 Copy需要的文件 displace.mod(不变), in.elastic(不变), , init.mod(结构信息), potential.mod(势文件), restart.Mgrelax, slurm 修改文件 init.mod,带温度的是不同的文件 Init.mod 区别：up 1.0e-3, 例子是1.0e-6，整体感觉和高温的差很多，还是用 1.0e-3 (需要多算几个测试一下 建议可以提高到 1.0e-2, castep 默认是0.003 即3.0e-3 ), 有温度就用2.0e-2 用read_restart restart.npt 替换read_data data.CdTeconcell，直接读取之前的结果 read_restart restart.Mgrelax cha...