应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 算DOS中非自洽计算应该怎样理解？
51/1返回列表
查看: 1722  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dontgogent

新虫 (小有名气)

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

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

» 猜你喜欢

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

查看全部散金贴
1楼 2020-03-31 17:40:05
已阅   回复此楼   关注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是直接自洽得到的，原理上允许这么做。
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

dontgogent
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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭