Квантово-химические расчеты в GAMESS

[NadezhdaKulikova]

Сейчас попробую прикрепить.

input-output.rar

[sanya1024]

Слишком много всего не нравится программе.
(1) Рекомендую задавать молекулу в декартовых координатах и автоматически генерировать внутренние (DLC):
$CONTRL NZVAR=<любое ненулевое число> $END
$ZMAT DLC=.TRUE. AUTO=.TRUE. $END
Возможно, это избавит Вас от первой ошибки (YOUR B MATRIX HAS A SMALL DETERMINANT), а также здорово облегчит Вам оптимизацию геометрии.

Еще очень рекомендую опцию
$SCF DIRSCF=.TRUE. $END
Это тоже существенно облегчит Вам жизнь и спасет от преждевременной смерти Ваш жесткий диск.

(2) Вы зачем-то задали нестандартный радиус молекулы растворителя -- вместо зашитого в программу DMSO со след. параметрами:
EPS= 46.700; EPSINF= 2.179; RSOLV= 2.455 A; VMOL= 70.940 ML/MOL;
TCE= .98200E-01 1/K; STEN= 42.860 DYN/CM; DSTEN= 0.0000; CMF= 0.0000
Вы задали RSOLV=2.3. Возможно, это вызвало вторую ошибку (TOO MANY SPHERES, -MXSP- MUST BE LARGER). Если это не спасет, попробуйте поиграть с параметром MXSP в группе $PCM.

Вся молекула в целом выглядит несколько странной. Вы уверены в том, что делаете?

[LKalvin]

Добрый день уважаемые форумчане!

Интересует следующее : с помощью какой команды либо как можно посмотреть зашитые в GAMESS растворители а также их параметры?

P.S. прошу простить за столь дилетантский вопрос и всё же ...

Заранее Вам очень благодарен и признателен!

[NadezhdaKulikova]

Прежде всего спасибо за развернутый ответ! Насчет всего из пункта (1) - попробую разобраться. О результатах напишу чуть позже. Насчет (2) - нестандартного радиуса - это мое "творчество". Когда на предыдущей серии моих кислот программа стала выдавать TOO MANY SPHERES, -MXSP- MUST BE LARGER. - я не нашла ничего лучшего, как поиграться с радиусом, так он и остался. Особенно же меня заинтересовало последнее Ваше наблюдение. Что именно Вам кажется странным? Сама молекула, или то, как выглядит ее геометрия в моем исполнении?

[o-oxhem]

Вы молекулу оптимизировали? Или нарисовали и ладно ... Кстати, базисный набор - 6-31G, нужно минимум - 6-31G(d), желательно 6-31G(d,p) (или иной соответствующий). RHF - это недостаточно. Имеет смысл начать с B3LYP.

[NadezhdaKulikova]

Обижаете! Конечно, оптимизировала! Да, я использую убогий базис, но мне в данной задаче главное - не найти абсолютную истину, а посмотреть общую тенденцию, а уж потом сделать вывод - надо ли искать абсолютную истину...

[NadezhdaKulikova]

О растворителях в мануале "Solvation models". Там где описывается задание простой модели влияния р-ля как среды - PCM - найдете список после SOLVNT=. Выбрав любой из них и вписав в инпут, на выходе найдете его базовые характеристики.

[annafizik]

Ребята помогите составить входной файл С2H6. Я в этом деле новичок - никак не могу разобраться!

[Fanfate]

Ребята помогите составить входной файл С2H6. Я в этом деле новичок - никак не могу разобраться!

Для какого расчета файл составить? И вы в Гамессе собираетесь считать?

В аттаче мой вариант файла для Firefly, для Гамесса почти всё тоже самое, только NZVAR и d5 нужно убрать