Gaussian: MP4 - ИК/КР частоты
-
- Сообщения: 426
- Зарегистрирован: Чт май 24, 2007 1:54 pm
Nord, динамикой-то занимаюсь (уж много там или нет), но немного в другом аспекте. Реализация довольно простая - выписываем каждые n шагов дипольный момент элементарной ячейки (в случае ab initio - надлежащим образом разобравшись с Berry phase с помощью MLWFC), открываем файлик в программе xmgrace и делаем преобразование Фурье
Что касается интегралов по траекториям - этим занимается в Бохуме Dominik Marx (CPMD), соответственно, Маркса и нужно читать, начиная примерно отсюда: http://www.theochem.ruhr-uni-bochum.de/ ... ex.en.html. Там же, полагаю, и вообще про ИК-спектры должно найтись. Насколько я знаю, пару лет назад они уже получали очень хорошие спектры для протонированного метана
Что касается интегралов по траекториям - этим занимается в Бохуме Dominik Marx (CPMD), соответственно, Маркса и нужно читать, начиная примерно отсюда: http://www.theochem.ruhr-uni-bochum.de/ ... ex.en.html. Там же, полагаю, и вообще про ИК-спектры должно найтись. Насколько я знаю, пару лет назад они уже получали очень хорошие спектры для протонированного метана
Плохо зная грамматику, сложные конструкции должны употребляться с осторожностью.
Это-то я и сам знаю. Более того, занимаюсь этим в квантовом случае. Но хотелось бы видеть последовательный вывод для классической механики, с учетом статистической природы состояния [/quote]Darth Vasya писал(а):Nord, динамикой-то занимаюсь (уж много там или нет), но немного в другом аспекте. Реализация довольно простая - выписываем каждые n шагов дипольный момент элементарной ячейки (в случае ab initio - надлежащим образом разобравшись с Berry phase с помощью MLWFC), открываем файлик в программе xmgrace и делаем преобразование Фурье
Спасибо!Darth Vasya писал(а):Что касается интегралов по траекториям - этим занимается в Бохуме Dominik Marx (CPMD), соответственно, Маркса и нужно читать, начиная примерно отсюда: http://www.theochem.ruhr-uni-bochum.de/ ... ex.en.html. Там же, полагаю, и вообще про ИК-спектры должно найтись. Насколько я знаю, пару лет назад они уже получали очень хорошие спектры для протонированного метана
Не важно, что о вас говорят современники, важно что о вас скажут потомки
Здесь слишком тесно. Гармонические частоты могут быть близкими к экспериментальным по величине, но иметь разный тип симметрии и отнесение - случайное совпадение по величине разных частот. Вообще желательно заранее прикидывать какие частоты должны получаться. Это легко сделать для характеристических частот - есть таблицы, для остальных нужно смотреть известные частоты гомологов.1177 писал(а):A 0-800 см-1?Сравнивать гармонический и экспериментальный спектры можно только в области 800 - 1500 см-1, здесь они близки (если базис с диффузными и поляризационными функциями и не чистый ХФ). С простым ХФ спектр везде плохой (есть случайности, но заранее не угадаешь).
Уважаемый 1177, заранее извиняюсь за чисто ликбезный вопрос. Вы учитываете, что 4,11-диметилтетрадекан в стереохимическом плане близкий родственник винной кислоты и имеет три изомера: мезо-форму и пару энантиомеров. Это может оказаться важно, поскольку экспериментальный спектр, скорее всего, снят (или будет снят) для смеси всех изомеров. Хотя, возможно, почти все различия и потонут в погрешностях расчёта и эксперимента.
Вообще же, IMHO, пример для заявленной цели
Вообще же, IMHO, пример для заявленной цели
как минимум весьма трудный. Поэтому, если задача носит существенно учебный характер, возможно, стоит подумать о других объектах, которые ранее подвергались детальному исследованию и расчёт которых не будет связан с большими затратами времени.1177 писал(а):Собственно, цель в том, чтобы сравнить результаты различных методов с экспериментом
Если люди не полагают, что математика проста, то только потому, что они не понимают, как на самом деле сложна жизнь (Джон фон Нейман)
1. GAMESS, не PC1177 писал(а):Перефразирую вопрос:
Есть ли стандартные программы (пакеты), способные провести расчет частот НЕ в гармоническом приближении? спасибо
2. GAUSSIAN
3. DALTON
http://www.kjemi.uio.no/software/dalton/dalton.html
4. MULTIMODE
http://www.chemistry.emory.edu/faculty/ ... multimode/
5. NWchem
http://www.emsl.pnl.gov/docs/nwchem/doc/user/index.html
и др.
Частоты, полученные в разных программах могут значительно отличаться.
Opt(tight) SCF(tight) Pop=none Freq(anharm,Noraman).1177 писал(а): А как это сделать в GAUSSIAN?
NoRaman - чтобы не считать КР интенсивности и поляризации. Если требуются, то нужно явно указать Freq(anharm,Raman).
GAUSSIAN какой версии? В ранних (и краденых) NoRaman не срабатывает и расчет удваивается по времени. Чем лучше соптимизируете геометрию, тем достоверней спектр.
Правильно говорят: "не говори ГОП..." При расчете n-пентана выдал ошибку:
Asymmetric Top with C1 Symmetry
Inertia Moments (amu-Ang**2) Ia: 30.30814 Ib: 251.57727 Ic: 265.45285
Rotational Constants (cm-1) a: 0.55621 b: 0.06701 c: 0.06351
Rotational Constants (GHz) a: 16.67469 b: 2.00884 c: 1.90384
VibMW allocation failure: iend,mxcore= 13533003 6251593
Error termination via Lnk1e in C:\G03W\l716.exe at Wed Nov 28 22:32:10 2007.
Геометрию чем оптимизировали?1177 писал(а):Правильно говорят: "не говори ГОП..." При расчете n-пентана выдал ошибку:
Asymmetric Top with C1 Symmetry
Inertia Moments (amu-Ang**2) Ia: 30.30814 Ib: 251.57727 Ic: 265.45285
Rotational Constants (cm-1) a: 0.55621 b: 0.06701 c: 0.06351
Rotational Constants (GHz) a: 16.67469 b: 2.00884 c: 1.90384
VibMW allocation failure: iend,mxcore= 13533003 6251593
Error termination via Lnk1e in C:\G03W\l716.exe at Wed Nov 28 22:32:10 2007.
Если каким-то боком задействую теорию возмущений, то к месту и не к месту ставлю Density=Current. Но здесь, думаю, не поможет.
Скинте *.gjf или лучше *.out файл.
У Вас, видимо, мало оперативки. Какова конфигурация компа? Всегда пишите в явном виде %mem= "max" и в # maxdisk="max".1177 писал(а):C5_anahrm.out # RHF/6-31G(d) Opt(tight) SCF(tight) Freq(anharm)
В RHF/6-31G(d) частоты будут плохие. Можно было изначально задать правильную симметрию, ответ чуть улучшиться, но не намного.
В B3LYP/6-311++G(3df,3pd) ангармонические частоты н-пентана (всё транс) я получил за полгода. Самое интересное, в низкочастотной области (< 1000 см-1), гармонические и ангармонические частоты довольно близки. Разность не превышает стандартных 25 см-1. Так, что если интересует эта область, то можно не мучиться и взять гармонические.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 56 гостей