Pauli-Villars Regularization Elucidated in Bopp-Podolsky's Generalized Electrodynamics

We discuss an inherent Pauli-Villars regularization in Bopp-Podolsky's generalized electrodynamics. Introducing gauge-fixing terms for Bopp-Podolsky's generalized electrodynamic action, we realize a unique feature for the corresponding photon propagator with a built-in Pauli-Villars regularization independent of the gauge choice made in Maxwell's usual electromagnetism. According to our realization, the length dimensional parameter $a$ associated with Bopp-Podolsky's higher order derivatives corresponds to the inverse of the Pauli-Villars regularization mass scale $\Lambda$, i.e. $a = 1/\Lambda$. Solving explicitly the classical static Bopp-Podolsky's equations of motion for a specific charge distribution, we explore the physical meaning of the parameter $a$ in terms of the size of the charge distribution. As an offspring of the generalized photon propagator analysis, we also discuss our findings regarding on the issue of the two-term vs. three-term photon propagator in light-front dynamics.


Introduction
Quantum electrodynamics (QED) may be regarded as a prototype of quantum field theories with well-established renormalization program which effectively regulates the infinities present in the local gauge field theory. Due to the infinities that cannot be gotten around, e.g. radiative corrections in QED, one needs to treat and tame such infinities taking a certain regularization procedure with the renormalization condition for physical amplitudes. The very impressive agreement between high precision measurements in accelerators and the predictions of quantum field theory in the presence of radiative corrections is the key for the indication of successful renormalization program. Phenomenological success of atomic model appears ultimately backed up by the successful QED renormalization program.
Historically, the problem of infinities first arose in the classical electrodynamics of point particles in the 19th and early 20th century. The well-known example is the mass of electron including the electromagnetic mass m em due to its own electrostatic field given by m em = e 2 8πre with the charge e and the radius r e of the electron, which becomes an infinity as r e → 0.
It may not be an overstatement that the early work of Lorentz and Abraham [1] including the bare mass of the spherical shell as well as m em to take a consistent point limit provided the inspiration for later development of the renormalization program in QED and other local field theories. Modifying the concept of point charge to an extended charge distribution lends also the physical meaning of charge renormalization as the charge screening due to the Dirac vacuum in QED.
In the same vein, Bopp [2] and Podolsky [3] attempted to remove infinities inherent in the usual treatment of point charges introducing higher order derivatives in the Lagrangian of electrodynamics while maintaining the equations of motion still linear in the fields and preserving gauge invariance. In particular, Podolsky discussed the classical aspects of his model presenting the equations of motion, energy-momentum tensor and plane wave field solutions [3]. Traditionally, however, it has become the case to view the model due to Bopp and Podolsky ("BP model") as a mechanism to describe massive photons without breaking gauge invariance as the propagating modes of the model comprise both massless photons as well as massive ones. In this work, we demonstrate that the BP model solution for a point charge in electrodynamics corresponds to an ordinary electrodynamic solution for a specific charge distribution. Motivated by this possible reinterpretation of BP model solution in electrodynamics, we further elucidate the BP model as a natural way of providing Pauli-Villars regularization in ordinary QED.
To discuss our work, it may be worthwhile to make a brief historical remark on previous works on the BP model. A few years later after the introduction of the BP Lagrangian, Podolsky and Kikuchi [4] went through the actual quantization of the model. They claimed that the usual quantization methods of the time could not be directly applied to the BP model of generalized electrodynamics due to the presence of higher order derivatives in the Lagrangian and therefore they needed to introduce extra auxiliary fields. The quantization was performed in an extended phase space and, to take gauge symmetry issues into account, a generalization of the Stueckelberg formalism was used. A review of BP model's original results by Podolsky and Schwed can be found in Ref. [5]. Some forty years later, the BP model was revisited when Galvão and Pimentel [6] first carefully analyzed its structure of constraints performing the instant-form canonical quantization with Dirac Brackets. In terms of the canonical Dirac-Bergmann formalism [7,8], BP model has three first-class constraints generating gauge symmetries [6]. The canonical quantization was performed after gauge fixing and promotion of Dirac Brackets into operator commutators. It is worth noting though that in Ref. [6] the term of BP model's with higher order derivatives was considered with the opposite sign. In fact, the electrodynamics of BP model with the opposite sign was then further investigated in a series of papers [6,9,10], leading to the discussion of tachionic propagating modes for the photon. With the advent of modern and more powerful quantization methods, Barcelos-Neto, Galvão and Natividade [9] performed the Batalin-Fradkin-Vilkovisky (BFV) quantization using two slightly different gauges, namely the usual Lorenz gauge and the other called generalized Lorenz gauge. More recently, Bufalo and Pimentel [11] extended the BFV analysis including matter fields. From the symplectic quantization point of view, it is also worth mentioning that an interesting duality connection between the BP model and the massive Proca model has been investigated in Ref. [12]. More recent further discussions on the BP model can be found in Refs. [13,14,15,16,17,18,19,20].
In the case of the Lorenz gauge, a natural gauge-fixing term originally introduced by Podolsky and Kikuchi [4] which considerably simplifies the calculation in the quantization process has long been passed without notice in the literature and has been only recently rescued by Bufalo, Pimentel and Soto [17]. The role of this term in obtaining a simple generalized photon propagator in a straightforward manner cannot be overemphasized. In particular, we show that this term permits a nice factorization of the generalized photon propagator in all gauges analyzed in this work. Thus, the BP model parameter dependent part of the propagator appears only as a global multiplicative factor turning its mathematical structure easier to analyze and interpret as a way of introducing Pauli-Villars regularization. We discuss that it is possible to split the propagator as a sum of two parts consisting of massive and massless modes which we interpret as the Pauli-Villars regularization.
Furthermore, to our knowledge neither the axial nor the light-front gauges have been discussed so far by the functional path-integral quantization point of view in the context of BP model in the literature. The canonical structure of BP model's generalized electrodynamics on the null-plane has been recently analyzed by Bertin, Pimentel and Zambrano [14] where the light-front Hamiltonian form evolution is considered. In Ref. [14], after unraveling the constraint structure in phase space, the generalized radiation gauge on the null-plane is adopted. In the present work, we follow a different approach considering the theory defined by the Lagrangian in the configuration space and introducing the gauge-fixing conditions in the integration measure of the generating functional via the Faddeev-Popov procedure generalized to the BP model case. The quantization is then performed in a covariant way and for instance the generalized Lorenz gauge can be achieved. The axial-gauges are obtained along the same lines, the breaking of relativistic covariance being only perceived by a particular choice of the axial direction vector n µ . In addition, we have the opportunity to extend the ideas introduced in Refs. [21,22] concerning the adoption of two simultaneous gaugefixing conditions leading to a so-called doubly transverse photon propagator in the light-front gauge. We show in this work explicitly how to handle the corresponding calculations in the BP model case.
Our work is organized as follows. In Section 2, we define our notation and conventions, review some basic facts about the BP model as a gauge field theory and physically interpret some of its classical properties. In particular, we demonstrate that the BP model solution for a point charge in electrodynamics corresponds to an ordinary electrodynamic solution for a specific charge distribution. This motivates us in Section 3 to elucidate the BP model as a natural way of providing Pauli-Villars regularization in ordinary QED. We discuss in Section 3 the covariant Lorenz gauge fixing obtaining the corresponding photon propagator and point out the necessity of the natural gauge-fixing term in the gauge fixing action. This term permits a natural factorization of the generalized photon propagator in all gauges analyzed in Section 3. The BP model parameter dependent part of the propagator appears only as a global multiplicative factor turning its mathematical structure easier to analyze and interpret as a way of introducing Pauli-Villars regularization. In Section 4, we provide an example of the BP model application discussing the second-order correction to the electron self-energy and show explicitly the consistency with the Pauli-Villars regularized result. We close in Section 5 with some final comments and concluding remarks. In Appendix A, we summarize the derivation of Eq.(31) used in Section 2.

Bopp-Podolsky's Generalized Electrostatics
The starting classical action in Bopp-Podolsky (BP)'s generalized electrodynamics containing second order space-time derivatives of the gauge field A µ (x) in Minkowski space is given by where a is a real number with physical dimension of length or inverse mass, known as Bopp-Podolsky's parameter [2,3]. We use Minkowski's coordinates with metric signature diag(η µν ) = (+1, −1, −1, −1) and the integration measure d 4 x in Eq.(1) runs throughout all space-time coordinates x µ . It is clear that BP's classical action S 0 is a natural higher derivatives Lorentz covariant generalization of ordinary Maxwell's electromagnetism -the latter being recovered for a = 0. Although with a different notation, we adopt the same original Bopp-Podolsky's [2,3,4,5] choice for the second term in Eq.(1). This choice is important if one wishes to interpret the extra degrees of freedom of BP model as physical massive excitations. Although the case of negative sign for the second term in S 0 above was considered in [6,9,10] describing tachionic mass excitations for the gauge field, we are not going to discuss it here but rather maintain the implementation consistent with the usual causality. Note that in Eq.(1), the short-hand F µν ≡ ∂ µ A ν − ∂ ν A µ stands for the ordinary electromagnetic field strength tensor which is naturally invariant under the gauge group U (1). That means the BP extra a-dependent term does not spoil the original gauge invariance of the action S 0 and, particularly, the propagator of the gauge field is not well defined before gauge fixing. We shall address the gauge fixing issue in Section 3 where the Lorenz and axial type gauge fixings will be discussed. In the following we briefly review a few immediate properties and consequences of action given by Eq.(1) and consider the static case obtaining the BP version of Poisson's equation as well as its general solution. For a point charge delta distribution, the BP model leads to a everywhere finite potential -we shall show that it is possible to generate this very same potential within the scope of ordinary electrodynamics using a suitable charge distribution.

Field Equations of Motion and General Static Solution
In order to understand the physical content encoded in Eq.(1) we couple the gauge field to an external source j µ (x) through and demand stationarity of the total action under functional variations of the gauge field A µ . As a result, the corresponding Euler-Lagrange equations of motion obtained from the minimal action principle read and can be understood in the present context as Maxwell's generalized equations. As the indices µ, ν run from 0 to 3, Eq.(3) represents a set of four fourth-order partial differential equations on the gauge field A µ . Conservation of the external current j µ comes straightforwardly from the antisymmetry of the field strength In terms of the measurable usual electric and magnetic fields respectively given by and the action given by Eq.(1) can also be expressed as while Eq. (3) gives rise to the four non homogeneous Maxwell equations, namely, one for the temporal component and three for the space components The homogeneous Maxwell equations, on the other hand, remain the same as they amount to the very identities which permit to describe the physical E and B fields, through Eqs. (5) and (6) above, in terms of the gauge potential field A µ . In particular the absence of magnetic monopoles still holds in BP's generalization as the magnetic field B remains divergenceless. Let's now focus on the static case where it is possible to obtain the general solution and discuss the physical meaning of the BP model. For a time-independent electromagnetic field, Eqs. (8) and (9) reduce respectively to and containing now only space derivatives, while the homogeneous ones simply state that B is divergenceless and E irrotational. In this case, inserting (5) disregarding time derivatives into Eq.(10) leads to a generalized Poisson equation where φ ≡ A 0 represents the electrostatic potential and we have written j 0 = 4πρ for the electric charge density. By defining the a-dependent fourth-order differential operator P a as we may rewrite Eq.(12) more compactly as Henceforth, we shall refer Eq. (14) as the Poisson-Bopp-Podoslky (PBP) equation.
To solve the PBP equation given by Eq. (14), equivalently Eq. (12), for an arbitrary given charge distribution ρ, let us first consider the case of an elementary distribution resulting from a fixed unity point charge localized at r 0 For this case the solution φ P,a (r) can be written as the difference between the usual Coulomb potential and the Yukawa potential evaluated at r = |r − r 0 |. In fact, the BP potential centered at r 0 , defined as satisfies Eq. (14) for the elementary Dirac delta charge distribution given by Eq. (15). This can be directly seen by applying the P a operator to each potential leading to Figure 1: Podolsky potential for a = 0.5 as a function of the distance. Although both Coulomb's and Yukawa's potentials blow up at the origin, their difference remains finite. In the limit a → 0 Podolsky's potential bends continually towards Coulomb's potential.
Applying the Green's function method, we then find that the general solution φ(r) of the PBP equation given by Eq. (14) subject to the boundary condition of vanishing potential at infinity can be written as a superposition of kernel elementary contributions (18) weighed by the given charge density ρ(r), that is where dτ denotes the integration volume element with respect to the dummy integration variable r . The ordinary electrostatic solution is recovered in the limit that the BP model parameter a goes to zero, i.e. a → 0, as expected.
For the point charge density given by Eq. (15), the ordinary electrostatic solution given by the Coulomb electrostatic potential φ C (r) = 1 r in Eq.(16) diverges at the local point r = r 0 , i.e. r = |r − r 0 | = 0, imposing the problem of infinities discussed in our Introduction, Sec.1, for the classical electrodynamics of point particles. On the other hand, for non-null a, the BP potential φ P,a (r) in Eq.(18) remains finite in the limit r → 0 approaching to the finite value φ P,a (0) = 1 a and reproduces back the Coulomb's characteristic 1/r behavior for large values of r compared to a. In Fig.1, we plot BP's potential as a function of r as well as its two constituent parts Coulomb's and minus Yukawa's for the numerical value a = 0.5 in the same unity as the measured distance variable r. Here, it may be worth noticing the important fact that Coulomb's potential doesn't have any length scale parameter while the BP's potential has a natural length scale provided by the real parameter a. What follows in the next subsection is that this parameter a introduced in the BP model action given by Eq.(1) can be equivalently reinterpreted as the length scale of a specific charge distribution removing the fiasco of divergence for a point charge particle in ordinary local gauge electrodynamics.
In particular, we demonstrate that the BP model solution for a point charge corresponds to an ordinary electrodynamic solution for a specific charge distribution. This motivates us in the subsequent section, Sec. 3, to elucidate the BP model as a natural way of providing Pauli-Villars regularization in ordinary QED.

The BP Potential from Ordinary Electrostatics
As discussed in our Introduction, Sec.1, a non-null BP a-parameter improves the model convergence properties by means of the higher-order derivatives term. As a matter of fact we have just seen Coulomb's potential given by Eq.(16) gets smoothed out becoming finite at the critical point r = r 0 when we generalize Poisson's equation including the convergence parameter a into the PBP equation given by Eq. (12) or Eq. (14). In other words, in the BP model, Eq.(16) generalizes to Eq. (18). However, we show below that it is also possible to obtain the same effect within the realm of ordinary electrostatics by using a specific charge distribution. Indeed, let's take a small positive length dimension parameter b and consider the normalized charge distribution Then, the general solution to Poisson's equation gives the corresponding potential Using the spherical coordinates and writing dτ = r 2 sin θ dφ dθ dr , one can see straightforwardly that the azimuthal part in φ brings a 2π contribution and the integral in Eq.(24) is directly worked out as Performing then the last two r integrations, we obtain finally Note the key fact that we dispensed completely the BP higher-derivative a-dependent term in obtaining Eq.(27) as we used only the standard ordinary electrodynamics. Nevertheless, it represents the very BP potential with b playing the role of the previous BP a-parameter.
To discuss this key point further, let's now address the following question: i.e., Which potential would be produced by the charge distribution given by Eq.(23) in the BP model context using the general solution given by Eq. (21)? The answer to this question turns out to be an (a, b)-dependent potential which we denote here by ψ (a, b, r). This sort of double BP potential can be explicitly calculated using the general solution given by (21) as The first part of Eq.(28) has been previously calculated and shown to represent φ P,b (r). Therefore, we write ψ(a, b, r) = φ P,b (r) −ψ(a, b, r) The integration of Eq.(30) is explicitly given in appendix A, which shows that it leads tō By plugging this result back into Eq. (29), we obtain the final result for the double BP potential as Specifically for r close to zero, we get a finite result without any divergence: We note here that in the a → 0 limit, Eq.(32) reproduces Eq. (27)

Gauge Fixing and Functional Quantization
As already stated in the last paragraph of our Introduction, Sec.1, due to gauge invariance, a direct propagator for the gauge field in BP model is ill-defined. That happens because we are working with a constrained system and to preserve explicit covariance we use more field variables than degrees of freedom [20]. In order to proceed with the quantization of the model, similarly to ordinary electrodynamics, we must choose a specific gauge suitable for perturbative calculations. In the following subsections, we show how to achieve the generalized Lorenz and axial gauges performing the functional quantization of BP model. In both cases, we shall obtain the Green functions generating functional by means of a suitable generalization of the Faddeev-Popov method.

Generalized Lorenz Gauge
Concerning the covariant Lorenz gauge we generalize the standard gauge-fixing term in the Maxwell Lagrangian to which is then added to the original Lagrangian in the integrand of action given by Eq.(1). Besides BP's a-parameter, the gauge-fixing given by Eq.(34) is allowed to depend on the additional free real gauge parameter ξ. For the particular case ξ = 1 the necessity of this natural term can already be seen in the original papers of Podolsky, Kikuchi and Schwed [4,5]. Although this isolated term cannot be found directly in Ref. [4], a careful reading shows that it is in fact inserted and summed up in their so called modified Lagrangian up to total divergences. However, the second part of Eq.(34) has been tacitly omitted in the more modern treatments of BP's model and only recently has it been reintroduced in BP's context by Bufalo, Pimentel and Soto [17].
With the addition of (34) the fixed action S LGF = S 0 + d 4 x L LGF now reads By demanding stationarity with respect to the gauge field, that is, by enforcing we obtain the equations of motion which generalize Eq.(3) in the case of null external source. In terms of the gauge field A µ the equations of motion given by Eq.(37) can also be rewritten as For the Feynman-t'Hooft gauge choice ξ = 1, it reduce to the simpler result This means that in the Lorenz gauge the free field equations of motion comprise a superposition of plane waves describing both massive and massless particles. As previously mentioned we remark the importance of the sign choice in Eq.(1) for this physical interpretation. In fact the general solution to Eq.(39) can be written as for arbitrary functions a µ (k) and b µ (k). The generalized photon propagator of the BP model in the current working Lorenz gauge can be read directly from the inverse of the differential operator acting on A µ in Eq.(38). Or equivalently, by performing integrations by parts and discarding boundary terms, the action S LGF in Eq.(35) may be rewritten as In momentum space, this gives rise to the operator and permits us to write the photon propagator as One can straightforwardly check that P µλ M λν = δ ν µ , confirming that Eq.(34) leads to a neat simple expression for the photon propagator above with the additional massive simple pole 1/a 2 . The Landau gauge result is obtained for ξ = 0.
In the following, we show that the gauge-fixing term given by Eq.(34) can be obtained from the original BP Lagrangian by imposing the condition in the Green's function generating functional by means of the well-known Faddeev-Popov procedure [26]. This can be done introducing Eq.(44) via a Dirac delta functional in the integration measure. Explicitly we can write the generating functional as where N is a normalization constant and ∆ F P represents the determinant which arises from the Jacobian of a gauge transformation in the condition given by Eq.(44), that is, By means of introducing a pair of anticommuting ghost fields (C,C) and an auxiliary Nakanishi-Lautrup [27,28] field B, we can rewrite the generating functional, after a convenient redefinition of the normalization factor, as Note that the action in the exponential argument is invariant under the BRS transformation where the BRS operator δ B has Grassmann parity one. A functional integration over the B field finally leads to with S LGF given by Eq.(35), and thus justifying the Lorenz gauge fixing term given by Eq. (34). We have explicitly shown how the condition given by Eq.(44) leads through the Faddeev-Popov procedure to the gauge fixing term given by Eq.(34) and calculated the corresponding propagator for BP's generalized electrodynamics. In the next subsection, we turn our attention to the axial and light-front gauges.

Axial and Light-Front Gauges
The axial gauge fixing has been originally introduced by Kummer [24] and since then has been studied in the literature for a long time. It is a noncovariant gauge in the sense that it relies on a choice of an arbitrary fixed direction in space-time n µ . It encompasses the light-front gauge as a special case when n µ is light-like. For an interesting and lively review of the axial, light-front as well as other noncovariant gauges in the context of non Abelian theories we cite Ref. [25]. In the present case of BP model, in order to implement the axial gauge fixing we pick up a specific constant non-null four-vector n µ in space-time and write where α stands for a free gauge parameter. In ordinary electrodynamics the axial gauge fixing term contains only the first part of Eq.(52). The second part constitutes the natural generalization for BP's electrodynamics. According to the nature of the directional vector n µ we have different types of axial gauges. Namely temporal axial, light-front axial and space or proper axial for respectively timelike, lightlike or spacelike n µ . The differences among them lead to important subtleties and turn out to be a key point in the canonical quantization when one needs clearly to pick up a time direction for the Hamiltonian evolution. In our present discussion however, we limit ourselves to the Lagrangian analysis proposing Eq.(52) for a general axial gauge. Addition of the space-time integral of L LF to BP's action given by Eq.(1) results in the axial gauge fixed action which, after discarding surface integration terms can be recast into By performing a Fourier transformation, this action can be rewritten in momentum space as with The photon propagator of the BP model in the light-front gauge is the inverse of M µν , reading explicitly Consistently, for the case a = 0 this result reduces to the usual photon propagator in the axial gauge [23]. Note that despite the last term in Eq.(56), a term proportional to n µ n ν does not show up in the propagator P µν . For the light-front gauge, we have n 2 = 0 and by choosing the gauge parameter α = 0 we get the simpler expression [29] P µν = −1 Similarly to the Lorenz gauge, here we can justify the term given by Eq.(52) by imposing the gauge condition (1 + a 2 )n µ A µ = 0 (59) in the generating functional by means of the Faddeev-Popov procedure. In fact, we have here the Faddeev-Popov determinant which can be exponentiated by means of the introduction of a pair of anticommuting ghost fields (C,C) leading to the generating functional where B is the Nakanishi-Lautrup field. Paralleling the Lorenz gauge case, here we also have the BRS symmetry Integration over the Nakanishi-Lautrup field leads to an effective action in the exponential argument as showing clearly the appearance of the proposed term given by Eq.(52).
In the usual Maxwell case, there is a well-known discussion in the literature regarding the propagator of the gauge field in the light-front. Recently it has been shown [22] that it is possible to consider a mixing of the Lorenz and light-front gauge fixings leading to the three-term photon propagator [21]. In the following we show that there exists a natural generalization of the ideas discussed in Ref. [22] to the current BP model. Specifically, in order to obtain a three-term propagator, we define the axial Lorenz gauge-fixing Lagrangian density as where now we denote by β the gauge free parameter. Proceeding analogously to the previous sections, we integrate Eq.(64) in space-time and add the result to the gauge action given by Eq.(1). After disregarding frontier surface terms, the total gauge fixing action now reads By inverting the differential operator defined in Eq.(65) in momentum space, similarly to the previous cases, we obtain the generalized photon propagator of the BP model in the axial Lorenz gauge as Going back to the particular light-front gauge case where n 2 = 0 and choosing the gauge parameter β = 0 we get which is the corresponding three-term generalized photon propagator in the light-front gauge for the BP model. Our result here with the three-term propagator is precisely consistent with the result that we obtained using the interpolation between the instant form dynamics and the light-front dynamics for the electromagnetic gauge field [30] and for the QED [31]. As discussed in Refs. [30] and [31], the last term in Eq. (67) is canceled by the instantaneous interaction in the light-front dynamics so that the two-term gauge propagator given by Eq.(58) provides effectively the same result for the physical amplitude without involving the instantaneous interaction. As usual, the propagator of the BP model exhibits an additional simple pole at 1/a 2 . Note further that the propagator given by Eq.(67) satisfies known as the double transverse property [21] in the Maxwell case. Curiously enough, the gauge-fixing term given by Eq.(64) can be justified by imposing the two gauge conditions given by Eqs.(44) and (59) simultaneously. In fact, for arbitrary r, those two conditions together imply and from which we can write the generating functional as justifying the gauge-fixing term L AL in Eq.(64) (for the case r = 1).

Application -The Electron Self-Energy
In this section, we discuss the second-order correction to the electron self-energy considering the propagator of the BP model for the photon. We specifically calculate the invariant amplitude for the Feynman diagram represented in Fig.2 which has one loop integration in the internal momentum k. This effectively illustrates how BP formulation parallels with the PV regularization. The invariant amplitude for the diagram in Fig.2 is given by where u(p) denotes a plane wave solution to Dirac's equation and with the photon propagator of the BP model in the Lorenz gauge given by Eq.(43). In the Feynman-t'Hooft gauge choice, ξ = 1, and we may write more directlŷ with the notational definitions and D P k ≡ k 2 (1 − a 2 k 2 ) .
Actually, even for an arbitrary value of ξ, the calculation of the simpler expression (76) is sufficient to obtain the full invariant amplitude in Eq.(74) for the originalΣ(p) in Eq.(75).
We have pursued the gauge fixing of BP's generalized electromagnetism in three distinct ways. Specifically we analyzed the standard Lorenz, axial and light-front gauge fixings. For all considered cases, we have achieved a clean and neat generalized photon propagator depending on two parameters, namely the free gauge parameter and BP's length dimensional parameter. As a general result, we have shown the propagator of the BP model can have the same structure of the usual Maxwell case with only an additional pole at k 2 = 1/a 2 . Note that the common multiplicative factor in the propagator can be split as Thus, it can be interpreted as the sum of two propagators describing a massless photon and a massive one corresponding to the Pauli-Villars regularization. We have shown that the different gauge fixings considered do not invalidate this appealing interpretation.
The integration in θ from 0 to π gives Cancelling the 2b 2 common factor and regrouping terms we obtain which represents the result used in the main text to be substituted into equation (29).