计算化学公社

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

[CP2K] 求助使用CP2K计算体系非键能随时间的变化曲线,能量怎么使用CP2K进行输出?

[复制链接 Copy URL]

53

帖子

0

威望

523

eV
积分
576

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师好,我想计算一个体系的非键能(Non-bonding energy)随时间的变化曲线,如下图所示:


这篇文章中非键能主要是Van der Waals and Coulomb interaction energies,请问这两个能量怎么使用CP2K进行输出。
我用GFN1-XTB的方法计算了一个MD任务,请问这个任务最后会输出这两个非键能吗?

6万

帖子

99

威望

5万

eV
积分
120158

管理员

公社社长

2#
发表于 Post on 2024-12-7 04:46:35 | 只看该作者 Only view this author
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并清楚、准确反映出帖子具体内容,避免有任何歧义和含糊性,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “使用CP2K计算体系非键能随时间的变化曲线” 改了,以后务必注意!
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120158

管理员

公社社长

3#
发表于 Post on 2024-12-7 04:48:11 | 只看该作者 Only view this author
那文章应当是用分子力场算的,GROMACS可以分别给出,CP2K不会单独给出不同物理成份
北京科音自然科学研究中心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 (御坂)

傻傻的木瓜

4#
发表于 Post on 2024-12-15 14:35:24 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2024-12-15 15:24 编辑
sobereva 发表于 2024-12-7 04:48
那文章应当是用分子力场算的,GROMACS可以分别给出,CP2K不会单独给出不同物理成份

社长您好,如果已用CP2K跑完分子动力学获得轨迹,原理上是否可以考虑采用您下面这两篇博文的方法,批量对各帧做能量分解并绘制成类似一楼的各成分对时间的曲线呢?
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685http://bbs.keinsci.com/thread-39446-1-1.html
使用Multiwfn做基于分子力场的能量分解分析
http://sobereva.com/442http://bbs.keinsci.com/thread-10907-1-1.html

编辑:使用场景举例,比如研究一个KCl+100H2O构成的水溶液体系中钾离子与水分子的相互作用(编辑2:应实验室要求确认此为合理的简化模型),有用CP2K跑AIMD的小几千帧轨迹也有用LAMMPS(实验室所限暂时没法用GROMACS)跑经典力场MD的大几万帧轨迹(转换成Multiwfn支持的文件格式),我在考虑的除了原理的问题还有批量计算的耗时。
√546=23.36664289109

6万

帖子

99

威望

5万

eV
积分
120158

管理员

公社社长

5#
发表于 Post on 2024-12-16 02:34:21 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2024-12-15 14:35
社长您好,如果已用CP2K跑完分子动力学获得轨迹,原理上是否可以考虑采用您下面这两篇博文的方法,批量对 ...

如果你是指自动抠出来团簇,结合Gaussian+Multiwfn自动做这样的能量分解分析(比如分析K+与第一配位层的水之间的作用),原理上可以
北京科音自然科学研究中心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 (御坂)

傻傻的木瓜

6#
发表于 Post on 2024-12-16 11:05:33 | 只看该作者 Only view this author
sobereva 发表于 2024-12-16 02:34
如果你是指自动抠出来团簇,结合Gaussian+Multiwfn自动做这样的能量分解分析(比如分析K+与第一配位层的 ...

需要抠团簇出来再计算,是因为上面那两种能量分解分析仅适用于孤立体系而不适用于周期性体系么?毕竟当钾离子浓度稍高一些时,无论盒子有多大都可能有配位层水分子在不同镜像间共享的情况,就是利用选择语句也不一定能取出明确的孤立团簇。

另一边,CP2K手册里也有个Reference指向J. Chem. Phys. 153, 214109 (2020)的
ENERGY_DECOMPOSITION_ANALYSIS字段,似乎可以兼顾周期性边界条件。把CP2K自带的benchmark中H2O-128.inp [注释为128 H2O (TIP5P,1bar,300K) a = 15.6404] 加载到Multiwfn生成单点能计算的输入文件,尝试用上述字段跑的结果如附件1所示。看起来CP2K先以STO-3G极小基为参照做IAO分析(这里跳过了默认开启、耗时超长的轨道定域化),然后分别用IAO和Mulliken布居分析产生的电荷来分解各原子能量,但是不完全对应前面提到的能量分解分析产生的成分。社长对此有何看法?

H2O-128.tar.xz (37.8 KB, 下载次数 Times of downloads: 2)
√546=23.36664289109

6万

帖子

99

威望

5万

eV
积分
120158

管理员

公社社长

7#
发表于 Post on 2024-12-17 03:00:51 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2024-12-16 11:05
需要抠团簇出来再计算,是因为上面那两种能量分解分析仅适用于孤立体系而不适用于周期性体系么?毕竟当钾 ...

抠的时候要考虑周期性,VMD里可以用pbwithin
抠团簇是为了用sobEDA等能量分解,只支持孤立体系

你说的CP2K里的那个能量分解没什么实际用处,那只是把总能量及其成份分解成原子贡献,并不给像常用的能量分解那样能给出片段间相互作用的物理成份

评分 Rate

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

查看全部评分 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

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 09:13 , Processed in 0.997856 second(s), 24 queries , Gzip On.

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