计算化学公社

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

[VASP] vasp中如何使用PBEsol泛函?

[复制链接 Copy URL]

256

帖子

0

威望

1124

eV
积分
1380

Level 4 (黑子)

在vasp中使用PBEsol泛函时,除了在INCAR文件中输入GGA=PS外,是否还要把POTCAR文件中的LEXCH=91改为LEXCH=PS(我的POTCAR文件是从vasp赝势中的PAW_GGA目录中拷贝的),谢谢解答!

471

帖子

0

威望

1795

eV
积分
2266

Level 5 (御坂)

2#
发表于 Post on 2018-12-17 18:54:45 | 只看该作者 Only view this author
PBE、RPBE、PBEsol都需要用PBE版的POTCAR,不可以用PW91的
用PBEsol,两个都要改,否则会有警告,说你强迫了VASP,输出的能量也差别很大
两种设置,其实你可以测试一下,再读一下OUTCAR,对比一下看

256

帖子

0

威望

1124

eV
积分
1380

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2018-12-17 19:30:13 | 只看该作者 Only view this author
本帖最后由 zhongyuabc 于 2018-12-17 19:37 编辑
hakuna 发表于 2018-12-17 18:54
PBE、RPBE、PBEsol都需要用PBE版的POTCAR,不可以用PW91的
用PBEsol,两个都要改,否则会有警告,说你强迫 ...

好的,谢谢!对于POTCAR中把PBE改为PBEsol只需要改LEXCH这一行就行了,是吗?POTCAR中的其他部分应该不用再改了吧?

34

帖子

0

威望

1003

eV
积分
1037

Level 4 (黑子)

4#
发表于 Post on 2022-3-26 10:45:15 | 只看该作者 Only view this author
您好,最近在尝试Cu(111)表面上吸附CO用不同泛函的测试。RPBE, revPBE和PBEsol是否只需要更改INCAR的GGA参数,POTCAR均用PAW_PBE的赝势即可呢?
越学越难顶。

81

帖子

1

威望

220

eV
积分
321

Level 3 能力者

5#
发表于 Post on 2022-5-14 22:16:13 | 只看该作者 Only view this author
hakuna 发表于 2018-12-17 18:54
PBE、RPBE、PBEsol都需要用PBE版的POTCAR,不可以用PW91的
用PBEsol,两个都要改,否则会有警告,说你强迫 ...

但这个修改LEXCH的做法是错误的。

参考vasp wiki的说明,见https://www.vasp.at/wiki/index.php/GGA

Important: VASP recalculates the exchange-correlation energy inside the PAW sphere and corrects the atomic energies given by the POTCAR file. For this to work, the original LEXCH tag must not be modified in the POTCAR file.

这个如果修改了LEXCH,不仅仅是能量会发生改变,甚至连lattce constant也会发生改变,具体原因不清楚。
虽然这个改变了的晶格常数似乎与实验值更加接近,但如果按vasp wiki的说明来看,这样做,应该是错误的。

如果需要修改,那么vasp wiki应该会明确说明的。没有说明,证明是不用修改。而且也特别指出不能修改。

9

帖子

0

威望

89

eV
积分
98

Level 2 能力者

6#
发表于 Post on 2022-7-18 21:16:11 | 只看该作者 Only view this author
Jiongci-21 发表于 2022-3-26 10:45
您好,最近在尝试Cu(111)表面上吸附CO用不同泛函的测试。RPBE, revPBE和PBEsol是否只需要更改INCAR的GGA ...

同问请问有答案了嘛

34

帖子

0

威望

1003

eV
积分
1037

Level 4 (黑子)

7#
发表于 Post on 2022-7-20 09:20:22 | 只看该作者 Only view this author
本帖最后由 Jiongci-21 于 2022-7-20 09:24 编辑
linbalala 发表于 2022-7-18 21:16
同问请问有答案了嘛

我测试了两种情况,1、INCAR按照manual设置相应泛函的GGA参数,POTCAR用PAW_PBE原文件不变;2、INCAR设置GGA以及改变POTCAR中LEXCH为相应参数。
两者都能正常进行计算,最终体系的绝对能量会有差别,但是相对能量比较接近。前者stdout输出显示,在计算前会提醒你INCAR的GGA参数变了,后者的warning消失了。
综合manual的说法,我觉得应该是只用改变INCAR的GGA参数,那个warning只是提醒你用了新的修正。



202207200922218976..png (19.59 KB, 下载次数 Times of downloads: 42)

上排为所述1;下排位所述2

上排为所述1;下排位所述2

评分 Rate

参与人数
Participants 2
eV +9 收起 理由
Reason
sunjk + 4 谢谢分享
tanglei + 5

查看全部评分 View all ratings

越学越难顶。

9

帖子

0

威望

89

eV
积分
98

Level 2 能力者

8#
发表于 Post on 2022-7-22 10:37:05 | 只看该作者 Only view this author
Jiongci-21 发表于 2022-7-20 09:20
我测试了两种情况,1、INCAR按照manual设置相应泛函的GGA参数,POTCAR用PAW_PBE原文件不变;2、INCAR设置 ...

十分感谢您详细的分享,辛苦您了!感谢!

9

帖子

0

威望

89

eV
积分
98

Level 2 能力者

9#
发表于 Post on 2022-7-24 10:22:42 | 只看该作者 Only view this author
Jiongci-21 发表于 2022-3-26 10:45
您好,最近在尝试Cu(111)表面上吸附CO用不同泛函的测试。RPBE, revPBE和PBEsol是否只需要更改INCAR的GGA ...

有没有试过BEEFVDW,RPEB+D3。另外请问有无计算时间的测试呢

34

帖子

0

威望

1003

eV
积分
1037

Level 4 (黑子)

10#
发表于 Post on 2022-7-26 20:42:37 | 只看该作者 Only view this author
linbalala 发表于 2022-7-24 10:22
有没有试过BEEFVDW,RPEB+D3。另外请问有无计算时间的测试呢

BEEFVDW没有用,组里没有编译。上述的计算全部都加了D3校正的。年初测试的,时间数据可能得找找,我晚点看看吧。
越学越难顶。

32

帖子

0

威望

253

eV
积分
285

Level 3 能力者

11#
发表于 Post on 2022-10-14 17:20:43 | 只看该作者 Only view this author
hgyhgy 发表于 2022-5-14 22:16
但这个修改LEXCH的做法是错误的。

参考vasp wiki的说明,见https://www.vasp.at/wiki/index.php/GGA

老师您好,那是不是只要在INCAR中输入GGA=PS就能用PBEsol泛函算吗?

81

帖子

1

威望

220

eV
积分
321

Level 3 能力者

12#
发表于 Post on 2022-10-26 10:04:50 | 只看该作者 Only view this author
本帖最后由 hgyhgy 于 2022-10-26 10:05 编辑
qq870624836 发表于 2022-10-14 17:20
老师您好,那是不是只要在INCAR中输入GGA=PS就能用PBEsol泛函算吗?

对。不用修改POTCAR的。修改了,结果是错误的。

这个既然vasp wiki没说,就不可能需要额外修改。

这个修改了,会导致晶格常数更接近实验,并不能成为有效的论据。

32

帖子

0

威望

253

eV
积分
285

Level 3 能力者

13#
发表于 Post on 2022-10-26 12:08:43 | 只看该作者 Only view this author
hgyhgy 发表于 2022-10-26 10:04
对。不用修改POTCAR的。修改了,结果是错误的。

这个既然vasp wiki没说,就不可能需要额外修改。

谢谢您

1

帖子

0

威望

18

eV
积分
19

Level 1 能力者

14#
发表于 Post on 2024-9-27 17:07:58 | 只看该作者 Only view this author
hgyhgy 发表于 2022-10-26 10:04
对。不用修改POTCAR的。修改了,结果是错误的。

这个既然vasp wiki没说,就不可能需要额外修改。

老师,我发现只在INCAR中修改GGA=PS,不修改POTCAR算出来的晶格常数也比直接用PBE算出来的晶格常数更接近实验值;所以有两个问题不太懂:(1)晶格常数更接近实验是不是PBEsol泛函特有的呀(2)PBEsol除了能算出更精准的晶格常数外,在计算数据上还有什么特点,在自洽scf计算时要继续使用PBEsol泛函吗

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 06:16 , Processed in 0.193345 second(s), 25 queries , Gzip On.

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