Материал: Нейтринные процессы во внешнем магнитном поле в технике матрицы плотности (90

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Дирака заряженного фермиона во внешнем постоянном однородном магнитном поле:

ψ(+)

 

 

(x) =

e−i(Ent−p2y−p3z)

U(+)

 

(η),

 

,p3

,s

 

 

 

 

 

n,p2

 

 

 

2EnLyLz

 

n,p2,p3,s

 

 

 

(+)

,p3

,s=ϱ(η) =

s

 

n(

)

 

−s

 

n−1( )

 

(3.25)

Un,p2

χ

V

χ

,

 

 

 

 

 

W

η

 

 

 

 

η

 

Un,p(+)2

,p3

,s=−ϱ(η) = V−s χn(η) + Ws χn−1(η),

 

где En = p23 + m2 + 2eBn – энергия частицы, индекс n=ν +(1−ϱs)/2 (ν = 0, 1, 2 . . . ) нумерует уровни Ландау частицы, xµ= (t, x, y, z), а Ly, Lz – нормировочные длины вдоль осей Oy и Oz. Здесь также введены функции одномерного гармонического осциллятора:

 

(eB)1/4e−η2/2

 

 

2 dk 2

 

χk(η) =

 

 

 

 

Hk(η),

Hk(η) = (1)keη

 

 

e−η ,

(3.26)

2kk!

 

 

k

π

 

где η = eB(x − ϱp2/eB), а Hk(η) – полиномы Эрмита, и биспиноры в следующем виде:

 

 

 

 

E + m

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

W+ =

1

 

n

0

,

W

 

=

1

 

En + m

 

, (3.27)

 

 

 

 

 

 

 

 

 

 

En + m

 

p3

 

En + m

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

p3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

V+ =

1

 

0

,

V

 

 

=

1

 

0

. (3.28)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

En + m

i2eBn

 

 

 

 

En + m

 

0

 

 

 

 

 

 

 

 

2eBn

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Каждый из определенных биспиноров непосредственно не имеет стандартной нормировки, однако такой нормировкой обладает их комбинация:

WsWs + VsVs = 2m,

(3.29)

где Ws, Vs – сопряженные по Дираку биспиноры. Отметим, что биспиноры Ws, Vs не ортогональны друг к другу.

Решение уравнения Дирака с отрицательной энергией может быть получено из выражения (3.25) формальной заменой:

En → −En, p2 → −p2, p3 → −p3,

16

что в биспиноре Un,p(+)2,p3,s(η) эквивалентно заменам

m → −m, Vs → −Vs.

Таким образом, отрицательно частотное решение может быть пред-

ставлено в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ()

 

 

 

(x) =

ei(Ent−p2y−p3z)

 

U()

 

 

 

η),

 

 

 

 

 

 

,p3

,s

 

 

 

 

 

 

 

 

 

 

 

 

 

n,p2

 

 

 

 

 

 

 

2EnLyLz

 

 

n,p2,p3,s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Un,p2

,p3

,s=ϱ

η) =s

χ

n(˜) +

V

−s

χ

n−1(˜)

,

(3.30)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

η

 

 

 

 

 

 

 

η

 

 

f

 

 

Un,p()2,p3

,s=−ϱη) = −V−s χnη) + Ws χn−1η),

 

 

 

 

(

 

 

 

 

 

)

e

 

(

 

e

 

 

 

 

)

, ˜f

 

 

 

где

s

= W

s

m

→ −

m

,

V

s =

V

s

 

→ −

m

 

 

 

 

2

, и

 

W

 

 

 

 

 

 

 

 

m

 

 

 

η =

 

 

eB(x + ϱp

/eB)

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

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

поляризации спина в виде:

[Σ × P ]3

 

 

µˆ3 = mΣ3 + ρ2

,

(3.31)

 

 

 

 

 

где Σ — оператор дираковского спина, P

= p − ϱeA, p — оператор

кинематического импульса,

()

ρ2 =

0

−iI .

(3.32)

 

iI

0

 

Положительно частотное решение, являющееся собственной функцией оператора µˆ3 может быть представлено в виде (3.25), где под биспинорами Ws и Vs следует понимать следующие [7]:

 

 

 

 

 

 

 

 

 

0

,

W+ = w

E

 

+ p˜

 

 

 

 

 

 

 

W

= w

En + p˜

 

n

0

 

,

(3.33)

 

 

 

 

 

 

 

0

 

 

 

 

 

 

p3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p3

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= v

0

 

 

 

 

 

p

 

 

,

 

 

 

 

 

 

V

−p3

,

V+ = v

 

03

 

 

 

(3.34)

 

 

 

 

 

 

 

0

 

 

 

En + p˜

 

 

 

 

 

 

 

 

 

 

 

 

En + p˜

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p˜ + m

 

 

 

 

 

p˜ − m

 

 

 

 

 

 

 

 

 

w =

 

 

 

, v =

 

 

 

,

p˜

 

=

 

 

2eBn + m2.

 

p (En + p˜ )

p (En + p˜ )

 

 

 

 

 

 

17

4.Матрица плотности заряженной частицы в постоянном однородном магнитном поле

Вданном разделе мы получим импульсное представление матрицы плотности частицы заряда ϱe, массы m в постоянном однородном магнитном поле. Покажем, что, просуммированное по поляризациям

и уровням Ландау, оно в пределе слабого поля B → 0 переходит в

b

известное выражение P ±mI, соответствующие положительно и отрицательно частотному решению уравнения Дирака в вакууме.

Чтобы получить матрицу плотности заряженной частицы в однородном магнитном поле, рассмотрим следующий интеграл:

 

 

 

 

 

 

 

 

 

(+)

 

(x) dp2 =

 

 

(4.1)

(+)

 

(+)

 

 

 

 

 

 

 

 

In,p3,s

(x, x) = ψn,p2,p3,s(x) ψn,p2,p3,s

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

[

 

 

 

]

 

 

 

e−i[En(t−t)−p3(z−z)]

 

 

 

 

 

 

 

 

(+)

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

−∞

e

2

 

 

Un,p2,p3,s(η) Un,p2,p3,s(η)

 

 

2EnLyLz

 

 

 

 

 

 

)

dp2,

 

 

 

 

 

 

 

 

ip (y

y

(+)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где η = eB(x−ϱp2/eB) и η=

 

eB(x−ϱp2/eB). Перейдём от интегри-

рования по переменной p2 к интегрированию по η. Выделив трансляционно неинвариантную фазу Φ(x, x)= eB(x + x)(y −y)/2, получим:

I(+)

(x, x) =

eB

 

eΦ(x,x) e−i[En(t−t)−p3(z−z)]F (ξ1, ξ2),

(4.2)

 

 

 

n,p3,s

 

 

 

 

L

 

 

 

 

 

 

 

e

 

e

 

 

 

2EnLyz

[

 

 

(+)

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

]

 

F (ξ1, ξ2) = eiϱ ξ1ξ2/2

 

(+)

(η) Un,p2,p3,s

(η − ξ1) dη,

(4.3)

e−iϱ ξ2η Un,p2,p3,s

 

 

 

 

 

 

 

 

 

 

 

 

где ξ1 =

 

eB(x − x) и ξ2 = eB(y − y). Далее, представив функцию

˜

в виде двумерного интеграла Фурье, мы можем образовать

F (ξ1, ξ2)

в интеграле In,p(+)3,s(x, x) фазовый множитель

e−ip (xx), где величина

pµ=(En, p1, p2, p3) может интерпретироваться как 4-импульс частицы.

e

Прямое и обратное Фурье преобразования функции F (ξ1, ξ2) удобно

18

выбрать в виде:

F (p1, p2) = (2π)2

e−i(ξ1p1+ξ2p2)/√eB F (ξ1, ξ2) 12,

(4.4)

 

1

−∞

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eB

.

 

 

(4.5)

F (ξ1, ξ2) =

ei(p1ξ1+p2ξ2)/eB

F (p1, p2)

 

 

e

−∞

 

 

 

 

 

 

 

 

 

 

 

 

dp1dp2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В итоге интеграл (4.2) запишется через Фурье-образ функции F (ξ

, ξ

)

следующим образом:

 

eΦ(x,x)

 

 

 

 

 

 

 

 

 

 

 

dp1dp2

e

1

2

 

In,p3,s(x, x) =

 

e

 

 

 

 

 

 

 

 

 

.

(4.6)

 

2EnLyLz

ip (x

)

F (p1, p2) eB

(+)

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

Так как интерес представляет не сам интеграл (4.1), а его подынтегральная функция, то для нее получим:

 

 

 

(+)

 

eΦ(x,x)

 

 

 

 

 

 

 

 

dp1

 

 

 

 

 

e

 

 

 

,

(4.7)

ψn,p2,p3,s(x) ψn,p2,p3,s(x) =

2EnLyLz

ip (x

)

F (p1, p2) eB

(+)

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (p1, p2) =e(ξ1/2−η−ϱp2/√

 

 

)ξ2 e−iξ1p1/

 

×

 

 

 

 

 

 

 

eB

eB

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

 

(+)

 

 

 

(+)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×[Un,p2,p3,s

(η) Un,p2,p3,s(η − ξ1)]

 

 

 

.

 

(4.8)

 

 

 

 

 

(2π)2

 

Таким образом, можно определить матрицу плотности ρ(+)n (p) дираковской частицы в импульсном представлении в следующем виде:

2π

(4.9)

ρ(+)n (p) = eB F (p1, p2).

Заметим, что интеграл по переменной ξ2 легко берется и дает сле-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дующий результат: 4π δ

 

ξ1 2η − 2ϱ p2/

eB , где δ(x) – δ-функция

легко вычислить интеграл по переменной ξ

.

Дирака, что позволяет (

 

 

 

 

 

)

 

 

 

 

1

 

Таким образом, для функции F (p1, p2) получим:

 

 

F (p1, p2) = e−i 2ϱ p1p2/eB e−i 2p1η/√

 

×

 

 

 

 

 

 

 

 

eB

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

(+)

 

 

 

(+)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×[Un,p2,p3,s

(η) Un,p2,p3,s(−η −

2ϱ p2/eB)]

 

.

(4.10)

π

19

Выпишем явно матрицы

Un,p(+)

 

 

 

 

(+)

 

 

 

 

для каждой из поляриза-

,p

,s

U

 

 

 

 

ций s по отдельности:

[

 

2

 

3

 

 

 

n,p2,p3,s]

 

 

 

 

 

F+ϱ(p1, p2) = F (p1, p2)

 

 

 

 

 

 

 

iab/2

 

 

 

 

 

 

= (1)n

e

 

 

 

dη e−iaη ×

(4.11)

π

 

 

 

 

 

 

 

 

 

s=ϱ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(η + b)

 

 

 

 

 

 

 

 

 

 

 

 

W

ϱ

W

 

χ

(η) χ

V

 

V

 

 

 

χ

n−1

(η) χ

n−1

(η + b) +

×[[

 

 

ϱ] n

 

n

 

 

 

[ −ϱ

 

 

−ϱ]

 

 

+ [Wϱ

 

−ϱ] χn(η) χn−1(η + b) [V−ϱ

 

ϱ] χn−1(η) χn(η + b)],

V

W

F−ϱ(p1, p2) = F (p1, p2)

 

 

 

 

 

 

 

 

iab/2

 

 

 

 

= (1)n eπ

 

 

dη e−iaη ×

(4.12)

 

 

 

 

 

 

 

s= ϱ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[[ ] [ ]

× Vϱ V ϱ χn(η) χn(η + b) − W−ϱ W −ϱ χn−1(η) χn−1(η + b) +

]

[ ] [ ]

+ W−ϱ V ϱ χn−1(η) χn(η + b) − Vϱ W −ϱ χn(η) χn−1(η + b) ,

где a=2p1/ eB, b=2ϱp2/ eB, и было использовано следующее свойство функций Эрмита: χn(−η) = (1)n χn(η). В общем случае интегралы такого типа сводятся к обобщенным полиномам Лагерра:

Lnk (x) =

x−kex dn

[xn+ke−x] ,

(4.13)

n! dxn

а именно,

e−iaη χn(η) χm(η + b) =

−∞

= (1)n−m eB 2 m−n m!/n! (b + ia)n−m e(iab−c2)/2 Lnm−m(c2) ,(4.14)

где n≥m и c2=(a2 + b2)/2. Поскольку в (4.11) и (4.12) функции χn(η) входят с индексами, различающимися не более чем на единицу, то более удобным оказывается использование следующих соотношений:

e−iaη χn(η) χn(η + b) =

 

e(iab−c2)/2 Ln c2

 

(4.15)

eB

,

−∞

(

)

 

20

Источник: https://studfile.net/preview/16723460/