我明白了，谢谢sob老师和zjxitcc老师百忙之中的指点！！！

如果用当前的簇模型是为了近似表现固体，这么算纯属瞎折腾。SCF总是很难收敛，而且还不能不考虑边界效应（需要加背景电荷、capped ECP，很折腾），而且还没法优化晶胞参数。

CP2K优化这种晶体又快又容易又不花钱

好的，因为我们只有Gaussian的版权，所以才没考虑其他，那我用CP2K试试，非常感谢老师的帮助

不用找报错、不用解决报错、不要试图解决SCF不收敛、不要再用孤立体系量化程序做这个体系计算，应立即改用周期性计算程序（例如QE, VASP, CP2K等）。