求助Ti内嵌C60（Ti@C60）开壳层单线态计算遇到的问题

JTDJZS

各位老师们好，我在计算Ti内嵌C60（Ti@C60）的结构优化，因为发现该结构的自旋多重度为1和3是的能量只相差了0.9kcal/mol，因此尝试计算开壳层下自旋多重度为1时的结构，阅读了谈谈片段组合波函数与自旋极化单重态 - 思想家公社的门口：量子化学·分子模拟·二次元帖子后，先是使用关键词：
# utpsstpss/genecp nosymm empiricaldispersion=gd3bj guess=mix stable
获得了稳定的波函数，并查看了自旋并不是处处为零，并且<S**2>为0.475，因此以此波函数进行结构优化：
# opt freq utpsstpss/genecp nosymm empiricaldispersion=gd3bj guess=read
得到的结构和闭壳层自旋多重度为1下优化得到的结构能量是完全一样的，<S**2>变为了0这是什么原因呢？
我还采用了# opt freq utpsstpss/genecp nosymm empiricaldispersion=gd3bj guess=(mix,always)关键词，也得到了相同的结构和能量，<S**2>也是0，这是说明并不存在开壳层单线态的情况吗？请求各位老师解答，并附上我的闭壳层单线态和开壳层单线态的输入输出文件。

sobereva

你这种体系根本不会有自旋极化单重态的电子态

LightRain

有人在C60发现过symmetry breaking solution
https://doi.org/10.1063/1.3661158
所以Ti内嵌也有可能发现大于0的<S^2>。虽然TPSS这种纯泛函比HF的symmetry breaking的程度可能少些

另外非限制性+guess = mix和stable = opt不是同一种方法，后者可能找到能量更低的解（虽然楼主的输入没用stable=opt)；以及非限制性+guess=mix+stable=opt和限制性+stable=opt也不相同

可能从非限制性+guess = mix得到的symmetry breaking的解，guess=read后回到了没有breaking的解