VASP计算时，模型建立完之后对盒子大小的优化

卡开发发

LDAUTYPE = 2只要任何U-J相同的情况U和J任意取值结果都一样。

zdt1991

好的，谢谢

zdt1991

卡开发发 发表于 2017-3-23 22:01
LDAUTYPE = 2只要任何U-J相同的情况U和J任意取值结果都一样。

好的，谢谢
非常感谢您的指导。

珊珊来迟

卡开发发 发表于 2017-3-4 16:00
虽2D系统但在vasp使用slab模型周期性仍然是3D的，c方向如果随之变化就会引进其他因素影响能量，所以只要a、 ...

老师关于优化晶格常数，我有几个问题想请教您，是一个三斜体系，有实验单晶（实验晶格常数为6.9988 7.2813 11.2715），以此为基础进行优化：1.若采用优化缩放系数的方法，我是这样设置的您看可以吗？
for a in  $(seq 0.98 0.01 1.07)                              
do                                                            
cat > POSCAR <<!                                             
crystal                                                      
$a                                                            
        6.9987998009         0.0000000000         0.0000000000
       -0.0369809416         7.2812061562         0.0000000000
       -2.3167224924        -2.7680103371        10.6779033364

。。。。。
用ISIF=2优化，取能量最低点对应的结构。相应的晶格常数为7.0688  7.3541  11.3842
2.若直接用ISIF=3的方法去优化晶格参数，
a）第一步已经用ISIF=3优化了晶胞形状和原子位置，还需要再用ISIF=2弛豫一下原子位置吗？

b）我见大师兄官网上的教程说用ISIF=3进行优化时需提高截断能，一般用600eV，我用这个方法优化出来a方向的晶格常数和实验差很多是为什么？（以实验得到的单晶作为POSCAR去优化的）
用ISIF=3优化得到的结果是7.2237  7.2351  11.2288
3.我知道的还有拟合状态方程的方法来获得晶格常数，所以我现在不太清楚怎么判断哪种方法得到的结构是最优的呢？
优化三斜晶系的晶格常数的话哪种方法比较常用呢？
期待您的回复！

卡开发发

珊珊来迟 发表于 2018-1-2 18:25
老师关于优化晶格常数，我有几个问题想请教您，是一个三斜体系，有实验单晶（实验晶格常数为6.9988 7.281 ...

1、调整缩放系数用isif=2优化的过程晶格的大小发生变化，但晶格的形状不会发生变化，这其实等同于引进了更多的约束条件，你可以考虑试试固定晶格体积而不固定形状的优化方式看看什么结果。
2、a)原则上不必重新用isif=2优化。
b)确实要用更高的动能截断和k点，最好测试到力和应力都收敛会比较稳妥。差很多的原因我暂时不好说，大概有可能：
1）自由度多了之后确实难优化一些。
2）可能是交换关联泛函选取的原因，比如体系是色散作用主导的可能一般的方法会很难算准确，需要特定的修正。
3）数值误差，比如动能截断和k点。
3、我不知道你所谓的拟合是哪种，拟合EOS就是调整缩放系数来拟合晶体状态方程。原则来说拟合EOS和优化晶格两种方法都是可以的。

珊珊来迟

卡开发发 发表于 2018-1-2 18:55
1、调整缩放系数用isif=2优化的过程晶格的大小发生变化，但晶格的形状不会发生变化，这其实等同于引进了 ...

谢谢您的回复！
老师用ISIF=4固定体积优化的时候需要保持晶格常数a b c的比例固定吗？
另外老师我不太清楚测试到应力收敛该怎么看？
我拟合EOS不是调整的缩放系数，是设置一系列不同的a，然后用isif=4优化，最后来拟合的。
那第一种调整缩放系数，然后用isif=4来优化，找能量最低点的方法其实和拟合EOS的方法是一样的了？

卡开发发

珊珊来迟 发表于 2018-1-2 20:24
谢谢您的回复！
老师用ISIF=4固定体积优化的时候需要保持晶格常数a b c的比例固定吗？
另外老师我不太 ...

不需要固定比例，a、b、c的比例会被优化。stress是一个3x3的矩阵，简单点你也可以看pressure。拟合EOS找结构可以看成就是找一个合适的函数来插值得到最低点。

珊珊来迟

卡开发发 发表于 2018-1-2 21:45
不需要固定比例，a、b、c的比例会被优化。stress是一个3x3的矩阵，简单点你也可以看pressure。拟合EOS找 ...

那用这种方法可以得到各个缩放系数下对应的晶格常数和能量，然后是可以直接取能量最低点的结构作为最终的晶格常数，还是要继续用已经得到的各个体积和能量接着去做EOS拟合再得到最终结构？

卡开发发

珊珊来迟 发表于 2018-1-2 22:33
那用这种方法可以得到各个缩放系数下对应的晶格常数和能量，然后是可以直接取能量最低点的结构作为最终的 ...

取能量最低的就行。