Молекула водорода

Моле́кула водоро́да — простейшая молекула, состоящая из двух атомов водорода. В её состав входят два ядра атомов водорода и два электрона. Вследствие взаимодействия между электронами и ядрами образуется ковалентная химическая связь. Кроме основной изотопной модификации H₂, существуют разновидности, в которых один или оба атома протия заменены другими изотопами водорода  — дейтерием или тритием: HD, HT, D₂, DT, T₂. Симметричность или несимметричность молекулы имеет значение при её вращении.

Модели молекулы водорода

Рис.1.Модель водорода по Бору и осевая проекция электронного кольца

Н.Бор ещё в 1913 году дал описание модели молекулы водорода^[1]. Согласно Бору, два внешних электрона, образующих молекулу, вращаются по одной и той же орбите вокруг линии, проходящей через ядра обоих атомов, и удерживают последние на определённом расстоянии друг от друга (рис.1). Боровская модель химической связи давала чёткую картину образования молекулярного водорода — динамическое равновесие системы, содержащей два протона, удерживающихся на определённом расстоянии друг от друга (d), притяжением кольца из двух электронов. При этом модель Н.Бора отражала кулоновскую электронную корреляцию — оба электрона за счёт электростатического отталкивания находятся в диаметрально противоположных точках электронного кольца вокруг линии, проходящей через ядра обоих атомов молекулы водорода.

Также как и боровская модель атома, боровская модель молекулы водорода не отражала и не могла отражать волновую природу электрона и статистическую интерпретацию волновой функции. Здесь уместно дополнить боровскую модель молекулы водорода достижениями квантовой механики, аналогично боровской модели атома. В этом случае гипотетическое кольцо электронов в модели молекулы водорода трансформируется в геометрическое место расположения плотности вероятности нахождения связывающих электронов. Электронную пару следует представлять электронным облаком в форме кольца, баранки, тора. При этом предполагается, что электроны в торе удалены друг от друга на максимально возможное расстояние, а электронная плотность в молекуле водорода имеет минимум — «кулоновскую дырку» на линии, соединяющей ядра молекулы водорода. Наиболее вероятный радиус электронного тора (r_e) определяется боровским радиусом a₀ и межъядерным расстоянием (d)

r_e² = a₀² — (d/2)²; r_e = 0,377 Å.

Рис.2.Контурная карта электронной плотности для молекулы водорода^[2]

В теории молекулярных орбиталей, разработанной в 1927—1932 г.г. Ф.Хундом, Р.Малликеном и др., модель молекулы водорода описывали контурными картами электронной плотности (рис.2).

Карта строится как разность между электронной плотностью молекулы и несвязанных атомов. Сплошные линии контурной карты отражают области увеличения электронной плотности (в единицах заряда электрона). Пунктирные — области её уменьшения в молекуле по сравнению с электронной плотностью несвязанных атомов водорода, сближенных на равновесное расстояние.^[2]

Поскольку в теории молекулярных орбиталей учитывается лишь фермиевская корреляция электронов, поэтому электронная плотность в молекуле водорода сконцентрирована на линии, соединяющей ядра молекулы (см. рис.2).

Рис.3.Модель молекулы водорода в теории валентных связей

Квантово-механическую теорию ковалентной связи в молекуле водорода разработали в 1927 году Вальтер Гайтлер и Фриц Лондон. Выполненный ими расчёт волновых функций явился основой теории валентных связей. В этой теории модель молекулы водорода описывали перекрыванием атомных орбиталей атомов водорода (1S типа). Схема перекрывания электронных облаков, образующих сигма-связь молекулы водорода, совпадает с линией соединения атомных центров молекулы (рис.3).

Расчёт молекулы водорода в основном сводится к определению интеграла перекрытия атомных орбиталей.

Гамильтониан

Гамильтониан молекулы водорода H₂ записывается в виде

$\hat{H} = - \frac{\hbar^2}{2M}\Delta_{R_1} - \frac{\hbar^2}{2M}\Delta_{R_2} - \frac{\hbar^2}{2m}\Delta_{r_1} - \frac{\hbar^2}{2m}\Delta_{r_2} + \frac{e^2}{| \mathbf{R}_1 - \mathbf{R}_2|} + \frac{e^2}{| \mathbf{r}_1 - \mathbf{r}_2|}$

$- \frac{e^2}{| \mathbf{R}_1 - \mathbf{r}_1|} - \frac{e^2}{| \mathbf{R}_1 - \mathbf{r}_2|} - \frac{e^2}{| \mathbf{R}_2 - \mathbf{r}_1|} - \frac{e^2}{| \mathbf{R}_2 - \mathbf{r}_2|}$,

где M — масса протона, m — масса электрона, $\mathbf{R}_i$ — координаты ядер, $\mathbf{r}_i$ — координаты электронов.

В адиабатическом приближении задача нахождения энергетических состояний молекулы водорода разбивается на два этапа. На первом этапе рассматривается только электронная подсистема, а ядра считаются зафиксированными в точках $\mathbf{R}_1$ и $\mathbf{R}_2$.

Гамильтониан электронной подсистемы имеет вид

$\hat{H}_{el} = - \frac{\hbar^2}{2m}\Delta_{r_1} - \frac{\hbar^2}{2m}\Delta_{r_2} + \frac{e^2}{| \mathbf{r}_1 - \mathbf{r}_2|} - \frac{e^2}{| \mathbf{R}_1 - \mathbf{r}_1|} - \frac{e^2}{| \mathbf{R}_1 - \mathbf{r}_2|} - \frac{e^2}{| \mathbf{R}_2 - \mathbf{r}_1|} - \frac{e^2}{| \mathbf{R}_2 - \mathbf{r}_2|}$,

Гамильтониан молекулы водорода симметричен относительно переменных $\mathbf{r}_1$ и $\mathbf{r}_2$, то есть не изменяется при смене нумерации электронов. Кроме того, он не зависит от спиновых переменных.

Волновые функции

Схематическое изображение симметричной волновой функции молекулы водорода и образования ковалентной связи

Схематическое изображение антисимметричной волновой функции молекулы водорода

Исходя из принципа тождественности частиц и учитывая то, что электроны являются фермионами, волновая функция молекулы водорода должна быть антисимметричной относительно перестановки электронов. Независимость гамильтониана от спиновых переменных позволяет выбрать её в виде произведения двух функций, одна из которых зависит только от координат и называется координатной волновой функцией, а вторая зависит только от спинов и называется спиновой функцией.

Существуют два возможных типа спиновых функций для двух электронов: с параллельными спинами и с антипараллельными спинами. Функции с параллельными спинами симметричны относительно перестановки электронов и образуют триплетные состояния, в которых суммарный спин равен единице. Функции с антипараллельными спинами образуют синглетные состояния с нулевым суммарным спином (см. Мультиплетность, Состояние (квантовая механика)). Они антисимметричны относительно перестановки электронов.

Поскольку волновая функция электрона должна быть антисимметричной относительно перестановки, она может быть произведением либо

• симметричной координатной функции на антисимметричную спиновую функцию, либо
• антисимметричной координатной функции на симметричную спиновую функцию.

Симметричную и антисимметричную координатные функции основного состояния в нулевом приближении теории возмущений можно построить исходя из волновых функций основного состояния атома водорода $\psi_{1s} \,$ .

• $\psi_{sym} = \frac{1}{\sqrt{2}} \left( \psi_{1s1}(\mathbf{r}_1)\psi_{1s2}(\mathbf{r}_2) + \psi_{1s1}(\mathbf{r}_2)\psi_{1s2}(\mathbf{r}_1) \right)$

  $\psi_{asym} = \frac{1}{\sqrt{2}} \left( \psi_{1s1}(\mathbf{r}_1)\psi_{1s2}(\mathbf{r}_2) - \psi_{1s1}(\mathbf{r}_2)\psi_{1s2}(\mathbf{r}_1) \right)$

Первый член каждой из приведенных волновых функций отвечает локализации электрона, пронумерованного индексом 1 на атоме с индексом 1 и электрона, пронумерованного индексом 2 на атоме с индексом 2. Второй член отвечает ситуации, когда первый электрон принадлежит второму атому, а второй — первому. При смене нумерации электронов первая (симметричная) функция не меняется, вторая (антисимметричная) — меняет знак.

Энергия электронных состояний

Энергия электронных состояний находится усреднением гамильтониана на волновых функциях

$\mathcal{E} = \int \psi^* \hat{H} \psi dV_1dV_2$.

Для двух возможных волновых функций основного состояния она приобретает вид

$\mathcal{E}_{sym/asym} = 2E_{1s} + \frac{e^2}{R} + \frac{Q \pm A}{1 \pm S^2 }$,

где $R = |\mathbf{R}_1 - \mathbf{R}_2|$ — расстояние между ядрами, $E_{1s}$ — значение энергии основного состояния атома водорода, $S$ — интеграл перекрытия, который определяется выражением

$S(R) = \int \psi_{1s}(|\mathbf{r} - \mathbf{R}_1|) \psi_{1s}(|\mathbf{r} - \mathbf{R}_2|) dV$,

$Q$ — кулоновский интеграл, определяемый как

$Q(R) = e^2 \int \psi^2_{1s}(|\mathbf{r}_1 - \mathbf{R}_1|) \psi^2_{1s}(|\mathbf{r}_2 - \mathbf{R}_2|) \left[ \frac{1}{|\mathbf{r}_1 -\mathbf{r}_2|} - \frac{1}{|\mathbf{r}_1 -\mathbf{R}_2|} - \frac{1}{|\mathbf{r}_2 -\mathbf{R}_1|} \right] dV_1 dV_2$,

$A$ — обменный интерграл, равный

$A(R) = e^2 \int \psi_{1s}(|\mathbf{r}_1 - \mathbf{R}_1|) \psi_{1s}(|\mathbf{r}_2 - \mathbf{R}_2|) \psi_{1s}(|\mathbf{r}_2 - \mathbf{R}_1|) \psi_{1s}(|\mathbf{r}_1 - \mathbf{R}_2|) \left[ \frac{1}{|\mathbf{r}_1 -\mathbf{r}_2|} \right.$

$\left. - \frac{1}{|\mathbf{r}_1 -\mathbf{R}_2|} - \frac{1}{|\mathbf{r}_2 -\mathbf{R}_1|} \right] dV_1 dV_2$.

Энергия молекулы водорода состоит из удвоенной энергии атома водорода, энергии кулоновского отталкивания между ядрами и членов, определяющих взаимодействие электронов между собой и с другим ядром. В зависимости от спинового состояния в формуле для $\mathcal{E}$ следует выбирать знак «+» или знак «-». Знак «+» отвечает синглетному состоянию, знак «-» — триплетному.

Величина $Q(R)$ — это средняя энергия кулоновского взаимодействия между двумя электронами и электронов с «чужим ядром». Величина $A(R)$ определяет вклад в энергию обменного взаимодействия, которая не имеет аналога в классической механике. Именно этот член определяет существенное различие между синглетными и триплетными состояниями.

Зависимость энергии основного состояния от расстояния между ядрами

Ковалентная связь

Энергия взаимодействия атомов в молекуле водорода для синглетного(сплошная линия) и триплетного (пунктирная линия) состояний

Обменный интеграл для молеклы водорода в определённом диапазоне расстояний между атомами отрицательный. Вследствие этого, для синглетных состояний он обеспечивает дополнительное притягивание между ядрами атомов, а для триплетных — дополнительное отталкивание. Дополнительное притягивание обусловливает появление минимума электронной энергии на расстоянии приблизительно в 1,5 радиуса Бора (что соответствует приблизительно 75 пм). Этот минимум возникает только для синглетного состояния, то есть для антипараллельных спинов. Таким образом устанавливается ковалентная связь между атомами. Электроны двух атомов водорода «спариваются».

Насыщение ковалентной связи

Для параллельных спинов, то есть для триплетного состояния, дополнительный положительный вклад в энергию приводит к тому, что атому отталкиваются на любом расстоянии. Этим объясняется насыщение ковалентной связи. Молекула водорода может состоять только из двух атомов. Третий атом водорода не может образовать связь с молекулой потому, что его спин непременно будет параллельным одному из спинов электронов в составе молекулы. Параллельность спинов приводит к отталкиванию между ядрами на любом расстоянии.

Необходимо отметить, что на больших расстояниях между атомами они притягиваются в любых состояниях благодаря Ван-дер-Ваальсовому взаимодействию, которое, однако, намного слабее ковалентной связи. Ван-дер-Ваальсово взаимодействие не учитывается в приведенной теории, поскольку требует рассмотрения возбуждённых состояний высших, чем 1s-орбиталь.

Движение ядер в молекуле водорода

Гамильтониан, описывающий относительное движение двух ядер в молекуле водорода имеет вид

$\hat{H}_{nucl} = - \frac{\hbar^2}{2\mu R^2} \frac{\partial }{\partial R} \left( R^2 \frac{\partial}{\partial R} \right) + \frac{\hat{L}^2}{2\mu R^2} + \mathcal{E}_{sym}(R)$,

где $\mu = M/2$ — приведённая масса двух ядер, $\mathcal{E}_{sym}(R)$ — энергия синглетного электронного состояния, зависящая от расстояния между ядрами (определена в предыдущем разделе), $\hat{L}$ — оператор углового момента, описывающий вращение молекулы.

Вращение

Оператор углового момента коммутирует с гамильтонианом, поэтому собственные функции гамильтониана можно выбрать общими с собственными функциями оператора углового момента, то есть сферические гармоники:

$\hat{L}^2 Y_{lm} = \hbar^2 l(l+1) Y_{lm}$.

Вклад вращения в энергию молекулы равен

$E_{rot} = \frac{\hbar^2 l(l+1)}{2\mu R^2}$.

Этот вклад мал, потому, что массы ядер намного больше, чем массы электронов. Он зависит от расстояния между ядрами R. При колебаниях ядер это расстояние меняется, поэтому вращательное движение, строго говоря, невозможно отделить от колебательного. Однако при малых l и при малых энергиях колебательного движения можно считать, что R равняется определённому среднему значению.

Вращательные уровни вырождены относительно магнитного квантового числа m_l. Кратность вырождения 2l+1.

Колебания атомов

Вблизи минимума функцию $\mathcal{E}_{sym}(R) + \frac{\hbar^2 l(l+1)}{2\mu R^2}$ можно разложить в ряд Тейлора, ограничившись только квадратичным членом.

$U(R) \approx \mathcal{E}_{sym}(R) + \frac{\hbar^2 l(l+1)}{2\mu R^2} = \mathcal{E}(R_0) + \frac{\hbar^2 l(l+1)}{2\mu R^2_0} + \frac {\mu \omega^2}{2}(R- R_0^2)$,

где

$\mu\omega^2 = \left. \frac{d^2U}{dR^2} \right|_{R=R_0}$

Таким образом, задача об энергетических состояниях колебания ядер сводится к задаче квантового гармонического осциллятора. Следовательно, для низших колебательных и вращательных уровней энергия молекулы водорода равна

$E_{s,n,l} = \mathcal{E}_{sym}(R_0) + \frac{\hbar^2 l(l+1)}{2\mu R^2} + \hbar \omega \left(v + \frac{1}{2}\right)$,

где $v = 0, 1 \ldots$ — колебательное квантовое число.

Низшая энергия колебания $\frac{\hbar \omega}{2}$ отвечает нулевым колебаниям.

При больших значениях квантовых чисел n и l колебательное и вращательное движение разделить невозможно.

Таким образом, энергию молекулы при малых квантовых числах можно записать в виде трёх слагаемых — электронного, колебательного и вращательного:

$E_{s,n,l} = E_{el} + E_{v} + E_r$

Малым параметром задачи является отношение массы электрона к массе протона. Электронное слагаемое не зависит от этого отношения, колебательное пропорционально $\sqrt{m/M}$, вращательное пропорционально $m/M$. Поэтому

$E_{el} \gg E_v \gg E_r$.

Следовательно, электронный уровень расщепляется на несколько колебательных, каждый из которых, в свою очередь, расщепляется на вращательные. Поскольку энергии вращательных уровней близки, то в оптических спектрах они сливаются в полосы — образуется полосатый спектр, который наблюдается для молекулы HD. Для симметричной молекулы H₂ переходы между орбитальными и колебательными уровнями запрещены правилами отбора, поэтому молекулярный протий не поглощает в инфракрасной области.

Пара- и ортоводород

Схематическое изображене спиновых изомеров водорода

Ядра атомов водорода, входящие в состав молекулы, являются протонами — фермионами со спином 1/2. По правилам сложения спинов ядерный спин молекулы может быть 0 или 1. Молекулы водорода с суммарным ядерным спином 0 называют параводородом, а молекулы с суммарным ядерным спином 1 и тремя возможными проекциями (−1, 0, 1) называют ортоводородом. В состоянии термодинамического равновесия при комнатной температуре отношение между ортоводородом и параводородом составляет 3:1. Однако при низких температурах в состоянии термодинамического равновесия молекулы параводорода доминируют благодаря разным вращательным свойствам: параводород является симметричным ротатором, ортоводород — несимметричным ротатором. Переход ортоводорода в параводород — медленный процесс, и состояние термодинамического равновесия при низких температурах устанавливается долго.

Характеристики молекулы водорода

Межъядерные расстояния и энергии диссоциации молекул водорода, его изотопов и молекулярного иона водорода^[3]

Молекула	Межъядерное расстояние,Å	Энергия диссоциации, Эв
H2	0,7416	4,477
HD	0,7414	4,512
D2	0,7416	4,555
DT	0,7416	4,570
T2	0,7416	4,588
HT	0,7416	4,524
H2+	1,06	2,648

Константы реакции диссоциации молекулярного водорода (K_p) и степень превращения (α) в зависимости от абсолютной температуры^[4]:

T, к	2000	3000	4000	5000	6000	8000
Кр	2,62 · 10-6	2,47 · 10-2	2,52	4,09 · 10	2,62 · 102	2,70 · 103
α	8,10 · 10-4	7,83 · 10-2	0,621	0,954	0,992	0,999

См. также

• Молекулярный ион водорода
• Рекомбинация (химия)
• Боровская модель химической связи
• Кулоновская электронная корреляция
• Статистическая интерпретация волновой функции
• Теория молекулярных орбиталей
• Теория валентных связей

Примечания

1. ↑ Бор Н. Избранные научные труды (статьи 1909-1925 г.г). — М.: Наука, 1973. — Т. 1. — С. 133. — 584 с.
2. ↑ ^{1 2} Химический энциклопедический словарь / гл. редактор И. Л. Кнунянц. — М.: Советская энциклопедия, 1983. — С. 646. — 792 с.
3. ↑ Справочник химика. — 2-е изд., перераб. и доп. — М.-Л.: ГНТИ Химической литературы, 1962. — Т. 1. — С. 338. — 1072 с.
4. ↑ Справочник химика. — 2-е изд., перераб. и доп. — М.-Л.: ГНТИ Химической литературы, 1964. — Т. 3. — С. 24. — 1008 с.

Литература

• Білий М. У. Атомна фізика. — Київ: Вища школа, 1973.  (укр.)
• Федорченко А. М. Теоретична фізика. Квантова механіка, термодинаміка і статистична фізика. Т.2. — Київ: Вища школа, 1993.
• Юхновський І. Р. Основи квантової механіки. — Київ: Либідь, 2002.
• Ландау Л. Д., Лившиц Е. М. Теоретическая физика. т. ІІІ. Квантовая механика. Нерелятивистская теория. — М.: Наука, 1974.