Добавил:
researchgate.net Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Yuriy Kruglyak. Quantum Chemistry_Kiev_1963-1991

.pdf
Скачиваний:
28
Добавлен:
28.02.2018
Размер:
16.43 Mб
Скачать

где скаляры

Q = b

 

a

ψ

 

ˆ

 

ϕ

,

P = ψ

 

ˆ

 

ϕ

, операторы

ˆ

 

t

b

 

,

ˆ

 

a s

 

, и

 

 

 

 

 

 

 

 

 

 

 

A

 

 

A

 

X =

 

 

Y =

 

 

очевидно, что

ˆ

эрмитов,

а

ˆ

– неэрмитов оператор. Цепочку равенств

A

B

(6) можно, естественно, продолжить.

Дираковский формализм особенно удобен при работе с проекционными операторами. По определению, проекция функции ψ(x) на нормированную функцию ϕ(x) задается интегральным оператором

 

ˆ

 

 

(x)ψ(x)dx′= Pϕ (x, x)ψ(x)dx

(7)

 

Pϕψ(x) =ϕ(x)ϕ

 

с ядром Pϕ (x, x ) =ϕ(x)ϕ

(x ) . Вместе с

тем

левый интеграл есть

скалярное

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

произведение этих двух функций, а именно:

 

 

 

 

тогда

 

ϕ (x)ψ(x)dx′=

ϕ

 

ψ ,

(8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

ψ =

 

ϕ

ϕ

 

ψ ,

(9)

 

 

 

 

 

 

 

 

Pϕ

 

 

 

или, иначе говоря, проекционный оператор есть

 

 

 

 

 

 

 

ˆ

 

ϕ

ϕ

 

.

 

 

 

(10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pϕ =

 

 

 

 

 

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

Пусть кет-вектор ϕ соответствует вектор-столбцу c , образованному из коэффициентов разложения по ортонормированному базису { χi }, а именно:

 

ϕ = ci

 

χi .

(11)

 

 

 

i

 

 

 

Тогда бра-вектор ϕ соответствует вектор-строке c, сопряженной к c . Скалярное произведение ϕ ψ соответствует матричному скалярному произведению cd , где вектор-столбец d соответствует кет-вектору ψ .

Операторы соответствуют матрицам в том же базисе, например,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ ϕ

 

cd= C

( d)

(12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

600

есть диадное произведение вектора-строки и вектора-столбца, т. е. матрица, а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

 

ϕ dc = ( d) C

(13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

есть скаляр.

Возьмем вектор-столбец (11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ = ci

 

χi =

ci χi

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и вектор-строку

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

 

= di χi

 

=

 

 

di χi

d.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда их оба произведения таковы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

 

ϕ = di χi

 

cj χj = di χi

 

χj cj = di δij cj

= di ci = dc,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

 

 

 

 

 

 

i, j

 

 

 

 

 

 

 

 

i, j

 

 

 

 

 

 

 

 

 

i

 

 

(14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ ϕ

 

=

 

ci χi

 

d j χj

=

ci

 

χi

χj

 

d j

= ciδij d j

= cidi = cd.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

 

 

 

 

 

 

i, j

 

 

 

 

 

 

 

 

 

 

 

 

 

i, j

 

 

 

 

 

 

 

 

 

i

 

 

 

 

Для проекционного оператора (10) имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ci χi

cj χj

 

= cicj χi

χj = cc

 

 

= Pϕ .

 

 

 

 

 

 

 

 

Pϕ = ϕ ϕ =

 

 

 

 

(15)

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

j

 

 

i, j

 

 

 

= d j

 

 

 

 

 

 

 

 

 

 

Действительно,

 

проекция

 

 

 

функции

 

ψ

 

 

χj

 

 

на

функцию

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ =

ci χi c , согласно (9), равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

ψ =

 

ϕ

 

ϕ

 

ψ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

ci χi

 

 

 

cj χj

dk χk = cicjdk

 

χi

χj

 

χk =

(16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

j

 

 

 

k

 

 

 

 

i, j,k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= cicjδ jkdk

 

χi

= cicjd j

 

χi

= (ccd)i

 

χi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i, j,k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i, j

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

601

в соответствии с тем, что диадное произведение ccесть матрица оператора проекции на подпространство, натянутое на вектор c ϕ, а матричные

элементы (Pϕ )k ,l проекционного оператора Pˆϕ в ортонормированном базисе {χk } равны матричным элементам этого диадного произведения:

(Pϕ )k,l = χk

 

ˆ

 

χl

= χk

 

ϕ ϕ

 

χl =

χk

ci χi

 

χj

χl

=

 

 

 

 

 

Pϕ

 

 

 

cj

 

 

 

 

 

 

 

 

 

 

 

i

j

 

 

(17)

 

 

 

 

 

 

 

 

 

 

 

 

 

=cicjδkiδ jl = ckcl = (cc)kl. i, j

 

 

Дираковский

 

формализм

 

легко

 

 

обобщается

на

случай

неортонормированного базиса {

 

χi

 

}. В этом случае в выражениях появляется

 

матрица перекрывания S с элементами Sij = χi

 

χj . Например,

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

 

ψ = ci χi

d j χj = ci χi

 

χj d j = ci

χi

 

χj d j = ci Sijd j = cS d . (18)

 

 

 

 

 

 

i

 

j

i, j

 

i, j

 

 

 

 

i, j

 

 

 

 

Работая с неортогональным базисом, следует различать матрицу

линейного оператора с элементами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aij

 

ˆ

 

 

 

 

 

(19)

 

 

 

 

 

 

 

 

= χi | A| χj

 

 

 

 

и матрицу линейного преобразования, соответствующего этому оператору. Пусть ϕc , ψ d и ϕ = Aˆ ψ . Тогда

 

 

 

ϕ =

ci χi

 

ˆ

d j χj

,

(20)

 

 

 

 

 

 

 

= A

 

 

 

 

 

 

i

 

 

j

 

 

 

умножая которое на χl

 

, получаем

 

ˆ

 

 

 

 

 

 

 

 

 

 

χl

 

ci χi = χl

d j χj ,

(21)

 

 

 

 

 

 

| A|

 

 

 

 

 

i

 

 

 

j

 

 

или иначе

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Slici

= Aljd j ,

 

(22)

 

 

 

 

 

 

i

j

 

 

 

 

 

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

Sc = Ad ,

(23)

602

откуда

 

 

 

 

 

 

 

 

c = S1Ad .

 

 

 

 

 

 

 

 

 

(24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак,

для

случая

неортогонального базиса

матрица

линейного

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

1

A , где

преобразования, описывающая действие линейного оператора A , есть S

 

элементы матрицы A задаются формулами (19). В случае проектора

 

ˆ

Pϕ это

 

 

 

 

 

 

 

 

ˆ

есть

 

S ,

 

тогда как матрица проекционного

значит, что матрица оператора Pϕ

Scc

 

оператора равна ccS .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

До сих

 

пор

рассматривался

оператор

 

проекции

(10) на

одномерное

пространство

ˆ

=

 

ϕ ϕ

 

,

натянутое

на

вектор

 

ϕ .

В случае

n-мерного

 

 

 

Pϕ

 

 

 

подпространства

 

 

ортонормированных

 

 

векторов

{

 

ϕi } (i =1,2,3,...,n) ,

 

 

 

 

проекционный оператор на это подпространство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

n

 

ϕi

ϕi

 

.

 

 

 

 

 

(25)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

Проекционные операторы

не

только

эрмитовые,

что следует

из его

определения (25) и правил (1) – (4) обращения с векторами бра и кет, но и

идемпотентны,

ˆ2

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

P

= P . Действительно,

 

 

 

 

 

 

ˆ

 

n

 

n

 

 

 

n

 

 

 

n

n

 

ˆ

2

=

ϕi

ϕi

ϕj

ϕj

=

ϕi

ϕi

ϕj

ϕj

=

ϕi δij ϕj

=

ϕi ϕi

P

 

= P . (26)

 

 

i=1

 

 

j=1

 

 

 

i, j

=1

ˆ

 

 

 

i, j=1

 

i=1

 

 

 

 

Проекционный оператор

 

остается инвариантным при замене векторов

 

 

P

{ϕi } на любой другой набор ортонормированных векторов {ψi }, лежащих в

том же самом подпространстве. Оба набора векторов {ϕi } и {ψi }

ортонормированны, что означает, что они связаны унитарным преобразованием U , а именно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ j

= Ukj

ϕk ,

 

 

 

 

 

 

 

 

 

 

(27)

так что оператор

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

k

=1

 

 

 

 

 

ψ j },

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P, построенный на векторах {

 

 

 

 

 

 

 

 

 

ˆ

 

n

 

 

 

 

 

n

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕj ϕj

=

Ukj

ϕk

 

 

 

 

 

ϕl

=

 

 

 

 

 

 

ϕk ϕl

=

P′=

 

Ulj

 

 

UkjUlj

 

 

 

j=1

 

 

 

 

 

j

=1

k=1

 

 

 

 

 

l=1

 

 

 

 

 

 

j,k,l=1

 

 

 

 

 

(28)

 

 

n

 

l

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

=

 

 

U

(U

)

 

 

ϕ

 

 

ϕ

 

=

δ

 

ϕ

 

 

ϕ

 

 

=

ϕ

 

ϕ

 

 

 

 

 

k

l

kl

k

l

 

k

k

= P,

 

 

 

kj

 

 

 

jl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k,l=1

 

j=1

 

 

 

 

 

 

 

 

 

 

 

 

k,l=1

 

 

 

 

 

 

 

 

 

 

k=1

 

 

 

 

 

 

 

что и доказывает сформулированное выше утверждение об инвариантности проекционного оператора.

Если векторы {ϕi } образуют базис (полный ортонормированный набор функций) в гильбертовом пространстве, то любую функцию ψ можно представить в виде разложения

 

ψ = ci

 

ϕi

(29)

 

 

 

i

 

 

 

603

с коэффициентами разложения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ci = ϕi

 

ψ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(30)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

подставляя которые в разложение (29), получаем

 

 

 

 

 

 

 

 

 

 

 

 

ψ = ci

 

ϕi

= ϕi

 

ψ

 

ϕi

=

 

 

 

ϕi

 

ϕi

 

 

ψ =

 

ϕi

ϕi

 

 

 

ψ Iˆ

 

ψ ,

(31)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

i

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Iˆ =

 

ϕi

 

ϕi

 

 

,

 

 

 

 

 

 

 

 

 

 

(32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

называют «разложением единицы» по функциям полного базиса, или иначе – «условием полноты» базиса.

Пусть нам дан эрмитов оператор Aˆ . Его собственные векторы ϕi образуют полный ортонормированный набор векторов. Следовательно, на этих

векторах можно построить разложение единицы (32).

 

 

 

 

ψ и вставим

Подействуем оператором A

 

на произвольную функцию

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

после него единичный оператор Iˆ :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

ψ

ˆ ˆ

 

ψ

ˆ

 

ϕi

ϕi

 

 

 

ψ

 

ˆ

 

ϕi

ϕi

 

 

 

ψ

 

ai

 

ϕi

ϕi

 

 

 

ψ

, (33)

 

 

 

 

 

 

 

 

 

 

 

A

 

= AI

 

= A

 

 

 

 

=

A

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

где учтено, что ϕi является собственным вектором оператора Aˆ :

Aˆ ϕi = ai ϕi

с собственным значением ai . Из (33) следует, что любой эрмитов оператор можно записать в виде его спектрального разложения

 

ˆ

 

ϕi ϕi

 

(34)

 

 

 

 

A = ai

 

 

 

i

 

 

 

 

по его собственным векторам и собственным значениям.

И в заключение, приведем доказательство утверждения о том, что след

ˆ

равен размерности n

подпространства, на которое

проекционного оператора P

он проектирует.

Пусть задан оператор

n

Pˆ = ϕi ϕi ,

i=1

который определен в некотором m-мерном линейном пространстве, заданном ортонормированным базисом {ψ j }, и поэтому векторы этого оператора Pˆ

можно разложить по базису {ψ j }:

m

ϕi = cij ψ j . (35)

j=1

604

След оператора P

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

m

ˆ

 

 

m

 

n

 

 

 

ψk

= ψk

 

 

ϕi ϕi

 

ψk .

(36)

 

 

 

 

 

TrP = ψk

P

 

 

 

 

 

k=1

 

 

 

k=1

 

i=1

 

 

 

 

 

Подставим разложение (35) в (36), а именно:

ˆ

m

n

m

 

 

∑∑cij

ψ j

TrP = ψk

 

k=1

i=1

j=1

 

m

cilψl ψk l=1

n m

=∑ ∑

i=1 j,k,l=1

n m

cijcilδkjδlk = ∑∑cikcik = n , (37) i=1 k=1

где мы воспользовались ортонормированностью базиса {ψk } и тем, что каждая из функций базиса {ϕi } (i =1,2,3,...,n) нормирована на единицу:

m

cikcik =1(i =1,2,3,...,n) .

k=1

Очевидно, что такой же результат имеет место для следа матрицы P , представляющей оператор Pˆ в ортонормированном базисе {ψ j }[2].

Литература

1.P. A. M. Dirac. A new notation for quantum mechanics,

Math. Proc. Cambridge Phil. Soc., v. 35: 3, 416 – 418 (1939).

2.И. Майер. Избранные главы квантовой химии. Доказательства теорем и вывод формул.

Москва: Бином, 2006.

Приложение П-3. Перестановки и их основные свойства

Под операцией перестановки Pˆ подразумевается такая операция, при которой изменяется порядок расположения неких объектов. Пусть исходная последовательность расположения объектов задается номерами 1,2,3,..., N , а после перестановки порядок номеров стал другим: P1, P2 , P3,..., PN . Набор объектов

{ P1, P2 , P3,..., PN } тот же самый, что и объектов {1,2,3,..., N }, однако объекты расположены в другом порядке.

Каждая перестановка Pˆ определяется последовательностью P1, P2 , P3,..., PN ,

что означает следующее: выполняя перестановку Pˆ , на первое место ставится объект, который был на месте P1 , на второе место – тот, который был на месте

P2 , и т. д ., на k -ое место – тот, который исходно был k -ым. Обычно некое расположение P1, P2 , P3,..., PN объектов 1,2,3,..., N также называется перестановкой.

Число различных перестановок N объектов, включая тождественную перестановку Iˆ , не изменяющую порядок расположения объектов, есть N ! Рассмотрим две перестановки Pˆ и Qˆ , которые характеризуются, соответственно, последовательностями P1, P2 , P3,..., PN и Q1,Q2 ,Q3,...,QN . Тогда произведение перестановок PQˆ ˆ = Rˆ есть новая перестановка Rˆ , которая

605

12B2 22B2

получается при выполнении сначала перестановки Qˆ , а потом – перестановки Pˆ , при этом в общем случае PQˆ ˆ QPˆ ˆ .

Перечислим три важных свойства перестановок:

1)Каждая перестановка Pˆ имеет единственную обратную к ней перестановку Pˆ 1 такую, что Pˆ 1Pˆ = PPˆ ˆ 1 = Iˆ .

2)Если к любой исходной последовательности объектов P1, P2 , P3,..., PN

применить все N ! возможных перестановок, то получится N ! разных последовательностей. Это означает, что множество полученных последовательностей не зависит от исходной последовательности P1, P2 , P3,..., PN .

От исходной последовательности зависит только их порядок. Таким образом,

множество всех перестановок R = PQ для

данной фиксированной перестановки

ˆ

ˆ ˆ

 

Q .

P то же, что и множество всех перестановок

ˆ

 

 

ˆ

3) Каждая перестановка характеризуется своей четностью. Четность

перестановки P может принимать

два

значения: +1 или –1. Если

ˆ

 

 

 

последовательность P1, P2 , P3,..., PN

можно получить четным числом транспозиций

(перестановок двух объектов), исходя из последовательности 1,2,3,..., N , то

четность перестановки

P равна +1. Если же число транспозиций нечетно, то

 

ˆ

 

P обозначают через

четность перестановки P равна –1. Четность перестановки

 

ˆ

 

ˆ

 

 

 

 

число транспозиций p

как (1)p . Четность произведения двух перестановок

R = PQ есть произведение четностей перестановок P

и Q : (1)

 

= (1)

 

.

ˆ ˆ ˆ

ˆ

ˆ

 

r

 

p+q

 

Приложение П-4. Интерпретация электронных возбуждений в молекулах посредством чисел заполнения естественных орбиталей на примере возбуждения 12 B2 22 B2 бензильного радикала [1, 2]

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

бензильном радикале с энергией E = 2.68 эв[3]. Квантовохимический расчет в приближении Паризера – Парра – Попла (ППП) приводит к близкому значению E = 2.64 эв[4]. Согласно [4] вклад основной конфигурации в терм 22B2 равен

3.7 %, а вклады одновозбужденных конфигураций (обозначение их см. в [5]),

2 4(4 6), 2 6, 3 5, 1 4(4 7), 1 7, 1 6(2 7) равны соответственно 57.2,

27.0, 6.0, 3.5, 2.4, 0.2 %. На этом основании говорят, что электронное

606

поглощение 12B2 22B2 обязано, в основном, одночастичным переходам 2 6 и

2 4(4 6) .

Можно предложить иной подход к отнесению возбужденных электронных состояний, основанный на использовании чисел заполнения niσ естественных

спин-орбиталей Niσ . Последние являются собственными функциями задачи на собственные значения для одночастичной матрицы плотности [6, 7]

 

Pσ Nσ = n

Nσ ,

(1)

 

i

iσ

i

 

где i – набор квантовых

чисел координатной части спин-орбитали, а

σ

спиновое квантовое число.

 

 

 

 

В ОХФ числа заполнения niσ строго равны 1 или 0. В МК-приближении

учитывается корреляция электронов с противоположными спинами, в результате чего числа заполнения становятся дробными и лежат в интервале 0 < niσ <1. Изменения чисел заполнения при возбуждении можно использовать

для отнесения соответствующего возбужденного терма. Покажем это на примере поглощения 12B2 22B2 в бензильном радикале [4,5].

Волновые функции вычислялись в приближении ППП. Соответствующий гамильтониан дается формулами (259 – 261/гл. 2), параметры гамильтониана – одноэлектронные hµνcore и двухэлектронные γµν подробно описаны в [8], где

также содержатся все необходимые далее сведения по основному терму 12B2 .

Расчеты выполнялись с помощью программы CI-3 [9].

В МК-приближении волновые функции обоих термов записывались в виде разложения по базисным векторам [8], каждый из которых есть собственная функция операторов Sˆz и Sˆ2 . Вычисления проводились с семью наборами

векторов – G, I, II, III, IV, V и F. Детерминант G включает только основную конфигурацию, энергия которой минимизируется подходящим выбором вариационных коэффициентов разложения в (259 – 261/гл. 2). Набор I наряду с основной конфигурацией содержит все одновозбужденные конфигурации, набор II дополнен всеми двухвозбужденными и т.д. Набор F содержит все конфигурации, относящиеся к представлению В2 группы C2v , которые можно

построить для семиуровневой π-модели бензильного радикала. Таким образом, набор конфигураций F реализует точное решение гамильтониана (259 – 261/гл. 2) для бензильного радикала. Размерности конфигурационных наборов таковы: 1, 15, 89, 239, 361, 399, 404. Более подробные сведения о построении базисных векторов и наборов G F даны в § 2.7.3.7.

Энергии термов 12B2 и 22B2 в зависимости от точности аппроксимации их волновых функций приведены в табл. 1.

607

 

 

Таблица 1

Понижение энергии термов

12B2

и 22B2 по мере уточнения их

пробных волновых функций, эв

 

 

 

Ψ

12B2

22B2

G

–210.827095

I

–210.996808

–208.359754

II

–211.698380

–208.823692

III

–211.735728

–209.150067

IV

–211.756423

–209.164167

V

–211.756735

F

–211.756817

–209.170150

Данные табл. 1 позволяют подсчитать энергию возбуждения E(12B2 22B2 ) в зависимости от точности аппроксимации волновых функций

начального и конечного состояний (табл. 2).

Таблица 2 Энергия возбуждения E(12B2 22B2 )

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

12 B2 \ 22 B2

I

II

III

IV

F

 

 

 

 

 

 

G

2.47

2.00

1.68

1.66

1.66

I

2.64

2.17

1.85

1.83

1.83

II

3.34

2.87

2.55

2.53

2.53

III

3.38

2.91

2.59

2.57

2.57

IV

3.40

2.93

2.61

2.59

2.59

V

3.40

2.93

2.61

2.59

2.59

F

3.40

2.93

2.61

2.59

2.59

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

На функциях от G до F далее вычислялись одночастичные матрицы плотности Pα и Pβ . Из расчетов основного терма 12B2 [8,9] следует, что

608

уточнение пробной волновой функции за счет включения четырехкратновозбужденных конфигураций (набор IV) уже несущественно для уточнения матриц плотности. Это видно также из таблиц 1 и 2. Поэтому далее в табл. 3 приводим элементы матриц плотности для обоих термов 12B2 и 22B2 ,

вычисленные с набором конфигураций III. В этой же таблице приводим аналогичные данные для основного терма 12B2 , соответствующие

одноконфигурационной аппроксимации волновой функции.

Таблица 3 Элементы матриц Pα и Pβ для термов 12B2 и 22B2 , вычисленные

на пробных функциях G, III и соответственно на III. Пропущенные элементы матриц совпадают с приведенными по соображениям симметрии.

 

 

 

12B

 

2

2B

 

 

 

 

 

 

2

 

 

2

 

µν

 

 

 

 

 

 

 

 

 

G

 

 

III

III

 

Pα

 

Pβ

 

Pα

Pβ

Pα

 

Pβ

 

 

µν

 

µν

 

µν

µν

µν

 

µν

 

11

0.500000

 

0.500000

 

0.444804

0.555196

0.732155

 

0.263845

 

22

0.537746

 

0.462254

 

0.591435

0.408565

0.495481

 

0.504510

 

33

0.500000

 

0.500000

 

0.465395

0.534605

0.562032

 

0.437968

 

44

0.521705

 

0.478296

 

0.580491

0.419509

0.660935

 

0.339065

 

77

0.902803

 

0.097197

 

0.861044

0.038956

0.487883

 

0.512117

 

12

0.300082

 

0.300082

 

0.268898

0.268898

0.162176

 

0.162176

 

13

0.000000

 

0.000000

 

0.010010

–0.010010

–0.079826

 

0.079826

 

14

–0.149104

 

–0.149104

 

–0.119800

–0.119800

0.040839

 

0.040839

 

17

0.218333

 

0.218333

 

0.260109

0.260109

0.243970

 

0.243970

 

23

0.339660

 

0.339660

 

0.332154

0.332154

0.379500

 

0.379500

 

24

–0.028623

 

0.028623

 

–0.050991

0.050991

–0.045547

 

0.045547

 

25

–0.160341

 

–0.160341

 

–0.135840

–0.135840

–0.069662

 

–0.069662

 

26

0.037746

 

–0.037746

 

0.040703

–0.040703

0.008972

 

–0.008972

 

27

–0.123306

 

0.123306

 

–0.109494

0.109494

0.031721

 

–0.031721

 

34

0.329329

 

0.329329

 

0.310294

0.310294

0.200558

 

0.200558

 

35

0.000000

 

0.000000

 

0.005482

–0.005482

0.021259

 

–0.021259

 

37

–0.021554

 

–0.021554

 

–0.032766

–0.032766

–0.035570

 

–0.035570

 

47

0.093502

 

–0.093502

 

0.035161

–0.035161

–0.040237

 

0.040237

 

Необходимые

нам числа

заполнения niσ для

термов 12B2 и 22B2 в

МК-приближении получаются диагонализацией приведенных в табл. 3 матриц плотности. Окончательные результаты суммированы в табл. 4.

609

Соседние файлы в предмете Квантовая химия