24小时热门版块排行榜    

查看: 1832  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 中青基了要发朋友圈吗? +6 349506619 2026-08-28 6/300 2026-08-31 08:43 by kingkocxr
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见