- •Компьютерная технология квантово-химических расчетов с помощью программного пакета “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. Заключение
Строение электронной оболочки атома
Волновая функция электрона и стационарное уравнение Шредингера
Современная
теория строения атома базируется на
квантово-механических представлениях.
В квантовой механике, в отличие от
классической, поведение электрона в
пространстве имеет вероятностный
характер и описывается с помощью так
называемой волновой
функции
,
зависящей от декартовых координат
электрона x,y,z.
Впервые статистическая интерпретация
волновой функции была дана Максом
Борном в
1926 г.: волновая функция
представляет собой некоторую (в общем
случае комплексную) функцию, квадрат
модуля которой равен плотности вероятности
пребывания электрона в заданной точке
пространства, а вероятность dW
обнаружить электрон в элементарном
объеме dV
( dV=dxdydz)
определяется соотношением:
dW = | |2dV (1.1.1)
Кроме того, волновая функция должна быть конечной, однозначной и непрерывной во всем пространстве, а также удовлетворять условию нормировки:
(1.1.2)
Последнее соотношение означает, что вероятность обнаружить электрон во всем пространстве равна единице (вероятность достоверного события).
В 1926 году австрийский физик Эрвин Шредингер, взяв за основу идею Луи де Бройля о волновом характере микрочастиц, воспользовался волновым уравнением классической механики для описания поведения электрона в атоме. Поскольку атом представляет собой замкнутую систему, то движение электрона в атоме как волны оказалось аналогичным распространению стоячей волны. С учетом соотношения де Бройля Шредингеру удалось получить уравнение для волновой функции электрона в следующем виде:
(q,t)
(q,t)
=i
(q,t),
(1.1.3)
где
-
оператор
полной энергии
(оператор
Гамильтона
или гамильтониан),
зависящий от координат частицы и времени
и представляющий собой сумму операторов
кинетической и потенциальной энергии
(
);
, где
-
постоянная
Планка,
равная 6.6210-34
Джс.
Уравнение
(1.1.3) называется нестационарным
уравнением Шредингера. В случае, когда
гамильтониан
явно
не зависит от времени, например, при
исследовании свойств атомов или молекул
в отсутствие внешних возмущений,
зависящих от времени, то уравнение
(1.1.3) можно преобразовать в стационарное
уравнение Шредингера:
= E, (1.1.4)
где и зависят только от координат q.
Оператором называется правило или закон, согласно которому каждой функции ставится в соответствие другая функция. Важное значение имеет последовательность действия операторов на функцию. Если результат действия двух операторов не зависит от последовательности их применения, то такие операторы называются коммутирующими.
В
квантовой механике каждой физической
величине ставится в соответствие
некоторый линейный
и самосопряженный (или эрмитов)
оператор, а все соотношения между
классическими величинами заменяются
на соотношения между соответствующими
им операторами. Оператор
называется линейным, если выполняется
соотношение:
[
1(х)+
2(х)]=
1(x)+
2(x)
, (1.1.5)
где и - произвольные константы.
Cамосопряженность или эрмитовость оператора означает, что для любых двух интегрируемых функций справедливо равенство:
1(х) ( 2(х))dx=2(х) ( 1(х))dx (1.1.6)
Операторы
декартовых координат x,
y
и z
частицы представляют собой те же самые
координаты, а их действие на любую
функцию сводится к простому умножению
функции на координату (
).
Операторы проекций импульса выглядят
более сложно и определяются производными
по координатам:
(1.1.7)
а
оператор кинетической энергии электрона
можно получить, если использовать
соотношение
=
,
(1.1.8)
где
называется
оператором
Лапласа.
Решение стационарного уравнения Шредингера (1.1.4) можно рассматривать как нахождение собственных функций оператора и соответствующих ему собственных значений энергии Е .
