beauyy 发表于 2024-6-1 01:26
是的，fchk文件里的数值和multiwfn里显示的一样，我刚才重算过还是这样的结果

不知道为什么这张图我这边加载不出来。如果是一堆无序的有正有负的数，应该是bug。

hebrewsnabla 发表于 2024-5-31 21:21
如果单看这个占据数，可能是bug。

fch里面的alpha orbital energy部分也是这样么？

是的，fchk文件里的数值和multiwfn里显示的一样，我刚才重算过还是这样的结果

如果单看这个占据数，可能是bug。

fch里面的alpha orbital energy部分也是这样么？

beauyy 发表于 2024-5-31 19:32
后续的UHF计算结果显示S2值是8.420901，比波函数稳定这一步的结果更大，也显示有多个自由基。使用的关键 ...

我没用高斯做过UNO的计算。不过我猜测是因为你写了hf关键词，程序做的是RHF，所以存的自然轨道也是RHF的。单重态的时候只有写UHF才会做UHF，不管读的轨道是RHF还是UHF。

如果还有问题可以上传这一步的输出文件。

如果说你下一步做UHF的时候没有发现<S**2>大于0，可以考虑检查波函数稳定性或（且）读取PBE0的轨道。并不是说guess=mix就一定能得到自旋极化的结果，这只是手段之一。

hebrewsnabla 发表于 2024-5-31 16:41
这个是哪一步计算的呢

是DFT做波函数优化这一步做的，关键词是stable=opt guess=(read,mix) def2tzvp em=gd3bj upbe1pbe

beauyy 发表于 2024-5-31 16:30
计算结果显示S2=2.38942，这个结果是不是说明有不止2个自由基？跟HF结果还是不一致，请问产生不一致的原 ...

这个是哪一步计算的呢

你这里提到的多个openshell singlet 的计算中，是否有任何一个的<S**2>大于0？

做openshell singlet 的计算，判断计算是否基本正确，首先要看的是<S**2>，因为这个不需要什么额外操作。<S**2>不对的话什么自由基数量就不谈了。