# Boundary divergences in vacuum self-energies and quantum field theory in curved spacetime

###### Abstract

It is well known that boundary conditions on quantum fields produce divergences in the renormalized energy-momentum tensor near the boundaries. Although irrelevant for the computation of Casimir forces between different bodies, the self-energy couples to gravity, and the divergences may, in principle, generate large gravitational effects. We present an analysis of the problem in the context of quantum field theory in curved spaces. Our model consists of a quantum scalar field coupled to a classical field that, in a certain limit, imposes Dirichlet boundary conditions on the quantum field. We show that the model is renormalizable and that the divergences in the renormalized energy-momentum tensor disappear for sufficiently smooth interfaces.

###### pacs:

03.70.+k;11.10.Gh## I Introduction

The vacuum energy produces measurable forces between neutral bodies reviews (). But even in the presence of a single body, the modes of the electromagnetic field are disturbed and produce a self-energy. Although irrelevant for the calculation of forces between different objects, this self-energy is in principle observable through its coupling to gravity.

In the case of perfect conductors, the vacuum energy density, or more generally the energy-momentum tensor, diverges near the boundaries, as noted for the first time by Deutsch and Candelas a long time ago DC (). These divergences are not the usual ones in quantum field theory, because they are present in the already renormalized energy-momentum tensor. The origin of the divergences is the unphysical assumption of perfect conductivity for all modes of the electromagnetic field since, on physical grounds, one expects any material to become transparent at high energies. More generally, divergences in the renormalized energy-momentum-tensor are present even for non-perfect conductors, as long as there is a sharp boundary between two media with different electromagnetic properties. In this case the origin of the divergences is that, for modes with extremely small wavelengths, it is unphysical to assume a sharp boundary, since the transition region between media becomes larger than the wavelength of the high frequency modes.

The problem of the divergences in the self-energies has been considered in flat spacetime in a series of works by Graham et al MIT (). It was pointed out that, if the boundary conditions are replaced by an interaction with a second field, the model is renormalizable and the self-energies are not divergent. In those works, the authors considered a toy model consisting of a vacuum scalar field coupled to a background (classical field) , with an interaction of the form . The sharp limit corresponds to consider a discontinuous , and the “perfect conductor” limit to take at a particular point, imposing Dirichlet boundary conditions on the vacuum field at that point. When these limits are taken, the self-energy depends on the ultraviolet cutoff, i.e. on the specific interaction between the high energy modes of the quantum field and the microscopic degrees of freedom of the body. Recently, Milton Milton () and Bouas et al Bouas () considered a similar problem, computing the energy density for a scalar vacuum field in some particular potentials (“soft walls”), showing that the energy density near a potential barrier is finite for sufficiently smooth potentials. For a discussion of some aspects of the coupling of the Casimir energy to gravity see Ref.revmilton ().

A complete analysis of the Casimir self-energies and their eventual gravitational implications must be performed in the context of quantum field theory in curved spacetimes. Following Ref.MIT (), one could replace the boundary conditions by interactions with a classical field, and show not only that the matter sector of the theory is renormalizable, but also that the usual divergences in the energy-momentum tensor can be absorbed in the coupling constants of the gravitational sector, resulting in finite and well defined semiclassical Einstein equations when the classical field is sufficiently smooth. The divergences associated to the unphysical limits will reappear when considering discontinuous classical fields. The aim of the present paper is to provide such analysis. We will consider a scalar vacuum field in curved spacetimes, coupled to a classical field that models the “mirror”. We will show that, for smooth , the model is renormalizable using the standard renormalization procedure for quantum fields in curved spacetimes. We will provide some examples to illustrate the appearance of divergences for sharp interfaces.

It is worth remarking that similar divergences in the renormalized energy-momentum tensor do appear even for free fields in curved spacetimes, if one assumes that the spacetime metric is not sufficiently smooth. A well known example, in cosmology, is the divergence in the energy density that appears when one considers an abrupt transition between two different epochs (e.g. De Sitter space to radiation dominated universe Ford ()), that involves a discontinuity in the curvature tensor. A less known example is the divergence of the vacuum polarization around spherically symmetric objects, at a sharp interface between the object and the vacuum stars ().

The paper is organized as follows. In the next section we introduce the model. In Section 3 we prove that the divergences in the energy-momentum tensor can be absorbed into the bare constants of the theory, yielding finite semiclassical Einstein equations. In Section 4 we present some explicit calculations for weakly coupled mirrors, showing that as long as the potential that models the mirror is smooth enough, there are no divergences in the renormalized energy-momentum tensor and in the renormalized coincidence limit of the two point function, . On the other hand, discontinuities in the potential or in its first two derivatives produce infinite answers. Section 5 contains a discussion of the main results. We will work with natural units and metric signature .

## Ii The Model

We consider a quantum vacuum field interacting with a background classical field on a curved spacetime. The action of the complete system is

(1) |

where

(2) |

and

(3) |

Here , , is the cosmological constant, and , i=1,2,3 are dimensionless parameters. The terms quadratic in the curvature are needed in order to renormalize the theory, as is the self-interaction term for . Note that the classical field provides a position-dependent mass for the field , and therefore the propagation of will be suppressed in regions where the mass is very high. In this sense models a “mirror”. For example, a thin mirror located at is described by the interaction

(4) |

and the perfect conductor limit corresponds to .

The classical field equations are

(5) |

(6) |

(7) |

where

(8) |

and

(9) |

The tensors come from the variation of the terms quadratic in the curvature contained in the gravitational action. The classical energy-momentum tensor for the field is

(10) |

while is the energy-momentum tensor for a free field with variable mass , that is,

(11) |

We consider now the semiclassical version of the theory, in which the field becomes a quantum field while and are treated classically. The Heisenberg equation for the quantum operator associated to is given by the classical Eq.(5). The evolution equations for the classical backgrounds are obtained by taking the mean value of the classical Eqs.(6) and (7):

(12) |

(13) |

We have therefore established a well defined model to study the divergences of the self-energy in Casimir calculations. The quantities and are formally divergent. The divergences must be absorbed into the bare constants of the theory. The renormalized version of contains information about the Casimir effect (force between different objects), as well as the self-energy. The gravitational effects produced by the vacuum energy can in principle be computed by considering it as a source of the semiclassical Einstein equations.

## Iii Renormalizability

The theory of interacting fields in curved spacetimes can be renormalized using a precise covariant procedure books qftcs (). In the present model it is necessary to introduce minor modifications to take into account that not only the metric but also one of the interacting fields is treated classically.

It will be particularly useful to adapt the renormalization method described in Ref.paz (). It is shown there that to analyze the renormalizability of theory in curved spaces at the level of the equations of motion, one can split the field as , where is the mean value of the field and is the quantum operator that describes the fluctuations around the mean value. To one loop order, satisfies a free field equation with a variable mass. The situation in our model is very similar, since the quantum field can be thought of a free field with variable mass. Therefore, we can follow closely Ref.paz ().

We define the renormalized quantities

(14) |

where and are constructed using the Schwinger DeWitt expansion up to fourth and second adiabatic order respectively adiab (). The divergences present in these quantities are to be absorbed into the bare constants of the theory.

The usual Schwinger DeWitt expansion for the propagator of a massive field reads SDW ()

(15) |

where is the Feynman propagator, is the van-Vleck determinant, is one half of the geodesic distance between and and is the number of spacetime dimensions, that acts as a regulator. The functions are defined by a set of recursive equations that follows from imposing the equation for the propagator.

The quantities to be subtracted to cancel the divergences of and are paz ()

(17) |

and

(18) |

We emphasize that we are not using a point-splitting regularization but a dimensional regularization. We will substitute and in the semiclassical equations by and , respectively, and absorb the infinities contained in and into the bare constants appearing in those equations.

### iii.1 Renormalization of the equation for the classical field

The coincidence limit of the two point function for a field with variable mass is given by paz ()

(19) |

where is an arbitrary constant with dimensions of mass. Inserting this expression into the semiclassical equation for and writing the bare constants in terms of renormalized ones

(20) |

we obtain

(21) |

Therefore, the divergences can be absorbed with the counterterms

(22) |

where and are finite contributions (they vanish in the minimal prescription scheme). Note that in the conformal case ( and ), only a counterterm for the self-coupling of is needed.

### iii.2 Renormalization of the semiclassical Einstein equation

Let us now consider the renormalization of the gravitational sector of the theory. As , the divergent part of is given by pazappendix ()

(23) |

This expression contains both geometric divergences and divergences dependent on the classical field (and its derivatives). The former should be absorbed into a redefinition of the gravitational constants appearing in the left hand side of the semiclassical Einstein equations, and the latter into a redefinition of the constants associated to the field, which appear in .

Inserting Eq.(23) into Eq.(12) one can show that the divergences dependent on can be absorbed using the same counterterms given in Eq.(22). This is a non trivial check of our calculations, and a necessary condition for the renormalizability of the theory. On the other hand, the geometric divergences can be absorbed into the gravitational constants, by choosing the following counterterns

(24a) | |||

(24b) | |||

(24c) |

(25a) | |||

(25b) |

This completes the proof of the renormalizability of the model, and is one of the main results of this paper. We have shown that, if the presence of a mirror is modeled by the interaction of the vacuum field with a classical background field, the divergences in the vacuum expectation value can be absorbed into the bare constants appearing in the semiclassical Einstein equations. Not only the gravitational constants are renormalized but also the bare constants associated to the Lagrangian of the classical background field. The renormalizability is valid for one or more mirrors of arbitrary shape, as long as they can be described by a smooth function . Note that the divergences in the energy-momentum tensor are considerably simpler in the conformal case and .

As a final remark, we stress that the quantities and will be finite for a sufficiently smooth background field and spacetime metric . We will illustrate this in the next section.

## Iv Explicit evaluations in the weak field approximation

In this section we present simple expressions for evaluating explicitly and for a given background potential , within a weak field approximation. For simplicity we will work in Minkowski space, and set so that the quantum field is massless. The sense in which we define the weak field approximation is that we require , or in other words . This can be achieved by having a weak coupling , or a weak and/or rapidly varying field . The calculation will stay at the lowest order in .

The procedure we follow is similar to the one used in stars (). It is based in solving the equation for the Feynman Green function,

(26) |

to first order in , and using the resulting to compute and in the coincidence limit. Expanding Eq.(26) by writing in a expansion and discarding the second-order term, we obtain

(27) |

This is solved explicitly using a momentum-space representation for :

(28) |

After switching variables to and , the -integrals can be done by passing to -dimensional space and employing standard dimensional regularization techniques, with the result in the coincidence limit being:

(29) |

Here is an arbitrary mass scale introduced in the regularization procedure. The divergence in the first term, which is purely local, can be absorbed in the coefficient of the action, as can be seen by comparison with the second line of Eq.(22). Hence in accordance with Eq.(17) we have

(30) |

For example, in the particular case in which the background potential depends only on the spatial coordinate , we shall have:

(31) |

This expression can be used to compute for the case of one or several parallel flat mirrors, which are “almost transparent” in the sense that the expression is valid to first order in the coupling .

We will prove now that according to this model is everywhere finite if the background potential is an integrable function. Under this assumption, it follows as a corollary of the Riemann-Lebesgue lemma that the Fourier transform of (which we note by ) falls off faster than . Hence is integrable, from which it follows that its Fourier transform is in turn well-defined and finite at all points. Since according to (31) it is precisely the Fourier transform of that gives the nonlocal part of , and the local part is the term involving , it follows that both the local and the nonlocal parts of are finite for all if is an integrable function.

We have therefore proved that the divergences in are removed if the background field modeling the mirrors is sufficiently well-behaved. We are not able to give a general proof using weaker assumptions than continuity. However, we will see that in simple concrete examples the assumption of continuity is sufficient to obtain finite results.

Let us first consider a discontinuous mirror of width , specified by the background potential:

(32) |

where is the unit step function. The parameter has length dimension and is introduced so that is normalized by . Introducing this expression into Eq.(31) we obtain as a result:

(33) |

where we omit the local contribution coming from the term. Notice that the result diverges logarithmically at the boundary of the mirror, where is discontinuous.

Consider next a continuous mirror given by

(34) |

wherein the discontinuities in Eq.(32) are resolved with linear interpolations of width . For this mirror the vacuum polarization is computed to be:

(35) |

where we have omitted again a local contribution proportional to . This result is everywhere finite; it can be compared to that of Eq.(33) in the plot in Figure 1.

We turn now to the study of to first order in . From the expression (11) it follows that (we drop the label from now on) is obtained from the coincidence limit of as:

(36) |

The coincidence limit of the last term is , and since is we may discard this term in our first-order approximation. The coincidence limit of the remaining term can be performed in an analogous way to that of , by passing to -dimensional space, using the momentum space representation (28) for the first-order , and using standard dimensional regularization techniques. The divergences appearing are proportional to , and can be absorbed into a redefinition of as we did before. After the renormalization is performed and we return to we are left with

(37) |

which is analogous to (30) for . For the case in which depends only on we get by analogy with (31):

(38) |

Comparing with our arguments showing finiteness for it is clear that the components of will not diverge as long as is integrable. To ensure this it is sufficient to require to be an integrable function, since then its Fourier transform falls off faster than . However, for the simple particular example of a mirror similar to (34) but with polynomials with higher degree of continuity interpolating between and , we have checked that continuity of is enough to render the result finite.

We stress at this point that the results obtained in this section have their gravitational counterparts: as shown in Ref.stars (), when one computes and for a free quantum field in a curved background, both quantities diverge at the points where the background is not sufficiently smooth. The examples described in Ref.stars () refer to the vacuum polarization around spherically symmetric objects, and divergences show up for sharp interfaces, where the matter density is discontinuous. The divergences disappear when the matter density and its first two derivatives are continuous across the interface. Similarly, in a cosmological context, divergences in the renormalized energy-momentum tensor are removed when one replaces an abrupt transition of the scale factor by a smoother transition in which the scale factor and its two first derivatives are continuous. The concrete example of the transition between the inflationary period and radiation domination is discussed in detail in Ref. Ford ().

High vacuum energy densities in the presence of boundaries may potentially produce large gravitational effects. In our model, the usual divergences are rendered finite when the sharp boundary is replaced by a sufficiently smooth background field. However, the quantum vacuum energy does attain high values if the background field varies over short distances. We can make a crude estimation of the maximum value of the energy density using dimensional analysis. From Eq.(IV), it is not difficult to see that if is varying over distances of order then we should expect the maximum value of to be of order , where is a representative value attained by . The quantity has units of (length), and determines the properties of the mirror. Therefore it will be identified with the square of the plasma frequency . Restoring and factors, the mass density associated to the vacuum fluctuations is of order . If we take to be in the range of the meters, representing a smearing of the conductor’s sharp boundary over atomic length scales, and to be in the range of Hz, a typical value of the plasma frequency, then we see that firstly, the first-order approximation assumed in this section is validated, and secondly, the peak values of the quantum energy densities at the boundary of the conductor are of order , too small to be gravitationally detectable.

## V Conclusions

In this paper we analyzed the coupling of the vacuum self-energy to gravity. We considered a model in which the presence of the bodies that disturb the modes of the quantum fields is described by the interaction with a background classical field. We have shown that the divergences in the energy-momentum tensor of the quantum fields are consistent with the semiclassical Einstein equations, that is, they can be absorbed into the bare constants of the theory. As expected, the divergences in the renormalized energy-momentum tensor noticed in previous works only appear when considering unphysical limits of perfect conductivity and/or sharp interfaces, that in our model would correspond to take a non-smooth background field. We have shown that modeling the mirrors in the Casimir effect with a background potential removes the divergences attached to the boundaries, as long as is a sufficiently smooth integrable function ( continuity is sufficient for and continuity for , though weaker assumptions are enough in simple examples). The proof is carried to first order in the coupling . These results are in tune with those found in Refs.Milton (); Bouas (), where finiteness of the results is shown with exact calculations for some particular examples of , and are analogous to those in Refs.Ford (); stars () for the vacuum polarization of free fields in curved backgrounds, when the metric and its first two derivatives are continuous. We expect that similar conditions for finiteness can be obtained beyond the perturbative approximation.

We have considered a toy model for a vacuum scalar field. We have presented some specific examples to illustrate the appearance of divergences for sharp interfaces. The examples suggest that the gravitational effects of the self-energy near mirrors are extremely small for realistic values of the relevant parameters. We expect similar results to be valid for more realistic models involving the electromagnetic field.

## Acknowledgements

This work was supported by UBA, CONICET and ANPCyT.

## References

- (1) P. W. Milonni, it The Quantum Vacuum, Academic Press, San Diego, 1994; M. Bordag, U. Mohideen, and V. M. Mostepanenko, Phys. Rep. 353, 1 (2001); K. A. Milton, The Casimir Effect: Physical Manifestations of the Zero-Point Energy (World Scientific, Singapore, 2001); S. Reynaud et al., C. R. Acad. Sci. Paris IV-2, 1287 (2001); K. A. Milton, J. Phys. A: Math. Gen. 37, R209 (2004); S.K. Lamoreaux, Rep. Prog. Phys. 68, 201 (2005); M. Bordag, G.L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Advances in the Casimir Effect, Oxford University Press, Oxford, 2009.
- (2) D. Deutsch, P. Candelas, Phys. Rev. D20, 3063 (1979).
- (3) N. Graham, R. L. Jaffe, V. Khemani, M. Quandt, M. Scandurra, H. Weigel, Nucl. Phys. B645, 49-84 (2002); ibidem Phys. Lett. B572, 196-201 (2003); ibidem Nucl. Phys. B677, 379-404 (2004).
- (4) K. Milton, arXiv:1107.4589.
- (5) J.D. Bouas, S.A. Fulling, F.D. Mera, K. Thapa, C.S. Trendafilova, and J. Wagner, arXiv:1106.1162.
- (6) K. A. Milton, in Lect. Notes Phys. 834 (2011). Edited by D. Dalvit, P. Milonni, D. Roberts, and F. da Rosa.
- (7) L. H. Ford, Phys. Rev. D35, 2955 (1987).
- (8) A. Satz, F. D. Mazzitelli, E. Alvarez, Phys. Rev. D71, 064001 (2005).
- (9) N. D. Birrell, P. C. W. Davies, “Quantum Fields In Curved Space,” Cambridge University Press (1982); S.A. Fulling, “Aspects of Quantum Field Theory in Curved Spacetime”, London Mathematical Society Student Texts (1989); L. Parker and D. Toms, “Quantum Field Theory in Curved Spacetime: Quantized Fields and Gravity”, Cambridge University Press (2009).
- (10) J. P. Paz, F. D. Mazzitelli, Phys. Rev. D37, 2170-2181 (1988).
- (11) The adiabatic order is usually defined as the number of derivatives of the metric tensor. In this case, and as discussed in Ref.paz (), the definition is generalized as follows: in terms containing derivatives of the classical field , the adiabatic order is the number of derivatives of plus the exponent of .
- (12) S. M. Christensen, Phys. Rev. D 14, 2490 (1976).
- (13) See Ref.paz (), Appendix B.