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

Аналитический метод решения многоэлектронного уравнения Дирака, основанный на многомерных спинорах. Монография

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

101

×eit(r1+r2 )ψ1* (1)ψ*2 (2) r1 α r2 ψ1 (1)ψ2 (2)dr1dr2.

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

 

1

 

1

 

eik (r1r2 )

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

dk ,

 

 

r1 r2

 

 

2π2

 

k2

 

 

 

 

 

 

то выражение (5.12) можно переписать в виде

Pэ1-э2 =

Γ (16)

+∞

 

eitρ

dt

 

1

 

 

α

I1 (k )I2 (k )dk ,

 

 

 

 

 

 

 

2π2

 

 

 

 

2πρ15 −∞ (it)10

 

 

 

 

k2

где

I1 (k ) = ψ1+ (r1 )eitr1+ikr1ψ1 (r1 )dr1 ,

I2 (k ) = ψ+2 (r2 )eitr2 +ikr2 ψ2 (r2 )dr2 .

(5.12)

(5.13)

(5.14)

(5.15)

(5.16)

При интегрировании по угловым переменным в (5.15) можно получить

 

 

 

 

 

 

 

1

1 e

ikr

e

ikr

 

 

Ω+1

 

(r1 )eikr1Ω 1

(r1 )dΩr

=

 

2π

 

eikr1xdx =

 

1

 

1

 

. (5.17)

 

 

 

 

 

 

 

 

 

 

ikr1

 

 

 

 

 

0

 

 

0

1

 

 

4π 1

2

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично для угловой части (5.16) можно записать

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

e

ikr

e

ikr

 

 

Ω+1

(r2 )eikr2 Ω 1

(r2 )dΩr

=

2π

eikr2xdx =

 

 

2

2

 

. (5.18)

 

 

 

 

 

ikr2

 

 

 

20+

20+

2

 

 

4π 1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дальнейшее интегрирование по координате r1 для интеграла I1 (k ) дает

 

 

 

1

 

 

eikr1 eikr1

 

 

 

 

1

 

(e

ir

k +t

 

 

ir

k t

) )r1dr1

 

I1 (k ) =

 

 

eitr1

 

 

 

 

r12dr1=

 

1(

 

)

e

1(

 

=

 

2Γ (3)

 

 

ikr

 

4ik

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

1

 

 

 

 

 

(it)

 

 

 

 

 

(it)

 

 

 

 

=

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

=

 

 

 

 

 

 

. (5.19)

 

(i(k + t))2

(i(k t))2

 

(k2 t2 )2

 

(k2 + (it )2 )2

 

4ik

 

 

 

 

 

 

 

 

 

Интегрирование по координате r2 , проведенное аналогичным образом, дает следующий вид интеграла I2 (k ):

102

I2 (k ) =

(it)

.

(5.20)

(k2 + (it)2 )2

 

 

 

Таким образом, в выражении (5.14) интеграл по переменной примет вид

 

(it)2

dk = (it)24π

 

 

1

 

 

 

dk =

 

k2

(k2 + (it)2 )4

(k2 + (it)2 )4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

7

 

 

 

 

 

 

 

 

 

 

1

 

{(it)

 

}

1

 

Γ

 

Γ

 

 

 

5π

2

 

= 4π (it)

2

 

2

2

 

2 2

 

=

 

 

. (5.21)

 

2 (it )2

4

 

 

 

 

 

Γ (4)

 

8(it)5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате для Pэ1-э2 можно получить

Pэ1-э

2

 

Γ (16) α 5π2 +∞

eitρ

Γ (16) 5α 2πρ14

75 α

 

=

 

 

 

 

 

 

dt =

 

 

 

 

 

=

 

ρ

. (5.22)

 

2πρ15

2π2

8

(it )15

2πρ15

16

Γ (15)

16

Рассмотрим второе слагаемое выражения (5.11), прямой член которого имеет вид

 

 

 

 

 

(Pэ1-э3 )

 

 

=

Γ (16)

+∞

 

eitρ

 

 

dt

×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пр

 

2πρ15 −∞ (it)10

 

 

 

 

 

×

eit(r1+r2 )ψ* (1)ψ* (2)

 

α

ψ

(1)ψ

3

(2)dr dr .

 

(5.23)

 

 

 

 

 

 

 

 

 

1

 

 

 

3

 

r1 r2

 

1

 

 

 

 

1

2

 

 

С использованием представления (5.13) можно записать

 

 

(Pэ1-э3 )

 

=

Γ (16)

+∞

 

eitρ

 

dt

1

 

α

I1 (k )I3 (k )dk ,

(5.24)

 

 

 

 

 

 

2π2

 

 

 

 

пр

 

2πρ15 −∞ (it)10

 

 

 

k2

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(it)

 

 

 

I1 (k ) = ψ1+ (r1 )eitr1+ikr1ψ1 (r1 )dr1

=

 

 

 

 

 

,

(5.25)

(k2 + (it)2 )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

103

 

 

 

 

I3 (k ) = ψ3+ (r2 )eitr2 +ikr2 ψ3 (r2 )dr2 .

(5.26)

Интеграл по угловым переменным (5.26)

имеет такой же вид,

как и для интеграла I1 (k ), и дает следующий результат:

+

 

(r2 )e

ikr

 

(r2 )dΩr =

2π 1

ikr

x

 

1 eikr2

eikr2

Ω 1

 

2

Ω 1

 

e

2

 

dx =

2

. (5.27)

 

 

 

2

0

 

 

20

2

4π 1

 

 

 

ikr2

Дальнейшее интегрирование по координате r2 для интеграла I3 (k ) дает

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

eikr2

eikr2

 

 

 

 

 

 

 

 

 

 

 

 

 

I3 (k ) =

 

 

 

 

 

L12 (itr2 )eitr2

 

 

 

 

 

 

 

r22 L12 (itr2 )dr2 =

 

 

 

 

2Γ (4)

 

 

 

ikr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

1

 

L12

(itr2 ) 2 (eir2(k +t) eir2(k t) )r2dr2.

 

 

(5.28)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12ik

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом явного вида полинома Лагерра можно записать

 

 

 

 

 

I

 

(k )

=

 

 

1

 

3 (itr

) 2

(

eir2(k +t) eir2(k t)

r dr

=

 

 

 

 

 

 

 

 

12ik

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

) 2 2

 

 

 

 

 

=

 

 

1

9 6(itr2 ) + (itr2 )2 (eir2(k +t) eir2(k t) )r2dr2 =

 

 

 

 

 

 

 

 

 

 

 

 

12ik

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k2 (it)3 (it)5 )

 

 

 

 

 

 

 

 

 

3(it )

 

 

 

 

2k2 (it) 6(it )3 4

 

 

 

=

 

+

 

 

 

 

 

 

 

 

 

 

=

 

(k2 + (it)2 )2

 

(k2 + (it)2 )3

 

 

(k2 + (it)2 )4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

((it)5 2k2 (it )3 + (it)5k4 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

(5.29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k2 + (it )2 )4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, в выражении (5.24)

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

примет вид

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

+ (it)5k

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

2

 

 

3

4

 

 

 

 

 

 

 

 

 

 

 

 

 

(it)

 

 

 

 

 

((it) 2k

 

it)

 

 

 

4π (it)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dk =

 

×

k2 (k2 + (it)2 )2

 

 

 

(k2 + (it )2 )4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

104

 

 

 

 

 

 

 

 

 

1 11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

7

 

 

 

 

5

 

Γ

 

 

Γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

Γ

 

 

 

Γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

Γ

 

Γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

×

(it )

 

 

 

2 2

 

2

 

 

(it)

 

 

 

2

 

 

+ 5

 

 

(it )

 

 

 

 

 

 

 

 

2

=

 

 

12

 

 

 

 

Γ (6)

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

Γ (6)

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

Γ (6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(it)

 

 

 

 

 

 

 

 

 

 

 

 

(it )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(it)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2π

 

 

 

Γ

 

 

 

 

 

Γ

 

 

 

7

 

 

9

 

 

 

 

 

 

7

 

 

1

 

 

 

 

 

 

 

3

 

 

 

1

 

 

 

 

 

 

π

2

 

 

1

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

2 2

 

 

 

 

2

 

 

 

 

+

5

 

 

 

=

 

 

 

 

 

 

 

 

.

(5.30)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(it )5

 

Γ (6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

it )5

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

2

 

 

 

2

 

 

2

 

 

 

 

 

 

 

В результате для прямого члена (Pэ1-э3 )

пр

 

можно получить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Pэ1-э3 )пр =

 

Γ (16) α π2 +∞ eitρ

 

 

 

 

 

 

 

 

 

Γ (16) α 2πρ14

 

 

 

15α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

dt =

 

 

 

 

 

 

 

 

4

 

 

=

 

 

 

.

(5.31)

 

2πρ15

2π2

(it)15

2πρ15

Γ (15)

 

4ρ

 

 

Для вычисления обменного члена необходимо вычислить ин-

теграл вида

 

 

 

 

 

 

 

 

 

I13 (k ) = ψ1+ (r1 )eitr1+ikr1ψ3 (r1 )dr1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.32)

Используя явный вид базисных функций, запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I13 (k ) =

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

L12 (itr1 )eitr1

eikr1 eikr1

r12dr1 =

 

 

 

 

 

 

 

 

 

 

 

 

2

3Γ (3)

 

 

 

 

 

 

ikr1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(it)

 

 

 

 

 

 

 

 

 

 

 

 

 

2 (it) 3(

it)3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4k2 (it )

 

 

 

 

 

 

 

 

=

1

 

 

 

 

3

 

 

 

 

 

 

 

 

 

+

 

k

 

=

1

 

 

 

 

 

 

 

 

 

 

.

(5.33)

3

 

(k2 + (it )2 )2

 

 

 

 

(k2 + (it)2 )3

 

 

 

 

 

3

 

 

(k2 + (it)2 )3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В интеграл по переменной k выражение I13 (k )

войдет в квадрате

 

 

(Pэ1-э3 )

 

=

 

Γ (16)

+∞

 

eitρ

 

 

dt

 

 

1

 

 

1

 

α

I13 (k )I13 (k )dk ,

(5.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обм

 

2πρ15 −∞ (it)10

 

 

2π2 3

 

 

k2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

I13 (k )I13 (k )dk =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

105

 

 

 

 

2

 

 

 

Γ

 

5

Γ

7

 

 

 

 

 

 

4k2 (it)

 

dk =

64π

 

 

2

 

2

3

π2

 

= 4π

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

.

(k2 + (it)2 )3

2(it)5

 

 

 

Γ (6)

 

8 (it)5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом для обменного члена можно получить

(Pэ1-э3 )

 

=

Γ (16)

+∞

eitρ

 

dt

 

1 α 3

π2

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(it)15

2π2 3 8

 

обм

 

 

2πρ15 −∞

 

 

 

 

 

 

 

 

=

Γ (16) 2πρ9 α

=

 

15

α .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

 

Γ 15

)

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16 ρ

 

 

 

 

 

2πρ

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Pэ24 )

 

= (Pэ1-э3 )

 

 

 

=

15α

,

 

 

 

пр

пр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4ρ

 

 

 

(Pэ24 )

 

 

= (Pэ1-э3 )

 

 

 

 

=

15

α .

 

 

 

обм

обм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16 ρ

 

 

Также можно показать, что

(Pэ1-э4 )пр = (Pэ23 )пр = 154ρα ,

(Pэ1-э4 )обм = (Pэ23 )обм = 0 .

(5.35)

(5.36)

(5.37)

(5.38)

(5.39)

(5.40)

При вычислении МЭ взаимодействия Pэ34 следует отметить, что

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

 

 

 

(Pэ34 )

 

=

Γ (16)

+∞

 

eitρ

 

dt ×

 

 

 

 

 

 

 

 

(it )10

 

 

 

 

 

 

 

 

пр

 

 

2πρ15 −∞

 

 

 

 

 

 

×

eit(r1

+r2 )ψ*

(1)ψ*

(2)

 

α

 

ψ

3

(1)

ψ

4

(2)dr dr .

(5.41)

 

 

 

 

 

 

3

 

 

4

 

 

r1

r2

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С использованием представления (5.13) можно записать

(Pэ34 )

 

=

Γ (16)

+∞

eitρ

dt

1

α

I3 (k )I4 (k )dk , (5.42)

 

 

(it)10

2π2

 

 

пр

 

2πρ15 −∞

 

 

k2

106

где

 

 

 

 

 

 

 

 

I3 (k ) = ψ3+ (r1 )eitr1+ikr1ψ3 (r1 )dr1 ,

 

 

 

 

 

 

 

(5.43)

 

 

 

 

 

 

 

 

I4 (k ) = ψ4+ (r2 )eitr2 +ikr2 ψ4 (r2 )dr2 .

 

 

 

 

 

 

 

(5.44)

Учитывая явный вид базисных функций ψ3

и ψ4 ,

можно показать,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

((it )5 2k2 (it )3 + (it)5k4 )

 

 

 

 

 

 

 

 

I3 (k ) = I4 (k ) =

 

.

 

 

 

(5.45)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k2 + (it )2 )4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

 

 

в выражении (5.42) интеграл по переменной k

примет вид

 

 

 

 

((it)5 2k2 (it)3 + (it)5k4 )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dk =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k2 (k2 + (it)2 )8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (it)2

4π

((it)8 4(it)6 k2 + 14(it)4 k4 20(it)2 k6 + 25k8 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dk =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k2 + (it)2 )8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Γ

 

1

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4π

 

 

 

Γ

2

 

 

9

 

11

 

 

13

 

 

 

 

 

 

1

 

 

 

9

 

 

 

11

 

 

 

1

 

 

3

 

9

 

 

=

 

 

 

 

2

 

 

 

 

4

 

 

 

 

+ 14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2(it)5

 

 

 

 

Γ (8)

 

 

 

 

2

2

 

 

 

 

2

 

2

 

2

2

 

 

2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

1

 

3

 

 

5

 

+ 25

1

 

 

3

 

5

 

=

 

 

 

7

 

 

 

π2

.

 

 

 

 

 

 

(5.46)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16 (it)5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

2

 

 

 

 

 

 

2 2 2

 

 

 

 

 

 

 

 

 

 

В результате

 

(Pэ34 )

 

 

 

 

Γ (16)

+∞

 

 

eitρ

 

 

 

 

 

 

π2 7α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

dt

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(it)15

 

 

2π2

 

 

64

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пр

 

 

 

 

2πρ15 −∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

Γ (16)

 

2πρ14

 

 

 

1

 

 

7α

 

=

 

105

α .

 

 

 

 

 

 

 

(5.47)

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

 

 

Γ 15

)

 

 

 

 

16

 

 

 

 

 

 

32 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πρ

 

 

 

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

107

Результаты расчета МЭ оператора электрон-электронного взаимодействия приведены в табл. 17.

Таблица 17

МЭ оператора электрон-электронного взаимодействия для базисных функций МУФ верхней компоненты

Номер

ψi (1)ψ j (2)

 

 

α

 

 

ψi (1)ψ j (2)

− ψi (1)ψ j (2)

 

 

 

α

 

ψ j (1)ψi (2)

 

 

 

 

 

 

 

 

оболочек

 

 

 

 

 

 

 

 

 

r r

 

 

r

r

 

 

 

 

i, j

 

 

1

 

2

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

1–2

 

 

 

75 α

 

 

 

 

 

0

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1–3

 

 

 

15 α

 

 

 

15 α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ρ

 

 

16 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2–3

 

 

 

15 α

 

 

 

 

 

0

 

 

 

 

 

4

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1–4

 

 

 

15 α

 

 

 

 

 

0

 

 

 

 

 

4

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2–4

 

 

 

15 α

 

 

 

15 α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ρ

 

 

16 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3–4

 

105 α

 

 

 

 

 

0

 

 

 

 

32

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Σ

 

 

 

 

 

 

 

 

675 α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

32 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогичные результаты можно получить для МЭ электронэлектронного взаимодействия для нижней компоненты. Угловые части вычисляются так же, как в главе 4 для литиеподобных ионов. Результаты вычисления МЭ оператора электрон-электронного взаимодействия для базисных функций МУФ нижней компоненты сведены в табл. 18.

108

Таблица 18

МЭ оператора электрон-электронного взаимодействия для базисных функций МУФ нижней компоненты

Номер

ξi (1)ξ j (2)

 

 

 

 

α

 

ξi (1)ξ j (2)

 

− ξi (1)ξ j (2)

 

 

 

α

 

ξ j (1)ξi (2)

 

 

 

 

 

 

 

 

 

 

оболочек

 

 

 

 

 

 

 

 

 

 

 

 

r

r

 

 

 

 

r

r

 

 

 

 

i, j

 

 

 

1

2

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1–2

 

 

 

 

75 α

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1–3

 

 

 

 

15 α

 

 

 

 

 

15 α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ρ

 

 

 

 

16 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2–3

 

 

 

 

15 α

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1–4

 

 

 

 

15 α

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2–4

 

 

 

 

15 α

 

 

 

 

 

15 α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ρ

 

 

 

 

16 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3–4

 

 

105 α

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

32 ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Σ

 

 

 

 

 

 

 

 

 

 

675 α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

32

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученные МЭ электрон-ядерного взаимодействия, опера-

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тора σi pi

и электрон-электронного взаимодействия определяют

i=1

константы в системе уравнений для амплитуд разложения верхней и нижней компонент ВФ.

Константы c1 и c2 характеризуют электрон-ядерное взаимодействие, для бериллиеподобных ионов c1 = c2 = 30α (см. табл. 16).

A

Константы a1 , a2 определяются МЭ оператора σi pi ; константы

i=1

d1 , d2 – межэлектронным взаимодействием. Для бериллиеподоб-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

109

ных

ионов

 

в

 

состоянии

 

 

0+

(см. (5.9),

 

(5.10)

и табл.

17, 18)

a = a = 15

,

d

 

= d

2

= 675

 

α .

 

 

 

 

 

 

 

 

 

 

 

 

1

2

2

 

1

 

 

 

32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С полученными константами уравнения (5.5) и (5.6) для ам-

плитуд разложения ВФ примут вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30αZ

675

 

 

 

 

 

 

 

 

 

′ −

15

 

 

1

 

 

 

 

32

 

α

 

 

 

 

 

 

M

 

 

M

 

 

E

+ 8

+

 

 

 

 

N = 0 ,

(5.48)

 

 

 

 

2ρ

4

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30αZ

675

α

 

 

 

 

 

 

 

 

N

 

15

 

 

1

 

 

32

 

 

 

 

 

 

 

 

 

+

 

 

N +

 

 

E +

 

 

 

 

 

M

= 0 .

(5.49)

 

 

 

 

 

2ρ

4

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.2. Волновая функция и энергия связи для 0+ состояния бериллиеподобных ионов с Z = 92–101

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

 

 

 

 

 

 

 

 

 

 

E = −2A

 

 

F 2

 

,

 

 

 

 

(5.50)

 

 

 

 

 

 

 

 

 

 

 

 

+ F 2

 

 

 

 

где A = 4 – число электронов,

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a a

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

a a

2

2

 

 

 

n

 

+ γ −1

1

 

 

 

 

 

 

 

 

 

A

 

n

+ γ −1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

2

 

 

 

 

P2

 

 

 

 

 

 

 

r

2

 

 

F = − A

 

 

 

 

 

 

 

 

 

+

 

+

 

 

 

 

 

 

 

 

 

, (5.51)

 

 

 

 

αP

 

 

 

 

P

 

 

 

 

 

 

α2P2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

a

a

 

 

(a + a

)2

 

P P

α2

 

 

 

 

 

 

 

 

и γ = 1+

 

1

 

2

+

 

 

1

2

 

 

1 2

 

 

 

 

 

,

a, c и d

коэффициенты,

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

A2

 

 

 

 

 

 

 

 

 

 

полученные при вычислении матричных элементов оператора

110

A

σi pi , электон-ядерного и электрон-электронного взаимодейст-

i=1

вия соответственно, nr – радиальное квантовое число.

A

Расчеты МЭ оператора σi pi , электрон-ядерного и элек-

i=1

трон-электронного взаимодействия показали, что для бериллиеподобных ионов

 

 

 

 

a1

= a2 = a,

 

 

 

(5.52)

 

 

 

 

P1 = P2 = P = (cZ d ),

 

 

 

 

 

 

 

 

 

 

 

поэтому энергия связи

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

E

= A

 

 

 

 

 

1 .

(5.53)

 

 

 

α2 (cZ d )2

 

 

 

 

 

 

 

 

 

 

 

 

 

1+

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

An +

A2a2 − α2 (cZ d )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

Для бериллиеподобных ионов в состоянии (1s2 2s2 )

получены ко-

эффициенты:

a = 15 , P = 30Z

675 .

 

 

 

 

 

 

2

 

 

32

 

 

 

 

Следует отметить, что формула (5.53) является общей для полной энергии связи иона с зарядом ядра Z и А электронами. Состояние иона определяется коэффициентами a, c и d, а также квантовым числом nr . В табл. 19 приведены значения энергии связи

бериллиеподобных ионов тяжелых элементов для различных значений квантового числа nr , в табл. 20 различными методами рас-

считана энергия связи бериллиеподобного иона урана. Из табл. 20 видно, что наиболее близкие результаты дают методы МУКФ и [20]. Авторы отдают предпочтение методу МУКФ, поскольку он позволяет учитывать все межэлектронные корреляции.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]