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

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

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

−∞

= eBn/2 e(iab−c2)/2 (b + ia)/c2 [Ln(c2) − Ln−1(c2)] , (4.16)

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

−∞

 

 

 

 

[

(

)

(

)]

 

=

0eBn/2 e(iab−c2)/2

(−b + ia)/c2

 

Ln c2

 

− Ln−1 c2

1

,(4.17)

где Ln(x) ≡ Ln(x) – полиномы Лагерра, и учтено свойство Ln−1(x) =

[

]

 

 

 

 

 

 

 

(n/x) Ln−1(x) − Ln(x) . Таким образом, после интегрирования по переменной η получаются следующие вклады в матрицу плотности от разных поляризаций:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

e−u/2{Ln(u) Wϱ W ϱ

− Ln−1(u) V−ϱ V −ϱ

}

 

 

 

 

π

 

 

 

e u/2

[

 

 

 

 

 

 

]

[

 

 

 

 

 

 

]

 

 

2n

 

 

 

 

 

 

 

 

 

 

+ (1)n

 

 

 

 

[Ln(u) − Ln−1(u)] ×

 

 

 

 

 

 

 

 

(4.18)

π

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[

 

 

 

]

 

 

 

 

 

[

 

 

 

 

]

 

 

 

 

×{(ϱp2 + ip1) Wϱ

V

−ϱ + (ϱp2 − ip1) V−ϱ

W

ϱ

 

 

},

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

e−u/2

{Ln(u) Vϱ V ϱ

− Ln−1

(u) W−ϱ

W −ϱ

}

 

 

π

 

 

 

e u/2

[

]

 

[

 

 

 

 

 

 

 

]

 

 

2n

 

 

 

 

 

 

 

 

 

+ (1)n

 

 

 

 

× [Ln(u) − Ln−1(u)] ×

 

 

 

(4.19)

π

 

u

 

 

 

 

 

 

 

 

 

 

 

[

 

 

 

 

]

 

 

 

 

[

 

 

 

 

 

 

 

]

 

×{(−ϱp2 + ip1) W−ϱ

V

ϱ (ϱp2 + ip1) Vϱ

W

−ϱ

},

 

где u = 2(p21 + p22)/eB. Вычисление билинейных комбинаций биспиноров W±ϱ и V±ϱ с использованием их явного вида (3.27) и (3.28) приво-

дит к результату:

 

p + m + 2eBn + m2

] Π±ϱ,

W±ϱW ±ϱ = 2

[(1 + 2eBn + m2 )

 

 

 

1

 

 

m

 

(b

)

 

2eBn

 

 

 

V±ϱV ±ϱ = 2

[(1 − √2eBn + m2 )

 

p

+ m − √2eBn + m2

]

Π±ϱ,

 

 

1

 

 

 

m

(b

)

 

2eBn

 

 

 

21

[

 

 

 

 

]

[

 

 

 

 

 

−ϱ

]

 

 

 

 

 

 

 

 

2

 

 

 

 

(ip1 − ϱp2) V−ϱ W ϱ

 

(ip1 + ϱp2) Wϱ V

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

b

 

 

 

 

 

 

 

=

 

 

 

 

eBn

[p − iϱ

p (pφγ)

] ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2eBn + m2

[

 

 

 

 

]

[

 

 

 

 

 

 

]

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

(ip1 + ϱp2) Vϱ W −ϱ

 

(ip1 − ϱp2) W−ϱ V ϱ

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[p +

 

 

b

] .

 

 

 

 

 

 

 

 

eBn

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

p (pφγ)

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2eBn + m2

Подстановка этих выражений в формулы

(4.18), (4.19) приводит к сле-

 

 

 

b

 

 

 

 

дующему выражению для матрицы плотности фермиона в импульсном представлении:

 

{ [(1 +

 

m

 

 

 

ρn,s(+)=ρ(p) = (1)n e−u/2

 

 

 

)p + p˜ + m] Πϱ Ln(u)

 

p˜

 

 

 

 

 

 

 

 

m b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

1

(u) +

 

 

 

 

 

 

 

 

 

 

 

[(1 p˜ )p − p˜ + m] Π−ϱ Ln−1

 

 

[

 

 

 

 

 

 

 

 

 

 

]

 

 

 

 

b

 

 

 

 

 

 

 

b

Ln−1(u)},

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

+ 2 p − iϱ

 

 

(pφγ)

 

(4.20)

 

 

p˜

 

 

 

{ [(1

 

 

 

 

m

 

 

 

ρn,s(+)=−ρ(p) = (1)n e−u/2

 

 

)p − p˜ + m] Πϱ Ln(u)

p˜

 

 

 

 

 

 

 

 

m b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

1

 

(u) +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[(1 + p˜ )p + p˜

+ m] Π−ϱ Ln−1

 

 

[

b

 

 

 

 

 

 

 

 

 

 

]

 

 

˜

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

+ 2

p + p˜ (pφγ) Ln−1(u)}.

 

(4.21)

Здесь pb = (pΛγ) = Enγ0 − p3γ3, pb = (pΛγ) = p1γ1 + p2γ2, (pφγ) =

˜

p2γ1 − p1γ2, φµν = Fµν/B и φ˜µν = Fµν/B — безразмерные тензор и ду-

˜

альный тензор электромагнитного поля, Λµν = (φφ)µν, Λµν = (φ˜φ˜)µν, Πϱ — оператор проекции спина частицы на направление магнитного поля (см. формулу (2.10)). Отметим, что входящие в (4.20) и (4.21) структуры pb , pb , (pφγ), p˜ , Πϱ, а также аргумент полиномов Лагерра u = 2p2 /eB = 2(p21 + p22)/eB являются инвариантами относительно преобразований Лоренца вдоль вектора напряженности магнитного поля. Эффективное разбиение 4-мерного пространства в постоянном однородном магнитном поле на два ортогональных подпространства —

22

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

После суммировании по поляризациям, матрица плотности приво-

дится к виду:

ρn(+)(p) =

ρn,s(+)(p) = (1)n 2 e−u/2 ×

 

 

s=±ρ

)(

)

 

[(

×pb + m Πϱ Ln(u) Π−ϱ Ln−1(u)

(4.22)

]

+ 2 pb L1n−1(u) .

Отметим, что выражение для ρ(+)n (p) может быть получено при использовании любого из проектирующих операторов, поскольку не содержит информации о спиновых свойствах заряженной дираковской частицы. Чтобы убедиться в этом, построим матрицу плотности, просуммированную по поляризациям частицы, на основе решения со спинорами (3.28), (3.27), отвечающими проекционному оператору Σ3. Для этого удобно воспользоваться следующими свойствами биспиноров, которые могут быть получены непосредственным вычислением:

 

 

 

 

 

 

 

(

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

(4.23)

 

 

Ws

W

s + Vs

V

s = p + mI Πs,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−s

 

 

 

 

 

 

 

 

 

 

 

(4.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(sp2 + ip1)

V W

W

V

−s

 

 

 

eBn p

 

,

 

 

 

 

sb

s

 

 

= 2

 

 

 

 

 

 

s=±1

[

 

 

 

 

 

 

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p b= p1γ1 + p2γ2.

 

 

Πs = (I + isγ1γ2) /2,

p

 

= Enγ0

p3γ3,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

I

– единичная матрица в пространстве Дирака, а γ

µ

– матри-

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

b

 

 

 

цы Дирака. Используя эти свойства и суммируя вклады в матрицу плотности от различных поляризаций, получим:

 

 

 

 

 

 

 

(+)

 

 

 

 

 

eΦ(x,x)

 

 

 

 

 

 

 

 

 

 

dp1

 

 

 

 

 

(+)

 

 

 

 

 

 

 

ip (x

x

)

(+)

 

 

 

 

ψn,p2,p3,s

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

 

e

 

 

 

ρn

(p)

 

,

s=

±

1

2EnLyLz

 

2π

 

 

 

 

 

{

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ(+)n (p) = (1)n 2 e−u/2

 

 

+ mI Πϱ

2 n

Ln(u)

 

(4.25)

 

 

p

 

 

 

p

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

[(b

 

 

 

)

 

 

 

2bn

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[(p + mI) Π−ϱ

 

 

 

 

p ] Ln−1(u) },

 

 

 

 

 

 

 

 

 

 

u

где

ϱ

– знак заряда частицы,

Φ(x, x

) = eB(x + x )(y

y

)/2

, и под-

 

 

 

b

 

 

 

 

 

 

b

 

 

 

разумевается, что L1(u) 0. Принимая во внимание рекуррентные соотношения на обобщенные полиномы Лагерра:

 

(u) (n + k) Lnk

d

(4.26)

n Lnk

1(u) = u du Lnk (u) = −u Lnk+11(u),

23

матрицу плотности (4.25) можно также представить в следующих фор-

мах:

[( )( )

ρ(+)n (p) = (1)n 2 e−u/2 p + mI Πϱ Ln(u) Π−ϱ Ln−1(u)

b

b

 

 

 

 

 

 

 

 

d Ln(u)

,

(4.27)

2 p

 

 

 

du

 

]

 

ρ(+)n (p) = (1)n 2 e−u/2[(p

+ mI1)(Πϱ Ln(u) Π−ϱ Ln−1(u) ) +

b

 

 

L

n−1

]

 

(4.28)

+ 2 p

 

 

 

(u) .

 

b

Как и следовало ожидать, это выражение совпадает с ранее полученным выражением (4.22), где был использован проекционный оператор µˆ3 (3.31).

При «наивном» суммировании этой матрицы плотности по n, то есть в предположении, что в пределе слабого поля (B → 0) дискретный спектр энергий заряженной частицы переходит в непрерывный

(En

E = p2 + m2), получаем стандартное вакуумное выражение:

(

 

 

 

) ∑

 

 

n=0 ρn(+)(p) = 2 e−u/2

p + m

n=0(1)n[Πϱ Ln(u) Π−ϱ Ln−1(u)] +

 

+ 4 e−u/2 p

b(1)nLn1

1

(u) = p + m − p = p + m, (4.29)

 

b

 

 

b

b b

n=1

где было использовано правило суммирования обобщенных полиномов

Лагерра:

 

 

 

 

 

 

(n + k)!

 

ex

 

 

 

 

 

 

(4.30)

(1)n

n! k!

Lnm+k(2x) =

2k+m+1

Lkm(x)

n=0

 

 

 

 

 

 

при значениях k = 0 и m = 0, 1.

В случае нерелятивистской частицы, пренебрегая всеми поперечными к полю компонентами импульса и полагая pˆ = mvˆ , получим из (4.20) и (4.21):

 

 

µ

ρn,s(+)=ϱ(p) = (1)n 2 e−u/2Ln(u) m 1 + v Πϱ,

)

 

 

 

 

2

(

1)+ v

 

 

 

ρn,s(+)=−ϱ(p) = (1)n+1 2 e−u/2Ln−1((u) mb

Π−ϱ, (4.31)

где

v

 

 

 

 

 

 

 

среды вдоль на-

 

= (1, 0, 0, v)/ 1

v

 

 

 

 

 

 

 

 

— 4-скорость движения b

 

правления поля. Нетрудно убедиться, что матрица плотности (4.31)

24

описывает состояние с определенной проекцией оператора дираковского спина Σ3 на направление магнитного поля.

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

 

p32

eBn

− gϱ˜

eBs

 

(4.32)

En,s = m +

 

+

 

 

,

2m

m

2m

где число n нумерует уровни Ландау, g˜ — аномальный магнитный момент в ядерных магнетонах для нуклонов и магнетонах Бора для электронов.

Просуммированная по s матрица плотности, соответствующая решению уравнения Дирака с отрицательной энергией (3.30), получается из (4.22) формальной заменой pµ → −pµ, что приводит к выражению:

 

( )

 

 

( )

 

 

e−iϱ Φ(x,x)

 

ip (x

x)

( )

 

dp1

 

 

 

 

 

 

 

 

 

ψn,p2,p3,s

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

 

e

 

ρn

(p)

 

,

s= 1

2EnLyLz

 

2π

±

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

{

p − m

 

 

 

 

 

 

 

 

 

ρn()(p) = (1)n 2 e−u/2

Πϱ Ln(u) Π−ϱ Ln−1(u) + (4.33)

 

 

 

 

 

 

(

 

)[

 

 

 

 

 

 

 

 

]

 

 

 

 

 

 

+b2 p

 

L1

 

(u) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n−1

 

}

 

 

 

 

 

 

 

b

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

Для полноты изложения, приведем известные матрицы плотности

для безмассового нейтрино с 4-импульсом k

µ

 

 

 

 

 

 

 

= (ω, k):

 

 

 

ψ(ν)

 

 

(ν)

 

e−ik(xx)

ρ(ν)

k ,

ρ(ν)

k

1 k

γ

,

 

 

ψ

 

k

(x)

 

k

(x ) =

 

 

(

)

 

( ) =

 

bµ(1

 

5)

(4.34)

 

2ωV

 

2

 

и для электронейтральной частицы с 4-импульсом P = (E, P )

и мас-

сой mN :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ(N)(x) ψ(N)(x) =

e−iP (xx)

ρ(N)(P ),

 

ρ(N)(P ) = P +mN ,

(4.35)

 

 

 

P,s

 

 

 

P,s

 

2EV

 

 

 

 

 

 

 

 

 

s=±1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25

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