计算化学公社

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

[CP2K] 在Multiwfn生成inp文件基础上进行基组快速设置的脚本(使用TZVP及以上基组时使用)

[复制链接 Copy URL]

205

帖子

3

威望

1894

eV
积分
2159

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
Multiwfn目前生成的cp2k文件时,使用较大基组(TZVP及以上)时,对所有元素使用TZVP-MOLOPT-GTH-qX基组(X为非赝势电子数),但是很多元素无该类型基组的定义,但是存在如TZVP-MOLOPT-PBE-GTH-qX类的基组,需要自行手动更改元素基组及Basis_File字段,在有大量inp文件时非常麻烦,该脚本能够帮助你快速修改所有元素的基组

在cp2k安装目录下的data文件夹内使用 grep "X" * > X.data获取基组索引文件用于该脚本,如grep "TZVP" * > TZVP.data,将TZVP.data文件拷贝至与该脚本的同一个目录下即可使用,之后脚本对inp文件进行更改时会读取TZVP.data文件,会根据inp文件中的第一个元素基组的第一个 “-” 之前的信息来决定读取哪个档次基组的data文件,大概率在读取 pob系列 等不把基组大小放在第一个“-”之前的基组时会出现bug,可以自行修改

双击.py文件以后根据提示自己选择基组即可使用,给两个data文件作为参考,是我在2023.1版本cp2k data目录下grep得到的
个人Python版本是3.7 64bit,可能需要注意

Note:
目前版本下只能把所有元素的基组改为同一类基组
目前脚本问题很多,暂时没工夫好好写一遍,但是对我的需求基本够用可能短时间内不太会进行修改了
Multiwfn默认设置下对于轻元素使用的是TZVP-MOLOPT-GTH-qX,暂时不确定与TZVP-MOLOPT-PBE-GTH-qX等基组存在多大区别

BasisSetForCp2k-v0.1.py

5.35 KB, 下载次数 Times of downloads: 10

TZV2P.data

23.32 KB, 下载次数 Times of downloads: 2

TZVP.data

33.36 KB, 下载次数 Times of downloads: 2

评分 Rate

参与人数
Participants 2
eV +12 收起 理由
Reason
hhqywj + 4 好物!
sobereva + 8

查看全部评分 View all ratings

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

GMT+8, 2024-11-25 17:40 , Processed in 0.248570 second(s), 25 queries , Gzip On.

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