Thesis - Report - Zayid Ahmed
Thesis - Report - Zayid Ahmed
Abstract
Charged particles are very sensitive to the magnetic field as it changes
their properties and particles behave differently according to the field.
Here we have studied how the Dirac solutions and the energy eigenval-
ues change in presence of uniform background magnetic field and their
gauge dependence. After that we have studied the method given by Ju-
lian Schwinger to calculate the propagator. The fermionic propagator
changes in presence of magnetic field. We have discussed it’s gauge de-
pendence and how it changes according to the strength of background
magnetic field.
2
Contents
1 Introduction 4
3 Schwinger’s method 9
3.1 photon propagator . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.2 fermionic propagator . . . . . . . . . . . . . . . . . . . . . . . . . 14
5 Conclusion 16
3
1 Introduction
As a charged fermion, we consider a electron being the particle having the
largest specific charge i.e, being the most sensitive to the external field influence.
The large hadron collider (LHC) and the relativistic heavy-ion collider (RHIC)
can produce a very strong magnetic field at the very early stages of collision. A
very strong magnetic field exists at the core of magnetars, and the same may
have existed at the early stages of universe. In light of this we study different
types of elementary particle decays and scattering processes. Theoretically the
solutions of Dirac equation in vaccum are well known, with the Dirac spinors
for positive and negative energy along with there corresponding completeness
relations.
In presence of background magnetic field the Dirac spinors are going
to change as they depend on the spatial coordinates due the gauge used to
solve the Dirac equation and also depend on the momentum. As a result, their
corresponding completeness relation are also going change. These realistic fields
are very complicated but for simplicity many times we assume these fields are
uniform. The Dirac equation can be exactly solved in presence of uniform
magnetic field, then we can proceed to quantize these solutions and calculate
elementary particle decays and scattering cross-section. We use the uniform
background magnetic field along z-axis. The dispersion relation we get after
solving the Dirac equation is En2 = p2z − m2 + 2n|eB|, and n = 0 state is non-
degenerate for both positive and negative energy solutions. Thus the ground
state is spin polarised. In the case of strong magnetic field (|eB| >> m2 ),
electron prefer to lie in the lowest Landau level, and electron will have only
longitudinal momentum. In the weak magnetic field the higher Landau levels
contribute. The dispersion relation is gauge independent, but the spinors and
completeness relation are gauge dependent for both positive and negative energy.
We use the Schwinger’s method to calculate the propagator for charged
scalar field and fermionic field. This method has some nice properties, as we are
not fixing the gauge in the beginning of calculations. As a result we end up with
an expression for the propagator with an arbitrary gauge. In the case strong
magnetic field it can be shown that one part of the propagator corresponds to the
contribution of the ground Landau level and the other part corresponds to the
contribution of higher Landau levels. We will see that propagator has a phase
factor which plays an important role in transforming the two-point function.
4
2.1 Choice of gauge
The background magnetic field of magnitude B along z-direction can be
fixed by choosing background gauge field AµB in many ways:
or
A0B = AyB = AzB = 0, AxB = −yB + b (2)
or
1 1
A0B = AzB = 0, AyB = xB + c, AxB = − yB + d (3)
2 2
Where a, b, c, and d are constants. AµB is classical background field and not a
quantized dynamical field. Where x and y are coordinates not 4-vectors. We
assume gauge configuration as given in Eq.(1) with b = 0 in our calculation.
We will use Dirac-Pauli representation of Dirac matrices.
Using the value of vector potential and evaluating RHS of the above equation,
∂
(E 2 − m2 )φ = [(−∇2 − (QB)2 y 2 − QB(2iy + σ3 )]φ (10)
∂x
5
As we can see that the variables x and z do not appear in the equation except
through the derivatives, by variable separable method we can write the solution
as,
φ = eip.Xy0 f (y) (11)
f (y) is a 2-component matrix which only depend on y, and possibly some mo-
mentum components. we have used p.Xy0 = px x + pz z, where px and Pz are the
eigenvalues of the momentum along x and z coordinates respectively.
σ3 appears in the equation which is a diagonal Pauli matrix. There will be two
independent solution of f (y), which can be taken to be the eigenstates of σ3
with eigenvalues s = −1, +1. We can choose two independent solution in the
form
F+ (y) 0
f+ (y) = , f− (y) = (12)
0 F− (y)
where σ3 fs = sfs , using this in Eq.(10) equation satisfied by Fs is,
d2 FS
+ (QBy + px )2 FS + (E 2 − m2 − p2z + QBs)FS = 0 (13)
dy 2
We will use a dimensionless variable,
p px
ξ = |Q|B(y + ) (14)
QB
which transform the equation in the form,
d2
[ − ξ 2 + as ]Fs = 0 (15)
dξ 2
where,
E 2 − m2 − p2z + QBs
as = (16)
|Q|B
Eq.(15) is Hermite’s equation and solution exists provided as = 2ν + 1 where
ν = 0, 1, 2, ....This gives us the energy eigenvalues
Nν are normalisation
p|Q|B 1/2
Nν = √ (19)
ν!2ν π
The functions Iν satisfy the completeness relation
X p
Iν (ξ)Iν (ξ∗ ) = |Q|Bδ(ξ − ξ∗ ) = δ(y − y∗ ) (20)
ν
6
Q was arbitrary.we will discuss two cases one for Q = −e (electron case) and
second Q = +e (positron case).
Q = −e Case
In−1 (ξ) 0
0 √ In (ξ)
U+ (y, n, py0 ) =
Enp−m
z
In−1 (ξ) , U− (y, n, py0 ) = 2neB In−1 (ξ)
(24)
√ En +m
− 2neB
En +m In (ξ) − Enp+m
z
In (ξ)
Q = e Case
The energy eigenvalues are,
7
With the same analogy as for the positive energy solution, the negative
energy solution for Dirac equation can be written as:
where X µ denotes the spacetime coordinate. And the spinors are given by,
pz √2nB
In−1 (ξ˜)
˜
En
√ +m En +m In−1 (ξ )
2nB In (ξ˜) pz ˜
, V+ (y, n, py0 ) = − En +m In (ξ )
V− (y, n, py0 ) = E +m
nI ˜
(27)
(ξ )
n−1
0
0 In (ξ˜)
8
solution s has a definite value and for all other higher Landau levels solution do
not have definite σz eigenvalue.
Considering the case of strong magnetic field we require to work n = 0
solutions and we can estimate the the strength of that field. Let say E is
the electron energy of the system and B is the strength of magnetic field. If
2neB > E 2 + m2 for a positive value of n then from the dispersion relation
E 2 = m2 + p2z + 2neB we can see that p2z becomes negative, that makes momen-
tum along z-axis imaginary, which is impossible. Thus for strong magnetic field
that makes 2eB greater than the square of electron energy minus the rest mass
square of the electron only n = 0 levels contributes to the energy levels and
only those corresponding wave functions should be used for the calculations. If
the electron energy of the system is of the order of 1M eV the for the magnetic
field magnitude greater than 1014 Gauss we should only have the n = 0 level
contributing towards the energy.
Thus, when the magnetic field magnitude is large only the lowest
Landau level will contribute to the electron energy but if the magnetic field
magnitude is low the higher Landau levels will start contributing.
3 Schwinger’s method
We will use the Schwinger’s method to calculate the propagator for scalar and
fermionic field in presence of uniform background magnetic filed.
We define the Green’s function as
H(x, ∂x )G(x, x0 ) = δ 4 (x − x0 ) (34)
The idea is to consider H as a hamiltonian that dscribes the proper time evolu-
tion of the system.
hx|x0 i = δ 4 (x − x0 )
[xα , pβ ] = −ιg αβ ,
pµ = ι∂µ
We introduce the unitary evolution operator which satisfy the Srondinger equa-
tion given by’
∂
i U (x, x0 : s) = H(x, p)U (x, x0 : s) (35)
∂x
9
where s is the proper time parameter. With boundary conditions;
Wew have
U (x, x0 ; s) = hx| e−ιHS |x0 i = hx| U |x0 i
(37)
Z 0
G(x, x0 ) = −ι dsU (x, x0 ; s) (38)
−∞
ι∂s = hx| H(x, p)U |x0 i = hx| U (s)U (s)H(x, p)U (s) |x0 i
.
Then the equation for U (x, x0 ; s) becomes an ordinary linear differential equation
which may be integrated as
h Z s i
U (x, x0 ; s) = C(x, x0 )exp − ι ds0 F (x, x0 ; s) (40)
where
H(x, p) = Πµ Πµ − m2
10
is the Hamiltonian and Πµ = pµ − eAµ (x)
To compute hx(s)| H(x(s), p(s)) |x0 (0)i we require to solve the equations motion
for xµ and Πµ , i.e. Hiesenberg
dxµ (s)
= ι[H, xµ ] (42)
ds
dΠµ (s)
= ι[H, Πµ ] (43)
ds
Let us compute the commutators of the above equations,
[H, xµ ] = [Πµ Πµ − m2 , xµ ]
= 2ιΠµ
[H, Πµ ] = [Πν Πν − m2 , Πµ ]
= Πν [Πν , Πµ ] + [Πν , Πµ ]Πν
= Πnν (−ιeFµν ) + (−ιeFνµ )Πν
= ιeΠν Fµν + ιeFµν Πµ
As we are working witht a constant uniform magnetic field along z-axis, in this
case ∂ ν Fµν = 0, i.e., ∂ ν (∂µ Aν − ∂ν Amu ) = 0 Equations of motion,
dxµ
= −2Πµ (44)
ds
dΠµ
= −2eFµν Πµ (45)
ds
Let us introduce some notations, since the space separates naturally in a parellel
and a perpendicular direction with regards to the magnetic field, for vectors aµ ,
bµ ,
aµ|| = (a0 , 0, a3 )
aµ⊥ = (0, a1 , a2 , 0)
Such that
(a · b)|| = a0 b0 − a3 b3
(a · b)⊥ = a1 b1 + a2 b2
(a · b) = (a · b)|| − (a · b)⊥
11
The equations of motion separates as;
dx⊥
= −2Π⊥
ds
dΠ⊥
= −2eF ij Π⊥
ds
dx||
= −2Π||
ds
dΠ||
=0
ds
where i, j = 1, 2
Let us look at the equations in the perpendicular direction,
x1 (s) Π1 (s)
x⊥ = , Π⊥ = (46)
x2 (s) Π2 (s)
dx⊥ (s)
= −2Π⊥
ds
dΠ⊥
= −2eBF Π⊥
ds
The solutions of these two equations are
12
therefore the Hamiltonian after doing the normal ordering i.e, we need to place
x(s) to the left and x(0) to the right, for that we will use the commutator
relation.
H⊥ = ΠT⊥ (s)Π⊥ (s) (51)
Let us look at the parallel part of the equations of motion, which are,
dx|| (s)
= −2Π||
ds
dΠ||
=0
ds
Thus the parallel part of the Hamiltonian is,
H|| = ΠT|| (s)Π|| (s) (52)
The Hamiltonian is;
H = H⊥ − H|| (53)
Since all the operators in H are in the correct order we can write,
2
−(eB)
hx0 (s)| H |x00 (0)i = (4 sin 02 0 00 00 1 02 0 00 002 2 0 00
2 (eBs)(x⊥ − 2x⊥ x⊥ + x⊥ ) − ιeB cot eBs + 4s2 (x|| − 2x|| x|| + x|| ) − m ) hx (s)|x (0)i
(54)
Thus the evaluation operator is,
C(x0 , x00 )
−ιeB 0 ι
hx0 (s)|x00 (0)i = exp (x − x00 )2⊥ cot eBs + (x0 − x00 )2|| + ιm2 s
s sin eBs 4 4s
C(x0 , x00 ) is the phase factor, we can compute this function. For that we need
to compute hx0 (s)| Πµ (s) |x00 (0)i and hx0 (s)| Πµ (0) |x00 (0)i, these two will give us
set of equations,
∂ e
(ι 0µ + eA(x0 ) + Eµν (x0 − x00 )2 )C(x0 , x00 ) = 0 (55)
2
and
∂ e
(ι 0µ + eA(x0 ) + Eµν (x0 − x00 )2 )C(x0 , x00 ) = 0 (56)
2
After solving these we can find function C(x0 , x00 ) , which is,
C(x0 , x00 ) = ιφ(x0 , x00 )
x0
( Z )
1
whereφ(x0 , x00 ) = exp ιe dxµ (Aµ (x) + Fµν (x0 − x00 )ν
x00 2
where φ(x0 , x00 ) is the phase factor and the other part is translationally invariant.
13
3.2 fermionic propagator
We will follow the same method here, the hamiltonian for a fermionic field is,
e
H = Πµ Πµ − m2 − σµν F µν (58)
2
The equations of motions are,
dxµ(s)
= −2Πµ
ds
dΠµ(s) e
= −2eFµν Πν − ι∂ ν Fµν − ∂µ Fµρ σ νρ
ds 2
For constant magnetic field, µ F µν = 0, we will end up with equations similar to
that for scalara field, which are;
dxµ(s)
= −2Πµ
ds
dΠµ(s)
= −2eFµν Πν
ds
Following the same notation as we did for scalar field propagator, we will sep-
arate the parallel and perpendicular parts of the equations of motion and the
calculate the Hamiltonian accordingly doing all the time ordering. As we can
see that this Hamiltonian has an extra term which is,
−e
σµν F µν = eBΣ3
2
where Σ3 = ιγ1 γ2 , B is the magnitude of magnetic field. The Hamiltonian is;
C(x0 , x00 )
eB 0 ι
hx0 (s) |x00 (0)i| = exp −ι (x − x00 )2⊥ cot eBs + (x0 − x00 )2|| −2 s + ιeBΣ3
s sin eBs 4 4s
(60)
where the phase factor is C(x0 , x00 ) = ιφ(x0 , x00 ) and
( Z 0 )
x
0 00 µ 1 0
φ(x , x ) = exp ιe dx (Aµ (x) + Fµν (x − x )6ν
x00 2
14
Propagator is,
Z ∞
0 00 0 00 1
exp (−2 s + BsΣ3
S(x , x ) = φ(x , x ) ds
0 s sin eBs
(62)
(" !#)
ι eBs
× exp − (x0 − x00 )2|| − (x0 − x00 )2⊥
4s tanh eBs
(63)
" #
1 eBs
× m+ (γ · x0|| − γ · x00|| ) − exp{−ιeBsΣ3 }(γ · x0⊥ − γ · x00⊥ )
2s sin eBs
(64)
In momentum space:
Z
ds exp{ιk · x}S(x0 − x00 )
S(k) =
Z ∞ (" #)
2 tanh eBs 2 2
= −ι ds exp ιs (γ · k|| ) − (γ · k⊥ ) − m
0 eBs
1
× (1 + γ1 γ2 tan eBs)(γ · k|| + m) − 2
γ · k⊥
cos (eBs)
15
particles and their two-point functions become background dependent when we
work in background field. That is why phase factor arises in the two-point func-
tion.
we need to consider that part of the propagator, s is time parameter. For the
weak field limit, we can demonstrate propagator into an infinite series in power
of B, and terms with lowest power of B are more important. And as B goes to
zero, we get back to the normal fermion propagator in absence of background
field.
5 Conclusion
The presence of background magnetic field effects the fermionic fields and the
propagator, this in turn effect the different types of scattering and decaying
processes. Presence of the magnetic field also helps us to understand the pro-
cesses more clearly as it separates the different type of particles. The processes
like γ− > νν − , e− e+ − > µ− µ+ , e− µ− − > e− µ− etc. can be studied using
the Dirac spinors and the fermionic propagator. The scattering cross-section
and the decay rates should be gauge independent when we calculate using these
calculations.
16
References
(1) K. Bhattacharya, Solution of the Dirac equation in presence of an uniform
magnetic field, 2007, arXiv:0705.4275 [hep-th].
(4) C. Itzykson and J. Zuber, Quantum Field Theory, Mc Graw Hill press,
International student edition, Second edition.
17