| 查看: 2656 | 回复: 8 | ||
jugengfans金虫 (正式写手)
|
[求助]
关于非自洽计算的几点疑惑。 已有1人参与
|
|
1.在计算band的非自洽计算过程中，为什么费米能级和自洽计算的有差异？ 我个人理解，因为在计算band的时候，K点在布里渊区中并不是MK那样遍历的，有的高对称K点可能会遗漏，所以在算占据数时会出现问题。 是这样吗。 2.为什么非自洽计算也会有迭代过程？ 你在静态计算的时候设置ICHARG=11，把电荷密度fix了，为什么vasp不是一步出结果，还在不停的迭代？ DFT的理论不是，你只要知道了CHGCAR就知道KS方程中的单电子势，就可以求解KS方程。 你非自洽计算的时候CHGCAR是固定的，那么应该带入KS方程直接就可以得到分子轨道，再从低到高算占据数，最后求和得到总波函数，事情不就完了吗？ 迭代的是什么呢？ 是不是因为，波函数其实有个mix在里面，也就是说VASP在迭代前预先搞了个波函数，然后和KS里的波函数线性组合了一下？ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
2026国自然函评费到账
已经有19人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【转】裸奔的Barabasi
已经有10人回复
- 宇宙的存在方式和世界的因果关系
已经有6人回复
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
2楼2016-05-08 23:48:47
|
本帖内容被屏蔽 |
3楼2016-05-09 00:30:45
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
4楼2016-05-11 10:48:26
【答案】应助回帖
★ ★
jugengfans: 金币+2, 问题是电子步的SCF 2016-05-12 00:44:51
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 |
5楼2016-05-11 20:29:41
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
6楼2016-05-12 00:11:03
7楼2016-05-12 03:18:00
8楼2020-07-21 11:02:08
9楼2021-08-16 21:04:49
20