计算化学公社

标题: 如何用gaussian计算氧气单重态的能量呀？ [打印本页]

作者
Author: xyyan1207 时间: 2018-11-16 19:38
标题: 如何用gaussian计算氧气单重态的能量呀？
我现在试过的下面这几种都不行啊，前三种正常计算完以后spin population还都是闭壳层的0和0，不是1和-1；最后一种有个报错。麻烦各路大神帮忙看看，这个氧气单重态的能量到底怎么算呢？？？

第一种：
%chk=o2-m1
#p  ub3lyp/def2TZVPP  stable=opt

Title Card Required

0 1
O             0.00000000 0.00000000 0.60729600
O             0.00000000 0.00000000 -0.60729600

第二种：
%chk=o2-m1
#p  ub3lyp/def2TZVPP  opt guess=mix

Title Card Required

0 1
O             0.00000000 0.00000000 0.60729600
O             0.00000000 0.00000000 -0.60729600

第三种：
%chk=o2-m1
#p  ub3lyp/def2TZVPP  opt

Title Card Required

0 1
O             0.00000000 0.00000000 0.60729600
O             0.00000000 0.00000000 -0.60729600

以上三种情况：计算完以后spin population还都是闭壳层的0和0，不是1和-1

第四种：
%chk=FragGuess
#p  ub3lyp/def2TZVPP  guess(Fragment=2,only) pop=none

Title Card Required

0 1 0 2 0 -2
O(Fragment=1)             0.00000000 0.00000000 0.60729600
O(Fragment=2)             0.00000000 0.00000000 -0.60729600

--Link1--
%chk=FragGuess
%mem=8GB
%nosave
#P ub3lyp/def2TZVPP/Auto opt Guess=Read Geom=AllCheck ...

这种情况会报错，报错是这样的  "The combination of multiplicity 2 and  8 electrons is impossible."

作者
Author: sobereva 时间: 2018-11-18 19:33
如果你说的1、-1是指的原子自旋布居，根本不可能出现那种状态
非要用对称破缺方式算氧气的单重态，那是两个整体轨道都是单占据的状态
然而靠UKS计算氧气的单重态能级是很不准的，应当用多参考方法

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