24小时热门版块排行榜    

查看: 1513  |  回复: 4

dontgogent

新虫 (小有名气)


[交流] 算DOS中非自洽计算应该怎样理解?

我的理解是读取CHGCAR以后就根据该CHGCAR作非自洽计算,不进行自洽迭代,但算DOS分两步计算时第二步读取第一步的CHGCAR以后为什么还要进行自洽迭代计算?到底该怎样理解非自洽计算?(难道说算DOS就是增大K点加上ISMEAR = -5再读取CHGCAR减少计算量做更精确的自洽计算而已,然后得到DOS图)
第二个问题只用一步计算小体系的DOS时到底是根据哪个参数决定它做非自洽计算?是ISMEAR = -5 吗?
如有解惑小弟拜谢
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

dontgogent

新虫 (小有名气)


2楼2020-03-31 17:40:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dontgogent

新虫 (小有名气)


3楼2020-03-31 17:40:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dontgogent(金币+1): 谢谢参与
1、非自洽指的是电荷密度不变来得到本征值和本征态(具体过程其实就是通过Harris-Foulkes泛函代替求解KS方程)。非自洽中迭代的原因有两个:
(1)正常自洽过程前n步其实做的也是非自洽,但是也要迭代,原因在于我们一般从原子密度做为初始密度。然而,空轨道没有定义,这样空轨道信息是采用随机数来初始化的,需要通过几步非自洽的迭代调整空轨道再进行后续自洽计算。
(2)DOS或者能带的非自洽是因为,你自洽采用的k网格和展宽可能不同于自洽过程中的,这些在自洽中不包含的k对应的本征值和本征态就得通过迭代得到。
你的理解是正确的,就是为了降低k点的使用,否则那些额外的k点对应的态都是需要求解KS方程来得到,还是非常昂贵的。杂化泛函例外,只能通过自洽的方式做。如果你觉得体系比较小,直接取比较密集的k点自洽做DOS是可以的;能带的情况也可以把0权重的k-path添加到MP网格后面,和杂化泛函一样的做法,这些在原理上都是没问题的。
2、理解错误,前面提到了,只用一步的过程其实得到的DOS是直接自洽得到的,原理上允许这么做。

» 本帖已获得的红花(最新10朵)

4楼2020-04-06 07:14:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dontgogent

新虫 (小有名气)


送红花一朵
引用回帖:
4楼: Originally posted by 卡开发发 at 2020-04-06 07:14:02
1、非自洽指的是电荷密度不变来得到本征值和本征态(具体过程其实就是通过Harris-Foulkes泛函代替求解KS方程)。非自洽中迭代的原因有两个:
(1)正常自洽过程前n步其实做的也是非自洽,但是也要迭代,原因在于我 ...

谢谢您解惑~

发自小木虫IOS客户端
5楼2020-04-14 11:56:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dontgogent 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见