SiO2表面吸附过渡金属配合物计算中SCF不收敛的问题

bogong

各位老师好：本人目前正在尝试计算SiO2表面对抗生素美他环素（MTC）以及对应的Cu2+配合物的吸附，前者计算一切正常，但后者含有金属配合物的体系在启动时即遇到了SCF不收敛的问题（并且完全没有收敛趋势）。经过个人的摸索以及尝试了社长博文《CP2K中遇到SCF难收敛时的解决方法》中仍然没有解决，目前确认的一些结果如下：
（1）本人去掉表面，单独尝试优化上方的(Cu2+-MTC)络合物本身，也是同样的问题——说明SCF不收敛并不是因为复合物中的吸附构型初猜不合理，而是络合物本身的问题
（2）关于络合物本身，目前该络合物中Cu2+虽然不是完全饱和的配位状态，但该结构在Gaussian中使用B3LYP/def2SVP这样的底线计算级别是完全可以正常优化的，因此结构上也具有一定合理性，且这种不饱和配位结构也的确是我们想研究的内容——说明SCF不收敛可能不是结构严重不合理，而是计算参数设置的问题
（3）目前在CP2K中的计算级别已经是DFT中的底线级别，即PBE泛函搭配DZVP-MOLOPT-SR-GTH-q6基组，根据对照组SiO2对MTC本身的吸附测试来看，若将计算级别降低至GFN1-xTB，则优化出的吸附构型与DFT级别下有明显区别——说明对于该体系GFN1-xTB可能精度过低，最好还是使用DFT级别的方法作为底线框架
（4）参考社长博文，我们尝试了使用smearing或者打开OT等方式，暂时都没有取得明显效果。由于本人对原理的理解尚且有限，为了避免盲目地继续尝试，在此附上我的结构文件和输入输出文件，谨向各位老师求教，感谢各位老师的帮助！