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

计算化学公社

 找回密码 Forget password
 注册 Register
计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 请教关于CP2K找不到势文件的解决方法

请教关于CP2K找不到势文件的解决方法

查看数: 4890 | 评论数: 10 | 收藏 Add to favorites 6
关灯 | 提示：支持键盘翻页<-左 右->
帖子模式
    组图打开中，请稍候......
白皛白
发布时间: 2022-11-28 22:49

正文摘要:

老师们好 前几天刚刚接触CP2K，想要尝试用CP2K计算一个含有Ce Cu O三种原子的体系，使用Multiwfn生成inp文件后发现出现了以下的报错信息 The requested basis set <DZVP-MOLOPT-SR-GTH-Q12> for element <CE ...

回复 Reply

honest 发表于 Post on 2023-2-9 09:06:25
reid 发表于 2023-2-8 19:43
例如：
&DFT
    BASIS_SET_FILE_NAME  /home/DFT/cp2k-2022.1/data/BASIS_MOLOPT_UZH

谢谢reid的指导，非常有用
honest 发表于 Post on 2023-2-9 09:05:35
Eudaimonia 发表于 2023-2-8 17:50
BASIS_SET_FILE_NAME  BASIS_MOLOPT
    BASIS_SET_FILE_NAME  BASIS_MOLOPT_UCL

这个是真没想到，多谢您的回复
reid 发表于 Post on 2023-2-8 19:43:30
honest 发表于 2023-2-8 10:42
卢老师，非常感谢您的解答。我在BASIS_SET_FILE_NAME后面加了Gd元素的基组文件路径后（/home/zheng/softw ...

例如：
&DFT
    BASIS_SET_FILE_NAME  /home/DFT/cp2k-2022.1/data/BASIS_MOLOPT_UZH
    BASIS_SET_FILE_NAME  /home/DFT/cp2k-2022.1/data/BASIS_MOLOPT
    POTENTIAL_FILE_NAME  /home/DFT/cp2k-2022.1/data/POTENTIAL_UZH

评分 Rate

参与人数
Participants 1		eV +2 收起 理由
Reason
honest + 2 谢谢

查看全部评分 View all ratings

reid 发表于 Post on 2023-2-8 19:40:20
honest 发表于 2023-2-8 10:42
卢老师，非常感谢您的解答。我在BASIS_SET_FILE_NAME后面加了Gd元素的基组文件路径后（/home/zheng/softw ...

BASIS SET可以多写几行，你用的那个可能含有Gd，但是没有FCHO，另起一行，把含FCHO的也写上
Eudaimonia 发表于 Post on 2023-2-8 17:50:19
honest 发表于 2023-2-8 10:42
卢老师，非常感谢您的解答。我在BASIS_SET_FILE_NAME后面加了Gd元素的基组文件路径后（/home/zheng/softw ...

    BASIS_SET_FILE_NAME  BASIS_MOLOPT
    BASIS_SET_FILE_NAME  BASIS_MOLOPT_UCL

没想到吧，这个关键字可以写好几个，在data文件夹下面的应该不需要把路径写完去

评分 Rate

参与人数
Participants 2		eV +7 收起 理由
Reason
HanC + 5 谢谢
honest + 2 谢谢

查看全部评分 View all ratings

honest 发表于 Post on 2023-2-8 10:42:41
sobereva 发表于 2023-2-8 03:11
分清楚基组名和基组定义文件
BASIS SET设的是基组文件里定义的基组名，不是让你在这里设基组文件路径
...

卢老师，非常感谢您的解答。我在BASIS_SET_FILE_NAME后面加了Gd元素的基组文件路径后（/home/zheng/software/cp2k-9.1/data/BASIS_MOLOPT_UCL），Gd元素的基组没有报错。但是output文件又显示其他元素（F，C，H，O）的基组文件找不到了，请问又该怎么设置呢？
我试了以下几种方式，都没有解决:
① /home/zheng/software/cp2k-9.1/data/BASIS_MOLOPT_UCL & BASIS_MOLOPT
②/home/zheng/software/cp2k-9.1/data/BASIS_MOLOPT_UCL & /home/zheng/software/cp2k-9.1/data/BASIS_MOLOPT
sobereva 发表于 Post on 2023-2-8 03:11:23
honest 发表于 2023-2-7 22:45
老师，您好！同样是CP2K的新人，计算一个MOF的能量时，遇到了基组中找不到Gd元素的问题。参考您给出的建 ...

分清楚基组名和基组定义文件
BASIS_SET设的是基组文件里定义的基组名，不是让你在这里设基组文件路径
基组文件路径是在BASIS_SET_FILE_NAME里设的
honest 发表于 Post on 2023-2-7 22:45:10
本帖最后由 honest 于 2023-2-7 22:47 编辑
乐平 发表于 2022-11-28 23:01
emmmmm……
如果你在 data/ 目录下 grep Ce *，应该能看到含有 Ce 的基组不止 BASIS_MOLOPT_LnPP1，还有以 ...

老师，您好！同样是CP2K的新人，计算一个MOF的能量时，遇到了基组中找不到Gd元素的问题。参考您给出的建议，我在后面加入了路径，但是依然显示找不到，该怎么处理呢？

inp文件中输入设置

output文件中提醒错误

Gd元素所在基组的路径查看

**本来基组文件是data文件下的，但是我按照路径输入后又显示路径大于80个字节，所以我将该文件copy到当前目前下，然后再把路径写进去的。

期待您的回复，非常感谢！


白皛白 发表于 Post on 2022-11-29 12:30:10
乐平 发表于 2022-11-28 23:01
emmmmm……
如果你在 data/ 目录下 grep Ce *，应该能看到含有 Ce 的基组不止 BASIS_MOLOPT_LnPP1，还有以 ...

我按照您说的方法尝试了一下，问题已经迎刃而解了，非常感谢您的指教。
乐平 发表于 Post on 2022-11-28 23:01:40
本帖最后由 乐平 于 2022-11-28 23:23 编辑

emmmmm……
如果你在 data/ 目录下 grep Ce *，应该能看到含有 Ce 的基组不止 BASIS_MOLOPT_LnPP1，还有以下内容
  1. [huan@master data]$ grep Ce *
  2. BASIS_ccGRB_UZH:# Cesium
  3. BASIS_MOLOPT_LnPP1:Ce DZV-MOLOPT-SR-GTH
  4. BASIS_MOLOPT_UCL:Ce SZV-MOLOPT-SR-GTH SZV-MOLOPT-SR-GTH-q12
  5. BASIS_MOLOPT_UCL:Ce DZVP-MOLOPT-SR-GTH DZVP-MOLOPT-SR-GTH-q12
  6. BASIS_MOLOPT_UCL:Ce TZVP-MOLOPT-SR-GTH TZVP-MOLOPT-SR-GTH-q12
  7. BASIS_MOLOPT_UCL:Ce TZV2P-MOLOPT-SR-GTH TZV2P-MOLOPT-SR-GTH-q12
  8. BASIS_MOLOPT_UCL:Ce SZV-MOLOPT-SR-GTH-q30
  9. BASIS_MOLOPT_UCL:Ce DZV-MOLOPT-SR-GTH-q30
  10. BASIS_MOLOPT_UCL:Ce DZVP-MOLOPT-SR-GTH-q30
  11. BASIS_MOLOPT_UCL:Ce TZVP-MOLOPT-SR-GTH-q30
  12. BASIS_MOLOPT_UCL:Ce TZV2P-MOLOPT-SR-GTH-q30
  13. BASIS_MOLOPT_UZH:# Cesium
  14. BASIS_MOLOPT_UZH:# Cesium
  15. grep: DFTB: Is a directory
  16. ECP_POTENTIALS:Ce nelec 28
  17. ECP_POTENTIALS:Ce ul
  18. ECP_POTENTIALS:Ce S
  19. ECP_POTENTIALS:Ce P
  20. ECP_POTENTIALS:Ce D
  21. ECP_POTENTIALS:Ce F
  22. ECP_POTENTIALS:Ce G
  23. ECP_POTENTIALS:# Ce - Lu: R.B. Ross, W.C. Ermler, S. Das, To be published
  24. ECP_POTENTIALS:Ce nelec 54
  25. ECP_POTENTIALS:Ce ul
  26. ECP_POTENTIALS:Ce S
  27. ECP_POTENTIALS:Ce P
  28. ECP_POTENTIALS:Ce D
  29. GTH_POTENTIALS:Ce GTH-BLYP-q12 GTH-BLYP
  30. GTH_POTENTIALS:Ce GTH-PADE-q12 GTH-LDA-q12 GTH-PADE GTH-LDA
  31. GTH_POTENTIALS:Ce GTH-PBE-q12 GTH-PBE
  32. GTH_POTENTIALS:Ce GTH-PBE-q30
  33. LnPP1_POTENTIALS:Ce GTH-PBE-q12
  34. grep: NNP: Is a directory
  35. POTENTIAL:Ce GTH-BLYP-q12 GTH-BLYP
  36. POTENTIAL:Ce GTH-PADE-q12 GTH-LDA-q12 GTH-PADE GTH-LDA
  37. POTENTIAL:Ce GTH-PBE-q12 GTH-PBE
  38. POTENTIAL:Ce GTH-PBE-q30
  39. POTENTIAL_UZH:#  58  Ce        54      4
  40. POTENTIAL_UZH:#  58  Ce  4    NA
  41. POTENTIAL_UZH:#  58  Ce 12    0.006572       -0.000986       -0.005229
  42. POTENTIAL_UZH:#  58  Ce 30    0.002963       -0.002831       -0.001768
  43. POTENTIAL_UZH:Ce GTH-PBE-q4 GTH-GGA-q4
  44. POTENTIAL_UZH:Ce GTH-PBE-q12 GTH-GGA-q12
  45. POTENTIAL_UZH:Ce GTH-PBE-q30 GTH-GGA-q30
  46. POTENTIAL_UZH:#  58  Ce        54      4
  47. POTENTIAL_UZH:#  58  Ce  4    NA
  48. POTENTIAL_UZH:#  58  Ce 12   -0.000363       -0.000431       -0.010483
  49. POTENTIAL_UZH:#  58  Ce 30   -0.002957       -0.004674       -0.008319
  50. POTENTIAL_UZH:Ce GTH-SCAN-q4 GTH-MGGA-q4
  51. POTENTIAL_UZH:Ce GTH-SCAN-q12 GTH-MGGA-q12
  52. POTENTIAL_UZH:Ce GTH-SCAN-q30 GTH-MGGA-q30
  53. POTENTIAL_UZH:#  58  Ce        54      4
  54. POTENTIAL_UZH:#  58  Ce  4    NA
  55. POTENTIAL_UZH:#  58  Ce 12   -0.002690        0.000477       -0.002189
  56. POTENTIAL_UZH:#  58  Ce 30   -0.004002       -0.003255       -0.001825
  57. POTENTIAL_UZH:Ce GTH-PBE0-q4 GTH-HYB-q4
  58. POTENTIAL_UZH:Ce GTH-PBE0-q12 GTH-HYB-q12
  59. POTENTIAL_UZH:Ce GTH-PBE0-q30 GTH-HYB-q30
  60. grep: xc_section: Is a directory
  61. xTB_parameters:                    Ce       0.350000   1.200000  0.839861  30.638143   5d   -36.440945 -5.245538   -5.032003 3.000000
  62. [huan@master data][        DISCUZ_CODE_82        ]nbsp;
复制代码

至少BASIS_MOLOPT_LnPP1， BASIS_MOLOPT_UCL， GTH_POTENTIALS， POTENTIAL， POTENTIAL_UZH 里有 Ce 元素的基组。



另外，你想调用上述基组的话，在 BASIS_SET 关键词后面把基组的完整的路径写上就行了。比如
  1. &KIND
  2.    ELEMENT Ce
  3.    BASIS_SET /YOUR_DIR/cp2k/data/BASIS_MOLOPT_LnPP1/DZV-MOLOPT-SR-GTH
  4.    POTENTIAL GTH-PBE
复制代码

不需要你自己复制 BASIS_MOLOPT_LnPP1 里的 Ce 信息然后粘贴到 BASIS_MOLOPT 里……




评分 Rate

参与人数
Participants 2		eV +8 收起 理由
Reason
stc980309 + 4 谢谢
白皛白 + 4 谢谢

查看全部评分 View all ratings

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

GMT+8, 2025-8-14 01:57 , Processed in 0.194979 second(s), 26 queries , Gzip On.

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