Yuriy Kruglyak. Quantum Chemistry_Kiev_1963-1991
.pdfгде скаляры |
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 . Скалярное произведение ϕ ψ соответствует матричному скалярному произведению c†d , где вектор-столбец d соответствует кет-вектору ψ .
Операторы соответствуют матрицам в том же базисе, например,
|
|
|
• |
|
|
|
|
|
|
|
• |
|
|
|
|
|
|
|
|
|
|
|
|
|
• |
|
|
|
ϕ ϕ |
|
↔ cd† = C |
( d† ) |
(12) |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
• |
|
|
|
|
|
|
|
|
|
|
|
|
|
• |
|
|
|
|
|
• |
|
|
600
есть диадное произведение вектора-строки и вектора-столбца, т. е. матрица, а
|
|
• |
|
|
|
|
|
• |
|
|
|
|
|
|
|
|
|
• |
|
ψ |
|
ϕ ↔ d†c = ( 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 = d†c, |
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
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 |
= ∑(cc†d)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 = c†S 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 = S−1Ad . |
|
|
|
|
|
|
|
|
|
(24) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Итак, |
для |
случая |
неортогонального базиса |
матрица |
линейного |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ |
|
−1 |
A , где |
преобразования, описывающая действие линейного оператора A , есть S |
|
||||||||||||||||||||||
элементы матрицы A задаются формулами (19). В случае проектора |
|
ˆ |
|||||||||||||||||||||
Pϕ это |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
ˆ |
есть |
|
† |
S , |
|
тогда как матрица проекционного |
|||||||||
значит, что матрица оператора Pϕ |
Scc |
|
|||||||||||||||||||||
оператора равна cc†S . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
До сих |
|
пор |
рассматривался |
оператор |
|
проекции |
(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
получается при выполнении сначала перестановки 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