Complex-mass shell renormalization of the higher-derivative electrodynamics

We consider a higher-derivative extension of QED modified by the addition of a gauge-invariant dimension-6 kinetic operator in the U(1) gauge sector. The Feynman diagrams at one-loop level are then computed. The modification in the spin-1 sector leads the electron self-energy and vertex corrections diagrams finite in the ultraviolet regime. Indeed, no regularization prescription is used to calculate these diagrams because the modified propagator always occurs coupled to conserved currents. Moreover, besides the usual massless pole in the spin-1 sector, there is the emergence of a massive one, which becomes complex when computing the radiative corrections at one-loop order. This imaginary part defines the finite decay width of the massive mode. To check consistency, we also derive the decay length using the electron--positron elastic scattering and show that both results are equivalent. Because the presence of this unstable mode, the standard renormalization procedures cannot be used and is necessary adopt an appropriate framework to perform the perturbative renormalization. For this purpose, we apply the complex-mass shell scheme (CMS) to renormalize the aforementioned model. As an application of the formalism developed, we estimate a quantum bound on the massive parameter using the measurement of the electron anomalous magnetic moment and compute the Uehling potential. At the end, the renormalization group is analyzed.

We consider a higher-derivative extension of QED modified by the addition of a gauge-invariant dimension-6 kinetic operator in the U (1) gauge sector. The Feynman diagrams at one-loop level are then computed. The modification in the spin-1 sector leads the electron self-energy and vertex corrections diagrams finite in the ultraviolet regime. Indeed, no regularization prescription is used to calculate these diagrams because the modified propagator always occurs coupled to conserved currents. Moreover, besides the usual massless pole in the spin-1 sector, there is the emergence of a massive one, which becomes complex when computing the radiative corrections at one-loop order. This imaginary part defines the finite decay width of the massive mode. To check consistency, we also derive the decay length using the electron-positron elastic scattering and show that both results are equivalent. Because the presence of this unstable mode, the standard renormalization procedures cannot be used and is necessary adopt an appropriate framework to perform the perturbative renormalization. For this purpose, we apply the complex-mass shell scheme (CMS) to renormalize the aforementioned model. As an application of the formalism developed, we estimate a quantum bound on the massive parameter using the measurement of the electron anomalous magnetic moment and compute the Uehling potential. At the end, the renormalization group is analyzed.

I. INTRODUCTION
Effective field theories (EFT) play a central role in modern physics. They cover almost all branches in physics such as nuclear systems derived from low-energy quantum chromodynamics [1], chiral perturbation theory [2][3][4][5], BCS theory formulated from conventional superconductivity [6], inflationary model in cosmology [7,8], gravitationally induced decoherence [9], and so on. Even our most fundamental theories, General Relativity and the Standard Model, are thought of as leading terms of some underlying theory [10,11].
The ideas concerning EFT have began with a nonlinear modification of Maxwell electromagnetism made in order to understand the photon-photon electrodynamical scattering process by Euler and Heisenberg [12] in the 1930s. At the same time, Fermi developed the theory of beta decay to describe the elementary process n → p+ e − +ν e in the framework of quantum field theory [13]. Even having some interesting features, Fermi and Euler-Heisenberg theories were not taken seriously since they were not renormalizable. Nevertheless, some years later, the development of the renormalization and the renormalization group techniques [14], along with the theorem derived by Appelquist and Carazzone [15]-which states that heavy mass particles can be decoupled from low energy dynamics under certain conditions-gave rise to the current EFT programme [16].
Effective theories allow us to simplify the description of a given physical process by taking into account the appropriate variables at a given energy scale, i.e., one can consider only the relevant degrees of freedom at a specific energy range. It is basically a low energy dynamics valid below some energy scale and which does not depend on the behavior in the ultraviolet regime. EFT have the advantage of reducing the number of degrees of freedom, turning the description of the physical system under consideration easier to deal with. An appropriate choice of degrees of freedom is a crucial point in the understanding of the problem.
One interesting set of EFT models are the so-called higher-order theories. This class of theories are characterized basically by the modification of the free propagator through the introduction of higher-order kinetic terms in the Lagrangian, which leads to a modified propagator that exhibits a better asymptotic behavior. However, these theories are usually plagued with ghosts states, giving rise to non-Hermitian interactions and scattering processes that do not preserve probabilities, violating the unitarity. There are many phenomena that are described in this framework such as dark energy [17][18][19], ultraviolet regulators [20][21][22], renormalizable gravity models [23], string theory [24], and supersymmetry [25], for instance. It also have been suggested that quantum gravity effects can give rise to partners to every field in the Standard Model [26,27], which can be included in a straightforward way by the aforementioned formalism, generating an extension of the SM. In addition, some work has recently demonstrated that higher-order derivative theories may emerge in the scope of noncommutative spacetimes [28][29][30][31][32]. In this framework, the introduction of a minimal length naturally generates higher-order kinetic terms, as in the case of electrodynamics and general relativity [33], for instance. In this sense, a better understanding of some basic features of higher-order theories would certainly help us to improve our knowledge of this set of models.
Concerning QED, a possible way to find an effective Lagrangian is through the addition of a gauge-invariant dimension-6 operator containing higher derivatives in the free Lagrangian of the U (1) sector, namely, where M is the only parameter added to theory and is responsible for introducing a cutoff in the ultraviolet regime of QED. This modification is very similar to Pauli-Villars regularization procedure [34], Lee-Wick model [35,36] at quantum level, and Podolsky electromagnetism [37,38] in the classical context. The new degree of freedom in QED changes dramatically the behavior of the theory at short distances. Nevertheless, unitarity is not preserved. [70] Many remarkable features are found in this QED extension such as the coexistence of Dirac magnetic monopoles and massive vector bosons [41], the presence of finite electromagnetic mass and self-energy of point particles [37,38]. Also, it provides an adequate scenario to solve the 4/3 problem [42][43][44] Our goal in this paper is precisely to address: (i) the one-loop radiative corrections, (ii) the perturbative renormalization of the HDQED in the complex-mass shell (CMS) scheme, (iii) a quantum bound on the Mparameter using the measurement of the electron anomalous magnetic moment, (iv) the decay width of the unstable mode taking into account the electron-positron scattering, (iv) the computation of the Uehling potential in the HDQED framework, and the (v) analysis of the renormalization group.
This work is organized as follows. In Sect. 2 we give a brief description of our model. In the following section, Sect. 3, we compute the second-order radiative corrections in the HDQED context. Also, the finite decay width of the unstable mode is explicity calculated. At the end of Sect. 3 we present the Uehling potential. The complex-mass shell renormalization scheme together with the renormalization group is presented in Sect. 4. We summarize our results in Sect. 5.

II. THE HIGHER-DERIVATIVE QED
The higher-derivative electromagnetism usually is defined by the U (1)-gauge-invariant Lagrangian density where F µν = ∂ µ A ν − ∂ ν A µ is the field strength and M is a mass parameter which introduces a length scale in the model. Using the Bianchi identity, the higher-order kinetic term can be rewritten as which implies that Lagrangian (2) takes the form The effective Lagrangian (4), besides being Lorentz invariant, gives origin to local field equations that are linear in the field quantities and are given by whereF µν = 1 2 ǫ µναβ F αβ . From the above fourth-order equations of motion, it is clear that the system will carry more degrees of freedom than Maxwell electromagnetism. In this sense, it is instructive understand what is the particle content of this model at the quantum level. To accomplish this, we first add the gauge-fixing term L ζ = − 1 2ζ (∂ µ A µ ) 2 to the Lagrangian (4), where ζ is the gaugefixing parameter, and then compute the propagator in momentum space, which assumes the form Contracting (7) with the external conserved currents J µ , one obtains The scattering amplitude (8) has poles at k 2 = 0 and k 2 = M 2 . Taking into account that J is space-like (J 2 < 0) [45][46][47], we have the residues From the above residues it is clear that the model under consideration carry two spin-1 modes, one massless and one massive. The positive sign residue corresponds to the spin-1 massless excitation, which can be related to the QED photon. The spin-1 massive mode possess a wrong sign residue which give rise to negative norm states (ghost states) and consequently causes unitarity violation. However, since we are treating the HDQED as an effective field theory, it will not concern us. Regarding to renormalizability, the modified spin-1 propagator always occurs coupled to conserved current, which means that terms proportional to the momenta k µ give vanish contributions. Hence, the modified propagator reduces to which is the difference between Maxwell and Proca propagators. The first term on the right of Eq. (10) is the QED propagator, which is renormalizable by power counting. The second corresponds to a massive spin-1 propagator. Massive vector theories have a bad behavior at high energies and do not go to zero asymptotically. However, there are two exceptions: (i) gauge theories with spontaneous symmetry breaking and (ii) neutral vector bosons coupled to conserved currents. The condition (ii) ensures the renormalizability of the Proca model, showing that the model is renormalizable by power counting.
According to Eq. (10), the nonrelativistic potential U (r) between two electrically charged particle can be expressed as which have the following properties: 1. The regularized potential is finite at the origin which is evidence that point-like particles has electromagnetic mass and self-energy finities.
2. The standard Coulomb potential is recovered when M r >> 1.
3. At short distances, i.e., M r << 1, U differs significantly from the Coulomb potential, as one can see from Fig. 1.

III. SECOND-ORDER RADIATIVE CORRECTIONS OF THE HDQED
The Lagrangian density characterizing the higherderivative spinor quantum electrodynamics can be written as where the first term on the right is the Dirac kinetic operator, m is the electron mass and the last term denotes the interacting term. In the lowest-order perturbation theory, the scattering of charged particles by each other is characterized by the matrix element of the operator in momentum space where J µ = eψγ µ ψ is the conserved charged fermionic current and D µν (k) is the Feynman propagator.
According to (15), the modified propagator always occurs coupled to conserved currents and the momenta dependence vanishes, reducing the propagator to which provides a better behavior at high energies since the propagator goes as k −4 at the asymptotic limit. By power counting, all higher-order process of the Lagrangian (2) are finite, except for the charge renormalization. In other words, the QED primary divergences related to the electron self-energy and vertex corrections are convergent. Nevertheless, the vacuum polarization diagram, as in QED case, is ultraviolet divergent. We will apply the dimensional regularization to cancel off the divergence arising from the vacuum polarization. Now, we shall undertake the computation of second-order radiative corrections of the HDQED.

A. Electron self-energy
The electron self-energy expression Σ(p) at one-loop approximation in the HDQED is given by where λ is an infrared regulator. A dimensional analysis shows that the modified propagator makes the integral (17) convergent at high energies. The explicit form of the electron self-energy integral is where α = e 2 /4π is the fine structure constant. Using the on-shell condition (p ′2 = p 2 = m 2 ) and remembering that the higher-order electrodynamics bring about small deviations from QED, which implies that the mass is very large compared to the electron mass (M ≫ m), the integral (18) reduces to which yields a finite value to the electron mass. Indeed, this result is not a surprise since the self-force acting on a point charge particle is finite and well defined in the HDQED scenario. As an aside, we remark that the 4/3 problem finds a natural explanation in the HDQED context since the electromagnetic mass enters in the Aharonov-Lorentz equation in a form consistent with special relativity [42][43][44]. Nevertheless, the infrared divergence remains.

B. Vertex correction
In the HDQED, the second-order vertex correction is given by As in the case of electron self-energy, the vertex correction (20) is ultraviolet convergent at short distances due to the modified propagator. Again, no regularization procedure is necessary. In order to obtain an explicit form of the vertex correction, we can use the Gordon identity of the Dirac current, which yields where the functions F 1 (q 2 ) and F 2 (q 2 ) are called the form factors and q ν := p ′ν − p ν is the four-momentum transfer at the vertex. At radiative corrections, the information concerning to vertex corrections are all contained in F 1 (q 2 ) and F 2 (q 2 ). So, computing the vertex correction at one-loop level provide us with the following form factors: From F 1 (q 2 ) and F 2 (q 2 ) it is clear that ultraviolet divergences do not affect the vertex corrections. The F 2 (q 2 ) form factor when q 2 → 0 is related to deviation from the electron magnetic moment standard prediction of the Dirac equation.

Anomalous magnetic moment of the electron
The electron anomalous magnetic moment is the most precise measurement in QED, having an accuracy up to 12 decimal places. This astonishing outcome can be used in the context of HDQED to put a quantum bound on the M -parameter and then compare with that of QED. To accomplish this, we first note that for an electron scattered by an external static magnetic field and at the limit q 2 → 0, the gyromagnetic ratio is [48] The form factor of the electron, F 2 (0), corresponds to a shift in the g-factor, usually quoted in the form F 2 (0) = g−2 2 . So, taking the limit q 2 → 0 in (23), it can be shown that The one-loop correction of the vertex diagram.
Integrating the above expression first with respect to x, gives Computing F 2 (0), we arrive at the conclusion that The first term of the above equation is equal to that calculated by Schwinger in 1948 [49]. Since then F 2 (0) has been calculated to order α 10 for QED. The second term of Eq. (27) is the most important correction related to the parameter M of the HDQED electrodynamics. Recent calculation concerning F 2 (0) in the framework of QED gives for the electron [50] where the uncertainty comes mostly from that of the best non-QED value of the fine structure constant α. The current experimental value for the anomalous magnetic moment is, in turn [51,52] Comparison of the theoretical value predicted by QED with the experimental one shows that these results agree in 1 part in 10 12 . As a consequence, Consequently, a lower limit on the M -parameter is M ≈ 438 GeV.

C. Vacuum polarization
As discussed previously, in the HDQED framework the modification introduced by the higher-order kinetic terms improves the behavior of the physical process in the high frequencies regime. As consequence, the above Feynman diagrams become finite at one-loop order. However, because the fermionic sector is unaltered in the HDQED, the vacuum polarization tensor will be identical to the one in QED, i.e., where the contribution at one-loop level of the polarization scalar Π(k 2 ) is Bearing in mind that the propagator always occurs coupled to conserved current, the dressed propagator in HDQED assumes the form which clearly reveals that the radiative corrections coming from the fermionic loops give contribution to Mparameter.
Performing the integration of (32) using dimensional regularization we find where the function f (k 2 ) is given by It is important to note that there is an imaginary part which emerges when k 2 > 4m 2 . This complex term is a clear indication that the massive mode develops a finite decay width, an aspect that will be discussed in detail in the next section. Also, the appearance of a complex massive pole implies that one cannot just apply the standard renormalization schemes. The appropriate framework to treat this issue needs to take into account the existence of unstable particles, a point that deserves some considerations and that will be done in Sect. 6.

Decay width of the unstable mode
Before proceeding with the renormalization procedure we will discuss the finite decay rate of the unstable mode. To begin with, we first introduce the auxiliary field formalism. This technique has the advantage of eliminate the higher-order kinetic terms through the introduction of auxiliary fields, reducing the Lagrangian to quadratic terms in the fields. So, defining the auxiliary field inserting in the Lagrangian (4), and defining A µ := B µ + Z µ , the Lagrangian (4) assumes the form where B µ and Z µ are related to the massless and massive fields, respectively. Here we promptly note that the introduction of auxiliary fields provide us a clear interpretation of the distinct modes at different energy levels.
According to (36), at low energy, the massless mode dominates over the massive one the description of the system, reproducing the results of QED. However, at high frequencies, the massive mode emerges and gives contributions to the physical observables as the anomalous magnetic moment of the electron, for instance. Thereafter, at short distances the unstable mode should develop a finite decay width. In the auxiliary field formalism the decay width description can be simplified by noting that at the ultraviolet regime only the massive propagator is relevant, namely, which is the Proca propagator with a minus sign. In electron-positron scattering, when the mass of the massive mode exceeds the electron-positron mass (M > 2m e ), the negative residue appears in the physical sheet. However, as we have discussed in the previous section, at this energy range the massive mode becomes complex and must decay in light particles. The decay width Γ can be obtained by the standard self-energy sum. Hence, at the narrow width approximation, the resummed massive propagator assumes the form Comparing the denominator of (38) with the imaginary part that arises in the self-energy sum when k 2 > m 2 in (34), we promptly find that the finite decay width of the M -parameter is To check the consistency of the previous decay width derivation, one can follow another route to obtain the decay rate of the M -parameter through the analysis of the electron-positron elastic scattering. The general expression for the decay rate of this process is where k is the massive spin-1 four-momentum and q and q ′ are related to the electron and positron momenta, respectively. The electron-positron elastic scattering amplitude is where ǫ µ is the polarization vector and u and v represent the fermions functions. Using the completeness relation, the decay factor Γ assumes the form Solving the above integral, we arrive at the decay rate (46) which is identical to the decay rate (39) found when considering the one-loop corrections in the spin-1 sector, so lending support to our derivation. It is worth to note that since M ≫ 2m, the Γ factor reduces to In general, there may be many decays modes, which means that the massive mode lifetime τ is given by where using our estimative for the mass value M ≈ 438 GeV gives a lifetime of τ 10 −24 s.

Uehling potential
One of the great triumphs of QED was the prediction of the small difference between the energy levels 2 S 1/2 and 2 P 1/2 in the hydrogen atom-the Lamb shift-which is related to the radiative corrections to the Coulomb potential coming from the vacuum polarization. Since in the HDQED framework the energy potential is not singular at the origin [44], we are interested in probe what sort of contribution will arise in this context. To start with, the general expression for the energy potential is given by Inserting the resummed propagator (32) in the energy potential expression (49), and neglecting terms of α 3 order, one obtains where the ξ-integral above is the integral representation of the Uehling potential in the HDQED framework. Since the above integral is hard to solve analytically, we will focus only on the asymptotic limit mr ≫ 1.
For mr ≫ 1, only the region 0 ≤ ξ − 1 ≪ (mr) −1 gives contributions to the integral. So, considering ξ ≃ 1, we obtain In the above limit, the ξ-integral becomes where Γ is the incomplete Gamma function and Φ is the confluent hypergeometric function of the first kind (Kummer's function). Since M > 2m in (50), the Uehling potential also gives rise to an imaginary part, which goes to zero when M goes to infinity. At the limit in which the massive mode is much bigger than the mass of the fermions, i.e., M ≫ m, the perturbative expansion at m/M gives for the real part of the Uehling potential, at the first order, the expression while for the imaginary part it gives which exponentially decay at the asymptotic limit.
To conclude we remark that in the limit that M goes to infinity, the ξ-integral becomes a Gamma function and the potential (50) reduces to which is the standard Uehling potential in QED.

A. Renormalized perturbation theory
The presence of a finite decay width implies that one needs to find an appropriate framework to describe unstable particles in the perturbation theory. It follows that applying directly the standard summation of self-energy violates the gauge invariance [53][54][55]. To circumvent this problem, many methods have been proposed over the last years [56][57][58][59][60][61][62][63][64][65][66]. Currently, the most general procedure to treat perturbative renormalization of unstable particles in quantum field theory is the so-called complex-mass shell (CMS) scheme. The CMS is an extension of the onshell procedure for unstable particles which is fully gauge invariant over all the phase space. In this formalism, the renormalization constant and the complex massive pole are defined at the location of the unstable propagator (for further information, see [67][68][69]).
Performing the perturbative renormalization one obtains the renormalized Lagrangian where the relations between the bare and renormalized quantities are and We call attention to the fact that, according to the usual renormalization procedure, one has to put counterterms into all diagrams, even the finite ones, since in an interacting theory the bare mass and the coupling constants are not equal to the physical parameters. A quick glance at the renormalized Lagrangian (56) shows us that five conditions are necessary in order to fix all the counterterms, which are given by The first constraint fixes the physical electron mass, the second and third fix the Dirac and sector-1 propagator to have a residue equal to 1, the fourth fixes the electron charge, and the last one defines the complex massive pole. Since the massive pole is complex, the renormalization constant will also be complex, which agrees with our claim that one needs to find a convenient framework to deal with the renormalization of unstable particles. In this sense, writing M 0 = √ Z M M , where Z M is the complex massive renormalization factor defined by one obtains the relation M 2 /M 2 0 := 1 − Π k 2 = M 2 . On the other hand, the M -parameter is renormalized by the renormalization function Z A . As well as QED, the Z A -factor is given by and the renormalization factors Z A and Z M are related by the expression We also point out that the Ward identity ensures that the factor Z ψ is equal to the vertex correction renormalization Z V , i.e., Z ψ = Z V . As a consequence, the HDQED does not give a further contribution to the charge renormalization.

B. Renormalization group
The Callan-Symanzik equation of the renormalization group is given by where Γ (n) is the one-particle irreducible Green function and µ is an arbitrary energy scale. The functions (β, γ M , γ A ) are related with the renormalization factors by Thus we may conclude that the beta function is unchanged in the HDQED framework. This happens because at high momenta the M -parameter is negligible, being given explicitly by β(e) = e 3 /12π 2 . Moreover, as well as in QED, the HDQED is not asymptotically free. Combining the functions (64) and Eqs. (58), it is easy to see that β(e) = e γ A (e), implying that γ A (e) = α/3π and The invariance of the Green function Γ (n) under a scale transformation, Γ (n) (e, M, µ) = Γ (n) (ē(t),M (t),μ(t) = µ e t ), leads to an effective constant couplingē(t) and an effectiveM (t)-parameter as a function of the dimensionless scale t, both satisfying the equations ∂ē ∂t = β(ē), whereē(t = 0) = e andM (t = 0) = M . The solution of (66) givesē Solving Eq. (67) we find which gives a running mass as a function of the energy scale.

V. FINAL REMARKS
In this paper we have investigated some issues related to one-loop corrections of HDQED. Since the higherorder kinetic terms in the spin-1 sector improves the behavior of the QED propagator at short distances, we have shown that the vertex correction and the electron self-energy are finite and depends explicitly on the Mparameter. In this sense, it became clear that the Mparameter acts as a natural regulator parameter which at the limit where M → ∞ recovers the QED divergences. Nevertheless, the vacuum polarization remains divergent. It is also worth to emphasize that we did not make use of any regularization procedure in the computations of these diagrams. The only assumption made was that the scattering processes should occur in the presence of external conserved currents. Despite the fact that the couplings with external conserved currents have eliminated the momenta dependence in the loop integrals, which simplified our task of solving the integrals, at the next leading order this prescription cannot be used anymore. Higherorder scattering processes induce the appearance of internal loops, which are not canceled when one takes into account the presence of conserved currents. We also call attention to the fact that our computations were made in the Lorenz gauge. However, any gauge condition is feasible [41].
From the analysis of the complex-mass shell conditions in the perturbative renormalization scheme, we have concluded that the HDQED enforces us to input another constraint in order to fix all the counterterms in the renormalized Lagrangian. Concerning the charge renormalization, the massless pole is the only one relevant while for the renormalization of the M -parameter one needs to take into account the complex massive pole. Analyzing the Callan-Symanzik equation, we have observed that the beta function and the anomalous dimension are unchanged by the presence of the M -parameter. Thus we may conclude that the M -parameter is negligible at the ultraviolet regime. Nevertheless, the gamma function associated to the M -parameter given by (65) and the effective M -parameter has an interesting behavior running with the energy scale. It is worth to remark that to the best of our knowledge it is the first time that the CMS scheme is applied to treat higher-order derivative electrodynamics theories.
As an interesting application we have found a quantum bound on the M -parameter using the value of the electron anomalous magnetic moment. This measurement has been used to estimate stringent constraints on possible theories beyond the standard model. Using the latest measurements of the mentioned phenomenon we found M = 438 GeV. Also, we have calculated the Uehling potential in the HDQED framework. Although the effect of Uehling potential to the Lamb shift is negligible, in muonic atoms the vacuum polarization is dominant, which in principle could be used to put a new bound in the mass parameter of HDQED. This new constraint on the M -parameter can be further investigated in the future. Regarding the modern particle physics experiments, we point out that the detection of this new degree of freedom could be related to physics beyond the Standard Model. In our prescription we have included the higher-order term after the spontaneous symmetry breaking has occurred in the electroweak theory (EW). However, this new degree of freedom could be added in the EW Lagrangian before the symmetry breakdown. In this way, the appearance of a mass −2 dimensional parameter at TeV scale could be a sign of a phenomenon arising from new physics, as effects of a large extra dimension, for instance.
In order to outline possible directions for future investigations, we would like to emphasize that the scenario developed here to treat unstable modes could be extended, at least in principle, to other spin fields. Also, issues concerning the unitarity of higher-order theories in the CMS formalism can be investigated in the future.