Arxiv:1509.02594V2 [Hep-Th] 30 Dec 2015 Hoyfo Lsia on Fve,A Swl Nw.Nevertheless Discu Supers Been Known

Home , Uehling potential

arXiv:1509.02594v2 [hep-th] 30 Dec 2015 hoyfo lsia on fve,a swl nw.Nevertheless discu supers been known. have in well electrodynamics contexts is nonlinear different of as many types in view, other arise theory, of actions point effective classical BI-type als a the Electrodynamics of th from (BI) dependence Born-Infeld Whereas theory space-time case, the Euler-Heisenbeg strengths. in the nonlinear field in (but the fields the of t in powers linear t and in (in correction corrections calculation Euler-Kockel Serber nonlinear wa the and as extends Incidentally, Uehling result Lagrangian. Euler-Heisenberg the Euler-Heisenberg polariza are the vacuum from contributions requir derived from of clearly contributions This class the a energy. incorporate interaction fac to the In Lagrangian of particles. form charged heavy structural between the energy interaction the in eedn sn h ag-nain u ahdpnetvariable path-dependent corrections but polarization gauge-invariant vacuum the using when done case were o the physic choice to holes the black electrodynamics from to independent extended is successfully wh physics been rule the has multiplication approach that new new shown a [29–3 later through q was product) implemented (or It (Moyal is theory which product field [35–37], quantum start noncommutative length) theorie a formulate field to using gravitat quantum way by novel of context studied theory this quantum been Within a has inten include not program. of does gravity focus SM quantum the the become that have fact the length fundamental and scenarios cnro uha onIfl hoy[1,mlihre atce [12 the particles is physics. millicharged this scattering illustrates photon-photon [11], to that theory contributions [6–10 example Born-Infeld significant confirmed An as been yet such [1–5]. not scenarios QED has of progresses days remarkable early despite interaction the nonlinear in to leading as and pairs electron-positron virtual ∗ lcrncades [email protected] address: Electronic nti otx,w lopitotta xesoso h Standard the of extensions that out point also we context, this In nprdb hs bevtos h ups fti ae st ex to is paper this of purpose the observations, these by Inspired h hnmnno aumplrzto nQatmElectrodynamics Quantum in polarization vacuum of phenomenon The neetnl,i hudb ealdhr httepyia ffc of effect physical the that here recalled be should it Interestingly, o asv ue-esne-ieeetoyais(Which varia (1 electrodynamics path-dependent Euler-Heisenbeg-like gauge-invariant massive the for fo of energy sh framework ph interaction analysis the the birefringence to Our vacuum modifications lowest-order the limit. the displays strong-field Electrodynamics the type of approximation the h edn oaihi re) ipa ogrne(1 long-range a display order), ap logarithmic the leading the in Electrodynamics Euler-Heisenberg-like However, ASnmes 47.e 26.n 13.40.Gp 12.60.Cn, 14.70.-e, numbers: PACS interactio the versions, non-commutative their for Again, /r ecnie ohmsieElrHiebr-ieadEuler- and Euler-Heisenberg-like massive both consider We 3 ye orcin pr rmaln-ag (1 long-range a from apart correction, type) - 1 eatet eFıiaadCnr Cient´ıfico-Tecnológic Centro F´ısica and de Departmento nvria ćiaFdrc at aá Valpara´ıso, Mar´ıa, Santa Técnica Federico Universidad eak nnniereetoyaisII electrodynamics nonlinear on Remarks Dtd eebr3,2015) 31, December (Dated: .INTRODUCTION I. arcoGaete Patricio /r /r 5 ye orcint h olm potential. Coulomb the to correction type) - .Aogtesm ie eas ealta alternative that recall also we line, same the Along ]. 5 1, ye orcint h olm potential. Coulomb the to correction type) - ,ti hscleetcagsbt h teghand strength the both changes effect physical this t, sdi h ieaue[18–23]. literature the in ssed omls,weeteitrcinptnilenergy potential interaction the where formalism, s nryi lrvoe finite. ultraviolet is energy n steadto fcreto em nteMaxwell the in terms correction of addition the es ecntn akrudfil ii) hc contains which limit), field background constant he ewe lcrmgei ed,rmi sexciting as remain fields, electromagnetic between s otissmlrnniercretost Maxwell to corrections nonlinear similar contains o anKol eoti e long-range new a obtain we mann-Kroll) xlie n[] ti fitrs ontc htthe that notice to interest of is it [5], in explained s ohcasso lcrdnmc,within Electrodynamics, of classes both r ersac ciiy[42] hshsisoii in origin its has This [24–28]. activity research se ∗ [40]. s rxmto ftesrn-edlmt(to limit strong-field the of proximation ] eto hudb ae tti on,t a to point, this at made, be should Mention 4]. o,a ocruvn iclistertcli the in theoretical difficulties circumvent to as ion, rn hoy[6 7.I diint Born-Infeld to addition In 17]. [16, theory tring corrections contains result Serber and Uehling e atmfil hoyi h rsneo minimal a of presence the in theory field uantum oe S)sc sLrnzivrac violating invariance Lorentz as such (SM) Model inpoes w motn xmlso such of examples important Two process. tion lsfraim neetnl enough, Interestingly formalism. bles rainlk atce 1–5 a aemore have may [13–15] particles axion-like or ] r ae noacut h rcdn studies preceding The account. into taken are edorpeiu tde 1,1]o nonlinear on 19] [18, studies previous our tend akrudfils.W ute eto ht as that, mention further We fields). background nmnn fewrs ecalculate we Afterwards, enomenon. aumplrzto per samodification a as appears polarization vacuum h yeo rdc 3] oercnl,this recently, More [39]. product of type the f w htmsieEuler-Heisenberg- massive that ows eWcmn-rl orcin hc a be can which correction, Wichmann-Kroll he loignncmuigpsto operators position non-commuting allowing s esneglk lcrdnmc in Electrodynamics Heisenberg-like IEetoyaisi itnuse,since distinguished, is Electrodynamics BI , cteigo htn ypoos which photons, by photons of scattering QD,aiigfo h oaiainof polarization the from arising (QED), c skona oo trpout[38]. star-product Voros as known is ich eValpara´ıso, de o Chile 2 between two static charges is determined by the geometrical condition of gauge invariance. One important advantage of this approach is that it provides a physically-based alternative to the usual Wilson loop approach. Accordingly, we shall work out the static potential for electrodynamics which include, apart from the Maxwell Lagrangian, additional terms corresponding to the Uehling, massive Euler-Heisenberg-like, Euler-Heisenberg Electrodynamics in the approximation of the strong-field limit (to the leading logarithmic order) and for their non-commutative versions. Our results show a 1 long-range L5-type correction to the Coulomb potential for both massive Euler-Heisenberg-like and Euler-Heisenberg Electrodynamics in the approximation of the strong-field limit (to the leading logarithmic order). Interestingly enough, 1 for massive Euler-Heisenbeg-like electrodynamics (Whichmann-Kroll) we obtain a new long-range L3 correction to the interaction energy. Nevertheless, for their non-commutative versions, the static potential becomes ultraviolet finite. The organization of the paper is as follows: In Section II, we reexamine Uehling Electrodynamics in order to establish a framework for the computation of the static potential. In Section III we consider Euler-Heisenberg-like (with a mass term) and show that it yields birefringence, compute the interaction energy for a fermion-antifermion pair and its version in the presence of a minimal length. In Section III, we repeat our analysis for Euler-Heisenberg Electrodynamics in the approximation of the strong-field limit. Finally, in Section IV, we cast our Final Remarks. In our conventions the signature of the metric is (+1, 1, 1, 1). − − −

II. BRIEF REVIEW ON THE UEHLING POTENTIAL

As already expressed, we now reexamine the interaction energy for Maxwell theory with an additional term corresponding to the Uehling correction (Uehling Electrodynamics). This would not only provide the setup theoretical for our subsequent work, but also fix the notation. To do that we will calculate the expectation value of the energy operator H in the physical state Φ , which we will denote by H Φ. We start off our analysis by considering the effective Lagrangian density [41]: | i h i 1 α = F 1 ∆ F µν , (1) L −4 µν − 3π M where 2 ∞ 1 2m 4m2 (m, x)= dτ 1+ 1 , (2) M 2 τ (τ + ∆) τ − τ Z4m r µ with ∆ ∂µ∂ . It should be noted that contains the effect of vacuum polarization to first order in the fine 2 ≡ e M structure constant, α = ~c , and, m, is the electron mass. To obtain the corresponding Hamiltonian, we will carry out the quantization of the theory using the Hamiltonian formalism. In this way, the canonical momenta are found to be Πµ = 1 α ∆ F µ0, and one immediately identifies − 3π M the primary constraint Π0 = 0, whereas the momenta are Πi = 1 α ∆ Ei (with E = F ). Accordingly, the − 3π M i i0 canonical Hamiltonian HC is 1 1 1 H = d3x Πi∂ A Π 1 α ∆ − Πi + F 1 α ∆ F ij . (3) C i 0 − 2 i − 3π M 4 ij − 3π M Z 0 i As usual, requiring the primary constraint Π to be stationary, leads to the secondary constraint Γ1 ∂iΠ = 0. It is straightforward to check that there are no further constraints in the theory. Hence, the extended≡ Hamiltonian 3 that generates translations in time then reads H = HC + d x (u0(x)Π0(x)+ u1(x)Γ1(x)), where uo(x) and u1(x) are arbitrary Lagrange multipliers. With the aid of equation (3) we find that A˙ (x) = [A (x) ,H]= u (x), which is R 0 0 0 an arbitrary function. Since Π0 = 0 always, neither A0 nor Π0 are of interest in describing the system and may be discarded from the theory. Therefore, the Hamiltonian becomes

1 1 1 H = d3x w(x)∂ Πi Π 1 α ∆ − Πi + F 1 α ∆ F ij , (4) i − 2 i − 3π M 4 ij − 3π M Z where w(x)= u1(x) A0(x). According to the usual− procedure, we impose a gauge condition such that the full set of constraints becomes second class. A convenient choice is 1 Γ (x) dzν A (z) dλxiA (λx)=0. (5) 2 ≡ ν ≡ i CZξx Z0 3 where λ (0 λ 1) is the parameter describing the spacelike straight path zi = ξi + λ (x ξ)i, and ξ is a ﬁxed point (reference point).≤ ≤ There is no essential loss of generality if we restrict our considerations to−ξi = 0. The Dirac brackets can now be determined and we simply note the only nontrivial Dirac bracket involving the canonical variables, that is,

1 A (x) , Πj (y) ∗ = δj δ(3) (x y) ∂x dλxiδ(3) (λx y) . (6) i i − − i − Z 0 In passing we also recall that the transition to a quantum theory is made by the replacement of the Dirac brackets by the operator commutation relations according to A, B ∗ ( i/~)[A, B]. With the foregoing information, we can now proceed{ to} obtain→ − the interaction energy. As already mentioned, in order to accomplish this purpose, we will calculate the expectation value of the energy operator H in the physical state Φ , where the physical states Φ are gauge-invariant ones. The physical state can be written as | i | i y iq i Φ Ψ(¯ y)Ψ(y′) = ψ¯(y)exp( dz A (z))ψ(y′) 0 , (7) | i≡ ~ i | i Zy′ where 0 is the physical vacuum state and the line integral appearing in the above expression is along a space-like path starting| i at y and ending at y, on a ﬁxed time slice. Making use of the′ above Hamiltonian structure [18], we ﬁnd that

y′ (3) Π (x) Ψ (y) Ψ (y′) = Ψ (y) Ψ (y′) Π (x) 0 + q dz δ (z x) Φ . (8) i i | i i − | i Zy

With the aid of equations (8) and (4), the lowest-order modiﬁcation in α of the interaction energy takes the form H = H + V + V , (9) h iΦ h i0 1 2 where H = 0 H 0 . The V , V terms are given by h i0 h | | i 1 2 2 y y q 3 i (3) i (3) V = d x dz δ (x z) dz′ δ (x z′) , (10) 1 2 − − Z Zy′ Zy′ and

2 y′ y′ q α ∞ ρ (τ) 3 (3) 2 1 i (3) V2 = dτ d x dzi′δ (z′ x) x dz δ (z x) , (11) 2 3π 2 τ − ∇ τ 2 − Z4m Z Zy − ∇x Zy

2 2 where ρ (τ)= 1+ 2m 1 4m . τ − τ We note that the term q (10) may look peculiar, but it is nothing but the familiar Coulomb interaction plus a self- √τ z z 1 e− | − ′| energy term [42]. Now making use of the Green function, G (z, z′)= , the term (11) can be rewritten in 4π z z′ the form | − |

2 y′ y′ 2 √τ y y′ q α ∞ ρ (τ) 2 i α q ∞ ρ (τ) e− | − | V2 = dτ dzi′ z dz G (z, z′)= dτ . (12) 2 3π 2 τ ∇ ′ −3π 2 2 τ y y Z4m Zy Zy Z4m | − ′| Since the second and third term on the right-hand side of Eq. (9) are clearly dependent on the distance between the external static fields, the potential for two opposite charges located at y and y reads ′ 2 q 1 α ∞ ρ (τ) √τL V = 1+ dτ e− , (13) −4π L 3π 2 τ Z4m where L = y y′ . Accordingly, one recovers the known Uehling potential, which finds here an entirely different derivation. | − | Before we proceed further, we wish to show that this result can be written alternatively in a more explicit form. Making use of [43]

tz ∞ e 1 1 χ (z)= dt − 1+ 1 , (14) n tn 2t2 − t2 Z1 r 4 we then get

q2 1 2α V = 1+ χ (2mL) . (15) −4π L 3π 1

By the transformation, t = cosh u [44], the functions χn can be reduced to the form [43]: 1 1 χn (z)= Kin 1 (z) Kin+1 (z) Kin+3 (z) , (16) − − 2 − 2 where the functions Ki denote Bessel function integrals. Hence we see that the interaction energy (with m = 1) becomes

2 2 3 3 q 1 2α L 5L 2L 3L 2L ∞ V = 1+ 1+ K (2L) + K (2L)+ + dtK (t) , (17) −4π L 3π 3 0 − 3 3 1 2 3 0 Z2L where K0(z) and K1(z) are modiﬁed Bessel functions. Finally, with the aid of asymptotic forms for Bessel functions, it is a simple matter to ﬁnd expressions for V for large and small L. Before concluding this subsection we discuss an alternative way of stating our previous result (13), which displays certain distinctive features of our methodology. We start by considering [42, 45]

V q ( (0) (L)) , (18) ≡ A0 − A0 where the physical scalar potential is given by

1 (t, r)= dλriE (t, λr). (19) A0 i Z0 This follows from the vector gauge-invariant ﬁeld expression

x (x) A (x)+ ∂ dzµA (z) , (20) Aµ ≡ µ µ − µ Zξ where the line integral is along a space-like path from the point ξ to x, on a fixed slice time. It is also important to observe that the gauge-invariant variables (19) commute with the sole first constraint (Gauss law), showing in this way that these fields are physical variables. Inasmuch as we are interested in estimating the lowest-order correction 1 i α − i to the Coulomb energy, we will retain only the leading term in expression E = 1 3π ∆ Π . Making use of this last expression, equation (19) gives − M 1 0 1 0 i J α ∞ dτ i J 0(t, r)= dλr ∂i + ρ (τ) dλr ∂i , (21) A − 2 3π 2 τ − 2 τ Z0 ∇ λr Z4m Z0 ∇ − λr i 0 to get the last line we used Gauss law for the present theory, that is, ∂iΠ = J (where we have included the external current J 0 to represent the presence of two opposite charges). Accordingly, for J 0(t, r) = qδ(3) (r), the potential for a pair of static point-like opposite charges located at 0 and L, is given by

2 q 1 α ∞ ρ (τ) √τL V = 1+ dτ e− , (22) −4π L 3π 2 τ Z4m after substracting a self-energy term.

III. EULER-HEISENBERG-LIKE MODEL

Proceeding in the same way as we did in the foregoing section, we shall now consider the interaction energy for Euler-Heisenberg-like electrodynamics. Nevertheless, in order to put our discussion into context it is useful to describe very brieﬂy the model under consideration. In such a case the Lagrangian density reads:

β2 1 1 p = 1 1+ 2 , (23) L 2 − β2 F − β2γ2 G 5

1 µν 1 ˜µν where have included two parameters β and γ. As usual, = 4 Fµν F , = 4 Fµν F , Fµν = ∂µAν ∂ν Aµ and ˜µν 1 µνρλ F G − F = 2 ε Fρλ. Let us also mention here that in our previous paper [19] we have studied the domain 0

µν ∂µF˜ =0, (25) where 2 Γ=1+ F G . (26) β2 − β2γ2 Also, it is straightforward to see that Gauss law becomes, D =0, (27) ∇ · where D is given by

E2 B2 (E B)2 2 D = 1 − · E + (E B) B . (28) − 2β2 − β2γ2 γ2 · " # 0 (3) Q e Again, from equation (27), for J (t, r) = eδ (r), we ﬁnd D = r2 rˆ, where Q = 4π . This then implies that for a point-like charge, e, at the origin, the expression Q E2 = 1 E , (29) r2 − 2β2 | | tells us that, for r 0, the electrostatic ﬁeld becomes singular at r = 0, in contrast to the 0

With the aid from equations (30) and (31) we find the electric permitivity, εij , and the inverse magnetic permeability, 1 µ− ij , tensors of the vacuum, that is, 2BiBj 1 2EiEj ε =Γ δ + , µ− =Γ δ , (34) ij ij γ2 ij ij − γ2 6 with Di = εij Ej and Bi = µij Hj . In accordance with our previous procedure [18, 19], we can now linearize the above equations. To do this, it is advantageous to introduce a weak electromagnetic wave (Ep, Bp) propagating in the presence of a strong constant external field (E0, B0). On these assumptions, we readily find that, for the case of a purely magnetic field (E0 = 0), the vectors D and H become B2 2 D = 1+ 0 E + (E B ) B , (35) 2β2 p γ2 p · 0 0 and

B2 1 H 0 B B B B = 1+ 2 p + 2 ( p 0) 0 , (36) 2β  2 B0 ·  β 1+ 2β2   where we have keep only linear terms in Ep, Bp. As before, we consider the z axis as the direction of the external magnetic ﬁeld (B0 = B0e3) and assuming that the light wave moves along the x axis, the decomposition into a plane wave for the ﬁelds Ep and Bp can be written as

i(wt k x) i(wt k x) Ep (x,t)= Ee− − · , Bp (x,t)= Be− − · . (37) In this case, it clearly follows that k2 ε µ E =0, (38) w2 − 22 33 2 and k2 ε µ E =0. (39) w2 − 33 22 3 As a consequence, we have two diﬀerent situations: First, if E B0 (perpendicular polarization), from (39) E3 = 0, 2 k ⊥ and from (38) we get w2 = ε22µ33. This then means that the dispersion relation of the photon takes the form

2 B0 1+ 2β2 n = 2 . (40) ⊥ v 3B0 u1+ 2 u 2β u t 2 k Second, if E B (parallel polarization), from (38) E = 0, and from (39) we get 2 = ε µ . This leads to || 0 2 w 33 22 2B2 n = 1+ 0 . (41) k s γ2 Thus we verify that in the case of a generalized Euler-Heisenberg electrodynamics the phenomenon of birefringence is present. We now pass to the calculation of the interaction energy between static point-like sources for a massive Whichman- Kroll-like model, our analysis follows closely that of references [18, 19]. The corresponding theory is governed by the Lagrangian density 1 1 m2 = F 2 + (F F µν )2 + A Aµ. (42) L −4 µν 32β2 µν 2 µ Next, in order to handle the second term on the right hand in (42), we introduce an auxiliary ﬁeld ξ such that its equation of motion gives back the original theory. This allows us to write the Lagrangian density as 1 ξ 1 m2 = F F µν + F F µν ξ2 + A Aµ. (43) L −4 µν 32β2 µν − 128β2 2 µ

ξ With the redeﬁnition η =1 2 , equation (43) becomes − 8β 1 1 m2 = ηF F µν (1 η)2 + A Aµ. (44) L −4 µν − 2 − 2 µ 7

Before we proceed to work out explicitly the interaction energy, we shall ﬁrst restore the gauge invariance in equation (44). Following an earlier procedure, we readily verify that the canonical momenta read Πµ = ηF 0µ, which results in the usual primary constraint Π0 =0, and Πi = ηF i0. In this way one obtains −

1 η m2 1 H = d3x Πi∂ A + Π2 + B2 A Aµ + (1 η)2 . (45) C i 0 2η 2 − 2 µ 2 − Z i 2 0 The consistency condition, Π˙ = 0, leads to the constraint Γ ∂iΠ + m A = 0. As a result, both constraints are second class. To convert the second class system into first class≡ we shall adopt the procedure described previously. Thus, we enlarge the original phase space by introducing a canonical pair of fields θ and Πθ. It follows, therefore, that 2 a new set of first class constraints can be defined in this extended space: Λ1 Π0 + m θ =0 and Λ2 Γ + Πθ = 0. This then shows that the new constraints are first class and, therefore, restore≡ the gauge symmetry. From≡ this, the new effective Lagrangian density, after integrating out the θ fields, becomes

1 m2 1 = F η + F µν (1 η)2. (46) L −4 µν ∆ − 2 − 2 m Now, writting σ = η + ∆ , the expression (46) can be brought to the form

2 1 k m2 = F σF µν 1 σ + , (47) L −4 µν − 128 − ∆ where k = 64β2. We are now ready to compute the interaction energy. In this case, the canonical momenta are Πµ = σF 0µ, with the usual primary constraint Π0 =0, and Πi = σF i0. Hence the canonical Hamiltonian is expressed as −

1 σ k m2 2 H = d3x Πi∂ A + Π2 + B2 + 1 σ + . (48) C i 0 2σ 2 128 − ∆ Z ( ) 0 i Time conservation of the primary constraint Π yields the secondary constraint Γ1 ∂iΠ = 0. Similarly for the σ constraint yields no further constraints and just determines the ﬁeld σ. In this case,≡ at leading order in β, the ﬁeldPσ is given by

m2 B2 3 1 σ = 1+ 1 Π2 , (49) 2  2 3  ∆ − 2β − 2β m2 B2 1+ ∆ 2β2  −    which will be used to eliminate σ. As before, the corresponding total (ﬁrst-class) Hamiltonian that generates the time 3 evolution of the dynamical variables is H = HC + d x (u0(x)Π0(x)+ u1(x)Γ1(x)), where uo(x) and u1(x) are the Lagrange multiplier utilized to implement the constraints. In the same way as was done in the previous subsection,R the expectation value of the energy operator H in the physical state Φ becomes | i 1 1 m2 − 15 15m2 1 H = Φ d3x Πi 1+ Πi + Π4 Π2 Π2 Φ , (50) h iΦ h | 2 ∆ 8β2 − 2β2 ∆ | i Z ( ) in this last line we have considered only quadratic terms in m2. In such a case, by employing (50), the lowest-order modiﬁcation in β2 and m2 of the interaction energy takes the form

H = H + V + V + V , (51) h iΦ h i0 1 2 3 where H = 0 H 0 . The V , V and V terms are given by h i0 h | | i 1 2 3 2 y 2 1 y q 3 i (3) m − i (3) V = d x dz′ δ (x z′) 1 dz δ (x z) , (52) 1 2 − − 2 − Z Zy′ ∇ Zy′ 8

4 y y y y 15q 3 (3) i (3) k (3) (3) V = d x dz δ (x z) dz′ δ (x z′) du δ (x u) dv δ (x v) , (53) 2 − 8β2 i − − − k − Z Zy′ Zy′ Zy′ Zy′ and 2 4 y y y y 15m q 3 (3) i (3) k (3) (3) V = d x dz δ (x z) dz′ δ (x z′) du δ (x u) dv δ (x v) . (54) 3 2β2 i − − − k − Z Zy′ Zy′ Zy′ Zy′ Finally, with the aid of the expressions (52), (53) and (54), the potential for a pair of static point-like opposite charges located at 0 and L, is given by q2 e mL q4 3 1 1 V = − + 5m2 . (55) −4π L 16πβ2 8π L2 − L3 observe that when m = 0, the profile (55) reduces to the known Whichmann-Kroll interaction energy. On the other hand, for m = 0, it should be noted the key role played by the mass term in transforming the Coulomb potential into the Yukawa6 one. Interestingly enough, an unexpected feature is found. In fact, the profile (55) displays a new 1 2 long-range L3 correction, where its strength is proportional to m . It is also important to observe that a analogous correction has been found in Born-Infeld electrodynamics in the context of very special relativity [47]. In this way we establish a new connection between nonlinear effectives theories. Before we proceed further, we should comment on our result. In the case of QED (Euler-Heisenberg Lagrangian 4 1 1 16 e ~ density) the parameter 2 is given by 2 = 4 7 , where me is the electron mass. In this context, we also recall β β 45 mec 16 the currently accepted upper limit for the photon mass, that is, mγ 2 10− eV. Thus, for the QED case, from equation (55) it follows that the second term on the right-hand side∼it would× be detectable in long-range distances ( 109 m). In other words, we see that detectable corrections induced by vacuum polarization with a mass term would∼ be present at low energy scales. From equation (55), it clearly follows that the interaction energy between heavy charged charges, at leading order in β, is not finite at the origin. Motivated by this, one may consider the above calculation in a non-commutative geometry, based in findings of our previous studies [18, 19]. In such a case, the electric field at leading order in β2 and m2, takes the form 1 2 m2 − 3 6m2 1 eθ δ(3) (x) E = 1+ + Π2 Π2 ∂ ∇ , (56) i ∆ 2β2 − β2 ∆ i − 2 " # ∇ ! 2 (3) θ (3) where it may be recalled that we are now replacing the source δ (x y) by the smeared source e ∇ δ (x y), with θ the parameter non-commutative. Now, making use of equation (19− ), we readily find that − (t, r)= (1) (t, r)+ (2) (t, r)+ (3) (t, r) . (57) A0 A0 A0 A0 The term (1) (t, r) was first calculated in [48], we can, therefore, write only the result: A0 2 m θ 2 2 2 (1) e 1 mr 1 ∞ 1 u m r m m θ (t, r)= q e− du e− − 4u q e . (58) A0 4π r − √π r2 √u − 4π " Z /4θ # While the terms (2) (t, r) and (2) (t, r), after some manipulation, can be brought to the form A0 A0 3 x 12q 1 2 (2) (t, r)= nî dui γ3 3 , u , (59) 0 2 3/ 6 2 4θ A β π 2 0 u Z 2 3 x 3m q 1 2 4θ 2 2 (3) (t, r)= nî dui γ2 3 , u γ 3 , u γ 1 , u , (60) 0 2 5/ 4 2 4θ 2 2 4θ 2 4θ A β π 2 0 u u − Z 2 a 3 r a x du /b u where γ /2, /4θ is the lower incomplete Gamma function defined by γ ( /b, x) 0 u u e− . Inserting these expressions in equation (18), we finally obtain the static potential≡ for two opposite charges q located at 0 andL as R 2 m θ 2 2 4 L 2 q e mL 1 ∞ 1 u m u 12q i i 1 3 3 u V = e− du e− − 4u nˆ du γ , 2 2 3/ 6 2 4θ −4π L − √π L √u − β π 2 0 u " Z /4θ # Z x 2 4 2 2 2 3m q i i 1 2 3 u 4θ 3 u 1 u 5 nˆ du 4 γ 2 , 4θ 2 γ 2 , 4θ γ 2 , 4θ , (61) − β2π /2 0 u u − Z which is finite for L 0. It is a simple matter to verify that in the limit θ 0 we recover our above result. → → 9

IV. LOGARITHMIC CORRECTION

We now want to extend what we have done to Euler-Heisenberg-like electrodynamics at strong fields. As already mentioned, such theories show a power behavior that is typical for critical phenomena [49]. In such a case the Lagrangian density reads: 1 b F F µν = F F µν 0 F F µν log µν , (62) L −4 µν − 8 µν 4λ2 2 2 3 e mec where b0 and λ are constants. In fact, by choosing b0 = 6π2 and λ = e~ , we recover the Euler-Heisenberg electrodynamics at strong fields [49]. In the same way as was done in the previous section one can introduce an auxiliary field, ξ, to handle the logarithm in (62). This leads to 1 = α F F µν α (F F µν )2, (63) L −4 1 µν − 2 µν

b0 b0ξ where α =1 (1 + log ξ) and α = 2 . 1 − 2 2 32λ A similar procedure can be used to manipulate the quadratic term in (63). Accordingly, by introducing a second auxiliary ﬁeld, η, one easily ﬁnds 1 η2 = F (α +4α η) F µν + α . (64) L −4 µν 1 2 4 2

By setting, σ = α1 +4α2η, we then have

1 µν 1 2 = σFµν F + (σ α1) . (65) L −4 64α2 − It is once again straightforward to apply the gauge-invariant formalism discussed in the foregoing section. The µ 0µ 0 ∂L canonical momenta read Π = σF , and at once we recognize the two primary constraints Π = 0 and σ ∂σ˙ = 0. The canonical Hamiltonian corresponding− to (65) is P ≡ 1 σ 1 H = d3x Π ∂iA + Π2 + B2 (σ α )2 . (66) C i 0 2σ 2 − 64α − 1 Z 2 0 i Requiring the primary constraint Π to be preserved in time, one obtains the secondary constraint Γ1 = ∂iΠ = 0. In the same way, for the constraint , we get the auxiliary field σ as Pσ b b B2 3b B2 ξ σ = 1 0 (1 + ln ξ)+ 0 ξ 1+ 0 . (67) 2  2 2 3  − 2 2λ 2λ b0 b0B 1 (1 + ln ξ)+ 2 ξ  − 2 2λ    Hence we obtain 1 b 3b ξ H = d3x Π ∂iA + Π2 + 0 (1 + log ξ) Π2 0 Π4 . (68) C i 0 2 4 − 2λ2 Z As before, requiring the primary constraint ξ to be preserved in time, one obtains the auxiliary field ξ. In this λ P case ξ = 6Π2 . Consequently, we get 1 6b H = d3x Π ∂iA + (1 + b ) Π2 0 Π4 . (69) C i 0 2 0 − λ2 Z Following the same steps that led to equation (50) we find that 1 3 H (1) = Φ d3x Π2 Π4 Φ . (70) h iΦ h | 2 − 8β2 | i Z It should be noted that this expression is similar to equation (50) in the limit m 0, except by the chaged sign in frot of the Π4 term. Hence we see that the potential for two opposite charges in 0→and L is given by q2 1 q4 1 V = . (71) −4π L − 6040β2π2 L5 10

V. FINAL REMARKS

Finally, within the gauge-invariant but path-dependent variables formalism, we have considered the confinement versus screening issue for both massive Euler-Heisenberg-like and Euler-Heisenberg Electrodynamics in the approximation of the strong-field limit. Once again, a correct identification of physical degrees of freedom has been fundamental for understanding the physics hidden in gauge theories. Interestingly enough, their non-commutative version displays an ultraviolet finite static potential. The above analysis reveals the key role played by the new quantum of length in our analysis. In a general perspective, the benefit of considering the present approach is to provide unifications among different models, as well as exploiting the equivalence in explicit calculations, as we have illustrated in the course of this work.

Acknowledgments

It is a pleasure to thank J. A. Helay¨el-Neto for helpful comments on the manuscript. This work was partially supported by Fondecyt (Chile) grant 1130426 and DGIP (UTFSM) internal project USM 111458.

[1] G. Breit and J. A. Wheeler, Phys. Rev. 46, 1087 (1934). [2] S. L. Adler, Ann. Phys. (N.Y.) 67, 599 (1971). [3] V. Costantini, B. De Tollis and G. Pistoni, Nuovo Cimento A 2, 733 (1971). [4] R. Ruffini, G. Vereshchagin and S-S. Xue, Phys. Rept. 487, 1-140 (2010). [5] G. V. Dunne, Int. J. Mod. Phys. Conf. Ser. 14, 42 (2012). [6] C. Bamber et al., Phys. Rev. D 60, 092004 (1999). [7] D. L. Burke et al., Phys. Rev. Lett. 79, 1626 (1997). [8] O. J. Pike, F. Mackenroth, E. G. Hill and S. J. Rose, Nature Photonics 8, 434 (2014). [9] D. Tommasini, A. Ferrando, H. Michinel and M. Seco, J. High Energy Phys. 0911, 043 (2009). [10] D.Tommasini, A. Ferrando, H.Michinel and M.Seco, Phys. Rev. A 77, 042101 (2008). [11] M. Born and L. Infeld, Proc. R. Soc. London, Ser. A 144, 425 (1934). [12] H. Gies, J. Jaeckel and A. Ringwald, Phys. Rev. Lett. 97, 140402 (2006). [13] E. Masso and R. Toldra, Phys. Rev. D 52, 1755 (1995). [14] P. Gaete and E. I. Guendelman, Mod. Phys. Lett. A 20, 319 (2005). [15] P. Gaete and E. Spallucci, J. Phys. A 39, 6021 (2006). [16] E.S. Fradkin and A.A. Tseytlin, Phys. Lett. B 163 (1985) 12. [17] E. Bergshoeff, E. Sezgin, C.N. Pope and P.K. Townsend, Phys. Lett. B 188 (1987) 70. [18] P. Gaete and J. Helayël-Neto, Eur. Phys. J. C 74, no. 3, 2816 (2014). [19] P. Gaete and J. Helayël-Neto, Eur. Phys. J. C 74, no. 11, 3182 (2014). [20] S. H. Hendi, Ann. Phys. 333, 282 (2013). [21] Z. Zhao, Q. Pan, S. Chen and J. Jing, Nucl. Phys. B 871, 98 (2013). [22] O. Miˇskovićand R. Olea, Phys. Rev. D 83, 024011 (2011). [23] S. Habib Mazharimousavi and M. Halilsoy, Phys. Lett. B 678, 407 (2009). [24] G. Amelino-Camelia, Nature 418, 34 (2002). [25] T. Jacobson, S. Liberati and D. Mattingly, Phys. Rev. D 67, 124011 (2003). [26] T. J. Konopka and S. A. Major, New J. Phys. 4, 57 (2002). [27] S. Hossenfelder, Phys. Rev. D 73, 105013 (2006). [28] P. Nicolini, Int. J. Mod. Phys. A 24, 1229 (2009). [29] E. Witten,Nucl. Phys. B 268, 253 (1986). [30] N. Seiberg and E. Witten, JHEP 9909, 032 (1999). [31] M. R. Douglas, N. A. Nekrasov, Rev. Mod. Phys. 73, 977-1029 (2001). [32] R. J. Szabo, Phys. Rept. 378, 207-299 (2003). [33] J. Gomis, K. Kamimura and T. Mateos, JHEP 0103, 010 (2001). [34] A. A. Bichl, J. M. Grimstrup, L. Popp, M. Schweda and R. Wulkenhaar, Int. J. Mod. Phys. A 17, 2219 (2002). [35] A. Smailagic and E. Spallucci, J. Phys. A 36, L517 (2003). [36] A. Smailagic and E. Spallucci, J. Phys. A 36, L467 (2003). [37] A. Smailagic and E. Spallucci, J. Phys. A 37, 1 (2004) [Erratum-ibid. A 37, 7169 (2004)]. [38] R. Banerjee, S. Gangopadhyay and S. K. Modak, Phys. Lett. B 686, 181 (2010). [39] A. B. Hammou, M. Lagraa and M. M. Sheikh-Jabbari, Phys. Rev. D 66, 025025 (2002). [40] P. Nicolini, A. Smailagic and E. Spallucci, Phys. Lett. B 632, 547 (2006). [41] W. Dittrich and M. Reuter, ”Effective Lagrangians in Quantum Electrodynamics”, Springer-Verlag (1985). [42] P. Gaete, Phys. Rev. D59, 127702 (1999). 11

[43] S. Klarsfield, Phys. Lett. B 66, 86 (1977) [44] A. M. Frolov and D. M. Wardlaw, Eur. Phys. J. B. 85, 348 (2012). [45] P. Gaete, Z. Phys. C76, 355 (1997). [46] S. Ferrara and A. Sagnotti, JHEP 1504, 032 (2015). [47] R. Bufalo, Phys. Lett. B 746, 251 (2015). [48] P. Gaete and E. Spallucci, J. Phys. A 45, 065401 (2012). [49] H. Kleinert, E. Strobel and S. S. Xue, Phys. Rev. D 88, 025049 (2013). arXiv:1509.02594v2 [hep-th] 30 Dec 2015 hoy te ye fnniereetoyaishv endiscu supers Nevertheles been have in well-known. electrodynamics contexts is nonlinear different of it many types as in other view, arise theory, of actions point effective classical BI-type als the a Electrodynamics of th (BI) from dependence Born-Infeld Whereas theory space-time case, the Euler-Heisenbeg strengths. in the nonlinear field in (but the fields the of t in powers linear t and in (in ex correction corrections calculation as Euler-Kockel Serber nonlinear Incidentally, the and extends Lagrangian. Uehling result Euler-Heisenberg Euler-Heisenberg the the polariza are from vacuum contributions requir derived from of clearly contributions This class the a energy. incorporate interaction fac to the In Lagrangian of particles. form charged heavy structural between the energy interaction the in un u ob needn rmtecoc ftetp fproduc of type smearin the a Vor of to as choice leads known the approach) is from state which independent (coherent alternativ rule non-commutativity be an multiplication to that new out out a turns point ex through also asymptotic We kernel implemented intro n heat to been obtained. formulate of is been to form development has way this modified divergences of novel a ingredient where a key of [41], that (Moya The operators theory recall product [38–40]. quantum qu star to in a a context particular, proposed using include this in by not Within interest, studied does of been program. SM has gravity operators the focus quantum position that the the fact in become the theoretical have in length origin fundamental its and scenarios violating cnro uha onIfl hoy[1,mlihre atce [12 the particles is physics. millicharged this scattering Alo illustrates photon-photon [11], [6–10]. to that theory contributions confirmed example Born-Infeld significant been An as yet such [1–5]. not scenarios QED has of progress days remarkable early despite the interaction in nonlinear to as leading and pairs electron-positron virtual ∗ lcrncades [email protected] address: Electronic nti esetv,w lopitotta xesoso h Stand the of extensions that out point also we perspective, this In h hnmnno aumplrzto nQatmElectrodynamics Quantum in polarization vacuum of phenomenon The neetnl,i hudb ealdhr httepyia ffc of effect physical the that here recalled be should it Interestingly, ineeg o ohcasso lcrdnmc.A result a (1 featu As unexpected (Whichmann-Kroll) Electrodynamics. electrodynamics of like classes both Euler-Heisenberg-type for massive the energy tion calculate that we shown Subsequently, is phenomenon. birefringence It limit. field Electr Euler-Heisenberg-like and Euler-Heisenberg-like t h edn oaihi re) ipasaln-ag ( long-range a displays order), logarithmic leading the t (to in Electrodynamics Euler-Heisenberg-like Furthermore, ASnmes 47.e 26.n 13.40.Gp 12.60.Cn, 14.70.-e, numbers: PACS t versions, non-commutative their for Besides, potential. /r yuigtegueivrat u ahdpnet variabl path-dependent, but gauge-invariant, the using By 3 ye orcin pr rmaln-ag (1 long-range a from apart correction, type) - 1 eatmnod ´sc n etoCient´ıfico-Tecnológi Centro F´ısica and de Departamento oermrso olna electrodynamics nonlinear on remarks Some nvria ćiaFdrc at aá Valpara´ıso, Mar´ıa, Santa Técnica Federico Universidad Dtd eebr3,2015) 31, December (Dated: .INTRODUCTION I. arcoGaete Patricio /r 5 dnmc nteapoiaino h strong- the of approximation the in odynamics ye orcint h olm potential. Coulomb the to correction type) - 1, ,ti hscleetcagsbt h teghand strength the both changes effect physical this t, sdi h ieaue[18–23]. literature the in ssed eitrcineeg sutailtfinite. ultraviolet is energy interaction he ewe lcrmgei ed,rmisa exciting as remains fields, electromagnetic between s eapoiaino h togfil limit strong-field the of approximation he ffc hc seuvln ota nonee in encountered that to equivalent is which effect g steadto fcreto em nteMaxwell the in terms correction of addition the es ecntn akrudfil ii) hc contains which limit), field background constant he otissmlrnniercretost Maxwell to corrections nonlinear similar contains o e r on.W bananwlong-range new a obtain We found. are res ∗ oetodrmdfiain oteinterac- the to modifications lowest-order o h aeo asv Euler-Heisenbeg- massive of case the for , aso hc osntsffrfo hr distance short from suffer not does which pansion rdc)[23] nti oncini becomes it connection this In [32–37]. product) l lie n[] ti fitrs ontc htthe that notice to interest of is it [5], in plained rn hoy[6 7.I diint Born-Infeld to addition In 17]. [16, theory tring corrections contains result Serber and Uehling e 4] ti otyntn eeta hstp of type this that here noting worthy is It [43]. t 1 inpoes w motn xmlso such of examples important Two process. tion ncmuaieqatmfil hoyhsbeen has theory field quantum on-commutative gtesm ie eas ealta alternative that recall also we line, same the ng rainlk atce 1–5 a aemore have may [13–15] particles axion-like or ] /r lcrdnmc ipastevacuum the displays Electrodynamics eiaino h oeetsaeapoc has approach state coherent the of derivation e akrudfils.W ute eto ht as that, mention further We fields). background aumplrzto per samodification a as appears polarization vacuum fitnersac ciiy[43] hshas This [24–31]. activity research intense of sfraimw td ohmassive both study we formalism es eWcmn-rl orcin hc a be can which correction, Wichmann-Kroll he nu edtere loignon-commuting allowing theories field antum 5 uechrn ttso h unu position quantum the of states coherent duce ,B lcrdnmc sdsigihd since distinguished, is Electrodynamics BI s, r oe S)sc sLrnzinvariance Lorentz as such (SM) Model ard ye orcint h Coulomb the to correction type) - rvtto,a ocruvn difficulties circumvent to as gravitation, cteigo htn ypoos which photons, by photons of scattering QD,aiigfo h oaiainof polarization the from arising (QED), od Valpara´ıso, de co ssa-rdc 4] nhw physics Anyhow, [42]. star-product os Chile 2 a class of non-local theory. In other words, non-commutativity is just a sub-class of possible non-local deformation [44, 45]. More recently, this new approach has been successfully extended to black holes physics [46], also in connection to holographic superconductors via AdS-CFT [47]. Inspired by these observations, the purpose of this paper is to extend our previous studies [18, 19] on nonlinear electrodynamics to the case when vacuum polarization corrections are taken into account. The preceding studies were done using the gauge-invariant but path-dependent variables formalism, where the interaction potential energy between two static charges is determined by the geometrical condition of gauge invariance. One important advantage of this approach is that it provides a physically-based alternative to the usual Wilson loop approach. Accordingly, we shall work out the static potential for electrodynamics which include, apart from the Maxwell Lagrangian, additional terms corresponding to the Uehling, massive Euler-Heisenberg-like, Euler-Heisenberg Electrodynamics in the approximation of the strong-field limit (to the leading logarithmic order) and for their non-commutative versions. Our results show a 1 long-range L5-type correction to the Coulomb potential for both massive Euler-Heisenberg-like and Euler-Heisenberg Electrodynamics in the approximation of the strong-field limit (to the leading logarithmic order). Interestingly enough, 1 for massive Euler-Heisenbeg-like electrodynamics (Whichmann-Kroll) we obtain a new long-range L3 correction to the interaction energy. Nevertheless, for their non-commutative versions, the static potential becomes ultraviolet finite. The organization of the paper is as follows: In Section II, we reexamine Uehling Electrodynamics in order to establish a framework for the computation of the static potential. In Section III we consider Euler-Heisenberg-like (with a mass term) and show that it yields birefringence, compute the interaction energy for a fermion-antifermion pair and its version in the presence of a minimal length. In Section III, we repeat our analysis for Euler-Heisenberg Electrodynamics in the approximation of the strong-field limit. Finally, in Section IV, we cast our Final Remarks. In our conventions the signature of the metric is (+1, 1, 1, 1). − − −

II. BRIEF REVIEW ON THE UEHLING POTENTIAL

2 ∞ 1 2m 4m2 (m, x)= dτ 1+ 1 , (2) M 2 τ (τ + ∆) τ − τ Z4m r µ with ∆ ∂µ∂ . It should be noted that contains the effect of vacuum polarization to first order in the fine 2 ≡ e M structure constant, α = ~c , and, m, is the electron mass. We may parenthetically note here that the presence of ∆ in equation (2) does not offer problems. On the one hand, as we shall explain below we restrict to the static case. On the other hand, in order to compute the interaction energy, we shall make a series expansion at leading order in α.

Before going on, two remarks are pertinent at this point. First, the modification of Coulomb’s law in equation (1) follows from the weak-field limit of the one-loop effective action of quantum electrodynamics (QED). Indeed, as was explained in [48], this modification can be written as 1 = Aµ(x)Π (x, y)Aν (y)d4xd4y, (3) Lwf 2 µν Z 2 where wf denotes weak field, and Πµν is the usual order-e polarization tensor of QED. As is well known, in momentum space, this tensor is given by Π (k)= k2g k k Π k2 , (4) µν µν − µ ν while the momentum space spectral representation of the polarization function Π k2 reads

α ∞ ρ (t) 1 Π k2 = k2 dt , (5) −3π 2 t k2 + t Z4m 3

2 2 with ρ (t) = 1+ 2m 1 4m . Next, due to the tensor structure of Π (k), can then be expressed in a t − t µν Lwf gauge invariant way, that q is, in terms of Fµν . As a consequence of this, equation (3) reduces to the modiﬁcation of Coulomb’s law appearing in equation (1). We mention in passing that in [48] the signature of the metric is diﬀerent from that used in this paper. Second, it should be noted that the theory described by equation (1) contains higher time derivatives, hence to construct the Hamiltonian one must use, for example, the Ostrogradsky method [49]. Accordingly, in the theory under consideration the velocities have to be taken as independent canonical variables. Let us also mention here that, in previous studies [50, 51], we have shown that although theories like to equation (1) contain higher derivatives, in the electrostatic case the canonical momentum conjugate to velocities disappears. Hence the new Legendre transformation to construct the Hamiltonian reduces to the standard Legendre transformation. It should, however, be emphasized here that the present paper is aimed at studying the static potential of the above theory, so that ∆ can be replaced by 2. Notice that, for notational convenience, we have maintained ∆ in equations (1) and (2), but it should be borne−∇ in mind that this paper essentially deals with the static case.

Now, we move on to compute the canonical Hamiltonian. For this end we perform a Hamiltonian constraint µ α µ0 0 analysis. The canonical momenta are found to be Π = 1 3π ∆ F . It is easy to see that Π vanishes, we then have the usual constraint equation, which according to Dirac− ’s theoryM is written as a weak ( ) equation, that is, Π0 0. It may be noted that the remaining non-zero momenta must also be written as weak equations.≈ This leads to Π≈i 1 α ∆ Ei (with E = F ). Accordingly, the canonical Hamiltonian H is ≈ − 3π M i i0 C

1 1 1 H d3x Πi∂ A Π 1 α ∆ − Πi + F 1 α ∆ F ij , (6) C ≈ i 0 − 2 i − 3π M 4 ij − 3π M Z which must also be written as a weak equation. Next, the primary constraint, Π0 0, must be satisfied for all ≈ times. An immediate consequence of this is that using the equation of motion, Z˙ [Z,HC], we obtain the secondary i ≈ constraint Γ1 ∂iΠ 0, which must also be true for all time. In passing we recall that we are considering the static case, hence this≡ new constraint≈ does not contain time derivatives. It is straightforward to check that there are no further constraints in the theory. Therefore, in the case under consideration, there are two constraints, which are first class. According to the general theory we obtain the extended Hamiltonian as an ordinary (or strong) equation by adding all 3 the first-class constraints with arbitrary constraints. We thus write H = HC + d x (u0(x)Π0(x)+ u1(x)Γ1(x)), where uo(x) and u1(x) are arbitrary Lagrange multipliers. It is also important to observe that when this new Hamiltonian is employed, the equation of motion of a dynamic variable may be writtenR as a strong equation. With the aid of 0 equation (6) we find that A˙ 0 (x) = [A0 (x) ,H]= u0 (x), which is an arbitrary function. Since Π 0 always, neither A0 nor Π0 are of interest in describing the system and may be discarded from the theory. In fact, the≈ term containing A0 is redundant, because it can be absorbed by redefining the function w(x). Therefore, the Hamiltonian is now given as

1 1 1 H = d3x w(x)∂ Πi Π 1 α ∆ − Πi + F 1 α ∆ F ij , (7) i − 2 i − 3π M 4 ij − 3π M Z where w(x)= u (x) A (x). 1 − 0 It must be clear from this discussion that the presence of the new arbitrary function, w(x), is undesirable since we have no way of giving it a meaning in a quantum theory. Hence, according to the usual procedure, we impose a gauge condition such that the full set of constraints becomes second class. A convenient choice is

1 Γ (x) dzν A (z) dλxiA (λx)=0. (8) 2 ≡ ν ≡ i CZξx Z0 where λ (0 λ 1) is the parameter describing the space-like straight path zi = ξi + λ (x ξ)i, and ξ is a ﬁxed point (reference point).≤ ≤ There is no essential loss of generality if we restrict our considerations to−ξi = 0. The Dirac brackets can now be determined and we simply note the only nontrivial Dirac bracket involving the canonical variables, that is,

1 A (x) , Πj (y) ∗ = δj δ(3) (x y) ∂x dλxiδ(3) (λx y) . (9) i i − − i − Z 0 4

In passing we also recall that the transition to a quantum theory is made by the replacement of the Dirac brackets by the operator commutation relations according to A, B ∗ ( i/~)[A, B]. With the foregoing information, we can now proceed{ to} obtain→ − the interaction energy. As already mentioned, in order to accomplish this purpose, we will calculate the expectation value of the energy operator H in the physical state Φ , where the physical states Φ are gauge-invariant ones. The physical state can be written as | i | i y iq i Φ Ψ(¯ y)Ψ(y′) = ψ¯(y)exp( dz A (z))ψ(y′) 0 , (10) | i≡ ~ i | i Zy′ where 0 is the physical vacuum state and the line integral appearing in the above expression is along a space-like path starting| i at y and ending at y, on a fixed time slice. The point we wish to emphasize, however, is that the physical fermion (Ψ′ (y)) is not the Lagrangian fermion (ψ (y)), which is neither gauge-invariant nor associate with an electric field. In fact, the physical fermion is the Lagrangian fermion together with a cloud (or dressing) of gauge fields. Making use of the above Hamiltonian structure [18], we find that

y′ (3) Π (x) Ψ (y) Ψ (y′) = Ψ (y) Ψ (y′) Π (x) 0 + q dz δ (z x) Φ . (11) i i | i i − | i Zy

With the aid of equations (11) and (7), the lowest-order modiﬁcation in α of the interaction energy takes the form

H = H + V + V , (12) h iΦ h i0 1 2 where H = 0 H 0 . The V , V terms are given by h i0 h | | i 1 2 2 y y q 3 i (3) i (3) V = d x dz δ (x z) dz′ δ (x z′) , (13) 1 2 − − Z Zy′ Zy′ and

2 y′ y′ q α ∞ ρ (τ) 3 (3) 2 1 i (3) V2 = dτ d x dzi′δ (z′ x) x dz δ (z x) , (14) 2 3π 2 τ − ∇ τ 2 − Z4m Z Zy − ∇x Zy

2 2 where ρ (τ)= 1+ 2m 1 4m . τ − τ We note that the term q (13) may look peculiar, but it is nothing but the familiar Coulomb interaction plus a self- √τ z z 1 e− | − ′| energy term [52]. Now making use of the Green function, G (z, z′)= , the term (14) can be rewritten in 4π z z′ the form | − |

2 y′ y′ 2 √τ y y′ q α ∞ ρ (τ) 2 i α q ∞ ρ (τ) e− | − | V2 = dτ dzi′ z dz G (z, z′)= dτ . (15) 2 3π 2 τ ∇ ′ −3π 2 2 τ y y Z4m Zy Zy Z4m | − ′| Since the second and third term on the right-hand side of Eq. (12) are clearly dependent on the distance between the external static fields, the potential for two opposite charges located at y and y reads ′ 2 q 1 α ∞ ρ (τ) √τL V = 1+ dτ e− , (16) −4π L 3π 2 τ Z4m where L = y y′ . Accordingly, one recovers the known Uehling potential, which finds here an entirely different derivation. | − | Before we proceed further, we wish to show that this result can be written alternatively in a more explicit form. Making use of [53]

tz ∞ e 1 1 χ (z)= dt − 1+ 1 , (17) n tn 2t2 − t2 Z1 r we then get

q2 1 2α V = 1+ χ (2mL) . (18) −4π L 3π 1 5

By the transformation, t = cosh u [54], the functions χn can be reduced to the form [53]: 1 1 χn (z)= Kin 1 (z) Kin+1 (z) Kin+3 (z) , (19) − − 2 − 2 where the functions Ki denote Bessel function integrals. Hence we see that the interaction energy (with m = 1) becomes

2 2 3 3 q 1 2α L 5L 2L 3L 2L ∞ V = 1+ 1+ K (2L) + K (2L)+ + dtK (t) , (20) −4π L 3π 3 0 − 3 3 1 2 3 0 Z2L where K0(z) and K1(z) are modiﬁed Bessel functions. Finally, with the aid of asymptotic forms for Bessel functions, it is a simple matter to ﬁnd expressions for V for large and small L. Before concluding this subsection we discuss an alternative way of stating our previous result (16), which displays certain distinctive features of our methodology. We start by considering [52, 55]

V q ( (0) (L)) , (21) ≡ A0 − A0 where the physical scalar potential is given by

1 (t, r)= dλriE (t, λr). (22) A0 i Z0 This follows from the vector gauge-invariant ﬁeld expression

x (x) A (x)+ ∂ dzµA (z) , (23) Aµ ≡ µ µ − µ Zξ where the line integral is along a space-like path from the point ξ to x, on a fixed slice time. It is also important to observe that the gauge-invariant variables (22) commute with the sole first constraint (Gauss law), showing in this way that these fields are physical variables. Inasmuch as we are interested in estimating the lowest-order correction 1 i α − i to the Coulomb energy, we will retain only the leading term in expression E = 1 3π ∆ Π . Making use of this last expression, equation (22) gives − M 1 0 1 0 i λr J (λr) α ∞ dτ i λr J (λr) 0 (t, r)= dλr ∂i 2 + ρ (τ) dλr ∂i 2 , (24) A − 3π 2 τ − τ Z0 ∇λr Z4m Z0 ∇λr − i 0 to get the last line we used Gauss law for the present theory, that is, ∂iΠ = J (where we have included the external current J 0 to represent the presence of two opposite charges). Accordingly, for J 0(t, r) = qδ(3) (r), the potential for a pair of static point-like opposite charges located at 0 and L, is given by

2 q 1 α ∞ ρ (τ) √τL V = 1+ dτ e− , (25) −4π L 3π 2 τ Z4m after subtracting a self-energy term.

III. EULER-HEISENBERG-LIKE MODEL

β2 1 1 p = 1 1+ 2 , (26) L 2 − β2 F − β2γ2 G 1 µν 1 ˜µν where have included two parameters β and γ. As usual, = 4 Fµν F , = 4 Fµν F , Fµν = ∂µAν ∂ν Aµ and ˜µν 1 µνρλ F G − F = 2 ε Fρλ. Let us also mention here that in our previous paper [19] we have studied the domain 0

µν ∂µF˜ =0, (28) where 2 Γ=1+ F G . (29) β2 − β2γ2 Also, it is straightforward to see that Gauss law becomes,

D =0, (30) ∇ · where D is given by

E2 B2 (E B)2 2 D = 1 − · E + (E B) B . (31) − 2β2 − β2γ2 γ2 · " # 0 (3) Q e Again, from equation (30), for J (t, r) = eδ (r), we ﬁnd D = r2 rˆ, where Q = 4π . This then implies that for a point-like charge, e, at the origin, the expression Q E2 = 1 E , (32) r2 − 2β2 | | tells us that, for r 0, the electrostatic ﬁeld becomes singular at r = 0, in contrast to the 0

With the aid from equations (33) and (34) we find the electric permittivity, εij , and the inverse magnetic permeability, 1 µ− ij , tensors of the vacuum, that is, 2BiBj 1 2EiEj ε =Γ δ + , µ− =Γ δ , (37) ij ij γ2 ij ij − γ2 7 with Di = εij Ej and Bi = µij Hj . In accordance with our previous procedure [18, 19], we can now linearize the above equations. To do this, it is advantageous to introduce a weak electromagnetic wave (Ep, Bp) propagating in the presence of a strong constant external field (E0, B0). On these assumptions, we readily find that, for the case of a purely magnetic field (E0 = 0), the vectors D and H become B2 2 D = 1+ 0 E + (E B ) B , (38) 2β2 p γ2 p · 0 0 and

B2 1 H 0 B B B B = 1+ 2 p + 2 ( p 0) 0 , (39) 2β  2 B0 ·  β 1+ 2β2   where we have keep only linear terms in Ep, Bp. As before, we consider the z axis as the direction of the external magnetic ﬁeld (B0 = B0e3) and assuming that the light wave moves along the x axis, the decomposition into a plane wave for the ﬁelds Ep and Bp can be written as

i(wt k x) i(wt k x) Ep (x,t)= Ee− − · , Bp (x,t)= Be− − · . (40) In this case, it clearly follows that k2 ε µ E =0, (41) w2 − 22 33 2 and k2 ε µ E =0. (42) w2 − 33 22 3 As a consequence, we have two diﬀerent situations: First, if E B0 (perpendicular polarization), from (42) E3 = 0, 2 k ⊥ and from (41) we get w2 = ε22µ33. This then means that the dispersion relation of the photon takes the form

2 B0 1+ 2β2 n = 2 . (43) ⊥ v 3B0 u1+ 2 u 2β u t 2 k Second, if E B (parallel polarization), from (41) E = 0, and from (42) we get 2 = ε µ . This leads to || 0 2 w 33 22 2B2 n = 1+ 0 . (44) k s γ2 Thus we verify that in the case of a generalized Euler-Heisenberg electrodynamics the phenomenon of birefringence is present. We now pass to the calculation of the interaction energy between static point-like sources for a massive Whichman- Kroll-like model, our analysis follows closely that of references [18, 19]. The corresponding theory is governed by the Lagrangian density 1 1 m2 = F 2 + (F F µν )2 + A Aµ. (45) L −4 µν 32β2 µν 2 µ Next, in order to handle the second term on the right hand in (45), we introduce an auxiliary ﬁeld ξ such that its equation of motion gives back the original theory. This allows us to write the Lagrangian density as 1 ξ 1 m2 = F F µν + F F µν ξ2 + A Aµ. (46) L −4 µν 32β2 µν − 128β2 2 µ

ξ With the redeﬁnition η =1 2 , equation (46) becomes − 8β 1 1 m2 = ηF F µν (1 η)2 + A Aµ. (47) L −4 µν − 2 − 2 µ 8

Before we proceed to work out explicitly the interaction energy, we shall ﬁrst restore the gauge invariance in equation (47). Following an earlier procedure, we readily verify that the canonical momenta read Πµ = ηF 0µ, which results in the usual primary constraint Π0 =0, and Πi = ηF i0. In this way one obtains −

1 η m2 1 H = d3x Πi∂ A + Π2 + B2 A Aµ + (1 η)2 . (48) C i 0 2η 2 − 2 µ 2 − Z i 2 0 The consistency condition, Π˙ = 0, leads to the constraint Γ ∂iΠ + m A = 0. As a result, both constraints are second class. To convert the second class system into first class≡ we shall adopt the procedure described previously. Thus, we enlarge the original phase space by introducing a canonical pair of fields θ and Πθ. It follows, therefore, that 2 a new set of first class constraints can be defined in this extended space: Λ1 Π0 + m θ =0 and Λ2 Γ + Πθ = 0. Notice that this new θ-field is not to be confused with the usual non-commutative≡ parameter. This then≡ shows that the new constraints are first class and, therefore, restore the gauge symmetry. As is well known, this procedure reproduces the usual Stückelberg formalism. From this, the new effective Lagrangian density, after integrating out the θ fields, becomes

1 m2 1 = F η + F µν (1 η)2. (49) L −4 µν ∆ − 2 − 2 m Now, writing σ = η + ∆ , the expression (49) can be brought to the form

1 k m2 2 = F σF µν 1 σ + , (50) L −4 µν − 128 − ∆ where k = 64β2. We are now ready to compute the interaction energy. In this case, the canonical momenta are Πµ = σF 0µ, with the usual primary constraint Π0 =0, and Πi = σF i0. Hence the canonical Hamiltonian is expressed as −

1 σ k m2 2 H = d3x Πi∂ A + Π2 + B2 + 1 σ + . (51) C i 0 2σ 2 128 − ∆ Z ( ) 0 i Time conservation of the primary constraint Π yields the secondary constraint Γ1 ∂iΠ = 0. Similarly for the σ constraint yields no further constraints and just determines the ﬁeld σ. In this case,≡ at leading order in β, the ﬁeldPσ is given by

m2 B2 3 1 σ = 1+ 1 Π2 , (52) 2  2 3  ∆ − 2β − 2β m2 B2 1+ 2  ∆ − 2β    which will be used to eliminate σ. As before, the corresponding total (ﬁrst-class) Hamiltonian that generates the time 3 evolution of the dynamical variables is H = HC + d x (u0(x)Π0(x)+ u1(x)Γ1(x)), where uo(x) and u1(x) are the Lagrange multiplier utilized to implement the constraints. In the same way as was done in the previous subsection,R the expectation value of the energy operator H in the physical state Φ becomes | i 1 1 m2 − 15 15m2 1 H = Φ d3x Πi 1+ Πi + Π4 Π2 Π2 Φ , (53) h iΦ h | 2 ∆ 8β2 − 2β2 ∆ | i Z ( ) in this last line we have considered only quadratic terms in m2. In such a case, by employing (53), the lowest-order modiﬁcation in β2 and m2 of the interaction energy takes the form

H = H + V + V + V , (54) h iΦ h i0 1 2 3 where H = 0 H 0 . The V , V and V terms are given by h i0 h | | i 1 2 3 2 y 2 1 y q 3 i (3) m − i (3) V = d x dz′ δ (x z′) 1 dz δ (x z) , (55) 1 2 − − 2 − Z Zy′ ∇ Zy′ 9

4 y y y y 15q 3 (3) i (3) k (3) (3) V = d x dz δ (x z) dz′ δ (x z′) du δ (x u) dv δ (x v) , (56) 2 − 8β2 i − − − k − Z Zy′ Zy′ Zy′ Zy′ and 2 4 y y y y 15m q 3 (3) i (3) k (3) (3) V = d x dz δ (x z) dz′ δ (x z′) du δ (x u) dv δ (x v) . (57) 3 2β2 i − − − k − Z Zy′ Zy′ Zy′ Zy′ Finally, with the aid of the expressions (55), (56) and (57), the potential for a pair of static point-like opposite charges located at 0 and L, is given by q2 e mL q4 3 1 1 V = − + 5m2 . (58) −4π L 16πβ2 8π L2 − L3 observe that when m = 0, the profile (58) reduces to the known Whichmann-Kroll interaction energy. On the other hand, for m = 0, it should be noted the key role played by the mass term in transforming the Coulomb potential into the Yukawa6 one. Interestingly enough, an unexpected feature is found. In fact, the profile (58) displays a new 1 2 long-range L3 correction, where its strength is proportional to m . It is also important to observe that a analogous correction has been found in Born-Infeld electrodynamics in the context of very special relativity [57]. In this way we establish a new connection between nonlinear effectives theories. Before we proceed further, we should comment on our result. In the case of QED (Euler-Heisenberg Lagrangian 4 1 1 16 e ~ density) the parameter 2 is given by 2 = 4 7 , where me is the electron mass. In this context, we also recall β β 45 mec 16 the currently accepted upper limit for the photon mass, that is, mγ 2 10− eV. Thus, for the QED case, from equation (58) it follows that the second term on the right-hand side∼it would× be detectable in long-range distances ( 109 m). In other words, we see that detectable corrections induced by vacuum polarization with a mass term would∼ be present at low energy scales. From equation (58), it clearly follows that the interaction energy between heavy charged charges, at leading order in β, is not finite at the origin. Motivated by this, one may consider the above calculation in a non-commutative geometry, based in findings of our previous studies [18, 19]. In such a case, the electric field at leading order in β2 and m2, takes the form 1 2 m2 − 3 6m2 1 eθ δ(3) (x) E = 1+ + Π2 Π2 ∂ ∇ , (59) i ∆ 2β2 − β2 ∆ i − 2 " # ∇ ! 2 (3) θ (3) where it may be recalled that we are now replacing the source δ (x y) by the smeared source e ∇ δ (x y), with θ the parameter non-commutative. Now, making use of equation (22− ), we readily find that − (t, r)= (1) (t, r)+ (2) (t, r)+ (3) (t, r) . (60) A0 A0 A0 A0 The term (1) (t, r) was first calculated in [50], we can, therefore, write only the result: A0 2 m θ 2 2 2 (1) e 1 mr 1 ∞ 1 u m r m m θ (t, r)= q e− du e− − 4u q e . (61) A0 4π r − √π r2 √u − 4π " Z /4θ # While the terms (2) (t, r) and (2) (t, r), after some manipulation, can be brought to the form A0 A0 3 x 12q 1 2 (2) (t, r)= nî dui γ3 3 , u , (62) 0 2 3/ 6 2 4θ A β π 2 0 u Z 2 3 x 3m q 1 2 4θ 2 2 (3) (t, r)= nî dui γ2 3 , u γ 3 , u γ 1 , u , (63) 0 2 5/ 4 2 4θ 2 2 4θ 2 4θ A β π 2 0 u u − Z 2 a 3 r a x du /b u where γ /2, /4θ is the lower incomplete Gamma function defined by γ ( /b, x) 0 u u e− . Inserting these expressions in equation (21), we finally obtain the static potential≡ for two opposite charges q located at 0 andL as R 2 m θ 2 2 4 L 2 q e mL 1 ∞ 1 u m u 12q i i 1 3 3 u V = e− du e− − 4u nˆ du γ , 2 2 3/ 6 2 4θ −4π L − √π L √u − β π 2 0 u " Z /4θ # Z x 2 4 2 2 2 3m q i i 1 2 3 u 4θ 3 u 1 u 5 nˆ du 4 γ 2 , 4θ 2 γ 2 , 4θ γ 2 , 4θ , (64) − β2π /2 0 u u − Z which is finite for L 0. It is a simple matter to verify that in the limit θ 0 we recover our above result. → → 10

IV. LOGARITHMIC CORRECTION

We now want to extend what we have done to Euler-Heisenberg-like electrodynamics at strong fields. As already mentioned, such theories show a power behavior that is typical for critical phenomena [58]. In such a case the Lagrangian density reads: 1 b F F µν = F F µν 0 F F µν log µν , (65) L −4 µν − 8 µν 4λ2 2 2 3 e mec where b0 and λ are constants. In fact, by choosing b0 = 6π2 and λ = e~ , we recover the Euler-Heisenberg electrodynamics at strong fields [58]. In the same way as was done in the previous section one can introduce an auxiliary field, ξ, to handle the logarithm in (65). This leads to 1 = α F F µν α (F F µν )2, (66) L −4 1 µν − 2 µν

b0 b0ξ where α =1 (1 + log ξ) and α = 2 . 1 − 2 2 32λ A similar procedure can be used to manipulate the quadratic term in (66). Accordingly, by introducing a second auxiliary ﬁeld, η, one easily ﬁnds 1 η2 = F (α +4α η) F µν + α . (67) L −4 µν 1 2 4 2

By setting, σ = α1 +4α2η, we then have

1 µν 1 2 = σFµν F + (σ α1) . (68) L −4 64α2 − It is once again straightforward to apply the gauge-invariant formalism discussed in the foregoing section. The µ 0µ 0 ∂L canonical momenta read Π = σF , and at once we recognize the two primary constraints Π = 0 and σ ∂σ˙ = 0. The canonical Hamiltonian corresponding− to (68) is P ≡ 1 σ 1 H = d3x Π ∂iA + Π2 + B2 (σ α )2 . (69) C i 0 2σ 2 − 64α − 1 Z 2 0 i Requiring the primary constraint Π to be preserved in time, one obtains the secondary constraint Γ1 = ∂iΠ = 0. In the same way, for the constraint , we get the auxiliary field σ as Pσ b b B2 3b B2 ξ σ = 1 0 (1 + ln ξ)+ 0 ξ 1+ 0 . (70) 2  2 2 3  − 2 2λ 2λ b0 b0B 1 (1 + ln ξ)+ 2 ξ  − 2 2λ    Hence we obtain 1 b 3b ξ H = d3x Π ∂iA + Π2 + 0 (1 + log ξ) Π2 0 Π4 . (71) C i 0 2 4 − 2λ2 Z As before, requiring the primary constraint ξ to be preserved in time, one obtains the auxiliary field ξ. In this λ P case ξ = 6Π2 . Consequently, we get 1 6b H = d3x Π ∂iA + (1 + b ) Π2 0 Π4 . (72) C i 0 2 0 − λ2 Z Following the same steps that led to equation (53) we find that 1 3 H (1) = Φ d3x Π2 Π4 Φ . (73) h iΦ h | 2 − 8β2 | i Z It should be noted that this expression is similar to equation (53) in the limit m 0, except by the changed sign in front of the Π4 term. Hence we see that the potential for two opposite charges in→0 and L is given by q2 1 q4 1 V = . (74) −4π L − 6040β2π2 L5 11

V. FINAL REMARKS

Finally, within the gauge-invariant but path-dependent variables formalism, we have considered the confinement versus screening issue for both massive Euler-Heisenberg-like and Euler-Heisenberg Electrodynamics in the approximation of the strong-field limit. Once again, a correct identification of physical degrees of freedom has been fundamental for understanding the physics hidden in gauge theories. Interestingly enough, their non-commutative version displays an ultraviolet finite static potential. The analysis above reveals the key role played by the new quantum of length in our analysis. In a general perspective, the benefit of considering the present approach is to provide a unification scenario among different models, as well as exploiting the equivalence in explicit calculations, as we have illustrated in the course of this work.

Acknowledgments

[1] G. Breit and J. A. Wheeler, Phys. Rev. 46, 1087 (1934). [2] S. L. Adler, Ann. Phys. (N.Y.) 67, 599 (1971). [3] V. Costantini, B. De Tollis and G. Pistoni, Nuovo Cimento A 2, 733 (1971). [4] R. Ruffini, G. Vereshchagin and S-S. Xue, Phys. Rept. 487, 1-140 (2010). [5] G. V. Dunne, Int. J. Mod. Phys. Conf. Ser. 14, 42 (2012). [6] C. Bamber et al., Phys. Rev. D 60, 092004 (1999). [7] D. L. Burke et al., Phys. Rev. Lett. 79, 1626 (1997). [8] O. J. Pike, F. Mackenroth, E. G. Hill and S. J. Rose, Nature Photonics 8, 434 (2014). [9] D. Tommasini, A. Ferrando, H. Michinel and M. Seco, J. High Energy Phys. 0911, 043 (2009). [10] D.Tommasini, A. Ferrando, H.Michinel and M.Seco, Phys. Rev. A 77, 042101 (2008). [11] M. Born and L. Infeld, Proc. R. Soc. London, Ser. A 144, 425 (1934). [12] H. Gies, J. Jaeckel and A. Ringwald, Phys. Rev. Lett. 97, 140402 (2006). [13] E. Masso and R. Toldra, Phys. Rev. D 52, 1755 (1995). [14] P. Gaete and E. I. Guendelman, Mod. Phys. Lett. A 20, 319 (2005). [15] P. Gaete and E. Spallucci, J. Phys. A 39, 6021 (2006). [16] E.S. Fradkin and A.A. Tseytlin, Phys. Lett. B 163 (1985) 12. [17] E. Bergshoeff, E. Sezgin, C.N. Pope and P.K. Townsend, Phys. Lett. B 188 (1987) 70. [18] P. Gaete and J. Helayël-Neto, Eur. Phys. J. C 74, no. 3, 2816 (2014). [19] P. Gaete and J. Helayël-Neto, Eur. Phys. J. C 74, no. 11, 3182 (2014). [20] S. H. Hendi, Ann. Phys. 333, 282 (2013). [21] Z. Zhao, Q. Pan, S. Chen and J. Jing, Nucl. Phys. B 871, 98 (2013). [22] O. Miˇskovićand R. Olea, Phys. Rev. D 83, 024011 (2011). [23] S. Habib Mazharimousavi and M. Halilsoy, Phys. Lett. B 678, 407 (2009). [24] G. Amelino-Camelia, Nature 418, 34 (2002). [25] T. Jacobson, S. Liberati and D. Mattingly, Phys. Rev. D 67, 124011 (2003). [26] T. J. Konopka and S. A. Major, New J. Phys. 4, 57 (2002). [27] S. Hossenfelder, Phys. Rev. D 73, 105013 (2006). [28] P. Nicolini, Int. J. Mod. Phys. A 24, 1229 (2009). [29] S. Pramanik, S. Ghosh and P. Pal, Annals Phys. 346, 113 (2014). [30] A. Farag Ali, S. Das and E. C. Vagenas, Phys. Lett. B 678, 497 (2009). [31] G. Amelino-Camelia, Living Rev. Relativity, 16, 5 (2013). [32] E. Witten,Nucl. Phys. B 268, 253 (1986). [33] N. Seiberg and E. Witten, JHEP 9909, 032 (1999). [34] M. R. Douglas, N. A. Nekrasov, Rev. Mod. Phys. 73, 977-1029 (2001). [35] R. J. Szabo, Phys. Rept. 378, 207-299 (2003). [36] J. Gomis, K. Kamimura and T. Mateos, JHEP 0103, 010 (2001). [37] A. A. Bichl, J. M. Grimstrup, L. Popp, M. Schweda and R. Wulkenhaar, Int. J. Mod. Phys. A 17, 2219 (2002). [38] A. Smailagic and E. Spallucci, J. Phys. A 36, L517 (2003). [39] A. Smailagic and E. Spallucci, J. Phys. A 36, L467 (2003). [40] A. Smailagic and E. Spallucci, J. Phys. A 37, 1 (2004) [Erratum-ibid. A 37, 7169 (2004)]. [41] E. Spallucci, A. Smailagic and P. Nicolini, Phys. Rev. D 73, 084004 (2006). [42] R. Banerjee, S. Gangopadhyay and S. K. Modak, Phys. Lett. B 686, 181 (2010). 12

[43] A. B. Hammou, M. Lagraa and M. M. Sheikh-Jabbari, Phys. Rev. D 66, 025025 (2002). [44] L. Modesto, J. W. Moffat and P. Nicolini, Phys. Lett. B 695, 397 (2011). [45] P. Nicolini, “Nonlocal and generalized uncertainty principle black holes,” arXiv:1202.2102 [hep-th]. [46] P. Nicolini, A. Smailagic and E. Spallucci, Phys. Lett. B 632, 547 (2006). [47] S. Pramanik, S. Das and S. Ghosh, Phys. Lett. B 742, 266 (2015). [48] W. Dittrich and M. Reuter, ”Effective Lagrangians in Quantum Electrodynamics”, Springer-Verlag (1985). [49] M. Ostrogradsky, Mem. Acad. St. Pet. VI 4, 385 (1850). [50] P. Gaete and E. Spallucci, J. Phys. A 45, 065401 (2012). [51] P. Gaete, Int. J. Mod. Phys. A 27, 1250061 (2012). [52] P. Gaete, Phys. Rev. D59, 127702 (1999). [53] S. Klarsfield, Phys. Lett. B 66, 86 (1977) [54] A. M. Frolov and D. M. Wardlaw, Eur. Phys. J. B. 85, 348 (2012). [55] P. Gaete, Z. Phys. C76, 355 (1997). [56] S. Ferrara and A. Sagnotti, JHEP 1504, 032 (2015). [57] R. Bufalo, Phys. Lett. B 746, 251 (2015). [58] H. Kleinert, E. Strobel and S. S. Xue, Phys. Rev. D 88, 025049 (2013).