Файл: Татевский В.М. Классическая теория строения молекул и квантовая механика.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 265

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

решений ф ( ,

q>NI2

этой

 

системы

и

несколько

разных

наборов

соответствующих чисел

е,,

B N

/ 2

и выбрав

тот

из

наборов

Ф ] (

ф ^ , которому соответствует

минимальное значение суммы

JV/2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Ни м ы не можем

быть

уверены,

что выбранный набор ф , , . . .

... ,

ф ^ д а е т функцию Ф,

соответствующую наинизшему значению £

среди других

функций Ф,

построенных на других

наборах ф ]

( . . .

 

 

 

 

 

 

N12

 

 

 

 

 

 

 

 

фд,/ 2 ,для

которых

суммы

2 %

больше,

так как в

выражение

 

Ё (XXVIII, 25) входит

 

kr=l

 

 

 

 

 

 

 

для

последний

член,

сложным

образом

за­

висящий от

функций

ф,

 

 

Фдг/ 2 .

Поэтому

для

выделения

того

набора решений системы независимых уравнений Хартри, которо­

му соответствует наинизшее значение Е (XXVIII, 25), мы

должны

будем непосредственно вычислить значения Е для

всех

наборов

решений указанной системы N/2 независимых уравнений Хартри.

Выделив таким путем набор решений этой системы

ф , ,

ф д ^

который дает наинизшее значение Е, мы все же не сможем утвер­ ждать, что это значение Е и соответствующая ему функция Ф яв­ ляются наилучшими при выбранной конфигурации ядер для опи­

сания наинизшего по энергии электронного состояния

системы

среди всех функций вида

(XXVIII, 4), так как мы заранее

строили

Ф из N/2

разных функций

ф А г

Е, кото­

Чтобы

действительно

определить наинизшее значение

рое можно получить с функцией

Ф общего

вида

(XXVIII, 4), мы

должны будем последовательно

повторить

все

решение задачи,

выбирая для построения функции

(XXVIII, 4) разные числа.и раз­

личных функций ф£ причем

 

 

 

 

JV

 

(XXVIII, 26)

Решая все соответствующие системы уравнений Хартри и вы­ числяя для всех наборов решений таких систем значения £, мы мо­ жем выбрать наинизшее значение Е, соответствующее выбранной конфигурации ядер. Наинизшее значение Е, полученное таким пу­ тем, и соответствующая ему функция Ф и будут давать оптималь­ ное приближение для наинизшего по энергии электронного состоя­ ния системы, которое может быть получено в варианте Хартри. На ряде других принципиальных и технических особенностей метода Хартри и возникающих при его применении специфических ослож­ нениях и трудностях мы останавливатьсяне будем.

Описанным выше путем может быть, в принципе, получено зна­ чение энергии для наинизшего по энергии состояния только для одной выбранной конфигурации ядер. Чтобы определить волновую функцию и энергию основного электронного состояния системы (если система вообще окажется в каком-либо из состояний единой


частицей), необходимо исследовать Е как функцию параметров Ri, кш-б, определяющих ядерную конфигурацию для всех воз­ можных функций Ф, которые можно получить вышеописанным ме­ тодом. Основному электронному состоянию системы будет соот­ ветствовать функция Ф, для которой £ ( # ь Язк-б) будет иметь наиболее глубокий минимум, лежащий ниже диссоциационных

пределов.

 

Поскольку

функция Ф в методе Хартри имеет тот же вид, что

и функция Ф

(XXVIII, 4) в методе Хюккеля, возможные спиновые

характеристики состояний, описываемых этой функцией, такие же,

как для функции Ф (XXVIII, 4)

в методе Хюккеля,

§

3. Вариант Фока

 

Волновая функция и энергия основного электронного состоя­

ния.

Вариант Фока отличается

от варианта Хартри видом зави­

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

выбора

функций

ури(х, у, z,

а),

когда

определитель

составляется

из N/2 разных функций щ(х,у,

z),

каждой

из

которых

соответ­

ствуют

две

функции

i|3ft(х, у,

z, а),

именно

а{а)щ{х,

 

у, z) и

В ( а ) ф ь ( л : ,

у,

z).

Функции

Ф

множества,

из

которого выбирается

функция,

 

наиболее

близкая

к функции х ¥ 0

основного

электронного

состояния

системы,

зависят

от

N/2

варьируемых

функций

ф й . От­

дельная функция Ф этого множества

будет иметь вид

 

 

 

 

 

а

Ы ф,

(1) В (а,) ф, (1) а (0,) ф 2

(1)

. . . а

(0,) ф ^

(1) В (0,)

ф ^ (1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

ф :

 

а

2 ) Фі (2) В (о») Ф, (2)

 

 

 

а

2 ) Ф_Д, (2) В (а2 ) Ф_д, (2)

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

« Ы

Ф, М Р Ы

Ф, (АО

. . . . .

а ( 0 „)Ф„(Л0 В (oN)

фд, (АО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( X X V I I I , 27)

где 1/ VА/!

 

— нормирующий множитель.

 

 

 

 

 

 

 

 

 

Варьируемые функции ф й предполагаются ортогональными и

нормированными, т. е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

Ф ;

(*,

у, г) Ф і

(х,

у,

z) dx

=

bkl

{ ZI

^

\ J

\

(XXVIII, 28)

Функция Фока (XXVIII, 27) удовлетворяет принципу антисим­ метрии. Действительно, перестановка пространственных и спино­ вых координат любых двух электронов приводит к перемене мес­ тами двух соответствующих строк в определителе (XXVIII, 27). При этом функция Ф меняет знак на обратный, но не изменяет своего значения. Следовательно, функция Фока будет давать


качественно правильное описание для вероятности разных распре­ делений электронов в пространстве вокруг ядер.

Поскольку мы ограничиваемся рассмотрением только таких

определителей, в которые входит N/2 разных

функций

<ph и в кото­

рых число столбцов, содержащих спиновые функции

а (а)

и В (а),

одинаково, все функции нашего множества

описывают

состояния

с Sz =

0. Можно показать, что для функций

нашего

множества и

квадрат

вектора суммарного спина S2

равен

нулю. Заранее

может

быть неизвестно, каково должно быть

значение Sz и S2

для

основ­

ного электронного состояния конкретной рассматриваемой системы. Поэтому, выбирая из множества функций Ф наиболее близкую

к функции Wo основного электронного состояния, мы решим

только

такой вопрос: какая из функций

Ф

вида

(XXVIII, 27), описываю­

щая состояние, для которого Sz

=

0 и

S2 = 0, наиболее

близка

к функции Wo- основного состояния системы по энергии.

 

Очевидно, что, рассматривая другие множества Ф, т. е. опреде­

лители вида (XXVIII, 27), но содержащие более чем N/2

разных

функций фА , и линейные комбинации таких определителей, которые могут соответствовать значениям Sz и S2, не равным нулю, мы, воз­ можно, смогли бы найти функцию Ф, более близкую к Wo (т. е. дающую более низкое значение энергии), чем оптимальная функ­ ция из множества (XXVIII,27). Поскольку мы не будем рассматри-

ривать

функций

Ф, содержащих

большее число различных функций

(fh, чем

N/2, а

также

и линейных комбинаций функции Ф, мы

очевидно, ограничимся

решением

вопроса о том, какая из функ­

ций Ф частного вида (XXVIII, 27)

наиболее близка к точной функ­

ции Wo основного электронного

состояния по энергии *. Энергия

состояния, которое приближенно описывается функцией Ф вида (XXVIII, 27), может быть записана в виде

( X X V I I I , 29)

где

оператор

Я

имеет

тот же

вид, что и в варианте Хартри,

т. е.

это оператор электронного уравнения

(XXIV, 2).

 

 

 

 

*

Функция

(XXVIII, 27)

может

описать

основные

состояния

только

таких

молекул, для которых эти состояния синглетны, т. е. для которых S2

= 0. Это

имеет место далеко не всегда. Например, все изученные двухатомные

молекулы

АВ,

где

А и В — атомы

элементов шестой группы

О, S, Se, Те,

имеют,

по имею­

щимся

данным,

триплетные

основные

состояния

(S

=

I ) . Если

квадрат суммар­

ного

вектора

спина

S2

для основного

состояния

молекулы неизвестен,

определе­

ние

приближенной волновой функции и энергии основного электронного состоя­

ния

сильно осложняется. Тогда нужно рассматривать определители типа

(XXVIII, 27)

со всеми

возможными

распределениями

электронов

по

спинфунк-

циям

а

и р,

т. е. определители с числами

п разных

функций

фи

в

интервале

N/2

^

п

=sj N,

и их определенные линейные комбинации, соответствующие

раз­

ным

значениям

квадрата полного спина S2 .

С

оптимальными

волновыми

функ­

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


Подставляя в это выражение для Е функцию

Ф в- виде

(XXVIII, 27) и оператор Н из электронного уравнения

(XXIV, 2),

учитывая условие ортонормированности функций <ph (XXVIII, 28), можно выразить энергию Е состояния, приближенно описывающе­

гося функцией Ф (XXVIII,

27), в виде

 

 

 

2

 

2

 

 

£ = 2 2 # « +

2 Vhi-Kki)

(XXVIII, ЗО)

где

fc=I

ft,

Z = l

 

 

 

 

 

 

 

 

Hkk =

J<Pft (О И (0 <fk

W d x i •

(XXVIII, 31)

=

Г Ф І ( / ) Ф > ( О Ф І ' С Д Ф І ( / )

 

^ =

г Ф ;

( О Ф > О ) Ф ; ( Д Ф ,

dXi

{XXVIIIi33)

Вывод выражений

(XXVIII,,30) — (XXVIII, 33)

дан в приложе­

нии 5 *.

 

 

 

 

 

Чтобы функция

Ф вида (XXVIII, 27)

давала

наинизшее значе­

ние энергии, необходимо, чтобы функции фь . . . , ф„, входящие в определитель (XXVIII, 27), удовлетворяли требованию

ЬЕ = 0 (XXVIII, 34)

(что следует из основной теоремы вариационного метода) и усло­ виям ортонормированности (XXVII Г, 28).

Как показано в приложении 6, из условий (XXVIII, 34) и (XXVIII,28) получаются уравнения Фока, которым должны удов­ летворять оптимальные функции, в виде

 

 

F (Фі

Ф Nft) Ф* =

е Л

(XXVIII, 35)

 

 

 

k =

1, 2,

. . . .

п

 

 

где

F — так называемый оператор Фока.

 

 

 

 

 

* Если

для N четного функция Фока в виде определителя выбрана так, что

все

ф А (* =

1, 2, . . . , Л/) различны, т. е. в

виде

 

 

 

 

а ( а , ) ф , (1)Р ( а , ) ф 2 (1) . . . В (<т,)Ф д г (1)

 

 

 

а (aN)

ф, (N) В (<т„) Ф2 (N)

 

( X X V I I I , 27а)

 

 

... В (о^) Ф„ (ЛГ)

то

путем, изложенным в приложении 5,

легко

показать, что

выражение для Е

будет:

 

N

 

N

 

'

 

 

 

 

 

 

 

 

 

Я =

2 Я ^ +

S

Vkt-Kki)

(XXVIII,30а)

 

 

 

fc=i

ft,

г=і

 

 

 

т. е. будет отличаться от (XXVIII, 30) только отсутствием коэффициента 2 перед Hkk и Jut и пределами сумм. Все дальнейшие результаты, полученные в гл. XXXI на основе выражения (XXVIII, 30), могут быть совершенно так же по­ лучены на основе формулы (XXVIII,30а), если соответствующие состояния рас­ сматриваемых молекул могут быть описаны функцией (XXVIII,27а).


В развернутом виде уравнения Фока (XXVIII, 35) будут

„,,,

 

,

r

f

Фг 0") Фг (/) .

ф^ ( * )

 

 

 

ffW

+

 

l

2

J

^

d r l

 

 

 

 

 

 

1=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

Ф/

(/) ФА О')

 

 

 

 

 

 

 

 

 

- 2 !

 

dx ^

fj'W = е Л

W

( X X V I I I , 36)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1-і

 

 

 

 

 

 

 

 

 

 

 

 

fc= 1, 2,

. . . . n

 

 

 

Уравнения

Фока

представляют

собой

систему

Л//2

уравнений

для Л72 искомых оптимальных функций <рь. Особенностью системы уравнений Фока является то, что оператор F зависит от всех иско­ мых функций щ. В отличие от уравнений Хартри во всех N/2 урав­ нениях Фока, определяющих N/2 искомых функций, оператор F один и тот же. В остальном система уравнений Фока аналогична системе уравнений Хартри. Так же, как в случае уравнений Хартри, система уравнений Фока может иметь несколько наборов решений, содержащих каждый N/2 функций щ и N/2 собственных значений ЄЙ, отвечающих этим функциям. Чтобы выбрать из этих возможных наборов решений уравнений Фока один набор, позволяющий опре­ делить приближенную функцию (XXVIII, 27), соответствующую наинизшей энергии, необходимо вычислить приближенное значение

энергии Е с каждым

из этих наборов функций щ и выбрать такой

из этих наборов, для

которого значение Е оказывается наименьшим.

Для выбора из полученных описанным путем функций Ф той, которая наиболее подходит для описания основного электронного состояния системы (если система в каких-либо из состояний, соот­ ветствующих полученным функциям Ф, является единой частицей),

нужно исследовать на

минимум величины E(Rr,

ЯзК-в), полу­

ченные с разными функциями Ф. За функцию Ф, наиболее

близкую

к W0, среди множества

функций (XXVIII, 27) должна

быть

принята

та, для которой E(RU

i?3 K -e) имеет наиболее глубокий мини­

мум, лежащий ниже диссоциационных пределов.

Распределение отрицательного заряда в пространстве вокруг ядер. Вычислим вероятность dWql) нахождения электрона с но­ мером і в элементе объема dx = dxdydz с координатами х, у, z при распределении остальных электронов в пространстве в соответ­ ствии с вероятностями отдельных конфигураций, используя функ­ цию Фока (XXVIII, 27). Из общего выражения для этой вероят­ ности (XXIV, 38) получим

dwf = dx J"

da

J"

Ф * Ф dx{...

dxt_x dxi+1

...dxN

( X X V I I I , 37)

a

1, 2

i-l, i+l

Ы

 

 

 

Интегрируя по всем спиновым координатам и по всем простран­ ственным координатам электронов, кроме пространственных коор-

4 Q 0