Q uantum field theory can be understood through gauge theories. It is already established that the gauge theories can be studied either perturbativel y or non-perturbatively. Perturbative means using Feynman diagrams and non-perturbative means using Path-integral method. Operator regularization (OR) is one of the exceptional methods to study gauge theories because of its two-fold prescriptions. That means in OR two types of prescriptions have been introduced , which gives us the opportunity to check the result in self consistent way. In an earlier paper , we have evaluated basic QED loop diagrams in (3 + 1) dimensions using the both methods of OR and Dimensional regularization (DR). Then all three results have been compared. It is seen that the finite part of the result is almost same. In this paper , we are interested to evaluate the same basic loop diagrams in (2 + 1) space-time dimensions, because of two reasons: the main reason in (2 + 1) space-time dimensions , these loops diagrams are finite, on other hand , there are divergences in (3 + 1) space-time dimensions and the other reason is to see validity of using OR to evaluate Feynman loop diagrams in all dimensions. Here we have used both prescriptions of OR and DR to evaluate the basic loop diagrams and compared the results. Interestingly the results are almost same in all cases.

Operator Regularization Dimensional Regularization Feynman Diagrams in QED Path-Integral Method Background Field Quantization and Generating Functional
1. Introduction

Gauge theories  -  describe the interactions of all known forces such as electromagnetic, week and strong interations. To get a clear picture of different features from these theories, renormalization is a must. These can be done in two ways. That means one can study these theories either perturbatively or non-perturbatively. In the perturbative method, we mainly get some Feynman diagrams from the theory. When we evaluate loop diagrams in some cases, the loop integrals are divergent. So we need to regularize the integrals. That is why perturbative method needs regularization prescription. There are many regularization methods available in the literature. However, the most popular and appropriately applicable methods are Dimensional regularization (DR)   , Pauli-Villars regularization  , Pre-regularization     and some others. In non-perturbative method, we expand the generating function in terms of path integrals, then using different techniques we try to renormalize the theory to find the different features of the particles involved in the interactions. Operator regularization (OR) method is one of the best non-perturbative methods to study gauge theories. The remarkable feature of this method is that it gives us the opportunity to study the theory both perturbatively and non-perturbatively. That means although OR is mainly a path integral method but at one stage there is an option to consider the term as a factor for operator for loop diagrams. Then we can evaluate Feynman diagrams using these operators. Operator regularization method was prescribed by D.G.C. McKeon et al.    to study gauge theories non-perturbatively. However, they mentioned there that at one stage one can also use this prescription perturbatively. That means at that point one can also draw Feynman diagrams. A.Y. Sheiek   has showed how one can use Feynman diagrammatic technique in OR method.

In an earlier paper  we have described OR method in both ways and evaluated one-loop Feynman diagrams in QED in (3 + 1) space-time dimensions. Also we have used DR method in evaluating these diagrams and compared the results with the results obtained in OR in both ways. We have seen all the results are exactly same, except a finite constant term which will not affect the renormalization of the theory. In this paper, we have used the same method to evaluate the basic QED one-loop Feynman diagrams in (2 + 1) space-time dimensions to see the basic difference between finite and infinite loop integrals. Because in (3 + 1) dimensions, the loop integrals are divergent, on the other hand, in (2 + 1) dimensions, the loop integrals are finite.

2. Operator Regularization Prescription

Operator regularization is a convenient method of computing quantum corrections in quantum field theory in the context of background-field quantization and using path integral method    , which were given by D.G.C. McKeon et al. In this method the Feynman diagrams of the usual perturbation series can be avoided. But at one stage there is an option to consider the factor used for operators and inverse operators for Feynman diagrams. That is from this prescription one can choose either path integral method or Feynman diagrammatic approach.

We have clearly explained how this OR can be applied for evaluating Feynman diagrams in ref.  . For self consistence let us write a few main steps of this prescription which has to be used in evaluating the one-loop Feynman diagrams in (2 + 1)-dimensions.

In gauge theories we mainly deal with the generating functions. Then after some simplification we end up with some types operator and inverse operators. Then how one can take care of these operators has been explained in this OR.

If we have an operator Ω then according to OR we can write

det Ω = exp ( t r ln Ω ) (2.1)

Let us regularize ln Ω in the following way:

ln Ω = − lim ε → 0 d n d ε n ( ε n − 1 n ! Ω − ε )                 ( n = 1 , 2 , 3 , ⋯ ) (2.2)

In facing no divergences we can always choose n to be greater than or equal to the number of “loop momentum integrals” or in other words order in ℏ .

Hence,

det Ω = exp [ t r { − lim ε → 0 d n d ε n ( ε n − 1 n ! Ω − ε ) } ] (2.3a)

and

Ω − m = ( − 1 ) m − 1 ( m − 1 ) ! d m d Ω m ( ln Ω ) = − lim ε → 0 d n d ε n ( ε n − 1 n ! Γ ( ε + m ) Γ ( m ) Γ ( ε ) Ω − ε − m ) (2.3b)

If we now rewrite Ω − ε as

Ω − ε = 1 Γ ( ε ) ∫ 0 ∞ d t   t ε − 1 exp ( − Ω t ) . (2.4)

in Equation (2.9) we arrive at the result

det Ω = exp [ − ξ ′ ( 0 ) ] (2.5a)

where we have defined the ξ-function

ξ ( ε ) = 1 Γ ( ε ) ∫ 0 ∞ d t   t ε − 1 t r exp ( − Ω t ) (2.5b)

This is the usual ξ-function regularization of the determinant of an operator   .

Equations (2.2) and (2.3) are the main steps of the Operator regularization which has to be used in evaluating the Green’s function of any problem. From this point we can divide the prescription in two fold way. That means if we want to follow the path integral method we have to use Schwinger expansion  for the operator like

det Ω = exp { − lim ε → 0 d d ε [ 1 Γ ( ε ) ∫ 0 ∞ d t   t ε − 1 t r ( e − Ω 0 t − t e Ω 0 t Ω I t 2 2 ∫ 0 1 d u   e − ( 1 − u ) Ω 0 t Ω I e − u Ω 0 t Ω I     − t 3 3 ∫ 0 1 d u   u ∫ 0 1 d v   e − ( 1 − u ) Ω 0 t Ω I e − u ( 1 − v ) Ω 0 t Ω I e − u v Ω 0 t Ω I + ⋯ ) ] } (2.6)

where, Ω = Ω 0 + Ω I with Ω 0 is independent of the background field f i and Ω I is at least linear in f i .

Then following the steps described in ref.  we can find the result of the problems in consideration.

On the other hand if we want to use perturbatuve method then we have to take n = 1 for one -loop, take n = 2 for two-loops in Equations (2.2) and (2.3) and so on.

From Equation (2.3b) we can write the general prescription of Operator regularization for the Feynman diagrams following  :

Ω − m = lim ε → 0 d n d ε n [ ( 1 + α 1 ε + α 2 ε 2 + ⋯ + α n ε n ) ε n n ! Ω − ε − m ] (2.7)

where the α n s are arbitrary. For one-loop diagrams it is enough to use n = 1. When m = 2 and n = 1, then Equation (2.7) taken the form

Ω − 2 = lim ε → 0 d d ε [ ε ( 1 + α ε ) Ω − ε − 2 ] (2.8)

In one loop calculations we can use (2.8) for operators. In the following sections we will use this prescription for evaluating the three basic one loop diagrams.

2.1. One Loop Fermions Correction in (2 + 1) Dim. Using Dimensional and Operator Regularizations [<xref ref-type="bibr" rid="scirp.87995-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.87995-ref27">27</xref>] [<xref ref-type="bibr" rid="scirp.87995-ref28">28</xref>] [<xref ref-type="bibr" rid="scirp.87995-ref29">29</xref>]

1) Dimensional Regularization Method:

Starting with the Feynman diagram for the one loop correction to the fermions line shown in Figure 1 which is represented by ( ∑ ( p ) ) :

Using the standard Feynman rules one can write ( ∑ ( p ) ) as,

( ∑ ( p ) ) = − i e 2 ∫ d 3 l ( 2 π ) 3 γ μ ( p − l − m ) [ ( p − l ) 2 − m 2 ] l 2 γ μ

Using the Feynman identity for combining the denominators, we can write

( ∑ ( p ) ) = − i e 2 ∫ 0 1 d x ∫ d 3 l ( 2 π ) 3 γ μ ( p − l − m ) γ μ [ ( p − l ) 2 x − m 2 x + l 2 ( 1 − x ) ] 2

Shifting the variable of integration as l ′ = l − p x and simplifying we get

( ∑ ( p ) ) = − i e 2 ∫ 0 1 d x ∫ d 3 l ′ ( 2 π ) 3 γ μ [ p ( 1 − x ) − m − l ′ ] γ μ [ l ′ 2 − m 2 x + p 2 x ( 1 − x ) ] 2

The term linear in l ′ integrates to zero because of symmetric integration, so

( ∑ ( p ) ) = − i e 2 ∫ 0 1 d x ∫ d 3 l ( 2 π ) 3 γ μ [ p ( 1 − x ) − m ] γ μ [ l 2 − m 2 x + p 2 x ( 1 − x ) ] 2 (2.1.1)

Which is taken as the common starting point for both Dimensional and Operator regularization.

Using Feynman identity and then γ-algebra, the above result becomes,

( ∑ ( p ) ) = e 2 μ π ( 4 π ) 2 ( 2 π μ ) ∫ 0 1 d x [ p ( 1 − x ) + 3 m ] [ m 2 − p 2 x ( 1 − x ) ] 1 / 2 = e 2 μ 2 ( 8 π ) ∫ 0 1 d x [ p ( 1 − x ) + 3 m ] [ m 2 − p 2 x ( 1 − x ) ] 1 / 2 (2.1.2)

Thus according to dimensional regularization, we see that there is no divergent part in (2 + 1)-dimensional space-time, because the integrals are finite in 3-dimensions.

2) Operator Regularization Method:

The same one-loop correction to fermion can be evaluated using OR, following the rule cited in Equations (2.5) and (2.6) in ref.  . The amplitude of the self-energy diagram as

( ∑ ( p ) ) = − i e 2 ∫ 0 1 d x lim ε → 0 d d ε ∫ d 3 l ( 2 π ) 3 ε ( 1 + α ε ) γ μ [ p ( 1 − x ) − m ] γ μ [ l 2 − m 2 + p 2 x ( 1 − x ) ] ε + 2

Using the standard integral

∫ d 2 w l ( 2 π ) 2 w 1 ( l 2 + M 2 ) A = 1 ( 4 π ) w Γ ( A ) Γ ( A − w ) ( M 2 ) A − w (2.1.3)

we get,

( ∑ ( p ) ) = − i e 2 ( 4 π ) 3 / 2 ∫ 0 1 d x   γ μ [ p ( 1 − x ) − m ] γ μ     ⋅ lim ε → 0 d d ε [ ε ( 1 + α ε ) Γ ( ε + 2 ) . Γ ( ε + 1 2 ) ( − m 2 x + p 2 x ( 1 − x ) ) ε + 1 2 ]

Here, lim ε → 0 d d ε [ ε ( 1 + α ε ) Γ ( ε + 2 ) . Γ ( ε + 1 2 ) ( − m 2 x + p 2 x ( 1 − x ) ) ε + 1 2 ] = lim ε → 0 d d ε [ ε ( 1 + α ε ) ( ε − 1 2 ) ! u − ( ε + 1 2 ) ( ε + 1 ) ! ]

where u = − m 2 x + p 2 x ( 1 − x )     = lim ε → 0 [ u − ( ε + 1 2 ) ( ε α ) ( − 1 2 + ε ) ( 1 + ε ) ! + u − ( ε + 1 2 ) ( 1 + α ε ) ( ε − 1 2 ) ! ( 1 + ε ) !         − u − ( ε + 1 2 ) ε ( 1 + α ε ) ( ε − 1 2 ) ! ln u ( 1 + ε ) ! + u − ( ε + 1 2 ) ε ( 1 + ε α ) Γ ( ε + 1 2 ) F ( ε + 1 2 ) ( 1 + ε ) !         − u − ( ε + 1 2 ) ε ( 1 + ε α ) ( ε − 1 2 ) ! Γ ( ε + 2 ) F ( ε + 2 ) ( 1 + ε ) ! ] = π u 1 / 2

Therefore,

lim ε → 0 d d ε [ ε ( 1 + α ε ) Γ ( ε + 2 ) . Γ ( ε + 1 2 ) ( − m 2 x + p 2 x ( 1 − x ) ) ε + 1 2 ] = π [ p 2 x ( 1 − x ) − m 2 x ] 1 / 2 (2.1.4)

Thus Equation (2.1.4) becomes,

( ∑ ( p ) ) = − i e 2 ( 4 π ) 3 / 2 ∫ 0 1 d x [ p ( 1 − x ) + 3 m ] π [ p 2 x ( 1 − x ) − m 2 x ] 1 / 2 = − e 2 ( 8 π ) ∫ 0 1 d x [ p ( 1 − x ) + 3 m ] [ m 2 − p 2 x ( 1 − x ) ] 1 / 2

This is the same form as like as obtained by dimensional regularization approach.

2.2. One Loop Photon Correction in (2 + 1) Dim. Using Dimensional and Operator Regularizations

Let us consider the Feynman diagram for the one loop correction to the photon line shown in Figure 2 which is represented by Π μ ν ( p ) :

The QED one loop correction to the photon line in (2 + 1)-dimensions is

Π μ ν ( p ) = e 2 ∫ d 3 l ( 2 π ) 3 T r [ γ μ ( l + p − m ) γ ν ( l − m ) [ ( l + p ) 2 − m 2 ] ( l 2 − m 2 ) ]

Combining the denominator using the Feynman identity and simplifying, we get

Π μ ν ( p ) = e 2 ∫ 0 1 d x ∫ d 3 l ( 2 π ) 3 T r [ γ μ ( l + p − m ) γ ν ( l − m ) ] [ ( l + p ) 2 x − m 2 x + ( l 2 − m 2 ) ( 1 − x ) ] 2 (2.2.1)

Now putting l ′ = l + p x in Equation (2.2.1), then we get,

Π μ ν ( p ) = e 2 ∫ 0 1 d x ∫ d 3 l ′ ( 2 π ) 3 T r [ γ μ ( l ′ − p x + p − m ) γ ν ( l ′ − p x − m ) ] [ ( l ′ + p ( 1 − x ) ) 2 x − m 2 x + ( ( l ′ − p x ) 2 − m 2 ) ( 1 − x ) ] 2

After simplification Equation (2.2.1) with l ′ → l becomes,

Π μ ν ( p ) = 4 e 2 ∫ 0 1 d x ∫ d 3 l ( 2 π ) 3 [ 2 l μ l ν { l 2 − m 2 + p 2 x ( 1 − x ) } 2     − 2 x ( 1 − x ) ( p μ p ν − p 2 δ μ ν ) { l 2 − m 2 + p 2 x ( 1 − x ) } 2 − δ μ ν { l 2 − m 2 + p 2 x ( 1 − x ) } ] (2.2.2)

If we apply the following integrals in the first and third terms in the integrand of Equation (2.2.2),

I) ∫ d d l l μ l ν ( l 2 + 2 l q − m 2 ) α = i π d / 2 α 1 ( − q 2 − m 2 ) α − d / 2 [ q μ q ν Γ α − d 2 + 1 2 g μ ν ( − q 2 − m 2 ) Γ ( α − 1 − d 2 ) ]

II) ∫ d d l 1 ( l 2 + 2 l q − m 2 ) α = ( − 1 ) α i π d / 2 Γ ( α − d / 2 ) Γ ( α ) ( − q 2 − m 2 ) α − d / 2

We arrive at,

Π μ ν ( p ) = 8 e 2 ( p μ p ν − p 2 δ μ ν ) ∫ 0 1 d x ∫ d 3 l ( 2 π ) 3 [ x ( 1 − x ) { l 2 − m 2 + p 2 x ( 1 − x ) } 2 ] (2.2.3)

Which is again taken as the common starting point for both Dimensional and Operator regularization for one loop correction to the photon lines.

Using Feynman identity-II, the above result becomes,

Π μ ν ( p ) = − 8 i e 2 ( p μ p ν − δ μ ν p 2 ) ( 2 π ) 3 / 2 ∫ 0 1 d x   x ( 1 − x ) [ p 2 x ( 1 − x ) − m 2 2 π μ 2 ] − ε / 2 = − 8 i e 2 μ ( p μ p ν − δ μ ν p 2 ) ( 2 π ) ∫ 0 1 d x x ( 1 − x ) [ p 2 x ( 1 − x ) − m 2 ] 1 / 2 (2.2.4)

Thus according to dimensional regularization, we see that there is no divergent part in (2 + 1)-dimensional space-time, because the integrals are finite in 3-dimensions.

Now proceeding with operator regularization and again following the same route, we get,

Π μ ν ( p ) = − 8 e 2 ( p μ p ν − p 2 δ μ ν ) ∫ 0 1 d x lim ε → 0 d d ε ∫ d 3 l ( 2 π ) 3 ε ( 1 + α ε ) x ( 1 − x ) [ l 2 − m 2 + p 2 x ( 1 − x ) ] ε + 2 (2.2.5)

Performing the momentum integral (2.1.3), so that relative to the Equation (2.2.5) we get, A = ε + 2 , w = 3 2 , M 2 = − m 2 + p 2 x ( 1 − x ) , then Equation (2.2.5) becomes,

Π μ ν ( p ) = − 8 e 2 1 ( 4 π ) 3 / 2 ( p μ p ν − p 2 δ μ ν ) ∫ 0 1 d x   x ( 1 − x )     ⋅ lim ε → 0 d d ε [ ε ( 1 + α ε ) Γ ( ε + 2 ) . Γ ( ε + 1 2 ) ( − m 2 + p 2 x ( 1 − x ) ) ε + 1 2 ]

Now using the Equation (2.1.4) then above equation reduces to,

Π μ ν ( p ) = − 8 e 2 1 ( 4 π ) 3 / 2 ( p μ p ν − p 2 δ μ ν ) x ( 1 − x ) π [ p 2 x ( 1 − x ) − m 2 x ] 1 / 2 = − 4 e 2 ( p μ p ν − δ μ ν p 2 ) ( 2 π ) ∫ 0 1 d x x ( 1 − x ) [ p 2 x ( 1 − x ) − m 2 ] 1 / 2

This is the same form as we obtained by dimensional regularization approach.

2.3. One Loop Vertex Correction in (2 + 1) Dim. Using Dimensional and Operator Regularizations

Let us now consider the Feynman diagram for the one loop correction to the vertex shown in Figure 3 which is represented by Γ ρ ( p , q ) .

The QED one loop correction to the vertex in (2 + 1)-dimensions is

Γ ρ ( p , q ) = ∫ d 3 l ( 2 π ) 3 [ ( − i e γ λ ) i ( p + l + m ) ( − i e γ ρ ) i q + l + m ( − i e γ σ ) δ τ σ l 2 ] = − i e 3 ∫ d 3 l ( 2 π ) 3 [ γ λ ( p + l − m ) γ ρ ( q + l − m ) γ λ ] l 2 [ ( p + l ) 2 − m 2 ] [ ( q + l ) 2 − m 2 ] (2.3.1)

Applying the 3-parameter Feynman formula for combining the denominator

and shifting the variable of integration l → l + p x + q y and simplify the denominator and numerator, we obtain,

Γ ρ ( p , q ) = − 2 i e 3 ∫ 0 1 d x ∫ 0 1 − x d y ∫ d 3 l ( 2 π ) 3     ⋅ [ γ λ ( l + p ( 1 + x ) + q y − m ) γ ρ ( l + q ( 1 + y ) + p x − m ) γ λ ] [ l 2 − m 2 ( x + y ) + p 2 x ( 1 − x ) + q 2 y ( 1 − y ) − 2 p ⋅ q x y ] 3

This integral contains convergent and divergent pieces. The part of the numerator quadratic in l is divergent, the rest convergent, so separating the divergent piece Γ ρ ( 1 ) ( p , q ) and convergent piece Γ ρ ( 2 ) ( p , q ) , i.e.

Γ ρ ( p , q ) = Γ ρ ( 1 ) ( p , q ) + Γ ρ ( 2 ) ( p , q )

Thus the divergent piece is,

Γ ρ ( 1 ) ( p , q ) = − 2 i e 3 ∫ 0 1 d x ∫ 0 1 − x d y ∫ d 3 l ( 2 π ) 3 γ σ l γ ρ l γ σ ( l 2 − M 2 ) 3 (2.3.2)

where, M 2 ≡ m 2 ( x + y ) − p 2 x ( 1 − x ) − q 2 y ( 1 − y ) + 2 p ⋅ q x y .

Which is taken as the common starting point for both Dimensional and Operator regularizations for one-loop correction to the vertex.

Using Feynman identity and then γ-algebra, the above result becomes,

Γ ρ ( 1 ) ( p , q ) = 2 e 3 μ Γ ( 1 2 ) ( 4 π ) 2 ∫ 0 1 d x ∫ 0 1 − x d y γ ρ ( M 2 ) 1 / 2 = 2 e 3 μ π γ ρ ( 4 π ) 2 ∫ 0 1 d x ∫ 0 1 − x d y 1 ( M 2 ) 1 / 2 = e 3 μ γ ρ ( 4 π ) ∫ 0 1 d x ∫ 0 1 − x d y 1 ( M 2 ) 1 / 2 (2.3.3)

where, M 2 ≡ m 2 ( x + y ) − p 2 x ( 1 − x ) − q 2 y ( 1 − y ) + 2 p ⋅ q x y .

Now proceeding with operator regularization through the rule (2.13) and (2.14), we get,

Γ ρ ( 1 ) ( p , q ) = − 2 i e 3 ∫ 0 1 d x ∫ 0 1 − x d y lim ε → 0 d d ε ∫ d 3 l ( 2 π ) 3     ε ( 1 + α ε ) γ σ l γ ρ l γ σ ( l 2 − M 2 ) ε + 3

Now performing the momentum integral -I from above, then we get,

Γ ρ ( 1 ) ( p , q ) = − i e e 2 ( 4 π ) 3 / 2 ∫ 0 1 d x ∫ 0 1 − x d y lim ε → 0 d d ε [ ε ( 1 + α ε ) Γ ( ε + 3 ) ⋅ Γ ( ε + 1 2 ) ( M 2 ) ε + 1 2 ] γ σ γ λ γ ρ γ τ γ σ

Now using the Equation (2.1.4) and γ-algebra, then above equation reduces to,

Γ ρ ( 1 ) ( p , q ) = − i e e 2 ( 4 π ) 3 / 2 ∫ 0 1 d x ∫ 0 1 − x d y π 2 [ M 2 ] 1 / 2 γ ρ = − i e e 2 γ ρ ( 16 π ) ∫ 0 1 d x ∫ 0 1 − x d y 1 [ M 2 ] 1 / 2 (2.3.4)

where, M 2 ≡ m 2 ( x + y ) − p 2 x ( 1 − x ) − q 2 y ( 1 − y ) + 2 p ⋅ q x y .

This is the same form as like as obtained by dimensional regularization approach.

3. Path Integral Form of Operator Regularization for One Loop Generating Functional in QED

The path integral form of OR for one-loop case is described in ref.  . That is if we consider the QED Lagrangian as,

L = − 1 4 ( ∂ μ Ω μ − ∂ ν Ω ν ) 2 + Ψ ¯ ( − i γ μ ∂ μ − e γ μ Ω μ − m ) Ψ − 1 2 α ( ∂ μ Ω μ ) 2 (3.1)

and let us expand this Lagrangian taking background field quantization of the fields in the following form, gauge field Ω μ and fermionic field ψ are respectively,

Ω μ = V μ + Q μ , for gauge field Ω μ

ψ = η + q , for fermionic field ψ ,

where V μ and η are the classical fields and Q μ and q are the quantum fields.

Therefore Equation (3.1) becomes,

L = − 1 4 ( ∂ μ V μ + ∂ μ Q μ − ∂ ν V ν − ∂ ν Q ν ) 2     + ( η ¯ + q ¯ ) ( − i γ μ ∂ μ − e γ μ ( V μ + Q μ ) − m ) ( η + q )     − 1 2 α ( ∂ μ V μ + ∂ μ Q μ ) 2

= 1 4 [ ( ∂ μ Q μ − ∂ ν Q ν ) + ( ∂ μ V μ − ∂ ν V ν ) ] 2 + η ¯ ( − i γ μ ∂ μ − m ) η         + η ¯ ( − i γ μ ∂ μ − m ) q + q ¯ ( − i γ μ ∂ μ − m ) η         + e ( η ¯ γ μ V μ η + η ¯ γ μ V μ q + η ¯ γ μ Q μ η + η ¯ γ μ Q μ q + q ¯ γ μ V μ η         + q ¯ γ μ V μ q + q ¯ γ μ Q μ η + q ¯ γ μ Q μ q ) − 1 2 α ( ∂ μ V μ + ∂ μ Q μ ) 2 (3.2)

Following ref.  the one-loop generating functional for Green’s functions Z 1 is

Z 1 = det D det 1 / 2 [ p 2 δ μ ν − ( 1 − 1 α ) p μ p ν − e 2 η ¯ γ μ D − 1 γ ν η − e 2 η ¯ γ ν D − 1 γ μ η ] = det A det B (say) (3.3)

where, D = γ μ ( − i ∂ μ − e V μ ) = p − e V

Here we see that Z 1 is the ratio of determinant of operators. Each of the determinants occurring in Equation (3.3) requires regularization and a corresponding ξ-function. The numerator and denominator separately contribute to Green’s functions with only external boson lines and with both external fermions lines and vertex function in massless QED respectively.

3.1. One-Loop Generating Functional and Loop Corrections for External Boson Lines

To find the loop corrections or to write the generating functional for external boson lines one has to make a close look at the numerator of Equation (3.3) and on the other hand for external fermion lines one has to take care of the denominator of Equation (3.3). So for bosonic case we have to regulate the det A through the use of ξ-function in Equation (2.5a) yielding

det A = Z 1 A = exp [ − lim ε → 0 ξ / A ( ε ) ] (3.1.1)

where,

ξ ( ε ) = 1 Γ ( ε ) ∫ 0 ∞ d t   t ε − 1 t r exp ( − Ω t ) with Ω = p − e V . (3.1.2)

As we mentioned in Section-2, after regularization we have to consider Schwinger expansion, to this view let us now identify the operator Ω 0 and Ω I with p and − e V , respectively, then by Equation (2.4), Equation (3.1.1) can be written as,

Z 1 A = exp [ − lim ε → 0 d d ε { 1 Γ ( ε ) ∫ 0 ∞ d t   t ε − 1 t r ( e − p t − t e − p t + ( − e V ) t 2 2       ⋅ ∫ 0 1 d u   e − ( 1 − u ) p t ( − e V ) e − u p t ( − e V )       − t 3 3 ∫ 0 1 d u   u ∫ 0 1 d v   e − ( 1 − u ) p t ( − e V ) e − u ( 1 − v ) p t ( − e V ) e − u v p t ( − e V ) + ⋯ ) } ] (3.1.3)

To one-loop order this series plays the same role as Feynman rules in the usual perturbation theory. Here we want to evaluate the one-loop correction to the two-point function for external photon in QED; we restrict our attention to the term bilinear in V μ on the right-hand side of Equation (3.1.3). This leaves us with

Z 1 V V A = exp [ − lim ε → 0 d d ε { e 2 Γ ( ε ) ∫ 0 ∞ d t t ε + 1 2     t r ( ∫ 0 1 d u   e − ( 1 − u ) p t V e − u p t V ) } ] (3.1.4)

Now let us complete the functional trace

T = t r ( ∫ 0 1 d u   e − ( 1 − u ) p t V e − u p t V ) (3.1.5)

Schwinger has pointed out that such traces are most easily evaluated in momentum space. We introduce a complete orthonormal set of states | p 〉 that are eigenstates of the operator p μ , where, in n dimensions,

〈 x | p 〉 = e i p ⋅ x ( 2 π ) n / 2 (3.1.6a)

and

〈 p | f | q 〉 = f ( p − q ) ( 2 π ) n / 2 (3.1.6b)

On the right-hand side of Equation (3.1.6b), f ( p − q ) is the Fourier transform of f ( x ) :

f ( p − q ) = ∫ d n x ( 2 π ) n / 2 f ( x ) e − i x ⋅ ( p − q ) (3.1.7)

Equation (3.1.5) takes the form,

T = ∫ d 3 p d 3 q d 3 r d 3 s 〈 p | e − ( 1 − u ) p t | q 〉 〈 q | V | r 〉 〈 r | e − u p t | s 〉 〈 s | V | p 〉 (3.1.8)

Upon inserting the complete set 1 = ∫ d 3 p | p 〉 〈 p | at the appropriate places, and using (3.1.6), we rewrite Equation (3.1.8) as,

T = ∫ d 3 p d 3 q d 3 r d 3 s V ( p − q ) ( 2 π ) 3 / 2 e − ( 1 − u ) p t ⋅ e i r . q ( 2 π ) 3 / 2 δ ( r − q )     ⋅ V ( r − s ) ( 2 π ) 3 / 2 e − u r t ⋅ e i s ⋅ p ( 2 π ) 3 / 2 δ ( s − p ) = ∫ d 3 p d 3 q ( 2 π ) 3     e − ( 1 − u ) p t − u q t V ( p − q ) V ( q − p ) (3.1.9)

After shifting the variable of integration p → p + q , Equation (3.1.9) becomes,

T = ∫ d 3 p d 3 q ( 2 π ) 3     e − [ q + ( 1 − u ) p ] t V ( p ) V ( − p ) (3.1.10)

Upon substituting Equation (3.1.10) into Equation (3.1.4), we find that

Z 1 V V A = exp [ − lim ε → 0 ξ ′ V V A ( ε ) ] (3.1.11a)

where,

ξ V V A ( ε ) = e 2 2 Γ ( ε ) ∫ 0 ∞ d t   t ε + 1 ∫ 0 1 d u ∫ d 3 p V ( p ) V ( − p ) ∫ d 3 q ( 2 π ) 3 e − [ q + ( 1 − u ) p ] t (3.1.11b)

We use Equation (2.4) to integrate over t, then (3.1.11b) becomes,

ξ V V A ( ε ) = e 2 Γ ( ε + 2 ) 2 Γ ( ε ) ∫ d 3 p V ( p ) V ( − p ) ∫ 0 1 d u ∫ d 3 q ( 2 π ) 3 [ q − ( 1 − u ) p ] ε + 2 [ q 2 − ( 1 − u ) 2 p 2 ] ε + 2 (3.1.12)

Now the last integral I 1 (say) of Equation (3.1.12) can be calculated as follows:

I 1 = ∫ d 3 q ( 2 π ) 3 [ q − ( 1 − u ) p ] ε + 2 [ q 2 − ( 1 − u ) 2 p 2 ] ε + 2 = ∫ d 3 q ( 2 π ) 3 [ q ∈ + 2 − ( ε + 2 ) ( q ) ε + 1 ( 1 − u ) p + ⋯ + ( 1 − u ) ε + 2 ( p ) ε + 2 ] [ q 2 − ( 1 − u ) 2 p 2 ] ε + 2

Differentiating Equation (3.1.12) with respect to ε and taking ε → 0 , we see that the product terms in ε will vanish. Hence in the numerator of I 1 only the first and last term will contribute.

∴ I 1 = ∫ d 3 q ( 2 π ) 3 [ ( q 2 ) ε + 2 2 ( q 2 − ( 1 − u ) 2 p 2 ) ε + 2 + ( 1 − u ) ε + 2 ( p ) ε + 2 ( q 2 − ( 1 − u ) 2 p 2 ) ε + 2 ] (3.1.13)

To evaluate this integral let us consider the standard integral,

∫ d n q ( 2 π ) n ( q 2 ) r ( q 2 + c 2 ) m = 1 ( 16 π 2 ) n / 4 ( c 2 ) ( n 2 ) + r − m Γ ( r + n 2 ) Γ ( m − r − n 2 ) Γ ( n 2 ) Γ ( m ) (3.1.14)

Using Equation (3.1.14) in Equation (3.1.13), we get,

I 1 = 1 ( 16 π 2 ) 3 / 4 [ { − ( 1 − u ) 2 p 2 } ( 3 2 + ε + 2 2 − ( ε + 2 ) ) Γ ( ε + 2 2 + 3 2 ) Γ ( ε + 2 − ε + 2 2 − 3 2 ) Γ ( 3 2 ) Γ ( ε + 2 )       + ( 1 − u ) ε + 2 ( p ) ε + 2 { − ( 1 − u ) 2 p 2 } ( 3 2 + 0 − ε − 2 ) ⋅ Γ ( 3 2 ) Γ ( ε + 2 − 3 2 ) Γ ( 3 2 ) Γ ( ε + 2 ) ] = 1 ( 4 π ) 3 / 2 1 Γ ( ε + 2 ) [ ( − 1 ) 1 − ε 2 ( 1 − u ) 1 − ε ( p ) 1 − ε 2 π Γ ( ε 2 + 5 2 ) Γ ( ε 2 − 1 2 )       + ( − 1 ) − ( ε + 1 2 ) ( 1 − u ) 1 − ε ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) Γ ( ε + 1 2 ) ]

Thus Equation (3.1.12) becomes,

ξ V V A ( ε ) = e 2 Γ ( ε + 2 ) 2 Γ ( ε ) 1 ( 4 π ) 3 / 2 1 Γ ( ε + 2 ) ∫ d 3 p V ( p ) V ( − p ) ∫ 0 1 d u [ ( − 1 ) 1 2 − ε 2     ⋅ ( 1 − u ) 1 − ε ( p ) 1 − ε 2 π Γ ( ε 2 + 5 2 ) Γ ( ε 2 − 1 2 )     + ( − 1 ) − ( ε + 1 2 ) ( 1 − u ) 1 − ε ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) Γ ( ε + 1 2 ) ]

= e 2 8 π 2 1 Γ ( ε ) ∫ d 3 p V ( p ) V ( − p ) [ ( − 1 ) 1 2 − ε 2 ( 1 2 − ε )     ⋅ ( p ) 1 − ε ( ε 2 + 3 2 ) ( ε 2 + 1 2 ) ( ε 2 − 1 2 ) { Γ ( ε 2 − 1 2 ) } 2     + π 2 ( − 1 ) − ( ε + 1 2 ) ( 1 2 − ε ) ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) Γ ( ε + 1 2 ) ]

= e 2 8 π 2 ∫ d 3 p V ( p ) V ( − p ) [ ( − 1 ) 1 2 − ε 2 ( 1 2 − ε )       ⋅ ( p ) 1 − ε ( ε 2 + 3 2 ) ( ε 2 + 1 2 ) ( ε 2 − 1 2 ) Γ ( ε ) 4 1 Γ ( ε )       + π 2 ( − 1 ) − ( ε 2 + 1 2 ) ( 1 2 − ε ) ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) Γ ( ε 2 + 1 2 ) 1 Γ ( ε ) ]

where we have used [ Γ ( ε − 1 2 ) ] 2 = Γ ( ε ) 4 .

= e 2 8 π 2 ∫ d 3 p V ( p ) V ( − p ) [ ( − 1 ) 1 2 − ε 2 ( p ) 1 − ε ( ε 3 + 3 ε 2 − ε − 3 ) 24 ( 2 − ε )       + π 2 ( − 1 ) − ( ε + 1 2 ) ( 1 2 − ε ) ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) 2 2 ε − 1 ( Γ ( ε + 1 2 ) ) 2 π Γ ( 2 ε ) ]

[Using Duplication formula]

= e 2 8 π 2 ∫ d 3 p V ( p ) V ( − p ) [ ( − 1 ) 1 2 − ε 2 ( p ) 1 − ε ( ε 3 + 3 ε 2 − ε − 3 ) 24 ( 2 − ε )       + ( − 1 ) − ( ε + 1 2 ) ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) 2 2 ε − 1 ( 2 ε + 1 ) ε ( 2 − ε ) ]

where we have used [ Γ ( 2 ε + 1 2 ) ] 2 = ( 2 ε + 1 4 ) Γ ( 2 ε + 1 ) .

To see sign of the term, let us put ε → 0 in the factor ( − 1 ) 1 2 − ε 2 and ( − 1 ) − ( ε + 1 2 ) , then these terms are equal to i and ?i. Thus the above equation becomes,

ξ V V A ( ε ) = i e 2 8 π 2 ∫ d 3 p V ( p ) V ( − p ) [ ( p ) 1 − ε ( ε 3 + 3 ε 2 − ε − 3 ) 24 ( 2 − ε )     − ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) 2 2 ε − 1 ( 2 ε + 1 ) ε ( 2 − ε ) ] (3.1.15)

Now differentiating Equation (3.1.15) w. r. to ε , we get

ξ ′ V V A ( ε ) = i e 2 8 π 2 ∫ d 3 p V ( p ) V ( − p ) [ 1 24 ( p ) 1 − ε ln ( p ) ( − 1 ) ( ε 3 + 3 ε 2 − ε − 3 ) ( 2 − ε )     + 1 24 ( p ) 1 − ε { ( 2 − ε ) ( 3 ε 2 + 6 ε − 1 ) − ( ε 3 + 3 ε 2 − ε − 3 ) ( − 1 ) ( 2 − ε ) 2 }     + ( p ) ε + 2 ln ( p ) ( p ) − ( 2 ε + 1 ) 2 2 ε − 1 ( 2 ε + 1 ) ε ( 2 − ε ) ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) ln ( p )     ⋅ ( − 2 ) 2 2 ε − 1 ( 2 ε + 1 ) ε ( 2 − ε ) + ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) 2 2 ε − 1 ln ( 2 ) ( 2 ) ( 2 ε + 1 ) ε ( 2 − ε )     + ( p ) ε + 2 ( p ) − ( 2 ε + 1 ) 2 2 ε − 1 ( 2 − ε ) ( 4 ε + 1 ) − ( 2 ε 2 + ε ) ( − 1 ) ( 2 − ε ) 2 ]

Hence

ξ ′ V V A ( 0 ) = lim ε → 0 ξ ′ V V A ( ε ) = i e 2 8 π 2 ∫ d 3 p V ( p ) V ( − p ) [ 1 24 p ⋅ ln ( p ) ( − 1 ) ( − 3 2 )     + 1 24 p ⋅ ( − 5 4 ) + ( p ) 2 ⋅ 1 p ⋅ 1 2 ⋅ 1 2 ] = i e 2 128 π 2 ∫ d 3 p V ( p ) V ( − p ) p [ ln ( p ) − 29 6 ] (3.1.16)

Substituting of Equation (3.1.16) into Equation (3.1.11a) yields our final expression for Z 1 V V A as,

Z 1 V V A = exp [ i e 2 128 π 2 ∫ d 3 p V ( p ) V ( − p ) p ( ln ( p ) − 29 6 ) ] (3.1.17)

This contributes to the to the one?loop generating functional for external bosons (photon) lines.

To find one-loop correction for external boson lines from above generating functional, we have to take logarithm on Equation (3.1.17) and then functional differentiation of the expansion with respect to momentum p.

Thus the one-loop correction for the external boson lines is,

= i e 2 128 π 2 p ( ln ( p ) − 29 6 ) (3.1.18)

Hence the result in (3.1.18) is finite and of the same form as we obtained by the diagrammatic form of Operator regularization and Dimensional regularization methods in Section-2. In this section we have shown and explained how one can choose the appropriate terms from the Schwinger expansion for the problem in hand.

3.2. One-Loop Generating Functional and Loop Corrections for External Fermion Lines and Vertex Function

In this case we focus on the denominator in Equation (3.3), so that let us regulate the det B through use of the ξ-function in Equation (3.11a) yielding

det B = Z 1 B = exp [ − 1 2 lim ε → 0 ξ ′ B ( ε ) ] , (3.2.1)

where,

ξ B ( ε ) = 1 Γ ( ε ) ∫ 0 ∞ d t   t ε − 1   t r exp { − t [ p 2 δ μ ν − ( 1 − 1 α ) p μ p ν     − e 2 η ¯ γ μ 1 p − e V γ ν η − e 2 η ¯ γ ν 1 p − e V γ μ η ] } (3.2.2)

In Equation (3.2.2) it is understood that the exponential is t r [ exp ( − t B ) ] , where

B μ ν ≡ p 2 δ μ ν − ( 1 − 1 α ) p μ p ν − e 2 η ¯ γ μ 1 p − e V γ ν η − e 2 η ¯ γ ν 1 p − e V γ μ η ≡ B 0 μ ν + B I μ ν (3.2.3)

where, B 0 μ ν is independent of the background field η and η ¯ , and B I μ ν is at least linear in η and η ¯ .

Now as before to use Schwinger expansion in this case let us use the Equation (2.4) and then taking bilinear in η and η ¯ on the on the right-hand side of Equation (3.2.1), we end up with

Z 1 B = exp [ − 1 2 { lim ε → 0 d d ε ( 1 Γ ( ε ) ∫ 0 ∞ d t   t ε   t r { exp ( p 2 δ μ ν − ( 1 − 1 α ) p μ p ν )             ⋅ t ( − 2 e 2 η ¯ γ μ D − 1 γ ν η ) } } ] (3.2.4)

The exponential factor in the trace of Equation (3.2.4) can be simplified using the complete set of orthonormal projections operators:

T μ ν ( p ) = ( δ μ ν − p μ p ν p 2 ) (3.2.5a)

L μ ν ( p ) = p μ p ν p 2 (3.2.5b)

These allows us to write ( e − t B 0 ) μ ν as

exp ( p 2 δ μ ν − ( 1 − 1 α ) p μ p ν ) t = ∑ n = 0 ∞ 1 n ! [ p 2 t ( T μ ν + 1 α L μ ν ) ] n = e − t p 2 T μ ν + e − t p 2 / α L μ ν (3.2.6)

and let us expand D − 1 in powers of the back-ground field in the ξ-function Equation (3.2.4):

D − 1 = 1 p − e V = 1 p ( 1 − e V 1 p ) − 1 = ( 1 p + 1 p e V 1 p + 1 p e V 1 p e V 1 p + ⋯ )

It is interesting to note that at this stage this is straightforward to apply the perturbative expansion of Equation (2.3) to this ξ-function and to select from the expansion those terms appropriate for any particular Greens function. This means that from the expansion we can choose appropriate terms that are associated with the related problems that we are interested in. Let us consider here the ξ-function for the fermion two-point function and the vertex function, we find

ξ B ( ε ) ≅ e 2 Γ ( ε ) ∫ 0 ∞ d t   t ε   t r [ η ¯ γ μ ( 1 p + 1 p e V 1 p ) γ ν η ( e − t p 2 T μ ν + e − t p 2 / α L μ ν ) ] (3.2.7)

Following ref.   in the Approach-A, we compute from Equation (3.2.7) the ξ-function in the limit of zero momentum transfer to the photon:

ξ B ( ε ) = e 2 8 π 2 ∫ d 3 p ( 1 p 2 ) Γ ( 1 − ε ) Γ ( 3 − ε ) [ 2 + ( α ε + 1 − 1 ) ( 2 − ε ) ( 1 + ε ) ]     ⋅ η ¯ ( p ) [ p + e 4 π 2 ( V ( 0 ) − 2 ε p p 2 p ⋅ V ( 0 ) ) ] η ( − p ) (3.2.8)

Therefore by Equation (3.2.1) the contributions to the one-loop generating functional is

Z 1 B ≅ exp { 1 2 [ e 2 8 π 2 ∫ d 3 p ( 3 2 + α ln α − α ln p 2 ) η ¯ ( p ) ( p + e 4 π 2 V ( 0 ) ) η ( − p )     − α e 3 16 π 4 ∫ d 3 p η ¯ ( p ) p η ( − p ) p ⋅ V ( 0 ) p 2 ] } (3.2.9)

This contributes to the one-loop generating functional for external fermion (electron) lines and vertex function in QED.

To obtain the one-loop correction for external fermion lines and vertex function, we have to take logarithm of Equation (3.2.9) and then functional differentiation with respect to momentum p.

Hence from Equation (3.2.9), we get

= e 2 16 π 2 ( 3 2 + α ln α − α ln p 2 ) η ¯ ( p ) p η ( − p )       + e 3 64 π 4 ( 3 2 + α ln α − α ln p 2 ) η ¯ ( p ) V ( 0 ) η ( − p )       − α e 3 32 π 4 η ¯ ( p ) γ μ p μ η ( − p ) p ν V ν ( 0 ) p 2 = e 2 16 π 2 ( 3 2 + α ln α − α ln p 2 ) η ¯ ( p ) p η ( − p )       + e 3 64 π 4 ( 3 2 + α ln α − α ln p 2 − 2 α ) η ¯ ( p ) V ( 0 ) η ( − p ) (3.2.10)

From the expansion (3.2.10) we can find the one-loop correction for the external fermion lines and one-loop vertex function by choosing the appropriate terms. This expression is of the same form as obtained by DR and OR methods with Feynman diagrams in Section-2.

Thus the one-loop correction for the external fermion lines is,

e 2 16 π 2 p ( 3 2 + α ln α − α ln p 2 ) (3.2.11)

and the one-loop correction to the vertex function is,

e 3 64 π 4 ( 3 2 + α ln α − α ln p 2 − 2 α ) (3.2.12)

The result in (3.2.11) and (3.2.12) is of the same form as we obtained by the diagrammatic form of Operator regularization and Dimensional regularization methods in Section-2.

4. Conclusion

Radiative corrections in quantum field theory are very important for renormalization of a theory. Because when we evaluate loop-diagrams in some cases we find both finite and infinite terms. To absorb the infinite parts is the renormalization. That is why best regularization method is needed to get the correct finite and infinite parts of the loop-diagrams. Dimensional regularization (DR) is one of the best methods to evaluate the loop diagrams, but it has also some limitations in some cases. Operator regularization (OR) method had been prescribed to overcome some of these problems. To compare the advantage of OR, we have evaluated basic loop diagrams in QED using OR and compared the result with DR. We found that the result is finite and comparable, which was expected. In an earlier paper also we have shown the comparable result between OR and DR where the result consists of finite and divergent parts. From these calculations it is clear that OR is an attractive method to use in evaluating loop diagrams.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Forkan, M. and Chowdhury, M.A.M. (2018) Divergence Free QED Lagrangian in (2 + 1)-Dimensional Space-Time with Three Different Regularization Prescriptions. Journal of Applied Mathematics and Physics, 6, 2067-2086. https://doi.org/10.4236/jamp.2018.610175

ReferencesYang, C.N. and Mills, R.L. (1954) Convervation of Isotopic Spin and Isotopic Gauge Invariance. Physical Review Journals Archive, 96, 191. https://doi.org/10.1103/PhysRev.96.191Glashow, S.L. (1961) Partial Symmetries of Weak Interaction. Nuclear Physics, 22, 579-588. https://doi.org/10.1016/0029-5582(61)90469-2Weinberg, S. (1967) A Model of Leptons. Nuclear Physics, 19, 1264.Salam, A. (1968) Elementary Particle Theory. In: Svartholm, N., Ed., Weak and Electromagnetic Interactions, Almquist Wiksell, Stockholm.Glashow, S.L., Illiopulos, J. and Maiani, I. (1970) Weak Interactions with Lepton-Hadron Symmetry. Physical Review D, 2, 1285. https://doi.org/10.1103/PhysRevD.2.1285Salam, A. and Strathdee, J. (1975) Transition Electromagnetic Fields in Particle Physics. Nuclear Physics B, 90, 203. https://doi.org/10.1016/0550-3213(75)90642-2Dowker, J. and Critchley, R. (1976) Effective Langrangian and Energy-Momentum Tensor in De Sitter Space. Physical Review D 13, 3224. https://doi.org/10.1103/PhysRevD.13.3224Hooft, G.’t and Veltman, M. (1972) Regularization and Renormalization of Gauge Fields. Nuclear Physics B, 44, 189. https://doi.org/10.1016/0550-3213(72)90279-9Garcial-Aguilar, J.D. and Perez-Lorenzana, A. (2015) Dimensional Regularization with Non Beta-Functions. Revista Mexicana de Fisica, 61, 149.Pauli, W. and Villars, F. (1949) On the Invariant Regularization in Relativistic Quantum Theory. Reviews of Modern Physics, 21, 434. https://doi.org/10.1103/RevModPhys.21.434Chowdhury, A.M., et al. (1985) Preregularization. Canadian Journal of Physics, 63, 1453. https://doi.org/10.1139/p85-244Chowdhury, A.M. and McKeon, D.G.C. (1986) Surface Terms and Radiative Corrections to the VVA Triangle Diagrams. Physical Review D, 33, 598. https://doi.org/10.1103/PhysRevD.33.598Chowdhury, A.M., McKeon, D.G.C. and Mann, R.B. (1986) Ambiguities of Chairal-Anomaly Graph in Higher Dimensions. Physical Review D, 33, 3090. https://doi.org/10.1103/PhysRevD.33.3090Chowdhury, A.M., et al. (1986) Preregularization and the Structure of Loop Momentum Ambiguities within Quantum Corrections to the Supercurrent. Physical Review D, 34, 619. https://doi.org/10.1103/PhysRevD.34.619McKeon, D.G.C. and Sherry, T.N. (1987) Operator Regularization of Green’s Functions. Physical Review Letters, 59, 532. https://doi.org/10.1103/PhysRevLett.59.532McKeon, D.G.C. and Sherry, T.N. (1987) Operator Regularization and One-Loop Green’s Functions. Physical Review D, 35, 3854-3872. https://doi.org/10.1103/PhysRevD.35.3854Chowdhury, A.M. and McKeon, D.G.C. (1987) Operator Regularization in Weinberg—Salam Model. Canadian Journal of Physics, 65, 1082-1089. https://doi.org/10.1139/p87-179Shiekh, A.Y. (2010) Operator Regularization of Feynman Diagrams in One Loop QED.Shiekh, A.Y. (2011) Operator Regularization of Feynman Diagrams at Multi-Loop Order.Forkan, M. and Abul Mansur Chowdhury, M. (2012) Operator Regularization in Evaluating Feynman Diagrams in QED in (3 + 1) Dimensional Space-Time. IOSR Journals, 3, 36-51.De Witt, B. (1967) Quantum Theory of Gravity. II. The Manifest Covariant Theory. Physical Review, 162, 1195-1239. https://doi.org/10.1103/PhysRev.162.1195Kallosh, R.E. (1974) The Renormalization in Non-Abelian Gauge Fields. Nuclear Physics B, 78, 293-312. https://doi.org/10.1016/0550-3213(74)90284-3Abbott, L.F. (1981) The Background Fields Method beyond One-Loop. Nuclear Physics B, 185, 189-203. https://doi.org/10.1016/0550-3213(81)90371-0Hawking, S. (1977) Zeta Function Regularization of Path Integral in Curved Space. Communications in Mathematical Physics, 55, 133-148. https://doi.org/10.1007/BF01626516Reuter, M. (1985) Chiral Anomalies and Zeta-Function Regularization. Physical Review D, 31, 1374-1385. https://doi.org/10.1103/PhysRevD.31.1374Schwinger, J. (1951) On Gauge Invariance and Vacuum Polarization. Physical Review, 82, 664-679. https://doi.org/10.1103/PhysRev.82.664Ryder, L.H. (1987) Quantum Field Theory. Cambridge University Press, Cambridge.Ramond, P. (1981) Field Theory: A Modern Primer. Benjamin/Cummings.Lee, T.D. (1981) Particle Physics and Introduction Field Theory. Harwood Academic Publisher, Amsterdam.