Это-то я и сам знаю. Более того, занимаюсь этим в квантовом случае. Но хотелось бы видеть последовательный вывод для классической механики, с учетом статистической природы состоянияDarth Vasya писал(а):Nord, динамикой-то занимаюсь (уж много там или нет), но немного в другом аспекте. Реализация довольно простая - выписываем каждые n шагов дипольный момент элементарной ячейки (в случае ab initio - надлежащим образом разобравшись с Berry phase с помощью MLWFC), открываем файлик в программе xmgrace и делаем преобразование Фурье
[/quote]Спасибо!Darth Vasya писал(а):Что касается интегралов по траекториям - этим занимается в Бохуме Dominik Marx (CPMD), соответственно, Маркса и нужно читать, начиная примерно отсюда: http://www.theochem.ruhr-uni-bochum.de/ ... ex.en.html. Там же, полагаю, и вообще про ИК-спектры должно найтись. Насколько я знаю, пару лет назад они уже получали очень хорошие спектры для протонированного метана
Здесь слишком тесно. Гармонические частоты могут быть близкими к экспериментальным по величине, но иметь разный тип симметрии и отнесение - случайное совпадение по величине разных частот. Вообще желательно заранее прикидывать какие частоты должны получаться. Это легко сделать для характеристических частот - есть таблицы, для остальных нужно смотреть известные частоты гомологов.1177 писал(а):A 0-800 см-1?Сравнивать гармонический и экспериментальный спектры можно только в области 800 - 1500 см-1, здесь они близки (если базис с диффузными и поляризационными функциями и не чистый ХФ). С простым ХФ спектр везде плохой (есть случайности, но заранее не угадаешь).
как минимум весьма трудный. Поэтому, если задача носит существенно учебный характер, возможно, стоит подумать о других объектах, которые ранее подвергались детальному исследованию и расчёт которых не будет связан с большими затратами времени.1177 писал(а):Собственно, цель в том, чтобы сравнить результаты различных методов с экспериментом
1. GAMESS, не PC1177 писал(а):Перефразирую вопрос:
Есть ли стандартные программы (пакеты), способные провести расчет частот НЕ в гармоническом приближении? спасибо
Opt(tight) SCF(tight) Pop=none Freq(anharm,Noraman).1177 писал(а): А как это сделать в GAUSSIAN?
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.
У Вас, видимо, мало оперативки. Какова конфигурация компа? Всегда пишите в явном виде %mem= "max" и в # maxdisk="max".1177 писал(а):C5_anahrm.out # RHF/6-31G(d) Opt(tight) SCF(tight) Freq(anharm)
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 23 гостя
