Dirac equation and the Melvin Metric

A relativistic wave equation for spin 1/2 particles in the Melvin space-time, a space-time where the metric is determined by a magnectic field, is obtained. The effects of very intense magnetic fields in the energy levels, as intense as the ones expected to be produced in ultra-relativistic heavy-ion collisions, are investigated.


I. INTRODUCTION
In the past few years systems with extreme magnetic fields (B > 10 15 G) have been proposed to exist. In the magnetar analysis, for example, magnetic fields of the order of 10 15 G at the surface [1], [2] and 10 17 G at the center are expected to exist. In ultra-relativistic heavy ion collisions, at √ s= 200 GeV (at RHIC), the magnetic field is expected to reach values as high as 10 19 G [3]- [5] and at LHC, at √ s= 7 TeV, 10 20 G.
When studying particles inside this systems (eletrons inside stars, and so on), the usual way is to solve the Dirac equation [6], [7], for example, and find the energy levels, and in general, the Minkowski space-time is considered. A question that is quite intriguing, is if when the magnetic field reach these values, some effect of the structure of the space-time may be observed. This is the purpose of this work.
Including the magnetic field in the metric is not a trivial question but some solutions exist, as the Melvin metric [8]- [12], where a magnetic universe, with a magnetic field in the z direction, is considered, or the Gutsunaev solution [13], [14], for a magnetic dipole.
Of course it is always possible to study an arbitrary shape of the magnetic field and solve (at least numerically) the resulting Einstein equations, but with the objective of finding analytical results, and with this procedure, exploring this effect in a first approximation, we will consider the Melvin metric and find the wave equations for Dirac particles subjected to a magnetic field inside this metric.
So, this paper will show the following contents: In section II, a brief review of the formulation of the Dirac equation in curved spaces will be made. In section III the wave equation in the Melvin metric will be worked out. The results and conclusons of this work will be shown in section IV.

II. DIRAC EQUATION IN CURVED SPACES
In this section a brief review about the wave equation for spin 1/2 particles in curved spaces will be made. The basic formulation and the equations that will be needed in the next sections will be shown.
A fundamental characteristic of the Dirac equation is its invariance under Lorentz transformations, so, when studying Dirac particles in curved spaces, it is interesting to preserve this aspect. One way to do that is by using the tetrads e (a) µ that may be defined in order to satisfy the equation where η (a)(b) is the Minkowski tensor, that represents a flat space-time, and g µν the metric tensor related to a space-time that possesses an arbitrary geometry [15]- [17].
Observing eq. (1), we may note that the tetrads may be used in order to project vectors from the curved space-time in the flat space-time with the expression that relates the form of a vector in different space-time geometries.
As it was said before, the Dirac equation when written in the Minkowski geometry posseses Lorentz symmetry, and we want that the equivalent equation, written in a curved space-time possesses the same characteristic. In fact, if we study the behavior of the elements of the Dirac equation under transformations that preserve the desired symmetries, we may understand which are the quantities that need to be added or modified in the equation.
We must note that a spinor transforms according to where ρ (Λ) = 1 + 1 2 iε (a)(b) Σ (a)(b) , and Σ (a)(b) is the spinorial representation of the generators of the Lorentz transformation, written in terms of the γ (c) matrices, Σ (a)(b) ≡ 1 4 i γ (a) , γ (b) [18]. The idea is to construct a covariant derivative ∇ (a) ψ that is locally Lorentz invariant, that means that we need to impose the transformation condition The usual way to obtain the form of the covariant derivative operator of the spinor is by supposing the combination and considering the form of the operator Ω µ transformation If we choose the combination of terms or in an equivalent form with the term Γ (a)µ(b) defined as where Γ ν µλ are the Christoffel symbols, eq. (6) is satisfied. Consequentlly eq. (5) satisfies the transformation condition demanded in (4) and then we get the final form of the covariant derivative operator We must remark that there are different definitions of covariant derivatives of tensors in the literature when studying the covariant Dirac equation, as for example in [19]- [21].
Considering the Dirac equation in a flat space-time and replacing the conventional derivative operator by the one obtained in (10) we obtain the desired wave equation for spin 1/2 particles in a curved space-time If the particle is submited to an external electromagnetic field, we may introduce this effect by a minimal coupling and the scalar action that leads to this equation is given by It is usual to define the term γ µ = e µ (a) γ (a) as a Dirac matrix in a given curved space-time and it is easy to verify that it satisfies the Clifford algebra

III. WAVE EQUATION IN THE MELVIN METRIC
If now we want to obtain the wave equations shown in the last section, in a spacetime that has its structure determined by a magnetic field, an useful idea wolud be to consider the Melvin metric [8]- [12]. The Melvin metric is a solution of the Einstein-Maxwell equations of the general relativity that represents a cylindrical magnetic universe. In his work [8], [9], Melvin considered a static magnetic field where its lines lie in cylindrical surfaces perpendicular to the radial direction, with intensity ∼ B 0 in the vicinity of the symmetry axis and falls as fast as B 0 /r 4 far away from the axis. This solution has wide-ranging applications in the literature, as, for example in the study of Kerr black holes [22], or in cosmology, where the possibility of the interaction of the magnetic field with the expansion of the universe may be considered, as for example in [23], [24].
The line element may be written in a cylindrically symmetric form (taking c=G= =1) where Λ (r) = 1 + 1 4 B 2 0 r 2 and B 0 is the magnetic field. This metric reflects the curvature of the space-time, determined by the existence of a magnetic field in the z direction. For this reason, a metric with axial symmetry is taken into account.
In the limit of vanishing the magnetic field we have Λ(r) = 1 and the equation (16) becomes that is the flat Minkowski space-time written in cylindrical coordinates. So, we may define and choose a diagonal tetrad basis e for which the equation (1) is satisfied, and then it is easy to determine its inverse form e µ (a) , We are interested in studying the effect of a magnetic field B 0 , that modifies the spacetime geometry, as it is shown in the Melvin metric. If we consider a Dirac particle inside this field, we may also investigate the effect of the minimal coupling in the wave equation for this particle considering the 4-potential A µ . In a flat space-time, a constant magnetic field in the z direction, B 0 , that may be related to an equivalent magnetic field in the Melvin metric by the equation (2) appears if the 4-potential has the only non-vanishing component given by Observing that the term γ µ Ω µ , in the curved space wave equation, relative to the tetrad (19) is given by equation (13) becomes In the limit Λ → 1, this equation reduces to the usual Dirac equation for a free particle in a flat space-time in a cylindrical coordinate system. So, we may interpret equation (23) as a generalization of a Dirac equation for a particle inside a magnetic field, an equation that has been extensively studied in the literature, in many contexts, as for example [25]- [29].
We are interested in observing the corrections in the energy spectrum, due to the alteration of the geometry of the space-time, determined in eq. (23) .

Making a transformation in equation (23)
we obtain a simplified form where the usual representation for the gamma matrices is considered that has no dependence in the z, t and φ variables. We will suppose a solution of the form where σ is the energy of the system, that assumes positive values for particles and negative values for antiparticles, p z is the momentum, m = ±1, ±2, ±3, ..., a quantum number, and M, the electron mass. So, equation (23) may be written in a explicit form In a first approximation we will solve the equation neglecting terms in higher orders of r, and then, the system of equations may be written as This set of equations may be decoupled by multiplying the first equation by the expression and then, using the second and the fourth to eliminate the terms containing the spinors R 1 , where b = mB 2 0 − qB 0 and m ′ = m − 2q/B 0 . In a similar way we derive equations for R 1 , R 2 , R 3 that may be resumed in the form where the positive sign refers to R 2 (r) and R 4 (r), and the negative one to R 1 (r) and R 3 (r) .
As we may observe, equation (38) is similar to the Schrödinger equation, and it may be written as may be identified as an effective potential and as we can see, the system has the form of an isotropic harmonic oscillator.
In fact, the solution of the equation (39) may be mapped into a 3-dimensional harmonic oscillator-like one in spherical coordinates. These solutions are given in terms of the associate Laguerre polynomials where S = br 2 2 and N i are normalization constants. The energy spectrum, relative to this solution is and using the definition of E given in eq. (39), we obtain the energy spectrum where s = ±1/2 is the spin quantum number of the particle and m = 1, 2, 3..., in order to keep the conservation of the parity. As it was pointed before, equation (41) is written with c = G = = 1. By making the conversion to the international system of units we have where ε 0 is the vacuum permissivity constant. In a system of units with = c = 1 we have were, in the last term inside the square root, the Planck mass m p appears scaling the magnetic field as (B 0 /m p ) 2 .

IV. RESULTS
Studying equation (41) we may recover some literature results. When the magnetic field and p z goes to zero, we obtain the expression for the rest energy of the particle Now, if the gravitational energy is neglected, the last term inside the square root vanishes and we have  Table 1. Observing the Table, we can see that the effect of the magnetic field appears when B increases and becomes important when these fields are intense.
The same calculations have been performed in a flat space-time with a minimal coupling.
The results are shown in Table 2. As we can see, the results are essentially the same, and, when comparing then with the ones found in Table 1, we can only find deviations when considering fields of the order of B ∼ 10 19 G and quantum numbers as large as n ∼ 10 30 .
These results are shown in Fig. 1. Observing these results, the conclusion that we may obtain, with a very good precision, is that the effect of the magnetic field in the metric is very small and may be neglected.
As it was explained before, eq. (41) is an approximation for small values of r, but when comparing the results with the exact numeric calculations that have been performed, up to the precision shown in Tables 1 and 2 the conclusions that we obtain are essentially the same, and then we may use this equation as a very good approximation. The corrections from terms with higher powers of r appear beyond the precision shown in these tables. This fact may be seen if one observes the approximations that have been made when eq.(32)- (35) have been obtained. Terms of the order (B 0 r) 2n , where n is an integer, appears. These terms increase with r 2 , and reaches the maximum values with some r M , that determines the size of the considered system. For heavy-ion collisions, for example, (B 0 r M ) 2 ∼ 10 −33 , for pulsars, 10 −11 . The greater value obtained was for magnetars, 10 −7 , and the effect of these corrections, as it was said before, are always beyond the precision shown in the tables, so, we may conclude that the accuracy of the approximation is good, even for extremely large magnetic fields and levels.
One remark that must be made is that when considering a covariant Dirac equation, a particle with gyromagnetic ratio g = 2 is taken into account. This is a good approximation for electrons inside stars, for example. But when considering higher energy processes, as high-energy collisions, deviations from this value are proposed to exist, and then, a way to study this question is to consider extensions of the covariant Dirac equation, as for example as it has been made in [21]. This kind of approach is left for future works.
Another aspect that may be taken into account is the Melvin metric. It is clear that (23) without a minimal coupling (solid line). Table 1: Energy levels in the Melvin metric with a minimal coupling, from the value of the magnetic field found in the Earth, up to the one expected to be produced in ultra-relativistic heavy-ion collisions.
many of the systems that have been studied does not have the magnetic field in the form of the one that determines this metric, but in some regions, with intense magnetic fields, at least as a first qualitative approximation, these results must be correct, and we expect that more careful calculations, with the metric determined by different shapes of the magnetic fields, confirm our results.