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

вопросы строения молекул и квантовой химии
Ответить
LKalvin
Сообщения: 389
Зарегистрирован: Ср ноя 17, 2010 10:01 pm

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

Сообщение LKalvin » Сб окт 13, 2012 8:10 am

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

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

Заранее Вам очень признателен и благодарен!
У вас нет необходимых прав для просмотра вложений в этом сообщении.

endlesslake
Сообщения: 136
Зарегистрирован: Чт авг 09, 2012 9:13 pm

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

Сообщение endlesslake » Сб окт 13, 2012 1:33 pm

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

Reizen
Сообщения: 69
Зарегистрирован: Пт ноя 21, 2008 12:54 pm

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

Сообщение Reizen » Сб окт 13, 2012 2:50 pm

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

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

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

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

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

Это трудоёмко, да - но достоверный результат можно получить только так.
...Reveal your face to me, guide me through the Stygian fields,
Enthral my soul to Sepedet's beams to serve your will!

(Nightwish - "The Pharaoh Sails To Orion")

Аватара пользователя
sanya1024
Сообщения: 1659
Зарегистрирован: Чт янв 20, 2011 3:24 pm

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

Сообщение sanya1024 » Сб окт 13, 2012 4:51 pm

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

Аватара пользователя
sanya1024
Сообщения: 1659
Зарегистрирован: Чт янв 20, 2011 3:24 pm

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

Сообщение sanya1024 » Сб окт 13, 2012 5:17 pm

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

LKalvin
Сообщения: 389
Зарегистрирован: Ср ноя 17, 2010 10:01 pm

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

Сообщение LKalvin » Сб окт 13, 2012 5:27 pm

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

Ответить

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

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

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