求助 如何用CP2K计算显示溶剂化表面能

chemistryMaga

各位老师好，我最近通过CP2K计算显示溶剂化表面能时遇到一些困惑，希望能得到老师们一些指点。首先我通过阅读文献发现显示水溶剂下矿物表面能的计算公式如下：
其次，仿照文献将碳酸钙体系中加入水并做了动力学模拟。
并以此作为输入文件通过Multiwfn来生成Cp2k的inp文件。选择了变胞优化并使Z轴方向尺寸不变，K点设置为 5 5 5，Cp2k顺利跑完程序。但带入公式计算时发现结果是一个很大的正值，不知道问题出在哪里，希望能得到老师们的一些帮助！