计算化学公社

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

[CP2K] cp2k跑分子动力学应该怎么分析结果?

[复制链接 Copy URL]

43

帖子

0

威望

549

eV
积分
592

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师,我使用cp2k跑了一个分子动力学,不知道怎么分析结果,目前只能分析距离随时间变化,我想得到一些关于电荷/电子转移还有成键的相关信息,应该如何做?
事以密成,语以泄败

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

2#
发表于 Post on 2024-9-24 21:42:06 | 只看该作者 Only view this author
结合Multiwfn做波函数分析,有的是方法
Multiwfn支持的分析化学键的方法一览
http://sobereva.com/471http://bbs.keinsci.com/thread-12488-1-1.html
谈谈怎么考察、计算、分析化学体系的电子密度
http://sobereva.com/715http://bbs.keinsci.com/thread-46380-1-1.html

可以每隔一定帧数导出一个molden文件,按下文批量用multiwfn进行分析
详谈Multiwfn的命令行方式运行和批量运行的方法
http://sobereva.com/612http://bbs.keinsci.com/thread-24929-1-1.html
上面的文章里就有现成的分析ORCA做AIMD过程中电子结构变化的例子


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

267

帖子

0

威望

941

eV
积分
1208

Level 4 (黑子)

3#
发表于 Post on 2024-11-13 13:57:07 | 只看该作者 Only view this author
sobereva 发表于 2024-9-24 21:42
结合Multiwfn做波函数分析,有的是方法
Multiwfn支持的分析化学键的方法一览
http://sobereva.com/471(h ...

sob老师请问一下,我已经用CP2K跑完了AIMD,现在应该已经不能再每隔一定帧数导出一个molden文件了。那我该怎么获得multiwfn需要的输入文件呢?是把xyz轨迹文件里的每一帧坐标都拆成单一文件计算单点能获得波函数吗?这样的话计算量似乎有点高

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

4#
发表于 Post on 2024-11-13 23:40:55 | 只看该作者 Only view this author
ljh123 发表于 2024-11-13 13:57
sob老师请问一下,我已经用CP2K跑完了AIMD,现在应该已经不能再每隔一定帧数导出一个molden文件了。那我 ...


没必要每步都产生,诸如可以每10帧产生一次
北京科音自然科学研究中心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

267

帖子

0

威望

941

eV
积分
1208

Level 4 (黑子)

5#
发表于 Post on 2024-11-14 13:56:31 | 只看该作者 Only view this author
sobereva 发表于 2024-11-13 23:40

没必要每步都产生,诸如可以每10帧产生一次

非常感谢sob老师非常细致认真的回答,总能解决我很多疑问。

894

帖子

3

威望

1685

eV
积分
2639

Level 5 (御坂)

傻傻的木瓜

6#
发表于 Post on 2024-11-14 15:43:19 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2024-12-15 15:24 编辑
sobereva 发表于 2024-11-13 23:40

没必要每步都产生,诸如可以每10帧产生一次

社长您好,我也在考虑利用之前您提到的REFTRAJ字段来这样重算部分MD轨迹,感觉似乎可以在重算中提高计算级别。比如一个KCl+100H2O构成的水溶液体系(编辑2:应实验室要求确认此为合理的简化模型),我打算先在便宜的PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH级别下用800 Ry截断能跑分子动力学,每5帧或10帧记录一个坐标的snapshot,然后在更贵的PBE0-D3(BJ)/TZVP-MOLOPT-PBE0-GTH级别下用稍高一些如1000 Ry截断能结合REFTRAJ重新计算每个snapshot的单点能和原子受力。坐标、单点能和原子受力是用于训练NNP的。

虽然高计算级别意味着时间加倍,但是保留用于高级别单点能计算的snapshot数量只有原先分子动力学的几分之一,如果估计同样计算资源(集群单个节点112核,单任务计算时长上限一周)单点能计算的总耗时和之前分子动力学模拟的总耗时相当的话还是可以接受的。请指点一下这个设计是否合理。

√546=23.36664289109

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

7#
发表于 Post on 2024-11-15 04:56:49 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2024-11-14 15:43
社长您好,我也在考虑利用之前您提到的REFTRAJ字段来这样重算部分MD轨迹,感觉似乎可以在重算中提高计算 ...

很合理
由于每一帧单独算单点+受力的时候没法像MD那样读取上一步收敛的波函数当初猜,所以应注意检查批量算单点+受力过程的SCF收敛情况

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
Uus/pMeC6H4-/キ + 1 感谢提醒

查看全部评分 View all ratings

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

894

帖子

3

威望

1685

eV
积分
2639

Level 5 (御坂)

傻傻的木瓜

8#
发表于 Post on 2024-11-15 11:41:42 | 只看该作者 Only view this author
sobereva 发表于 2024-11-15 04:56
很合理
由于每一帧单独算单点+受力的时候没法像MD那样读取上一步收敛的波函数当初猜,所以应注意检查批 ...

波函数收敛的确很重要呢。

其实从已有的水溶液研究文献J. Phys. Chem. Lett. 2017, 8, 7, 1545–1551Nature Chemistry 16, 1838–1844 (2024)来看,在CP2K中利用ADMM技术时计算级别甚至可达revPBE0-D3(BJ)方法结合TZV2P档次基组。不过还有几个小细节要确认一下是否正确:

1.我看到博文《CP2K做杂化泛函计算的关键要点和简单例子》写了PBE0泛函的定义方法,用了手册所述&PBE字段书写。我理解的revPBE0和revPBE的关系也类似于PBE0与PBE,revPBE0由75%的revPBE交换项+25%的HF交换项+100%的revPBE相关项(实与原PBE一致)构成,那么书写时只要加一行PARAMETRIZATION REVPBE就行。

2.在CP2K的data/BASIS_MOLOPT_UZH中有专为PBE0优化的TZV2P-MOLOPT-PBE0-GTH系列基组,但是用revPBE0时应该用回data/BASIS_MOLOPT_UCL里的TZV2P-MOLOPT(-SR)-GTH。(毕竟BASIS_MOLOPT_UZH没写原始文献)

3.上述JPCL论文采用的ADMM辅助基组是cpFIT3,但如今最好改成上述博文推荐的admm-tzp。
√546=23.36664289109

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

9#
发表于 Post on 2024-11-15 11:48:32 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2024-11-15 11:41
波函数收敛的确很重要呢。

其实从已有的水溶液研究文献J. Phys. Chem. Lett. 2017, 8, 7, 1545–1551 ...

1 或许正确,但我没具体了解过revPBE0。若无特殊必要用,还是建议用流行得多的PBE0

2 只要是杂化泛函,最好都用给PBE0优化的基组,原理上比用对纯泛函优化的基组更理想

3 对

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
Uus/pMeC6H4-/キ + 1 感谢

查看全部评分 View all ratings

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

894

帖子

3

威望

1685

eV
积分
2639

Level 5 (御坂)

傻傻的木瓜

10#
发表于 Post on 2024-11-15 13:23:07 | 只看该作者 Only view this author
sobereva 发表于 2024-11-15 11:48
1 或许正确,但我没具体了解过revPBE0。若无特殊必要用,还是建议用流行得多的PBE0

2 只要是杂化泛函 ...

用revPBE-D3做液态水体系AIMD的文章能看到好几篇(没准今年第一性原理泛函流行程度投票里面有因此投票给这个泛函的),感觉可能跟J. Phys. Chem. C 2014, 118, 50, 29401–29411比较若干泛函的结论有关:
It is found that, among these seven functionals, BLYP, BLYP-D2, and PBE-D3 functionals produce slower dynamics of spectral diffusion, while BLYP-D3, revPBE-D3, DRSLL-PBE, and DRSLL-optB88 functionals produce faster dynamics. The latter results are comparable to each other, and they are in good agreement with available experimental results.

想想也是,如果其他类型的体系也要低级别跑完MD后高级别重算单点能和受力,通用性更强的PBE0会更合适。
√546=23.36664289109

267

帖子

0

威望

941

eV
积分
1208

Level 4 (黑子)

11#
发表于 Post on 2024-12-13 20:52:27 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2024-11-14 15:43
社长您好,我也在考虑利用之前您提到的REFTRAJ字段来这样重算部分MD轨迹,感觉似乎可以在重算中提高计算 ...

你好,我在用REFTRAJ重算能量的时候遇到一些问题。因为我是用NPT跑的轨迹,因此每一帧的盒子尺寸都在波动,而使用REFTRAJ是似乎只能用一个盒子尺寸,轨迹只能用xyz格式,这样重跑的话对每一帧的能量似乎有一定影响。请问这种情况你遇到过吗,是不是只能单独提取结构算单点能了呢

894

帖子

3

威望

1685

eV
积分
2639

Level 5 (御坂)

傻傻的木瓜

12#
发表于 Post on 2024-12-13 21:02:48 | 只看该作者 Only view this author
ljh123 发表于 2024-12-13 20:52
你好,我在用REFTRAJ重算能量的时候遇到一些问题。因为我是用NPT跑的轨迹,因此每一帧的盒子尺寸都在波动 ...

我自己没做过这种情况的重跑,但手册说可以添加CELL_FILE_NAME字段指定记录盒子信息的.cell文件来重跑,应该用先前&MOTION/&PRINT/&CELL控制输出的文件就行。当然要保证先前&MOTION/&PRINT/&CELL/EACH和&MOTION/&PRINT/&TRAJECTORY/EACH的记录频率一致,盒子和轨迹的记录能对上。

在&FORCE_EVAL/&SUBSYS那里也需要定义一个盒子和结构,我感觉用重跑第一帧的就行。
√546=23.36664289109

267

帖子

0

威望

941

eV
积分
1208

Level 4 (黑子)

13#
发表于 Post on 2024-12-13 21:13:23 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2024-12-13 21:02
我自己没做过这种情况的重跑,但手册说可以添加CELL_FILE_NAME字段指定记录盒子信息的.cell文件来重跑, ...

我尝试过CELL_FILE_NAME字段,结果他的盒子信息默认用.cell文件最后一帧的盒子尺寸,还是无法实现NPT的重算

894

帖子

3

威望

1685

eV
积分
2639

Level 5 (御坂)

傻傻的木瓜

14#
发表于 Post on 2024-12-13 21:30:09 | 只看该作者 Only view this author
ljh123 发表于 2024-12-13 21:13
我尝试过CELL_FILE_NAME字段,结果他的盒子信息默认用.cell文件最后一帧的盒子尺寸,还是无法实现NPT的重 ...

读取盒子信息还需要把默认为F的VARIABLE_VOLUME设成T哦。
√546=23.36664289109

267

帖子

0

威望

941

eV
积分
1208

Level 4 (黑子)

15#
发表于 Post on 2024-12-13 21:34:09 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2024-12-13 21:30
读取盒子信息还需要把默认为F的VARIABLE_VOLUME设成T哦。

原来如此,十分感谢!还有就是想请教下,通过REFTRAJ如何输出动能呢?我用REFTRAJ得到的.ener文件里只有势能,其他全为0

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 15:40 , Processed in 0.529495 second(s), 21 queries , Gzip On.

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