[求助] 求助：用vasp算B4结构的E－V线的时候，出现莫名拐点，困扰了很久，急求解决！！

计算过程如下： INCAR.relax: SYSTEM=HfN ENCUT=600 ISTART=0 ICHARG=2 ISMEAR=1;SIGMA=0.18 IBRION=2 NSW=100 POTIM=0.2 ISIF=4 EDIFF=1E-5 EDIFFG=-1E-3 PREC=Accurate INCAR.static: SYSTEM=HfN ENCUT=600 ISTART=0 ICHARG=2 ISMEAR=-5 PREC=Accurate   K点取19*19*15，以Gamma点为中心。 运行脚本为： for i in `seq -w -53 -1 -70`; do cat > POSCAR <<!    hcp                                   $i 0.5000000000000000        -0.86600000000000000     0.000000000000000       0.5000000000000       0.86600000000000000    0.000000000000000        0.000000000000000       0.000000000000000   1.65715000000000000   2  2   Direct 0.33333333  0.66666666   0.00000000 0.66666666  0.33333333   0.50000000 0.33333333  0.66666666   0.37500000 0.66666666  0.33333333   0.87500000 ! cd /public/users/jy002/chenlongqing/clq/B4/Vother cp INCAR.relax INCAR vasp（运行命令） cp CONTCAR POSCAR cp INCAR.static INCAR vasp（运行命令） 也就是，根据参考文献取一系列体积，然后对每个体积进行基矢和原子坐标优化，初始的c/a和u因子都是从参考文献查找确定。 但是结果出现了很头痛的问题： 体积从57到70的优化和能量都很不错很平滑，但是体积56却是一个突变点，能量开始骤降，得到的EV图如下 我查看了体积为56的优化过程，发现vasp把它优化成了正交晶系(base centered orthorhombic)，并且c/a也突变到了1.32（本来应该是1.65左右）。我用了很多方法(比如加入SYMPREC=1E-4，或者调小EDIFFG)，但是都不能让vasp明白我计算的是六方晶系，这点大家能给点意见吗？ 另外，我怀疑是不是在体积为56时对应压强下发生了相变，所以会出现这样的情况，所以我看了下我的相图(相图我是只取体积为57之后的正常的EV线拟合得出的)，发现当B4结构体积为56时，对应的压强为23.2GPa，而相变图上在18GPa的时候B1和B4会发声相变，不知道这点有没有联系，还望大家帮我分析一下，非常感谢啊！十万火急啊！！！拜谢！！！

回复此楼