关于vaspkit201计算弹性常数

Xr.Ang

这是一个磁性结构为<100>AFMD的FCC结构Fe
我正在复现这个结构的文献上的弹性常数，文献里就是通过能量法计算的

我通过EOS+isif=2来优化晶胞
使用vaspkit201计算晶格常熟的时候还是会出现以下warning，而且出现了弹性常数为负的情况
  | Warning: The unstrained structure has non-zero elastic energy.|
| I strongly suggest you reperform structure optimization again.|

我去查看了STRAIN_ENERGY.dat文件，确实某些文件夹下的0结构弹性能不为0
但是如果我用isif=3优化原胞，那结构就和文献结构差太多了

我该怎么解决这个问题呢

有可能是对称性影响了弹性系数的计算吗vaspkit自动识别的空间群是129，由于磁性结构的影响，我的结构并不是严格的面心立方，我该如何去修改这个值呢