Fluctuations of Quantum Fields via Zeta Function Regularization
Guido Cognola1, ? , Emilio Elizalde2, ? and Sergio Zerbini1, ?
Dipartimento di Fisica, Universit` a di Trento and Ist
ituto Nazionale di Fisica Nucleare Gruppo Collegato di Trento, Italia Instituto de Ciencias del Espacio, IEEC/CSIC, IFAE, and University of Barcelona, Spain
arXiv:hep-th/0201152v2 5 Aug 2002
Explicit expressions for the expectation values and the variances of some observables, which are bilinear quantities in the quantum ?elds on a D dimensional manifold, are derived making use of zeta function regularization. It is found that the variance, related to the second functional variation of the e?ective action, requires a further regularization and that the relative reg2 ularized variance turns out to be N , where N is the number of the ?elds, thus being independent on the dimension D . Some illustrating examples are worked through. The issue of the stress tensor is also brie?y addressed. 05.30.-d,05.30.Jp,11.10.Wx,11.15.Ex
Typeset using REVTEX
firstname.lastname@example.org, email@example.com firstname.lastname@example.org
? email:email@example.com, ? email:
Vacuum ?uctuations play an important role in many physical processes. The Casimir e?ect is one of its most interesting physical manifestations and it has been experimentally veri?ed. It is also well known that the Casimir e?ect is related to the presence of a non vanishing vacuum energy (see for example [1,2]). This fact mainly occurs when one is dealing with non trivial space-times, where the spatial sections are topologically non trivial spaces or manifolds with boundaries. Another interesting issue where quantum vacuum e?ects are present is the Quantum Field Theory in curved space-time [1,3,4]. Recall that within the semi-classic approach to quantum gravity, the basic equation reads G?ν + g?ν Λ = 8πG < T?ν > , (1)
where G?ν is the Einstein tensor, Λ the cosmological constant and T?ν is the vacuum expectation value of the matter stress tensor (we use h ? = c = 1). As a consequence, ?uctuations of the stress tensor can induce ?uctuations of the classical gravitational ?eld and, in order to justify the semiclassical approximation, it appears very important to have reliable a priori estimates of these ?uctuations. Fluctuations of the stress tensor were studied in [5,6] making use of canonical methods. Fluctuations of Casimir forces were investigated in [7,8]. Alternatively, other authors  have investigated the same problem by making use of zeta function regularization. With regard to this issue, it is well known that the notion of e?ective action (or e?ective potential) plays an important role as a powerful tool in relativistic quantum ?eld theory. This quantity, however, is ill de?ned, since, within the Euclidean formulation and in the one-loop approximation, the one-loop e?ective action contains functional determinants of elliptic operators, which have to be regularized. Zeta-function regularization [10–12], (see also  for the generalization to elliptic pseudo-di?erential operators and [15–17] for physical applications) was introduced by a number of authors as a convenient tool in order to deal with the evaluation of functional determinants. It permits to give a meaning in the sense of analytic continuation —a mathematically very precise procedure— to quantities that are formally divergent. In this paper, we would like to revisit the zeta function regularization approach for the evaluation of expectation values < O > and their quantum ?uctuations. It is our opinion that zeta function regularization is a very powerful tool at our disposal, as compared with other methods and that this issue deserves a careful investigation. ? We will mainly consider two quantities: O = φ2 and the stress tensor trace O = T? 2 2 2 and their corresponding variances ? O =< O > ? < O > . Within our formalism, it is convenient to introduce the relative variance  ?2 rO = ?2 O . < O >2 (2)
A di?erent relative variance, though directly related to the previous one, has been introduced by Kuo and Ford . It reads
? 2O = ? r
?2 ?2 O rO = . 2 <O > 1 + ?2 rO
? Considering now the operator O = T? , we observe that the trace of the Einstein equations reads ? G? ? + D Λ = 8πG < T? > .
Thus, ?uctuations of the stress tensor trace induce classical ?uctuations of the scalar curvature. Furthermore, this trace ?uctuation may be used to have an estimate of the ?uctuations of the whole stress tensor, since for conformally invariant quantum ?elds in homogeneous space-times, one has g?ν α < Tα > . (5) < T?ν >= D The issue concerning the validity of semiclassical gravity has been discussed in [6,18]. We also recall that the ?rst variation of the e?ective action is related to the vacuum expectation value of physical quantities, while the second variation of the e?ective action is associated with the variance. Within the zeta function regularization procedure, one has to deal with traces of complex powers of elliptic operators. The ?rst variation of the e?ective action is well de?ned by the use of zeta function regularization, while the second variation is intrinsically ill de?ned, unless one makes use of suitable variations with disjoint supports. In the coincidence limit, the physically interesting case, one has to make use of a further regularization . Our main result is the following: modulo regularization problems, the relative variance 2 2 ?2 turns out to be ?2 r O = N , (thus ?r O = N +2 ), where N is the number of scalar ?elds in some multiplet. This result seems to be general, that is independent on the quantity one is dealing with, for example the stress tensor trace. Our results are compatible and should be compared with the ones recently obtained, for N = 1, regarding the vacuum energy density ?uctuations via smeared ?elds and point separation , which give relative variances of the order of unity albeit dimensionally dependent. On the other hand, coming back to the case of N neutral scalar non-interacting ?elds, we recover a well known criterion for neglecting the quantum gravity ?uctuations in the large N limit . The content of the paper is as follows. In Sect. II, zeta function regularization and heatkernel techniques are brie?y summarized. In Sect. III, the ?rst variation of the e?ective action computed is shown to be related to the vacuum expectation values of observables. In Sect. IV, the second variation is shown to be associated with the variance and the ?nal result is presented with the help of a further analytic and ad hoc regularization. In Sect. V, some examples are presented. The paper ends with some concluding remarks and an Appendix, where the ?rst and second variations of the trace of the complex power of an elliptic operator are explicitly computed.
II. ZETA FUNCTION REGULARIZATION OF THE EFFECTIVE ACTION
In this section, we will summarize some basic aspects of the heat-kernel and zeta function regularization methods. For the sake of simplicity we shall here restrict ourselves to scalar 3
?elds, but the method is also valid in more general situations. In the case of a neutral scalar ?eld, the one-loop Euclidean partition function, reads  Γ(1) = ? ln Z = 1 A ln det 2 , 2 ? (6)
where ?2 is a renormalization parameter. To begin with, recall the de?nition of the zeta function related to an elliptic operator A Γ(s)ζ (s|A) = 1 Γ(s)
ts?1 Tr e?tA dt ,
which is valid for large values of ?s. For small t, and for second order elliptic di?erential operators, it will be assumed that the heat-trace has the following asymptotics
r ?D 2
where Kr are the integrated heat-kernel coe?cients. In principle they may be computed (see, for example ). We will also assume the validity of the local heat-kernel expansion, namely
r ?D 2
If the asymptotic expansions (8) and (9) hold, a standard argument leads to the following meromorphic extension of the zeta function and its local counterpart (J (s) being an analityc function) Γ(s)ζ (s|A) = Kr r ? D + J (s ) . r =0 s + 2 Kr (x) r ?D + J (s, x) . r =0 s + 2
Γ(s)ζ (s|A)(x) =
As a result, at s = 0, the global and local zeta functions are regular and their derivatives exist. We also have ζ (0|A)(x) = KD (x). As is well known, KD (x) = 0 in any odd dimensional manifold without boundary. Under hypothesis above, zeta function regularization corresponds to doing the following: ln det A d A ′ = ? ζ (0 | ) = ? lim Tr(?2s A?s ) . s→0 ds ?2 ?2 (12)
However we would like just to recall that, in some situations, it might be necessary to generalize the above de?nition in the form ln det 1 d2 A = ? lim Tr(s?2s A?s ) . s→0 2 ds2 ?2 4 (13)
When Tr(A?s ) is regular at s = 0, the two de?nitions (12) and (13) coincide, but, in some cases, ζ (s|A) has a simple pole at s = 0 and so Eq. (12) is no more valid. This fact may well be present at the level of the e?ective action (see for example [22–24]). For the sake of simplicity, in the sequel we will assume the validity of Eqs. (10) and (11). The latter gives ζ (0|A)(x) = KD (x) . (14)
For future use, we also observe that when KD?2 (x) = 0, the local zeta function has a simple pole at s = 1. In fact we may write ζ (1 + s|A)(x) = 1 KD?2 (x) + G(1 + s, x) , Γ(1 + s) s (15)
where G(1, x) is a regular function. It is given by G(1, x) = PPζ (1|A)(x) ? γ KD?2 (x) , γ being the Euler’s constant and PPζ (1|A)(x) = lim ζ (1 + s|A)(x) ?
KD?2 (x) s
the ?nite part of ζ (s|A) at s = 1. When we consider an odd dimensional manifold without a boundary, this singularity is absent.
III. THE FIRST VARIATION OF THE EFFECTIVE ACTION: VACUUM EXPECTATION VALUES
In this section, we will evaluate vacuum expectation values < O > of some speci?c quantities O , such as the stress tensor trace or conformal anomaly and the square of the ?eld (?uctuation). These quantities involve the product of two quantum ?elds at the same point, x, and are therefore ill de?ned. They require a regularization. We shall consider a multiplet of N scalar ?elds denoted by φ in an external ?eld, described through a classical action given by S= 1 2 dx φLφ , (18)
where L is a suitable (matrix valued) di?erential operator de?ned on a D -dimensional smooth manifold. To begin with, let us recall the formal trick that allows one to get the vacuum expectation value of the bilinear O = φKφ within a path integral approach . We will consider two cases: K = I ,the identity matrix, in the case of the ?eld ?uctuation O = φ2 , and ? K = c1 + c2 L, in the case of stress tensor trace O = T? . Here, c1 and c2 are constants and moreover, in the conformally invariant case, c1 = 0. If we denote by α(x) a suitable classical source, we may consider the Euclidean partition function 5
Z (α ) =
L(α) = det 2 ?
Here we are assuming that, in the massive case, the multiplet has the same common mass. In this way, there is no multiplicative anomaly (see, for example, ) and L(α) = L + αK may be regarded as a simple di?erential operator. A formal functional derivation leads to < O (x) >= ?2 δ ln Z (α) δα .
Making use of (19) and zeta function regularization (see Eq. (59) in the Appendix), we may give a meaning to the above formal expression, namely δ ln Z (α) = ? N d 2s ? sTr L?s?1 (α)δL(α) lim s → 0 2 ds . (21)
In the case O = φ2 , δL(α) = δα. Thus < φ2 (x) >= N lim d s?2s ζ (s + 1|L)(x) . s→0 ds (22)
In the equations above, ζ (s|L)(x) is the local zeta function. As a result, making use of the meromorphic expansion (15), one gets (see [27–29]) < φ2 (x) >= N PPζ (1|L)(x) + NKD?2 (x) ln ?2 . When D is odd and the manifold is without boundary, we simply have < φ2 (x) >= Nζ (1|L)(x) . (24) (23)
In the other case, namely when one is dealing with the stress tensor trace, δL(α) = δα(c1 + c2 L). As a consequence,
? < T? (x) >= N lim s→0
d s?2s (c1 ζ (s + 1|L)(x) + c2 ζ (s|L)(x) , ds
and, as a result,
? < T? (x) >= c1 < φ2 (x) > +c2 N KD (x) .
In the conformally coupled case c1 = 0 and one has the usual conformal anomaly, due only to quantum e?ects.
IV. THE SECOND VARIATION OF THE EFFECTIVE ACTION: THE VARIANCE
We have seen that the ?rst variation of the e?ective action is associated with the vacuum expectation value < O > of bilinear quantities in quantum ?elds. Let us show now that the second variation of the e?ective action is related to the variance ?2 O =< O 2 > ? < O >2 . To begin with, the second variation of the partition function (19) gives < O (x)O (y ) >= 4 On the other hand, we have the identity δ 2 ln Z (α) 1 1 δZ (α) δZ (α) δ 2 Z (α ) = ? 2 . δα(x)δα(y ) Z (α) δα(x)δα(y ) Z (α) δα(x) δα(y ) As a consequence, the variance is given by ?2 O (x, y ) =< O (x)O (y ) > ? < O (x) >< O (y ) >= 4 Now, it appears convenient to introduce the relative variance ?2 r O (x, y ) = ?2 O (x, y ) . < O (x) >< O (y ) > (30) δ 2 ln Z (α) δα(x)δα(y ) .
1 δ 2 Z (α ) Z (α) δα(x)δα(y )
The coincidence case, x = y , is particularly interesting from the physical point of view. It is formally given by
2 ?2 r O (x) = ?r O (x, y )|x=y =
δ 2 ln Z (α) δ 2 α(x)
δ ln Z (α) δα(x)
Let us evaluate the second variation of ln Z (α). In general, in the coincidence limit one gets an ill de?ned quantity and a further regularization is required, as explained in the Appendix. Making use of (21), (13) and (63) in the Appendix, one has δ2 δ1 ln Z (α)(ε) =
N?2ε d ?2s lim × 2Γ2 (1 + ε) s→0 ds Γ(s) . (32)
dv v ε (u + v )s Tr e?uA δ2 Ae?vA δ1 A
For ε > 0 and su?ciently large, the integrand is regular at s = 0, and we have δ2 δ1 ln Z (α)(ε) =
∞ ∞ N?2ε ε dv v ε Tr e?uA δ2 Ae?vA δ1 A du u 2Γ2 (1 + ε) 0 0 ?2ε Tr A?1?ε δ2 AA?1?ε δ1 A . =N 2
This is our general formula for the second variation of the functional determinant. It is in agreement with a result obtained in . Let us consider the speci?c variation related to our observables. First, for O = φ2 , δL(α) = δα, and in the coincidence limit we have < ?2 φ2 > (ε) = 2N?2? ζ 2 (1 + ε|L)(x). (34)
For x = y , ζ (1 + ε|L)(x, y ) is regular at ε = 0, but in the coincidence limit y → x, ζ (1|L)(x) is well de?ned only for odd D and boundary free manifolds. ? In the conformally coupled case and for O = T? , one has δL(α) = c2 δαL. As a result
? 2 2 2 < ?2 T? > (ε) = 2Nc2 2 ζ (ε|L)(x) = 2Nc2 ζ (0|L)(x) + O (ε)
and we can remove the regularization parameter because ζ (s|L)(x) is regular at the origin. Some remarks are in order. If we limit ourselves to the odd dimensional case, in a manifold without boundary, the analytic continuation works in a simple manner for both ? the quantity O = T? and O = φ2 and no scale dependence appears in the ?nal expressions. As a consequence, it turns out that the relative variance, ?2 rO = is always equal to ?2 rO = and ? 2O = ? r 2 . N +2 (38) 2 , N (37) < O 2 > ? < O >2 , < O >2 (36)
The expressions above are valid for a multiplet of scalar N ?elds. In the case of other quantum ?elds, the formula for the relative variance could be a little bit di?erent (see for example ref.  for the treatment of gauge ?elds). In the even dimensional case, or/and in the presence of boundary, the situation is more complicated and a further renormalization seems unavoidable.
V. EXAMPLES A. D ?dimensional torus
In this Section we will consider, as ?rst example, the D ?dimensional torus. This is a symmetric ?at manifold with ?nite volume and the local zeta function is simply the ratio of the global zeta function and the torus volume. Thus, we may limit ourselves to the discussion involving the zeta function. This is a general conclusion valid for every symmetric space. The zeta function for the case of a massive (N = 1) scalar ?eld is given by 8
ζD (s|L) = with – 
4π 2 R2
ZD (s; (mR)2 )
ZD (s; (mR)2 ) = (2π )D/2 (mR)D?2s
Γ(s ? D/2) 2s/2+D/4+1 π s (mR)?s+D/2 + × Γ(s) Γ(s) √ s/2?D/4 n2 KD/2?s 23/2 πmR n2 .
Here R is the radius of the torus and m the mass of the ?eld. This expansion is exponentially convergent. It is to be seen not just as a big mass expansion (the convergence is then extremely fast), for it is valid in a very wide range of values of mR: 1 ? mR < ∞. In the case of a massless ?eld, the convenient expression to be used is quite di?erent (this is explained in detail in ): 21+s D?1? j/2 Z D (s ) = π Γ(s ? j/2)ζR (2s ? j ) Γ(s) j =0
n=1 nj ∈Z
On the other hand, for the sake of completeness, it is interesting to have a perturbative expression for the case when the mass of the ?eld is very small but di?erent from zero. This is obtained by means of binomial expansion in the equation de?ning the zeta function. The result is:
Kj/2?s 2πn n2 j ?.
ZD (s; (mR)2 ) =
Γ(s + k ) (mR)2k ZD (s + k ). k ! Γ(s)
This expression, combined with the preceding one, Eq. (41), yields the desired low mass expansion. Such expansions for the zeta function are not to be found in the literature. In fact, explicit expressions of the type (41) have appeared in the seminal paper  for the ?rst time. They are convenient, on the one hand, because they exhibit the pole structure of the zeta function explicitly. On the other, they are useful from the computational point of view, because they consist of a term including the main contribution, together with a series that converges extremely fast: only a few ?rst terms need to be computed in order to obtain results as accurate as desired. From these expressions, the ?rst and second variations of the e?ective action are immediate to compute. Essentially, what we get are expressions of the following kind: (i) in the massless case (p is a natural number, a small one for any of the operators here considered) Z D (p ) Γ(p ? 1)2 = × [ZD (p ? 1)]2 2p?2Γ(p)
D ?1 j =0 D ?1 j =0
π j/2 Γ(p ? j/2)ζR (2p ? j ) + 4π 2 SD,j (p)
[π j/2 Γ(p ? 1 ? j/2)ζR (2p ? 2 ? j ) + 4π 2 SD,j (p ? 1)] 9
where the SD,j (p) are fast convergent series providing only corrections to the main terms, (ii) ZD (p; (mR)2 ) = [ZD (p ? 1; (mR)2 )]2 Z D (p ) ZD (p + 1) Z D (p ) 1 + 2p ? (mR)2 + · · · , 2 2 [ZD (p ? 1)] [ZD (p ? 1)] [ZD (p)]2 (44)
in the case of a ?eld of very small mass (after doing a small-mass expansion as described above), and (iii) (mR)2p?4 Γ(p ? 1)2 ZD (p; (mR)2 ) = × [ZD (p ? 1; (mR)2 )]2 (2π )D/2 Γ(p) (45)
[Γ(p ? 1 ? D/2) + 2p/2?D/4+1/2 π p?1?D/2 (mR)p?1?D/2 + ΣD (p ? 1)]
Γ(p ? D/2) + 2p/2?D/4+1 π p?D/2 (mR)p?D/2 + ΣD (p)
in the massive case, where again ΣD (p) is a very fast convergent series. These expressions simplify very much when poles of the gamma function appear (for even dimension D ). They are quite easy to deal with. As advanced before, no singularity appears for D odd, and the whole expression remains, in that case (the series being well approximated by a couple of terms).
B. The Casimir slab
As a second example, we will revisit the computation of the local zeta function related to the Casimir slab, namely a massless scalar quantum ?eld con?ned in the x direction between in?nite parallel planar Dirichlet boundaries, located at x = 0 and x = a. In this case, the local zeta function is not trivial, due to the presence of the planar boundaries (see, for example [2,3,15,16]). To start with, recall the form of the local diagonal heat-kernel, which depends only on the con?ning coordinate x, i.e. < x|e?tL |x >= 2 a(4πt)
D ?1 2
nπx n2 π 2 t exp ? 2 ) . a a
Here L = ??2 is the Dirichlet Laplacian in the slab. Making use of the trigonometric duplication formula for the sine and Poisson-Jacobi re-summation formula, we may rewrite the above expression in the form < x|e?tL |x >= 1 (4πt) 2
∞ (na + x)2 a2 n2 )? exp (? ) . t t n=?∞
All the Seeley-DeWitt Kr (x) coe?cients vanish, but the ?rst one K0 (x) = 1. As a conse? quence, we may anticipate that < T? (x) > is vanishing. Making use of the Mellin transform and the above expression for the heat-kernel, and analytically regularizing the integral  10
dt tz = 0 ,
the analytic continuation for the local zeta function may be obtained. We present the result in a simple and symmetric form (for another equivalent form corresponding to D = 4, see ) ζ (s|L)(x) = × D (4π ) 2 Γ(s) [2ζR (D ? 2s) ? ζH (D ? 2s, x/a) ? ζH (D ? 2s, 1 ? x/a)] . ? s)a2s?D Γ( D 2 (49)
In the above expression, ζR (z ) and ζH (z, q ) are the Riemann and Hurwitz zeta functions respectively. This result can also be obtained by making use of the re-summation techniques explained in [35,15]. Some comments are in order. The local expression we have obtained is already in the renormalized form. Furthermore, the related trace involves a (in?nite, but trivial) volume VT in the transverse directions and an integration over x. Since
dq ζH (z, q ) = 0 ,
the resulting zeta function reads ζ (s|L) = ? s)a2s?D+1 2VT Γ( D 2 (4π ) 2 Γ(s)
ζ R (D ? 2 s ) ,
and this is exactly the zeta function associated with the Casimir slab con?guration. As far as the application to the evaluation of the vacuum expectation value is concerned, ?rst, let us consider D > 2. The zeta function is regular at s = 1 and the result is < φ (x) >= × D (4π ) 2 [2ζR (D ? 2) ? ζH (D ? 2, x/a) ? ζH (D ? 2, 1 ? x/a)] .
Γ( D ? 1)a2?D 2
The expression is ?nite everywhere and gives a vanishing result at x = 0 and x = a. Due to the simple geometry of the planar boundaries, the boundary divergences are not present. Things are di?erent in D = 2. In this case, the zeta function has a pole at s = 1 and the dependence on the scale ? appears. For the sake of completeness, we give the result. It reads < φ2 (x) >= a? πx 1 γ + ln . + ln sin 2π π a (53)
In this case, the boundary divergences are present. A few other exact analytic results obtained via zeta function regularization can be found in .
VI. CONCLUDING REMARKS
In this paper, we have revisited the use of zeta function regularization approach to the evaluation of expectation values of physical quantities and their related quantum variances. The former are associated with the evaluation of the ?rst variation of the e?ective action, while the latter are related to the second variation. We have shown that the ?rst variation can be regularized by use of zeta function techniques, and explicit expressions for the vacuum expectation values have been exhibited. For issues concerning the second variation, in general, zeta function regularization works well only when one is dealing with o? diagonal terms, since the coincident limit is highly singular. A further analytic regularization has to be introduced to treat the coincidence limit. ? We may summarize our results as follows. For a multiplet of N scalar ?elds and O = T? , in the conformally coupled case, analytic regularization gives a relative variance exactly 2 equal to N , independently on the dimension of the boundary free manifold. For O = φ2 , 2 , but now limited to the odd we have obtained again a relative variance exactly equal to N dimensional case without boundary. These restrictions can be removed in some particular situations, like the case of the Casimir slab, where one is dealing with a ?at manifold with ?at boundaries. The even dimensional case seems to require however further renormalization. The formalism can be directly applied to the expectation value of the stress tensor (see ), but only after much more e?ort and work. Again, the ?rst variation of the e?ective action with respect to the metric tensor gives < T?ν >. With regard to this issue, for an evaluation using local zeta function regularization see [28,36]. The second variation is related to the variance, and for the o? diagonal case, there is again factorization. In the special case of homogeneous space-times and for a conformally ? coupled scalar ?eld, one recovers the O = T? case.
In order to compute the variation of the trace of an elliptic invertible operator A, one has to take into account the fact that the variation (deformation) of A does not commute with A. From A?1 A = I , we have δA?1 = ?A?1 δAA?1 . 1 2πi (54)
For the calculation of complex powers of A we can use the Cauchy-Dunford representation A?s = ? 1 2πi
dzz ?s (A ? z )?1 ,
in which C is a suitable contour on the complex z ?plane. As a consequence, we have δA?s =
dzz ?s (A ? z )?1 δA(A ? z )?1 .
Making use of the two representations (for Re z > 0 and Re s > 0): Γ(s)(A ? z )?s = dtts?1 etz e?tA , (57)
and 1 1 =? Γ(z ) 2πi one gets δ TrA?s = ?sTr(A?s?1 δA) . (59)
dw (?w )?z e?w ,
If the operator A is self-adjoint, then there exist eigenvalues and eigenvectors, λn and Ψn , such that AΨn = λn Ψn , and we have δ TrA?s = 1 2πi 1 = 2πi = ?s
dzz ?s Tr((A ? z )?2 δA) dzz ?s
n ?s ?1 λn
< Ψn , δAΨn >= ?sTr(A?s?1 δA) .
1 < Ψn , δAΨn > (z ? λ n )2
For the second variation, making use again of (57) and (58), one gets δ2 δ1 TrA?s = ?sTr(δ2 A?s?1 δ1 A) 1 = ?s dzz ?s?1 Tr((A ? z )?1 δ2 A(A ? z )?1 δ1 A) . 2πi C
This expression is valid only if Supp δ1 has void intersection with Supp δ2 . Since we are interested in the coincidence limit y = x, we have to introduce an additional analytic regularization. The simplest one is to replace (A ? z )?1 with (A ? z )?1?ε and employ the representation in terms of the Mellin transform: Γ(1 + ε)(A ? z )?1?ε =
du uε euz e?uA ,
in Eq. (63) and then use again (58). As a result, we arrive to the formula δ2 δ1 TrA?s (ε) = ? 1 1 Γ2 (1 + ε) Γ(s)
dvv ε (u + v )s Tr e?uA δ2 Ae?vA δ1 A .
This representation for the second variation of the trace of a complex power of an elliptic operator A is in agreement with the one obtained in , making use of a di?erent method based on Schwinger’s perturbative expansion.
We are grateful to Valter Moretti for a useful discussion. This investigation have been supported by the cooperation agreement INFN (Italy)–DGICYT (Spain). EE has been ?nanced also by DGI/SGPI (Spain), project BFM2000-0810, and by CIRIT (Generalitat de Catalunya), contract 1999SGR-00257. 13
 N.D. Birrell and P.C.W. Davies, Quantum Fields in Curved Space (Cambridge University Press, Cambridge, 1982).  A.A. Grib, S.G. Mamayev and V.M. Mostepanenko, Vacuum quantum E?ects in Strong Fields (Friedmann Laboratory Publishing, St. Petersburg, 1994).  S.A. Fulling, Aspects of Quantum Field Theory in Curved Space-time (Cambridge University Press, Cambridge, 1989).  R.M. Wald, Quantum Fields in Curved Spacetime and Black Hole Thermodynamics (Chicago University Press, Chicago, 1994).  L.H. Ford, Phys. Rev. 47, 1 (1993); Phys. Rev. 51, 1962 (1995);  C.I. Kuo and L.H. Ford, Phys. Rev. 47, 4510 (1993).  G. Barton, J. Phys. A: Math. Gen. 24, 991 (1991); J. Phys. A: Math. Gen. 24, 5533 (1991).  C. Eberlein, J. Phys. A: Math. Gen. 25, 3015 (1992); J. Phys. A: Math. Gen. 25, 3038 (1992).  N. G. Phillips and B.L. Hu, Phys. Rev. 55, 6123 (1997).  D.B. Ray and I.M. Singer, Advances in Math. 7, 145 (1971).  J.S. Dowker and R. Critchley, Phys. Rev. D13, 3224 (1976).  S.W. Hawking, Commun. Math. Phys. 55, 133 (1977).  M. Kontsevich and S. Vishik. Functional Analysis on the Eve of the 21st Century. volume 1, 173–197, (1995).  Kirill A. Kazakow. Quantum ?uctuations of e?ective ?elds and the correspondence principle, hep-th/0207218 (2002).  E. Elizalde, S. D. Odintsov, A. Romeo, A.A. Bytsenko and S. Zerbini, Zeta Regularization Techniques with Applications (World Scienti?c, Singapore, 1994).  E. Elizalde, Ten Physical Applications of Spectral Zeta Functions (Springer Verlag, Berlin, 1995).  A.A. Bytsenko, G. Cognola, L. Vanzo and S. Zerbini. Phys. Rep. 266, 1 (1996).  B.L. Hu and N. G. Phillips, Int. J. Theor. Phys. 39, 1817 (2000).  N.G. Phillips and B.L. Hu, Phys. Rev. 62, 084017 (2000).  E. Toumboulis, Phys. Letters B70, 361 (1977); C.G. Callan, S.B. Giddings, J.A. Harvey, and A. Strominger, Phys. Rev. 45, 1005 (1992); D. Binosi, V. Moretti, L. Vanzo and S. Zerbini, Phys. Rev. 59, 104017 (1999).  T.P. Branson, P.B. Gilkey, K. Kirsten and D.V. Vassilevich, Nucl. Phys. B563, 603 (1999).  A.A. Bytsenko, G. Cognola, and S. Zerbini, Class. Quant. Grav. 14, 615 (1997).  A.A. Bytsenko, A.E. Goncalves, and S. Zerbini, Mod. Phys. Lett. A16, 1479 (2001).  A.A. Bytsenko, E. Elizalde, and S. Zerbini, Phys. Rev. D64, 105024 (2001).  G. Cognola and S. Zerbini, Mod. Phys. Letters 3, 599 (1988).  E. Elizalde, L. Vanzo and S. Zerbini, Comm. Math. Phys. 194, 613 (1998); E. Elizalde, G. Cognola and S. Zerbini, Nucl. Phys. B 532, 407 (1998); E. Elizalde, A. Filippi, L. Vanzo and S. Zerbini, Phys. Rev. D 57, 7430 (1998);  D. Iellici and V. Moretti, Phys. Lett. B425, 33 (1998).  V. Moretti, Commun. Math. Phys. 201, 327 (1999). 14
       
D. Binosi and S. Zerbini, J. Math. Phys. 40, 5106 (1999). E. Elizalde, J. Phys. A34, 3025 (2001). E. Elizalde, J. Comput. Appl. Math. 118, 125 (2000). E. Elizalde, Commun. Math. Phys. 198, 83 (1998). M. Bordag, E. Elizalde and K. Kirsten, J. Math. Phys. 37, 895 (1996). I.M. Gelfand and G.E. Chilov, Les Distributions, tome 1 (Dunod, Paris, 1972). A. Actor, Fortschr. Phys. 35, 793 (1987); Fortschr. Phys. 43, 141 (1989). V. Moretti, Journ. Math. Phys. 40, 33843 (1999).