24小时热门版块排行榜    

查看: 3052  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dyc_2008

捐助贵宾 (著名写手)

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

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

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

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

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

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

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

一布福田,二修菩提
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
fzx2008: 金币+3, 谢谢指教! 2012-11-01 10:40:52
dyc_2008: 回帖置顶 2012-11-02 13:00:06

dyc_2008

捐助贵宾 (著名写手)

★ ★ ★ ★ ★
fzx2008: 金币+5, 有见解! 2012-10-31 18:39:56

stefsui

木虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖

真情无价108

木虫 (小有名气)

普通表情 高级回复(可上传附件)
信息提示
请填处理意见