- •Компьютерная технология квантово-химических расчетов с помощью программного пакета “gaussian” Методическое пособие
- •Введение
- •Строение электронной оболочки атома
- •Волновая функция электрона и стационарное уравнение Шредингера
- •1.2 Решение уравнения Шредингера для атома водорода и водородоподобных атомов
- •1.3 Многоэлектронные атомы
- •Аналитические аппроксимации атомных орбиталей и атомные базисные наборы
- •2.1 Метод молекулярных орбиталей (мо лкао)
- •Атомные псевдопотенциалы
- •II. Теория электронного строения молекул
- •В заключение этого раздела необходимо обратить внимание еще на одно весьма важное обстоятельство.
- •2.2 Теория возмущений Мёллера-Плессе
- •2.3 Теория функционала плотности
- •2.4 Полуэмпирические методы квантовой химии
- •2.5 Модель реакционного поля для учета эффектов сольватации
- •III. Порядок работы с программным пакетом квантово-химических расчетов “ гауссиан”
- •3.1 Структура файла исходных данных
- •0 1 (Спецификация молекулы)
- •Variables:
- •3.2 Структура выходного файла
- •Variables:
- •Initialization pass.
- •Isotopes: n-14,h-1,h-1,h-1
- •Item Value Threshold Converged? Maximum Force .043925 .000450 no
- •Virtual (a1) (e) (e) (e) (e) (a1) (a1) (e) (e) (a1)
- •Ir Inten -- 381.2474 22.7331 22.7331
- •22.61138 (Kcal/Mol)
- •Vibrational temperatures: 1230.35 2673.15 2673.15 5243.52 5468.34
- •Vibrational 22.653 .605 .172
- •3.3 Алгоритм составления z-матрицы
- •3.3.1. Использование внутренних координат атомов.
- •Variables:
- •Variables:
- •3.3.2 Использование декартовых координат атомов.
- •3.3.3 Смешанное использование внутренних и декартовых координат.
- •Variables:
- •3.3.4 Использование псевдоатомов.
- •Variables:
- •IV. Оптимизация молекулярной геометрии
- •Расчет молекулы при фиксированной геометрии
- •Оптимизация геометрии
- •V. Расчет путей химических реакций
- •5.1 Поверхность потенциальной энергии реакций
- •5.2 Поиск переходного состояния.
- •Variables:
- •Variable Old X -de/dx Delta X Delta X Delta X New X
- •Item Value Threshold Converged?
- •Variable Old X -de/dx Delta X Delta X Delta X New X
- •Item Value Threshold Converged?
- •Ir Inten -- 18.4533 18.9748 128.0921
- •1 Imaginary frequencies ignored.
- •Variables:
- •Variables:
- •Iteration 8
- •I X gradient displacement newx
- •Item value threshold converged?
- •Ir Inten -- 18.5781 18.9228 128.1958
- •1 Imaginary frequencies ignored.
- •5.3 Расчет пути химической реакции.
- •Irc calculation in mass-weighted internal coordinates
- •Variables:
- •VI. Учет влияния эффектов сольватации
- •Version: matrix inversion
- •174 Tesserae over a maximum of 1500
- •VII. Дополнительные сведения
- •7.1. Порядок следования разделов в исходном файле
- •7.2. Методы расчета, реализованные в Гауссиане
- •7.3. Встроенные базисные наборы и псевдопотенциалы атомов
- •7.4. Наиболее часто употребляемые директивы
- •7.4.1 Опции директивы Opt
- •7.4.2 Опции директивы Freq
- •7.4.3 Опции директивы scrf
- •7.4.4 Опции директивы scf
- •7.4.5 Опции директивы Pop
- •7.4.6 Опции директивы Guess
- •7.4.7 Опции директивы Geom
- •7.4.8 Опции директивы irc
- •7.5. Задание базисных наборов и псевдопотенциалов во входном файле
- •VIII. Заключение
Variables:
R1 3.06
R2 0.96
U1 90.
U2 127.0
D1 269.9
D2 4.4
D3 -184.4
Положения шестнадцати атомов Hg здесь описаны их декартовыми координатами x, y и z, причем начало системы координат совпадает с положением первого атома Hg. При этом в каждой строке после символов Hg ставится ноль (0 без точки), чтобы показать, что далее следуют декартовы координаты атома. Для атома О и двух атомов Н использованы внутренние координаты, описывающие ориентацию молекулы Н2О относительно поверхности. Так для атома О указано, что его расстояние от первого атома (атома Hg) R1 составляет 3.06 Å, угол Hg(2)-Hg(1)-O(17) или U1 равен 90о (здесь цифрами в скобках указаны номера атомов в порядке их следования в Z-матрице) и диэдральный угол Hg(3)-Hg(2)-Hg(1)-O(17) или D1 равен 269.9о.
3.3.4 Использование псевдоатомов.
Многие молекулы характеризуются довольно высокой симметрией, вследствие чего некоторые внутренние координаты молекулы имеют одинаковые значения. Например, в тригонально-пирамидальной молекуле аммиака NH3 (точечная группа симметрии C3v) все три связи NH имеют одинаковую длину, одинаковое значение имеют также три валентных угла H-N-H. Чтобы при квантово-химических расчетах учесть свойства симметрии молекулы, в описании ее геометрии с помощью Z-матрицы часто бывает удобно наряду с реальными атомами молекулы использовать так называемые псевдоатомы. Псевдоатомы обозначаются символами Х или ХХ, их может быть любое количество, а располагаться они могут в любой строке Z-матрицы. Псевдоатомы используются квантово-химической программой только для описания молекулярной геометрии, а на стадии выполнения расчетов в квантово-химической процедуре игнорируются, т.е. их введение не оказывает влияния на расчетные результаты.
В случае молекулы аммиака удобно ввести псеводоатом Х, расположив его на оси симметрии молекулы, как это показано на рис.3.3.4.1. При этом удобно пронумеровать атомы, взяв в качестве первого (исходного) атома псевдоатом Х. Положение атома N(2) можно описать, задав расстояние N(2)-X(1). Значение этого расстояния может быть произвольным (например, 2.0 Å). Положение атома H(3) можно описать заданием длины связи N(2)-H(3) и плоскостного угла Х(1)-N(2)-H(3). Очевидно, что этот угол будет несколько больше тетраэдрического, поскольку валентный угол Н-N-H в молекуле аммиака несколько меньше тетраэдрического (около 107.5о). Поэтому для определенности в качестве исходного значения угла Х(1)-N(2)-H(3) можно взять угол, равный около 110.0о. Положения атомов Н(4) и Н(5) можно описать с помощью их расстояний до атома азота N(2)-H(4) и N(2)-H(5), плоскостных углов Х(1)-N(2)-H(4) и Х(1)-N(2)-H(5), а также диэдральных углов H(3)-X(1)-N(2)-H(4) и H(3)-X(1)-N(2)-H(5). Диэдральный угол H(3)-X(1)-N(2)-H(4) это угол между плоскостью, в которой расположены атомы H(3), X(1) и N(2), и плоскостью, в которой расположены атомы X(1)-N(2)-H(4). Очевидно, что эти плоскости пересекают основание правильной пирамиды (правильный треугольник) по биссектрисам углов с вершинами H(3) и Н(4). Поскольку угол между биссектрисами равен 120о, то такое же значение имеет и диэдральный угол H(3)-X(1)-N(2)-H(4). Тогда диэдральный угол H(3)-X(1)-N(2)-H(5) будет равен 240о (или -120о). В результате искомую Z-матрицу можно записать в следующем виде:
X
N 1 2.0
H 2 nh 1 xnh
H 2 nh 1 xnh 3 120.0
H 2 nh 1 xnh 3 -120.0
