Effect of vacuum polarization of charged massive fermions in an Aharonov--Bohm field

The effect of vacuum polarization of charged massive fermions in an Aharonov-Bohm (AB) potential in 2+1 dimensions is investigated. The causal Green's function of the Dirac equation with the AB potential is represented via the regular and irregular solutions of the two-dimensional radial Dirac equation. It is shown that the vacuum current density contains the contribution from free filled states of the negative energy continuum as well as that from a bound unfilled state, which can emerge in the above background due to the interaction of the fermion spin magnetic moment with the AB magnetic field while the induced charge density contains only the contribution from the bound state. The expressions for the vacuum charge and induced current densities are obtained (recovered for massless fermions) for the graphene in the field of infinitesimally thin solenoid perpendicular to the plane of a sample. We also find the bound state energy as a function of magnetic flux, fermion spin and the radius of solenoid as well as discuss the role of the so-called self-adjoint extension parameter and determine it in terms of the physics of the problem.


I. INTRODUCTION
Quantum systems of relativistic fermions in external fields in 2+1 dimensions attract considerable interest related to the possibility of applying the results obtained for simplified models to study fundamental physical phenomena such as the AharonovBohm and quantum Hall effects [1,2], as well as hightemperature superconductivity [3]. When the external field configuration has the cylindrical symmetry, a natural assumption is that the relevant quantum mechanical system is invariant along the symmetry (z) axis and the system then becomes essentially two-dimensional in the xy plane. Solutions of the Dirac equation with the AB potential in 2+1 dimensions have been used to describe scattering of spin-polarized electrons by infinitesimally thin solenoid [4] as well as the behavior of fermions in the field of a cosmic string [5].
Important physical quantities are the vacuum charge and current densities induced by the background field. Vacuum polarization effects in the (2+1)-dimensional quantum electrodynamics with homogeneous magnetic field and with a nonzero fermion density were studied in [6,7].
New interest to two-dimensional fermion systems with the energy spectrum governed by Dirac Hamiltonian was revived in connection with problems of graphene (see, e.g., [8] and [9][10][11]). In graphene, the electrons near the Fermi surface can be described in terms of an effective Lorentz-invariant theory, with their kinetic energy determined by Dirac's dispersion law, and the dynamics of the electron at low energies is described by the Dirac equation in 2+1 dimensions for a zero-mass fermion [9]. It should be noted that while a description of electron states in graphene in [12][13][14] were based on the Dirac equation for massless fermions, work [15] has shown that the massive case can also be created.
At the same time, the effective fine structure constant in graphene is large, which leads to a new possibility for studying quantum electrodynamics in the strong-coupling regime. Charged impurity screening in graphene in terms of vacuum polarization were investigated in [14,[16][17][18][19]. The induced vacuum current in graphene in the field of a solenoid was investigated in [20] and a wonderful phenomenon was revealed: the induced current turns out to be a finite periodic function of the magnetic flux. Induced vacuum condensates such as the induced charge density and current (and other ones) for massless fermions in the background of a singular magnetic vortex in (2+1)-dimensional space-time were investigated in [21,22].
The the Dirac equation with the AB potential is governed by a (singular) Hamiltonian that requires the supplementary definition in order for it to be treated as a self-adjoint quantum-mechanical operator. Based on the form asymmetry self-adjoint extension method [23,24], the most relevant physical quantities, such as energy spectrum, wave functions and the self-adjoint extension parameter, were determined for some quantum systems with the AB potential by applying the so-called self-adjoint boundary conditions in [25][26][27][28]. It is useful to note that the self-adjoint extension regularization were applied for two-dimensional models in quantum field theory (see [29,30]). As a result some uncertainty appears in the prediction of physical quantities [29,30]. Essentially, a physical interpretation of self-adjoint extensions is a purely physical problem and each extension can (and must) be understood through an appropriate physical regularization [24].
In the present paper, we use the physical regularization method for constructing the fermion wave function and determining the energy in a bound state as a function of magnetic flux, fermion spin. We introduce a small parameter R, which, if we take in mind a real solenoid, is the finite radius of this solenoid (see, for instance, [4,31,32].
We shall adopt the units where c = = 1.

II. SOLUTIONS TO THE DIRAC EQUATION IN AN AHARONOV-BOHM FIELD
The space of particle quantum states in two spatial dimensions is the Hilbert space H = L 2 (R 2 ) of square-integrable functions Ψ(r), r = (x, y) with the scalar product The Dirac equation for a fermion in a given external field can be obtained just as in 3+1 dimensions. The Dirac γ µ -matrix algebra is known to be represented in terms of the two-dimensional Pauli matrices σ j and the parameter s = ±1 can be introduced to label two types of fermions in accordance with the signature of the two-dimensional Dirac matrices [33] and is applied to characterize two states of the fermion spin (spin "up" and "down") [4]. Then, the Dirac Hamiltonian for a fermion of the mass m and charge e = −e 0 < 0 in an Aharonov-Bohm potential A 0 = 0, A r = 0, A ϕ = B/r, r = x 2 + y 2 , ϕ = arctan(y/x) is where P µ = −i∂ µ −eA µ is the generalized fermion momentum operator (a three-vector). The Hamiltonian (3) should be defined as a self-adjoint operator in the Hilbert space of square-integrable two-spinors Ψ(r) with the scalar product (1). The total angular momentum J ≡ L z + sσ 3 /2, where L z ≡ −i∂/∂ϕ, commutes with H D .
In real physical space, because of the existence of finite magnetic flux inside solenoid Φ = 2πB the term including the spin parameter appears in the form of an additional delta-function interaction of spin with magnetic field of solenoid in the Dirac equation squared. The additional potential − seB δ(r) r will be taken into account by boundary conditions. Such kind of point interaction also appears in several Aharonov-Bohm-like problems [34][35][36][37].
Eigenfunctions of the Hamiltonian (3) are (see, [38]) where E is the fermion energy, l is an integer. The wave function Ψ(t, r) is an eigenfunction of the operator J with eigenvalue j = l + s/2. Taking into account the easily checked relations where µ ≡ e 0 B, we reduce the problem to that for the radial Hamiltonianȟ in the Hilbert space of doublets F (r) square-integrable on the half-line: whereȟ in the range r > 0. Eliminating, for instance, g l (r) we derive the differential equation for f l (r) and then the lower ("small") component of doublet is found from relation Then, for E 2 − m 2 > 0 the radial solutions can be written via the Bessel functions: Here A is a constant, ν = |l + µ|, p = √ E 2 − m 2 and the upper (lower) signs should be taken for l + µ > 0 (l + µ < 0).
Let us represent where n denotes the largest integer ≤ µ, i.e. the integer part of µ, and 0 < β < 1 is the fractional part of µ. Hence n = 0, 1, 2, ... for µ > 0 and n = −1, −2, −3, ... for µ < 0. Note that signs of e and B can be fixed and the potential is invariant under the changes e → −e, s → −s, it hence suffices to consider only the case e < 0, µ > 0. One can suppose that a bound state exists due to the interaction of the fermion spin magnetic moment with AB magnetic field, which must be attractive. In the case µ > 0, the potential is attractive for s = −1 and repulsive for s = 1, in the case µ < 0, it is attractive for s = 1 and repulsive for s = −1. Then, it is seen that when β = 0 the upper and lower components of (11) are integrable near r = 0 only for ν ± s > −1.
Linearly independent at ν = 0 solution U (r; E) is determined as follows: Here B is a constant. We shall need the irregular (i.e. integrable at r → ∞) solutions for E 2 − m 2 < 0. Such irregular solutions, nontrivial at ν = n/2, n = 1, 2, ... are the linear combination of F and U and also can be represented via the MacDonald functions: where C is a constant and λ = √ m 2 − E 2 . The irregular solutions are integrable at r → ∞.

III. A BOUND FERMION STATE ENERGY. PHYSICAL REGULARIZATION
A δ(x) potential is convenient to take into account artificially by means of the continuity relations. For this we replace the spin term (5) by (see, [4]) and take account of it by means of the continuity relations in R. The quantity Bδ(r − R)/R should not be considered as the real magnetic field inside a flux tube but only as a field model allowing us to take it into account by means of the continuity relations for solutions in two ranges. Although the functional structure of Eqs. (5) and (15) are quite different, as discussed in [32], we are free to use any form of potential provided that only the contribution of the form (5) is excluded. Now we can find the wave function and energy of bound state by means of solutions to the Dirac equation in ranges r < R and r > R and the potential (15) can be taken into account by means of the continuity relations. Obviously, for the model with the zero AB potential in the range r < R the radial solutions must satisfy (9) at µ = 0.They are written via the Bessel functions integrable near r = 0: where C is a constant, |l| ± s ≥ 0 and the upper (lower) signs should be taken for l > 0 (l < 0). The continuity relations can be written as The left-and right-hand sides of (17) are calculated using the asymptotic representation for the Bessel functions in the limit z ≪ 1: Here Γ(x) is the Euler gamma function of argument x.
The energy spectrum is determined by As a result, we obtain It follows from Eq. (20) that the bound state energy is determined by the poles of gamma function Γ(1 + |l| ± s) at l = 0: where we must put n = 0 and take the upper (lower) signs for s = −1, β > 0 (s = 1, β < 0). In result, we obtain the same equation for these two cases in the form: We define particle (antiparticle) bound states as states that tend to the upper (lower) continuous boundary m (−m) upon adiabatically slow switching of the external field. Only one particle (antiparticle) bound state with s = −1 (s = 1) emerges. It is seen that these equations keep for the case l + n = ∓1 (s = ±1). One can find from (22) that an adiabatic increase of the magnetic flux (β) from 0 to 1 lifts a particle (antiparticle) energy level E = m → E = −m (E = −m → E = m) and for |β| = 1/2 the particle and antiparticle energies are equal E = ± √ m 2 − R −2 , where the upper (lower) signs should be taken for particle (antiparticle). We believe that doublet (N is a normalization factor) represents right the particle radial wave function (with s = −1) in bound state. It is evident, the wave function (23) is singular at r = 0 but square integrable on the half-line [0, ∞) with the measure rdr.
The bound state energy for the considered system was derived in [28] by means of the self-adjoint extension method in the form: Here ξ = tan(θ/2) (2π ≥ θ ≥ 0) parameterizes the self-adjoint extensions the radial Dirac Hamiltonian, which are different for various θ except for two equivalent cases θ = 0, 2π (or ξ = ±∞). We can determine the self-adjoint extension parameter in terms of the physics of the problem, i.e. the parameter R. By comparing equations (22) and (22) we arrive, for example, for β > 1/2

IV. VACUUM CHARGE AND CURRENT DENSITIES
Now we consider the densities of vacuum electric charge and vacuum electric current due to to vacuum polarization. These quantities are determined by the three-vector j µ (r), which is expressed via the Green's function of the Dirac equation as follows where C is the path in the complex plane of E enclosing all the singularities along the real axis E depending upon the choice of the Fermi surface (we chose E F = −m). The singularities of G(r, r ′ ; E) can be simple poles associated with the discrete spectrum (in the range −m < E < m), and two cuts (−∞, −m] and [m, ∞) associated with the continuum spectrum in the ranges |E| ≥ m. As was shown in [39]), for the partial Green's function in a Coulomb field in 3+1 dimensions, the path C may be deformed to run along the singularities on the real E axis as follows: C = C − + C p + C + , where C − is the path along the negative real E axis (ReE < 0) from −∞ to 0 turning around at E = 0 with positive orientation, C p is a circle around the bound states' singularities with −m < E < 0, and C + is the path along the positive real E axis (ReE > 0) from ∞ to 0 but with negative orientation (i.e. clockwise path) turning around at E = 0. We note that due to the Furry theorem, the spatial component of induced vacuum current in an AB potential in 2+1 dimensions should be an odd function of µ. In the considered case this is j ϕ -component of induced vacuum current. It will be recalled that the j ϕ -component is the vector product of vectors j and n = r/r and has the only component [j × n] = j x n y − j y n x ; it is the so-called pseudoscalar. Thus, the j ϕ -component of induced current is determined as follows For the Dirac equation in a Coulomb field in 3+1 dimensions, the radial partial Green's function is given by [39] G l (r, where W(E) is the Wronskian and ψ R (r) and ψ I (r) are the regular and irregular solutions of the radial Dirac equation (Ȟ−E)ψ(r) = 0. One can show that the Green's function in our case can be represented via the regular and irregular solutions of the two-dimensional radial Dirac equation just as in 3+1 dimensions. It is convenient to apply doublet (11) as the regular solutions and (14) as the irregular ones. In such a way, we first construct the Green's function for the model with the nonzero (at r < R) AB potential by analogy with the Coulomb case in 3+1 dimensions. The (r-independent) Wronskian, defined by two doublets (11) and (14) as Here the upper (lower) signs should be taken for l + µ > 0 (l + µ < 0). For the induced vacuum charge and current in the AB potential, we obtain and In (30)  ; x = λ(β, R)r, where λ(β, R) is determined by (22) so x ∼ c(r/R), c ∼ 1.
Using recurrent relations for the Bessel functions and summing over s, we find and Integral over E from the first term of integrand (32) gives 0, because integrals taken along the paths C − and C + cancel in pairs. Now, it is possible to deform the paths C − and C + to the imaginary E axis: and Here, symbols b and v characterize the contributions from bound and free states, respectively. In (34) the upper (lower) signs should be taken for l + µ > 0 (l + µ < 0), I ν (z) = e −iπν/2 J ν (iz) is the modified Bessel function of first kind of argument z = √ m 2 + E 2 r. An essential detail is that the vacuum charge (j b 0 (r)) and current (j b ϕ (r)) are singular at the origin, localized near the point r = 0 and exponentially small at r ≫ R: Summation over l in (34) gives j v 0 (r) = 0. For massless fermions, this result was obtained and explained in [20] as follows: as far as j v 0 (r) should be the odd function of µ due to the Furry theorem, therefore, it must be pseudoscalar that would contradict to the parity conservation of the two-dimensional Dirac equation for massless particle. In the two-dimensional model with mass term, this term is not invariant with respect to the operations of time inversion and of spatial parity. Nevertheless, there appears the nonzero vacuum charge density j b 0 (r) due to the vacuum charge of bound (empty!) state; the vacuum charge spatial distribution is defined by the modulus squared of the fermion wave function in the bound state (with β > 1/2). The vacuum current density contain the vacuum currents j b ϕ (r) and j v ϕ (r) that are, respectively, due to the vacuum current of bound state and the vacuum current of free (filled!) states j v ϕ (r), which is spread out over an energy range of the negative energy continuum. It is worth to note that the vacuum charge density j v 0 (r) is induced by the homogeneous background magnetic field in the massive QED 2+1 [6,7], but this is not so in the QED 3+1 .
Using the easily checked representation and replacing variable x by y according to sinh x = sinh −1 y, we obtain At m = 0 (37) takes the form which coincides exactly with that for the induced vacuum current of massless fermions for the first time obtained in [20]. The integrals over y from any summand of (38) diverge, so some quantity δ ≪ 1 should be introduced as a lower limit of integration over y in order to we could change the order of summation and integration [20]. After that we first take integral over E by means of formula [39] As a result, we obtain Taking the sum over l by means of and then integral over y with using formula [39] we finally arrive at We emphasize that j b ϕ (r) = 0 at m = 0 because the normalization constant N ∼ √ m. A massless charged fermion can not be bound with an AB potential. As far as the main contribution to the integral over E at r ≫ 1/m is given by the region E ∼ 1/r, one can expect that Eq. (43) can be used for estimation of massive case if we replace r 2 by r 2 1 + (mr) 2 in the denominator of (43).
One can see that the induced current depends only on the fractional part β of µ and is finite periodical function of the magnetic flux µ. For the first time these results were obtained in [23].
It should be noted that the nonzero probability current arises under the effect of a constant uniform magnetic field on an electron bound by an attractive delta-function potential [40,41] The spatial distribution of probability current density resembles the spatial induced vacuum current j b ϕ (r) but the j b ϕ (r) current arises due to vacuum polarization, i.e. when bound state is not filled. We also note that vacuum polarization must manifest itself in such a way so as to modify (change) the external potential.
If we take in mind a real solenoid, then to learn the role of finite small radius R we consider again the semirealistic model with the nonzero AB potential at r < R. For such a model we can use the linear combination of the J ν (pr) and J −ν (pr) Bessel functions but not the Bessel (J |l|±s (pr)) and Neumann (N |l|±s (pr)) functions as the ψ R , ψ I solutions for r < R. For r > R the K ν , I ν modified Bessel functions apply. Then, continuity at r = R determines the linear combination of the J ν and J −ν solutions that joins the K ν function giving the irregular solution. Continuity at r = R also can determine the particular linear combination of the K ν and I ν that joins the regular (i.e., integrable near r = 0!). With such constructed solutions the Greens function will gain finite-size correction functions at r < R and r > R (see, [39]). It is wonderful that for the considered model we do not need to construct the regular solutions in the range r < R because we can apply regular solutions constructed by means of the self-adjoint adjoint method in [42]. Regular solutions (doublets) must satisfy the so-called self-adjoint boundary conditions [16] Physically, the self-adjoint boundary conditions show that the radial component of probability current density is equal to zero at the origin (the "origin" does not produce particles). For our problem the needed regular solution has the form where D is the normalization factor and, as above, 1 > β > 0, ξ = tan(θ/2). Note that the correct values of the self-adjoint parameter determine the behavior of the upper (lower) component of doublet (45) at the origin. Particularly, the case θ = 0 (θ = π) is equivalent to insisting that the upper (lower) component stays regular at the origin for any s = ±1 and, generally speaking, for µ > 0. If θ = 0, π both components of the doublet contain singular terms at the origin. We need to estimate the contribution of the correction functions to the induced vacuum current. To estimate the Green's correction function at r ≫ R we do not need to sum over l and s, it hence suffices to consider only the case, for example, s = −1. First, from simple system we obtain At r ≫ R the contribution of the Green's (singular) correction function to the integrand in Eq. (33) is determined with C 1 (ER). Estimating C 1 with using Eq. (18), one can obtain: C 1 (ER) ∼ (ER) 2 cos(θ/2), (ER) 2β sin(θ/2) The main contribution to the integral over E at r ≫ R is given by E ∼ 1/r so that the contribution of the Green's correction function to the induced current is suppressed by the factors (R/r) 2 cos(θ/2), (R/r) 2β sin(θ/2). These results are in agreement with estimations made in [20]. The contribution of the Green's (singular) correction function to the vacuum charge density j v 0 (r) may be nonzero and contain the same factors.

V. SUMMARY
We have investigated the effect of vacuum polarization for charged massive fermions in an AB potential in 2+1 dimensions using the causal Green's function of the Dirac equation with the AB potential represented via the regular and irregular solutions of the two-dimensional radial Dirac equation, which takes into account the fermion spin. It is shown that the vacuum current density contains the contribution from free filled states of the negative energy continuum as well as that from a bound (empty!) state, which can emerge in the above background due to the interaction of the fermion spin magnetic moment with the AB magnetic field while the nonzero vacuum charge density appears only due to the contribution from the bound state. We have derived (recovered for massless fermions) expressions for the vacuum charge and vacuum current densities for graphene in the field of infinitesimally thin solenoid perpendicular to the plane of a sample.
"A quantum-tunneling system using two-dimensional ionic structures in a linear Paul trap" was realized recently in [43], in which "the charged quantum-tunneling particles should be affected by the vector potential of a magnetic field throughout the entire process, even during quantum tunneling, and, thus, the AB effect should occur for tunneling particles. Authors "were successful in observing the AB effect of tunneling particles using this system". It was revealed a wonderful effect: "the tunneling rate of the structure periodically depends on the strength of the magnetic field, whose period is the same as the magnetic-flux quantum φ 0 . It seems that this result is in agreement with ours described, for example, by Eq. (43).