Синглет с 4мя расспареными электронами

[Anna_und_Koenig]

Добрый день!
Буду очень благодарна, если кто-нибудь поможет мне разобраться с таким вопросом:
есть ли такая возможность посчитать в GAMESS синглет (или хотя бы триплет), но с распаренными по орбиталям электронами? В частности меня интересует железо.
Параметры NAOS и NBOS почему-то игнорируются. Я расчитываю MCSCF. Видимо не хватает корреляции...

[sanya1024]

Инпут в студию. И переходите на FireFly -- для MCSCF нет ничего лучше.

[Anna_und_Koenig]

Инпут.
$CONTRL SCFTYP=mcscf RUNTYP=energy MULT=1 MAXIT=500 icharg=0 $end
$GUESS GUESS=moread norb=146 $END
$SYSTEM TIMLIM=8355 MEMORY=450000000 $END
$STATPT OPTTOL=5.0E-4 $END
$SCF nconv=4 DAMP=.t. DMPCUT=1.0 vtscal=.f. $END
$MCSCF cistep=GUGA MAXIT=500 acurcy=0.0001 engtol=0.00001 $END
$GUGDIA ITERMX=800 NSTATE=4 $END
$GUGEM FASTCI=.T. PACK2=.T. $END
$DRT GROUP=C1 NMCC=135 NDOC=4 NAOS=3 NBOS=3 NVAL=1 FORS=.t. $END
$DATA
mc
C1
FE 26.0 -0.0117812333 0.0542689226 -0.1777223941
N31 6
D 1
1 1.6000000000 1.00000000
D 1
1 0.4000000000 1.00000000

O 8.0 0.0454774960 -0.1376348540 1.9080496610
N31 6
D 1
1 1.6000000000 1.00000000
D 1
1 0.4000000000 1.00000000

O 8.0 -0.6621916723 0.6391803545 2.6946055351
N31 6
D 1
1 1.6000000000 1.00000000
D 1
1 0.4000000000 1.00000000

N 7.0 -2.0644375060 0.1876099030 -0.0524799370
N31 6

N 7.0 0.1250026390 2.0904763730 0.0663038750
N31 6

N 7.0 2.0681333710 -0.1100483540 -0.0501845190
N31 6

N 7.0 -0.1302601270 -2.0344264300 -0.1712562650
N31 6

H 1.0 -3.4247354531 2.2353132384 -8.3671324355
N31 6

N 7.0 0.0361127809 0.1753203162 -2.3929934346
N31 6

O 8.0 -2.9779030854 1.7480771451 -7.6778929285
N31 6

H 1.0 -0.1215685601 -1.0422243447 -7.7995672659
N31 6

N 7.0 0.0834374257 -0.1870584527 -7.3269392151
N31 6

H 1.0 0.7188690499 0.3645297940 -7.8671347135
N31 6

C 6.0 -1.0901801295 0.5454990095 -6.8815270562
N31 6

H 1.0 -1.8161656382 -0.1564895930 -6.4936095869
N31 6

C 6.0 -1.7536758135 1.2958635524 -8.0153636754
N31 6

O 8.0 -1.2639844444 1.4965841691 -9.1040860435
N31 6

C 6.0 -0.7181113304 1.5357359398 -5.7553665242
N31 6

H 1.0 0.0879762260 2.1711688314 -6.1147282797
N31 6

H 1.0 -1.5606984251 2.1791106841 -5.5435505643
N31 6

C 6.0 -2.7812374810 1.3241562320 0.0441834260
N31 6

H 1.0 -0.9826598120 5.1941152400 0.2958723450
N31 6

H 1.0 3.3754110140 2.9918439990 0.0947184900
N31 6

H 1.0 5.1562404480 1.0235165520 0.0334830580
N31 6

H 1.0 4.9998855320 -1.6333384510 -0.0419345890
N31 6

C 6.0 4.2738921210 0.4216899880 -0.0005954620
N31 6

H 1.0 1.6661067450 5.0136255600 0.2308475060
N31 6

C 6.0 0.9538715640 4.2182805700 0.1894802050
N31 6

C 6.0 1.2952629840 2.8036089730 0.0956933890
N31 6

C 6.0 -0.3872611180 4.3099513190 0.2223214950
N31 6

C 6.0 2.5649561720 2.2896180560 0.0563395430
N31 6

C 6.0 2.9236697750 0.9283239920 -0.0011353350
N31 6

C 6.0 4.1931567540 -0.9319126280 -0.0381252230
N31 6

C 6.0 -0.9202764580 2.9535281190 0.1462399050
N31 6

C 6.0 2.7958830710 -1.2776646060 -0.0673731190
N31 6

C 6.0 2.2674900630 -2.5515420260 -0.1182409440
N31 6

C 6.0 0.9184801060 -2.9103001240 -0.1701889580
N31 6

C 6.0 0.3984988410 -4.2413004960 -0.2411285860
N31 6

C 6.0 -0.9612219510 -4.1401541680 -0.2783645270
N31 6

C 6.0 -1.2806720640 -2.7454167320 -0.2272374040
N31 6

C 6.0 -2.5807987320 -2.2013437020 -0.2170817590
N31 6

C 6.0 -2.9308038880 -0.8806705470 -0.1320899050
N31 6

C 6.0 -4.2830789760 -0.3559614430 -0.0908772650
N31 6

C 6.0 -4.1945247460 0.9855177560 0.0172590920
N31 6

C 6.0 -2.2607139470 2.6181144340 0.1324183880
N31 6

H 1.0 -2.9677272190 3.4210046490 0.2086571270
N31 6

H 1.0 -4.9929499060 1.6926613070 0.0810833890
N31 6

H 1.0 -5.1697675330 -0.9510343710 -0.1335047090
N31 6

H 1.0 -3.3856294920 -2.9105704930 -0.2627430400
N31 6

H 1.0 -1.6730222200 -4.9358578960 -0.3263881190
N31 6

H 1.0 0.9859459090 -5.1340109690 -0.2544410500
N31 6

H 1.0 2.9725461250 -3.3619034900 -0.1176189670
N31 6

H 1.0 -1.3017192717 1.8062397149 -2.7871298863
N31 6

C 6.0 -0.6445730814 1.0780585478 -3.1992781952
N31 6

C 6.0 0.7547440617 -0.5803873487 -3.1866813182
N31 6

H 1.0 1.3885736129 -1.3756798267 -2.8755136160
N31 6

C 6.0 -0.3139236390 0.8500644454 -4.4930606424
N31 6

N 7.0 0.5771669360 -0.2086791817 -4.4652000500
N31 6

H 1.0 0.9521533818 -0.6210703293 -5.2893136821
N31 6

$END
$VEC
......


Выдачи, к сожалению, при себе нет. Выдает 1ое состояние дубл.-дубл. пару, 2 неспаренных спина.

[sanya1024]

При fors=.t. включается CASSCF, т.е., генерируются все правильные по спину комбинации на данном кол-ве электронов и орбиталей.
Ваш инпут заказывает расчет только одного состояния, низшего с mult=1. А Вам надо считать сразу несколько состояний и делать усреднение эл. плотности по нужным из них (обычно это несколько низших подряд). Для этого нужно добавить
$GUGDIA NSTATE=<кол-во заказываемых состояний с некоторым запасом> prttol=0.01 $END
$GUGDM2 wstate(1)=<массив из весов усредняемых состояний> $END

Обычно в wstate(1) все веса берут равными 1.0, т.е., пишут столько единичек, по скольким состояниям усредняем.

Сразу говорю: в GAMESS-US Ваша задача гарантированно не влезет, а понадобится делать много проб, чтобы определиться с кол-вом состояний и усреднением (а возможно, и с акт. пространством). Переход на FireFly для этой задачи -- жизненная необходимость.

[Anna_und_Koenig]

Я считаю в FireFly, написала GAMESS автоматически.
Я правильно поняла, вы предлагаете считать aldet? Я пробовала, он, к сожалению, совсем уводит от СПЗ и дает состояние 5Fe-3O2-

[sanya1024]

Совершенно не обязательно. То, что я написала -- это как раз для guga. Я предлагаю считать больше одного состояния и делать усреднение по состояниям с таким расчетом, чтобы в усреднение попали состояния, где доминируют нужные Вам конфигурации.

Итак, Вы хотите синглет, но со всеми распаренными электронами. Строго говоря, это невозможно, поскольку такой слейтеровский детерминант не является чистым по спину, а нужна некая линейная комбинация детерминантов такого типа -- конфигурационная функция состояния (configuration state function, CSF). Поэтому здесь лучше вариант guga, работающий с конфигурациями. Однако состояние в многоэлектронной системе -- штука еще более сложная, и на выходе CASSCF Вы получите состояния, являющиеся линейной комбинацией конфигураций. Единственное ограничение -- ограничение по спину (guga).
Судя по всему, желаемое Вами заполнение орбиталей достигается не в низшем синглетном состоянии. Поэтому и нужно заказать расчет нескольких состояний. А чтобы не только нижнее, но и следующие состояния имели оптимальные для себя орбитали, нужно делать усреднение эл. плотности по этим состояниям.

Думаю, надо посмотреть на Ваши стартовые орбитали (покажите выдачу, из к-рой Вы вставляли орбитали в этот инпут в виде $VEC) и прикинуть, что тут можно сделать.