如何用gaussian计算氧气单重态的能量呀？

xyyan1207

我现在试过的下面这几种都不行啊，前三种正常计算完以后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."

sobereva

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