计算化学公社

标题: 求助：VASP+phonopy DFPT方法算二维材料声子谱总有一个超大虚频 [打印本页]

作者
Author: helpcal 时间: 2023-10-20 03:20
标题: 求助：VASP+phonopy DFPT方法算二维材料声子谱总有一个超大虚频
本帖最后由 helpcal 于 2023-12-30 23:58 编辑

材料是稳定的结构，不同文献都有报道且有用2种方法得到的声子谱。最初的文献是用661超胞的finite difference方法计算，我用此方法能得到一样结果。但是DFPT始终失败，算完之后OUTCAR里最后一个虚频总是很大。我算两者的INCAR里区别就是NSW和IBRION，其他输入参数一样。也有文献是用DFPT算的。算声子的INCAR已贴。单独的那个声子谱是提高精度之后的结果，即力收敛是EDIFFG=-1E-08.但是结果只是细节上的收敛更好，但那个超大虚频始终存在并且没什么变化。

想请教一下各位这个问题如何解决。基于以上个人认为精度方便的提高只是在细节方面改善，而不会产生巨大改变比如去除这个超大虚频。

补充一下，对于221原胞，我加上了finite difference的图（第一张），其gamma点频率几乎是0，对比DFPT的超大虚频，感觉是不是DFPT这里方法哪有问题？

作者
Author: waitingseven 时间: 2023-10-20 08:58
DFPT方法用了什么版本的VASP，其次是gam还算standard，k点用的啥

作者
Author: helpcal 时间: 2023-10-20 09:04

waitingseven 发表于 2023-10-20 08:58
DFPT方法用了什么版本的VASP，其次是gam还算standard，k点用的啥

5.4.4的vasp。vasp_std. 原胞k点是12 12 1. 221超胞K点是661(刚试了下12 12 1的结果一样），以此类推，到661胞k点是2 2 1。

作者
Author: ahxb 时间: 2023-10-20 09:29
ENCUT取250是否足够？

作者
Author: helpcal 时间: 2023-10-20 09:39

ahxb 发表于 2023-10-20 09:29
ENCUT取250是否足够？

文献就是250。当然我也测试过400的结果也是一样。这些精度方面参数各种都调过，对结果影响没有很大。

作者
Author: helpcal 时间: 2023-10-21 21:58
发现问题是不能加dipol correction。请问有人知道原因么？ vdw修正不能加，加了直接报错算不了这个是知道的。VASP的DFPT限制挺大啊，优势看起来只有用起来方便一点。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)