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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于非自洽计算的几点疑惑。
91/1返回列表
查看: 2667  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

jugengfans

金虫 (正式写手)

[求助] 关于非自洽计算的几点疑惑。 已有1人参与

1.在计算band的非自洽计算过程中，为什么费米能级和自洽计算的有差异？
  我个人理解，因为在计算band的时候，K点在布里渊区中并不是MK那样遍历的，有的高对称K点可能会遗漏，所以在算占据数时会出现问题。
是这样吗。
2.为什么非自洽计算也会有迭代过程？
你在静态计算的时候设置ICHARG=11，把电荷密度fix了，为什么vasp不是一步出结果，还在不停的迭代？
      DFT的理论不是，你只要知道了CHGCAR就知道KS方程中的单电子势，就可以求解KS方程。
你非自洽计算的时候CHGCAR是固定的，那么应该带入KS方程直接就可以得到分子轨道，再从低到高算占据数，最后求和得到总波函数，事情不就完了吗？
迭代的是什么呢？
    是不是因为，波函数其实有个mix在里面，也就是说VASP在迭代前预先搞了个波函数，然后和KS里的波函数线性组合了一下？
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-05-06 22:27:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)
顶起来。 应该是和mix有关系 个人认为
赞一下 回复此楼
2楼2016-05-08 23:48:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

obaica

禁言 (著名写手)
本帖内容被屏蔽
3楼2016-05-09 00:30:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)
谢谢楼上顶贴。。
回复此楼
4楼2016-05-11 10:48:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★
jugengfans: 金币+2, 问题是电子步的SCF 2016-05-12 00:44:51
1.在计算band的非自洽计算过程中，为什么费米能级和自洽计算的有差异？
  我个人理解，因为在计算band的时候，K点在布里渊区中并不是MK那样遍历的，有的高对称K点可能会遗漏，所以在算占据数时会出现问题。
是这样吗。

> 是的

2.为什么非自洽计算也会有迭代过程？
你在静态计算的时候设置ICHARG=11，把电荷密度fix了，为什么vasp不是一步出结果，还在不停的迭代？
      DFT的理论不是，你只要知道了CHGCAR就知道KS方程中的单电子势，就可以求解KS方程。
你非自洽计算的时候CHGCAR是固定的，那么应该带入KS方程直接就可以得到分子轨道，再从低到高算占据数，最后求和得到总波函数，事情不就完了吗？
迭代的是什么呢？
    是不是因为，波函数其实有个mix在里面，也就是说VASP在迭代前预先搞了个波函数，然后和KS里的波函数线性组合了一下？

> 原因很單純~ 程式是死的，但~人是活的
    簡單的說，"非自洽" 顧名思意就是 "沒有疊代過程"
    VASP 計算要如何 "沒有疊代過程"?
    就是 NSW=0 或 1
    ICHARG = 11 這個參數?K不會影響到 NSW
    ICHARG 只是告訴 VASP 如何處理 initial charge
    所以如果你有給 NSW 大於 1，又給了 ICHARG = 11
    VASP 還是會去做 self-consistent calculation
赞一下(2人) 回复此楼
5楼2016-05-11 20:29:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jugengfans

金虫 (正式写手)
引用回帖:
5楼: Originally posted by jpchou at 2016-05-11 20:29:41
1.在计算band的非自洽计算过程中，为什么费米能级和自洽计算的有差异？
  我个人理解，因为在计算band的时候，K点在布里渊区中并不是MK那样遍历的，有的高对称K点可能会遗漏，所以在算占据数时会出现问题。
是这样 ...

谢谢解答，不过你说的是离子步自洽，而我想知道的是设置NSW=0的时候，电子步为什么同样会进行SCF.
赞一下 回复此楼
6楼2016-05-12 00:11:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)
喔~ 抱歉~ 誤會了
回复此楼
7楼2016-05-12 03:18:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buddy_84

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by jugengfans at 2016-05-08 23:48:47
顶起来。 应该是和mix有关系 个人认为

楼主，请问关于“非自洽计算，也有scf迭代的问题”解决了吗?
赞一下 回复此楼
8楼2020-07-21 11:02:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大杰大烤鸭

新虫 (小有名气)
请问这问题里的MK是个啥？
赞一下 回复此楼
9楼2021-08-16 21:04:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jugengfans 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 AVl0JAsg9hyo 2026-08-29 5/250 2026-08-30 19:57 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:46 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 18:51 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 18:29 by hZiFeudZyoGR
[博后之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 18:28 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 18:18 by hZiFeudZyoGR
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[基金申请] 能否申诉？ +4 echo8914667 2026-08-30 4/200 2026-08-30 16:03 by sundiv
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 看板上这么多中的，有点像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 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
关闭
请填处理意见
关闭