标题: 关于CP2K中自相互作用校正（SIC）的参数设置和结果的问题 [打印本页] 作者Author: gingaorb 时间: 2023-6-20 02:17 标题: 关于CP2K中自相互作用校正（SIC）的参数设置和结果的问题 各位小哥哥，各位大佬们，请教一个问题呀~
是这样的，我研究的是一个丙酮酸分子在一个含有水合电子的水团簇中的体系，为了比较准确的描述水合电子和丙酮酸之间的相互作用，我参考了文献“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