|
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等基组存在多大区别
|
评分 Rate
-
查看全部评分 View all ratings
|