ПС или локальний/глобальный минимум (FF)?

[LKalvin]

Доброго времени суток уважаемые форумчане!

Производил оптимизацию геометрических параметров в FF получил результат с одной мнимой частотой, значит ли это что я получил ПС вместо локального/глобального минимума?

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

[endlesslake]

Возможно, но не обязательно. Еще это может быть артефакт из-за вычисления гессиана численным дифференцированием, особенно с использованием PCM.

[Reizen]

Самый вероятный вариант, что - да, всё так и есть, как говорит программа.

Но если структура выглядит сомнительно (у меня нет под рукой ChemCraft-а, и я не вижу, что за ПС там получилось), то можете провести дополнительные проверки:

а) поставьте PURIFY=.t. в группе $FORCE и пересчитайте гессиан (можно ради экономии времени не пересчитывать, а просто скопировать всю группу $VIB из irc-файла в inp-файл - программа сама её "подхватит"). Мнимая частота может исчезнуть, а может, впрочем, и появиться вторая... В любом случае, лучше не снимать флажок PURIFY=.t. на всех последующих шагах.

б) попробуйте ужесточить критерий поиска экстремума на ППЭ (поставьте, например, OPTTOL=0.000001 или 0.0000005 в группе $STATPT и запустите оптимизацию из вашей финальной точки). Тому количеству мнимых частот, которое получится после этого расчёта, уже можно доверять с намного большей уверенностью.

в) настоящая достоверная проверка (то, что было выше - это так, экспресс-проверка в надежде, что мнимая частота исчезнет сама собой ) - это выполнить спуск при помощи расчёта IRC (см. $IRC-группу в хэлп-файле) из получившегося ПС в обе стороны вдоль переходного вектора. Если оба спуска будут действительно спусками и приведут в те или иные минимумы - значит, строение ППЭ именно таково, каким вы его видите, а ваша структура - действительно ПС.

Это трудоёмко, да - но достоверный результат можно получить только так.

[sanya1024]

Структура в выдаче выглядит как ПС... но все не так просто.
Во-первых, очень не рекомендуется считать численные гессианы в PCM.
Во-вторых, дефолтный критерий сходимости геометрии 1e-4, а для получения корректных численных гессианов его лучше ужесточить на порядок (как выше уже сказали). Ну, и почистить (PURIFY=.t.) не повредит.
С другой стороны, мнимая частота в 117 см-1 -- это немало. Много шансов за то, что она не исчезнет.

[sanya1024]

Да, забыла добавить: обязательно ставьте $FORCE nvib=2 $end -- считаться будет вдвое дольше, зато гораздо лучше точность. Опять-таки можно подсунуть $VIB из расчета с nvib=1.

[LKalvin]

Спасибо всем кто откликнулся