有机分子Annihilation of the first spin contaminant报错求助

bean233

论坛里的各位老师大佬们好，我在使用Gaussian16在对一个有机分子进行结构优化时，使用的B3LYP/6-31G+(d,p)的计算方法结构优化以及频率计算，总是Annihilation of the first spin contaminant的报错，使用了论坛里解决SCF不收敛的方法贴中的绝大多数方法，发现还是会出现报错，请问是哪里出了问题？报错文件以及gjf文件均放在了压缩包中，麻烦各位老师帮忙答疑解惑，感谢！

zjxitcc

（1）Annihilation of the first spin contaminant不是报错，是正常输出，与你的SCF不收敛问题无关。当前结构似乎不太合理，硝基-NO2的朝向看起来非常扭曲，最好是先用GFN2-xTB或PM7方法先优化一下，提取优化好的结构 再用DFT优化。
（2）在（1）的基础上，如果还有SCF不收敛、难收敛问题，可以加关键词scf(xqc,maxcycle=300)，不用添加其他你认为能帮助SCF收敛的关键词。优化完成后，需读取波函数，检查UB3LYP波函数稳定性，如果波函数稳定，那么频率分析有无虚频 这个结论才有意义。注意，请勿跳过（1）直接进行（2），那样会造成对我回答的曲解。
（3）你上传的文件里写的基组是6-31g(d,p)，不是你文字里说的6-31G+(d,p)，二者应当对应，不要随便凭记忆敲。
（4）你当前算的是有机分子结合两个Li原子。如果你脑袋里想算的是有机分子结合两个Li+离子，则说明你的想法与你的计算文件不一致，注意检查。

bean233

zjxitcc 发表于 2025-9-15 17:10
（1）Annihilation of the first spin contaminant不是报错，是正常输出，与你的SCF不收敛问题无关。当前结 ...

好的好的，感谢老师解答，我现在就去尝试一下。
您提到的（4），我目前计算的这个有机分子嵌入锂离子要去质子化，整体电荷数还是为0，您说的错误是这个吗？还是说还存在其他问题

zjxitcc

bean233 发表于 2025-9-15 17:13
好的好的，感谢老师解答，我现在就去尝试一下。
您提到的（4），我目前计算的这个有机分子嵌入锂离子要 ...

那你记得删掉两个H，我不太清楚你的结构，所以看不出你是否删了两个H