计算化学公社

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

[CP2K] 请问cp2k怎样计算吉布斯自由能

[复制链接 Copy URL]

117

帖子

0

威望

1223

eV
积分
1340

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 reid 于 2022-11-24 19:42 编辑

使用cp2k 2022.1版,拿氧气分子练习计算自由能。先优化后计算频率,输入文件见O2-opt.inp和O2-freq.inp。均顺利计算完毕,但查看频率结果时(O2-freq.out)发现“Gibbs energy correction [kJ/mol]: Infinity”,请问这是什么原因?请教如何用cp2k计算自由能,谢谢!

O2-opt.inp

3.69 KB, 下载次数 Times of downloads: 76

O2-freq.inp

3.48 KB, 下载次数 Times of downloads: 83

O2-freq.out

55.61 KB, 下载次数 Times of downloads: 55

255

帖子

1

威望

2666

eV
积分
2941

Level 5 (御坂)

2#
发表于 Post on 2022-11-24 20:11:08 | 只看该作者 Only view this author
最简单的方法,使用社长的Shermo处理你的振动分析输出文件,直接就按照Gaussian输出的样式给你了
http://sobereva.com/soft/shermo/
我唯一知道的就是我一无所知,但我是化学小迷弟

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

3#
发表于 Post on 2022-11-24 20:20:34 | 只看该作者 Only view this author
你加了三个方向的PBC,你把PBC改成NONE试试。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

117

帖子

0

威望

1223

eV
积分
1340

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2022-11-24 21:06:40 | 只看该作者 Only view this author
chands 发表于 2022-11-24 20:20
你加了三个方向的PBC,你把PBC改成NONE试试。

您的意思是都改成非周期?

117

帖子

0

威望

1223

eV
积分
1340

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-11-24 21:10:16 | 只看该作者 Only view this author
ChemG 发表于 2022-11-24 20:11
最简单的方法,使用社长的Shermo处理你的振动分析输出文件,直接就按照Gaussian输出的样式给你了
http://s ...

您好,我知道社长这个软件。但是输出文件里如果没有自由能校正因子也能用软件导出来?谢谢

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

6#
发表于 Post on 2022-11-24 21:19:48 | 只看该作者 Only view this author
reid 发表于 2022-11-24 21:06
您的意思是都改成非周期?

对。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

117

帖子

0

威望

1223

eV
积分
1340

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-11-24 21:29:41 | 只看该作者 Only view this author
reid 发表于 2022-11-24 21:06
您的意思是都改成非周期?

PBC那里改成NONE,但Gibbs Engergy correction 仍旧是 InfinityF:\serve\常用分子热力学能\1.bmp

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

8#
发表于 Post on 2022-11-24 21:35:15 | 只看该作者 Only view this author
reid 发表于 2022-11-24 21:29
PBC那里改成NONE,但Gibbs Engergy correction 仍旧是 Infinity

文件没传上来。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

9#
发表于 Post on 2022-11-25 00:50:35 | 只看该作者 Only view this author
reid 发表于 2022-11-24 21:10
您好,我知道社长这个软件。但是输出文件里如果没有自由能校正因子也能用软件导出来?谢谢

哪来的什么所谓的“自由能校正因子”,根本没这种东西。那叫自由能热校正量。
CP2K只要能正常输出振动频率,Shermo就能算出热力学量
北京科音自然科学研究中心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

117

帖子

0

威望

1223

eV
积分
1340

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2022-11-25 08:50:06 | 只看该作者 Only view this author
非常感谢几位的回答,问题已解决,就是用sob大神的Shermo,流程很简单,先结构优化,再跑freq,最后跑单点。把单点能填到ini中,再把freq的out文件拖入Shermo窗口,马上就出来结果。再一次佩服大神!

14

帖子

0

威望

210

eV
积分
224

Level 3 能力者

11#
发表于 Post on 2024-9-26 22:17:04 | 只看该作者 Only view this author
reid 发表于 2022-11-25 08:50
非常感谢几位的回答,问题已解决,就是用sob大神的Shermo,流程很简单,先结构优化,再跑freq,最后跑单点 ...

请问吧单点能填到ini中是具体指什么,直接用freq的out拖入Shermo窗口里面的总校正能量可以代表该状态下的总能量吗

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

12#
发表于 Post on 2024-9-27 07:54:22 | 只看该作者 Only view this author
HAC 发表于 2024-9-26 22:17
请问吧单点能填到ini中是具体指什么,直接用freq的out拖入Shermo窗口里面的总校正能量可以代表该状态下的 ...

先把此文仔细看了
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552http://bbs.keinsci.com/thread-17494-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

14

帖子

0

威望

210

eV
积分
224

Level 3 能力者

13#
发表于 Post on 2024-9-27 09:30:18 | 只看该作者 Only view this author
sobereva 发表于 2024-9-27 07:54
先把此文仔细看了
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/55 ...

嗯嗯,谢谢sob老师

25

帖子

0

威望

239

eV
积分
264

Level 3 能力者

14#
发表于 Post on 2024-10-29 09:14:18 | 只看该作者 Only view this author
reid 发表于 2022-11-25 08:50
非常感谢几位的回答,问题已解决,就是用sob大神的Shermo,流程很简单,先结构优化,再跑freq,最后跑单点 ...

说到底就是最后也没用CP2K,直接用shermo算自由能就行对吧?

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

15#
发表于 Post on 2025-2-15 22:02:11 | 只看该作者 Only view this author
sobereva 发表于 2024-9-27 07:54
先把此文仔细看了
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/55 ...

卢老师,shermo一般是针对静态的分子或者周期性构型进行自由能计算的。如果是用的cp2k中的fpmd功能,这种任务能否像gmx一样使用US伞状采样做自由能计算呢
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 17:05 , Processed in 0.194696 second(s), 23 queries , Gzip On.

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