24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

miracle3130

铁杆木虫 (小有名气)

应助: 0 (幼儿园)
金币: 5932.2
帖子: 182
在线: 20.7小时
虫号: 716393
注册: 2009-03-06

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

如题，对于共线或者是非共线的态密度计算，到底是自洽计算（ISTART=0, ICHARG=2或1），还是非自洽计算（ISTART=1,ICHARG=11）。
望指教！

回复此楼

» 猜你喜欢

1楼2009-09-26 10:09:36

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

bingmou

金虫 (著名写手)

1ST强帖: 5
应助: 2 (幼儿园)
贵宾: 0.446
金币: 431.9
散金: 171
红花: 17
帖子: 1207
在线: 114.8小时
虫号: 487243
注册: 2007-12-28
专业: 材料科学/半导体
管辖: 第一性原理

★ ★
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很大，那是不是要修改波函数呢？所以，这么多步都只是在求波函数而已，不是在寻找最优的电荷分布。