| 查看: 1864 | 回复: 4 | |||
[交流]
算DOS中非自洽计算应该怎样理解？
|
|
我的理解是读取CHGCAR以后就根据该CHGCAR作非自洽计算，不进行自洽迭代，但算DOS分两步计算时第二步读取第一步的CHGCAR以后为什么还要进行自洽迭代计算？到底该怎样理解非自洽计算？（难道说算DOS就是增大K点加上ISMEAR = -5再读取CHGCAR减少计算量做更精确的自洽计算而已，然后得到DOS图） 第二个问题只用一步计算小体系的DOS时到底是根据哪个参数决定它做非自洽计算？是ISMEAR = -5 吗？ 如有解惑小弟拜谢 |
回复此楼» 猜你喜欢
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
上海工程技术大学激光智能制造课题组招收博士研究生
已经有10人回复
-
申博发邮件
已经有11人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有12人回复
-
我的奶奶
已经有3人回复
2楼2020-03-31 17:40:14
3楼2020-03-31 17:40:25
★
dontgogent(金币+1): 谢谢参与
dontgogent(金币+1): 谢谢参与
|
1、非自洽指的是电荷密度不变来得到本征值和本征态（具体过程其实就是通过Harris-Foulkes泛函代替求解KS方程）。非自洽中迭代的原因有两个： （1）正常自洽过程前n步其实做的也是非自洽，但是也要迭代，原因在于我们一般从原子密度做为初始密度。然而，空轨道没有定义，这样空轨道信息是采用随机数来初始化的，需要通过几步非自洽的迭代调整空轨道再进行后续自洽计算。 （2）DOS或者能带的非自洽是因为，你自洽采用的k网格和展宽可能不同于自洽过程中的，这些在自洽中不包含的k对应的本征值和本征态就得通过迭代得到。 你的理解是正确的，就是为了降低k点的使用，否则那些额外的k点对应的态都是需要求解KS方程来得到，还是非常昂贵的。杂化泛函例外，只能通过自洽的方式做。如果你觉得体系比较小，直接取比较密集的k点自洽做DOS是可以的；能带的情况也可以把0权重的k-path添加到MP网格后面，和杂化泛函一样的做法，这些在原理上都是没问题的。 2、理解错误，前面提到了，只用一步的过程其实得到的DOS是直接自洽得到的，原理上允许这么做。 |
» 本帖已获得的红花（最新10朵）
dontgogent4楼2020-04-06 07:14:02
5楼2020-04-14 11:56:36