High-temperature expansion of the one-loop effective action induced by scalar and Dirac particles

The complete nonperturbative expressions for the high-temperature expansion of the one-loop effective action induced by the charged scalar and the charged Dirac particles both at zero and finite temperatures are derived with account for possible nontrivial boundary conditions. The background electromagnetic field is assumed to be stationary and such that the corresponding Klein-Gordon operator or the Dirac Hamiltonian are self-adjoint. The contributions of particles and antiparticles are obtained separately. The explicit expressions for the $C$-symmetric and the non $C$-symmetric vacuum energies of the Dirac fermions are derived. The leading corrections to the high-temperature expansion due to the nontrivial boundary conditions are explicitly found. The corrections to the logarithmic divergence of the effective action that come from the boundary conditions are derived. The high-temperature expansion of the naive one-loop effective action induced by charged fermions turns out to be divergent in the limit of a zero fermion mass.


Introduction
The evaluation of the high-temperature expansions of the one-loop effective actions is a well elaborated procedure at the present moment. As for quantum electrodynamics (QED), the first results in this field date back to the papers [1,2]. Surprisingly, as far as we know, a uniquely defined (in the nonperturbative sense) expression for the complete high-temperature expansion of the one-loop effective action in QED in the Minkowski spacetime is not derived for a general stationary electromagnetic background. There are two issues in obtaining such a hightemperature expansion: i) To treat properly and nonperturbatively the external electric field; ii) To organize the high-temperature expansion in a manageable way that allows one to deduce its complete form. As for the latter problem, we may distinguish the works [3][4][5][6][7]. In [5][6][7], the high-temperature expansion was obtained but for the case A 0 = 0 and using the Euclidean approach. As for the former problem, we mention the papers and the book [8][9][10][11][12][13], where the leading terms of the high-temperature expansion for charged particles in the electric field were found. In the present paper, we shall derive the complete nonperturbative high-temperature expansion of the one-loop effective action induced by the charged scalars and the Dirac fermions in a stationary electromagnetic field of a general form in the Minkowski spacetime. We shall also find the contribution to this expansion coming from the nontrivial (MIT bag) boundary condition [14,15] imposed on the massive Dirac fermions.
When the external stationary electromagnetic field is purely magnetic, the Dirac equation can be squared to the equation of Klein-Gordon (KG) type with a self-adjoint operator. This facilitates the problem of finding the one-loop thermodynamic potential since the powerful methods based on the use of the spectral zeta-function of the Laplace type operators can be employed [3][4][5][6][7][16][17][18][19]. However, in the presence of the electric field, these methods cannot be immediately applied in the case of Dirac fermions because the naive squaring of the Dirac equation leads to the KG type equation with a non-Hermitian operator (see, e.g., [20][21][22]). Furthermore, this KG type operator possesses a kernel larger than the initial stationary Dirac equation. Therefore, we have to modify the approaches of [5-7, 13, 16-19] to cope this problem. The main idea we shall rest on consists in using the square of the Dirac Hamiltonian instead of the "square" of the Dirac equation operator (see (96)). This idea is not new and was employed, for example, in [23][24][25]. However, as for deriving the one-loop effective action, spectrum ε k (ω) of such an operator is discrete with an accumulation point at infinity [35][36][37]. Let us introduce the operator where the contour C runs from below upwards slightly to the left and parallel to the imaginary axis, and τ ν := |τ | ν e iν arg τ , arg τ ∈ [0, 2π). The operator (1) is trace-class when Re ν < 0. Let ω α k , α = 1, n(k), n(k) < ∞, be the real solutions to the equation ε k (ω α k ) = 0, and ε ′ k (ω α k ) = 0.
Then, the function ζ + (ν, ω) = Tr K −ν understood as a generalized function of ω, is analytic for Re ν < 1 and admits an analytic continuation to the whole complex ν plane. Let ϕ(ω) be the trial function, then where Tr θ(K(ω)) is a generalized function counting the number of positive eigenvalues of K(ω). Suppose that ε k (ω) is infinitely differentiable at the points ω α k , and the condition (3) holds true. Then, there is an asymptotic expansion when ε → +0. Consequently, according to the theorem 1, the generalized function ζ + (ν, ω) possesses singularities in a form of simple poles at the points ν ∈ N, and the function ζ + (ν, ω)/Γ(1 − ν) is an entire function of ν.
The first term in this expression is a contribution from particles, and the second one corresponds to antiparticles. The theorem 1 allows us to give a simple proof of the general formula for the high-temperature expansion of the Ω-potential [17,18]. Recall the main assumptions made about the zeta-function: 2. ζ + (ν, ω) is absolutely locally integrable on the ray [ω c , +∞) (when Re ν < 1); 3. There is an asymptotic expansion when ω → +∞: As was shown in [17,18] by a direct calculation, the third property holds for the operator K(ω) being a Fourier transform of a KG type operator. The coefficients ζ + k (ν) are expressed through the expansion coefficients of the heat kernel appearing in (1) (see, e.g., [13,17,18,24,25]). The second property follows from the definition (4).
Let us start with the case of the Bose-Einstein distribution and introduce the notation where µ ∈ (0, ω c ). Using the Mellin transform, we replace a nontrivial functional dependence on ω in the integrand by a simple power law where the contour C 1 runs parallel to the imaginary axis from below upwards slightly to the left of the point s = −1. Then, where Formula (12) is valid only when the contour C 1 is shifted to the domain of convergence of the integral representation for σ s ν (µ), i.e., when Re s < 2 Re ν − d − 1. In order to derive the high-temperature expansion of (12), we have to shift the contour C 1 to the domain of positive Re s and take into account the pole structure of the integrand. The location of the poles of the function Γ(−s)ζ(−s) is evident, and the structure of singularities of the function σ s ν (µ) can be determined by means of the theorem 1. The change of the variable, x = (ω − µ) −1 , transforms the integral to the desired form It is not difficult to find the asymptotic expansion of the integrand, where D := d + 1. All the theorem 1 conditions are satisfied. Therefore, the function σ s ν (µ) can be continued analytically to the domain Re s < 2 Re ν − d + N , where It possesses simple poles at the points s = 2ν − d − 2 + k, k = 1, N + 1, with the residues −a k−1 (µ). The implication of the theorem 1, guarantees the convergence of (12) and allows us to shift the contour C 1 to the right. Hence, we obtain where l 0 = ⌊2 Re ν − d + N ⌋. Notice that the term with l = −1 is understood as a limit, and the contour C 1 now runs slightly to the left of the line Re s = N + 2 Re ν − d.
If N → ∞, then setting aside an exponentially suppressed, as β → +0, term, we deduce from (18) the asymptotic expansion As for fermions, the considerations are quite analogous, but with the difference that the Riemann zeta-function ζ(z) should be replaced by the Dirichlet eta-function η(z) := (1 − 2 1−z )ζ(z) in all the appearances. Therefore, The term with l = −1 is absent in the expansion since the eta-function does not possess a singularity at z = 1 unlike the zeta-function. The contributions at the even positive l in the second term in (18), (19), and (20) vanish. With the aid of formulas (15), (16) for the singularities of σ s ν (µ), one can also check directly that expressions (19), (20) are finite when ν → 0.
The contribution from antiparticles to the high-temperature expansion of the Ω-potential is derived similarly. Let If the operator K(ω) is invariant under the simultaneous substitution ω → −ω, A 0 → −A 0 , then and Introduce the following functions understood in the sense of analytical continuation when Re s ≥ 2 Re ν − d − 1. Then, The proof of this formula is absolutely analogous to the one given above for particles.
3 Omega-potential of fermions at T = 0 In this section, we derive formula (9) of [19] for the nonrenormalized one-loop Ω-potential of fermions at zero temperature and nonzero chemical potential with the contribution of the energy of zero-point fluctuations. The formula may be of use when the solution of the Dirac equation reduces to the solution of KG equation with a self-adjoint operator.
There are two definitions of the energy of vacuum fluctuations for the Dirac fermions in the literature: the C-symmetric and the non C-symmetric one (see [29] for details). Let us start with the nonsymmetric one. According to the nonsymmetric definition (see, e.g., [30,31]), the energy of vacuum fluctuations of fermions is the energy of the "Dirac sea". Then, for fermions at zero temperature but nonzero chemical potential, where Λ is a cut-off parameter, E vac is a nonrenormalized vacuum energy, E part is an average energy of particles, and is the spectral density of energies of antiparticles. Therefore, where β 0 is some regularization parameter. The last equality is accurate within the renormalization ambiguity. Thus, where it is assumed that Tr θ(K(0)) = 0. In this case, Tr θ(K(µ)) for µ > 0 gives the number of states with a frequency 0 < ω < µ and, for µ < 0, it gives the number of states with a frequency µ < ω < 0. Assigning +1 to the charge of particles and −1 to the charge of antiparticles, we obtain i.e., we arrive at formula (9) of [19].
Applying the formulas of Sec. 2 to the expression (30), we deduce To justify the applicability of the formulas of Sec. 2, we have to assume that ζ + (ν, ω) = 0 in some neighborhood of ω = 0 and partition the integral in (32) into two: The consideration of Sec. 2 is to be applied to the second integral, while the first integral is an entire function of s tending to zero at | Im s| → ∞.
If we are interested only in the logarithmic divergence and the finite part of (30) as β 0 → 0, then where β 0 tends to zero after the limit ν → 0. Evaluating this limit, we come to where f.p. denotes the finite part of the expression as ν → 0. The expression (35) needs to be renormalized as, e.g., in (92).
When the C-symmetric definition for the vacuum energy is used (see, for example, [13,29]), Such an expression for the vacuum energy originates from the transformation of the initially Weyl-ordered Hamiltonian of fermionic fields to the normal form. Then, the vacuum energy of the Dirac fermions reads as where ν → 0. Taking the limit ν → 0 and keeping only the finite and divergent terms when β 0 → 0, we arrive at the nonrenormalized C-symmetric vacuum energy Notice that However, where T is an observation period, for the operator under the trace sign on the right-hand side is not trace-class. The left-hand side of (40) can be regarded as the definition for the right-hand side of (40) in the case of stationary background fields. The one-loop correction to the total nonrenormalized Ω-potential is written as The one-loop correction to the effective Lagrangian, L ef f = −Ω c tot , is real and corresponds to the standard choice of the vacuum state for stationary background fields.

High-temperature expansion for the Dirac fermions
In the previous sections, we have derived the formulas for the high-temperature expansion of the one-loop Ωpotential induced by fermions in the case when the spectral problem for the Dirac equation can be reduced to solving the KG type equation with a self-adjoint operator. This situation is realized in the absence of the external electric field (see, e.g., [19,33,[40][41][42] and references therein) and for certain special configurations of the electromagnetic fields (see [22]).
In this section, we shall derive the formula for the high-temperature expansion of the Ω-potential induced by fermions in the external stationary electromagnetic fields of a general configuration subject to the two restrictions: 1. The external electromagnetic field is such that the Dirac Hamiltonian is self-adjoint; 2. The Dirac Hamiltonian does not possess zero modes.
The first requirement says that the electromagnetic field does not have too strong singularities, and the electromagnetic potentials do not grow too fast at spatial infinity (if the problem is posed in the entire space). The second requirement can be violated in superstrong electromagnetic fields (see, e.g., [29] for details). This condition is rather technical and can be relaxed because the zero modes do not contribute to the one-loop thermodynamic potential of fermions (see (52)).
The eigenvalue eigenvector problem takes the form where ψ k is a Dirac bispinor and and We assume that the chemical potential µ conjugate to the electric charge is included into the definition of the potential A 0 (see, e.g., [43]). With this definition, the sign of the chemical potential is opposite to that used in Secs. 2, 3 and in [19]. The bispinor ψ k obeys the boundary condition [14,15,44] where Γ is a smooth boundary of the domain where the problem (42) is posed. The matrix Π is a projector, and we choose it in the form corresponding to the so-called MIT bag boundary condition [14,15] Π : where n µ is the inward unit normal to the surface Γ.
The spectral functions introduced allow one to evaluate the integrals of the form where ρ(ω) is the spectral density of the Hamiltonian H D , provided that f (ω) can be expressed in the form of the inverse Mellin transform with the contour C running downwards parallel to the imaginary axis at Re s < 2 − d. The integral (47) can be written in the form where we have substituted the Mellin transform of the spectral density It is assumed that the integral over s on the right-hand side of formula (49) converges. The representation (49) makes it possible to find the high-temperature expansion of the Ω-potential in terms of the spectral functions ζ s (ν) and η s (ν). Indeed, in virtue of the formula it follows from (49) that the contribution of particles to the one-loop Ω-potential can be written as where Re s < 2 − d.
In order to find the expansion in the rising powers of β, we need to move the contour C to the right up to the required power of β and take into account the singularities of the integrand in the complex s plane. In general, the integral along C moved to +∞ does not tend to zero. Therefore, when the contour C is moved to +∞, the resulting series in the rising powers of β is only asymptotic. The integral along C gives the remainder of this expansion and is exponentially suppressed when β → +0.
Let us investigate the singularities of ζ s (ν) and η s (ν) in the ν plane and their behavior when | Im ν| → ∞ (see, e.g., [35,45]). These spectral functions can be expressed in terms of the heat kernel trace It is useful to introduce [45] Then Hence, it is sufficient to study the analytic properties of ζ s (ν, ε) at ε → 0.
Since (H D − ε) 2 is a Laplace type operator (see (81)), the following asymptotic expansion holds for τ → +0 (see [44] for a review) where a k/2 (ε) are the heat kernel expansion coefficients. The integrand of (54) is absolutely integrable. On the upper limit, the integral over τ in (54) is convergent for any ν because ker H 2 D = 0. Therefore, introducing some cutoff Λ on the upper integration limit in (54), we can apply the theorem 1. As a result, we obtain the structure of singularities of ζ s (ν): where a k/2 := a k/2 (0). From (55) we have where a ′ k/2 := ∂ ε a k/2 (ε) at ε = 0. It follows from these formulas that The theorem 1 also implies Consequently, when | Im ν| → ∞, the functions |ζ s (ν)| and |η s (ν)| tend to infinity not faster than e | Im ν|(π/2+0) . Notice that, in the case we consider, a stronger estimate even holds (see the remark after theorem 5.5.2 in [35]). This means that the integral over s on the right-hand side of (52) converges, and the contour C can be moved to the right. Now we see that the integrand of (52) has the pole singularities at the points s = −d, ∞. The poles at Re s ≤ 0 and s = 2l + 2, l = 1, ∞, are simple. The other ones are of the second order. Let us introduce the notationζ The functionsζ s (ν),η s (ν) are the finite part of the Laurent series (i.e., the coefficient c 0 ) of the functions ζ s (ν), η s (ν) in the vicinity of the point ν. Then the high-temperature expansion of the contribution of particles to the one-loop Ω-potential is where D := d + 1, and we have taken into account that a ′ 0 = 0. As long as the contribution of antiparticles has the same form as (62) Therefore, the total Ω-potential does not contain the functions η s (ν),η s (ν), and the heat kernel expansion coefficients a ′ k . Bearing in mind (59), we see that the asymptotic expansion in β of the total Ω-potential for fermions without the vacuum energy contribution is expressed through the heat kernel expansion coefficients save the term proportional toζ s (−1/2).
Let us step back for a while and consider (52) once again. If we substitute the representation (53) into (52) and change the order of integration over s and τ , the following integral arises This integral is reduced to the Jacobi theta-function. Using the representation [46] where H is the Hankel contour and the principal branch of the power function is taken, it is not difficult to obtain that Developing the logarithm as a series and integrating it term by term, we arrive at where ω n := π(2n + 1)/β are the Matsubara frequencies for fermions. Hence, the contribution of ζ s (1 − s)/2 to (52) is This result is somewhat expected as the evaluation of the path-integral for fermions over ψ † and ψ instead ofψ and ψ gives at finite temperature (see, e.g., [47,48]) rather than the logarithm of the determinant of the square of the Dirac equation operator. We shall not use the representation (69) below, but it can be useful for the numerical evaluation of (52).
As we have already noted in Sec. 3, one can introduce the two definitions for the energy of vacuum fluctuations of the Dirac fermions: the C-symmetric (36) and the non C-symmetric (26) one. In accordance with the non C-symmetric definition, where β 0 is the regularization parameter and E are the antiparticle energies. Using (59) and (62), we obtain the nonrenormalized expression for the vacuum energy where one should set µ = 0. The first line contains the power divergencies only. The logarithmic divergence responsible for the conformal anomaly and the finite part are presented in the last line. The C-symmetric vacuum energy E c vac is defined in (36). The explicit expression for E c vac is obtained from (72) by throwing out all the terms containing a ′ k/2 . The symmetric definition of the vacuum energy seems to be a more preferable one. It is symmetric with resect to the exchange of notions of a particle and an antiparticle. Besides, this definition implies that the finite part of the vacuum contribution is exactly canceled by the corresponding contribution from the thermal part of the Ω-potential in the high-temperature limit and at zero chemical potential. We see from (72) that, in the four-dimensional spacetime, the finite parts of the symmetric and nonsymmetric vacuum energies differ only by the surface term proportional to a ′ 5/2 . As it should be, the energy of vacuum fluctuations (26), (36) and the corresponding effective action are real-valued for the fields satisfying the conditions mentioned at the beginning of this section. The vacuum state is standardly defined (see, e.g., [29] for details) as the ground state of the Hamiltonian of quantum Dirac fields. The particles are defined as the modes with ω k > 0 (see (42)), whereas the antiparticles are the modes with ω k < 0. The total Ω-potential with the contribution of the vacuum fluctuations becomes The high-temperature expansion of Ω th is presented in (62), where the contribution from antiparticles should be taken into account, and E ren vac is the renormalized vacuum energy (see (92)). The termζ s (−1/2) in Ω th at zero chemical potential is exactly canceled by the analogous term in E ren vac [9,10,[26][27][28]. Hence, the total asymptotic expansion in β of the one-loop effective Lagrangian L (1) ef f induced by the Dirac fermions is expressed only in terms of the heat kernel expansion coefficients, when µ = 0. This property is an analog of the fact that the terms with σ l ν , for l even, vanish in the high-temperature expansion (18). Notice that, for a nonzero chemical potential, the cancelation ofζ s (−1/2) in L (1) ef f is not complete (see [27] and (93)). The contributions ofζ s (−1/2),η s (1/2 − l), l = 0, ∞, in (62) are not expressed through the heat kernel expansion coefficients. Nevertheless, as is well known, the heat kernel expansion allows one to find the asymptotic expansion of these contributions in the inverse powers of a large mass m. Here we shall find such an expansion only forζ s (−1/2) andη s (1/2) appearing in the finite at β → 0 part of the high-temperature expansion. Let and Then, substituting (74) into this expression and developing it as a series in τ with the aid of (56), we deduce Substituting (74) into (54) and using the expansion (75), we come tō where f.p. means the finite part with respect to ǫ → 0, the relation (76) has been used in the last line, and the prime at the sum sign recalls that all the terms that are singular due to the gamma function must be omitted. Employing the relation [46] ψ(1 + n) = −γ + H n , H n := n k=1 k −1 , we have eventuallȳ As far asη s (1/2) is concerned, it is obtained fromζ s (−1/2) by means of (55). Substituting (79) into (72), we derive the asymptotic expansion where µ = 0. The last term in this expression looks as the standard expansion of the one-loop effective action with respect to a large mass (see, e.g., [21,49]). However, one should bear in mind thatā k/2 depend on m. The maximal power of m appearing inā k/2 will be found later in Sec. 5. The self-adjoint positive-definite Laplace type operator defining the heat kernel (54) takes the form where P 0 := ε − A 0 and D i = P i − α i P 0 . The operator iγ 0 γ 5 is self-adjoint in the Hilbert space specified by the condition (44). Inserting (iγ 0 γ 5 ) 2 under the trace sign in (54), we infer that for ε = 0. Remark that the derivative with respect to ε in (55) and in a ′ k/2 at ε = 0 can be replaced by the minus derivative with respect to the chemical potential µ entering into the definition of A 0 . Therefore, we put henceforth P 0 = −A 0 in (81).
The operator γ 5 is not compatible with the boundary condition (44) and does not map the Hilbert space into itself. Nevertheless, if we consider the spectral problem (42) in the entire space, i.e., if we do not impose the boundary condition (44), then γ 5 becomes a self-adjoint operator in the respective Hilbert space. In this case, inserting (γ 5 ) 2 under the trace sign in (54), we see that the spectral density ρ(ω) is symmetric with respect to the replacement m → −m. In the next section, we shall explicitly show that the boundary condition (44) violates the symmetry of the spectral density under m → −m.

Leading terms of the high-temperature expansion
Making use of the general formulas of the previous section, let us obtain the leading terms of the high-temperature expansion of the Ω-potential of Dirac fermions in a bag. The first few coefficients of the heat kernel expansion associated with the Laplace type operator with nontrivial boundary conditions are given in [44,50]. The higher coefficients without the surface terms are presented in [51]. The requirement of self-adjointness of the Dirac Hamiltonian with the boundary condition (44) leads to whereP = γ µ P µ and We have introduced here the extrinsic curvature of the hypersurface Γ where the boundary condition is defined (see, for details, [13]): L ab := −e µ a e ν b ∂ µ n ν + n λ Γ λ µν e µ a e ν b = −e i a e j b ∂ i n j , a, b = 1, 2, where e µ a are orthonormal vectors tangent to the surface, n µ = (0, n) is the inward unit normal, and Γ λ µν is the Levi-Civita connection (in our case it is equal to zero). Hence, employing the notation from [44,50], we have from (81) and (83): The expressions for the traces arising in evaluating the heat kernel expansion coefficients are given in Appendix B. Substituting the traces found in (109)-(112) to the general formulas 1 of [44,50], we obtain 1 Notice that we did not check these formulas.
where τ are the coordinates on the boundary surface, and h is the determinant of the metric induced on Γ. Notice that, in a 2 , the last term in the first line is canceled by the first term in the second line. The last term in a 2 can be omitted. As for a 5/2 , the terms that do not vanish under differentiation with respect to µ are only presented. The volume terms in the coefficients a k , k = 0, 2, coincide with the "pseudo-trace" expansion coefficients ([13], problem 7.21, see also (100)). The derivatives with respect to −µ are written as As we have discussed above, the last expression gives the difference between the finite parts of the C-symmetric and the non C-symmetric vacuum energies of Dirac fermions. Upon renormalization, this contribution has to be completely canceled out by the counterterms since the integrand in a ′ 5/2 has the dimension not exceeding 3 and is not Lorentz-invariant, even if one restores the dependence on the normal vector n µ . Thus, after renormalization, the C-symmetric and the non C-symmetric vacuum energies coincide.
The leading correction to the high-temperature expansion due to the nontrivial boundary condition (44) comes from the coefficient a 1 . When the boundary Γ is a sphere of the radius R, the extrinsic curvature L ab = δ ab /R. In this case, substituting (87) into (62) doubled and setting m = 0, we reproduce the leading correction in (48) of [32] to the free energy due to the boundary condition (44). As we see, the coefficient a 1 and the corresponding contribution to (62) are not symmetric with respect to m → −m. Consequently, the spectral density ρ(ω) is not symmetric under this replacement either. The volume contributions to a k/2 are invariant with respect to m → −m.
The coefficient a 2 is prefixed to ln β 0 in the vacuum energy (72) and, as is well known, related to the conformal anomaly. Formula (87) provides the generalization of the standard expression for a 2 to the case of the electrically charged Dirac fields obeying the boundary condition (44). In the case of the neutral massive Dirac field confined to a sphere, the surface term coming from a 2 coincides with (3.1) of [23]. The surface term in a 2 is not symmetric under m → −m. The term mA 2 0 in the surface contribution to a 2 should be completely canceled out by the counterterm inasmuch as it is not Lorentz-invariant.
In order to renormalize the other contributions to the vacuum energy, we need to find the terms inā k with the mass dimension less than or equal to 4, the dimension of the coefficients depending on m being not taken into account. So we have to determine the maximal power of m entering intoā k/2 . As for the volume terms in a k/2 , such contributions come from the second term in E. The maximal power of m in the surface terms stems from the first term in S. In this case, Therefore, if one takes into account the surface terms, the expansion in m −2 in the last contribution in (79) and (80) is, in fact, not an expansion in the rising powers of m −1 . In order to obtain the large mass expansion with the surface contributions, one has to resum the expansion in (79) and (80). We leave this issue for a future research and, in considering the large mass expansion, shall cast out the surface terms. The expansion obtained thereby is valid for the problem (42) posed in the whole space without the boundary condition (44). The necessary for the renormalization procedure terms arē The other volume terms inā k/2 , k ≥ 3, are suppressed in (79), (80) by a power of a mass. The asymptotic large mass expansion of the vacuum energy (80) without the surface terms takes the form whereÃ 0 := A 0 | µ=0 . In accordance with the standard renormalization rules (see, e.g, [52]), all the written terms in (91) should be completely canceled out by the counterterms provided one normalizes the effective action at a zero photon momentum and a zero chemical potential. As a result, we have from (72): Substituting this expression into (73), we find that the complete asymptotic high-temperature expansion of the total one-loop Ω-potential of the Dirac fermions with account for the contribution of antiparticles and the vacuum energy (but without the surface terms) reads as . (93) Let us note once again that, at µ = 0, this expansion is expressed solely in terms of the heat kernel expansion coefficients a k , and this is not a large mass expansion. Because of the logarithmic contribution, this expression is not defined when m → 0, i.e., there exists an infrared divergence in this limit, at least at zero chemical potential.
One can try to solve this problem by introducing an effective mass for the fermions and sum thereby an infinite number of the ring diagrams (see, e.g., [48,53]). The contribution of the photons to (93) with their effective mass should be also taken into account. The leading terms of the high-temperature expansion of the Ω-potential (62) expanded in a large mass m without the contribution of antiparticles and the surface terms are written as For E = 0, this expression coincides with the leading terms of the expansion (105) of [19] provided one changes the sign of the chemical potential. If we take into account the contribution of antiparticles, i.e., if we take the doubled symmetric in A 0 part of (94), then we obtain formula (4.2) of [10] up to a renormalization ambiguity. The nonsymmetric in A 0 part of (94) can be used, for example, to find the number of electron-positron pairs in the system (see, for details, [19,33]) at finite temperature.

Comparison with other approaches
Let us compare the method we have developed with the approaches presented in the literature. Of course, we do not pretend here to a complete description of the whole literature devoted to this subject, but mention only the main approaches that we know and that are close to the method used by us. Usually, the one-loop correction to the effective action is expressed in terms of the trace [21,49,54] Tr 4 e −τ K , where The quantity (95) cannot be directly used for a nonperturbative evaluation of the effective action since the operator under the trace sign is not trace-class [55,56]. Moreover, any operator function f (K) different from zero is not trace-class. In order to see this, it is sufficient to consider the trace in the momentum representation in the domain of large |p µ |. Then where K(p, x) is the principal symbol of K, which is obtained from K by the replacement of i∂ µ by p µ and neglecting the terms of the lower order in p µ . The integrand of (97) is invariant under the Lorentz transformations of the vector p µ preserving the spacetime metric g µν (x) at the fixed point x. The action of the Lorentz group foliates the momentum space onto the noncompact orbits where the integrand is constant. Therefore, the integral (97) is not absolutely convergent and can converge only conditionally what implies that the operator f (K) is not trace-class. Perturbatively, this property is revealed as the conditional convergence of the renormalized Feynman diagrams [57,58] in the Minkowski spacetime. For the electromagnetic field of a general configuration, the expression (95) seems to be understood only perturbatively, and the trace should be evaluated in a certain basis (see, for example, (40)).
In the case of the stationary electromagnetic potential A µ , one can consider the operator and pose the so-called nonlinear spectral problem for it (see the book [13] and references therein). Then, formally, the evaluation of the effective action at finite and zero temperatures is reduced to the evaluation of the expression (see Secs. 2, 3 above) Tr e −τ K(ω) .
In [13], Sec. 6, the pseudo-trace, is introduced. Here ω k are the solutions to the equation and λ k (ω) are the eigenvalues of the operator ω 2 − K(ω). The high-temperature expansion of the Ω-potential is expressed through the expansion of (100) for small t. As long as the operator (98) is not self-adjoint, equation (101) can have complex roots. The methods presented in Sec. 6 and problem 7.21 of [13] (see also [39]) are not applicable, at least immediately, in this case since they are essentially based on the self-adjointness of the operator ω 2 − K(ω), i.e., on the real-valuedness of λ k (ω) for real ω (see the remark after (6.6) of [13] and the requirements ii), iii) after (2.7) of [39]). Besides, the physical interpretation of the complex solutions of the nonlinear spectral problem for the operator (98) is unclear. The kernel of the operator (98) can be larger than the set of stationary solutions of the initial Dirac equation. Unfortunately, this issue is poorly discussed in the literature. In [54], this problem is not touched at all. In the standard textbook [29], it is wrongly assumed that K is a self-adjoint operator (see the remark after (10.108b)). In [59], Sec. 6.2.3, it is supposed that K possesses a real-valued spectrum and a set of eigenvectors that is complete and orthonormal with respect to the Lorentz-invariant scalar product. In [60], it is assumed, in considering the thermodynamic properties of fermions in external electromagnetic fields, that the electric field is completely screened, and actually the case of E = 0 is investigated. A more detailed discussion of this issue is presented in [21], p. 551, but its solution is not given there. Besides, it is supposed that the imaginary parts of the eigenvalues of the operator K are of the same sign. However, in the case A = 0, it is easy to show using the Majorana representation of γ-matrices (see, e.g., [22]) that if ε k and ψ k are the eigenvalue and the eigenvector of (98), then ε * k and γ 0 ψ * k are also the eigenvalue and the eigenvector of (98) (we do not impose the boundary condition (44)). In this representation, equation (98) with the potential A 0 (x 2 ) becomes where E 2 = −∂ 2 A 0 . We see that, for the function A 0 (x 2 ) of a general form, the operator (98) possesses the complex eigenvalues ε k . Indeed, let ε k ∈ R for some A 0 (x 2 ). Then, perturbing the potential by some δA 0 (x 2 ), we obtain with the help of the standard formula of perturbation theory that δε k ∈ C. Nevertheless, the methods based on the use of the expression (99) are nonperturbatively applicable for some classes of the electromagnetic fields. For example, such electromagnetic fields are 1. The fields satisfying the condition F µν n ν = 0, where n µ is a constant isotropic vector. The solution of the Dirac equation with these fields is reduced to the solution of the KG type equation with a Hermitian operator (see, for details, [22]); 2. The fields with A 0 = 0. In this case, the operator (98) is self-adjoint, and (81) is reduced to ω 2 − K(ω) when ε = 0. Therefore, the approaches expounded in Secs. 2 and 4 are equivalent. The vacuum energy and the Ω-potential written in terms of the spectral zeta-function associated with the squared Dirac equation operator are presented in this case in many papers and books (see, e.g., [19,33,[40][41][42] and references therein); 3. The constant electromagnetic field. In this case, the last term in (98) can be factored out from the exponent in (99), and the problem is reduced to the respective problem for the KG equation with a Hermitian operator.
There is also a third approach to the problem based on the Wick rotation prescription (see, e.g., [5][6][7][61][62][63] for details). According to this method, the one-loop correction to the effective action induced by Dirac fermions formally reads as The gauge connection A E µ included into P µ is real. The chemical potential enters (103) as the imaginary additive iµ to A E 0 . The spectrum ofP E is real-valued for µ = 0 and, in the even-dimensional space, is symmetric with respect to zero. This follows from the fact that there exists the matrix γ 5 in the even-dimensional space that maps the eigenvectors associated with the positive eigenvalues to the eigenvectors associated with the negative eigenvalues. It is assumed in (103) that the determinant is regularized in such a way that this symmetry of the spectrum ofP E is preserved. The determinant on the right-hand side of (103) can be evaluated, for example, with the help of the spectral zeta-function. Then, in the even-dimensional space, the multiplicative anomaly arises (see, for details, [7,64]) which leads to the additional corrections to the expression (103). Unfortunately, in the case A 0 = 0, we do not know a rigorous proof of the equivalence (in the nonperturbative sense) of this approach to the Minkowski spacetime methods presented, for example, in Secs. 2, 4 above, or in [13]. For A 0 = 0, it was shown in [7] that these approaches are equivalent for ultrastatic spacetimes. Notice that the Euclidean approach developed in [13] is distinct from that is described here. In [13], having performed the Wick rotation, the background fields become complex (see Sec. 7.4, [13]), and the corresponding wave operator becomes non self-adjoint even in the massless case.
All the methods mentioned above must be perturbatively equivalent to the approach developed in Sec. 4 even in the case A 0 = 0. In other words, the asymptotic expansions of the effective action in the coupling constant or, what is the same in the case at issue, in the inverse powers of a large mass must coincide up to a renormalization ambiguity. This is fulfilled, in particular, if a k/2 coincide with the expansion coefficients of the pseudo-trace (100) associated with the operator (98) (see (80)). In Sec. 5, we showed that this relation holds for k = 0, 2, 4 without the surface terms. However, we have not succeeded in proving this relation for an arbitrary k. The perturbative equivalence to the Euclidean approach entails that the coefficients of the large mass expansion in the last line of (80) multiplied by T , where T is the observation period, are equal to the heat kernel expansion coefficients constructed for the operatorP 2 E , where the inverse Wick rotation and renormalization are performed. The direct proof of this statement deserves a separate study. Let us stress that, at first,ā k/2 in (80) must be expanded in m −1 .
Notice that the method of deriving the high-temperature expansion that employs the squaring of the Dirac Hamiltonian was used in [24,25] in the case of a stationary gravitational background. The background connection associated with the internal gauge symmetries was not considered in [24,25]. Besides, in [13,24,25], the expression for the complete high-temperature expansion was not derived. The method we have developed in the previous sections allows us to obtain the complete high-temperature expansion of the contributions of particles and, separately, of antiparticles to the Ω-potential, the nontrivial boundary conditions being also taken into account. This expansion is expressed in terms of the spectral zeta-function of the positive-definite selfadjoint operator of Laplace type. In [13,24,25,39], the contributions from the boundary conditions and the separate contributions from particles and antiparticles were not considered. In (1.66) of [7], the high-temperature expansion was obtained in terms of pseudo-differential operators in the case A 0 = 0.
In [9,10], the approximate expression for the one-loop thermodynamic potential of the Dirac particles was found in a stationary external electromagnetic field. In this expression, the zeroth and the first derivatives of A µ are only taken into account. Several leading terms of the high-temperature expansion were also explicitly found in [9,10]. Unfortunately, the rigorous nonperturbative derivation of formula (2.1) in [10] is absent. In [9], the reference is made to the paper [54], so the approach of [9,10] appears to be equivalent to the method based on the expression (99) that we have discussed above. The expression (4.2) of [10] for the first leading terms of the high-temperature expansion coincides with (94) up to renormalization provided that the doubled symmetric in A 0 part of (94) is taken. In [11,12,65] (see also [8]), a systematic method of deriving the gradient corrections to the effective action both at zero and finite temperatures was elaborated, the Schwinger's answer for the constant electromagnetic field generalized to a finite temperature being taken as the zeroth order approximation. According to this approach, the vacuum energy contains an imaginary contribution when the electric field is present. It is not clear whether the expression (1) of [12] is well-defined nonperturbatively or not.

Conclusion
To sum up, we have obtained an explicit well-defined nonperturbative expression for the complete high-temperature expansion of the one-loop Ω-potential induced by scalars and Dirac fermions with nontrivial boundary conditions. The contributions of particles and antiparticles were treated separately that allows one to find the number of pairs of particles and antiparticles in the system (see, for details, [19,33]). The high-temperature expansion can be used to obtain the one-loop energy of vacuum fluctuations at zero temperature both for bosons and fermions. So we have obtained the vacuum energies for bosons [17][18][19] and fermions. In particular, we have found the explicit expression for the difference between the C-symmetric and the non C-symmetric vacuum energies (see, for details, [29]) for the Dirac fermions. It turns out that the finite parts of these energies differ only by a surface term, which is absent if the Dirac fields obey the standard asymptotic conditions at spatial infinity, i.e., the condition (44) is not imposed. Thus, in this case, we may conclude that both definitions of the vacuum energy give the same result upon renormalization of the divergent part. In a general case, we have shown that this finite surface term should be completely canceled out by the respective finite counterterm. We have also found the leading correction to the high-temperature expansion due to the nontrivial boundary condition (44). Besides, the correction to the coefficient a 2 responsible for the conformal anomaly has been derived in the case of the massive charged Dirac fields obeying the MIT bag boundary condition (44). Finally, we have proved that the complete asymptotic in β high-temperature expansion of the one-loop effective action induced by the Dirac fermions at zero chemical potential is expressed solely in terms of the heat kernel expansion coefficients even in the case when the mass of the field is not large.
The difference between the C-symmetric and non C-symmetric definitions of the vacuum energies can be found for quarks in a bag. However, at first, one should generalize the above considerations to the case of non-Abelian background gauge fields. Another possibility to examine experimentally these two definitions can be realized in a graphene sheet with boundaries. In this case, the non C-symmetric definition is even more preferable from the physical viewpoint as the "Dirac sea" (the electrons in the valence band) is really present. The Lorentz symmetry is not a fundamental symmetry in this case either. Therefore, one may impose different normalization conditions than that we have used in Sec. 5.
Several questions remain unanswered and are left for a future research. First, we did not prove by an explicit calculation that, up to a renormalization ambiguity, the large mass expansion of the vacuum energy (80) coincides with the standard answer given, for example, by the Euclidean approach. Second, we did not obtain the large mass expansion of the vacuum energy with the surface contributions taken into account. The naive heat kernel expansion is to be resummed in this case. The third problem, which seems to be the most difficult one, is to prove (or disprove) the nonperturbative equivalence of the different approaches to the evaluation of the oneloop effective action at zero and finite temperatures. Of course, the procedure developed admits immediate generalizations to the curved spacetime, the chiral gauge connections, and other space dimensions. We leave these generalizations for a separate study.

A Some analytic properties of the Mellin transform
Let us consider the analytic properties of the integrals of the form as functions of the complex variable λ. The following theorem holds for such integrals (see, e.g., [34,66]). Proof. The integrals I λ and ∂ λ I λ are absolutely convergent for Re λ > −1, and, consequently, I λ is an analytic function of λ in this domain. In this domain, we have The right-hand side of this equality is defined for Re λ > −2, λ = −1, and provides an analytic continuation of I λ to this domain. The function I λ has a simple pole at the point λ = −1 with the residue a 0 . Continuing this process by subtracting the terms of the asymptotic series (106) from ϕ(x), we arrive at the first assertion of the theorem. For Re λ > −2 − N , we have Making the substitution y = ln x, we reduce the integral in this expression to the Fourier transform. In the domain of the λ plane under consideration, the integrand is absolutely integrable, and, consequently, according to the Riemann-Lebesgue lemma (see, e.g., [67]), the integral in (108) tends to zero for | Im λ| → ∞. The out of the integral terms in (108) also tend to zero when | Im λ| → ∞ and Λ > 0.
In many cases, this theorem allows one to investigate the analytic properties of the function I λ without an explicit evaluation of the integral (105). If the upper limit of the integral (105) is Λ = +∞, then one can split the integral into two ones and reduce them to (105) by a change of the integration variable, x → x −1 , in the integral with the infinite integration limit. If the asymptotic expansion (106) is performed in the fractional powers of x, then one can bring the integral at hand to the form (105) by a suitable change of the integration variable x and redefinition of λ. Differentiating I λ with respect to λ, one can generalize the theorem to the case when the asymptotic expansion of ϕ(x) contains the terms x k ln l x. In this case, I λ possesses multiple poles in the complex λ plane. Theorem 1 applied to the spectral density of the self-adjoint operator A > 0 gives the structure of singularities of the spectral zeta-function ζ(ν, A) in the complex ν plane and relates the residues of the spectral zeta-function to the heat kernel expansion coefficients.
(112) In the expressions (110)-(112), we retain only the terms that give nonzero contributions to the derivative with respect to the constant chemical potential µ included into A 0 .