请问cp2k是否支持在k点下进行sccs计算？

kogareru1z

很多人都提到过现在的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

sobereva

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

UW_0728.

起码在开发者方面未经充分验证。你贴出来的那个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上给开发者看一看。