计算化学公社

标题: 请问cp2k是否支持在k点下进行sccs计算？ [打印本页]

作者
Author: kogareru1z 时间: 2026-3-17 16:58
标题: 请问cp2k是否支持在k点下进行sccs计算？
很多人都提到过现在的cp2k不支持在k点下进行隐式溶剂sccs计算，但我在cp2k的github的issue里看到了相关讨论，有人测试发现在k点下cp2k也能运行而不会报错https://github.com/cp2k/cp2k/issues/4854
，于是我也用简单的体系进行了测试，大家能否分析下结果如何？

kpoint  sccs Total energy          G_pol             G_cav             G_disp             G_rep
gamma none -364.88870917641327
221    none -364.85357375844427
441    none -364.85173616660711
881    none -364.85280089293377
12121 none -364.85283231450165
gamma on    -364.89078089295714 -0.00017897341641 0.00000000000000 -0.07607055191891 0.02229308507076
221    on    -364.76088528382962    0.00142187395651 0.00000000000000 -0.07432420500056 0.02231312531325
441    on    -364.76023696627271    0.00142420698242 0.00000000000000 -0.07432435369119 0.02231342131549
881    on    -364.76033126723001    0.00142647071695 0.00000000000000 -0.07432430992082 0.02231341765116
12121 on    -364.76034800553265    0.00142503612488 0.00000000000000 -0.07432429303536 0.02231342278912

作者
Author: sobereva 时间: 2026-3-18 05:29
官方明确说过不支持，至少开发者没测试、验证过合理性。即便考虑k点时结果在数值上看起来没问题、能随k点增加收敛，也不代表结果就正确，使用此数据需要后果自负

作者
Author: UW_0728. 时间: 2026-3-18 06:13
本帖最后由 UW_0728. 于 2026-3-18 06:38 编辑

起码在开发者方面未经充分验证。你贴出来的那个issue链接里面就讨论得很明白了：添加SCCS功能的那位老师说他写这段代码时根本没有考虑k点下的可行性，所以对SCCS在k点下能正常工作表示很惊讶；不过也表明了也许SCCS真的就天然兼容k点（或者只需要非常简单的改进）。按我个人对源代码的理解（可能有误），SCCS直接考虑k点用问题不大，因为代码里面是接收现成的通过每步SCF计算好的电子密度（无需区分是Gamma-only还是多k点）然后再进行SCCS迭代。

如果时间精力足够充裕，想要充分探索和验证SCCS与k点的兼容性，你需要做更多、做很多不同角度的测试，不限于收敛性这一方面。最简单的，多k点（a*a*a）下溶剂化能数据的a*a*a倍与相应a*a*a超胞的数据作比较；此外还要进行与受力（静态的受力计算，以及不变胞的几何优化，等等）等属性相关的测试，甚至不同类的体系都要测一遍。完后你也可以把你充分测试后的结果分享到那则issue上给开发者看一看。

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