Оптимизация в C1, отрицательная частота

[Vit Nhoc]

Собственно сабж. Почему так бывает? Молекула CH3COOH. Выходной файл:

http://chemcraftprog.com/tmp/PkCORR.out

[Гесс]

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

А, некоторые люди непереносят "отрицательные" частоты, более корректно они все же "мнимые".

[Vit Nhoc]

Собственно сабж. Почему так бывает? Молекула CH3COOH, расчёт в wB97XD/6-31++G(D,P). Выходной файл:

http://chemcraftprog.com/tmp/PkCORR.out

Судя по координатам последних шагов оптимизации, эта молекула должна быть симметричной (Cs).

[VTur]

Используйте другой потенциал

[Vit Nhoc]

В смысле, другой функционал?

[VTur]

wB97XD измените не другой.

Вы прокалибровали метод?

[Гесс]

Зачем??? За что? Стартовав с конечной геометрии этого расчета гауссиан с очень хорошими шансами приведет оптимизацию в другом функционале в то же ПС. А если еще и guess прочитать - то шансы еще повысятся.

[VTur]

Вы сначала сделайте. Потом обсуждайте.

Функционалы бывают разные. Особенно с приставкой w.

[Гесс]

Гм. Следует признавать свои ошибки.
Я проверил эту молекулу с рядом функционалов (PBE, B3LYP, CAM-B3LYP, ωPBEhPBE, LC-ωPBE, ωB97 и ωB97XD) с тем же базисом, в D релизе, начиная с финальной геометрии расчета Vit Nhocа.
1) ωB97XD единственный функционал находящий мнимую частоту.
2) при повороте на 30° оптимизация выкручивает молекулу обратно в то же TS.
3) При запуске IRC из этого TS процедура заканчивается через один шаг,

--------------------------------------------------------------------------
Summary of reaction path following
--------------------------------------------------------------------------
Energy RxCoord
1 0.00000 0.00000
2 0.00000 0.10242
--------------------------------------------------------------------------

По видимому мы в каком то очень плоском месте ППЭ, тогда странно что такая большая частота.

[VTur]

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