使用stable=opt之后，log文件最后一段出现ub3lyp。

倪宇晴

输入文件是
%mem=6GB
%nprocshared=14
%chk=Me+O2_stable.chk
# stable=opt b3lyp/6-311+g(d,p) empiricaldispersion=gd3bj

Title Card Required

0 1
O                 -2.67394900   -0.32488400    0.56784500
O                 -2.79527300   -0.66071200   -0.58348000
O                  2.72235600   -0.18395800   -0.53913800
O                  2.74747200   -0.80577600    0.49324300
C                 -0.00049900    1.57987700    0.04918400
H                 -0.86460700    2.23766000   -0.05307400
H                  0.10916700    0.97966300   -0.85554000
H                  0.89584100    2.18357100    0.19706900
H                 -0.14225300    0.92248600    0.90867800


log文件最后一段为

为什么log文件最后一段会出现UB3LYP。请各位老师指点迷津。

sobereva

虽然一开始是默认用RB3LYP，但由于stable=opt会考察放宽限制时能量是否降低，所以考虑了非限制性开壳层计算的情况

倪宇晴

sobereva 发表于 2018-1-18 22:46
虽然一开始是默认用RB3LYP，但由于stable=opt会考察放宽限制时能量是否降低，所以考虑了非限制性开壳层计算 ...

谢谢老师答复

倪宇晴

sobereva 发表于 2018-1-18 22:46
虽然一开始是默认用RB3LYP，但由于stable=opt会考察放宽限制时能量是否降低，所以考虑了非限制性开壳层计算 ...

老师，我优化的是一个甲烷和两个氧气分子。
stable=opt的log文件里有The wavefunction has an RHF -> UHF instability.的语句出现
读取stable=opt的chk继续优化的时候发现变成了UB3LYP。
我的目的是优化的是单重态的氧和甲烷。
这是不是意味着优化单重态的氧需要用开壳层。
还有我如何确定我优化的结果里的氧是单重态的。

sobereva

倪宇晴 发表于 2018-1-19 18:44
老师，我优化的是一个甲烷和两个氧气分子。
stable=opt的log文件里有The wavefunction has an RHF -> UH ...

当前已经暗示，基态就是开壳层了。你应当绘制自旋密度，看看到底单电子是怎么分布的，搞清楚基态时候体系真正的电子结构状态
如果你要把体系当成闭壳层单重态来算，得到的就是激发态了

倪宇晴

sobereva 发表于 2018-1-19 19:32
当前已经暗示，基态就是开壳层了。你应当绘制自旋密度，看看到底单电子是怎么分布的，搞清楚基态时候体系 ...

感谢老师，谢谢您的回复。