Спин-орбитальное взаимодействие для O2

[Shurick]

Привет посетителям форума!

Наверное мне следует обратиться на форум Грановского, но спрошу пока здесь, по-русски.
Мне нужно посчитать матричный элемент спин-орбитального взаимодействия между низшими синглетным и триплетным термами системы двух атомов кислорода, разведённых на произвольное расстояние.

Вопрос 1. Можно ли это посчитать в GAUSSIANe? Как я понял, в Гауссиане можно посчитать SOC только между состояниями одной мультиплетности. Прав ли я?

Вопрос 2. Возможно ли это правильно сделать в Firefly? Причина моих сомнений вот в чём. Я взял MO из усреднённого по состояниям MCSCF расчёта в базисе ACC-PVDZ и подставил их по $VEC1 при расчёте RUNTYP=SPINORBT в том же базисе. Входной файл выглядит у меня при этом что-то навроде

$CONTRL SCFTYP=MCSCF RUNTYP=SPINORBT CITYP=GUGA MULT=3 $END
$TRANST METHOD=BREIT NFZC=3 NOCC=9 NUMVEC=1 NUMCI=2 IROOTS(1)=3,1 $END
$DRT GROUP=C1 NMCC=3 NDOC=5 NVAL=1 FORS=.TRUE. $END
$CIDRT1 GROUP=C1 FORS=.T. NFZC=3 NDOC=5 NVAL=1 $END
$CIDRT2 GROUP=C1 FORS=.T. NFZC=3 NDOC=4 NALP=2 $END
$BASIS EXTFIL=.T. GBASIS=ACC-PVDZ NDFUNC=3 DIFFSP=.T. $END
$GUESS GUESS=MOREAD $END
$DATA
O2
C1
O 8 -0.677128746 0.042301867 0
O 8 0.595939753 -0.120058286 0
$END

Однако результаты меня удивили. Как их интерпретировать?

SPIN-ORBIT COUPLING CONSTANTS (CM-1)

----------------------------------------------------
SINGLET STATE | TRIPLET STATE | SOCC(1E) | SOCC(2E)

----------------------------------------------------
1 1 0.001 0.000
2 1 0.001 0.000
3 1 252.582 160.704

Неужели матричный эл-т СОВ для 1-го синглетного (a1Dg?) и 1-го триплетного (Sigma_g-) терма равен нулю?

Зато межу 3-им синглетом и 1-м триплетом аж 160 cm-1 (целых 0.02 эВ!)!
Наверное, я что-то недопонял.
Пожалуйста, помогите.

[sanya1024]

Я бы советовала делать спин-орбитальные расчеты в GAMESS-US -- там гораздо более развитые средства для этого. Инпут придется слегка модифицировать (см. мануал Гамесса).
Что касается Ваших результатов -- какие щели между состояниями? если состояния разделены большой щелью, SOC будет маленьким, и наоборот.

Можно на выдачу посмотреть?

[Shurick]

Во вложении out-файл.

[sanya1024]

Действительно странно: единственный триплет лежит ниже первых двух вырожденных синглетов на 0.68 эВ (с такой щелью взаимодействие будет заметным) и ниже третьего синглета на 1.13 эВ. Почему-то в выдаче вижу матричные элементы HSO типа <Singlet1,MS= 0|HSOTriplet1,MS= 1>, <Singlet2,MS= 0|HSOTriplet1,MS= 1>, <Singlet3,MS= 0|HSOTriplet1,MS= 1> и <Singlet1,MS= 0|HSOTriplet1,MS= 0>, <Singlet2,MS= 0|HSOTriplet1,MS= 0>, <Singlet3,MS= 0|HSOTriplet1,MS= 0>, но не вижу <Singlet1,MS= 0|HSOTriplet1,MS= -1>, <Singlet2,MS= 0|HSOTriplet1,MS= -1>, <Singlet3,MS= 0|HSOTriplet1,MS= -1> -- странно.
Может, выложите заодно еще инпут (с $VEC1 группой)? я его Гамессом прогоню, минутное дело. Мне кажется, что в Вашем инпуте чего-то не хватает, но в мануале FF эта часть описана довольно фигово, а в Гамессе со времен "разветвления" прибавилось много чего в т.ч. в части спин-орбиты.

[Shurick]

Выкладываю инпут-файл

[sanya1024]

Ловите готовенького. Анализировать потом будем.