去离子水 发表于 2015-10-28 09:54
拟合eos得到平衡体积后，在POSCAR的第二行写上-V0，然后做isif=3。最后从CONTCAT里计算a，c的值，这种做 ...

可以，你也可以改变形状和改变体积交替进行，不过工作量相应大一些。非立方晶格的情形fitting EOS没准不同的c或是a/c的情形对应的v0可能会不一样，这点非常不容易处理。

卡开发发 发表于 2015-10-27 23:44
差不多，一般这样的精度也就行了。其实EOS也只是个近似，晶体内部的实际变量应该不止三个晶格常数这样，说 ...

拟合eos得到平衡体积后，在POSCAR的第二行写上-V0，然后做isif=3。最后从CONTCAT里计算a，c的值，这种做法对吗？还有没有其他方法在知道平衡体积后得到a，c的值

差不多，一般这样的精度也就行了。其实EOS也只是个近似，晶体内部的实际变量应该不止三个晶格常数这样，说白了就是做不同的晶格常数下晶格的能量，然后找能量极小点，fitting EOS相当于对这些点内插得到极值点而已（dE/dΩ=0）。

原则上来说isif当中isif=3涉及到的优化方式可靠的话能够得到准确的结构而不必fitting EOS，前提是stress和force在几个参数的收敛性测试条件下（通常参考能量，但实际远远不够），当然这有可能也只是理论上讲讲，实际可能影响的因素很复杂。如果再考虑这么复杂的收敛性测试，看起来反而fitting EOS更加省力一些。