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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于原子自洽，电子自洽，非自洽三步运算的思考。
131/2返回列表12下一页
查看: 3050  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dyc_2008

捐助贵宾 (著名写手)

[交流] 关于原子自洽，电子自洽，非自洽三步运算的思考。已有10人参与

先向各位大牛道歉，
因为上次我发了个帖子，被某牛人认为是侮辱了物理。
所以，这次要先道歉，如果我侮辱了物理，请各位严厉的指出。
我一定会进行严格的自我批评。

1、原子自洽，也就是结构优化，其实vasp也写了一个CHGCAR文件。但是，我们为什么非要用，电子自洽得出的CHGCAR呢？
小虫认为，结构优化的时候，每个原子步，都是基于当前的POSCAR，进行一次电子自洽，自洽完后，计算各原子所受到的力，判断原子要不要移动移动。如此反复多次，最后一次，发现达到推出的要求了（步数达到了，或者精度达到了）。这时，我们可以让VASP写一个CHGCAR，而且，这个CHGCAR，就是依据最后的原子坐标‘电子自洽’出来的电荷密度。但是我们偏偏不用它，而要另外进行一次电子自洽。各位看官，这是否多余了呢？不是。因为：结构优化后，胞发生了变化，布里渊区发生了变化，所以截断动能对应的K，在K空间中区域不再是一个球，可能变成了椭球，也就是说，计算所用之基函数，是基于初始坐标的，对结构优化后的坐标，这个基函数有点不合时宜。所以，我们要进行一次电子自洽，得到我们要的电荷密度。

结论：不用结构优化的CHGCAR，不是因为这个CHGCAR没有自洽，而是因为基函数与原子坐标不匹配得不怎么好。
2、所以，计算态密度一定要以‘电子自洽’得到的CHGCAR为输入文件。
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理 VASP 量子理论见解 科研

» 本帖已获得的红花（最新10朵）

liangqi05

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
一布福田，二修菩提
1楼2012-10-31 18:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dyc_2008

捐助贵宾 (著名写手)
★ ★ ★ ★ ★
fzx2008: 金币+5, 有见解！ 2012-10-31 18:39:56
自己顶自己，乱谈自洽非自洽

结构自洽就是我们照合影时，优化每个人所站的位置
静态自洽，就是你站的位置定了。大家再来优化一下各自的pose
好把整张照片弄的河蟹一点。在微观世界，可以认为 ： 电荷密度 就是这张我们需要的合影

至于态密度，就是来分析下，共有多少个pose，没个pose有多少人在摆。
V手势3人
吐舌头 5人
。
。
。
高难度的一字马， 0 人
更到难度的倒立，0人

难度越高，在微观世界可认为是能级越高。虽然有这个态的存在，但是这个态上却没有电子。
赞一下(5人) 回复此楼
一布福田，二修菩提
2楼2012-10-31 18:28:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

liangqi05

铜虫 (小有名气)
送鲜花一朵
高手呀高手
回复此楼
把握机遇
3楼2012-10-31 19:08:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★ ★ ★ ★
小木虫: 金币+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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stefsui

木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
学习下
对于scf只知道个大概～～惭愧
赞一下 回复此楼
5楼2012-11-01 23:49:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

真情无价108

木虫 (小有名气)
初学者，感觉有帮助
赞一下 回复此楼
6楼2012-11-02 11:21:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

T_liyi_T

铜虫 (小有名气)
受教了
回复此楼
7楼2012-11-03 09:39:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贝芷苡

铜虫 (小有名气)
明白好多
回复此楼
坚持每天进步一点点
8楼2013-01-03 10:44:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paopaotu326

金虫 (正式写手)
比喻的真好，呵呵
回复此楼
9楼2013-01-11 11:07:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

拨云见日no1

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by dyc_2008 at 2012-10-31 18:28:19
自己顶自己，乱谈自洽非自洽

结构自洽就是我们照合影时，优化每个人所站的位置
静态自洽，就是你站的位置定了。大家再来优化一下各自的pose
好把整张照片弄的河蟹一点。在微观世界，可以认为 ： 电荷密度 就是 ...

很形象，受教了
回复此楼
10楼2013-05-20 10:10:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dyc_2008 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭