Gluing Formula for Casimir Energies

We provide a completely new perspective for the analysis of Casimir forces in very general piston configurations. To this end, in order to be self-contained, we prove a “gluing formula” well known in mathematics and relate it with Casimir forces in piston configurations. At the center of our description is the Dirichlet-to-Neumann operator, which encodes all the information about those forces. As an application, the results for previously considered piston configurations are reproduced in a streamlined fashion.


Introduction
Casimir energy and force computations are often plagued by divergences.In order to control these divergencies, different regularization techniques have been applied [1][2][3][4][5][6][7][8][9].However, these typically lead to different divergencies, which raises the question of their interpretation.In order to avoid this type of problem, configurations for which there are unambiguous finite answers have received considerable interest.Among these are the "piston configurations" introduced by Cavalcanti [10].
A piston configuration consists of two chambers M 1 and M 2 divided by a movable partition N into two compartments; see Figure 1.

Introduction
Casimir energy and force computations are often plagued by divergences.In order to control these divergencies, different regularization techniques have been applied [4,5,8,14,16,17,29,35,40].However, they typically lead to different divergencies which raises the question of their interpretation.In order to avoid this type of problem, configurations for which there are unambiguous finite answers have received considerable interest.Among them are the so-called piston configurations introduced by Cavalcanti [10].
A piston configuration consists of two chambers M 1 and M 2 divided by a movable partition N into two compartments; see Figure 1.An interesting question then is how the geometry and topology of the chambers influence the force on the piston, and how that force depends on the boundary condition imposed.Many configurations like flat Date: January 15, 2018.2000 Mathematics Subject Classification.Primary: 58J20; Secondary: 14F40.Key words and phrases.Casimir energies, Casimir forces, BFK-gluing formula, zeta-determinant, Dirichlet-to-Neumann operator.
The first author was supported by the Baylor University Summer Sabbatical Programme.The second author was supported by the National Research Foundation of Korea with the Grant number 2016R1D1A1B01008091.  Interesting questions then are how the geometry and topology of the chambers influence the force on the piston, and how that force depends on the boundary condition imposed.Many configurations, such as flat pistons at zero temperature [11][12][13][14][15][16][17][18][19][20] or finite temperatures [21][22][23][24], as well as curved pistons [25][26][27][28][29], have been analyzed on the basis of the spectrum of a Laplace-type operator associated with M 1 and M 2 .It is the aim of this article to introduce a completely new perspective on this type of analysis.
The starting point for developing this new perspective is the gluing formula on a compact oriented Riemannian manifold.Although this formula was proven much more generally in [30] (see also [31,32]), in our context, using it for Laplacians is all we need, and in Section 2, we provide an independent proof and an elementary example.In Section 3, we use the gluing formula to derive a gluing formula for Casimir energies.This can be used to reformulate the Casimir force on a piston in terms of the Dirichlet-to-Neumann map, which is done in Section 4. Examples are considered to indicate how streamlined this new perspective is.The conclusions summarize the most important results and the plan of where we go from here.

Proof of the BFK-Gluing Formula for Zeta-Determinants of Laplacians
In this section, we introduce the Burghelea-Friedlander-Kappeler gluing formula, which we call the BFK-gluing formula from now on.We give an independent proof for the case of relevance for our paper.
Let (M, g) be an (m + 1)-dimensional compact oriented Riemannian manifold.Furthermore, let N be a compact hypersurface of M, so that the closure of M − N has two components, M 1 and M 2 .The manifold M therefore consists of the manifolds M 1 and M 2 glued together; see Figure 1.It is a natural question to ask if there are relationships between spectral quantities on M 1 , M 2 , and the manifold M = M 1 ∪ N M 2 obtained when gluing them together.For a particular spectral quantity, the BFK-gluing formula sheds some light on this question.

Remark 1.
The assumption for the gluing formula to hold is that ∂N = ∅.In Figure 1, this means M is a two-dimensional surface and N is a circle.We come back to this assumption in the conclusions, providing an example in which the BFK-gluing formula holds despite that ∂N = ∅.
In order to formulate the BFK-gluing formula, let ∆ M denote the Laplacian acting on smooth functions, and denote by ∆ M 1 ,D and ∆ M 2 ,D the restriction of ∆ M to M 1 , M 2 with Dirichlet boundary conditions imposed on N. In the case that ∂M = ∅, we also impose Dirichlet boundary conditions on ∂M.Let α j , α 1,j and α 2,j be the eigenvalues of (0 ≤ λ ∈ R): We denote the associated zeta functions by ζ M (s), ζ M 1 ,D (s) and ζ M 2 ,D (s) and have the standard definitions: In case there is a zero mode on M, for example α 0 = 0, it is excluded in the summation.These zeta functions are meromorphic functions in the complex plane, and in particular they are analytic about s = 0) are well-defined quantities [33] that one defines the associated determinants through: ln Det P = −ζ P (0) These quantities make up the left-hand side (LHS) of the BFK-gluing formula.A relation to quantum field theory is established by choosing the manifolds M, M 1 , and M 2 to have the product structures M = M × S 1 , M 1 = M 1 × S 1 , and M 2 = M 2 × S 1 , where S 1 is a circle of perimeter β.In this case, these quantities are the partition sums for a finite-temperature quantum field theory at temperature 1/β of a non-interacting massive scalar field living on M, M 1 , and M 2 ; see Section 3. The remaining ingredient for the BFK-gluing formula is the Dirichlet-to-Neumann operator.Let ∂ N be a unit normal vector field to N such that it points outward on M 1 and inward on M 2 .The Dirichlet-to-Neumann operator R(λ) : We then define This specifies an elliptic pseudodifferential operator of order 1; the associated zeta function ζ DN (s) is regular at s = 0, and ζ DN (0) is again well defined [31,33].
The BFK-gluing formula relates these quantities through a polynomial: which is determined as an integral of some local density on N and which is completely determined by data on a collared neighborhood of N. In detail, we have or using zeta functions instead, If zero modes are present, slight modifications to Equation (3) occur [30], which will, however, be irrelevant for our application.
Given this is the basis of everything that follows, for a self-contained presentation, we include the proof of this statement for zeta-determinants of Laplacians on a compact or complete Riemannian manifold given in [30].
The underlying idea to prove Equation ( 3) is that by taking [m/2] + 1 derivatives on both sides, equality of both sides is obtained.For the formulation of the proof, we introduce some notation.As before, let (M, g) be a complete oriented (m + 1)-dimensional Riemannian manifold with a compact hypersurface N. We denote by M 0 the closure of M − N so that (M 0 , g 0 ) is a Riemannian manifold with boundary N + ∪ N − , where g 0 is a metric induced from g and N + = N − = N.We here note that M 0 may or may not be connected.Furthermore, as before, we denote by ∆ M a Laplacian acting on C ∞ (M), and by ∆ M 0 ,D , the extension of ∆ M to M 0 with the Dirichlet boundary condition on N + ∪ N − .By a Laplacian, we mean a symmetric non-negative differential operator of order 2 with the principal symbol σ L (∆ M )(x, ξ) = ξ 2 satisfying the property of the finite propagation speed (FPS).It is shown in [34] that ∆ M is essentially self-adjoint.Furthermore, we assume that ∆ M 0 ,D is a non-negative invertible operator.These assumptions show that for λ ∈ R + ∪ {0}, (∆ M + λ) −1 and (∆ M 0 ,D + λ) −1 are bounded operators defined on L 2 (M).In many cases, operators arising from geometry satisfy these assumptions [35][36][37].
It is known (p 327 of [31]; see also [38]) that e −t∆ M − e −t∆ M 0 ,D is a trace-class operator and that for where a j ∈ R. Because Tr e −t∆ M − e −t∆ M 0 ,D is a trace-class operator and ∆ M 0 ,D is an invertible operator, the dimension of the kernel of ∆ M is finite.Hence Theorem 2.2 of [39] tells that for some c > 0 as t → ∞, we have This shows that for λ ∈ R + , is exponentially decreasing for t → ∞.We note that This shows that ) is a trace-class operator and where If we repeat this argument, after ν = [ m 2 ] + 1 integrations, we have for some c > 0 that Similarly, it can be shown that for t → 0 + , and hence is integrable.
Taking the limit t → 0 + in Equation (10), this leads to the following result.
We define a relative zeta function by 2 and has an analytic continuation to the whole complex plane, having a regular value at s = 0. We define the relative zeta-determinant by For general facts about the relative zeta-determinant, we refer to [39].
Lemma 2. For λ > 0, we have the following equality: Proof.We note that for Re s > m 2 , Using Equations ( 11) and ( 12) together with Equation ( 14), we have We next analyze the Dirichlet-to-Neumann operator For this purpose, we need to define some auxiliary operators.We recall that the boundary of M 0 is N + ∪ N − .We choose a unit normal vector field ∂ N along N, which points outward on N + and inward on N − .We define For later use, we define the restriction maps γ 0 and γ 0 as follows: Finally, we define the Poisson operator P D (λ) : In fact, ψ is given by where f is an arbitrary extension of ( f , g) to M 0 .Because ∆ M 0 ,D + λ is an invertible operator, ψ exists uniquely.We then define the Poisson operator P D (λ) by Using the above operators, the Dirichlet-to-Neumann operator R(λ) defined in Equations ( 1) and ( 2) is then reformulated as follows.
where "•" means the composition of operators.We need the following two auxiliary lemmas.
Proof.From the definition of the Poisson operator, P D (λ) satisfies the following two equalities: Taking the derivative with respect to λ, from here we have are bounded operators and L 2 (M 0 ) can be identified with L 2 (M).Hence, ∆ M 0 ,D + λ −1 can be regarded as an operator acting on L 2 (M).The next lemma shows the relation between the two operators.
Proof.Because the operators on the LHS and right-hand side (RHS) are bounded, it is enough to check the equality on From Definition 1 and Lemmas 3 and 4, we then have the following equalities: From the above lemma, we have the following equalities: which shows that where we have used the fact that γ 0 • P D (λ) • δ diag = Id.Lemmas 2 and 4 then finally yield the following equality: where in the third step, cyclicity of the trace has been used.Summarizing the above argument, we have shown the BFK-gluing formula in the context of Laplace-type operators.
Theorem 1.There exists a real polynomial P(λ) of degree less than or equal to [ m 2 ] such that This is equivalent to Equations ( 3) and ( 4).The fact that p j is given as an integral of some density on N follows from the λ → ∞ behavior of ln Det ∆ M + λ, ∆ M 0 ,D + λ and ln Det R(λ) [30].

Example 1.
As an illustration, we consider a simple one-dimensional example.Let 0 < c < 1 and I c = [0, c].The setting for the gluing formula is I 1 = I c ∪ {c} [c, 1].We consider the following eigenvalue problem: The eigenfunctions and eigenvalues then are and the associated zeta function is given by the zeta function ζ R (s) of Riemann.In detail, and thus Similarly the intervals [c, 1] and I 1 can be treated and the relevant combination on the LHS of Theorem 1 then is On the other side, the Dirichlet-to-Neumann map involves the solutions to the following problems: In terms of these, one defines x=c Explicitly, it is easy to see that For this example, we therefore verify that Theorem 1 is satisfied with P(0) = − ln 2

Gluing Formula for Casimir Energies
We next relate the gluing formula to a finite-temperature quantum field theory of a non-self-interacting massive scalar field in a piston geometry.Relevant one-particle energy spectra then follow from eigenvalues of suitable Laplace-type operators.As indicated earlier, to relate this quantum field theory to the considerations of Section 2, let M i = M i × S 1 , where S 1 denotes a circle of perimeter β, and M i denotes a smooth compact Riemannian manifold with smooth boundary N .This manifold represents the spatial dimensions of the space-time where the quantum field lives.As far as Figure 1 goes, what one sees there is really M 1 and M 2 , although for the application to the gluing formula, M i = M i × S 1 is needed.
Let τ ∈ [0, β] parameterize S 1 and impose periodic boundary conditions along that circle.Furthermore, let y be coordinates on M i .For the analysis of Casimir energies, the relevant eigenvalue problems then are, with λ being the mass of the quantum field, More explicitly, because of the product structure, where E 2 i,j are the eigenvalues of the "spatial part" of the Laplacian, namely we have The associated zeta functions, encode the finite-temperature Casimir energy via Leaving a discussion of finite ambiguities aside, when discussing forces, there will be none; it is an interesting question to ask how the Casimir energy behaves when gluing together two manifolds with identical boundaries.More precisely, the question is what the relation is between E Cas , the Casimir energy on M = M × S 1 , and E (1) Cas .An answer is provided by the gluing formula; namely from Equation ( 4), one has This is the gluing formula for Casimir energies.The change of the Casimir energy of the manifolds M 1 and M 2 when glued together along their boundary N × S 1 is described by the zeta function associated with the Dirichlet-to-Neumann map and a (partially known) expression in terms of geometric tensors of their boundary N [30,40].

Dirichlet-to-Neumann Map and Casimir Forces on Pistons
The gluing formula for Casimir energies allows immediately a reformulation of the Casimir force in piston configurations.Let a denote the position of N within the "cylindrical part" of the configuration in Figure 1.Then the force on the piston is Cas + E Cas − E We note that P(λ) is completely independent of the location of N within the cylindrical part, and the force is encoded completely by the Dirichlet-to-Neumann map.Thus the force is fully determined by the "harmonic functions" of M 1 and M 2 , that is, by solutions of (∆ M j + λ)φ j = 0.
In order to observe an immediate advantage of this formulation, we recall that the Dirichlet-to-Neumann map is a map from functions on N to functions on N.That is, the force is determined from a map between an m-dimensional manifold instead of the (m + 1)-dimensional manifold we started with.Thus this formulation achieves a dimensional reduction without any effort.
We now use this new formulation to obtain the Casimir force on pistons for configurations that have been dealt with before.It will become apparent that the dimensional reduction for explicit examples technically means that one summation is explicitly performed without ever worrying about it.
We consider the configuration analyzed in, for example, [18,23].Each chamber is assumed to have the structure Typically, N = C × K, where C is the (visible) cross-section of the cylinder and K denotes a manifold representing additional Kaluza-Klein dimensions.Finally, I 1 = [0, a] and I 2 = [a, L] are the intervals giving the extension of each cylindrical part of the chambers.Assuming Dirichlet boundary conditions at x = 0 and x = L, we next formulate the relevant boundary value problems to obtain the Dirichlet-to-Neumann operator for this setting.In the left chamber we need to consider By choosing ϕ(τ, y) = e 2πinτ/β ϕ j (y), the problem of Equation ( 31) turns into the following easily solvable problem: where µ 2 j are eigenvalues of the following eigenvalue problem: The general solution of the ordinary differential equation in Equation ( 32) is where ρ n,j = 2πn Imposing the boundary condition at x = 0 shows The boundary conditions at x = a give ψ 1,j (a) = 1 = A e ρ n,j a − e −ρ n,j a =⇒ A = 1 e ρ n,j a − e −ρ n,j a and thus ψ 1,j (x) = e ρ n,j x − e −ρ n,j x e ρ n,j a − e −ρ n,j a Similarly, for the right chamber we need to solve Along the same lines, one obtains e ρ n,j (L−a) − e −ρ n,j (L−a) For the Dirichlet-to-Neumann map, we need the first derivative at x = a: and thus These are the eigenvalues of the Dirichlet-to-Neumann map with eigenfunctions e 2πinτ/β ϕ j (y)!It is clear that eigenvalues of the Dirichlet-to-Neumann map are fairly complicated functions of the eigenvalues of the "underlying" Laplacian on N.
Equation (35) shows that the zeta function associated with the Dirichlet-to-Neumann map is cosh(ρ n,j a) sinh(ρ n,j a) + cosh(ρ n,j (L − a)) sinh(ρ n,j (L − a)) −s (36) The analysis of this zeta function is simplified by observing that cosh(ρ n,j a) sinh(ρ n,j a) + cosh(ρ n,j (L − a)) sinh(ρ n,j (L − a)) = sinh(ρ n,j L) sinh(ρ n,j a) sinh(ρ n,j (L − a)) which allows us to write ζ DN (s) as From the large eigenvalue behavior of ρ n,j (Equation ( 34)), this series representation is clearly convergent for s > m.This form is very suitable for finding ζ DN (0), and we obtain where All series in Equation (38) are clearly absolutely convergent.
For the Casimir force on the piston, the only relevant pieces are a-dependent pieces, and we have where the β-differentiation has been performed.We compare this result with known answers.First, we consider the zero temperature limit, that is, β → ∞.In this limit, the n-summation turns into an integral, and we have Assuming for the moment that µ 2 j > 0, the integral can be performed using [41], 8.432.9: such that the Casimir force on the piston reads which is in exact agreement with Equation (4.3) in [18].The result clearly shows that the piston is pulled towards the closer wall.Also at a finite temperature, the result of Equation ( 39) is equivalent to known results presented in [18,23], although these appear to be different.The reason is that typically in finite-temperature computations, the Matsubara sum involving the n-summation above is manipulated.Here instead, the sum from the intervals I 1 and I 2 is automatically manipulated (performed), and thus final answers appear to be different, as effectively resummations for different summations have been performed.
The case for which µ j = 0 is possible, for example, µ 0 = 0 with a degeneracy g 0 ; for n = 0, this leads to zero modes on N. The contribution to ζ DN (0) can be computed along the lines of Example 1; see Equation (1).This is independent of β and does not enter the Casimir energy or force.Thus Equation (39) remains valid with the summation starting at j = 0, repeated according to its multiplicity.When considering the zero temperature limit, the required integral this time is of the k-sum, we first rewrite the k-sum and apply a Mellin transform to obtain Thus a gluing formula may exist also in the more general context of ∂N = ∅.
In order to more explicitly see the polynomial structure of P(λ), we note the following (for small enough λ): where ζ ∆ N is the zeta function associated with the Laplacian on N. From here, one finds where for the case ∂N = ∅, exact agreement with [31] is found (once possible zero modes on N are taken into account).We note that because of boundary contributions in heat kernel coefficients, values and residues of ζ ∆ N do not vanish generically, and thus the polynomial will not vanish for m odd, as it does for ∂N = ∅.
In addition to the described generalization to ∂N = ∅, in order to be applicable to the electromagnetic field, we plan to generalize the BFK-gluing formula to different boundary conditions.Thus, instead of imposing Dirichlet boundary conditions on the piston, Robin boundary conditions will be allowed, and the corresponding relevant Robin-to-Robin map needs to be found.
It is hoped that using the vast amount of literature on the Dirichlet-to-Neumann map, (see, e.g., [43][44][45][46][47][48][49]), new insights can be gained on what the decisive topological and geometrical factors are that lead to attractive or repulsive Casimir forces.

2 + 2 [ln 2 − 1 ln 1 −
respect to s at s = 0 can easily be taken, and one findsζ a (0) = −a FP ζ N − 1 2a Res ζ N − 1 e −2aρ jwhere FP and Res denotes the finite part and the residue, respectively, of the zeta function at the argument indicated.Adding up the relevant combination, the BFK-gluing formula is verified withP(λ) = − ln 2 ζ N (0)