计算化学公社

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

[CP2K] 读赝势文件给cp2k molden file添加Cell和Nval section

[复制链接 Copy URL]

14

帖子

1

威望

144

eV
积分
178

Level 3 能力者

本帖最后由 ykhuang 于 2023-3-14 19:07 编辑

读你本地的cp2k/data/文件夹下赝势文件,之后给cp2k输出的molden文件加[Cell]、[Nval] section,然后顺便还可以调起Multiwfn,当然了你得设置好...(啊突然想起来我一开始只是为了想要做一下双正交化的...又写小程序玩了一下午
调用方式有很多,你可以
1. 在IDE里面打开主动模块,设置moldenFileName = [你的molden文件(含路径)],也可以
2. 命令行里面> python cp2kMoldenPostProcessingToMultiwfn.py -i [你的molden文件(含路径)](更多的选项请“--help”)。
程序会读赝势文件然后提取出来价电子处理的电子数量,存成字典,导出json,这样下回就可以再用了。

注意哈,斜杠同时还是个转义符号,所以得转义掉,就像pseudoPotentialPath最后两个“\\”一样(啰嗦给不写程序的人)

# Section: pseudopotentials I/O
pseudoPotentialFile = 'GTH_POTENTIALS' # only GTH_POTENTIALS is supported, for supporting more kinds of pseudopotential files,
                                       # modify _cp2kImportNumOfValenceElectrons.py
pseudoPotentialPath = 'D:\cp2k-2023.1\data\\'
functional = 'PBE'
boolSmallCorePseudoPotential = True
pseudoPotentialDatabaseGenerator = '_cp2kImportNumOfValenceElectrons.py' # what contents you must not change
pPDGPath = 'E:\Multiwfn_3.8_dev_bin_Win64\\' # absolute path of _cp2kImportNumOfValenceElectrons.py file
# Section: molden file I/O
moldenFileName = 'F:\***\***\***.molden'
# should include absolute path of your molden file
# however if you also supply another molden file by stdin, then that will be read instead
outMoldenPath = '.\\'
# Section: Multiwfn
multiwfnPath = 'E:\Multiwfn_3.8_dev_bin_Win64\\'
boolCallMultiwfn = True


PseudoPotentialFile:读取的赝势文件名
pseudoPotentialPath:赝势文件文件夹路径,推荐是写绝对路径
functional:你DFT用的泛函
boolSmallCorePseudoPotential:考虑到有的元素的赝势分小核和大核,True和False自行选择
pseudoPotentialDatabaseGenerator:读取赝势文件的那个程序,名字没事的话一般不用动,里面内容要动的话一般也是为了支持更多格式的赝势文件
驼峰命名法的好处就是可以肆意在后文缩写,然后缩写pseudoPotentialDatabaseGenerator -> pPDG,
pPDGPath:读赝势文件程序的文件夹路径,同样推荐写绝对路径
moldenFileName:你的molden文件路径,如果你还在命令行调用的时候又提供了新的,那就会默认读那个新的,即设置stdin优先级更高
outMoldenPath:添加[Nval]信息之后新的molden文件放在哪里
multiwfnPath:你安装Multiwfn的地方
boolCallMultiwfn:决定要不要唤起Multiwfn

展示几种调用方式:
1. 目录下没有json文件&以设置moldenFileName变量方式传入molden&不唤起Multiwfn

2. (刚刚跑了上一run所以)目录下有了json文件&stdin方式传入molden&不唤起Multiwfn

3. 目录下有了json文件&stdin方式传入molden&唤起Multiwfn

那就这样,886

_cp2kImportNumOfValenceElectrons.py

4.65 KB, 下载次数 Times of downloads: 8

被动模块,正常运行的必需文件

cp2kMoldenPostProcessingToMultiwfn.py

9.66 KB, 下载次数 Times of downloads: 12

主动模块,可以在IDE里面运行,也可以使用命令行

评分 Rate

参与人数
Participants 2
威望 +1 eV +4 收起 理由
Reason
JamesBourbon + 4 好物!
sobereva + 1

查看全部评分 View all ratings

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 00:25 , Processed in 0.194364 second(s), 26 queries , Gzip On.

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