计算化学公社

 找回密码 Forget password
 注册 Register
Views: 11202|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 改变k点后还能继续使用之前的波函数和电荷密度么?

[复制链接 Copy URL]

63

帖子

0

威望

213

eV
积分
276

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
我先用一个较低的k点优化结构并得到了WAVECAR和CHGCAR,然后我将k点调高并设置ISTART=1, ICHARG=1。
我自己试了一下,发现k点调高后ISTART=1, ICHARG=1得到的结果和ISTART=0, ICHARG=2得到的结果结构上完全一样,能量之差1E-5个eV。
能不能确定这样做是没问题的?

3753

帖子

3

威望

1万

eV
积分
19671

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2017-9-14 14:38:35 | 只看该作者 Only view this author
如果提升k点后ISTART=0和ISTART=1做构型优化完全一样,能量差别也很小,那么没有关系。但一般情况下不推荐这样做,有时候有可能导致构型稍有差异(肉眼未必能看出来)。一般而言,ISTART=1和ICHARG=11会通过波函数构造密度或直接读取密度来进行非自洽迭代(一般通过Harris functional的方法之类的),如果收敛标准也不足够高,有可能会收敛不到真正的基态上。

不太建议一开始随便去取一个低的k点,最好做一下收敛性测试,k点的收敛一般很快。如果体系需要的k很大,可以加一些展宽来加速k点收敛以减小k点使用,但不能太大,自洽算能量的时候减小展宽增大k点即可。

评分 Rate

参与人数
Participants 2
eV +5 收起 理由
Reason
sobereva + 1
ra2ghgzh + 4 谢谢

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

471

帖子

0

威望

1795

eV
积分
2266

Level 5 (御坂)

3#
发表于 Post on 2017-9-15 14:12:00 | 只看该作者 Only view this author
其实吧,这样做也没什么不可以,俺在旧社会时曾经这么做过,测试好K点后,接下来的计算,先用较少的k点,快速获得波函数,然后在此基础上正常K点下计算。当时条件差,优化个结构都是以天为单位计算,先少量k点先算得波函数,再正常k点计算,目的就是省点时间。
现在大家都富豪了,没人愿意这么“多此一举“了

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 15:44 , Processed in 1.350541 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list