Квантово-химические расчеты в программе Gaussian. Серба П.В - 21 стр.

UptoLike

ВУЗ: 

ТТИ ЮФУ | Таганрог

Составители: 

Рубрика: 

80
Рис. 3.4. Сечения поверхностей потенциальной энергии вдоль координаты
реакции для реакций распада (а), обмена (б) и изомеризации (в).
В Gaussian сканирование потенциальной энергии производится
методами, изложенными в п. 3.3.
Равновесные структуры соответствующие исходным реагентам и
продуктам реакции находятся с использованием команды Opt (см. 3.1.).
Поиск седловых точек производится с использованием команды Opt
с
опциями TS или Saddle=N (N — порядок седловой точки). Эти команды
выполняются совместно с опцией CalcFC. Опция CalcFC используется для
вычисления начальных силовых констант, которые необходимы для
нахождения седловой точки. Ниже представлен пример расчета структуры
исходных реагентов, переходного состояния и продуктов реакции
.
41
H 0.51 -0.88 1.92
H 0.51 0.88 1.92
H7,C2,1.1,C1,111.2,H6,120.
H8,C2,1.1,C1,111.2,H6,-120.
В левой колонке использованы Декартовы координаты, данные справа
представляют типовую Z-матрицу. Следует заметить, что первые три атома в
пределах Z-матрицы не используют полное число параметров; только в
четвертом и последующих атомах используется полное число параметров.
Другая форма Z-матрицы для этой же молекулы:
0 1
C1
C2 C1 RCC
H3 C1 RCH C2 ACCH
H4 C1 RCH C2 ACCH H3 120.
H5 C1 RCH C2 ACCH H3 -120.
H6 C2 RCH C1 ACCH H3 180.
H7 C2 RCH C1 ACCH H6 120.
H8 C2 RCH C1 ACCH H6 -120.
Variables:
RCH = 1.5
RCC = 1.1
ACCH = 111.2
В этой Z-матрице численные значения длины связи и углов были
заменены переменными. Сами же численные значения переменных даны в
отдельном разделе после спецификации последнего атома. Разделы отделены
пустой строкой или строкой Va ri ab le s:
Периодические системы определяют с обычной спецификацией
молекулы для элементарной ячейки. Только дополнительно необходимо
ввести один, два или три
вектора трансляции, приложенные к спецификации
молекулы (без разделения пустой строкой), указывая направление
репликации.
Например, следующий спецификация определяет одномерное
вычисление для неопрена:
# PBEPBE/6-31g(d,p)/Auto SCF=Tight
neoprene, [-CH2-CH=C(Cl)-CH2-] optimized geometry
0 1
C,-1.9267226529,0.4060180273,0.0316702826
H,-2.3523143977,0.9206168644,0.9131400756
H,-1.8372739404,1.1548899113,-0.770750797
C,-0.5737182157,-0.1434584477,0.3762843235
H,-0.5015912465,-0.7653394047,1.2791284293
C,0.5790889876,0.0220081655,-0.3005160849
C,1.9237098673,-0.5258773194,0.0966261209
H,1.772234452,-1.2511397907,0.915962512
H,2.3627869487,-1.0792380182,-0.752511583

Страницы