| 查看: 3061 | 回复: 12 | |||||||||||||||||||
[交流]
关于原子自洽，电子自洽，非自洽三步运算的思考。已有10人参与
|
|||||||||||||||||||
|
先向各位大牛道歉， 因为上次我发了个帖子，被某牛人认为是侮辱了物理。 所以，这次要先道歉，如果我侮辱了物理，请各位严厉的指出。 我一定会进行严格的自我批评。 1、原子自洽，也就是结构优化，其实vasp也写了一个CHGCAR文件。但是，我们为什么非要用，电子自洽得出的CHGCAR呢？ 小虫认为，结构优化的时候，每个原子步，都是基于当前的POSCAR，进行一次电子自洽，自洽完后，计算各原子所受到的力，判断原子要不要移动移动。如此反复多次，最后一次，发现达到推出的要求了（步数达到了，或者精度达到了）。这时，我们可以让VASP写一个CHGCAR，而且，这个CHGCAR，就是依据最后的原子坐标‘电子自洽’出来的电荷密度。但是我们偏偏不用它，而要另外进行一次电子自洽。各位看官，这是否多余了呢？不是。因为：结构优化后，胞发生了变化，布里渊区发生了变化，所以截断动能对应的K，在K空间中区域不再是一个球，可能变成了椭球，也就是说，计算所用之基函数，是基于初始坐标的，对结构优化后的坐标，这个基函数有点不合时宜。所以，我们要进行一次电子自洽，得到我们要的电荷密度。 结论：不用结构优化的CHGCAR，不是因为这个CHGCAR没有自洽，而是因为基函数与原子坐标不匹配得不怎么好。 2、所以，计算态密度一定要以‘电子自洽’得到的CHGCAR为输入文件。 |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理 | VASP | 量子理论见解 | 科研 |
» 本帖已获得的红花（最新10朵）
» 猜你喜欢
谈谈两天一夜的“延安行”
已经有14人回复
-
职称评审没过，求安慰
已经有14人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有15人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
EST投稿状态问题
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
垃圾破二本职称评审标准
已经有10人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP静态自洽计算，能量升高7eV，然后局域磁矩消失，请问是什么问题？
已经有8人回复
- 关于赝势与自洽计算
已经有15人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- 静态自洽计算遇到的问题
已经有5人回复
- 自洽场的fortran程序
已经有11人回复
- 静态自洽计算能量时EDIFFG的问题
已经有14人回复
- 【求助】我在用ATK做自洽计算的时候，出现如下的错误：
已经有3人回复
- 【求助】VASP中，“结构优化”和“电子自洽计算”的关系？
已经有8人回复
- 【求助】结构优化之后的静态自洽计算
已经有12人回复
- 【讨论】自洽计算与非自洽计算的含义、所做的工作
已经有10人回复
- 【求助】SCF(自洽场理论)是怎么回事
已经有3人回复
- 【原创】并行自洽计算后得到和节点数一样多的波函数文件wfc
已经有4人回复
- 【求助】vasp的电子自洽和原子核优化过程是否可以分开？
已经有4人回复
- 【求助】静态自洽计算读WAVECAR有用吗？
已经有9人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
- 【求助】态密度计算究竟是自洽计算，还是非自洽计算？
已经有24人回复
- 【求助】ALGaN/GaN 薛定谔泊松方程自洽求解（插AlN层）
已经有5人回复
★ ★ ★ ★ ★
fzx2008: 金币+5, 有见解！ 2012-10-31 18:39:56
fzx2008: 金币+5, 有见解！ 2012-10-31 18:39:56
送鲜花一朵gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+3, 谢谢指教！ 2012-11-01 10:40:52
dyc_2008: 回帖置顶 2012-11-02 13:00:06
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+3, 谢谢指教！ 2012-11-01 10:40:52
dyc_2008: 回帖置顶 2012-11-02 13:00:06
stefsui
木虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 2252.3
- 散金: 114
- 红花: 1
- 帖子: 155
- 在线: 472.7小时
- 虫号: 872985
- 注册: 2009-10-15
- 专业: 理论和计算化学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
真情无价108
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1594.9
- 散金: 90
- 红花: 2
- 帖子: 296
- 在线: 179.2小时
- 虫号: 1516328
- 注册: 2011-11-29
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
贝芷苡
铜虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 262.1
- 散金: 1
- 帖子: 146
- 在线: 49.2小时
- 虫号: 2016934
- 注册: 2012-09-20
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
paopaotu326
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1023
- 散金: 5
- 红花: 5
- 帖子: 675
- 在线: 178.6小时
- 虫号: 951482
- 注册: 2010-02-02
- 专业: 金属材料
拨云见日no1
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 262.3
- 帖子: 38
- 在线: 69.8小时
- 虫号: 2463781
- 注册: 2013-05-14
- 专业: 凝聚态物性I:结构、力学和