计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3960|回复 Reply: 15
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 使用multiwfn生成cp2k输入文件后,发现没有对应基组该如何计算呢?

[复制链接 Copy URL]

156

帖子

0

威望

3822

eV
积分
3978

Level 5 (御坂)

最近计算Yb元素合金,在multiwfn程序中生成的输入文件中,Yb的基组为 TZV2PX-MOLOPT-GTH-q34 但是在对应的BASIS_MOLOPT 或者BASIS_MOLOPT_LnPP2  、BASIS_MOLOPT_UCL 都没有发现有对应的基组,想问下,这种基组文件我该去哪里找呢?

此外,我发现基组的价电子是要与赝势匹配的,POTENTIALS赝势文件里面只有q24,要想得到q34 就只有选择 LnPP1_POTENTIALS,但这个关键词不能重复,LnPP1_POTENTIALS对于其他非镧系金属又不能描述,请问这个问题该如何解决呢?

最后,关于这类镧系金属元素,我是否需要加U呢,之前测试U值发现晶格参数越来越大,反而偏离了实验值(还不如不加U,但又怕审稿人问起),也没有找到对应的文献,不知该如何是好。

182

帖子

0

威望

1486

eV
积分
1668

Level 5 (御坂)

16#
发表于 Post on 3 day ago | 只看该作者 Only view this author
sobereva 发表于 2026-8-10 18:46
北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)培训里明确讲了:

感谢社长的耐心回复,我知道我的问题所在了
另外发现超算安装CP2K没写CP2K_DATA_DIR的环境变量

6万

帖子

99

威望

6万

eV
积分
128902

管理员

公社社长

15#
发表于 Post on 3 day ago | 只看该作者 Only view this author
杨小狗 发表于 2026-8-10 16:31
DFT模块的基组和赝势也必须给出全路径吗?

北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)培训里明确讲了:

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
128902

管理员

公社社长

14#
发表于 Post on 3 day ago | 只看该作者 Only view this author
杨小狗 发表于 2026-8-10 16:29
您好,关于您(2)中的回答,我用的是CP2K-2025.1版本想要在输入文件中指定Yb基组的文件,但我用您的方法 ...

北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)培训里明确讲了怎么确认自带的基组/赝势文件中支持某种元素的基组/赝势有哪些



之后基组名、赝势名、基组文件名、赝势文件名都指定对了自然就不会报错
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

182

帖子

0

威望

1486

eV
积分
1668

Level 5 (御坂)

13#
发表于 Post on 3 day ago | 只看该作者 Only view this author
杨小狗 发表于 2026-8-10 16:29
您好,关于您(2)中的回答,我用的是CP2K-2025.1版本想要在输入文件中指定Yb基组的文件,但我用您的方法 ...

DFT模块的基组和赝势也必须给出全路径吗?

182

帖子

0

威望

1486

eV
积分
1668

Level 5 (御坂)

12#
发表于 Post on 3 day ago | 只看该作者 Only view this author
本帖最后由 杨小狗 于 2026-8-10 16:30 编辑
乐平 发表于 2023-8-3 11:27
(1) CP2K 的基组和赝势都在安装目录的 XXXX/cp2k-yy/data 路径里,这里的 XXXX 是你自己的安装路径,yy 代 ...

您好,关于您(2)中的回答,我用的是CP2K-2025.1版本想要在输入文件中指定Yb基组的文件,但我用您的方法书写了inp文件,但输出的out文件一直报错找不到Yb基组定义
……
&END KIND
    &KIND Yb
      ELEMENT Yb
      BASIS_SET /public3/home/sc54133/POT/BASIS_MOLOPT_LnPP2/DZVP-MOLOPT-SR-GTH-q11
      POTENTIAL /public3/home/sc54133/POT/LnPP2_POTENTIALS
    &END KIND
    &KIND Y
      ELEMENT Y
      BASIS_SET DZVP-MOLOPT-SR-GTH-q11
      POTENTIAL GTH-PBE
    &END KIND
    &KIND Er
      ELEMENT Er
      BASIS_SET /public3/home/sc54133/POT/BASIS_MOLOPT_LnPP2/DZVP-MOLOPT-SR-GTH-q11
      POTENTIAL /public3/home/sc54133/POT/LnPP2_POTENTIALS
    &END KIND
    &KIND Cl
      ELEMENT Cl
      BASIS_SET DZVP-MOLOPT-SR-GTH-q7
      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
#   WFN_RESTART_FILE_NAME xitu_opt-RESTART.wfn

……
*******************************************************************************
*   ___                                                                       *
*  /   \                                                                      *
* [ABORT]                        The requested basis set                      *
*  \___/   </public3/home/sc54133/POT/BASIS_MOLOPT_LnPP2/DZVP-MOLOPT-SR-GTH-q *
*    |          11> for element <Yb> was not found in the basis set files     *
*  O/|                     <BASIS_MOLOPT><BASIS_MOLOPT_LnPP2>                 *
* /| |                                                                        *
* / \                                          aobasis/basis_set_types.F:1463 *
*******************************************************************************


===== Routine Calling Stack =====

            3 read_qs_kind
            2 create_qs_kind_set
            1 CP2K

……

35

帖子

1

威望

458

eV
积分
513

Level 4 (黑子)

11#
发表于 Post on 2023-8-19 01:33:05 | 只看该作者 Only view this author
ms860309 发表于 2023-8-4 00:17
感謝回覆

並沒有看到有特殊的符號 or 換行 如\n這種

更新一下
使用CP2K github上面新的BASIS_pob已經有包含rev2
且能正常使用
https://github.com/cp2k/cp2k/blob/master/data/BASIS_pob

35

帖子

1

威望

458

eV
积分
513

Level 4 (黑子)

10#
发表于 Post on 2023-8-4 00:17:46 | 只看该作者 Only view this author
乐平 发表于 2023-8-3 23:45
从报错的信息看,程序希望获得你的 pob-DZVP-rev2 文件中第三行,第一列是整数,但是发现非整数的字符。
...

感謝回覆

並沒有看到有特殊的符號 or 換行 如\n這種
cluster上沒有dos2unix,但是裝了win版轉換後丟上去狀況一樣
似乎不是這方面問題

1160

帖子

0

威望

4779

eV
积分
5939

Level 6 (一方通行)

9#
发表于 Post on 2023-8-3 23:45:40 | 只看该作者 Only view this author
ms860309 发表于 2023-8-3 23:25
借題詢問一下,我想使用PBE/pob-DZVP-rev2,如下做了定義
basis set是在這邊下載的 https://www.chemiebn. ...

从报错的信息看,程序希望获得你的 pob-DZVP-rev2 文件中第三行,第一列是整数,但是发现非整数的字符。

可能你下载的文件格式不是 Linux 格式?  换行符是 Windows 的?

你可以在命令行输入
  1. vi -b pob-DZVP-rev2
复制代码

查看一下文件,如果出现特殊符号,说明是格式问题。

如果的集群里有 dos2unix 工具,可以直接在命令行输入
  1. dos2unix pob-DZVP-rev2
复制代码

转换格式

35

帖子

1

威望

458

eV
积分
513

Level 4 (黑子)

8#
发表于 Post on 2023-8-3 23:25:43 | 只看该作者 Only view this author
本帖最后由 ms860309 于 2023-8-3 23:29 编辑

借題詢問一下,我想使用PBE/pob-DZVP-rev2,如下做了定義
basis set是在這邊下載的 https://www.chemiebn.uni-bonn.de ... /software/ssc_basis



input file則是作如下修改



但執行時卻會出現這個錯誤


想問問我哪邊弄錯了

156

帖子

0

威望

3822

eV
积分
3978

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2023-8-3 17:30:42 | 只看该作者 Only view this author
sobereva 发表于 2023-8-3 17:06
Multiwfn根据你选的基组的名字和元素自动产生基组名,实际基组文件里没定义的就自己改

好的,我以后多留意一下这个,谢谢社长~

156

帖子

0

威望

3822

eV
积分
3978

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2023-8-3 17:30:08 | 只看该作者 Only view this author
乐平 发表于 2023-8-3 16:11
(1) 自己确认一下是保险的,Multiwfn 有 bug 的话,跟社长反馈,他会很快更新。

(2) 看 CP2K 官方文档 ...

好的好的,谢谢您~

6万

帖子

99

威望

6万

eV
积分
128902

管理员

公社社长

5#
发表于 Post on 2023-8-3 17:06:01 | 只看该作者 Only view this author
1138711019 发表于 2023-8-3 13:28
谢谢大佬详细的回答,对于第一个问题,我采用的是BASIS_MOLOPT_UCL : Yb  TZV2P-MOLOPT-SR-GTH-q24,比较 ...

Multiwfn根据你选的基组的名字和元素自动产生基组名,实际基组文件里没定义的就自己改。

不能简单说加不加U,要看具体加多少,加0.001和没加一样。
如果没有适合的、专门用于你研究的问题对应的已知的Ueff参数就甭加,否则可能有害无益。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1160

帖子

0

威望

4779

eV
积分
5939

Level 6 (一方通行)

4#
发表于 Post on 2023-8-3 16:11:39 | 只看该作者 Only view this author
1138711019 发表于 2023-8-3 13:28
谢谢大佬详细的回答,对于第一个问题,我采用的是BASIS_MOLOPT_UCL : Yb  TZV2P-MOLOPT-SR-GTH-q24,比较 ...

(1) 自己确认一下是保险的,Multiwfn 有 bug 的话,跟社长反馈,他会很快更新。

(2) 看 CP2K 官方文档,https://manual.cp2k.org/cp2k-202 ... BASIS_SET_FILE_NAME

&DFT
&BASIS_SET_FILE_NAME 可以重复
&POTENTIAL_FILE_NAME 不能重复,需要仔细挑选

156

帖子

0

威望

3822

eV
积分
3978

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2023-8-3 13:28:01 | 只看该作者 Only view this author
乐平 发表于 2023-8-3 11:27
(1) CP2K 的基组和赝势都在安装目录的 XXXX/cp2k-yy/data 路径里,这里的 XXXX 是你自己的安装路径,yy 代 ...

谢谢大佬详细的回答,对于第一个问题,我采用的是BASIS_MOLOPT_UCL : Yb  TZV2P-MOLOPT-SR-GTH-q24,比较好奇multiwfn是怎么写出了 TZV2PX-MOLOPT-GTH-q34 的,我以为有,但自己又一直没找到对应的基组文件。

第二个问题就是我在元素部分定义了基组和赝势,还需要再&DFT部分再定义吗,二者会不会冲突呢

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

GMT+8, 2026-8-13 11:42 , Processed in 0.248165 second(s), 24 queries , Gzip On.

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