Теоретический ИК спектр

[asp]

Вопрос по ИК-спектру рассчитанному в Firefly.
Почему частоты идут не по порядку?
И первые 6 надо отбрасывать?

[Гесс]

Ну судя по виду этих первых 6 частот они реально смахивают на вращение и смещение.
но если честно то я удивлен что фаерфлай не выкидывает их самостоятельно а оставляет решение пользователю.
К тому же разница между 111 и реальной 127 весьма невелика. Гамесс насколько я видел вгоняет эти 6 частот в практически полный ноль и в окончательном перечне частот их уже не показывает.
Вообщем
Пусть это смотрят люди работающие в фаерфлае, их тут много.

[Lantano]

=> Гесс.
Очень хорошо, что не выкидывает. Вы правильно отметили, что одна из "нулевых" частот близка к фундаментальной. Остальные тоже не маленькие.

С этим расчетом что-то не так. Недостаточно точно локалилован минимум, недостаточная точность решения электронной задачи, большой шаг при расчете гессиана, может быть еще что-нибудь. Да, при численном расчете частот (а у Вас он именно такой) нулевые частоты могут заметно отличаться, но эти отличия не должны быть большими.

[Гесс]

asp - чисто на всякий случай - а геометрия оптимизирована?

Lantano - ну не выкидывает так не выкидывает, мне стоит принять это во внимание когда я доберусь до Фаерфлая, просто я несколько удивлен что фаерфлай с гамессом отличаются в этом пункте (да и вообще не так много программ кажется считают пользователя настолько умным).
- а увеличение сетки интегрирования не позволит определиться четче с тем что это за частоты?
- а что Firefly не имеет аналитических частот для ДФТ?

[Lantano]

Я Firefly давно ничего серьезного не считал, а уж DFT и подавно. Но любое увеличение точности электронной задачи при численном расчете гессиана можно только приветствовать. А что тут важней - увеличить точность расчета интегралов или сетку увеличить - не знаю. Я бы сделал обе вещи.
Обычно при численном расчете гессиана пишу что-то такое:
$CONTRL INTTYP=HONDO ICUT=15 ITOL=30 $END
$FORCE VIBANL=.T. NVIB=2 VIBSIZ=0.005 PURIFY=.F. RDHESS=.F. PRTIFC=.T. DECOMP=.T. $END
Самое главное во второй строке NVIB=2. А вот VIBSIZ=0.005, может быть, понадобится подбирать.

[asp]

Геометрия оптимизирована.
NVIB=2 уменьшило эти частоты, но не исключило их полностью.

[VTur]

Не понятно, что значит полностью исключить

Расчет идет в масвзвешенных декартовых координатах. Их 3N штук. Колебательных степеней свободы 3N-6. Поэтому 6 вращательно-поступательных степеней будет всегда. Главное сделать эти величины много меньше частот.
Можно и занулить, но исключить?

Ну, поставьте PURIFY=.TRUE.

[ Post made via Android ]

[Гесс]

дайте новый аутпут.