g16计算簇模型含有两个C自由基的体系自旋极化单重态，自旋密度处处为0问题求助

sun877469558

各位老师好，本研究是簇模型含有两个C自由基的体系，自由基进行耦合，

计算自旋极化单重态时遇到了问题，
使用  ①  #p ub3lyp/6-31G* em=gd3bj guess=mix nosymm stable=opt   （结果：自旋密度处处为0）
        ②  #p ub3lyp/6-31G* em=gd3bj guess=mix nosymm stable （结果：波函数不稳定）
        ③  #p ub3lyp def2svp em=gd3bj guess=mix nosymm stable=opt（结果：自旋密度处处为0）
        ④  #p ub3lyp def2svp  em=gd3bj guess=mix nosymm stable（结果：自旋密度处处为0）

使用四种关键词去跑都无法到最稳定基态对称破缺波函数

附件是两个碳自由基不同距离的输入与输出文件，分别是1.9埃与3.8埃

请问各位老师如何解决这个问题？

Uus/pMeC6H4-/キ

这团簇是什么原始结构按什么方法截取的，确认周围一堆甲烷乙烷之类的小分子都有必要出现么，是否算完单点能还需要几何优化？除了C193-C195-N190所在分子有点像自由基之外，没看出来第二个自由基在哪、距离是量取的哪里，有没有尝试过取出最核心的两个碳自由基单独这样算？

guess=mix只是一种通过混合HOMO与LUMO来准备波函数初猜的方式，我觉得对如此复杂的团簇体系不如指定原子的片段初猜合适。

sobereva

没法简单看出单电子位点在哪。做个三重态计算或许还能了解一下单电子出现位点的倾向性
明显是闭壳层的分子别加进去，纯属添乱、浪费计算时间。至少等你确定了单电子位点后再考虑
另外，注意HF成份高的泛函出现对称破缺的倾向性比B3LYP更强，必要时考虑用wB97XD