Effects of the fermionic vacuum polarization in QED

Some effects of vacuum polarization in QED due to the presence of field sources are investigated. We focus on effects with no counter-part in Maxwell electrodynamics. The the Uehling interaction energy between two stationary point-like charges is calculated exactly in terms of Meijer-G functions. Effects induced on a hydrogen atom by the vacuum polarization in the vicinity of a Dirac string are considered. We also calculate the interaction between two parallel Dirac strings and corrections to the energy levels of a quantum particle constrained to move on a ring circumventing a solenoid.

The vacuum polarization was also studied in coordinate space in reference [22], where the Green's function was calculated in higher orders beyond the Uehling term.
In this paper we study some effects produced by the vacuum polarization of the fermionic field. We focus on situations with no counterpart in classical electrodynamics, studying setups where sources for the electromagnetic field can interact via the vacuum polarization. In section (II) we start by studying the standard interaction between two point-like steady charges corrected by the vacuum polarization, in lowest order in the fine structure constant. We calculate exactly the Uehling interaction between them. Our result has two main advantages: it is an exact result (in lowest order in the fine structure constant) valid for any distance between the charges and the result is given by a simple expression written in terms of a K-Bessel functions and Meijer-G functions, what makes it easier to be plotted, once those functions are well known in the literature.
In section (III) we find out the field produced outside a Dirac string due to the vacuum polarization, in lowest order in the fine structure constant. We show that we have a magnetic field outside the string anti-parallel to the internal magnetic flux and we calculate the corrections in order α (the fine structure constant) in the energy levels of a quantum particle constrained to move on a ring (2-D quantum rigid rotor). In section (IV) we show that a hydrogen atom, in its ground state, interact with the string via a kind of Zeeman effect. This interaction falls down very quickly when the distance between the atom and the string increases. This force is attractive when the total angular momentum of the electron is parallel to the internal magnetic flux of the string, and repulsive in the opposite case.
In section (V) we show that it emerges an interaction between two Dirac strings due to the vacuum polarization of the fermionic field. We calculate this interacting force exactly, up to order α, for any distance between the strings when they are parallel to each other. Section (VI) is devoted to some comments and final remarks.

II. UEHLING INTERACTION
In this section we calculate exactly the Uehling interaction between two point-like steady charges.
The gauge sector of the Classical Electrodynamics can be described by the Lagrangian where the last term is a gauge fixing one, A µ is the electromagnetic field, F µν = ∂ µ A ν − ∂ ν A µ is the field strength and J µ is the external source.
From the Lagrangian (1) one obtain the dynamical equation for which the corresponding propagator D M (x − y) satisfies the differential equation The sub-index M in (3) means that we have the quantities calculated for the Maxwell theory.
The solution for (3) is given by the Fourier integral in the four momentum p µ , whereD is the Fourier transform of the propagator D M (x − y) The solution for the field equation (2) is given by the integral Substituting the Fourier integrals for the field configuration and for the external source, and the integral (4) in Eq. (6) and using the fact that With the aid of Eq's (8) and (5) we can write the Fourier transform of the external source as a function of the Fourier transform of the gauge field obtained from the Maxwell electrodynamics, For a steady external source, J µ (x) = J µ (x), one can show that the energy stored in the electromagnetic field is given by It is well known in the literature that the net QED effects of the fermionic vacuum bubbles can be taken into account by a correction in the gauge field propagator [2,5], as follows where with α standing for the fine structure constant and m, the mass of the electron.
The corrected propagator for the electromagnetic field is given by the Fourier integral In the presence of an external source, the corrected field configuration is Substituting the second Eq. (7) and Eq. (13) in Eq. (14), we can show that where, in the second line, we used Eq. (11).
With the aid of Eq. (8), we can rewrite Eq. (15) in the form In the second line of Eq. (16), the second term on the right hand side can be interpreted as a correction due to the vacuum polarization for the classical field configuration A µ M (x). The energy stored in the electromagnetic field due to the presence of a static external source J µ (x), taking into account the corrections imposed by the fermionic vacuum bubbles, is given by Substituting Eq. (13) in (17), using expression (11), integrating out in dy 0 and in dp 0 and using the fact that dy 0 e ip 0 y 0 = 2πδ(p 0 ), we have Notice that Eq. (18) is the energy obtained from the Maxwell electrodynamics, E M , with a correction term added.
Now, let us consider the external source produced by two stationary charges, q 1 and q 2 , placed at positions a 1 and a 2 , respectively, Substituting Eq's (19), (12) and (5) in (18), discarding the terms of self energies (the ones corresponding to the interactions of a given charge with itself) and performing some simple manipulations, we can write where we defined the vector a = a 1 − a 2 and its corresponding modulus a = |a|.
The right and side of Eq. (20) is the same as the one found in the calculations of the Uehling potential [2]. We shall calculate it exactly in this section.
The first term on the right hand side of Eq. (20) is the coulombian interaction between the charges. This well known result [3,[23][24][25][26][27] is obtained from E M (the interaction energy given by the Maxwell electrodynamics) and the propagator defined in (5) and (4). The integral in the second term can be calculated by changing the integration Now we perform the integral in d 3 q with the results of reference [23]. The result is Substituting the result (22) in Eq. (20), performing the change of integration variable u = 1/ √ 1 − v 2 and making some simple manipulations, we have All integrals in Eq. (23) can be performed exactly. The first one is given by where K 0 (2ma) is the K-Bessel function of second kind [28], and the other two ones are with MeijerG standing for the G-Meijer functions [29].
Substituting the result (24) in Eq. (23) we have finally the Uehling interaction Notice that expression (26) is exact, (up to order α) valid for any distance a.
It is usual to interpret the result (26) as a coulombian interaction between a test charge q 2 and an effective one given by where we defined the function In the graphic (1) we can see a plot for the function (28). It gives the charge distribution induced in the fermionic vacuum by the presence of a point-like stationary electric charge. For ma → 0 the function g(ma) diverges as ln(ma) [2,5]. When ma → ∞ the function g(ma) goes to zero as e −2ma /(ma) 3/2 [2,5].

III. FIELD OF A DIRAC STRING
In this section we discuss some effects of the vacuum polarization in the vicinity of a Dirac string. We choose a coordinate system where the Dirac string lies along the z-axis with internal magnetic flux Φ. Its corresponding source is given by [30][31][32][33][34] If Φ > 0 we have the internal magnetic field pointing in theẑ direction. For Φ < 0, the internal magnetic field points in the opposite direction. The sub-index ⊥ means we are taking just the components of a given vector perpendicular to the string. For instance, p ⊥ = (p x , p y , 0) is the momentum perpendicular to the string.
From Eq. (29) we can identify the Fourier transform of the Dirac string sourcẽ Substituting the source (29) in expression (6), with the Maxwell propagator (4), we have the four-potential as expected for a Dirac string. The potential (31) produces a vanishing electromagnetic field outside the z axis [38] The vacuum polarization effects in the vicinity of a Dirac string can be obtained by substituting, in the solution for field configuration (15), the definitions (5) and (12) and the source for the Dirac string (30). The result is where A µ M(D) (x) is given by (31) and we defined Integrating out expression (33) in p 0 and p 3 and noticing that just its 1 and 2 components are non-vanishing, we can write The integral above is calculated in reference [23] where K 0 (x) stands for the Bessel function of second kind. So, acting with the derivatives, performing the change of integration variable and using the cylindrical coordinates, with ρ = |x ⊥ | = x 2 + y 2 , we can write the vector potential (34) in the form withφ standing for the unitary vector for the azimuthal coordinate.
The potential (37) is static, has vanishing zero component and does not produce any electric field. Its rotational gives a magnetic field outside the solenoid. It is simpler, first, to calculate the relevant derivatives for the rotational and, after, perform the integration over the ξ variable. The result is, where we defined the function For small and long distances from the string, the magnetic field (38) reads A. The 2-D quantum rigid rotor The vector potential (37) induces a modification in the Aharonov-Bohm bound states. It can be seen by taking a very simple example of a 2-dimensional quantum rigid rotor composed by a non-relativistic quantum particle of mass M constrained to move along a circular ring surrounding the Dirac string. Let us take the ring on the plane z = 0 centered at the origin with radius b.
As stated in Eq. (32), the total vector potential produced by the string is composed by the one obtained from Maxwell theory (31) added by the correction (37). Performing the integral in Eq. (37), we can write where we defined the function which is always positive, as one can see in the graphic (3).
It is well known in the literature [35,36] that the energy levels of a two dimensional quantum rigid rotor are modified when it circumvents an infinite solenoid. In this case we have a very simplified version of the so called Aharonov-Bohm bound states [36]. In this section we consider the modification introduced by the vacuum polarization in the energy levels of a quantum rigid rotor. For this task, we take a quantum rigid rotor composed by a particle with mass M and electric charge q, restricted no move along a ring of radius b. We adopt a coordinate system where the ring lies on the plane z = 0, centered at the origin. We also consider a Dirac string placed along the z axis, with internal magnetic flux Φ. In this case, the Hamiltonian for the charged particle reads where F (mb) is defined by (42).
The energy eigenfunctions of the hamiltonian (43) are given by where n = 0, ±1, ±2, · · · is any integer and A is a normalization constant. Up to order α, the corresponding energy levels are The first term on the right hand side of (45) is the well known Aharonov-Bohm energy [35] and the second term is a correction due to the vacuum polarization.

IV. STRING-ATOM INTERACTION
The external magnetic field created by a Dirac string can lead to physical phenomena. Let us consider some of its effects produced on a single hydrogen atom. For this task, we take a coordinates system where the atom is placed at the origin and the Dirac string, parallel to the z axis, along the line (d, 0, z). In this setup, the Dirac string is placed at a distance d from the atom. We shall restrict to the situation where d is much higher in comparison with the atomic distances. The magnetic field produced in this case can be written by shifting the coordinate x in Eq. (38), as follows where, in the last line, we used the fact that the coordinates x and y are evaluated in the atomic distances and d >> x, y.
In a typical experiment d is a macroscopic distance in order of centimeters, what makes the values of the product md very large. So it is legitimate to approximate expression (47) for md >> 1, as follows So, in this setup, the Dirac string produces a kind of Zeeman effect on the atom, once the field (48) is, approximately, constant and uniform along the atom. In this regime, the field produced by the Dirac string is lower than the atomic magnetic field, so we shall use the results of the Zeeman effect (with external magnetic field pointing along the −ẑ direction) for weak external field approximation to study the system, where the fine structure constant dominates the energy corrections [37]. We shall take the hydrogen atom on its ground state (n = 1, ℓ = 0, m ℓ = 0, j = s = 1/2). In this case, the degeneracy is broken as follows where E 0 is the non-perturbed ground state energy of the hydrogen atom, µ B is the Bohr magneton and m j = ±1/2 is the azimuthal quantum number for the total angular momentum (for the ground state, m j = m s ). The minus signal for the Zeeman contribution is due to the fact that the external magnetic field (48) points along the −ẑ direction.
As usual, in this case, the contribution of the nuclear magnetic moment (proton) is not relevant [37].
Substituting Eq. (48) in (49) we have The energy (50) exhibits a dependence on the distance d between the Dirac string and the atom and produces a force on the atom (taking the string as fixed) given by where, in the last line, we discarded a term of order md.
When m j = 1/2, the projection of the total electronic angular momentum of the electron points in the same direction of the internal magnetic flux of the string and the force (51) becomes negative, what means that it exhibits an attractive nature. When the the total angular electronic momentum points in the opposite direction, with respect to the internal magnetic flux, the force (51) is positive, what means a repulsive behavior.
In a scattering experiment with a beam of unpolarized atoms propagating in the vicinity of a solenoid, we must have a bifurcation of the beam, according to the total angular momentum of each atom in the beam.
From the graphic (4) we can see the behavior of the magnitude of the force (51) as a function of md.

V. THE INTERACTION BETWEEN TWO STRINGS
In this section we study the interaction between two Dirac strings due to the vacuum polarization. We take two Dirac strings parallel to each other with the first one lying along the z axis and the second one lying along the line (a 1 , a 2 , z). Defining the vector a = (a 1 , a 2 , 0), we can identify the distance between the strings by the modulus of a, | a |= (a 1 ) 2 + (a 2 ) 2 . From Eq. (29), we can write the sources for the two strings, as follows Substituting the sources (52) in Eq. (18), discarding the self energies of each string, noticing that E M = 0 (there is no interaction energy between two Dirac strings in the Maxwell Electrodynamics), integrating out in dy 0 , dp 0 , dk 0 , dk ′0 , dk 3 , dk ′3 , dx 3 and dp 3 , identifying the length of a Dirac string L = dy 3 and performing some simple manipulations, we have the interaction energy between the strings per unit of length Now we define the differential operator use the fact that p 2 ⊥ e ip ⊥ ·a = −∇ 2 a e ip ⊥ ·a and substitute the definitions (5) and (12) in Eq. (53), what leads to Using the result (35), acting with the operator ∇ 2 a and performing the change of integration variables (36), the energy per unit length (55) reads The integral in Eq. (56) can be calculated. The final result for the energy is The gradient of the energy (57) with respect to a (with a minus signal) gives the force between the strings per unity of length which is repulsive when Φ 1 and Φ 2 have the same sign and attractive when the magnetic fluxes have opposite signs.
In Eq. (58),â stands for the unit vector in the direction of a and the function inside brackets is positive along 0 ≤ ma ≤ ∞.
As a last comment, we highlight that in some other gauge theories, even without radiative corrections, we can have interactions between Dirac strings and other sources for the gauge field, what is the case of Lee-Wick Electrodynamics [30] and theories with explicit Lorentz symmetry breaking [31,32].

VI. CONCLUSIONS AND FINAL REMARKS
In this paper we investigated some effects of the vacuum polarization, in lowest order in the fine structure constant, due to the presence of field sources for the electromagnetic field. All these effects were obtained from the the vacuum polarization tensor of the QED. We obtained an exact expression for the Uehling interaction (Uehling potential) between two point-like steady charges in terms of a K-Bessel function and MeijerG functions. Our results are compatible with that ones obtained in the literature, approximately, for long and small distances between the charges.
We also investigated some phenomena produced outside a Dirac string. One of these effects are the modifications induced in the energy levels of a quantum rigid rotor which circumvents a Dirac string. We have calculated these modifications exactly (for any radius of the quantum rigid rotor) up to order α. We have shown that a hydrogen atom interact with a Dirac string via a kind of Zeeman effect. The nature (attractive or repulsive) of this interaction depends on the orientation of the total magnetic moment of the electron with respect to the internal magnetic flux of the string. We restricted to the case where the atom is in its ground state and far away from the string. When the total magnetic moment of the electron is parallel to the internal magnetic flux of the string, the interaction is attractive, on the contrary, it is repulsive. As expected, this interaction is very small and falls down very fast when the distance between the string and the atom increases.
We have also investigated the interaction which emerges between two Dirac strings due to the vacuum polarization.
We restricted to the case where the strings are parallel or anti-parallel to each other. We showed that the strings attract each other when their internal magnetic fluxes are anti-parallel. This interaction is repulsive when the internal magnetic fluxes are parallel to each other. We have computed this interacting force exactly (up to order α), showing that it falls down vary fast when the distance between the strings increases.