| 查看: 2409 | 回复: 12 | |||||||
[交流]
关于原子自洽，电子自洽，非自洽三步运算的思考。已有10人参与
|
|||||||
|
先向各位大牛道歉， 因为上次我发了个帖子，被某牛人认为是侮辱了物理。 所以，这次要先道歉，如果我侮辱了物理，请各位严厉的指出。 我一定会进行严格的自我批评。 1、原子自洽，也就是结构优化，其实vasp也写了一个CHGCAR文件。但是，我们为什么非要用，电子自洽得出的CHGCAR呢？ 小虫认为，结构优化的时候，每个原子步，都是基于当前的POSCAR，进行一次电子自洽，自洽完后，计算各原子所受到的力，判断原子要不要移动移动。如此反复多次，最后一次，发现达到推出的要求了（步数达到了，或者精度达到了）。这时，我们可以让VASP写一个CHGCAR，而且，这个CHGCAR，就是依据最后的原子坐标‘电子自洽’出来的电荷密度。但是我们偏偏不用它，而要另外进行一次电子自洽。各位看官，这是否多余了呢？不是。因为：结构优化后，胞发生了变化，布里渊区发生了变化，所以截断动能对应的K，在K空间中区域不再是一个球，可能变成了椭球，也就是说，计算所用之基函数，是基于初始坐标的，对结构优化后的坐标，这个基函数有点不合时宜。所以，我们要进行一次电子自洽，得到我们要的电荷密度。 结论：不用结构优化的CHGCAR，不是因为这个CHGCAR没有自洽，而是因为基函数与原子坐标不匹配得不怎么好。 2、所以，计算态密度一定要以‘电子自洽’得到的CHGCAR为输入文件。 |
回复此楼» 收录本帖的淘贴专辑推荐
 第一性原理 | VASP | 量子理论见解 | 科研 |
» 本帖已获得的红花（最新10朵）
liangqi05» 猜你喜欢
2024国社科通讯评审
已经有17人回复
-
34岁读博士晚吗
已经有27人回复
-
最新消息：2024国自然下载文件名变了
已经有5人回复
-
我标书代码变了
已经有7人回复
-
Nature 11日发文，中国著名学者们称造假迫不得已
已经有6人回复
-
关于2023的收入
已经有35人回复
-
请问事业编制和年薪制冲突吗？
已经有8人回复
-
硕博巨婴，也许才刚刚开始
已经有74人回复
-
初始合伙人来啦！（生物试剂耗材标准品）
已经有14人回复
-
面青地会评时间
已经有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人回复
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
|
1, 这个CHGCAR是相当有用的。 其实在做scf的时候，初始电荷密度是有一定的选择的。如果初始电荷密度不太合理的话，scf需要很多次才能收敛。就是因为开始的时候我们不知道什么样的电荷密度是收敛的，所以一般我们都选用原子的电荷密度（基于原子轨道）。的确，在relax原子的时候，每一次由于原子位置或晶格的变化，会导致电荷密度或波函数有一个小的偏离，但这种偏离不大。所以我们在优化的时候还是用了上一次的CHGCAR。 验证这一点实际很简单，你可以做两个计算，一个是直接优化。 另一个是做一步优化，然后把CONTCAR拷贝成POSCAR，接着优化。如此反复。第二种办法的时间绝对比第一种时间长。 2, 尽管所有的教程都是推荐把CHGCAR作为输入。 但本人认为， 如果你不计较是否多花时间的话，可以不这么做。直接用非常稠密的k点做dos就可以。至于为什么，它的思想跟第一类似。先用少量的k点做scf，然后用非常dense的k点，这样时间会少一些。 但我认为，以现在机器的计算能力，完全没必要。 |
4楼2012-11-01 08:45:46
2楼2012-10-31 18:28:19
3楼2012-10-31 19:08:48
stefsui
木虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 2252.3
- 散金: 114
- 红花: 1
- 帖子: 155
- 在线: 472.7小时
- 虫号: 872985
- 注册: 2009-10-15
- 专业: 理论和计算化学
5楼2012-11-01 23:49:00
真情无价108
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1594.9
- 散金: 90
- 红花: 2
- 帖子: 296
- 在线: 179.2小时
- 虫号: 1516328
- 注册: 2011-11-29
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
6楼2012-11-02 11:21:56
7楼2012-11-03 09:39:38
贝芷苡
铜虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 262.1
- 散金: 1
- 帖子: 146
- 在线: 49.2小时
- 虫号: 2016934
- 注册: 2012-09-20
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
8楼2013-01-03 10:44:27
paopaotu326
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1075
- 散金: 5
- 红花: 5
- 帖子: 670
- 在线: 177.4小时
- 虫号: 951482
- 注册: 2010-02-02
- 专业: 金属材料
9楼2013-01-11 11:07:12
拨云见日no1
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 262.3
- 帖子: 38
- 在线: 69.8小时
- 虫号: 2463781
- 注册: 2013-05-14
- 专业: 凝聚态物性I:结构、力学和
10楼2013-05-20 10:10:41
15