Пространственная симметрия

Перестановочная симметрия волновой функции

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

если два любые электрона обмениваются состояниями (АСО), то глобальная волновая функция умножается на число –1 (перестановочная антисимметрия).

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

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

Р ₁₂ [ Ф₁ ]= Р ₁₂ [(Аa × Вb)] = Вb × Аa = Ф₃ (¹ k × Ф₁)

Следовательно, функция Ф₁ не является собственной для оператора перестановки и не обладает требуемой антисимметричностью. Аналогичный вывод получим для функции Ф₃: она переходит не сама в себя, а в функцию Ф₁. Для преодоления указанной трудности можно воспользоваться специальным математическим приемом — симметризацией. Всякий раз, когда имеются несимметричные объекты, переходящие не сами в себя, а друг в друга, из них можно построить линейную комбинацию, обладающую симметрией. Например, построим разность:

Ф₁₃ = Ф₁ – Ф₃ = Аa × Вb – Вb × Аa

Такая комбинация уже является антисимметричной:

Р ₁₂ [ Ф₁₃ ] = Р ₁₂ [(Аa × Вb – Вb × Аa) = (Вb × Аa – Аa × Вb) =

= (–1) × (Аa × Вb – Вb × Аa) = (–1) × Ф₁₃

Отсюда следует вывод: коэффициенты при первой и третьей базисных функциях должны быть одинаковы по величине и противоположны по знаку: С ₁ = – С ₃. Проделав аналогичный анализ для остальных базисных функций, мы обнаружим еще три антисимметричные комбинации:

Ф₂₄ = Ф₂ – Ф₄ = Аb × Вa – Вa × Аb (С ₂ = – С ₄)

Ф₅₇ = Ф₅ – Ф₇ = Аa × Вa – Вa × Аa (С ₅ = – С ₇)

Ф₆₈ = Ф₆ – Ф₈ = Аb × Вb – Вb × Аb (С ₆ = – С ₈)

Функции ионных РФ, построенные в виде определителя, уже сами по себе являются антисимметричными. В результате, из 10 неизвестных коэффициентов у нас осталось только 6.

Молекула водорода имеет весьма высокую пространственную симметрию. Очевидно, что электронная оболочка также должна быть симметричной, поскольку поведение электронов определяется природой ядер.

Плотность электронного облака определяется квадратом волновой функции. Ясно, что в пространственно симметричных (относительно ядерного остова) точках эти квадраты должны быть равны между собой, а значения самой волновой функции в симметричных точках могут отличаться только знаком. Таким образом, при действии любого из операторов пространственной симметрии на волновую функцию она должна переходить сама в себя; при этом волновая функция может вся целиком умножаться на +1 или –1. Следовательно, для молекулы водорода возможно всего два типа электронных состояний:

пространственно симметричные F [ Ф_s ] = (+1) × Ф_s

пространственно антисимметричные F [ Ф_as ] = (–1) × Ф_as

Для молекулы водорода в качестве операции симметрии чаще всего используют инверсию (i). Соответственно, волновые функции подразделяются на:

четные (симметричные), для которых i [ Ф_g ] = (+1) × Ф_g

нечетные (антисимметричные), для которых i [ Ф_u ] = (–1) × Ф_u

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

Проанализируем пространственную симметрию антисимметричных функций, найденных выше. Оператор инверсии меняет местами ядра, что выражается заменой друг на друга символов А и В в выражениях для волновых функций. (На спиновые состояния электронов и описывающие их спиновые функции a и× b оператор инверсии не влияет.) Например, для функции Ф₁₃ найдем:

i [ Ф₁₃ ] = i [ Аa × Вb – Вb × Аa ] =

= [ Вa × Аb – Аb × Вa ] = (– 1) Ф₂₄

Следовательно, функции Ф₁₃ и Ф₂₄ не обладают необходимой пространственной симметрией, они не являются ни четными, ни нечетными. Однако, посредством процедуры симметризации из них можно построить как четную, так и нечетную комбинации.

Рассмотрим сумму:

Ф⁺ = Ф₁₃ + Ф₂₄ = Аa × Вb – Вb × Аa + Аb × Вa – Вa × Аb =

= АВ (ab + ba) – ВА (ab + ba) = [ AB – BA ](ab + ba)

Действуя оператором инверсии на пространственную часть, нетрудно определить, что она является нечетной:

i [ AB – BA ] = [ BA – AB ] = (–1)[ AB – BA ]

Теперь построим разность:

Ф^– = Ф₁₃ – Ф₂₄ = Аa × Вb – Вb × Аa – Аb × Вa + Вa × Аb =

= АВ (ab – ba) + ВА (ab – ba) = [ AB + BA ](ab – ba)

Пространственная часть этой функции является четной:

i [ AB + BA ] = BA + AB = (+1)[ AB + BA ]

Исследуем функции Ф₅₇ и Ф₆₈.

Ф₅₇ = Аa × Вa – Вa × Аa = [ АВ – BA ] (aa)

Ф₆₈ = Аb × Вb – Вb × Аb = [ АВ – BA ] (bb)

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

Наконец, проанализируем функции Ф₉ и Ф₁₀.

Ф₉ = Аa × Аb – Аb × Аa = [ AA ] (ab – ba)

Ф₁₀ = Вa × Вb – Вb × Вa = [ BB ] (ab – ba)

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

Ф'⁺ = [ AA + BB ] (ab – ba) (четная)

Ф'^– = [ AA – BB ] (ab – ba) (нечетная)

В результате получили 6 функций, обладающих как перестановочной, так и пространственной симметрией. В перестановочном отношении все они антисимметричны, а в пространственном — распадаются на два класса: четные и нечетные. Суммируем результаты в виде таблицы.

Функ–ция	Выражение через атомные орбитали	Перестановочная симметрия	Пространственная симметрия	S	MS
(Ф^+)u	[AB – BA] (ab + ba)	антисимм.	нечетная
(Ф^–)g	[AB + BA] (ab – ba)	антисимм.	четная
Ф57	[АВ – BA] (aa)	антисимм.	нечетная		+1
Ф68	[АВ – BA] (bb)	антисимм.	нечетная		–1
(Ф'^+)g	[AA + BB] (ab – ba)	антисимм.	четная
(Ф'^–)u	[AA – BB] (ab – ba)	антисимм.	нечетная

В последних двух столбцах указаны значения квантовых чисел — спинового (S) и магнитного спинового (M_S), определяющих величину и ориентацию вектора глобального спина молекулы.

Теперь можно сконструировать три варианта глобальной волновой функции, включающих только однотипные слагаемые.

1. Четное синглетное состояние (спин S = 0)

Ф_g = D ₁(Ф^–) + D ₂(Ф'⁺) = { D ₁[ AB + BA ] + D ₂[ AA + BB ]}(ab – ba)

Эта функция содержит два коэффициента (D ₁ и D ₂), из которых фактически только один неизвестен (условие нормировки). Пространственная часть является безузловой и электронное облако не имеет разрывов.

2. Нечетное триплетное состояние спин (S = 1)

Ф_u = F ₁(Ф₅₇) + F ₂ (Ф⁺) + F ₃ (Ф₆₈) =

= [ AB – BA ] { F ₁(aa) + F ₂(ab + ba) + F ₃(bb)}

Данная функция содержит три коэффициента (F ₁, F ₂ и F ₃), из которых только два неизвестны. Эти коэффициенты влияют только на ориентацию суммарного вектора спина молекулы. Накладывая магнитное поле, можно сделать два из этих коэффициентов нулевыми. Пространственная форма электронного облака не зависит от ориентации спина.

3. Нечетное синглетное состояние (спин S = 0)

Ф'_u = Ф'^– = [ AA – BB ] (ab – ba)

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

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

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

Энергетические характеристики молекулы водорода.

Располагая выражениями для волновых функций, мы можем рассчитать энергию молекулы в каждом из трех возможных состояний по стандартной формуле: Е = ò Ф* Н Ф dv. Для этого нужно знать явный вид оператора Гамильтона, который для молекулы водорода записывается в следующем виде:

Н = Н _a + H _b + H _ab

где Н _a и Н _b — операторы Гамильтона для отдельных атомов водорода, содержащих электроны № 1 и № 2, соответственно.

Н _a = (–h²/2 m ₁)Ñ₁² – e ²/ r _1a и Н _b = (–h²/2 m ₂)Ñ₂² – e ²/ r _2b

Оператор Н _ab содержит члены, характеризующие взаимодействие между атомами (энергии притяжения электронов к "чужим" ядрам (1 ® a и 2 ® b), а также энергии межэлектронного и межъядерного отталкивания) и относится сразу к обоим атомам.

Н _ab = – e ²/ r _1b – e ²/ r _2а + е ²/ r ₁₂ + e ²/ r _ab

Рассмотрим для примера волновую функцию основного состояния

Ф = D ₁(АВ + ВА) + D ₂(AA + BB).

Вклад ионных РФ существенно меньше, чем ковалентных (D ₂ << D₁). Поэтому, в первом приближении, можно ограничиться учетом только ковалентных форм и взять приближенное выражение для волновой функции:

Ф» D (АВ + ВА).

Вначале найдем величину коэффициента D. Это можно сделать из условия нормировки.

òФ*Ф dv = D ² ò (AB + BA)*(AB + BA) dv = D ² [ ò (A*B* × AB) dv +

+ ò (A*B* × BA) dv + ò (B*A* × AB) dv + ò (B*A* × BA) dv ] = 1

В этом выражении имеется четыре интеграла. Рассмотрим первый из них. В подынтегральном выражении два сомножителя относятся к первому электрону, а два — ко второму. Такие интегралы можно представить в виде произведения более простых интегралов, каждый из которых включает только координаты только первого или только второго электрона. В результате шестимерный интеграл распадается на два трехмерных, каждый из которых представляет собой стандартное условие нормировки для атомных орбиталей. Поскольку мы принимаем атомные орбитали А и В нормированными, оба интеграла-сомножителя оказываются равными 1.

ò (A*B* × AB) dv = ò [ A*A × B*B dv _a × dv _b =

= òA*A dv _a × òB*B dv _b = 1 × 1 = 1

Проанализируем второй интеграл:

ò (A*B* × BA) dv = ò [ A* ₍₁₎ B ₍₁₎× B* ₍₂₎ A ₍₂₎ dv ₁ × dv ₂ =

= òA* ₍₁₎ B ₍₁₎ dv ₁ × òB* ₍₂₎ A ₍₂₎ dv ₂ = S × S = S ²

Трехмерные интегралы-сомножители в этом выражении называются интегралами перекрывания (S). Величина такого интеграла определяется величиной области перекрывания, в которой обе атомные орбитали (А и В) заметно отличаются от нуля.

Очевидно, что S ® 0, когда ядра а и b удаляются друг от друга, и S ® 1, когда ядра сближаются.

Третий интеграл построен точно так же, как и второй, а четвертый — так же, как первый. Следовательно, выражение для нормировки приобретает вид:

òФ*Ф dv = D ² [ 1 + S ² + S ² + 1] = 1

Отсюда получаем выражение для нормировочного множителя и волновой функции:

D = 1/(2 + 2 S ²)^1/2 и Ф» (АВ + ВА)/(2 + 2 S ²)^1/2

Подставляя это выражение в формулу для вычисления энергии, получим:

Е = ò Ф* Н Ф dv = [1/(2 + 2 S ²)] ò(АВ + ВА) * Н (АВ + ВА) dv

Раскрыв скобки в подынтегральном выражении, получим вместо одного сложного интеграла четыре более простых.

ò(АВ + ВА) * Н (АВ + ВА) dv = ò(АВ) * Н (АВ) dv + ò(АВ) * Н (ВА) dv + ò(ВА) * Н (АВ) dv + ò (ВА) * Н (ВА) dv = I + II + III + IV

Проанализируем эти четыре интеграла. Для начала можно отметить, что по структуре интегралы попарно равны друг другу(переводятся друг в друга операцией симметрии — инверсии): I = IV и II = III, и поэтому достаточно детально рассмотреть только два из них, например I и II.

Подставляя в интеграл I выражение для гамильтониана, получим сумму трех еще более простых интегралов:

ò(АВ) * Н (АВ) dv = ò(АВ) * (Н _a + H _b + H _ab)(АВ) dv =

= ò(АВ) *(Н _a)(АВ) dv + ò(АВ) *(Н _b)(АВ) dv + ò(АВ) *(Н _ab)(АВ) dv

Оператор Гамильтона в первом интеграл действует только на волновые функции типа А, описывающие состояние электрона № 1. Волновые функции типа В, описывающие состояние электрона № 2, для оператора Н _а выступают как простые числовые множители и их можно выносить за знак оператора. Это дает возможность разбить первый из трех интегралов в произведение двух трехмерных интегралов-сомножителей:

ò(АВ) *(Н _а)(АВ) dv = ò В*B dv _b • ò А* (Н _а) А dv _a

Первый сомножитель представляет собой условие нормировки для атомной орбитали В и поэтому равен 1. Второй сомножитель представляет собой энергию электрона, принадлежащего атому водорода с ядром а, т.е. Е _Нa.

Второй интеграл построен аналогично первому и равен энергии, которую имеет электрон, принадлежащий атому водорода с ядром b (Е _Нb).

Третий интеграл содержит гамильтониан, описывающий оба электрона. Поэтому этот интеграл нельзя разложить в произведение двух одноэлектронных и его нужно рассматривать как меру энергии, характеризующую интенсивность взаимодействия между двумя атомами водорода. Эта энергия называется кулоновским интегралом (J _ab). Ее можно представить и в более детальной записи:

ò(АВ) *(Н _ab)(АВ) dv = ò(АВ) * (– e ²/ r _1b – e ²/ r _2а + е ²/ r ₁₂ + e ²/ r _ab)(АВ) dv =

= ò(АВ) * (– e ²/ r _1b)(АВ) dv + ò(АВ) * (– e ²/ r _2а)(АВ) dv +

+ ò(АВ) * (е ²/ r ₁₂)(АВ) dv + ò(АВ) * (e ²/ r _ab)(АВ) dv =

= ò(В*B) dv ₂ × ò A* (– e ²/ r _1b)(А) dv ₁ + ò(A*A) dv ₁ × ò B* (– e ²/ r _1b)(B) dv ₂ +

+ ò(АВ) * (е ²/ r ₁₂)(АВ) dv + (e ²/ r _ab)ò(А*A) dv ₁ × ò(B*В) dv ₂ =

Произведение e × (A*A) dv _a представляет собой электрический заряд dq _a, локализованный в малом объеме dv _a, принадлежащем атому Н_а. Произведение dq _a × e / r _1b, стоящее под знаком первого интеграла, — это потенциальная энергия взаимодействия заряженного объема dv _a электронного облака, созданного электроном № 1 с ядром атома Н_b. Поэтому первый интеграл представляет собой полную энергию кулоновского притяжения электрона № 1, расположенного на орбитали A, к ядру b; второй — энергию кулоновского притяжения электрона, расположенного на орбитали B, к ядру a; третий — энергию кулоновского отталкивания двух электронных облаков, а четвертый — энергию кулоновского отталкивания двух ядер. Кулоновский интеграл, следовательно, отражает все кулоновские взаимодействия между составными частями двух атомов, ядра которых сближены на расстояние r _ab.

Таким образом, для интеграла I получаем следующую оценку:

I = Е _Нa+ Е _Нb+ J _ab = 2 E _H + J

Перейдем к анализу интеграла II. Подставляя в этот интеграл выражение для гамильтониана, получим, как и в предыдущем случае, сумму трех простых интегралов:

ò(АВ) * Н (ВА) dv = ò(АВ) * (Н _a + H _b + H _ab)(ВА) dv =

= ò(АВ) *(Н _a)(ВА) dv + ò(АВ) *(Н _b)(ВА) dv + ò(АВ) *(Н _ab)(ВА) dv

Первый из этих интегралов может быть разложен в произведение двух одноэлектронных:

ò(АВ) *(Н _a)(ВА) dv = ò В*А dv _b • ò А* (Н _a) В dv _a

Первый сомножитель в этом выражении представляет собой интеграл перекрывания (S). Второй сомножитель можно преобразовать, с учетом того, что гамильтониан является самосопряженным оператором, для которых выполняется правило:

ò А* Н В dv = ò В* Н * А dv

Атомная орбиталь А является собственной функцией для оператора Н _a и для нее выполняется соотношение:

Н _AА = Н _A * А = Е _На × А

Следовательно,

ò А* (Н _a) В dv _a = ò В* Н _a * А dv _a = ò В* Е _На А dv _a =

= Е _На × ò В*А dv _1a = Е _На × S

С учетом первого множителя получаем величину Е _На × S ².

Второй интеграл устроен аналогично первому и поэтому равен Е _Нb × S ².

Третий интеграл не поддается разложению на одноэлектронные сомножители и должен рассматриваться как глобальная характеристика молекулы, точнее как характеристика межатомного взаимодействия. Этот интеграл называется обменным (K). Его, как и кулоновский, можно представить в более детальной записи:

ò(АВ) *(Н _ab)(ВA) dv = ò(АВ) * (– e ²/ r _1b – e ²/ r _2а + е ²/ r ₁₂ + e ²/ r _ab)(ВA) dv =

= ò(АВ) * (– e ²/ r _1b)(ВA) dv + ò(АВ) * (– e ²/ r _2а)(ВA) dv +

+ ò(АВ) * (е ²/ r ₁₂)(ВA) dv + ò(АВ) * (e ²/ r _ab)(ВA) dv =

= ò(В*A) dv ₂ × ò A* (– e ²/ r _1b)(B) dv ₁ + ò(A*B) dv ₁ × ò B* (– e ²/ r _1b)(A) dv ₂ +

+ ò(АВ) * (е ²/ r ₁₂)(ВA) dv + (e ²/ r _ab)ò(А*B) dv ₁ × ò(B*A) dv ₂ =

Из формы интегралов в этом выражении видно, что они представляют собой особые кулоновские поправки к полной энергии молекулы. Они возникают вследствие того, что в области перекрывания плотность электронного облака молекулы не равна сумме плотностей двух исходных атомных электронных облаков. К этой сумме добавляется интерференционный вклад. Другими словами, в области перекрывания возникает дополнительный электрический заряд, величина которого в малом объеме dv равна eS (A*B) dv ₁ = eS (В*А) dv ₂. Этот дополнительный заряд взаимодействует, с одной стороны, с обеими ядрами, а с другой стороны — его части взаимодействуют между собой (уплотнение электронного облака в межъядерном пространстве приводит к повышению потенциальной кулоновской энергии).

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

В итоге получаем следующую оценку интеграла II:

II = Е _На × s ² + Е _Нb × s² + K = 2 E _H × s ² + K

Теперь, с учетом равенств I = IV и II = III, запишем полное выражение для энергии молекулы:

E_g = (I + II + III + IV)/(2 + 2 s ²) = (2 E _H + J + 2 E _H × s ² + K +

+ 2 E _H × s ² + K + 2 E _H + J)/(2 + 2 s ²) =

= [4 E _H (1 + s ²) + 2 (J + K)]/2(1 + s²) = 2 E _H + (J + K)/(1 + s²)

Таким образом, энергия молекулы водорода представлена в виде суммы двух вкладов:

1) энергии двух атомов, из которых молекула образовалась, и

2) энергетической поправки на межатомные взаимодействия.

Из формулы видно, что эта поправка состоит из двух частей — кулоновской и обменной. Их величины можно вычислить, зная явный вид атомных орбиталей А и В.

Подчеркнем, что найденное выражение для энергии носит общий характер в методе ВС. Независимо от размера молекулы, энергия всегда представляется в виде суммы энергий изолированных атомов, из которых молекула образована, и поправок (кулоновских и обменных) на межатомные взаимодействия.

Аналогичные вычисления можно выполнить и для других состояний. Например, для триплетного возбужденного состояния с волновой функцией Ф_u получим:

E_u = 2 E _H + (J – K)/(1 – s ²)

Следует подчеркнуть, что значение обменного интеграла отрицательно. Поэтому энергия триплетного состояния выше, чем основного (E_u > E_g).

Дата добавления: 2015-06-04