CASSCF в GAMESS. Бен ай нид хелп.

[Vinianel]

Уважаемые коллеги, добрый день!

Имею проблему с CASSCF в Gamess, хотелось бы разобраться, что я делаю не так. Зараннее хочу предупредить что пока очень плохо разбираюсь в теме CASSCF как таковой, в голове каша.
Очень хотелось бы услышать мнение работающих непосредственно в Gamess, а так же мнение уважаемой sanya1024.

Итак, пытаюсь сделать CASSCF расчет для комплекса с кобальтом. Имеется геометрия, оптимизированная в DFT.
Сперва получила набор стартовых орбиталей из ХФ расчета. Посмотрела на них в Chemcraft-e, определилась, какие хочу включить в активное пространство. Количество активных электронов определила по заселенностям орбиталей.

Дальше задаю CASSCF расчет, шапочка инпута:

Код: Выделить всё

$CONTRL SCFTYP=MCSCF MAXIT=20 COORD=ZMT
ICHARG=1 MULT=1 RUNTYP=ENERGY $END
$SYSTEM TIMLIM=10000 MWORDS=1000 $END
$BASIS EXTFIL=.T. GBASIS=TZVP $END

 $MCSCF CISTEP=ALDET $END
 $DET GROUP=C1 ISTSYM=1 NCORE=264 NACT=12
 NELS=8 ITERMX=5 $END
 $TRANS DIRTRF=.T. $END
 $ciinp castrf=.t. $end
 $GUESS GUESS=MOREAD NORDER=1 IORDER(272)=263
 IORDER(263)=272 NORB=2950 $END

Потом идет геометрия молекулы и массив $VEC из вышеупомянутого ХФ расчета.

Программа выдает ошибку:

Код: Выделить всё

 ***** ERROR *****
 THIS RUN FOUND 0 CI EIGENVECTORS WITH S= 0.00,
 BUT YOU REQUESTED STATE AVERAGING OF 1 ROOTS.
 PLEASE EXAMINE YOUR CHOICE OF -NSTATE- INPUT DATA.

В чем может быть проблема?
Я наткнулась на похожее обсуждение вот в этой ветке: viewtopic.php?f=71&t=101684&sid=a1e6558 ... 169f4c3feb

На всякий случай перепроверила, является ли синглет в моем случае более стабильным. DFT рассчет показывает, что да.

Может быть я неправильно задаю активное пространство?
Буду признательна за любой совет.

По поводу того, почему считаю GAMESS-ом, а не Firefly (так как скорее всего тема в любом случае всплывет при обсуждении). Дело в том, что наш cluster support отказался иметь дело с Firefly по причине его 32-битности. Мне удалось установить Fireflу на старенький домашний кластер, однако результатов по моим комплексам там не дождаться.

[Tonnie Fox]

у Вас не хватает NSTATE=, WSTATE(1)=, iroot= (в случае SS-CASSCF расчетов) в блоке $DET
и ntrack=, istate=(в случае SA-CASSCF расчетов) в блоке $MCSCF

именно на это и ругается программа

[Vinianel]

Хм. Поправьте меня, если я не права, но я думала, что в SS-CASSCF в случае, если указанные вами параметры не заданы, по умолчанию присваетвается NSTATE=1, WSTATE(1)=1,0 и IROOT=1. Т.е. то что и нужно для синглета.

[Tonnie Fox]

Хм, да. Честно говоря, никогда не пытался не задавать нужные мне параметры, поэтому это первое что пришло в голову.
Попробуйте в $DET добавить PURES=.f. и сгененировать стартовый набор состояний при NSTATE=9 для вашей системы (где будут и синглеты и триплеты). Там уже найдете нужное Вам.

[Vinianel]

О, спасибо большое, попробую!

[Vinianel]

Эм. Попробовала.
Программа начала считать 1 electron integrals и выдала ошибку:

Код: Выделить всё

INITIAL GUESS ORBITALS GENERATED BY MOREAD   ROUTINE.
 *** ERROR IN READMO:   PROBLEM READING ORBITALS!
 POSSIBLY A DAMAGED OR MANGLED ORBITAL INPUT GROUP?
 ERROR OCCURED AT ORBITAL=     2 (MODULUS 100=   2),
          ITS LINE NUMBER=     1
 DATA READ FROM INPUT WAS ORBITAL=     1 LINE=   288
 EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Tue Jun  6 08:09:14 2017
             10263420  WORDS OF DYNAMIC MEMORY USED
 CPU     0: STEP CPU TIME=    34.78 TOTAL CPU TIME=      169.9 (    2.8 MIN)
 TOTAL WALL CLOCK TIME=      172.1 SECONDS, CPU UTILIZATION IS  98.76%
 DDI Process 0: error code 911
[cli_0]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 911) - process 0

[Tonnie Fox]

дайте аут файл dft расчетов комплекса и punch файл от него. так проще будет, чем угадывать проблему

[Vinianel]

Аутпут и punch файл прикрепила.

[Tonnie Fox]

Да вроде всё нормально с Вашим комплексом. Попробовал в убогом базисе на домашнем пылесосе, орбитали из hf расчета, остальное почти всё как у Вас (акт. пространство поменьше).
Перепроверьте стартовые орбитали, чтобы соответствовали оптимизированной геометрии, и геометрию в дефолтных для гамеса координатах (COORD=UNIQUE). Должно сработать.

[Vinianel]

Спасибо большое, попробую.