几个可能原因，截断能和K点太小了没收敛，应变步长太大了，结构优化没达到平衡点。

Sihan Zhang 感谢您的回答，在计算之前我已进行过截断能和K点的收敛测试，几何优化依旧会出现这样的问题，但我将原本的BCC结构换成FCC之后反而没有出现这个问题，我百思不得其解。

Nana Chen 如果你很确认你的材料不特殊的话，可以试一下调应变步长，因为一般软件用的都是差分算二阶导数。MS上我记得比较好用的是CASTEP，Dmol3别人一般说有点问题，也可能是软件的问题。不过具体这个方便我计算做的不多，你可以试试在科音的论坛上问这个方面的问题，更容易得到有效的解答。