捐赠本论坛 (Donate this site)
开启辅助访问

计算化学公社

 找回密码 Forget password
 注册 Register
在11月21-24日期间，本论坛仅限等级≥level 4的成员发贴，带来的不便请谅解！其他成员若有学术问题请等到能发贴的时候再发帖。此期间不要给管理员sobereva发私信问学术问题
计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First Principle) › 用VASP计算U值，具体怎么算？

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

查看数: 1105 | 评论数: 8 | 收藏 Add to favorites 2
关灯 | 提示：支持键盘翻页<-左 右->
帖子模式
    组图打开中，请稍候......
2572094680
发布时间: 2023-8-8 22:17

正文摘要:

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

回复 Reply

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
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基组元素混用吗？

  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
复制代码
sobereva 发表于 Post on 2023-8-10 08:56:35
1138711019 发表于 2023-8-10 08:47
尝试了PBEsol，结果也不是很好，偏离实验值太多，不知道是不是我其他设置有问题，
提升基组的影响确实不 ...

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


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

借楼问下卢老师，我计算稀土体系时，用不加U的PBE就比实验值晶格常数高，然后再加U值就更高了，这种情况该如何处理呢？软件用的cp2k，而且改变基组从DZVP到TZVP或者更高精度的时候，晶格常数也在增加，不知道这样合理不，我该用哪种基组呢？
阿司匹林泡腾片 发表于 Post on 2023-8-9 09:24:36
我之前看过的一个帖子：点这
sobereva 发表于 Post on 2023-8-8 22:25:52
我不清楚那个脚本，只能说如果有实验值可参考，自己弄Ueff值也不是什么费劲的事，参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里讲DFT+U部分的这个典型文献里的做法：


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

GMT+8, 2024-11-23 23:07 , Processed in 0.263841 second(s), 26 queries , Gzip On.

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