Effects of the fermionic vacuum polarization in QED

TL;DR

This paper investigates unique effects of vacuum polarization in QED caused by field sources, including exact calculations of interaction energies and their influence on atomic and string configurations, highlighting phenomena absent in classical Maxwell theory.

Contribution

It provides exact Meijer-G function expressions for Uehling interaction energy and explores vacuum polarization effects on hydrogen atoms and Dirac strings, revealing novel quantum phenomena.

Findings

Exact Meijer-G function expression for Uehling energy

Vacuum polarization effects on hydrogen atom energy levels

Interaction between Dirac strings and quantum particles

Abstract

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.

Equations

$$\mathcal{L}=-\frac{1}{16\pi }{F}_{\mu \nu }{F}^{\mu \nu }-J_{\mu }{A}^{\mu }-\frac{1}{8\pi }({\partial }_{\mu }{A}^{\mu }{)}^2$$

$${\partial }^{\mu }{\partial }_{\mu }{A}_M^{\nu }=4\pi J^{\nu }$$

$${\partial }_{\mu }{\partial }^{\mu }{D}_M(x-y)=4\pi {\delta }^4(x-y).$$

$$D_M(x-y)=\int \frac{d^{4}p}{(2\pi {)}^4}{\tilde{D}}_M(p)e^{-ip(x-y)},$$

$${\tilde{D}}_M(p)=-\frac{4\pi }{p^2}$$

$$A_M^{\mu }(x)=\int d^{4}x{D}_M(x-y)J^{\mu }(y).$$

$$A_M^{\mu }(x)=\int \frac{d^{4}p}{(2\pi {)}^4}{{\tilde{A}}_M}^{\mu }(p)e^{-ipx},$$

$$J^{\mu }(y)=\int \frac{d^4{p}^{\prime }}{(2\pi {)}^4}{\tilde{J}}^{\mu }(p^{\prime })e^{-i{p}^{\prime }y},$$

$${{\tilde{A}}_M}^{\mu }(p)={\tilde{D}}_M(p){\tilde{J}}^{\mu }(p).$$

$${\tilde{J}}^{\mu }(p)=-\frac{p^2}{4\pi }{{\tilde{A}}_M}^{\mu }(p).$$

$$E_M=\int d^3\mathbf{x}{d}^{4}y\frac{1}{2}{J}_{\mu }(\mathbf{x})D_M(x-y)J^{\mu }(\mathbf{y}).$$

$${\tilde{D}}(p)={\tilde{D}_{M}}(p)\big{[}1+\Pi(p)\big{]}=-\frac{4\pi}{p^{2}}\big{[}1+\Pi(p)\big{]}\ .$$

$$\Pi (p)=-\frac{\alpha }{\pi }{\int }_0^{1}dv\frac{v^2\left(1-\frac{1}{3}{v}^2\right)}{v^2+\frac{4m^2}{p^2}-1},$$

$$D(x-y)=\int \frac{d^{4}p}{(2\pi {)}^4}\tilde{D}(p)e^{-ip(x-y)}.$$

$$A^{\mu }(x)=\int d^{4}xD(x-y)J^{\mu }(y).$$

$$A^{\mu }(x)$$

$$=$$

$$A^{\mu }(x)$$

$$=$$

$$E=\int d^3\mathbf{x}{d}^{4}y\frac{1}{2}{J}_{\mu }(\mathbf{x})D(x-y)J^{\mu }(\mathbf{y})$$

$$E=E_M+\frac{1}{2}\int d^3\mathbf{x}{d}^3\mathbf{y}\int \frac{d^3\mathbf{p}}{(2\pi {)}^3}{J}_{\mu }(\mathbf{x})J^{\mu }(\mathbf{y}){\tilde{D}}_M(p^0=0,\mathbf{p})\Pi (p^0=0,\mathbf{p})e^{i\mathbf{p}(\mathbf{x}-\mathbf{y})}.$$

$$J^{\mu }(\mathbf{x})=q_1{\eta }^{\mu 0}{\delta }^3(\mathbf{x}-{\mathbf{a}}_1)+q_2{\eta }^{\mu 0}{\delta }^3(\mathbf{x}-{\mathbf{a}}_2).$$

$$E=\frac{q_{1}q_{2}}{a}+4q_{1}q_{2}\alpha\int_{0}^{1}dvv^{2}\Bigl{(}1-\frac{1}{3}v^{2}\Bigr{)}\int\frac{d^{3}{\bf p}}{(2\pi)^{3}}\frac{e^{i{\bf p}\cdot{\bf a}}}{{\bf p}^{2}(1-v^{2})+4m^{2}}\ ,$$

$$\int \frac{d^3\mathbf{p}}{(2\pi {)}^3}\frac{e^{i\mathbf{p}\cdot \mathbf{a}}}{{\mathbf{p}}^2(1-v^2)+4m^2}=\frac{1}{(1-v^2{)}^{3/2}}\int \frac{d^3\mathbf{q}}{(2\pi {)}^3}\frac{e^{i\mathbf{q}\cdot \mathbf{a}/\sqrt{1-v^2}}}{{\mathbf{q}}^2+4m^2}.$$

$$\int\frac{d^{3}{\bf p}}{(2\pi)^{3}}\frac{e^{i{\bf p}\cdot{\bf a}}}{{\bf p}^{2}(1-v^{2})+4m^{2}}=\frac{1}{4\pi}\frac{1}{1-v^{2}}\frac{1}{a}\exp\bigg{(}-\frac{2ma}{\sqrt{1-v^{2}}}\biggr{)}$$

$$E$$

$$=$$

$${\int }_1^{\infty }du\frac{1}{\sqrt{u^2-1}}{e}^{-2mau}=K_0(2ma),$$

$${\int }_1^{\infty }du\frac{1}{\sqrt{u^2-1}}\frac{1}{u^2}{e}^{-2mau}=\frac{1}{2}(ma{)}^3\text{mbox}MeijerG([[],[0]],[[-1/2,-1,-3/2],[]],(ma{)}^2),$$

$${\int }_1^{\infty }du\frac{1}{\sqrt{u^2-1}}\frac{1}{u^4}{e}^{-2mau}=\frac{1}{2}(ma{)}^5\text{mbox}MeijerG([[],[0]],[[-1/2,-2,-5/2],[]],(ma{)}^2),$$

$$\displaystyle E=\frac{q_{1}q_{2}}{a}\Biggl{[}1+\frac{\alpha}{3\pi}\Biggr{(}2K_{0}(2ma)-\frac{1}{2}(ma)^{3}\mbox{MeijerG}([[],[0]],[[-1/2,-1,-3/2],[]],(ma)^{2})$$

$$\displaystyle-\frac{1}{2}(ma)^{5}\mbox{MeijerG}([[],[0]],[[-1/2,-2,-5/2],[]],(ma)^{2})\Biggr{)}\Biggr{]}\ .$$

$$\displaystyle q_{eff}=q_{1}\biggl{[}1+g(ma)\biggr{]}\ .$$

Full text

I Introduction

Since the establishment of the QED, the effects regarding the the vacuum polarization had been drawing attention mainly in situations with no counterpart in classical Electrodynamics. In this context, we can mention the Uehling potential Uehling , which is the electromagnetic potential associated with a single point-like stationary charge corrected in lowest order in the fine structure constant. The calculation of this potential is usually found in the literature perturbatively in the momentum space for a point-like source GreinerQED ; Zee ; Kaku ; Itzykson and was also generalized for a charge distribution with finite radius JPA271994 . Recently, the Uehling potential was calculated exactly in terms of Bessel Integral functions Frolov . In the context of non relativistic quantum mechanics, the Uehling potential can lead, for instance, to effects on hydrogen-like atoms JPB492016 . Atomic effects of QED not accounted by Uehling potential were also studied in the literature PRA481993 .

Another interesting scenario created by the vacuum polarization, with no counterpart in classical Electrodynamics, is related to the effects which emerge around solenoids and Dirac strings. In this context we can mention, for instance, the radiative corrections to the Aharonov-Bohm scattering PRD611999 , the interaction between two solenoids FoundPhys231993 , vacuum currents produced around a Dirac string PRD471993 ; IJMP61991 , Bremsstrahlung and pair production in the Aharonov-Bohm potential IJTP442005 and so on.

The vacuum polarization can be modified by the presence of an external field. The most common situations studied in this context are the effects produced by the presence of an external magnetic field in the Coulomb interaction and in the hydrogen-like atoms PR1011956 ; PRA381988 ; JPB271994 ; NPB5852000I ; NPB5852000II ; NPB6092001 . Modifications in the nuclear Coulomb field induced by a strong laser field were also considered in the literature PRA722005 .

The vacuum polarization was also studied in coordinate space in reference PRA892014 , 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 $\alpha$ (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 $\alpha$, 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

[TABLE]

where the last term is a gauge fixing one, $A^{\mu}$ is the electromagnetic field, $F^{\mu\nu}=\partial^{\mu}A^{\nu}-\partial^{\nu}A^{\mu}$ is the field strength and $J^{\mu}$ is the external source.

From the Lagrangian (1) one obtain the dynamical equation

[TABLE]

for which the corresponding propagator $D_{M}(x-y)$ satisfies the differential equation

[TABLE]

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^{\mu}$,

[TABLE]

where

[TABLE]

is the Fourier transform of the propagator $D_{M}(x-y)$

The solution for the field equation (2) is given by the integral

[TABLE]

Substituting the Fourier integrals for the field configuration and for the external source,

[TABLE]

and the integral (4) in Eq. (6) and using the fact that $\int d^{4}ye^{i(p-p^{\prime})y}=(2\pi)^{4}\delta^{4}(p-p^{\prime})$, we can show that

[TABLE]

With the aid of Eq’s (10) 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,

[TABLE]

For a steady external source, $J^{\mu}(x)=J^{\mu}({\bf x})$, one can show that the energy stored in the electromagnetic field is given by

[TABLE]

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 GreinerQED ; Itzykson , as follows

[TABLE]

where

[TABLE]

with $\alpha$ 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

[TABLE]

In the presence of an external source, the corrected field configuration is

[TABLE]

Substituting the second Eq. (7) and Eq. (15) in Eq. (16), we can show that

[TABLE]

where, in the second line, we used Eq. (13).

With the aid of Eq. (10), we can rewrite Eq. (17) in the form

[TABLE]

In the second line of Eq. (20), 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^{\mu}_{M}(x)$.

The energy stored in the electromagnetic field due to the presence of a static external source $J^{\mu}({\bf x})$, taking into account the corrections imposed by the fermionic vacuum bubbles, is given by

[TABLE]

Substituting Eq. (15) in (23), using expression (13), integrating out in $dy^{0}$ and in $dp^{0}$ and using the fact that $\int dy^{0}e^{ip^{0}y^{0}}=2\pi\delta(p^{0})$, we have

[TABLE]

Notice that Eq. (24) 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 ${\bf a}_{1}$ and ${\bf a}_{2}$, respectively,

[TABLE]

Substituting Eq’s (25), (14) and (5) in (24), 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

[TABLE]

where we defined the vector ${\bf a}={\bf a}_{1}-{\bf a}_{2}$ and its corresponding modulus $a=|{\bf a}|$.

The right and side of Eq. (26) is the same as the one found in the calculations of the Uehling potential GreinerQED . We shall calculate it exactly in this section.

The first term on the right hand side of Eq. (26) is the coulombian interaction between the charges. This well known result Zee ; BaroneHidalgo1 ; BaroneHidalgo2 ; HBB ; CBB ; BB 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 variables ${\bf q}={\bf p}\sqrt{1-v^{2}}$, as follows

[TABLE]

Now we perform the integral in $d^{3}{\bf q}$ with the results of reference BaroneHidalgo1 . The result is

[TABLE]

Substituting the result (28) in Eq. (26), performing the change of integration variable $u=1/\sqrt{1-v^{2}}$ and making some simple manipulations, we have

[TABLE]

All integrals in Eq. (29) can be performed exactly. The first one is given by

[TABLE]

where $K_{0}(2ma)$ is the K-Bessel function of second kind Arfken , and the other two ones are

[TABLE]

with MeijerG standing for the G-Meijer functions Gradshestein .

Substituting the result (32) in Eq. (29) we have finally the Uehling interaction

[TABLE]

Notice that expression (36) is exact, (up to order $\alpha$) valid for any distance $a$.

It is usual to interpret the result (36) as a coulombian interaction between a test charge $q_{2}$ and an effective one given by

[TABLE]

where we defined the function

[TABLE]

In the graphic (1) we can see a plot for the function (40). It gives the charge distribution induced in the fermionic vacuum by the presence of a point-like stationary electric charge. For $ma\to 0$ the function $g(ma)$ diverges as $\ln(ma)$ Itzykson ; GreinerQED . When $ma\to\infty$ the function $g(ma)$ goes to zero as $e^{-2ma}/(ma)^{3/2}$ Itzykson ; GreinerQED .

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 $\Phi$. Its corresponding source is given by BHN ; LB1 ; LB2 ; Fernanda ; Anderson

[TABLE]

If $\Phi>0$ we have the internal magnetic field pointing in the $\hat{z}$ direction. For $\Phi<0$, the internal magnetic field points in the opposite direction. The sub-index $\perp$ means we are taking just the components of a given vector perpendicular to the string. For instance, ${\bf p}_{\perp}=(p_{x},p_{y},0)$ is the momentum perpendicular to the string.

From Eq. (43) we can identify the Fourier transform of the Dirac string source

[TABLE]

Substituting the source (43) in expression (6), with the Maxwell propagator (4), we have the four-potential

[TABLE]

as expected for a Dirac string. The potential (45) produces a vanishing electromagnetic field outside the $z$ axis111Along the $z$ axis, the magnetic field in infinity.

The vacuum polarization effects in the vicinity of a Dirac string can be obtained by substituting, in the solution for field configuration (17), the definitions (5) and (14) and the source for the Dirac string (44). The result is

[TABLE]

where $A^{\mu}_{M(D)}(x)$ is given by (45) and we defined

[TABLE]

Integrating out expression (47) in $p^{0}$ and $p^{3}$ and noticing that just its $1$ and $2$ components are non-vanishing, we can write

[TABLE]

The integral above is calculated in reference BaroneHidalgo1

[TABLE]

where $K_{0}(x)$ stands for the Bessel function of second kind. So, acting with the derivatives, performing the change of integration variable

[TABLE]

and using the cylindrical coordinates, with $\rho=|{\bf x}_{\perp}|=\sqrt{x^{2}+y^{2}}$, we can write the vector potential (48) in the form

[TABLE]

with $\hat{\phi}$ standing for the unitary vector for the azimuthal coordinate.

The potential (53) 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 $\xi$ variable. The result is,

[TABLE]

where we defined the function

[TABLE]

The function (55) is always positive, as one can see in the graphic (2), so the magnetic field (54) points in the opposite direction in comparison with the magnetic flux of the string.

For small and long distances from the string, the magnetic field (54) reads

[TABLE]

III.1 The 2-D quantum rigid rotor

The vector potential (53) 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. (46), the total vector potential produced by the string is composed by the one obtained from Maxwell theory (45) added by the correction (53). Performing the integral in Eq. (53), we can write

[TABLE]

where we defined the function

[TABLE]

which is always positive, as one can see in the graphic (3).

It is well known in the literature GriffithsQM ; Sakurai 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 Sakurai . 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 $\Phi$. In this case, the Hamiltonian for the charged particle reads

[TABLE]

where $F(mb)$ is defined by (60).

The energy eigenfunctions of the hamiltonian (61) are given by

[TABLE]

where $n=0,\pm 1,\pm 2,\cdots$ is any integer and $A$ is a normalization constant. Up to order $\alpha$, the corresponding energy levels are

[TABLE]

The first term on the right hand side of (63) is the well known Aharonov-Bohm energy GriffithsQM and the second term is a correction due to the vacuum polarization.

For small and large values of $mb$, the energy (63) reads

[TABLE]

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. (54), as follows

[TABLE]

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 (67) for $md>>1$, as follows

[TABLE]

So, in this setup, the Dirac string produces a kind of Zeeman effect on the atom, once the field (70) 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 $-{\hat{z}}$ direction) for weak external field approximation to study the system, where the fine structure constant dominates the energy corrections Cohen . We shall take the hydrogen atom on its ground state ($n=1$, $\ell=0$, $m_{\ell}=0$, $j=s=1/2$). In this case, the degeneracy is broken as follows

[TABLE]

where $E_{0}$ is the non-perturbed ground state energy of the hydrogen atom, $\mu_{B}$ is the Bohr magneton and $m_{j}=\pm 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 (70) points along the $-{\hat{z}}$ direction.

As usual, in this case, the contribution of the nuclear magnetic moment (proton) is not relevant Cohen .

Substituting Eq. (70) in (73) we have

[TABLE]

The energy (74) 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

[TABLE]

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 (75) 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 (75) 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 (75) 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 ${\bf a}=(a^{1},a^{2},0)$, we can identify the distance between the strings by the modulus of ${\bf a}$, $\mid{\bf a}\mid=\sqrt{(a^{1})^{2}+(a^{2})^{2}}$. From Eq. (43), we can write the sources for the two strings, as follows

[TABLE]

Substituting the sources (78) in Eq. (24), 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^{\prime 0},\ dk^{3},\ dk^{\prime 3},\ dx^{3}$ and $dp^{3}$, identifying the length of a Dirac string $L=\int dy^{3}$ and performing some simple manipulations, we have the interaction energy between the strings per unit of length

[TABLE]

Now we define the differential operator

[TABLE]

use the fact that ${\bf p}_{\perp}^{2}e^{i{\bf p}_{\perp}\cdot{\bf a}}=-\nabla_{a}^{2}e^{i{\bf p}_{\perp}\cdot{\bf a}}$ and substitute the definitions (5) and (14) in Eq. (81), what leads to

[TABLE]

Using the result (51), acting with the operator $\nabla_{a}^{2}$ and performing the change of integration variables (52), the energy per unit length (85) reads

[TABLE]

The integral in Eq. (86) can be calculated. The final result for the energy is

[TABLE]

The gradient of the energy (87) with respect to ${\bf a}$ (with a minus signal) gives the force between the strings per unity of length

[TABLE]

which is repulsive when $\Phi_{1}$ and $\Phi_{2}$ have the same sign and attractive when the magnetic fluxes have opposite signs. In Eq. (88), ${\hat{a}}$ stands for the unit vector in the direction of ${\bf a}$ and the function inside brackets is positive along $0\leq ma\leq\infty$.

Figure (5) shows the graphic for the modulus of the force (88) divided by $\frac{8}{3\pi}\alpha\Phi_{1}\Phi_{2}m^{3}$ as a function of $ma$. For small and large values of $ma$, the behavior of the force (88) is

[TABLE]

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 BHN and theories with explicit Lorentz symmetry breaking LB1 ; LB2 .

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 $\alpha$. 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 $\alpha$), showing that it falls down vary fast when the distance between the strings increases.

Acknowledgments. F.A. Barone thanks to CNPq (Brazilian agency) under the grant 311514/2015-4 for financial support.