Делалю Гугой CASSCF(8,8) расчет геометрии возбужденного состояния в FF.
Код: Выделить всё
$CONTRL SCFTYP=MCSCF EXETYP=RUN RUNTYP=OPTIMIZE COORD=UNIQUE
MAXIT=200000 INTTYP=HONDO
ICHARG=-2 MULT=3 UNITS=ANGS ECP=READ $END
$MCSCF CISTEP=GUGA MAXIT=20000 FORS=.T. SOSCF=.T. $END
$DRT GROUP=C1 NMCC=60 NDOC=3 NALP=2 NVAL=3
FORS=.T. $END
$GUGDIA NSTATE=5 ITERMX=2000 $END
$GUGDM IROOT=1 NFLGDM(1)=1,0,0,0,0 $END
$GUGDM2 WSTATE(1)=1,0,0,0,0 $END
$GUESS GUESS=MOREAD NORB=142 $END
При расчете CASSCF(50,25) геометрии основного (и возбужденных тоже) состояния в M=1, порядок МО сохраняется (т.е. такой же как в RHF), а вот при переходе к M=3, МО поплыли.
Вот такой проблемс я думал обычно стартовые МО нужны, чтобы порядок был таким каким его задали, а тут оказывается они могут еще перемешиваться.