| 查看: 1831 | 回复: 4 | |||
[交流]
算DOS中非自洽计算应该怎样理解？
|
|
我的理解是读取CHGCAR以后就根据该CHGCAR作非自洽计算，不进行自洽迭代，但算DOS分两步计算时第二步读取第一步的CHGCAR以后为什么还要进行自洽迭代计算？到底该怎样理解非自洽计算？（难道说算DOS就是增大K点加上ISMEAR = -5再读取CHGCAR减少计算量做更精确的自洽计算而已，然后得到DOS图） 第二个问题只用一步计算小体系的DOS时到底是根据哪个参数决定它做非自洽计算？是ISMEAR = -5 吗？ 如有解惑小弟拜谢 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
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
投票:
