计算化学公社

标题: GGA=B3时，对应的POTCAR应该用哪个？ [打印本页]

作者
Author: SUN12 时间: 2023-7-27 16:40
标题: GGA=B3时，对应的POTCAR应该用哪个？
请问INCAR中的GGA选用B3LYP时，对应的POTCAR应该选择哪个啊，这几种的POTCAR中的LEXCH都不是B3，，，[attach]72332[/attach]

作者
Author: wypkdhd 时间: 2023-7-28 15:05
本帖最后由 wypkdhd 于 2023-7-28 15:08 编辑

B3 B3LYP[20] with VWN3[8] for LDA correlation
上方文字出自https://www.vasp.at/wiki/index.php/GGA
你应该用的赝势，为LDA。
参考https://www.vasp.at/wiki/index.php/Available_PAW_potentials

PS：如果你使用的vasp版本是vasp5.4或者vasp6，请使用potpaw_LDA.54.tar，赝势文件在http://bbs.keinsci.com/thread-16082-1-1.html的第四章。

作者
Author: SUN12 时间: 2023-7-28 22:41

wypkdhd 发表于 2023-7-28 15:05
B3 B3LYP[20] with VWN3[8] for LDA correlation
上方文字出自https://www.vasp.at/wiki/index.php ...

好嘞，非常感谢

作者
Author: 卡开发发 时间: 2023-7-29 10:34

wypkdhd 发表于 2023-7-28 15:05
B3 B3LYP[20] with VWN3[8] for LDA correlation
上方文字出自https://www.vasp.at/wiki/index.php ...

B3LYP with VWN3 for LDA correlation的意思是，B3打头的B3LYP的关联泛函部分使用了LDA-VWN3的形式，其中另外一个B5则是结合VWN5，但并不是说这个B3LYP适合配合LDA的赝势。而其中LDA的赝势其实对应的也不是VWN，而是CA-PZ的形式，当然，如果真的进行LDA计算影响不会很大，大部分的LDA形式结果差异不会很大。

LZ的情况我猜也许可以用potPAW_GGA（应该对应PW91）或者PBE，原理上说，应该选用各方面复现B3LYP全电子计算结果比较好的那一个，选用哪种其实目前应该没啥定论因为没有足够的数据作为支撑。

作者
Author: wypkdhd 时间: 2023-7-29 10:54

卡开发发发表于 2023-7-29 10:34
B3LYP with VWN3 for LDA correlation的意思是，B3打头的B3LYP的关联泛函部分使用了LDA-VWN3的形式，其中 ...

你说的对，我的说法有误。其实搞了这几年计算，我现在用的还是pbe，顶多有beefvdw。看文献也是，对于vasp，不管啥泛函，都是统一的vasp5.4 pbe赝势。

作者
Author: SUN12 时间: 2023-7-29 21:11

卡开发发发表于 2023-7-29 10:34
B3LYP with VWN3 for LDA correlation的意思是，B3打头的B3LYP的关联泛函部分使用了LDA-VWN3的形式，其中 ...

卡神又尝试了一下PW91的赝势,，但是计算基本上很快结束，  错误提示如下，这是为啥呢----------------------------------------------------------------------------
|                                                                            |
|  ADVICE TO THIS USER RUNNING 'VASP/VAMP' (HEAR YOUR MASTER'S VOICE ...):  |
|                                                                            |
|    You enforced a specific xc-type in the INCAR file,                   |
|    a different type was found on the POTCAR file                         |
|       I HOPE YOU KNOW, WHAT YOU ARE  DOING                            |
|                                                                            |
-----------------------------------------------------------------------------

LDA part: xc-table for VWN3       , standard interpolation
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...
WAVECAR not read
entering main loop
   N    E                   dE          d eps    ncg    rms       rms(c)

===================================================================================
= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
= PID 24495 RUNNING AT cu09
= EXIT CODE: 9
= CLEANING UP REMAINING PROCESSES
= YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================
Intel(R) MPI Library troubleshooting guide:
   https://software.intel.com/node/561764
===================================================================================

作者
Author: 卡开发发 时间: 2023-7-30 11:25
本帖最后由卡开发发于 2023-7-30 11:33 编辑

SUN12 发表于 2023-7-29 21:11
卡神又尝试了一下PW91的赝势,，但是计算基本上很快结束，错误提示如下，这是为啥呢------------------- ...

你实际的错误我不清楚，也许要去OUTCAR看看，也许是编译问题或者其他兼容性方面的问题（这点我暂时判断不了）。上面赝势的Warning与此无关，目前VASP官方提供的赝势应该也没有和GGA=B3所对应的。Warning的情况绝大多数能运行下去，但是通常都是有用的建议，赝势产生与泛函不匹配的时候其实确实该慎重处理。

作者
Author: SUN12 时间: 2023-7-30 12:24

卡开发发发表于 2023-7-30 11:25
你实际的错误我不清楚，也许要去OUTCAR看看，也许是编译问题或者其他兼容性方面的问题（这点我暂时判断不 ...

嗯嗯，谢谢卡神

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