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

вопросы строения молекул и квантовой химии
asp
Сообщения: 122
Зарегистрирован: Вт ноя 11, 2003 11:34 am

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

Сообщение asp » Сб мар 07, 2015 8:42 am

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

Аватара пользователя
Гесс
Сообщения: 9856
Зарегистрирован: Ср фев 15, 2012 11:19 pm

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

Сообщение Гесс » Сб мар 07, 2015 2:09 pm

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

Аватара пользователя
Lantano
Сообщения: 162
Зарегистрирован: Чт окт 09, 2008 10:35 pm

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

Сообщение Lantano » Сб мар 07, 2015 7:18 pm

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

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

Аватара пользователя
Гесс
Сообщения: 9856
Зарегистрирован: Ср фев 15, 2012 11:19 pm

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

Сообщение Гесс » Сб мар 07, 2015 7:47 pm

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

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

Аватара пользователя
Lantano
Сообщения: 162
Зарегистрирован: Чт окт 09, 2008 10:35 pm

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

Сообщение Lantano » Сб мар 07, 2015 8:38 pm

Я 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
Сообщения: 122
Зарегистрирован: Вт ноя 11, 2003 11:34 am

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

Сообщение asp » Пн мар 09, 2015 8:57 pm

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

VTur
Сообщения: 7360
Зарегистрирован: Пт авг 31, 2007 1:36 pm

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

Сообщение VTur » Пн мар 09, 2015 9:43 pm

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

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

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

[ Post made via Android ] Изображение
После отстоя требуйте долива

Аватара пользователя
Гесс
Сообщения: 9856
Зарегистрирован: Ср фев 15, 2012 11:19 pm

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

Сообщение Гесс » Вт мар 10, 2015 2:30 am

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

Ответить

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

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

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