RUNTYP=TRANSITN

вопросы строения молекул и квантовой химии
Ответить
PTornado
Сообщения: 25
Зарегистрирован: Сб дек 23, 2006 10:38 am

RUNTYP=TRANSITN

Сообщение PTornado » Вт дек 11, 2007 1:08 pm

Здравствуйте!
Я пытаюсь воспроизвести расчет характеристик разрешенных излучательных синглет-синглетных электро-дипольных моментов переходов в молекуле О2, как это сделано у Г.И. Кобзева в "ПРИМЕНЕНИЕ НЕЭМПИРИЧЕСКИХ И ПОЛУЭМПИРИЧЕСКИХ МЕТОДОВ В КВАНТОВО-ХИМИЧЕСКИХ РАСЧЕТАХ", чисто с обучающей целью и у меня появляется ошибка:

ERROR...$CIDRT MUST USE -FZC- FOR -CI- WAVEFUNCTIONS.
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Tue Dec 11 15:00:48 2007

Хотя я все сделал вроде бы правильно. Вот сам инпут файл:
$CONTRL SCFTYP=NONE RUNTYP=TRANSITN MULT=1 cityp=guga $END
$CONTRL MAXIT=200 NZVAR=0 $END
$SYSTEM MEMORY=10000000 $END
!$BASIS GBASIS=N31 NGAUSS=6 $END
$BASIS GBASIS=TZV $END
$GUESS GUESS=moread norb=9 $END
$CIDRT1 GROUP=C1 NFZC=4 NDOC=4 NVAL=1 NALP=0 IEXCIT=2 $END
$DRT1 GROUP=C1 NMCC=4 NDOC=4 NALP=1 NVAL=0 FORS=.TRUE. $END
$TRANST NOCC=9 NFZC=4 IROOTS(1)=3 $END
$GUGDIA ITERMX=200 NSTATE=3 $END
$DATA
O2...RHF/TZV
C1
O 8.0 -0.6003326977 0.000000000 0.0000000000
O 8.0 0.6003326977 0.0000000000 0.0000000000
$END
$VEC1
1 1 3.88675988E-01 3.38755995E-01 3.83850001E-02-5.05999988E-03 6.36000012E-04
1 2 2.22199992E-03 0.00000000E+00 0.00000000E+00-1.02400000E-03 0.00000000E+00
........
$END

Вот не получается и все. Подскажите пожалуйста, где здесь ошибка.
Спасибо.

VTur
Сообщения: 7357
Зарегистрирован: Пт авг 31, 2007 1:36 pm

Сообщение VTur » Чт дек 13, 2007 1:07 am

Вы считаете возбужденное состояние бирадикала?

Nord
Сообщения: 2227
Зарегистрирован: Сб фев 14, 2004 5:36 pm

Сообщение Nord » Чт дек 13, 2007 8:45 am

Группа $DRT1 не нужна -- она нужна для runtyp=mcscf, а у Вас написано mcscf=none. Скорее все, Вы забыли "!" в самом первом символе строки
Не важно, что о вас говорят современники, важно что о вас скажут потомки

PTornado
Сообщения: 25
Зарегистрирован: Сб дек 23, 2006 10:38 am

Сообщение PTornado » Чт дек 13, 2007 12:51 pm

VTur:
Нет. Я рассматриваю возбужденное синглетное состояние молекулы O2. Все электроны спарены.

Nord:
Когда я убираю группу $DRT1 у меня также появляется ошибка:
THIS -CI- OR -MCSCF- RUN REQUIRES A FREE FORMAT $DRT1 INPUT GROUP.
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Thu Dec 13 14:41:38 2007

Как я только уже не пробовал. Кстати, у самого Кобзева, в группе $GUESS стоит norb=10. Хотя так не должно быть. Я уже изменил на norb=9.

Nord
Сообщения: 2227
Зарегистрирован: Сб фев 14, 2004 5:36 pm

Сообщение Nord » Чт дек 13, 2007 2:08 pm

Костыль: ну поставьте тогда scftyp=mcscf
Не важно, что о вас говорят современники, важно что о вас скажут потомки

PTornado
Сообщения: 25
Зарегистрирован: Сб дек 23, 2006 10:38 am

Сообщение PTornado » Чт дек 13, 2007 5:57 pm

Спасибо, но все равно не получается. Та же ошибка:

ERROR...$CIDRT MUST USE -FZC- FOR -CI- WAVEFUNCTIONS.

Буду разбираться.

Ответить

Вернуться в «квантовая химия и моделирование»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 27 гостей