We give a comprehensive presentation of methods for calculating the Casimir force to arbitrary accuracy, for any number of objects, arbitrary shapes, susceptibility functions, and separations. The technique is applicable to objects immersed in media other than vacuum, nonzero temperatures, and spatial arrangements in which one object is enclosed in another. Our method combines each object's classical electromagnetic scattering amplitude with universal translation matrices, which convert between the bases used to calculate scattering for each object, but are otherwise independent of the details of the individual objects. The method is illustrated by re-deriving the Lifshitz formula for infinite half spaces, by demonstrating the Casimir-Polder to van der Waals cross-over, and by computing the Casimir interaction energy of two infinite, parallel, perfect metal cylinders either inside or outside one another. Furthermore, it is used to obtain new results, namely the Casimir energies of a sphere or a cylinder opposite a plate, all with finite permittivity and permeability, to leading order at large separation.

Abstract. We develop an exact method for computing the Casimir energy between arbitrary compact objects, both with boundary conditions for a scalar field and dielectrics or perfect conductors for the electromagnetic field. The energy is obtained as an interaction between multipoles, generated by quantum source or current fluctuations. The objects' shape and composition enter only through their scattering matrices. The result is exact when all multipoles are included, and converges rapidly. A low frequency expansion yields the energy as a series in the ratio of the objects' size to their separation. As examples, we obtain this series for two spheres with Robin boundary conditions for a scalar field and dielectric spheres for the electromagnetic field. The full interaction at all separations is obtained for spheres with Robin boundary conditions and for perfectly conducting spheres.

The Casimir problem is usually posed as the response of a fluctuating quantum field to externally imposed boundary conditions. In reality, however, no interaction is strong enough to enforce a boundary condition on all frequencies of a fluctuating field. We construct a more physical model of the situation by coupling the fluctuating field to a smooth background potential that implements the boundary condition in a certain limit. To study this problem, we develop general new methods to compute renormalized one-loop quantum energies and energy densities. We use analytic properties of scattering data to compute Green's functions in time-independent background fields at imaginary momenta. Our calculational method is particularly useful for numerical studies of singular limits because it avoids terms that oscillate or require cancellation of exponentially growing and decaying factors. To renormalize, we identify potentially divergent contributions to the Casimir energy with low orders in the Born series to the Green's function. We subtract these contributions and add back the corresponding Feynman diagrams, which we combine with counterterms fixed by imposing standard renormalization conditions on low-order Green's functions. The resulting Casimir energy and energy density are finite functionals for smooth background potentials. In general, however, the Casimir energy diverges in the boundary condition limit. This divergence is real and reflects the infinite energy needed to constrain a fluctuating field on all energy scales; renormalizable quantum field theories have no place for ad hoc surface counterterms. We apply our methods to simple examples to illustrate cases where these subtleties invalidate the conclusions of the boundary condition approach.

We have developed an exact, general method to compute Casimir interactions between a finite number of compact objects of arbitrary shape and separation. Here, we present details of the method for a scalar field to illustrate our approach in its most simple form; the generalization to electromagnetic fields is outlined in Ref. [1]. The interaction between the objects is attributed to quantum fluctuations of source distributions on their surfaces, which we decompose in terms of multipoles. A functional integral over the effective action of multipoles gives the resulting interaction. Each object's shape and boundary conditions enter the effective action only through its scattering matrix. Their relative positions enter through universal translation matrices that depend only on field type and spatial dimension. The distinction of our method from the pairwise summation of two-body potentials is elucidated in terms of the scattering processes between three objects. To illustrate the power of the technique, we consider Robin boundary conditions φ − λ∂nφ = 0, which interpolate between Dirichlet and Neumann cases as λ is varied. We obtain the interaction between two such spheres analytically in a large separation expansion, and numerically for all separations. The cases of unequal radii and unequal λ are studied. We find sign changes in the force as a function of separation in certain ranges of λ and see deviations from the proximity force approximation even at short separations, most notably for Neumann boundary conditions.