计算化学公社

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

[CP2K] CP2K SCCS单点时Solvation energy是否能直接读out文件中的Solvation free energy

[复制链接 Copy URL]

4

帖子

0

威望

29

eV
积分
33

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 heiheida 于 2026-1-25 11:53 编辑

各位老师们好,我目前使用CP2K的SCCS进行单点,水溶剂条件,RHO_MIN 和 RHO_MAX分别设置为0.0001以及RHO_MAX 0.005,截断能为600 Ry。
根据输出文件的  SCCS| Solvation free energy [Hartree]  -0.07312178911993部分,Solvation energy应该是-1.88 eV,但是液相中总能量-1445.436192490559506 a.u. ,气相中-1445.499683535024133 a.u. ,根据溶剂中总能量和气相中总能量能量之差,Solvation energy却是+1.72 eV 。

6万

帖子

99

威望

6万

eV
积分
125115

管理员

公社社长

2#
发表于 Post on 2026-1-24 21:53:54 | 只看该作者 Only view this author
说明液相中总能量和气相中的能量都是哪读的
北京科音自然科学研究中心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

4

帖子

0

威望

29

eV
积分
33

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2026-1-24 22:26:58 | 只看该作者 Only view this author
sobereva 发表于 2026-1-24 21:53
说明液相中总能量和气相中的能量都是哪读的

好的老师,我是在真空下进行结构优化,然后把优化好的结构在溶剂以及真空下进行单点,液相和气相中的能量度的就是溶剂以及真空下的单点能

6万

帖子

99

威望

6万

eV
积分
125115

管理员

公社社长

4#
发表于 Post on 2026-1-25 15:52:45 | 只看该作者 Only view this author
heiheida 发表于 2026-1-24 22:26
好的老师,我是在真空下进行结构优化,然后把优化好的结构在溶剂以及真空下进行单点,液相和气相中的能量 ...

以你这么自己算的为准
SCCS模型下输出文件里直接给你的Solvation energy只是个近似,因为算其中极性部分贡献的时候,并没有用到真正在真空下迭代求解出的波函数对应的能量。
北京科音自然科学研究中心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

4

帖子

0

威望

29

eV
积分
33

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2026-1-25 23:36:39 | 只看该作者 Only view this author
sobereva 发表于 2026-1-25 15:52
以你这么自己算的为准
SCCS模型下输出文件里直接给你的Solvation energy只是个近似,因为算其中极性部分 ...

好的,谢谢老师!请问老师我这个溶剂下的能量是不是有很大问题,需要继续调整RHO_MIN 和 RHO_MAX使溶剂下能量更低?
我测过截断能500 Ry到750 Ry,能量变化不大。然后也试过截断能为600 Ry下把RHO_MIN 和 RHO_MAX调整为0.0001 以及 0.001,Solvation energy从+1.72 eV变为+0.25 eV左右。
我用VASP在真空以及溶剂条件下对结构进行单点,通过能量差算出来的Solvation energy是-0.66 eV左右。

6万

帖子

99

威望

6万

eV
积分
125115

管理员

公社社长

6#
发表于 Post on 2026-1-25 23:50:33 | 只看该作者 Only view this author
heiheida 发表于 2026-1-25 23:36
好的,谢谢老师!请问老师我这个溶剂下的能量是不是有很大问题,需要继续调整RHO_MIN 和 RHO_MAX使溶剂下 ...

不知道你算的是什么东西,无法评论。有的体系在有的溶剂里的溶解自由能就是正的。如果是孤立体系,可以和Gaussian结合SMD算的溶解自由能对比判断合理性。如果确信溶解自由能理应是负的,可尝试进一步调整RHO_MIN 和 RHO_MAX(但如果是之后要能量横向求差的话,由于误差抵消效应,溶解自由能为正不一定是必须解决的问题)

另外,也要看溶解自由能的极性和非极性部分各自是多少。如果连前者都为正,那说明问题就很严重了。
北京科音自然科学研究中心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

4

帖子

0

威望

29

eV
积分
33

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2026-1-31 20:17:07 | 只看该作者 Only view this author
本帖最后由 heiheida 于 2026-1-31 20:20 编辑
sobereva 发表于 2026-1-25 23:50
不知道你算的是什么东西,无法评论。有的体系在有的溶剂里的溶解自由能就是正的。如果是孤立体系,可以和 ...

谢谢老师!我在算溶剂化下N2的吸附能,我看了我的输出Polarisation energy是正的,所以我的这个能量有很大问题。我一直在尝试改变sccs相关参数,但是不管如何Polarisation energy都是正的。最后将TZV2P-MOLOPT-PBE-GTH换成了DZVP-MOLOPT-SR-GTH之后Polarisation energy以及计算的Solvation energy就都对了。
我想请教一下,我之前看算能量最好用TZV2P级别,那么我是否可以用DZVP-MOLOPT-SR-GTH下算出来的Solvation energy加上TZV2P下的单点能(无sccs)作为体系的能量呢?

6万

帖子

99

威望

6万

eV
积分
125115

管理员

公社社长

8#
发表于 Post on 2026-2-1 00:42:09 | 只看该作者 Only view this author
heiheida 发表于 2026-1-31 20:17
谢谢老师!我在算溶剂化下N2的吸附能,我看了我的输出Polarisation energy是正的,所以我的这个能量有很 ...

原理上没什么问题
北京科音自然科学研究中心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, 2026-2-16 16:56 , Processed in 0.186368 second(s), 20 queries , Gzip On.

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