Добрый день!
Пример: 2-Ме-1,3-диоксолан.
ПРИРОДА PBE/3z:
eng>$Energy
eng> E=-3.074205298469e+02
eng> D= 3.88164395e-01 9.96774290e-02 2.48075939e-01
или же
X = 0.9866 Y = 0.2534 Z = 0.6305 (Debye)
Gaussian B3LYP/6-31G(d,p):
X = 1.0694 Y = 0.0001 Z = 0.5670 (Debye).
--
То, что выдаёт Gaussian - это компоненты дипольного момента вдоль главных осей инерции. Это величины, которые, в принципе, можно измерить с помощью микроволновой спектроскопии. Так для второго примера (2-Ме-1,3-диоксолан) экспериментально получены следующие компоненты:
X = 1.055 Y = 0.000 Z = 0.624
А теперь вопросы: в какой системе координат выдаёт компоненты дипольного момента программа ПРИРОДА? Если это система главных осей инерции, то почему такая разница с Gaussian? Если это не система главных осей, то как перевести посчитанные компоненты дипольного момента в эту систему?
Компоненты дипольного момента при расчёте в ПРИРОДА
Re: Компоненты дипольного момента при расчёте в ПРИРОДА
Ответ от автора: "Дипольный момент выдаётся в той же системе координат, как и всё остальное (G=... H=... ), то есть исходной, неповёрнутой."
Кто знает основы механики (что такое тензор инерции) и линейной алгебры (собственные векторы матрицы 3x3, линейное преобразование)?
Кто знает основы механики (что такое тензор инерции) и линейной алгебры (собственные векторы матрицы 3x3, линейное преобразование)?
Re: Компоненты дипольного момента при расчёте в ПРИРОДА
У меня есть скрипт, который может переориентировать молекулу по главным осям эллипсоида инерции.
xyz -i -ellips=2 file.xyz
См. также комментарии к подпрограмме ellips этого скрипта.
xyz -i -ellips=2 file.xyz
См. также комментарии к подпрограмме ellips этого скрипта.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 19 гостей