Seeley-DeWitt coefficients in N = 2 Einstein-Maxwell supergravity theory and logarithmic corrections to N = 2 extremal black hole entropy

: We investigate the heat kernel method for one-loop effective action following the Seeley-DeWitt expansion technique of heat kernel with Seeley-DeWitt coefficients. We also review a general approach of computing the Seeley-DeWitt coefficients in terms of background or geometric invariants. We, then consider the Einstein-Maxwell theory embedded in minimal N = 2 supergravity in four dimensions and compute the first three Seeley-DeWitt coefficients of the kinetic operator of the bosonic and the fermionic fields in an arbitrary background field configuration. We find the applications of these results in the computation of logarithmic corrections to Bekenstein-Hawking entropy of the extremal Kerr-Newman, Kerr and Reissner-Nordström black holes in minimal N = 2 Einstein-Maxwell supergravity theory following the quantum entropy function formalism.


Introduction
Bekenstein-Hawking entropy formula is an approximated formula in two-derivative classical general relativity where one ignores the higher derivative terms in the action.Remarkably, this formula is very simple and universal in the sense that it is completely governed by horizon area, does not bother about the past history or detailed matter content of the particular black hole.In string theory, the Bekenstein-Hawking formula is subjected to two types of corrections -higher curvature corrections and quantum corrections [1].The higher curvature corrections arise due to considering higher derivative terms in the classical action, are captured by the Wald's generalization of Bekenstein-Hawking entropy.On the other hand, the quantum corrections are important for black holes of a size comparable to Planck size.The Bekenstein-Hawking formula works perfectly for macroscopically large black holes (A H ≫ 1) 1 and is expected to receive quantum corrections as one lowers the size of the black hole.In this work, we are interested in a certain class of quantum corrections as discussed below.
Recent years there have been remarkable research activities2 in order to find quantum corrections to black hole entropy formula and the quantum-corrected entropy received the following form, Here c ln A H G N is the logarithmic correction, proportional to the logarithmic of the size of the black hole, with a proportionality constant c and O(A −1 ) are power-law corrections to the leading order term in the black hole entropy formula.Although the total quantum corrections will depend on the details of the UV completion of the low energy theory, the logarithmic corrections can be determined from the knowledge of only low energy modes.These corrections can be computed from one-loop effective action of the massless fields in the particular theory.These logarithmic corrections are universal in this sense that they appear in any generic theory of gravity.But c does not have the same sort of universality: it can be determined only from low energy modes, i.e., the spectrum of massless fields present in the theory, and are unconcerned about UV completion of the theory [3][4][5][6][7][8][9][10][11][12][13].These correction terms are important as they provide us with the non-trivial information about the microstates of the black holes and hence providing a valuable testing ground for any candidate theory of quantum gravity.We still consider macroscopically large black holes to make the logarithmic terms as leading order corrections [3,4].
Quantum entropy function formalism [14][15][16] is used in [3][4][5][6][7][8] to compute the logarithmic corrections to the entropy of supersymmetric extremal black holes.This procedure uses the presence of AdS 2 factors in the near horizon geometry of extremal black holes and rules of AdS 2 /CFT 1 correspondence.However, it does not depend on supersymmetry and can be applied to any extremal black holes [9].Computation of the logarithmic corrections involves finding the one-loop effective action of the massless fields in the near horizon geometry of black hole.Evaluation of effective action was investigated by many different methods, but we consider the most potent and traditional heat kernel method (discussed in section 2.3).The main advantage of this method is that it is explicitly covariant, as well as the most suitable tool to regularize the UV divergence of the one-loop effective action.Further, in order to compute the heat kernel for one-loop effective action we pursue the Seeley-DeWitt expansion technique [17][18][19][20][21][22][23][24][25][26][27][28][29], followed by a general approach to compute the Seeley-DeWitt coefficients for fluctuations of gauge, graviton and matter fields around an arbitrary background field configuration demonstrated in [30].In this approach, the Seeley-DeWitt coefficients are computed for an arbitrary field configuration without imposing any limitations on the background geometry.Moreover, they are expressed in terms of local invariants calculated from background fields and their covariant derivatives.So one can employ these results for any arbitrary background field configuration.
In comparison, an incomplete list of many other approaches to compute heat kernel includes some perturbative expansion methods, 3 eigenfunction expansion method [3][4][5], a group theoretic approach [31][32][33], quasinormal mode approach [34], etc.But most of these approaches are either background geometry dependent or model dependent or associated with some complex computations.Of these the methods followed in [3][4][5] and [11][12][13] are closest to the one advocated in this work; so below we have given a detailed comparison between the methods.The advantages of the current approach are also outlined below eq.(3.8).The eigenfunction method to calculate the Seeley-DeWitt coefficients followed in [3][4][5] involves explicit knowledge of all the eigenfunctions and eigenvalues of the kinetic operator operating on the bosonic and the fermionic fields of the corresponding matter multiplet in the AdS 2 × S 2 background geometry.This is quite complicated and restricted to this type of geometry only.In the absence of rotational symmetry of the background geometry, the computation of eigenfunctions of the kinetic operator is indeed a difficult task.Whereas the Seeley-DeWitt expansion approach involves expansion of the heat kernel in powers of small proper time and is applicable for any arbitrary background geometry beyond AdS 2 × S 2 .Here the Seeley-DeWitt coefficients are expressed in terms of geometric invariants and can be applied for any arbitrary background.Again, the technique followed by Larsen et al. in [11][12][13] is kind of hybrid -it is an extension of the general Seeley-DeWitt approach, where computation has been done by redefining fields present in the theory which induces loss of generality within the approach.
The technical aim of this paper is twofold -first, following the general Seeley-DeWitt approach [30] (reviewed in section 3) we compute necessary Seeley-DeWitt coefficients of the quadratic order field fluctuations in the minimal N = 2 Einstein-Maxwell supergravity theory (EMSGT) and then find their useful applications in the logarithmic correction to the Bekenstein-Hawking entropy formula for extremal black holes.The results for the fermionic sector of the minimal N = 2 EMSGT are newly computed for an arbitrary background, while the bosonic sector results are reviewed and recalculated from [9].These Seeley-DeWitt coefficient results for the minimal N = 2 EMSGT in four dimensions agree with eigenfunction approach results of [5] for AdS 2 × S 2 background, as well as with the hybrid approach results of [11].Then the logarithmic corrections to the entropy of extremal black holes in minimal N = 2 EMSGT are computed using the quantum entropy function formalism.The result is stated in eq.(6.15).The result is new in the literature for extremal Kerr-Newman black holes.Any microscopic theory explaining the entropy of these black holes should reproduce the result specified in eq.(6.15).With appropriate limits eq.(6.15) produces the Kerr and Reissner-Nordström results which are consistent with [6] and [5], respectively.
The rest of this paper is structured as follows.In section 2, we will discuss the heat kernel method in four-dimensional Euclidean gravity.In this context, we will express the one-loop effective action of our theory in terms of heat kernels of the corresponding kinetic operator, and also introduce Seeley-DeWitt coefficients in the short-time expansion of these heat kernels.In section 3, we will review the general approach, necessary formulae and steps for calculating the Seeley-DeWitt coefficients.In section 4, we consider Einstein-Maxwell theory embedded in minimally coupled N = 2 supergravity in 4D and briefly explore the basic structure, equations of motion and important identities of the resulting minimally coupled N = 2, d = 4 EMSGT.In section 5, we divide the total field contents of the minimally coupled N = 2, d = 4 EMSGT into bosonic and fermionic sectors, then employ the general approach to calculate the first three Seeley-DeWitt coefficients for the two sectors separately. 4In section 6, we will use the results for Seeley-DeWitt coefficients to compute logarithmic corrections to the entropy of extremal Kerr-Newman, Reissner-Nordström and Kerr black holes in minimally coupled N = 2 EMSGT following quantum entropy function formalism, and comment on the consistency of the results.

Heat kernel in 4D Euclidean gravity
In this section, we briefly summarize the heat kernel approach in the context of 4D Euclidean gravity theory.Effective action plays a crucial theoretical role in quantum field theory and quantum gravity since it encodes almost all information about a quantum system.Unfortunately, it is associated with some mathematical pathologies and hence very difficult to compute.Heat kernel method is found to be a useful tool in this purpose.

Euclidean path integrals: the set-up
We start with the partition function Z in Euclidean path integral representation [35], is the Euclidean action of the matter fields ξ propagating through a geometry described by the metric g over an arbitrary manifold M; D[g, ξ] suggests that we have to carry out the functional integration (2.1) over all possible matter fields ξ and metric g.The leading contribution to the path integral (2.1) is expected to arise from the metric and the matter fields that are near the background fields ḡ and ξ.Here (ḡ, ξ) extremize our action (2.2), i.e. they are the solution of the classical field equations with proper boundary conditions and periodicity.At this stage, one can express g and ξ in terms of their fluctuations around the background as and expand the action in Taylor series around the classical solution (ḡ, ξ) , where S[ḡ, ξ] is the action of the classical background fields ḡ and ξ; S 2 [g, ξ] is the quadratic order fluctuated action in the fluctuations g and ξ.Note that there is no linear term in eq.(2.4) as we consider the action is not coupled to any external sources of ξ.Replacing eq.(2.4) into eq.(2.1) and ignoring the higher order terms, we have is the background field contribution to the partition function while the second term of eq.(2.5) is regarded as the contribution of the matter quanta and thermal gravitons (fluctuation of metric) on the background geometry [36,37].

The one-loop effective action
The quantum corrected one-loop effective action W [38,39] is expressed as The quadratic term of action S 2 in eq.(2.4) can be expressed as where Λ is a differential operator which controls the quadratic order quantum fluctuations around the classical background.For bosonic fields, Λ is a second-order differential operator.Again, it is real, elliptic and self-adjoint since our background metric ḡ has Euclidean signature (i.e.positive and real).But for fermionic fields, Λ is of first-order.We will derive all the important results for bosonic fields first and then extend them for fermionic fields.
If ξ is a graviton or gauge field then Λ will receive a complicated tensor structure, which has been suppressed here for simplicity.
Considering Λ has a discrete spectrum of eigenvalues λ i , the one-loop effective action given eq.(2.6) will take the following form [40] where χ = ±1 for bosons and fermions respectively.Note that the spectrum of Λ may in practice also be continuous, and our basic structure of the heat kernel method and the general Seeley-DeWitt computation approach (sections 2.3 and 3) are valid for both the discrete and continuous spectrums.

The heat kernel method
Practically, one finds it difficult to compute W directly, especially when fluctuations are gauge fields and gravitons.The most straightforward method is to compute the summation in eq.(2.8) provided we have the knowledge about the complete basis of normalized eigenfunctions {f i } and eigenvalues {λ i } obeying Λf i = λ i f i .But this is a challenging task.Also, from the series in eq.(2.8) it is evident that the one loop determinant detΛ diverges (UV-divergence in quantum field theory) and traditional way to regulate this is the heat kernel method : expressing W in terms of diagonal matrix element of heat kernel via a special integral representation, called proper time representation (2.15).
In terms of eigenfunctions {f m i } and eigenvalues {λ i } of the kinetic operator Λ, heat kernel (of the operator e −sΛ ) can be defined as [5,6] where m, n are set of indices of the eigenfunctions at two points x, y respectively on M; s is an auxiliary time variable, called proper time.The eigenfunctions {f m i } have been normalized so that where G mn is a metric in the space of eigenfunctions (fields) induced by ḡ.Contracting over indices m, n, let us define the heat kernel K(x, y; s) K(x, y; s) = G mn K mn (x, y; s), (2.11) which satisfies the diffusion equation with boundary conditions K(x, y; 0) = δ(x, y), where Λ x describes the fact that in eq.(2.12) the operator Λ is supposed to act only on the first argument of K(x, y; s).The trace of heat kernel D(s), called heat trace, can be obtained by setting x = y (for avoiding the dependence of eigenfunctions) in eq.(2.11) and then integrating over all the M: Heat trace D(s) is important because it stores information about the spectrum of the kinetic operator Λ.For example, one may evaluate 5 This is called Schwinger-DeWitt proper time representation [41,42] of the effective action, where ϵ is an ultraviolet cutoff, naturally restricted by the Planck length, i.e. ϵ ∼ l p 2 ∼ G N .Also, one may identify the one-loop corrected effective Lagrangian density ∆L eff from eq. (2.15) as The main benefit of this approach lies in the fact that heat kernels are ultraviolet finite for all positive proper time (s). 5 In deriving eq.(2.15), we have used the identity limϵ→0

General approach for computing Seeley-DeWitt coefficients
The most common method for calculating Seeley-DeWitt coefficients in physics is the De-Witt iterative technique (reviewed in sec.4.3 of [30]).But, we will follow the method of Gilkey [43] (explained in sec.4.1 of [30]) which seems to be more efficient for a curved background with arbitrary spin.In this section, we will review the main ideas, the basic structure and useful formulae to compute the Seeley-DeWitt coefficients based on the manual [30].We now start by considering the quadratic order fluctuated action (2.7) controlled by the kinetic differential operator Λ and rewrite its schematic form in terms of field indices as where { ξm } is the set of all fluctuating fields.In order to successfully apply the procedure, Λ is constrained by the following conditions: ■ Λ must be of a second order differential operator.It is useful to note that the quadratic order fluctuations of the fermions are governed by first order operators / D, i.e. δ 2 L = ξ / D ξ.There we first restructure them as Dirac-type and Hermitian one, and define [5] det where This will include an additional 1/2 factor in the description of heat kernel (or in the formula of Seeley-DeWitt coefficients), which is compensated by a factor of 2 arises due to the fact that squaring the complex operator doubles the integration variables.
■ Λ must be made up of pseudo-differential operators which excludes negative and fractional powers of the derivatives.
■ Λ should be in a minimal form of type g µν D µ D ν + Q, where it has a scalar principal part, i.e. second order derivatives are contracted with the metric.
In summary, we conclude that the present method demands Λ to be a Laplace-type operator having the following schematic form where G mn is some arbitrary matrix, 7 D ρ is the ordinary covariant derivative with connections determined by the background metric, and N ρ , P are arbitrary matrices constructed from the background fluctuations.Let us define a new covariant derivative D ρ including the field connection ω ρ , and rewrite eq. ( 3.3) by replacing D ρ with D ρ as We aim to express eq.(3.5) in the standard form and extract the following necessary matrices, where is the curvature associated with D ρ .So, one needs to identify the matrices G, P and N from the relation (3.5) and use those results in the formula (3.7) for obtaining ω, E and Ω.With these definitions, the formulae for the first three Seeley-DeWitt coefficients [30,44] are summarized as following 7 Precisely, G is the effective metric tensor which contracts indices of an operator acting on ξm and we can also use G to find trace values of these operators.For example, if G mn be the metric tensor related to an arbitrary two rank tensor B mn , then The formulae given in eq.(3.8) are applicable only for manifolds having no boundary. 8 Hence, all the total derivative terms have been ignored as they will appear as boundary terms in the integral (2.17) and hence will disappear after integration.Again for manifolds satisfying the condition R = 0 one can neglect terms proportional to R in this formulae.We consider free 9 quantum fields on a Riemannian manifold having no boundaries so that the Seeley-DeWitt coefficients can be described by local invariants induced from the background fields and metric.The main advantage of calculating the Seeley-DeWitt coefficients following this general Seeley-DeWitt approach is -(a) after considering quadratic fluctuations for different fields in a given theory we have simple and straightforward formulae to compute the Seeley-DeWitt coefficients, (b) the coefficients are expressed in terms of geometric invariants, hence one can employ these results for any arbitrary background field configuration without imposing any limitation on the background geometry.It is important to mention that for the fermionic case, we need to multiply an additional minus sign (correspond to χ = −1) in the above formulae of Seeley-DeWitt coefficients.We also note that the present method assumes the fermionic field fluctuations as Dirac spinors, and hence for casting Majorana and Weyl (having half the degrees of freedom of Dirac spinor) fermion contribution, one needs to multiply a factor of −1/2 in the formulae (3.8).Finally, we summarize the key steps of computing Seeley-DeWitt coefficients as follows: 10   1.Consider fluctuations for all possible fields present in the theory around a classical background and express the action in quadratic order of fluctuations.
2. Gauge fix the theory of choice by adding appropriate gauge fixing term to the quadratic action and find the quadratic fluctuation controlling operator Λ (defined in the eq.(3.1)).
Identify the appropriate ghosts and keep them separate.
3. With all the restrictions mentioned, restructure Λ so that we can have the canonical form (3.3).Extract the results for the matrices G, P and N .also.But in this paper, our main motive is to compute logarithmic corrections of black hole entropy which demands to choose the no boundary formulae (discussed in section 6). 9"Free" means we are considering only interaction with background geometry, all other interactions are turned off. 10 For a better understanding of this method and some elementary examples, readers may refer [11,44].

Physical Field Content
Bosonic sector: one graviton (g µν ), one graviphoton (A µ ) Fermionic sector: two gravitino (ψ µ , ϕ µ ) 4 Minimal N = 2, d = 4 Einstein-Maxwell supergravity theory In this section, our motive is to briefly summarize a minimally coupled N = 2 supergravity embedded with the solutions of Einstein-Maxwell theory in four dimensions.There are a class of N = 2 supergravity solutions to the equations of motion, determined by certain conditions [45].These conditions reduce the full N = 2 supergravity equations of motion to a much simpler form of equations of motion for Einstein-Maxwell theory.The resultant minimal N = 2, d = 4 Einstein-Maxwell supergravity theory (EMSGT) plays a significant role in string phenomenology, black holes, and holography.In general, N = 2 supergravity theories are coupled to additional matter fields like n V vector multiplets and n H hyper multiplets [45,46].But in this work, we will focus on the simpler case of minimally coupled N = 2 supergravity where n V = n H = 0.However, both the theory have the same type of classical or black hole solutions.At any point, one can compute the Seeley-DeWitt coefficients and the logarithmic contributions for n V and n H fluctuations separately, and can extend the minimal supergravity results to these matter coupled theories.The minimal N = 2 EMSGT is associated with the field content provided in table 1.Both the bosonic and fermionic sectors of the field content are described by the Lagrangians L EM and L gravitini respectively in the action (4.1).Note that all these fields are free in the sense that they are only coupled to the metric (or graviton).We consider standard d = 4 Einstein-Maxwell theory embedded in minimal N = 2 supergravity described by the following action, with where R is the Ricci scalar induced from metric g µν , F µν the electromagnetic field strength of gauge field A µ in Einstein-Maxwell theory, ψ µ , ϕ µ are two gravitini species and H µν = −ε µνρσ F ρσ is the dual field strength.11ψ µ , ϕ µ are Majorana spinors in 4D and the action (4.1) describes black holes in pure N = 2 supergravity.It is useful to note that we have set 12 G N = 1/16π.

Equations of motion and identities
Let us denote (ḡ µν , Āµ ) as an arbitrary solution to the classical equations of motion of our theory and define corresponding background field strength tensor Fµν ≡ ∂ µ Āν − ∂ ν Āµ .
Then the geometry and matter fields are related by Einstein equation, The background gauge field strength Fµν satisfies Maxwell equations and Bianchi identity, Eq.( 4.3) and eq.( 4.4) represents the basic equations of motion of N = 2, d = 4 EMSGT.Again, the trace of Eq.( 4.3) yields R = 0, ( Hence for the remaining parts of this paper, we will use this condition and ignore the terms proportional to R. Again using (4.5), (4.3) reduces to, The field content is divided into two sectors -fermions and bosons, for which the contributions will be calculated separately.We actually compute results for the fermionic sector, while the bosonic sector results are reviewed and recalculated from [9].Finally, we will sum up these results in order to find the complete Seeley-DeWitt coefficients of the theory.

N = 2 SUGRA: Bosons
Let's now consider the Einstein-Maxwell theory described by the Lagrangian L EM in (4.1) with the action, We aim to calculate the Seeley-DeWitt coefficients for this theory by considering quadratic order fluctuations around a classical background.This methodology has been developed in [9], where the authors have computed only the a 4 coefficient for calculating the logarithmic correction to extremal Kerr-Newman black hole entropy.But we find it worthful to briefly review this work and calculate all the first three Seeley-DeWitt coefficients for the theory, which may be useful in other applications.Let (ḡ µν , Āα ) be any arbitrary classical solution to the equations of motion of Einstein-Maxwell theory, and around this classical background we consider following fluctuations of metric and gauge field 13 where h µν and a α are the graviton and graviphoton respectively.Moreover, the theory is gauge fixed by adding the following gauge-fixing term, and the corresponding ghost action is added to the action (5.1).Up to quadratic order in the fluctuations the gauge fixed action excluding the ghost part is given by [9], where we denote □ ≡ D µ D µ , h = h µ µ = tr(h µν ), Fµν ≡ ∂ µ Āν − ∂ ν Āµ and ξm = {h µν , a α }.Kindly note that Λ is not in Laplace-type form due to the presence of the kinetic term h□h.We need to restructure Λ and express it in the canonical form.In order to proceed further, we introduce the operator G µν αβ 14 G µν αβ = 1 2 (ḡ µα ḡνβ + ḡµβ ḡνα − ḡµν ḡαβ ), (5.5) 13 The graviton and graviphoton have been scaled in such a manner that the kinetic terms for them have the same normalization in the quadratic order fluctuated action (5.4). 14The general form of G is ξmG mn ξn = hµν G µν αβ h αβ + aαG αβ a β , where the form of G µν αβ is given in eq. ( 5.5) and G αβ = ḡαβ .
and rewrite the kinetic part as h µν □h µν − 1 2 h□h + a α □a α = G µν αβ h µν □h αβ + ḡαβ a α □a β .With necessary modifications of the total derivative terms to make Λ Hermitian, the resulting form of Λ becomes where We identify the two matrices P and N (defined in (3.3)) as (5.8) By comparing Λ from eq. (5.6) with the standard form (3.6), one can extract the useful matrices I, ω, E and Ω, and find their traces in order to compute the Seeley-DeWitt coefficients.It is now straightforward to express the identity matrix as ξm I mn ξn = h µν G µν αβ h αβ + a α ḡαβ a β . (5.9) Making use of the results of N and P from eq. (5.8) in the formulae (3.7), we express Note that in writing the expression for Ω, we have used the following commutation relation, (5.11) The useful traces for the matrices are finally written down in (5.12).We skip the details of the intermediate steps as the calculation is straightforward although tedious.
tr(I) = 10 + 4 = 14, (5.12) Substituting these trace values in the formulae (3.8), we find the Seeley-DeWitt coefficients of Einstein-Maxwell theory (without ghost).The results are (5.13) The ghost fields are non-interacting and hence can be treated separately.The corresponding ghost action is given by, where b µ , c µ (vector fields) are diffeomorphism ghosts related to the graviton; b, c (scalar fields) are the ghosts associated with the graviphoton.Adjusting this action up to total derivative, we express the corresponding Λ as, and comparing with (3.3) identify the matrices N and P as We follow the general approach described in section 3 and making use of formulae (3.7) we extract the following necessary matrices (5.17) Using (5.17), one can compute the following traces for the ghost fields, tr(I) = 4 + 4 + 1 + 1 = 10, and can use them in order to find Seeley-DeWitt coefficients for ghost fields corresponding to the bosonic action of the theory.This yields Note that an overall minus sign is put on each coefficient since ghosts fields follow the reverse of usual spin-statistics.

N = 2 SUGRA: Fermions
The fermionic sector of N = 2 includes two Majorana gravitini fields described by the Lagrangian L gravitini expressed in (4.2).To gauge fix the fermionic action we add the following term to L gravitini , and we can express the quadratic part of the corresponding gauge fixed action in Dirac form as Now, we can rewrite the action for gravitini fluctuations ξm = {ψ µ , ϕ ν } as (5.22) The action contains first-order Dirac-type differential operator / D, but for computing the Seeley-DeWitt coefficients following general approach section 3 we need to make it quadratic in order.With our choice of four-dimensional Euclidean space-time, gamma matrices are Hermitian, the gravitino conjugates are ψµ = ψ † µ , φµ = ϕ † µ , and (γ ν γ ρ γ µ D ρ ) † = −γ ν γ ρ γ µ D ρ .We will consider F µν to be real.Hence / D turns out to be Hermitian, i.e. / D † = / D. Now one can follow the method mentioned in section 3 to obtain the necessary Laplace-type operator Λ of the fermionic sector as (5.23) The second derivative terms of eq. ( 5.23) can be simplified to the following form [44], (5.24) Using (5.24) and re-arranging the terms (5.23) can be expressed as, (5.25) So we have adjusted Λ as the prescribed form (3.5), where the matrices N and P hold the following form As mentioned in section 3, now our task is to express Λ as the standard form (3.6).Then, using the formulae (3.7) we find and also identify the corresponding I, E and Ω as the following compact form (5.28) We then computed all the necessary trace values using the information (5.28) to obtain the Seeley-DeWitt coefficients for the gravitini fields, (5.29) Although these computations are extremely lengthy, one can proceed in a systematic way (refer Appendix A for details of this computation).Making use of these trace results (5.29), one computes the first three Seeley-DeWitt coefficients (3.8) for the fermions (two gravitini) in our N = 2 SUGRA.The results are (5.30) Note that we have put an overall −1/2 factor by hand on each coefficient to justify the fermion and Majorana condition, as discussed in section 3.
Next, the choice of gauge fixing in eq.(5.21) introduces six bosonic ghosts in our system with the action, where bA , cA , and ẽA are spin half ghosts with A = 1, 2 for the two gravitini species.These Majorana ghosts are minimally coupled, hence one can obtain the Seeley-DeWitt coefficients as 6 times the spin-1 2 Seeley-DeWitt coefficients with an opposite sign, i.e. a 2n ghost,f = −6a 1/2 2n .The detailed calculations for a 1/2 2n have been described in [44] 15 for a general background, and it is very much straightforward.This led us to express the results for ghost part of fermionic fields as (5.32) 15 See § 2.2.1 of [44] and finally, use its eq.2.19 with R = 0 for a 1/2 2n .

Combined Seeley-DeWitt coefficients in minimal N = 2 EMSGT
In previous subsections, we have calculated Seeley-DeWitt coefficients for the field fluctuations of the bosonic and fermionic part of N = 2 SUGRA and also computed the same for corresponding ghost parts separately.These ghost fields have no interaction with any of graviton, graviphoton and gravitini fields.Therefore, the contributions from the ghost fields can be directly added to the results from the gauge fixed part.The total Seeley-DeWitt coefficients are obtained by adding the bosonic and fermionic contributions.
Combining the results (5.13) and (5.19) we get the Seeley-DeWitt coefficients for the bosonic sector of N = 2 SUGRA as This clearly reflects the result derived in [9].Similarly, for the fermionic part, we sum up the results (5.30) and (5.32).This yields Our results for the fermionic sector also agrees with eigenfunction approach results [5] for AdS 2 × S 2 background, as well as with the hybrid approach results [11]. 16We then find the full N = 2 SUGRA Seeley-DeWitt coefficients with the contribution graviton, graviphoton, and gravitini fluctuations as a N =2 2n = a b 2n + a f 2n .The results are (5.35) The following section deals with the use of these coefficients in the logarithmic correction to the entropy of extremal black holes.
6 Application of Seeley-DeWitt coefficients: logarithmic corrections to the black hole entropy In this section, we will summarize the general procedure to find the logarithmic correction to extremal back holes in an arbitrary theory by making use of the Seeley-DeWitt coefficient a 4 (x) of the same theory [14][15][16]35].Then we will use this procedure for extremal black holes in minimal N =2 Einstein-Maxwell supergravity theory.It has been observed that Euclidean quantum gravity [35] gives a more detailed picture of the quantum corrections [3,4,6] to the black hole entropy.We will follow the quantum entropy function formalism [14][15][16], in order to calculate the logarithmic corrections to the entropy of string theoretic black holes.For extremal black holes, AdS 2 /CFT 1 correspondence yields an auxiliary way of describing the horizon degeneracy: they are equivalently given by finite (boundary independent) part of the partition function Z AdS 2 ×K of the string theory in the near horizon geometry structured as AdS 2 ×K, where K is compact space fibered over AdS 2 .The classical limit provides us with the relation: Z finite AdS 2 ×K = e S BH , where S BH is the quantum corrected black hole entropy.The partition function Z AdS 2 ×K computed by carrying out Euclidean path integral in AdS 2 × K weighted by e −S (S is Euclidean action), over all string fields that asymptotically approaches to the near horizon geometry of black hole and for Z finite AdS 2 ×K one needs to consider only one-loop term which is finite and boundary independent [3][4][5][6].Logarithmic corrections (∆S BH ) are special type of one-loop quantum corrections, requires integrating out non-zero modes of massless fields17 in Z AdS 2 ×K , and then find one-loop effective action W or effective Lagrangian density ∆L eff in quadratic order fluctuation of S (see section 2) in the background AdS 2 × K.We find logarithmic corrections as Here L (nz) eff is the one-loop corrected effective Lagrangian density with the contribution of only non-zero modes of massless fields, i.e. one needs to throw away zero-modes contribution from the overall result.Then we express [5,6] where K(x, x; 0) is the zero mode contribution in K(x, x; s). 18Note that logarithmic part arises only from the s independent part of the expansion (2.17) of K(x, x; s) after integrating over the range ϵ ≪ s ≪ A H , where ϵ ∼ G N = 1/16π.Finally, we express the key formula for finding logarithmic correction to extremal black hole entropy as where the non-zero mode and zero mode contributions are respectively, where β is the scaling dimension of fields arises for incorporating corrections due to zero mode contribution in ∆S BH [5,6].One can evaluate C local using by computing a 4 of massless fields in near horizon geometry of concerned black hole, while C zm can be obtained from the following special combined formula [11], Here K (number of rotational isometries) is 3 if for black holes having J = 0, otherwise it is 1; N SUSY (number of preserved supercharges) is 4 for BPS black holes, 0 otherwise.In order to compute the logarithmic correction to extremal black hole entropy in any arbitrary theory, one needs to first compute a 4 (x) for the quadratic fluctuations of massless fields in the theory.Then find the results for C local (using the a 4 (x) result) and C zm respectively, and finally use these results in the formula (6.4).With cognizance of the above structure, our next aim to compute logarithmic corrections to the entropy of extremal Kerr-Newman, Reissner-Nordström and Kerr black holes19 in minimal N = 2 EMSGT by making use of a 4 (x) result (5.35) of the theory.

Logarithmic corrections to N = 2 extremal black hole entropy
We will start with a Kerr-Newman black hole [9] to calculate associated logarithmic correction.In appropriate units, the metric of a general Kerr-Newman black hole characterized by mass M , charge Q, and angular momentum J is expressed as where a = J/M .The horizon is situated at and the corresponding classical Bekenstein-Hawking entropy of black hole having a horizon area (6.10) The extremal limit M → a 2 + Q 2 yields .11)Now for the metric (6.8), we have [47,48]  4 , det ḡ = (r 2 + a 2 cos 2 θ) 2 sin 2 θ. (6.12) In the extremal limit, one can achieve [9] horizon where b = a/Q = J/M Q.
■ Using the a 4 (x) result (5.35) for minimal N = 2 EMSGT and the limits (6.13) in the relation (6.5), and also making use of the formula (6.7), for extremal Kerr-Newman black hole (J ̸ = 0, Q ̸ = 0) we have This is one of our central results.Similarly, we can get the results for Kerr and Reissner-Nordström black holes by applying appropriate limits as discussed below.
■ For an extremal Kerr black hole (J ̸ = 0, Q = 0), we have to take b → ∞ in the limits of background invariants given in (6.13).This yields hence the logarithmic correction becomes This is in accordance with the result in [6] (set n V = 1, n 3/2 = 2, n S = n F = 0 in eq.(3.12) and include an additional −4 ln A H term for the correction due to zero-modes of gravitini).
■ For the case of an extremal Reissner-Nordström black hole (J = 0, Q ̸ = 0), one has to put the limit b → 0 in the limits of background invariants given in (6.13).This gives hence the logarithmic correction becomes This agrees with the result in [5] (see eq.(5.32)).
To conclude, we have computed the first three Seeley-DeWitt coefficients for minimal N = 2 Einstein-Maxwell supergravity theory following the general Seeley-DeWitt approach [30] (reviewed in section 3).The fermionic sector results are newly calculated in this approach, while the bosonic sector results are reviewed and recalculated from [9].This computation does not impose any limitations on the background geometry.The coefficients are expressed in terms of the local invariants of the theory and can be evaluated for any arbitrary background.These Seeley-DeWitt coefficient results for the minimal N = 2 Einstein-Maxwell supergravity theory in four dimensions agree with eigenfunction approach results [5] for AdS 2 ×S 2 background, as well as with the hybrid approach results [11].As an application of these results, we have computed the logarithmic corrections to Bekenstein-Hawking entropy of the extremal Kerr-Newman, Kerr and Reissner-Nordström black holes in minimal N = 2 Einstein-Maxwell supergravity theory by properly integrating these coefficients in the near horizon geometry of these black holes.The Kerr and Reissner-Nordström results agree with [6] and [5] respectively, while the Kerr-Newman result is new to the literature.Any microscopic theory explaining the entropy of such black holes is expected to reproduce the subleading logarithmic correction terms specified in eq.(6.15) in addition to the leading Bekenstein-Hawking term.
Making use of expressions of ω and P from eq. (5.27) and eq.(5.26) respectively, one can simplify E to the following form  (A.9) Also, Therefore, . Then, Y U i V j = 0, (A. Then, L 2 P 2 = 0, So, Hence, we have

4 .. 7 .. 8
Define a new effective covariant derivative D ρ with field connection ω as eq.(3.4) and making use of it, modify Λ further to have the form (3.5). 5. Comparing the new form of Λ with standard form (3.6) read off the results for matrices I, E and Ω using eq.(3.7).6. Compute tr (I), tr (E), tr (E 2 ) and tr (Ω 2 ).Use this data in the formulae (3.8) to find the result for the first three Seeley-DeWitt coefficients (without ghosts) a no ghost 2n Making use of background equations of motion, gamma and tensor identities, simplify a no ghost 2n There are Seeley-DeWitt coefficient formulae for manifold with boundary (see eqs. (5.29)-(5.33) of [30])

)
In addition, the dual gauge field strength Hµν satisfies the Maxwell-Bianchi equations of the form, D µ Hµν = 0, D [µ Hνρ] = 0. (4.7)It is also important to mention the gravitational Bianchi identity R µ[νρσ] = 0. (4.8)We have used eqs.(4.4) and (4.6) to (4.8) to simplify our calculations and express them in closed forms.5 Seeley-DeWitt coefficients in minimal N = 2, d = 4 Einstein-Maxwell supergravity theory In this section, we will calculate the first three Seeley-DeWitt coefficients of the quadratic order fluctuations of minimally coupled N = 2, d = 4 EMSGT fields discussed in section 4.