计算化学公社

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

作者
Author: Eudaimonia 时间: 2023-11-20 18:01
标题: 在Multiwfn生成inp文件基础上进行基组快速设置的脚本(使用TZVP及以上基组时使用)
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等基组存在多大区别

欢迎光临计算化学公社 (http://bbs.keinsci.com/)