计算化学公社

标题: vasp中如何使用PBEsol泛函？ [打印本页]

作者
Author: zhongyuabc 时间: 2018-12-17 13:36
标题: vasp中如何使用PBEsol泛函？
在vasp中使用PBEsol泛函时，除了在INCAR文件中输入GGA=PS外，是否还要把POTCAR文件中的LEXCH=91改为LEXCH=PS(我的POTCAR文件是从vasp赝势中的PAW_GGA目录中拷贝的），谢谢解答！

作者
Author: hakuna 时间: 2018-12-17 18:54
PBE、RPBE、PBEsol都需要用PBE版的POTCAR，不可以用PW91的
用PBEsol，两个都要改，否则会有警告，说你强迫了VASP，输出的能量也差别很大
两种设置，其实你可以测试一下，再读一下OUTCAR，对比一下看

作者
Author: zhongyuabc 时间: 2018-12-17 19:30
本帖最后由 zhongyuabc 于 2018-12-17 19:37 编辑

hakuna 发表于 2018-12-17 18:54
PBE、RPBE、PBEsol都需要用PBE版的POTCAR，不可以用PW91的
用PBEsol，两个都要改，否则会有警告，说你强迫 ...

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

作者
Author: Jiongci-21 时间: 2022-3-26 10:45
您好，最近在尝试Cu（111）表面上吸附CO用不同泛函的测试。RPBE, revPBE和PBEsol是否只需要更改INCAR的GGA参数，POTCAR均用PAW_PBE的赝势即可呢？

作者
Author: hgyhgy 时间: 2022-5-14 22:16

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应该会明确说明的。没有说明，证明是不用修改。而且也特别指出不能修改。

作者
Author: linbalala 时间: 2022-7-18 21:16

Jiongci-21 发表于 2022-3-26 10:45
您好，最近在尝试Cu（111）表面上吸附CO用不同泛函的测试。RPBE, revPBE和PBEsol是否只需要更改INCAR的GGA ...

同问请问有答案了嘛

作者
Author: Jiongci-21 时间: 2022-7-20 09:20
本帖最后由 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只是提醒你用了新的修正。

作者
Author: linbalala 时间: 2022-7-22 10:37

Jiongci-21 发表于 2022-7-20 09:20
我测试了两种情况，1、INCAR按照manual设置相应泛函的GGA参数，POTCAR用PAW_PBE原文件不变；2、INCAR设置 ...

十分感谢您详细的分享，辛苦您了！感谢！

作者
Author: linbalala 时间: 2022-7-24 10:22

Jiongci-21 发表于 2022-3-26 10:45
您好，最近在尝试Cu（111）表面上吸附CO用不同泛函的测试。RPBE, revPBE和PBEsol是否只需要更改INCAR的GGA ...

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

作者
Author: Jiongci-21 时间: 2022-7-26 20:42

linbalala 发表于 2022-7-24 10:22
有没有试过BEEFVDW，RPEB+D3。另外请问有无计算时间的测试呢

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

作者
Author: qq870624836 时间: 2022-10-14 17:20

hgyhgy 发表于 2022-5-14 22:16
但这个修改LEXCH的做法是错误的。

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

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

作者
Author: hgyhgy 时间: 2022-10-26 10:04
本帖最后由 hgyhgy 于 2022-10-26 10:05 编辑

qq870624836 发表于 2022-10-14 17:20
老师您好，那是不是只要在INCAR中输入GGA=PS就能用PBEsol泛函算吗？

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

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

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

作者
Author: qq870624836 时间: 2022-10-26 12:08

hgyhgy 发表于 2022-10-26 10:04
对。不用修改POTCAR的。修改了，结果是错误的。

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

谢谢您

作者
Author: ZYLin 时间: 2024-9-27 17:07

hgyhgy 发表于 2022-10-26 10:04
对。不用修改POTCAR的。修改了，结果是错误的。

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

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

作者
Author: 星星太 时间: 2025-4-3 15:36

Jiongci-21 发表于 2022-7-20 09:20
我测试了两种情况，1、INCAR按照manual设置相应泛函的GGA参数，POTCAR用PAW_PBE原文件不变；2、INCAR设置 ...

同学你好，请问pbesol的赝势库如何下载呢

作者
Author: tangfengzhang 时间: 2025-5-17 14:45

星星太发表于 2025-4-3 15:36
同学你好，请问pbesol的赝势库如何下载呢

同问，怎么下载pbesol赝势

欢迎光临计算化学公社 (http://ccc.keinsci.com/)