Перейти к содержимому
Метод Хюккеля. Аллильный радикал и бутадиен

Метод Хюккеля. Аллильный радикал и бутадиен

Второй простой метод квантовых расчетов — метод Хюккеля для сопряженных π-систем. Это метод МО ЛКАО, в котором рассматриваются только pzp_z-орбитали атомов углерода сопряженной цепи, а матричные элементы не вычисляются, а задаются:

  • Hνν=αH_{\nu\nu} = \alphaкулоновский интеграл, одинаковый для всех атомов углерода;
  • Hνμ=βH_{\nu\mu} = \betaрезонансный интеграл, для соседних атомов (отрицательная величина); для несоседних Hνμ=0H_{\nu\mu} = 0;
  • Sνν=1S_{\nu\nu} = 1, Sνμ=0S_{\nu\mu} = 0 — перекрыванием разных орбиталей пренебрегаем.

Аллильный радикал

Три атома углерода, три pzp_z-орбитали χ1,χ2,χ3\chi_1, \chi_2, \chi_3, три π-электрона:

Аллильный радикал CH₂=CH–CH₂·: три p_z-орбитали на атомах 1, 2, 3

ψМО=c1χ1+c2χ2+c3χ3 \psi_{МО} = c_1\chi_1 + c_2\chi_2 + c_3\chi_3

Вариационная система с тремя уравнениями; подставляем хюккелевские значения матричных элементов (H11=αH_{11} = \alpha, H12=βH_{12} = \beta, H13=0H_{13} = 0, S11=1S_{11} = 1, S12=0S_{12} = 0, …):

{c1(H11εS11)+c2(H12εS12)+c3(H13εS13)=0c1(H21εS21)+c2(H22εS22)+c3(H23εS23)=0c1(H31εS31)+c2(H32εS32)+c3(H33εS33)=0 \begin{cases} c_1\left( H_{11} - \varepsilon S_{11} \right) + c_2\left( H_{12} - \varepsilon S_{12} \right) + c_3\left( H_{13} - \varepsilon S_{13} \right) = 0 \\ c_1\left( H_{21} - \varepsilon S_{21} \right) + c_2\left( H_{22} - \varepsilon S_{22} \right) + c_3\left( H_{23} - \varepsilon S_{23} \right) = 0 \\ c_1\left( H_{31} - \varepsilon S_{31} \right) + c_2\left( H_{32} - \varepsilon S_{32} \right) + c_3\left( H_{33} - \varepsilon S_{33} \right) = 0 \end{cases} {c1(αε)+c2β=0c1β+c2(αε)+c3β=0c2β+c3(αε)=0 \begin{cases} c_1\left( \alpha - \varepsilon \right) + c_2\beta = 0 \\ c_1\beta + c_2\left( \alpha - \varepsilon \right) + c_3\beta = 0 \\ c_2\beta + c_3\left( \alpha - \varepsilon \right) = 0 \end{cases}

Делим каждое уравнение на β\beta и вводим безразмерную переменную

x=αεβ x = \frac{\alpha - \varepsilon}{\beta} {c1x+c2=0c1+c2x+c3=0c2+c3x=0 \begin{cases} c_1 x + c_2 = 0 \\ c_1 + c_2 x + c_3 = 0 \\ c_2 + c_3 x = 0 \end{cases}

Секулярное уравнение

x101x101x=0x(x21)x=0x(x22)=0 \begin{vmatrix} x & 1 & 0 \\ 1 & x & 1 \\ 0 & 1 & x \end{vmatrix} = 0 \qquad \Rightarrow \qquad x\left( x^2 - 1 \right) - x = 0 \qquad \Rightarrow \qquad x\left( x^2 - 2 \right) = 0

Это секулярное уравнение. Корни:

x=0,x=±2 x = 0, \qquad x = \pm\sqrt{2}

Энергии — из ε=αβx\varepsilon = \alpha - \beta x:

ε=α+2β,ε=α,ε=α2β \varepsilon = \alpha + \sqrt{2}\beta, \qquad \varepsilon = \alpha, \qquad \varepsilon = \alpha - \sqrt{2}\beta

Так как β<0\beta < 0, нижняя орбиталь — α+2β\alpha + \sqrt{2}\beta (связывающая), верхняя — α2β\alpha - \sqrt{2}\beta (разрыхляющая), а орбиталь с энергией α\alphaнесвязывающая: ее энергия равна энергии исходной pzp_z-орбитали. Три электрона: два на связывающей, один на несвязывающей.

Диаграмма МО аллильного радикала: α + √2β (два электрона), α — несвязывающая (один электрон), α − √2β

Коэффициенты

x=0x = 0. Система дает

{c2=0c1+c3=0c2=0{c2=0c1=c3 \begin{cases} c_2 = 0 \\ c_1 + c_3 = 0 \\ c_2 = 0 \end{cases} \qquad \Rightarrow \qquad \begin{cases} c_2 = 0 \\ c_1 = -c_3 \end{cases}

Условие нормировки:

c12+c22+c32=1c12+c32=12c12=1,c1=12 c_1^2 + c_2^2 + c_3^2 = 1 \qquad \Rightarrow \qquad c_1^2 + c_3^2 = 1 \qquad \Rightarrow \qquad 2c_1^2 = 1, \quad c_1 = \frac{1}{\sqrt{2}} ψ0=12χ112χ3 \psi_0 = \frac{1}{\sqrt{2}}\chi_1 - \frac{1}{\sqrt{2}}\chi_3

x=2x = -\sqrt{2}.

{2c1+c2=0c12c2+c3=0c22c3=0{c2=2c1c2=2c3{c1=c3c2=2c1 \begin{cases} -\sqrt{2}c_1 + c_2 = 0 \\ c_1 - \sqrt{2}c_2 + c_3 = 0 \\ c_2 - \sqrt{2}c_3 = 0 \end{cases} \qquad \Rightarrow \qquad \begin{cases} c_2 = \sqrt{2}c_1 \\ c_2 = \sqrt{2}c_3 \end{cases} \qquad \Rightarrow \qquad \begin{cases} c_1 = c_3 \\ c_2 = \sqrt{2}c_1 \end{cases} c12+c22+c32=1c12+2c12+c12=14c12=1,c1=12 c_1^2 + c_2^2 + c_3^2 = 1 \qquad \Rightarrow \qquad c_1^2 + 2c_1^2 + c_1^2 = 1 \qquad \Rightarrow \qquad 4c_1^2 = 1, \quad c_1 = \frac{1}{2}

x=2x = \sqrt{2} — так же, только c2=2c1c_2 = -\sqrt{2}c_1. Для удобства 222=22\dfrac{\sqrt{2}}{\sqrt{2}\cdot\sqrt{2}} = \dfrac{\sqrt{2}}{2}:

ψ2=12χ1+22χ2+12χ3 \psi_{-\sqrt{2}} = \frac{1}{2}\chi_1 + \frac{\sqrt{2}}{2}\chi_2 + \frac{1}{2}\chi_3 ψ2=12χ122χ2+12χ3 \psi_{\sqrt{2}} = \frac{1}{2}\chi_1 - \frac{\sqrt{2}}{2}\chi_2 + \frac{1}{2}\chi_3

Вид орбиталей

Орбитали аллильного радикала вдоль цепи: связывающая (плотность на всех трех атомах, максимум в центре), несвязывающая ψ₀ (узел на среднем атоме), разрыхляющая (два узла)

На связывающей орбитали электронная плотность распределена по всем трем атомам с максимумом на среднем. Несвязывающая ψ0\psi_0 имеет узел на среднем атоме: ее электрон сидит только на крайних атомах. Разрыхляющая — с двумя узлами.

Не видим двойной связи. π-Связь охватывает всю структуру в целом. Так же мы не видим отдельного электрона. Мы видим его одновременно распределенным между 1 и 3 атомом с равной вероятностью.

Электронная плотность и порядок связи

Матрица связей, известная из метода МО,

Pλσ=2jcjλcjσ P_{\lambda\sigma} = 2\sum_j c_{j\lambda}c_{j\sigma}

очень информативна — несет информацию об электронной плотности на атомах и о порядке связи. Ее диагональный элемент — электронная плотность qνq_\nu на атоме ν\nu, сумма по занятым орбиталям с числом электронов nin_i на каждой:

qν=iзанcνi2ni q_\nu = \sum_i^{зан} c_{\nu i}^2\,n_i

Недиагональный элемент — порядок π-связи между атомами ν\nu и μ\mu:

Pνμ=iзанciνciμni P_{\nu\mu} = \sum_i^{зан} c_{i\nu}c_{i\mu}\,n_i

Для аллильного радикала (два электрона на ψ2\psi_{-\sqrt{2}}, один на ψ0\psi_0):

q1=(12)22+(12)21=1,q2=(22)22=1,q3=1 q_1 = \left( \frac{1}{2} \right)^2 \cdot 2 + \left( \frac{1}{\sqrt{2}} \right)^2 \cdot 1 = 1, \qquad q_2 = \left( \frac{\sqrt{2}}{2} \right)^2 \cdot 2 = 1, \qquad q_3 = 1

Электронная плотность в единицах заряда электрона. Избыточный заряд на атоме — разность между числом электронов, которое атом отдал в π-систему, и плотностью:

ην=nνqν \eta_\nu = n_\nu - q_\nu

Если на атоме выше, чем единица, то избыточный заряд — отрицательный; если ниже — положительный. В радикале все q=1q = 1, зарядов нет.

Порядок π-связи между 1 и 2 атомом: произведение коэффициентов первого и второго на каждой занятой орбитали, умноженное на число электронов:

P12=21222+1120=22,P23=22212=22 P_{12} = 2 \cdot \frac{1}{2} \cdot \frac{\sqrt{2}}{2} + 1 \cdot \frac{1}{\sqrt{2}} \cdot 0 = \frac{\sqrt{2}}{2}, \qquad P_{23} = 2 \cdot \frac{\sqrt{2}}{2} \cdot \frac{1}{2} = \frac{\sqrt{2}}{2}

Обе связи одинаковы, порядок каждой 0,71\approx 0{,}71. Аллильный радикал — полностью делокализованная электронная структура:

Аллильный радикал и аллил-анион с пунктирными делокализованными связями

Аллил-анион

Из рассмотрения аллильного радикала методом Хюккеля легко получить анион: орбитали те же, добавляется один электрон на несвязывающую орбиталь.

Диаграмма МО аллил-аниона: два электрона на α + √2β и два на несвязывающей α

q1=(12)22+(12)22=12+1=1,5,q2=(22)22=1,q3=1,5 q_1 = \left( \frac{1}{2} \right)^2 \cdot 2 + \left( \frac{1}{\sqrt{2}} \right)^2 \cdot 2 = \frac{1}{2} + 1 = 1{,}5, \qquad q_2 = \left( \frac{\sqrt{2}}{2} \right)^2 \cdot 2 = 1, \qquad q_3 = 1{,}5 η1=η3=11,5=0,5 \eta_1 = \eta_3 = 1 - 1{,}5 = -0{,}5

Распределение избыточного заряда в аллил-анионе: −½ на концах, 0 в центре

Избыточный отрицательный заряд — на концах.

Итак, метод Хюккеля решает три типа задач:

  1. Дана молекула, надо найти энергии орбиталей.
  2. Уже даны энергии орбиталей, и надо найти волновые функции.
  3. Даны волновые функции, и надо найти распределение электронов.

Бутадиен

Четыре атома, четыре pzp_z-орбитали, четыре π-электрона:

Бутадиен CH₂=CH–CH=CH₂: четыре p_z-орбитали

ψ=c1χ1+c2χ2+c3χ3+c4χ4,x=αεβ \psi = c_1\chi_1 + c_2\chi_2 + c_3\chi_3 + c_4\chi_4, \qquad x = \frac{\alpha - \varepsilon}{\beta} x1001x1001x1001x=0{c1x+c2=0c1+c2x+c3=0c2+c3x+c4=0c3+c4x=0 \begin{vmatrix} x & 1 & 0 & 0 \\ 1 & x & 1 & 0 \\ 0 & 1 & x & 1 \\ 0 & 0 & 1 & x \end{vmatrix} = 0 \qquad\qquad \begin{cases} c_1 x + c_2 = 0 \\ c_1 + c_2 x + c_3 = 0 \\ c_2 + c_3 x + c_4 = 0 \\ c_3 + c_4 x = 0 \end{cases}

Учет симметрии

Определитель четвертого порядка раскрывать долго. Если провести учет свойств симметрии, задача упрощается. Элементы симметрии, которые имеются в молекулярной системе, одновременно действуют и на молекулярные орбитали. Если атомы являются симметричными относительно некоторой операции симметрии, то квадраты их коэффициентов в молекулярных орбиталях будут равны.

Бутадиен и ось симметрии C₂, переводящая атом 1 в 4 и 2 в 3

Ось C2C_2 переводит атом 1 в 4, а 2 в 3:

c12=c42,c22=c32 c_1^2 = c_4^2, \qquad c_2^2 = c_3^2

Значит, каждая орбиталь либо симметрична (SS), либо антисимметрична (AA) относительно этой оси — дополнительные условия связи коэффициентов.

Симметричные орбитали: c1=c4c_1 = c_4, c2=c3c_2 = c_3. Подставляем в систему — четыре уравнения сводятся к двум:

{c1x+c2=0c1+c2x+c2=0c2+c2x+c1=0c2+c1x=0x11x+1=0 \begin{cases} c_1 x + c_2 = 0 \\ c_1 + c_2 x + c_2 = 0 \\ c_2 + c_2 x + c_1 = 0 \\ c_2 + c_1 x = 0 \end{cases} \qquad \Rightarrow \qquad \begin{vmatrix} x & 1 \\ 1 & x + 1 \end{vmatrix} = 0 x2+x1=0x=1,62,x=0,62 x^2 + x - 1 = 0 \qquad \Rightarrow \qquad x = -1{,}62, \quad x = 0{,}62

Антисимметричные орбитали: c1=c4c_1 = -c_4, c2=c3c_2 = -c_3:

{c1x+c2=0c1+c2(x1)=0x11x1=0 \begin{cases} c_1 x + c_2 = 0 \\ c_1 + c_2\left( x - 1 \right) = 0 \end{cases} \qquad \Rightarrow \qquad \begin{vmatrix} x & 1 \\ 1 & x - 1 \end{vmatrix} = 0 x2x1=0x=0,62,x=1,62 x^2 - x - 1 = 0 \qquad \Rightarrow \qquad x = -0{,}62, \quad x = 1{,}62

Четыре корня — четыре орбитали. Четыре электрона занимают две нижние, связывающие:

Диаграмма МО бутадиена: α + 1,62β и α + 0,62β (связывающие, по два электрона), α − 0,62β и α − 1,62β (разрыхляющие)

Энергия делокализации

Энергия связывания сопряженной системы — энергия π-электронов на молекулярных орбиталях минус их энергия на изолированных pzp_z-орбиталях:

Eсв=2(α+1,62β)+2(α+0,62β)4α=4,48β E_{св} = 2\left( \alpha + 1{,}62\beta \right) + 2\left( \alpha + 0{,}62\beta \right) - 4\alpha = 4{,}48\beta

Это отрицательная величина (β<0\beta < 0). Принято сравнивать с молекулой этилена: в молекуле этилена два электрона на орбитали α+β\alpha + \beta, и энергия связывания EF=2βE_F = 2\beta. В бутадиене две «двойные связи» дали бы 2EF=4β2E_F = 4\beta. Разность — есть энергия резонанса, или энергия делокализации:

Eрез=Eсв2EF=0,48β E_{рез} = E_{св} - 2E_F = 0{,}48\beta

Общая формула для линейной цепи

Секулярные определители аллила и бутадиена устроены одинаково: xx на диагонали, единицы рядом с ней, нули везде еще:

x101x11x101x=0 \begin{vmatrix} x & 1 & & & 0 \\ 1 & x & 1 & & \\ & 1 & x & \ddots & \\ & & \ddots & \ddots & 1 \\ 0 & & & 1 & x \end{vmatrix} = 0

Эта общность наталкивает на мысль, что существует какая-то общая формула. Для линейной цепи из NN атомов корни выражаются явно:

xk=2cosπkN+1,k=1,2,,N x_k = -2\cos\frac{\pi k}{N + 1}, \qquad k = 1, 2, \dots, N

, где NN — порядок определителя (число атомов в цепи сопряжения), kk — номер уровня. Проверка: при N=3N = 3 получаем 2,0,2-\sqrt{2}, 0, \sqrt{2}, при N=4N = 4±1,62\pm 1{,}62 и ±0,62\pm 0{,}62.