计算化学公社

标题: CP2K SCCS单点时Solvation energy是否能直接读out文件中的Solvation free energy [打印本页]

作者
Author: heiheida 时间: 2026-1-24 17:09
标题: CP2K SCCS单点时Solvation energy是否能直接读out文件中的Solvation free energy
本帖最后由 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 。

作者
Author: sobereva 时间: 2026-1-24 21:53
说明液相中总能量和气相中的能量都是哪读的

作者
Author: heiheida 时间: 2026-1-24 22:26

sobereva 发表于 2026-1-24 21:53
说明液相中总能量和气相中的能量都是哪读的

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

作者
Author: sobereva 时间: 2026-1-25 15:52

heiheida 发表于 2026-1-24 22:26
好的老师，我是在真空下进行结构优化，然后把优化好的结构在溶剂以及真空下进行单点，液相和气相中的能量 ...

以你这么自己算的为准
SCCS模型下输出文件里直接给你的Solvation energy只是个近似，因为算其中极性部分贡献的时候，并没有用到真正在真空下迭代求解出的波函数对应的能量。

作者
Author: heiheida 时间: 2026-1-25 23:36

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左右。

作者
Author: sobereva 时间: 2026-1-25 23:50

heiheida 发表于 2026-1-25 23:36
好的，谢谢老师！请问老师我这个溶剂下的能量是不是有很大问题，需要继续调整RHO_MIN 和 RHO_MAX使溶剂下 ...

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

另外，也要看溶解自由能的极性和非极性部分各自是多少。如果连前者都为正，那说明问题就很严重了。

作者
Author: heiheida 时间: 2026-1-31 20:17
本帖最后由 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）作为体系的能量呢？

作者
Author: sobereva 时间: 2026-2-1 00:42

heiheida 发表于 2026-1-31 20:17
谢谢老师！我在算溶剂化下N2的吸附能，我看了我的输出Polarisation energy是正的，所以我的这个能量有很 ...

原理上没什么问题

欢迎光临计算化学公社 (http://bbs.keinsci.com/)