计算化学公社

标题: 用VASP计算U值，具体怎么算？ [打印本页]
作者
Author: 2572094680    时间: 2023-8-8 22:17
标题: 用VASP计算U值，具体怎么算？
请问一下大家有没有用VASP计算过U值啊，具体怎么算呀，找了一下官网上的没看太懂
(, 下载次数 Times of downloads: 8)

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

(, 下载次数 Times of downloads: 8)
作者
Author: 阿司匹林泡腾片    时间: 2023-8-9 09:24
我之前看过的一个帖子：点这
作者
Author: 1138711019    时间: 2023-8-9 10:24
sobereva 发表于 2023-8-8 22:25
我不清楚那个脚本，只能说如果有实验值可参考，自己弄Ueff值也不是什么费劲的事，参考北京科音CP2K第一性原 ...

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

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

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

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

(, 下载次数 Times of downloads: 9)
作者
Author: 1138711019    时间: 2023-8-10 10:29
sobereva 发表于 2023-8-10 08:56
CP2K有对镧系用的大核赝势，培训里我也提了。LnPP1那种中/小核赝势算不动或难收敛时可以试试大核。

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

  1.   &KIND Yb
  2.       ELEMENT Yb
  3.       BASIS_SET DZVP-MOLOPT-GTH-q11
  4.       POTENTIAL_FILE_NAME /home/XXX/Downloads/cp2k-2023.1/data/LnPP2_POTENTIALS
  5.             POTENTIAL GTH-PBE-q11
  6.     &END KIND
  7.     &KIND Al
  8.       ELEMENT Al
  9.       BASIS_SET DZVP-MOLOPT-SR-GTH-q3
  10.       POTENTIAL GTH-PBE
  11.     &END KIND
  12.     &KIND Li
  13.       ELEMENT Li
  14.       BASIS_SET DZVP-MOLOPT-SR-GTH-q3
  15.       POTENTIAL GTH-PBE
  16.     &END KIND
  17.   &END SUBSYS

  19.   &DFT
  20.     BASIS_SET_FILE_NAME  BASIS_MOLOPT
  21.     BASIS_SET_FILE_NAME  BASIS_MOLOPT_LnPP2
  22.     POTENTIAL_FILE_NAME  POTENTIAL
复制代码

不过测试LnPP1的时候却能够顺利运行，但是LnPP1中Yb的描述是GTH-PBE-q34，这个可比GTH-PBE-q11价电子更多，会不会运算量更大呢？我怕消耗资源更多。
  1.   &KIND Yb
  2.       ELEMENT Yb
  3.       BASIS_SET DZV-MOLOPT-SR-GTH
  4.       POTENTIAL_FILE_NAME /home/XXX/Downloads/cp2k-2023.1/data/LnPP1_POTENTIALS
  5.             POTENTIAL GTH-PBE-q34
  6.     &END KIND
  7.     &KIND Al
  8.       ELEMENT Al
  9.       BASIS_SET DZVP-MOLOPT-SR-GTH-q3
  10.       POTENTIAL GTH-PBE
  11.     &END KIND
  12.     &KIND Li
  13.       ELEMENT Li
  14.       BASIS_SET DZVP-MOLOPT-SR-GTH-q3
  15.       POTENTIAL GTH-PBE
  16.     &END KIND
  17.   &END SUBSYS

  19.   &DFT
  20.     BASIS_SET_FILE_NAME  BASIS_MOLOPT
  21.     BASIS_SET_FILE_NAME  BASIS_MOLOPT_LnPP1
  22.     POTENTIAL_FILE_NAME  POTENTIAL
复制代码
作者
Author: 1138711019    时间: 2023-8-10 13:58
1138711019 发表于 2023-8-10 10:29
感谢卢老师的建议，我按照你的说法，如下写了Yb的描述，但是计算报错，报错为Incompatible Atomic Occupa ...

找到报错Incompatible Atomic Occupations Detected  的原因了，
https://groups.google.com/g/cp2k ... p;utm_source=footer



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3