Теоретический ИК спектр
Re: Теоретический ИК спектр
Вопрос по ИК-спектру рассчитанному в Firefly.
Почему частоты идут не по порядку?
И первые 6 надо отбрасывать?
Почему частоты идут не по порядку?
И первые 6 надо отбрасывать?
У вас нет необходимых прав для просмотра вложений в этом сообщении.
Re: Теоретический ИК спектр
Ну судя по виду этих первых 6 частот они реально смахивают на вращение и смещение.
но если честно то я удивлен что фаерфлай не выкидывает их самостоятельно а оставляет решение пользователю.
К тому же разница между 111 и реальной 127 весьма невелика. Гамесс насколько я видел вгоняет эти 6 частот в практически полный ноль и в окончательном перечне частот их уже не показывает.
Вообщем
Пусть это смотрят люди работающие в фаерфлае, их тут много.
но если честно то я удивлен что фаерфлай не выкидывает их самостоятельно а оставляет решение пользователю.
К тому же разница между 111 и реальной 127 весьма невелика. Гамесс насколько я видел вгоняет эти 6 частот в практически полный ноль и в окончательном перечне частот их уже не показывает.
Вообщем
Пусть это смотрят люди работающие в фаерфлае, их тут много.
Re: Теоретический ИК спектр
=> Гесс.
Очень хорошо, что не выкидывает. Вы правильно отметили, что одна из "нулевых" частот близка к фундаментальной. Остальные тоже не маленькие.
С этим расчетом что-то не так. Недостаточно точно локалилован минимум, недостаточная точность решения электронной задачи, большой шаг при расчете гессиана, может быть еще что-нибудь. Да, при численном расчете частот (а у Вас он именно такой) нулевые частоты могут заметно отличаться, но эти отличия не должны быть большими.
Очень хорошо, что не выкидывает. Вы правильно отметили, что одна из "нулевых" частот близка к фундаментальной. Остальные тоже не маленькие.
С этим расчетом что-то не так. Недостаточно точно локалилован минимум, недостаточная точность решения электронной задачи, большой шаг при расчете гессиана, может быть еще что-нибудь. Да, при численном расчете частот (а у Вас он именно такой) нулевые частоты могут заметно отличаться, но эти отличия не должны быть большими.
"Если миру не нужен свет - Не беда, проживёт и так.
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Re: Теоретический ИК спектр
asp - чисто на всякий случай - а геометрия оптимизирована?
Lantano - ну не выкидывает так не выкидывает, мне стоит принять это во внимание когда я доберусь до Фаерфлая, просто я несколько удивлен что фаерфлай с гамессом отличаются в этом пункте (да и вообще не так много программ кажется считают пользователя настолько умным).
- а увеличение сетки интегрирования не позволит определиться четче с тем что это за частоты?
- а что Firefly не имеет аналитических частот для ДФТ?
Lantano - ну не выкидывает так не выкидывает, мне стоит принять это во внимание когда я доберусь до Фаерфлая, просто я несколько удивлен что фаерфлай с гамессом отличаются в этом пункте (да и вообще не так много программ кажется считают пользователя настолько умным).
- а увеличение сетки интегрирования не позволит определиться четче с тем что это за частоты?
- а что Firefly не имеет аналитических частот для ДФТ?
Re: Теоретический ИК спектр
Я 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, может быть, понадобится подбирать.
Обычно при численном расчете гессиана пишу что-то такое:
$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, может быть, понадобится подбирать.
"Если миру не нужен свет - Не беда, проживёт и так.
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Если стал не нужен поэт - Не беда, он наймётся в кабак..."
Re: Теоретический ИК спектр
Геометрия оптимизирована.
NVIB=2 уменьшило эти частоты, но не исключило их полностью.
NVIB=2 уменьшило эти частоты, но не исключило их полностью.
Re: Теоретический ИК спектр
Не понятно, что значит полностью исключить
Расчет идет в масвзвешенных декартовых координатах. Их 3N штук. Колебательных степеней свободы 3N-6. Поэтому 6 вращательно-поступательных степеней будет всегда. Главное сделать эти величины много меньше частот.
Можно и занулить, но исключить?
Ну, поставьте PURIFY=.TRUE.
[ Post made via Android ]
Расчет идет в масвзвешенных декартовых координатах. Их 3N штук. Колебательных степеней свободы 3N-6. Поэтому 6 вращательно-поступательных степеней будет всегда. Главное сделать эти величины много меньше частот.
Можно и занулить, но исключить?
Ну, поставьте PURIFY=.TRUE.
[ Post made via Android ]
После отстоя требуйте долива
Re: Теоретический ИК спектр
дайте новый аутпут.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 21 гость