用VASP计算U值，具体怎么算？

2572094680 · 发表于 Post on 2023-8-8 22:17:04

请问一下大家有没有用VASP计算过U值啊，具体怎么算呀，找了一下官网上的没看太懂

sobereva · 发表于 Post on 2023-8-8 22:25:52

我不清楚那个脚本，只能说如果有实验值可参考，自己弄Ueff值也不是什么费劲的事，参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里讲DFT+U部分的这个典型文献里的做法：

阿司匹林泡腾片 · 发表于 Post on 2023-8-9 09:24:36

我之前看过的一个帖子：点这

1138711019 · 发表于 Post on 2023-8-9 10:24:19

sobereva 发表于 2023-8-8 22:25
我不清楚那个脚本，只能说如果有实验值可参考，自己弄Ueff值也不是什么费劲的事，参考北京科音CP2K第一性原 ...

借楼问下卢老师，我计算稀土体系时，用不加U的PBE就比实验值晶格常数高，然后再加U值就更高了，这种情况该如何处理呢？软件用的cp2k，而且改变基组从DZVP到TZVP或者更高精度的时候，晶格常数也在增加，不知道这样合理不，我该用哪种基组呢？

sobereva · 发表于 Post on 2023-8-9 20:32:36

1138711019 发表于 2023-8-9 10:24
借楼问下卢老师，我计算稀土体系时，用不加U的PBE就比实验值晶格常数高，然后再加U值就更高了，这种情况 ...

尝试其它泛函，比如PBEsol等
既然加了有害（而且无论Ueff的大小都令结果变差），显然就不该加。
要么就尝试杂化泛函。
另外考虑是否其它设置不当导致结果烂，比如k点考虑不当、磁性设置不当等
一般来说DZVP-MOLOPT-SR-GTH档次就够用了，提升基组一般不至于有多大影响。

1138711019 · 发表于 Post on 2023-8-10 08:47:21

sobereva 发表于 2023-8-9 20:32
尝试其它泛函，比如PBEsol等
既然加了有害（而且无论Ueff的大小都令结果变差），显然就不该加。
要么就 ...

尝试了PBEsol，结果也不是很好，偏离实验值太多，不知道是不是我其他设置有问题，
提升基组的影响确实不大，所以最后选择DZVP-MOLOPT-SR-GTH，没加U。
另外想问一下，cp2k是否能像vasp一样，将稀土元素的部分f电子当成内壳层处理（比如Yb_2和Yb_3），因为后面体系建了表面还扩胞了，cp2k计算老不收敛，所以想着能否像vasp这样简化一下有利于计算。
谢谢卢老师~

sobereva · 发表于 Post on 2023-8-10 08:56:35

1138711019 发表于 2023-8-10 08:47
尝试了PBEsol，结果也不是很好，偏离实验值太多，不知道是不是我其他设置有问题，
提升基组的影响确实不 ...

CP2K有对镧系用的大核赝势，培训里我也提了。LnPP1那种中/小核赝势算不动或难收敛时可以试试大核。

1138711019 · 发表于 Post on 2023-8-10 10:29:22

sobereva 发表于 2023-8-10 08:56
CP2K有对镧系用的大核赝势，培训里我也提了。LnPP1那种中/小核赝势算不动或难收敛时可以试试大核。

感谢卢老师的建议，我按照你的说法，如下写了Yb的描述，但是计算报错，报错为Incompatible Atomic Occupations Detected qs_kind_types.F:3405。
搜索一下，说是原子占据不兼容，但是Yb元素对应的基组赝势都是对的啊，难道说SR基组的元素不能和非SR基组元素混用吗？

&KIND Yb
ELEMENT Yb
BASIS_SET DZVP-MOLOPT-GTH-q11
POTENTIAL_FILE_NAME /home/XXX/Downloads/cp2k-2023.1/data/LnPP2_POTENTIALS
POTENTIAL GTH-PBE-q11
&END KIND
&KIND Al
ELEMENT Al
BASIS_SET DZVP-MOLOPT-SR-GTH-q3
POTENTIAL GTH-PBE
&END KIND
&KIND Li
ELEMENT Li
BASIS_SET DZVP-MOLOPT-SR-GTH-q3
POTENTIAL GTH-PBE
&END KIND
&END SUBSYS
&DFT
BASIS_SET_FILE_NAME BASIS_MOLOPT
BASIS_SET_FILE_NAME BASIS_MOLOPT_LnPP2
POTENTIAL_FILE_NAME POTENTIAL

复制代码

不过测试LnPP1的时候却能够顺利运行，但是LnPP1中Yb的描述是GTH-PBE-q34，这个可比GTH-PBE-q11价电子更多，会不会运算量更大呢？我怕消耗资源更多。

&KIND Yb
ELEMENT Yb
BASIS_SET DZV-MOLOPT-SR-GTH
POTENTIAL_FILE_NAME /home/XXX/Downloads/cp2k-2023.1/data/LnPP1_POTENTIALS
POTENTIAL GTH-PBE-q34
&END KIND
&KIND Al
ELEMENT Al
BASIS_SET DZVP-MOLOPT-SR-GTH-q3
POTENTIAL GTH-PBE
&END KIND
&KIND Li
ELEMENT Li
BASIS_SET DZVP-MOLOPT-SR-GTH-q3
POTENTIAL GTH-PBE
&END KIND
&END SUBSYS
&DFT
BASIS_SET_FILE_NAME BASIS_MOLOPT
BASIS_SET_FILE_NAME BASIS_MOLOPT_LnPP1
POTENTIAL_FILE_NAME POTENTIAL

复制代码

1138711019 · 发表于 Post on 2023-8-10 13:58:13

1138711019 发表于 2023-8-10 10:29
感谢卢老师的建议，我按照你的说法，如下写了Yb的描述，但是计算报错，报错为Incompatible Atomic Occupa ...

找到报错Incompatible Atomic Occupations Detected 的原因了，

https://groups.google.com/g/cp2k ... p;utm_source=footer

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] 用VASP计算U值，具体怎么算？