- •Компьютерная технология квантово-химических расчетов с помощью программного пакета “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=0.95
R2=1.48
А=105.0
D=120.0
В этой записи учтена симметрия молекулы, а именно, две длины связи О-Н обозначены одной и той же переменной, поскольку в реальной молекуле они одинаковы. То же самое справедливо и для двух валентных углов Н-О-О. Если такая Z-матрица будет использована для оптимизации геометрии молекулы (в этом случае в строке опций после символа # в любом месте нужно указать Opt=Z-matrix), то значения всех переменных будут подвергнуты оптимизации для того, чтобы определить равновесную геометрию молекулы, которой соответствует наименьшая полная энергия системы.
В случае, когда необходимо выполнить частичную оптимизацию геометрии молекулы, т.е. оставить фиксированными какие-либо внутренние координаты, то в Z-матрице вслед за разделом переменных можно ввести дополнительный раздел констант (constants), например:
H
O 1 R1
O 2 R2 1 A
H 3 R1 2 A 1 D
Variables:
R1=0.95
R2=1.48
А=105.0
Constants:
D=120.0
В этом случае оптимизироваться будут переменные R1, R2 и A, а диэдральный угол D=120o останется неизменным.
Примечание: Обычно в программе Гауссиан один раздел отделяется от другого пустой строкой при использовании приведенной выше формы Z-матрицы вместо обозначений разделов Variables и Constants можно использовать пустую строку. Вместо знака равенства можно использовать один или несколько пробелов.
3.3.2 Использование декартовых координат атомов.
Несмотря на то, что использование внутренних координат атомов в общем случае является более удобным для составления Z-матрицы, иногда возникает необходимость использовать для описания положения атомов в молекуле их декартовы координаты x, y и z. В качестве примера приведем Z-матрицу молекулы Н2О, составленную в декартовых координатах для экспериментальных значений длины связи R(O-H)=0.957 Å и валентного угла Н-О-Н, равного 104.5о:
O 0. 0. 0.
H 0.7567 0.5859 0.
H -0.7567 0.5859 0.
Примечание: В этом примере за начало декартовой системы координат принят атом кислорода, координаты x, y и z которого равны 0, хотя выбор точки отсчета является произвольным (в качестве таковой может быть любой атом молекулы или совершенно произвольная точка пространства).
3.3.3 Смешанное использование внутренних и декартовых координат.
При проведении квантово-химических расчетов иногда возникает необходимость при составлении Z-матрицы использовать для описания положения атомов как внутренние, так и декартовы координаты атомов. Ниже приведена Z-матрица для оптимизации геометрии молекулы воды на поверхности ртути. В этом примере поверхность ртути моделируется небольшим фрагментом (кластером), состоящим из шестнадцати атомов Hg (рис.3.3.3.1).
Рис.3.3.3.1
Ориентация молекулы Н2О, адсорбированной на поверхности кластера Hg16
Ниже приведена Z-матрица для этой системы:
Hg 0 0. 0. 0.
Hg 0 3. 0. 0.
Hg 0 1.5 2.59808 0.
Hg 0 -1.5 2.59808 0.
Hg 0 -3. 0. 0.
Hg 0 -1.5 -2.59808 0.
Hg 0 1.5 -2.59808 0.
Hg 0 4.5 -2.59808 0.
Hg 0 0. 5.19615 0.
Hg 0 -4.5 -2.59808 0.
Hg 0 -1.5 0.86603 -2.44949
Hg 0 0. -1.73205 -2.44949
Hg 0 1.5 0.86603 -2.44949
Hg 0 0. 3.4641 -2.44949
Hg 0 -3. -1.73205 -2.44949
Hg 0 3. -1.73205 -2.44949
O 1 R1 2 U1 3 D1
H 17 R2 1 U2 2 D2
H 17 R2 1 U2 2 D3
