关于CP2K中自相互作用校正（SIC）的参数设置和结果的问题

gingaorb

各位小哥哥，各位大佬们，请教一个问题呀~
是这样的，我研究的是一个丙酮酸分子在一个含有水合电子的水团簇中的体系，为了比较准确的描述水合电子和丙酮酸之间的相互作用，我参考了文献“A molecular dynamics study of the hydroxyl radical in solution applying self-interaction-corrected density functional methods”中的研究使用一个自旋密度自相互作用校正（SD-SIC）方法。我用的泛函是BLYP-D3(BJ)，基组是TZV2P-GTH。
根据参考文献“Degradation of Per- and Polyfluoroalkyl Substances with Hydrated Electrons: A New Mechanism from First-Principles Calculations”中的方法，确定SD-SIC方法中SIC_SCALING_A和SIC_SCALING_B参数的数值，要通过所研究体系关键物种的垂直电子亲和能（VEA）来确定。
因此我的研究思路是：先通过Gaussian优化丙酮酸分子结构（b3lyp/def2tzvp+D3(BJ)），在分别计算丙酮酸的单点能和带-1电荷丙酮酸的单点能（CCSD(T)/AUG-cc-pVTZ），得到丙酮酸单点能-341.9286534 a.u.以及带-1电荷丙酮酸的单点能-341.9261644 a.u.（三者文件名分别为B3LYP-opt.gif、CCSD-sp.gif、CCSD-1-sp.gif）。
然后我在CP2K中用Gaussian优化好的结构分别计算丙酮酸的单点能（-67.3296633424807 a.u.）和带-1电荷丙酮酸的单点能（两者者文件名分别为cp2k-Pyruvate_SP.inp、cp2k-Pyruvate-1_SD-SIC_SP.inp）,其中设置的SIC部分的关键词和参数为：
    &TDDFPT
      &SIC
        SIC_METHOD AD
        SIC_SCALING_A 0.2
        SIC_SCALING_B 0
      &END SIC
    &END TDDFPT

然后，奇怪的情况发生了，无论我如何改变SIC_METHOD、SIC_SCALING_A和SIC_SCALING_B，所得到的结果都完全相同！！！为-67.3497019097407 a.u.，就很离谱，我不明白会出现这种情况的原因是什么？为什么我无论如何改变SIC参数设置计算得到的单点能结果都相同？是我输入文件哪里写错了还是我对该方法的理解有错误？请各位小哥哥，各位大佬们指点指点迷津好吗？o(╥﹏╥)o

sobereva

如果模拟总时间不很长、体系不很大，与其折腾这个，不如直接用杂化泛函跑

gingaorb

sobereva 发表于 2023-6-20 03:08
如果模拟总时间不很长、体系不很大，与其折腾这个，不如直接用杂化泛函跑

体系191个水分子模拟时间估计在20 ps左右，担心用B3LYP或者PBE0可能跑不动~o(╥﹏╥)o，请问站长大大有没有什么好的建议吗？

sobereva

gingaorb 发表于 2023-6-21 10:11
体系191个水分子模拟时间估计在20 ps左右，担心用B3LYP或者PBE0可能跑不动~o(╥﹏╥)o，请问站长大大有没 ...

先跑起来试试，根据跑10步以后的平均每步耗时估计总耗时是否能接受。
要么联系文章作者要输入文件，并且表示会引用他们的文章，很可能直接给你