求助Fe3+配合物的优化中出现挤在一起的不合理构型

suosuosky

求助各位老师，我在优化还有Fe3+d的配合物时，经常出现拥挤在一起的情况（很不合理），不知道是什么原因造成的，请各位老师指导一下。

（1）Fe3+结合6个H2O进行优化，使用关键词：
# opt freq b3lyp/genecp scrf=(PCM,solvent=water) SCF=xqc
其中加SCF=xqc关键词，是一直不收敛尝试了很多次后才加的，结果优化很多步后构型如下图所示（没算完，但是貌似构型很不合理了）。
输出文件已上传（Fe3+_6w.log）

（2）Fe3+结合另外一个体系优化，使用关键词：
#p opt(gdiis,maxstep=5,notrust,recalc=3) freq b3lyp/genecp scrf=(solve nt=h2o,pcm)
结果还是出现了挤在一起的情况，输出文件已上传（M+Fe3+_2w.log）。
是不是构型搭建或者关键词使用有问题，因为在搭建的其他Fe3+配合物时也有类似的问题新出现。

请各位老师指导，谢谢大家！

zjxitcc

输出文件里很多警告 Warning!  O  atom    2 has  8 valence electrons but only   0 basis functions.
如果你有经验的话，其实看到“拥挤在一起”就会反应出是输入文件里基组格式有误，有些原子没有指定基组。

tianmafei

原因应该是没有给O原子指定基组。把输入文件贴出来看一下。

suosuosky

zjxitcc 发表于 2022-3-14 20:18
输出文件里很多警告 Warning!  O  atom    2 has  8 valence electrons but only   0 basis functions.
如 ...

谢谢老师的回复和提醒

ionexchangeC

这种体系的SCF收敛难度通常不是很高，所以遇到这种只能用二次收敛方法才能收敛的情况，可以先考虑一下是不是关键词写错了（除了未给原子指定基组之外，还有一个常见的就是自定义赝势基组时写了gen而不是genecp）