计算化学公社

标题: 使用stable=opt之后，log文件最后一段出现ub3lyp。 [打印本页]

作者
Author: 倪宇晴 时间: 2018-1-18 22:41
标题: 使用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。请各位老师指点迷津。

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

作者
Author: 倪宇晴 时间: 2018-1-18 23:00

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

谢谢老师答复

作者
Author: 倪宇晴 时间: 2018-1-19 18:44

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

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

作者
Author: sobereva 时间: 2018-1-19 19:32

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

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

作者
Author: 倪宇晴 时间: 2018-1-19 19:54

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)