24小时热门版块排行榜    

CyRhmU.jpeg
查看: 5140  |  回复: 24
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

miracle3130

铁杆木虫 (小有名气)

[交流] 【求助】态密度计算究竟是自洽计算,还是非自洽计算?已有9人参与

如题,对于共线或者是非共线的态密度计算,到底是自洽计算(ISTART=0, ICHARG=2或1),还是非自洽计算(ISTART=1,ICHARG=11)。
望指教!
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 第一原理资料汇编 vasp计算资料 第一性原理杂
VASP vasp 学习贴

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

★ ★
gavinliu7390(金币+2):谢谢解答! 2010-04-14 14:53
非自恰的。所谓自恰计算,一般是电子步计算的时候不断的迭代电子密度,得到能量最低的状态。
所以在计算态密度的时候,都是先进行一次ISMEAR=0的自恰计算,得到体系的最低能量状态。然后进行一次ISMEAR=-5(这个计算能量比较准确)非自恰计算得到体系的态密度,因为计算态密度不需要迭代电子步,只需要根据已有的CHARGE文件就可以得到特征值了,然后就可以计算态密度了。
所以,计算态密度至少在VASP中中是一次自恰计算,再计算一次非自恰计算。
10楼2010-04-11 15:28:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

★ ★
zzy870720z(金币+2):谢谢 2010-05-21 12:28:57
引用回帖:
Originally posted by xiao72379 at 2010-04-14 14:48:47:

比较肯定bingmou 兄的说法,不过我本人有2个疑惑请教:
1、先进行的自洽计算是不是ISMEAR一定是0呢?因为我用过-5,我师兄用过1
2、我查资料,包括问高年级师兄姐,大家都说第二步是非自洽的,不过就像上面的 ...

1,之所以要用ISMEAR=0是因为这样做相对计算量比较小,而且对于一些体系比较大的体系K点很小,不能用ISMEAR。而ISMEAR=-5的计算总能很准确,但是并不使用于弛豫,在手册的ISMEAR部分说道:The only drawback is that the method is not variational with respect to the partial occupancies. Therefore the calculated forces and the stress tensor can be WRONG by up to 5 to 10% for metals.
2. 如果你仔细看一下输出文件,会发现和自恰计算有一个不同,就是rms(c)一直是空的,而rms(c)表示的 the difference between input and output charge density. 也就是说电荷密度一直没有变化。其实在这里并没有电子电荷分布的mix。那为什么有那么多能量的变化呢?给VASP一个电荷分布,系统猜一个波函数,得到能带,得到总能,发现猜的波函数太次,rms很大,那是不是要修改波函数呢?所以,这么多步都只是在求波函数而已,不是在寻找最优的电荷分布。
14楼2010-05-20 23:47:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

引用回帖:
Originally posted by uuv2010 at 2010-04-11 17:23:45:

为什么非自洽计算的输出文件中还有电子步的迭代呢?

同上
15楼2010-05-20 23:47:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

引用回帖:
Originally posted by wei.ji at 2010-05-21 10:54:27:
算态密度的话,精度高的电荷密度貌似更重要些。

没看出你想说什么。计算DOS的时候必然用的是CHACAR文件而不是CHA
21楼2010-05-21 11:47:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 miracle3130 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见