Heun functions describing bosons and fermions on Melvin's spacetime

Employing a pseudo-orthonormal coordinate-free approach, the solutions to the Klein--Gordon and Dirac equations for particles in Melvin spacetime are derived in terms of Heun's biconfluent functions.


Introduction
The study of relativistic particles in static magnetic fields has a long history and is still attracting considerable attention, especially for cases where someone deals with curved manifolds.
Even though on Minkowski spacetime the relativistic behavior of an electron in various magnetostatic configurations is well understood (see, for example, the Johnson and Lippmann's paper [1]), a weakness on curved spacetime regards the explicit gauge-covariant formulation.
Recently, when dealing with slowly rotating neutron stars which have been termed as magnetars [2], it has been assumed that their huge magnetic induction in the core and crust, B ∼ 10 14 − 10 15 (G), is affecting the spacetime geometry. A way out could be the search for general relativistic solutions with the magnetic field considered as a perturbation of the spherically symmetric background [3]. Another way is to assume that magnetized metrics, as the one belonging to the Melvin class [4,5] may be reliable candidates for describing these highly compact astrophysical objects with a dominant axial magnetic field [6].
Within a coordinate-dependent formulation, switching between canonical and pseudo-orthonormal basis, the above mentioned authors are integrating the system of four coupled first-order differential equations, in the first approximation, neglecting the terms in higher orders of the polar radial coordinate ρ. Their solutions are expressed in terms of generalized Laguerre polynomials, similarly to the case of the Dirac equation in cylindrical coordinates on a flat manifold [7].
In the present work, we are applying a free of coordinates method to analyze the Klein-Gordon and Dirac equations describing particles evolving in a Melvin's spacetime. Employing the Cartan's formalism, we are computing all the essential geometrical objects for writing down the corresponding matter fields and Einstein's equations.
It turns out that the SO(3, 1) × U (1)−gauge covariant Klein-Gordon equation can be exactly solved, its solutions being given by the Heun biconfluent functions [8][9][10]. The same happens with the approximate expression of the second-order differential system derived from the Dirac equation. The Heun functions, either general or confluent, are main targets of recent investigations and have been obtained for massless particles evolving in an Universe described by the metric function written as a non-linear mixture of Schwarzschild, Melvine and Bertotti-Robinson solutions [11].

The geometry
Recently, in [12], it was presented the procedure of transforming a known static symmetric solution to Einstein-hydrodynamic equations into a magnetized metric, by (nonlinearly) adding the magnetic field. In spherical coordinates, this has the general form with the metric functions g 11 and g 00 depending only on r and where, for the moment, B * is a parameter related to the magnetic field intensity.
In the pseudo-orthonormal Cartan frame corresponding to the metric (1), for the potential where B 0 is the strength of the magnetic field on the axis, the Maxwell tensor components, corresponding to a poloidal magnetic field with B ρ and B θ , are given by the relations pointing out a prolate (in shape) star. Once we assume g 00 = g 11 = 1, we can switch to cylindric coordinates {ρ, ϕ, z, t}, by ρ = r sin θ , z = r cos θ , so that the magnetized metric (1) turns into the simple Melvin expression with Within an SO(3, 1)−gauge covariant formulation, we introduce the pseudoorthonormal frame whose corresponding dual base is so that the metric (6) gets the Minkowskian form ds 2 = η ab ω a ω b , with η ab = [1, 1, 1, −1]. The first Cartan equation, with 1 ≤ b < c ≤ 4 and Γ a .
[bc] = Γ a .bc − Γ a .cb , leads to the following connection one-forms: where Λ is the derivative of Λ with respect to ρ.
Employing the second Cartan equation one derives the curvature two-forms R ab = R abcd ω c ∧ ω d , with 1 ≤ c < d ≤ 4, leading to the curvature components pointing out the special radius value ρ * = 2/B * , for which only the components R 1212 = B 2 * /Λ 4 = R 3434 are surviving and the Weyl tensor vanishes. Since the scalar curvature is zero, the Einstein tensor components are given by the Ricci tensor components, as In the pseudo-orthonormal frame whose dual bases is (8), it turns out that the potential (4), generating the magnetic induction along Oz, gets the familiar expression and the essential component of the Maxwell tensor reads Using the energy-momentum tensor components in the Einstein equations G ab = κ 0 T ab , one gets the following relation between the parameters B * and B 0 , with κ 0 = 8πG/c 4 .

Exactly solvable Klein-Gordon equation
In this section, we are going to construct the wave function of the charged bosons, considered as test particles evolving in the crust of a relativistic magnetar. The complex scalar field of mass µ, minimally coupled to gravity, is described by the SO(3, 1) × U (1) gauge-covariant Klein-Gordon equation which, in the pseudo-orthonormal frame with the dual bases (8), reads The above form suggests the variables separation which leads to the following differential equation for the un-known function φ, with Λ defined in (7).
This can be exactly integrated, its solution being expressed in terms of the Heun biconfluent function as [9,10] where the variable and the parameters are respectively given by and (20) Let us point out that the Heun biconfluent equation has one regular singularity at the origin and one irregular at ∞ and can be obtained, from the Heun general equation, by a process of successive confluences [10].
Regarding the asymptotic behavior of the function (18), solution to the equation (17), that has a singularity in ρ → 0, due to the exponential term, this is vanishing for large y−values. On the other hand, for a regular solution at the origin y = 0 (where HeunB(0) = 1), one has to choose the plus sign of α in (20).

The SO(3, 1) × U (1)−gauge covariant Dirac Equation
For relativistic fermions of mass M , coupled to the external magnetic field generated by (13), the Dirac equation has the SO(3, 1) × U (1)−gauge covariant expression γ a Ψ ;a + M Ψ = 0 , where ";" stands for the covariant derivative In view of the relations (10), the term expressing the Ricci spin-connection in the equation (21) reads where we have introduced the function The explicit form of the Dirac equation (21) being one may use the variables separation to derive the differential equation satisfied by the part depending on ρ, i.e.
With the following function substitution the above equation becomes where and we are going to use the Dirac representation for the γ i matrices, where σ µ denotes the usual Pauli matrices.
In the followings, we are assuming that the particle is not moving along the magnetic field direction, i.e. p z = 0, and the bi-spinor χ is of the form so that the equation (28) decouples in two equations for the (two-component) spinors ζ and η i.e.
Applying the usual procedure, one gets the following differential equations and which cannot be analytically solved. However, by imposing the condition B 2 * qB 0 and neglecting the powers of ρ larger than 3, the equations (33) and (34) get the simpler forms and The corresponding solutions, i.e.
are expressed in terms of Heun's biconfluent functions [8][9][10] of variables and parameters and therefore the components of the bi-spinor ψ in (27) are given by Using the expressions (40) in (25), one may compute the radial current density, meaning particles per unit time and per unit covariant 2-surface and the corresponding (radial) current, represented in the figure 1, as a function of .
One may notice that, for x 1, the current is suddenly increasing from zero to a maximum value, which depends on the ratio M/ω and on the magnetic field intensity. The case corresponding to massless fermions is significantly less complicated. Thus, the equation with the variables separation (25), leads to the following differential equation satisfied by the part depending on ρ, i.e.
Within the framework of the gauge-invariant geometry, based on the semidirect product of the local groups SO(3, 1) and U (1), the present paper is focusing on the Klein-Gordon and Dirac equations describing particles evolving in a background endowed with the Melvin's metric. By making use of the Cartan's formalism, we have derived the corresponding Einstein-Melvin equations leading to the essential relation between the model's parameters (14). As a remark, switching to the canonical bases, the third covariant induction component is given by the expression B z = |g|F 12 (c) = |g| e 1 a e 2 b F ab ≡ B 0 .
In case of bosons, the Klein-Gordon equation can be integrated exactly, its solution being given by the Heun biconfluent functions of parameters (20).
The equations (33) and (34) coming from the Dirac equation (24) have rather complicated expressions, containing several additional terms which were neglected in [6]. In the assumption B 2 * qB 0 , we have been able to find solutions for the equations (35) and (36), expressed in terms of Heun's biconfluent functions (38). The corresponding radial particle-current (41), represented in the figure 1, is starting from zero, at the origin ρ = 0 and exhibits a rather non-trivial behavior, characterized by a sudden growth to a maximum value. This one is followed by a local minimum which might signal, in the approximation we have used, the presence of a plateau. The seemingly far away unlimited increasing is a result of the violation of the approximation holding condition which demands x < 1.