преждевременное завершение работы GAMESS-US

[petrovps]

Добрый день, коллеги!
Нужна Ваша помощь. Столкнулись с проблемой: при расчете радикала и катион-радикала вылетает ошибка (файлы прилагаю). Молекулы считает нормально. Используем GAMESS-US, компилировал с библиотеками Intel, Linux Mint. Подскажите, пожалуйста, решение этой проблемы.
Спасибо.

[VTur]

Попробуйте убрать карту hess=guess

[petrovps]

Не помогло...

[VTur]

Пишет то же самое? Тогда геометрия слишком кривая

[sanya1024]

Прежде чем писать "кривая геометрия", было бы неплохо открыть файлы и посмотреть. Нормальная геометрия.
ТС, Вы забыли задать оптимизацию во внутренних координатах: $contrl nzvar=<любое ненулевое число> и $ZMAT DLC=.t. AUTO=.t. $END
Есть и еще огрехи, не имеющие прямого отношения к этой проблеме, но чреватые другими проблемами. Например, в $contrl надо бы добавить ispher=1 -- использование сферических d функций вместо декартовых (придется пересчитывать всю серию -- это равносильно смене базиса). Впрочем, базис тоже не повредит сменить: Вы уверены, что Вам нужен базис с d функциями на C и O, но без p функций на H?

[VTur]

Ну, аще. Специалист пришел ...

[petrovps]

Прежде чем писать "кривая геометрия", было бы неплохо открыть файлы и посмотреть. Нормальная геометрия.
ТС, Вы забыли задать оптимизацию во внутренних координатах: $contrl nzvar=<любое ненулевое число> и $ZMAT DLC=.t. AUTO=.t. $END

Спасибо. Попробую.

Есть и еще огрехи, не имеющие прямого отношения к этой проблеме, но чреватые другими проблемами. Например, в $contrl надо бы добавить ispher=1 -- использование сферических d функций вместо декартовых (придется пересчитывать всю серию -- это равносильно смене базиса). Впрочем, базис тоже не повредит сменить: Вы уверены, что Вам нужен базис с d функциями на C и O, но без p функций на H?

Насчет базиса: такой посоветовал мне коллега, а он ссылается на статью со схожим объектом. Задача простая-определить потенциал ионизации ряда родственных структур. А какой базис на Ваш взгляд лучше для данной задачи?

[amge]

Ну, аще. Специалист пришел ...

Вообще-то, sanya права, с виду геометрия хорошая. Я ее даже на всякий случай заоптимизировал Природой - сошлась за несколько шагов, почти не изменившись. Действительно хорошая, причина ошибки явно не в ней.

Молекула плоская. Может, оптимизация с Cs поможет?

[VTur]

Больше 6 собственных чисел гессиана близки к 0. Можно сменить метод расчета начального гессиана, можно расчитать начальный гессиан честно, можно чуть изменить геометрию. Автор сам решит, что проще.
Расчет в декартовых функциях лучше описывает валентные электроны. Это можно посмотреть по частотам - сравнить расчет с 5d 7f и 6d 10 f базисными функциями - в хороших сетках декартов базис лучше воспроизводит колебания водородов

[ Post made via Android ]

[sanya1024]

Ну, аще. Специалист пришел ...

Хамите.

Мой совет вполне рабочий: оптимизация в автоматически сгенерированных внутренних координатах сошлась за 5 шагов (аттач).

Что касается базиса... Чтобы лучше описать валентные электроны, нужно брать адекватный базис. Например, трехэкспонентный (хотя бы 6-311G(d,p)). p функции на водородах тоже, как ни странно, нужны -- хотя бы потому, что там возможно образование внутримолекулярной водородной связи (между двумя OH группами на фениле). Сейчас проверяю, насколько раздуется/сдуется задача при добавке p функций на водороды и при замене декартовых d функций на сферические.

Да, еще одну полезную вещь забыли: $SCF dirscf=.t. $END

[sanya1024]

Ну что: 6-31G(d) декартовы d функции -- 277 функций в базисе
6-31G(d,p) декартовы d функции -- 310 функций в базисе (тут уже dirscf оказалось критично)
6-31G(d,p) сферические d функции -- 293 функции.
Итого: с добавкой p функций и выкидыванием лишних d базис расширился, причем расширился в правильную сторону (не за счет избыточных и линейно зависимых поляризационных функций на углеродах, а за счет недостающих поляризационных функций на водородах).

VTur, а Ваш любимый Гауссиан всегда по дефолту использует сферические функции (проверьте сами). Да и большинство других программ тоже. Поэтому лучше всегда включать сферические функции, чтобы сравнение с другими расчетами было корректно.

[sanya1024]

А вот и выдачи с добавкой p функций и со сферическими d (аттач).

[petrovps]

Огромное спасибо, sanya1024! Вы очень помогли.

[VTur]

VTur, а Ваш любимый Гауссиан всегда по дефолту использует сферические функции (проверьте сами). Да и большинство других программ тоже. Поэтому лучше всегда включать сферические функции, чтобы сравнение с другими расчетами было корректно.

Да, я всегда пишу в нем 6d 10f
И в ГАМЕССЕ тоже пишу в явном виде ISPHER=.
Что бы потом, возвратившись к расчету, не просматривать весь вывод,а сориентироваться по первой странице выдачи.

[petrovps]

Добрый день!

Снова столкнулся с проблемой, не хочет считать нормально, обрывается. Файлы во вложении. Подскажите, пожалуйста, что делать.

Спасибо.

[VTur]

У Вас осциллирует роцедура SCF. Нужно изменить условия сходимости.

[petrovps]

Подскажите, пожалуйста, как изменить условия сходимости?

[sanya1024]

Итак, Вы считаете нейтральный радикал. Вы уверены, что Вам нужен именно UHF, а не более легкий вариант ROHF? Кстати, поставьте inttyp=hondo вместо pople.
Настройки сходимости SCF находятся в группе $SCF, смотрите мануал. Но то, что у Вас получилось, никакими настройками $SCF не лечится.