计算化学公社

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

[CP2K] 关于cp2k生成molden文件及Multiwfn后处理作图的一些请教

[复制链接 Copy URL]

113

帖子

0

威望

1548

eV
积分
1661

Level 5 (御坂)

本帖最后由 nusiew 于 2022-7-9 15:36 编辑

1. 我看之前论坛上说 “cp2k生成molden文件需要用对角化而不要用OT的,否则里面没轨道能级信息”,我用OT方法和对角化方法分别输出的molden文件对比发现,二者形式上似乎没有区别,请问molden轨道能级信息从哪里可以看到,OT方法输出的molden文件是否缺少轨道能级信息?两种方法输出的molden文件放蓝奏云盘里了(https://wwd.lanzouy.com/iESbF07lk64j

2. 请问cp2k输出的molden文件是否是标准的molden文件?我参考Multiwfn手册附录5部分基于基函数和GTF分别进行了check,信息如下图。手册中写的是 If both the two maximum deviations are evidently larger than zero, then there must be some severe problems in the input file; if any one of them is very close to zero, the input file should be fully compatible with Multiwfn. 测试结果大概在0.5左右,这个结果是否意味着molden文件可用?

3. 如下,请问用Multiwfn绘制的ELF图的横纵坐标的范围是否可以修改?由于输出的图的横纵坐标值很模糊,请问大家是怎么解决这个问题的?能否把数据导出来绘图呢?谢谢大家!








两个对易的厄米算子可以有共同本征函数集

464

帖子

11

威望

3948

eV
积分
4632

Level 6 (一方通行)

2#
发表于 Post on 2022-7-9 20:40:50 | 只看该作者 Only view this author
1. molden格式文件里面有个Ene字段专门记录了能级信息,Google搜molden format
2. 是
3. 把数据导出,放到Origin里作图
自由发挥,野蛮生长

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

3#
发表于 Post on 2022-7-10 05:09:07 | 只看该作者 Only view this author
1 搜Ene=看后面的数值

2 怎么用CP2K产生能给Multiwfn做波函数分析用的molden文件在下文写明了,手册2.9.2.1节也写明了
使用Multiwfn结合CP2K通过NCI和IGM方法图形化考察固体和表面的弱相互作用
http://sobereva.com/588http://bbs.keinsci.com/thread-21742-1-1.html
严格照着做就不可能有问题

3 绝对不要直接截图,要在菜单里选择保存图片
并且搞清楚Multiwfn怎么控制输出的图片尺寸,看
谈谈与计算化学有关的作图的图像清晰问题
http://sobereva.com/607http://bbs.keinsci.com/thread-24177-1-1.html

保存图片得到的图一点都不模糊



北京科音自然科学研究中心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

113

帖子

0

威望

1548

eV
积分
1661

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2022-7-10 10:24:46 | 只看该作者 Only view this author
本帖最后由 nusiew 于 2022-7-10 10:41 编辑

谢谢sobereva老师和丁越老师的解答,对于cp2k默认输出的molden文件修改的地方总结为两个,一个是添加晶格常数信息,还有一个是如果使用了赝势,要将element index替换成实际的价电子数,是不是直接使用赝势中的被赝势赝化的电子数就可以了?比如说对于V,直接将原子序数23替换成13?
        &KIND V
            BASIS_SET DZVP-MOLOPT-SR-GTH
            POTENTIAL GTH-PBE-q13
        &END KIND


202207101016111500..png (159.21 KB, 下载次数 Times of downloads: 12)

202207101016111500..png

202207101017102778..png (27.38 KB, 下载次数 Times of downloads: 13)

202207101017102778..png
两个对易的厄米算子可以有共同本征函数集

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

5#
发表于 Post on 2022-7-26 06:50:07 | 只看该作者 Only view this author
nusiew 发表于 2022-7-10 10:24
谢谢sobereva老师和丁越老师的解答,对于cp2k默认输出的molden文件修改的地方总结为两个,一个是添加晶格常 ...

北京科音自然科学研究中心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

176

帖子

0

威望

1452

eV
积分
1628

Level 5 (御坂)

6#
发表于 Post on 2025-7-26 23:23:19 | 只看该作者 Only view this author
本帖最后由 杨小狗 于 2025-7-26 23:24 编辑

借贴寻求一个问题的答案:请问,cp2k在计算Cu(001)d-center的时候,将胞拓展8x8x1之后,产生的molden文件太大了,载入multiwfn之后加载完成自动闪退,使用[CELL].txt让molden文件载入cell信息之后,就不会闪退了,但重新排布电子占据的时候,排好后返回主界面依然会闪退,请问是molden文件太大造成的吗?
2.如果molden文件太大,ultraedit打开的时候也提示文件太大应该怎么解决?谢谢

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

7#
发表于 Post on 2025-7-27 01:59:39 | 只看该作者 Only view this author
杨小狗 发表于 2025-7-26 23:23
借贴寻求一个问题的答案:请问,cp2k在计算Cu(001)d-center的时候,将胞拓展8x8x1之后,产生的molden文件 ...

1 可能是。找个大内存机子
北京科音自然科学研究中心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

176

帖子

0

威望

1452

eV
积分
1628

Level 5 (御坂)

8#
发表于 Post on 2025-7-28 08:54:02 | 只看该作者 Only view this author
sobereva 发表于 2025-7-27 01:59
1 可能是。找个大内存机子

收到,感谢社长解答

176

帖子

0

威望

1452

eV
积分
1628

Level 5 (御坂)

9#
发表于 Post on 2025-7-28 17:47:16 | 只看该作者 Only view this author
杨小狗 发表于 2025-7-28 08:54
收到,感谢社长解答

加大机子内存,molden文件可以载入multiwfn了,原因出在内存不够,保存不下molden的数据

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 05:43 , Processed in 0.172909 second(s), 23 queries , Gzip On.

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