the interaction between neutral objects is dominated by fluctuation forces due to the coordinated behavior of fluctuating charges or collective modes inside the objects . at zero temperature or at sufficiently small distances ,
the interactions result from quantum fluctuations .
they are important on atomic scales as well as between macroscopic bodies . a prominent example for the latter is the casimir force between parallel metallic plates due to current or electromagnetic field fluctuations @xcite .
during the last decade experimental verifications of this effect have been performed with increasing precision .
these high precision measurements were enabled by the use of curved surfaces instead of parallel planar mirrors in order to avoid the problem of parallelism .
the most commonly employed geometry is a sphere - plate setup that was used in the first high - precision tests of the casimir effect @xcite .
this geometry has been successfully used ever since in most of the recent experimental studies of casimir forces between metallic surfaces @xcite . in order to keep the deviations from two parallel plates sufficiently small , spheres with a radius much larger than the surface distance
have been used .
the effect of curvature has been accounted for by the `` proximity force approximation '' ( pfa ) @xcite .
this scheme is assumed to describe the interaction for sufficiently small ratios of radius of curvature to distance .
however , this an uncontrolled assumption since pfa becomes exact only for infinitesimal separations , and corrections to pfa are generally unknown .    at the other extreme , the interaction between a planar surface and objects that are either very small or at asymptotically large distance is governed by the casimir - polder potential that was derived for the case of an atom and a perfectly conducting plane @xcite .
this limit has been probed experimentally with high precision for a bose - einstein condensate that was trapped close to a planar surface @xcite .
recently , there have been first attempts to go beyond the two extreme limits of asymptotically large and small separations by measuring the casimir force between a sphere and a plane over a larger range of ratios of sphere radius to distance @xcite .
so far , no theoretical prediction is available that can describe the electromagnetic casimir interaction between a compact object and a planar surface at all distances , including the important sphere - plate geometry . until recently , progress in understanding the geometry dependence of fluctuation forces was hampered by the lack of practical methods that are applicable at all separations .
conceptually , the effect of geometry and shape is difficult to study due to the non - additivity of fluctuation forces .
explicit consequences of this non - additivity and also non - monotonic changes in the force have been recently predicted for a pair of cylinders next to planar walls @xcite .
this behavior has been interpreted in terms of collective charge fluctuations inside the bodies @xcite .
for some decades , there has been considerable interest in the theory of casimir forces between objects with curved surfaces .
two types of approaches have been pursued .
attempts to compute the force explicitly in particular geometries and efforts to develop a general framework which yields the interaction in terms of characteristics of the objects like polarizability or curvature . within
the second type of approach , balian and duplantier studied the electromagnetic casimir interaction between compact perfect metals in terms of a multiple reflection expansion and derived also explicit results to leading order at asymptotically large separations @xcite . for parallel and partially transmitting plates
a connection to scattering theory has been established which yields the casimir interaction of the plates as a determinant of a diagonal matrix of reflection amplitudes @xcite . for non - planar ,
deformed plates , a general representation of the casimir energy as a functional determinant of a matrix that describes reflections at the surfaces and free propagation between them has been developed in ref .
later on , an equivalent representation has been applied to perturbative computations in the case of rough and corrugated plates with finite conductivity @xcite .
functional determinant formulas have been used also for open geometries that do not fall into the class of parallel plates with deformations . for the electromagnetic casimir interaction between a planar plate and infinitely long cylinders , a partial wave expansion of the functional determinant
has been developed @xcite .
the same results have been reobtained and used to compute corrections to the pfa for the cylinder - plate geometry in ref .
kenneth and klich identified the inverted green s function in the functional determinant as a t - matrix and derived a formal result for the casimir interaction in terms of this matrix in the case of scalar fields in a medium with a space and frequency dependent speed of light @xcite .
recently we described a new method based on a multipole expansion of fluctuating charges that makes possible accurate and efficient calculations of casimir forces and torques between any number of compact objects @xcite .
the method applies to electromagnetic fields and dielectrics as well as perfect conductors .
it also applies to other fields , such as scalar and dirac , and to any boundary conditions . in this approach ,
the casimir energy is given in terms of the fluctuating field s scattering amplitudes from the individual objects , which encode the effects of the shape and boundary conditions .
an equivalent partial wave expansion has been applied to a scalar model for dielectrics @xcite .
casimir interactions due to scalar field fluctuations serve as a simplified model for the full electromagnetic interaction .
this model is usually easier to analyze and provides an important tool in developing conceptually new approaches and in estimating geometry dependencies . for the sphere - plate geometry with dirichlet boundary conditions bulgac _
_ obtained the casimir interaction over a wide range of separations from a modified krein trace formula @xcite .
most notably in the context of scalar fields , a versatile numerical world - line algorithm based on monte carlo methods has been developed and applied to a number of interesting geometries , including the here studied sphere - plate interaction @xcite .    in this work ,
we extend our approach developed in refs .
@xcite to describe the interaction of compact objects in the presence of a plane mirror .
our general result holds for scalar and electromagnetic fields . in the electromagnetic case the mirror
is assumed to be perfectly conducting but the compact objects can be dielectrics or perfect conductors .
the derivation of the general result for the casimir interaction , see eqs .   and ,
is given in section [ sec : derivatiobn ] by combining a functional integral approach and the method of images . in section [ sec
: plate - sphere ] we apply our approach to compute the interaction of a sphere with a plane mirror over a wide range of separations for the scalar dirichlet and neumann problem and for the electromagnetic field that is most relevant to applications .
we provide a large distance expansion of the interaction , generalizing the casimir - polder potential to include higher order multipole polarizabilities . for small separations , we compute the two leading correction terms to the pfa .
in this section we first review the functional integral formulation of casimir interactions between compact objects for a scalar field and the electromagnetic field .
the fluctuating field is integrated out in order to obtain an effective action for the fluctuating sources on the objects .
then we show that the interaction of compact objects with an infinite plane mirror can be described by the equivalent problem of the interaction of compact objects with their mirror images in otherwise empty space , i.e. , without the plane mirror .
this equivalence holds for dirichlet or neumann boundary conditions ( at the mirror and the objects ) for a scalar field and for a perfectly conducting mirror plane and arbitrary dielectric objects in the case of electromagnetic fluctuations .
we first consider a real quantum field @xmath0 in order to introduce the method of images in the path integral formulation of casimir interactions .
we assume that the space is divided into two half spaces by a mirror plane at @xmath1 . @xmath2
fixed closed surfaces @xmath3 , @xmath4 , are located in the right half space ( @xmath5 ) .
these surfaces are regarded as the boundaries of objects on which either dirichlet or neumann boundary conditions are imposed .
the action of the unconstrained field in minkowski space is @xmath6 = \frac{1}{2}\int dt \int _
> d{{\bf x}}\left\ { \frac{1}{c^2 } ( \partial_t \phi)^2   - ( \nabla\phi)^2\right\ } \ , .\ ] ] here the @xmath7-integration runs over the right half space , indicated by @xmath8 .
the free energy @xmath9 of the constraint field @xmath0 at inverse temperature @xmath10 is represented by the euclidean functional integral @xmath11_{{\mathcal c}}\exp\left ( -s_e[\phi]/\hbar \right ) \
, , \ ] ] where the euclidean action @xmath12 follows from eq .   after the wick rotation @xmath13 , @xmath14 with @xmath15 .
the functional integral is over all fields that are periodic in the @xmath16-interval from @xmath17 to @xmath18 and that obey the boundary conditions at the surfaces , indicated by the subscript @xmath19 .
the surfaces are fixed and the boundary conditions are time independent
. hence each fourier component of the field with respect to @xmath20 obeys the constraints at the surfaces separately . to make use of this property ,
we expand @xmath21 as @xmath22 with matsubara frequencies @xmath23 and @xmath24 . the functional integral splits into independent functional integrals over the @xmath25 and the logarithm of @xmath26 can be written as the sum @xmath27_{{\mathcal c}}\exp \left[-\frac{\lambda}{2\hbar c } \int _
x}}\left\ { \kappa_n^2 |\phi_n|^2 + |\nabla \phi_n|^2 \right\}\right ] \nonumber\ , .\end{aligned}\ ] ] in the following we are interested in the limit of zero temperature . then @xmath28 , and the sum over @xmath29 can be replaced by the integral @xmath30 and @xmath25 is replaced by @xmath31 .
combining positive and negative @xmath32 , we get @xmath33 with    @xmath34_{{\mathcal c}}[{{\mathcal d}}\phi^*({{\bf
x}},\kappa)]_{{\mathcal c}}\exp \left[-\frac{\lambda}{\hbar c } \int _
x}}\left\ { \kappa^2 |\phi|^2 + |\nabla \phi|^2 \right\}\right ] \ , .\ ] ]    the ground state energy is obtained from @xmath26 as @xmath35 the casimir energy is the difference between the ground state energy of the surfaces at their actual distance and the sum of the ground state energies of the separate surfaces which is obtained by removing the surfaces to infinite separation , @xmath36 where @xmath37 is the partition function for the surfaces with infinite separation . in the following
we suppress the label @xmath32 .
next , the constraints at the objects are implemented by functional @xmath38-functions @xcite . for dirichlet boundary conditions , @xmath39 , on the surfaces @xmath3 ,
the constraint functional integral can be expressed in terms of an unconstrained integral by using @xmath40_{{\mathcal c}}[{{\mathcal d}}\phi^*]_{{\mathcal c}}= \\ & & \!\!\!\!\!\int \!\!{{\mathcal d}}\phi { { \mathcal d}}\phi^ * \prod_{\alpha=1}^n \int \!\!{{\mathcal d}}\varrho_\alpha { { \mathcal d}}\varrho^*_\alpha \exp\left [ i\!\!\int_{\sigma_\alpha}\!\!\!\!\ ! d{{\bf x}}\,\left\{\varrho_\alpha^*({{\bf x } } ) \phi({{\bf x } } ) + { { \mbox{c.c.\,}}}\right\ } \right ]   , \nonumber\end{aligned}\ ] ] where the @xmath38-functions at each position of the surfaces have been written as an integral over a source field @xmath41 that is non - zero on the surfaces @xmath3 only .
when we use this representation of the constraints in eq .  , the now unconstrained functional integral over @xmath42 is gaussian and yields    @xmath43 \nonumber\\ & \equiv & { \mathcal z}_0 \prod_{\alpha=1}^n\int { { \mathcal d}}\varrho_\alpha{{\mathcal d}}\varrho^*_\alpha e^{-(\hbar c/\lambda ) \tilde s[\varrho]}\ , , \end{aligned}\ ] ]    where @xmath44 is the partition function of the field in the right half space without the surfaces @xmath3 .
for neumann boundary conditions at the surfaces @xmath3 the field @xmath42 in the exponential of eq .
is replaced by @xmath45 and @xmath46 in eq .
is replaced by @xmath47 , where @xmath48 is the normal derivative pointing out of the objects . here
@xmath46 is the free green s function in the right half space which is given by @xmath49 where @xmath50 is the mirror image of @xmath51 and the @xmath52(@xmath53 ) sign applies to dirichlet ( neumann ) boundary conditions at the mirror plane at @xmath1 .
the green s function of free unbounded space is given by @xmath54 using eq .
the action @xmath55 $ ] defined by eq .
can be expressed in terms of the original sources @xmath41 and their mirror sources @xmath56 , @xmath57 = \frac{1}{2 } \sum_{\alpha\beta } \left\ {   \int_{\sigma_\alpha } d{{\bf x}}\int_{\sigma_\beta } d{{\bf x } } ' \varrho^*_\alpha({{\bf x } } ) g_0({{\bf x}},{{\bf x}}',\kappa ) \varrho_\beta({{\bf x } } ' ) \right .
\nonumber\\ & + & \!\!\ ! \left .
\int_{\sigma_\alpha } d{{\bf x}}\int_{\sigma^r_\beta } d{{\bf x } } ' \varrho^*_\alpha({{\bf x } } ) g_0({{\bf x}},{{\bf x}}',\kappa ) \varrho^r_\beta({{\bf x } } ' ) + { { \mbox{c.c.\,}}}\right\ } \end{aligned}\ ] ] for dirichlet boundary conditions at the surfaces @xmath3 and with @xmath58 replaced by @xmath59 for neumann boundary conditions at the surfaces @xmath3 .
here we have introduced the mirror sources @xmath60 where the @xmath52 ( @xmath53 ) sign applies to a dirichlet ( neumann ) mirror .
the mirror sources are located on the mirror surfaces @xmath61 that are obtained from the @xmath3 by @xmath62 for all surface positions , see fig .
[ fig : objects ] .
the first term of the action of eq .   describes the interaction of the surface sources in the absence of the mirror plane .
the second term couples each surface source to all mirror sources .
since the mirror problem is now described by an action in _ free _ space with sources and mirror sources , we can apply the concepts of the previously developed approach for casimir interactions between compact objects in unbounded space @xcite .
below , we provide an explicit derivation for the case of dirichlet boundary conditions at the surfaces @xmath3 but we shall also indicate how the derivation has to be modified for neumann boundary conditions .
( objects @xmath63 ) in the right half space ( @xmath5 ) and mirror surfaces @xmath61 ( objects @xmath64 ) .
@xmath65 , @xmath66 are local coordinate vectors measured from the object s origins , @xmath67 is the center - to - mirror distance of surface @xmath68 ( object @xmath69 ) . ]
the action of eq .
is composed of two qualitatively different terms that we will now consider separately .
firstly , there are terms that couple sources on different surfaces ( @xmath70 ) where we use the term `` surface '' in the following for the original and the mirror surfaces .
we shall call these terms _ off - diagonal_. as _ diagonal _ terms we shall denote those which couple sources on the same surface ( @xmath71 possible ) . both type of terms can be expressed in terms of the multipole moments of the sources .
_ off - diagonal terms _
 these terms couple sources on different objects , @xmath72 for @xmath73 and the original sources to the mirror sources , @xmath74 for _ all _ @xmath75 , @xmath10 .
here we have introduced local coordinates @xmath65 that are measured relative to an arbitrarily chosen origin that is located inside the surface @xmath75 , see fig .
[ fig : objects ] .
we have also defined the fields @xmath76 which are the classical fields generated by the sources .
for neumann boundary conditions on the surfaces @xmath3 , @xmath77 and @xmath78 in eqs .   and have to be replaced by @xmath79 and @xmath80 , respectively . also , @xmath58 in eq .   has to be replaced by @xmath81 .
since we can assume that every position on @xmath3 is outside a sphere enclosing @xmath68 or @xmath82 , we can use the partial wave expansion of @xmath58 for @xmath83 when we consider eq .
with coordinates @xmath84 relative to the origin inside surface @xmath68 , the field that is generated by the source @xmath85 can be written as @xmath86 where we have defined the multipole moments of the source @xmath85 as @xmath87 the field of the mirror sources can be expressed in the same form , @xmath88 where @xmath89 denotes the local coordinates that are measured relative to the origin inside surface @xmath90 .
the multipole moments of the mirror source are given by @xmath91 due to eq .   and
for neumann boundary conditions at the surfaces @xmath3 the above expressions remain valid if the multipole moments @xmath93 are replaced by neumann multipoles which have the form of eq .   but with @xmath94 replaced by @xmath95 $ ] .    in order to express the action of eq .   in terms of multipole moments
, we have to write the field generated by the surface source @xmath85 as a function of the local coordinate @xmath65 that is regular at the origin inside @xmath3 .
this can be done using translation matrices @xmath96 which relate outgoing ( @xmath97 ) and regular ( @xmath98 ) spherical bessel functions , @xmath99 the matrix elements @xmath100 of @xmath96 depend on the vector @xmath101 from the origin inside @xmath3 to the origin inside @xmath68 , see fig .
[ fig : objects ] , and are given by @xcite    @xmath102    where we have assumed that the cartesian coordinate frames associated with the two origins have identical orientation , i.e. , they are related by a translation .
the summation over @xmath103 involves only a finite number of terms since the 3-@xmath104 symbols vanish for @xmath105 and @xmath106 .    using the translation formula of eq .  , the field generated by the source on surface @xmath68 , given by eq .
, can be written as function of the coordinate @xmath65 as @xmath107 where we have introduced @xmath108 . similarly we obtain for the field generated by the mirror sources of eq .
, now expressed as function of @xmath65 , @xmath109 where we defined @xmath110 and used eq .  .
notice that the latter formula applies also to the case @xmath111 which describes the translation between the surface @xmath3 and its mirror image so that the argument of the translation matrix becomes @xmath112 where @xmath113 is the normal distance between the origin of surface @xmath3 and the mirror plane .
for this case of translations along the @xmath114-direction the translation matrix simplifies to @xcite    @xmath115    when we substitute the result for the fields of eqs .  , into eqs .
, , we obtain the action in terms of the original multipole moments , @xmath116 which applies to dirichlet as well as neumann boundary conditions at the surfaces @xmath3 .
to simplify notation , we have defined the modified translation matrix elements @xmath117,\ , ( \alpha\neq\beta ) \nonumber\\\end{aligned}\ ] ] here we have used the definition of eq .   and @xmath118^*=(-1)^lj_l(i\kappa r)$ ] .
in addition , we applied the symmetry relations @xmath119 and @xmath120 which follow from eqs .
, , symmetry properties of 3-j symbols and spherical harmonics , and @xmath121 where @xmath122 is the ( real valued ) modified bessel function of second kind .
notice that the actions of eqs .  , couple sources on different surfaces and hide the dependence on the particular boundary conditions and shape of the surfaces .
_ diagonal terms _
 these are the self - action terms @xmath123 in eq .  , which we have expressed here in terms of the classical field defined in eq .
for neumann boundary conditions on the surface @xmath3 we have to replace again @xmath124 by @xmath125 in eq .  .
here and in the following we only use the coordinate system associated with the origin inside @xmath3 , and hence drop the label @xmath75 on the coordinates . the classical field generated by the source @xmath41 on surface @xmath3 obeys the helmholtz equation @xmath126 for positions @xmath7 that are located on the surface @xmath3 the _ total _ field @xmath42 generated by all sources must obey the same helmholtz equation .
the part of the total field that is generated by sources other than @xmath41 can be regarded as incident field @xmath127 at the surface @xmath3 , which obeys in region around @xmath3 , that is free of sources other than @xmath41 , the homogeneous helmholtz equation .
hence the total field can be written as @xmath128 for all @xmath7 located on @xmath3 . for neumann boundary conditions at @xmath3 the green s function in eq .
is replaced again by @xmath129 .
we would like to evaluate the action of eq .   in terms of multipole moments .
hence , we must consider field configurations with a fixed source on surface @xmath3 that is characterized by its multipole moments .
this implies that we have to find the incident field @xmath127 that induces a prescribed set of multipoles @xmath130 on @xmath3 .
the multipole moments can be identified as the amplitudes of the scattered field which is given by @xmath124 with @xmath7 located outside of the surface @xmath3 . using the partial wave expansion of eq .
, we get @xmath131 from scattering theory we know that the amplitudes of the scattered field are related to the amplitudes of the regular incident field by the transition matrix @xmath132 which is related to the scattering matrix @xmath133 by @xmath134 . if we expand the incident field as @xmath135 the amplitudes of the scattered field are given by @xmath136 where the @xmath137 denote the matrix elements of the @xmath132-matrix of the surface @xmath3 .
hence , the amplitudes of the incident field have to be given by @xmath138^{-1}_{lml'm ' } q_{\alpha , l'm ' } \ , .\ ] ] for dirichlet boundary conditions , the total field @xmath42 of eq .   has to vanish on @xmath3 so that @xmath139 on the surface .
hence , using eqs .   and , the action of eq .
can be expressed in terms of the multipole moments of @xmath41 , @xmath140^{-1}_{lml'm ' } \,q_{\alpha , l'm'}\ ] ] with the modified @xmath141-matrix defined by @xmath142 here we have used the definition of the multipoles in eq .   to integrate over the surface and applied the relation @xmath143 .
for neumann boundary conditions on @xmath3 , the normal derivative of the total field of eq .   has to vanish on @xmath3 so that @xmath144 on the surface .
when we use the definition of the multipole moments for neumann boundary conditions , we obtain again eq .   but with the matrix @xmath145 for neumann boundary conditions .
now we can combine all results for off - diagonal and diagonal terms and express the total action of eq .   in terms of the original multipole moments .
since we have @xmath146 we obtain from eqs .  ,
and the total action of the multipole moments @xmath147 = \sum_{\alpha\beta } q^*_\alpha { \mathbb m}^{\alpha\beta } q_\beta \ , , \ ] ] where we have suppressed the sum over the indices @xmath148 , @xmath149 , @xmath150 , @xmath151 and defined the matrix @xmath152^{-1 }   \delta_{\alpha\beta } - \tilde { \mathbb u}^{\alpha\beta } ( 1-\delta_{\alpha\beta } ) \pm \tilde { \mathbb u}^{r,\alpha\beta } \right\}\ ] ] for dirichlet ( @xmath53 ) or neumann ( @xmath52 ) boundary conditions at the mirror plane . the partition function of eq .
is then obtained by integrating over all multipole moments , @xmath153 \ , .\ ] ] the gaussian integral over the multipoles is proportional to the inverse determinant of @xmath154 . finally , we substitute into eq .   to obtain the casimir energy , @xmath155 where the determinant is taken with respect to the partial wave indices @xmath148 , @xmath149 and the surface indices @xmath75 , @xmath10 .
the matrix @xmath156 is the result of moving the surfaces to infinite separation , where the translation matrices @xmath157 , @xmath158 vanish so that @xmath159^{-1 }   \delta_{\alpha\beta } \ , .\ ] ] in the special case of one compact surface in front of the mirror plane eq .
simplifies to @xmath160 for dirichlet ( + ) or neumann ( - ) boundary conditions at the mirror plane .
this expression applies to dirichlet , neumann and even more general boundary conditions at the compact surface which enter only via the matrix @xmath161 .
notice that this result depends on the original matrix elements ( without tilde ) since the phase factors of eqs .   and
drop out when taking the matrix product of @xmath162 and @xmath163 .
this general result shows that the casimir interaction between a mirror and an object with arbitrary shape and boundary condition can be obtained from the transition matrix @xmath161 of the object and the translation matrix @xmath164 that describes the ( classical ) interaction between the induced source and its mirror image .
the derivation of the casimir energy for a scalar field can be extended to electromagnetic field fluctuations in the presence of dielectric objects @xcite .
the result will have the form of an effective action for electric and magnetic multipoles of the current densities @xmath165 inside the objects .
we consider again @xmath2 objects that are located in the right half space that is bounded by a perfectly conducting plane at @xmath1 . at his plane
the tangential electric field and the normal magnetic field vanish , @xmath166 , @xmath167 .
the material objects are assumed to be dielectrics that are characterized by a frequency dependent dielectric function @xmath168 and permeability function @xmath169 .
the partition function can be factorized again into a product of factors @xmath170 at a fixed wick rotated frequency @xmath32 .
hence , in the following we consider all expressions at fixed @xmath32 and suppress the label @xmath32 .
the euclidean action for the electromagnetic field in the presence of macroscopic media without external sources can be expressed as @xmath171\ ] ] in terms of the macroscopic fields @xmath172 , @xmath173 . here
the _ energy _ density of the field is integrated since under a wick rotation to imaginary time the lagrangian in real time is generally transformed to the hamiltonian in imaginary time . in this description
the induced ( bound ) currents inside the material objects have been absorbed into the definition of the macroscopic fields .
the partition function for this action is given by a functional integral over the vector potential @xmath174 and the scalar potential @xmath175 ( after introducing a faddeev - popov gauge fixing term ) where the fields are expressed in terms of the potentials as @xmath176    an alternative description in terms of the fields @xmath177 , @xmath178 only is obtained if the bound charges ( @xmath179 ) and currents ( @xmath180 ) density inside the objects are not substituted by the macroscopic fields but considered explicitly . then the wick rotated action can be written in terms of the potentials as @xmath181 & = &   \int _
d{{\bf x}}\bigg [ \kappa^2 ( |{\mathbf a}|^2+|\phi|^2 )   + \sum_{j=1}^3 |\nabla a_j|^2 + |\nabla\phi|^2 \nonumber\\ & + & \left({\mathbf a}{\mathbf j}^ * + \phi \rho^ * + { { \mbox{c.c.\,}}}\right)\bigg ] \ , , \end{aligned}\ ] ] where we have chosen the feynman gauge .
the partition function is obtained by integrating over both potentials @xmath174 , @xmath175 and sources @xmath180 , @xmath179 . however , in the latter integration the currents and charges must be weighted according to the energy cost for inducing them on the objects .
this energy cost must depend on shape and material of the objects .
we will see below that this can be achieved by rewriting the self - energies of the separate objects in terms of the incident field that generates the polarizations and magnetizations which give rise to the induced current .
hence the action of eq .
is independent of material and shape of the objects and these properties enter the partition function through proper weights on the currents that measure the susceptibility of the objects to current fluctuations .
we proceed by integrating out the unconstrained fluctuations of the potentials @xmath174 , @xmath175 in the action of eq .  .
this integration yields the partition function as a _ weighted _ functional integral over sources which we indicate at this stage by a subscript @xmath182 on the integration variable .
when we denote the current density in the interior @xmath63 of object @xmath75 by @xmath165 , we get    @xmath183_w [ { { \mathcal d}}{\mathbf j}^*_\alpha]_w   \exp\left[-\frac{1}{2}\sum_{\alpha\beta } \int_{d_\alpha } \!\!d{{\bf x}}\int_{d_\beta } \!\!d{{\bf x } } ' \left\ { { \mathbf j}_\alpha^*({{\bf x } } ) { { \mathcal g}}_>({{\bf x}},{{\bf x}}',\kappa ) { \mathbf j}_\beta({{\bf x } } ' ) + { { \mbox{c.c.\,}}}\right\}\right ]   \nonumber\\ & = & { \mathcal z}_0 \prod_{\alpha=1}^n \int    [ { { \mathcal d}}{\mathbf j}_\alpha]_w   [ { { \mathcal d}}{\mathbf j}^*_\alpha]_w\ , e^{-\tilde s[{\mathbf j } ] } \ , , \end{aligned}\ ] ]    where we have used the continuity equation @xmath184 to eliminate the charge density by introducing the tensor green s function for the half space @xmath185 \ , , \ ] ]    where @xmath58 is the free scalar green s function of eq .  , @xmath186 is the identity matrix and @xmath187 is defined below eq .  .
the action of eq .   can be expressed in terms of the original current densities @xmath165 and the mirror current densities . their components parallel and perpendicular to the mirror plane
are @xmath188 the action then reads @xmath189 & = & \frac{1}{2 } \sum_{\alpha\beta } \left\ {   \int_{d_\alpha } d{{\bf x}}\int_{d_\beta } d{{\bf x } } ' { \mathbf j}^*_\alpha({{\bf x } } ) { { \mathcal g}}_0({{\bf x}},{{\bf x}}',\kappa ) { \mathbf j}_\beta({{\bf x } } ' ) \right . \nonumber\\ & + & \left .
\int_{d_\alpha } \!\!d{{\bf x}}\int_{d^r_\beta } \!\!d{{\bf x } } ' { \mathbf j}^*_\alpha({{\bf x } } ) { { \mathcal g}}_0({{\bf x}},{{\bf x}}',\kappa ) { \mathbf j}^r_\beta({{\bf x } } ' ) + { { \mbox{c.c.\,}}}\right\ } \ , .\nonumber\\\end{aligned}\ ] ] now the sources are coupled by the free , infinite space green s function @xmath190 the mirror currents are located on the mirror objects @xmath64 that are obtained from the original objects by reflection at the mirror plane .
this action has the same structure as in the case of scalar fields .
it can be expressed in terms of multipole moments of the current densities very similarly to the scalar case .
again , we consider diagonal and off - diagonal terms separately .    _ off - diagonal terms _  we introduce again local coordinates @xmath65 that are measured relative to an origin inside @xmath63 .
the terms that couple the original sources on different objects can be written as @xmath191 for @xmath73 and the terms involving mirror sources become @xmath192 for all @xmath75 , @xmath10 .
here we have introduced the electric fields @xmath193 which are generated by the current densities . for positions @xmath84
that are located outside a sphere that encloses the object @xmath69 , the electric field can be expressed in terms of the electric and magnetic multipole moments of the current density @xmath194 for @xmath195 , @xmath196 , @xmath197 where @xmath198 , @xmath199 are the regular , divergence - less solutions @xmath200 \\ { \mathbf n}_{lm}^{j}({{\bf x } } ) & = & \frac{1}{\lambda}\frac{1}{i\kappa } \nabla\times \nabla\times [ { { \bf x}}j_l(i\kappa r)y_{lm}(\hat{{\bf x}})]\end{aligned}\ ] ] of the vector helmholtz equation and @xmath201 . the electric fields of the original sources can then be written as @xmath202 \ , , \ ] ] which is an expansion in the outgoing solutions @xmath203 \\ { \mathbf n}_{lm}^{h}({{\bf x } } ) \!\!&= & \!\!\frac{1}{\lambda}\frac{1}{i\kappa } \nabla\times \nabla\times [ { { \bf x}}h_l^{(1)}(i\kappa r)y_{lm}(\hat{{\bf x}})]\end{aligned}\ ] ] of the vector helmholtz equation . the electric fields of the mirror sources can be also expressed in terms of the multipole moments of the original sources , @xmath204 \ , , \ ] ] where @xmath205 denotes the local coordinates of the mirror object @xmath206 . using eq .  , the definition of the vector solutions of eq .   and
@xmath207 , the moments of the mirror currents are given by @xmath208 the expansion of the electric fields into outgoing vector waves with respect to the origin of the object that generates the field does not allow us to perform the integrations in eqs .  , .
we would like to expand the electric field generated by object @xmath69 in terms of vector waves that are regular at the origin of object @xmath63 .
when the coordinate systems associated with the two objects have identical orientation , this can be done by relating the outgoing and regular vector waves by a translation matrix @xmath209 . generalizing the result for scalar fields , the translation matrix couples both types of vector solutions ,    @xmath210    with @xmath211 .
the matrix elements are given by @xcite , @xmath212 \sqrt{\frac{(2l+1)(2l'+1)(2l''+1)}{l(l+1)l'(l'+1 ) } } \nonumber\\ & & \times\begin{pmatrix}l&l'&l''\\0&0&0\end{pmatrix } \begin{pmatrix}l&l'&l''\\m&-m'&m'-m\end{pmatrix } h^{(1)}_{l''}(i\kappa|{{\bf x}}_{\alpha\beta}| ) y_{l''m - m'}(\hat{{\bf x}}_{\alpha\beta } ) \ , , \\
\label{eq : em - trans - mat - elements - gen - c } c_{l'm'lm}({{\bf x}}_{\alpha\beta } ) & = &   \frac{\kappa}{\sqrt{l(l+1)l'(l'+1 ) } } { { \bf x}}_{\alpha\beta } \cdot \bigg [ \hat{{\bf x}}\frac{1}{2 } \left(\lambda^+_{lm } a_{l'm'lm+1}({{\bf x}}_{\alpha\beta } ) + \lambda^-_{lm } a_{l'm'lm-1}({{\bf x}}_{\alpha\beta } ) \right)\nonumber\\ & & \hspace*{3.6cm}+\ , \hat{{\bf y}}\frac{1}{2i } \left(\lambda^+_{lm }   a_{l'm'lm+1}({{\bf x}}_{\alpha\beta } ) - \lambda^-_{lm } a_{l'm'lm-1}({{\bf x}}_{\alpha\beta } ) \right)\nonumber\\ & & \hspace*{3.6cm}+\ , \hat{{\bf z}}\ , m \,a_{l'm'lm}({{\bf x}}_{\alpha\beta } ) \bigg]\end{aligned}\ ] ] with @xmath213    and @xmath214 .
it is useful to combine the translation matrix elements to the matrix @xmath215 that acts on the vector @xmath216 of multipole moments .    using the translation formula of eqs .
, , we can express the field @xmath217 in eq .   in terms of regular vector solutions and perform the integration . this integration yields the multipole moments defined by eqs .  , .
the action of eq .   can then be written as @xmath218 where we have defined the modified translation matrix @xmath219 \ , , \ ] ] where @xmath220 denotes the conjugate transpose of the matrix of eq .  .
for translations along the z - axis , @xmath221 , the latter expression is diagonal in @xmath149 and simplifies to @xmath222 due to the symmetries @xmath223 and @xmath224 .
the action of eq .   that couples the original and mirror sources can be cast into a form similar to eq .
with a translation matrix @xmath225 that is defined by @xmath226 where the phase factors relating the original and mirror multipole moments in eqs .  , have been absorbed in the translation matrix . with this definition eq .
becomes @xmath227 with the modified translation matrix @xmath228 \ , .\ ] ] for @xmath111 , the action of eq .
describes the interaction between a source and its mirror image . in this case
the translation vector is @xmath229 where @xmath113 is the normal distance between the origin of object @xmath63 and the mirror plane .
hence , the translation is along the z - axis and the translation matrix elements of eqs .
, simplify to    @xmath230 \sqrt{\frac{(2l+1)(2l'+1)}{l(l+1)l'(l'+1 ) } } ( 2l''+1)\nonumber\\ & & \times\begin{pmatrix}l&l'&l''\\0&0&0\end{pmatrix } \begin{pmatrix}l&l'&l''\\m&-m&0\end{pmatrix } h^{(1)}_{l''}(i\kappa 2l_\alpha )   \ , , \\    \label{eq : em - trans - mat - elements - z - c } c_{l'm'lm}(-2l_\alpha\hat{{\bf z } } ) & = & -\delta_{m'm } \kappa 2l_\alpha m ( -1)^m   i^{l - l ' } \sum_{l '' } i^{-l '' }   \sqrt{\frac{(2l+1)(2l'+1)}{l(l+1)l'(l'+1 ) } } ( 2l''+1)\nonumber\\ & & \times\begin{pmatrix}l&l'&l''\\0&0&0\end{pmatrix } \begin{pmatrix}l&l'&l''\\m&-m&0\end{pmatrix } h^{(1)}_{l''}(i\kappa 2l_\alpha ) \ , .\end{aligned}\ ] ]    these matrix elements obey the symmetry relations @xmath231 , @xmath232 so that the matrix of eq .
simplifies for @xmath111 to @xmath233 .    _ diagonal terms _
 so far we have described the interaction between _ arbitrary _ multipoles inside the material objects .
but the functional integral over currents of eq
.   contains weights that measure the energy cost for inducing currents on an object with particular shape and material composition . in the following
we will show that these weights can be implemented straightforwardly when we express the diagonal terms of the action of eq .
, @xmath234 in terms of the incident field which generates a prescribed current that corresponds to the induced polarization and magnetization inside the object .
then the relation between the induced current and the incident field ensures that the currents are weighted properly . to see this we employ the macroscopic formulation of the electromagnetic action , see eq .  .
the energy of eq .   is associated with the process of building up the induced current @xmath165 inside the object .
this energy must therefore be equal to the change in the total macroscopic field energy that results when the object is placed into an external field .
thus we obtain @xmath235 \ , , \ ] ] where integration extends over all space .
field vectors with subscript @xmath17 represent the incident field that is generated by some fixed external sources in otherwise empty space and field vectors without label stand for the total field from the external and induced sources after adding the object .
@xmath236 can be also written as @xmath237 \nonumber\\ & + &   \frac{1}{2 } \int \!d{{\bf x}}\left [ ( { \mathbf e}+{\mathbf e}_0 ) ( { \mathbf d}^*-{\mathbf d}^*_0 ) \right . \nonumber\\ & & \quad \quad\,\,\,\left
. + ( { \mathbf b}+{\mathbf b}_0 ) ( { \mathbf h}^ * -{\mathbf h}_0^ * ) + { { \mbox{c.c.\,}}}\right ] \ , . \end{aligned}\ ] ] the second integral of this expression vanishes .
this can be seen by setting @xmath238 , @xmath239 and using that @xmath240 , @xmath241 since the external sources @xmath242 , @xmath243 are fixed when the object is added .
notice that the bound ( induced ) sources do not appear explicitly since they are included in the macroscopic fields . since @xmath244 , @xmath245 and @xmath246 , @xmath247 with @xmath248 , @xmath249 inside the objects and @xmath250 , @xmath251 outside the objects , we get @xmath252 \ , , \ ] ] where integration runs only over the interior of the object .
the material dependent functions @xmath253 and @xmath254 can be expressed in terms of the polarization @xmath255 and magnetization @xmath256 of the object .
the relations between macroscopic fields yield @xmath257 and @xmath258 . when we substitute @xmath259 in eq .
( [ eq : em - mac - action-2 ] ) we can integrate by parts to obtain @xmath260 \
, , \ ] ] where we have combined the polarization and magnetization to yield the induced current density @xmath261 notice that the incident field @xmath262 in eq .
( [ eq : em - mac - action - final ] ) depends on the current density @xmath165 since @xmath262 has to induce the prescribed current density .
hence , the problem of expressing the diagonal part of the action in terms of multipole moments has been reduced to computing the incident field that has to act on the object to generate a given current density . in scattering theory one
usually encounters the opposite problem . for an incident field
one would like to compute the scattered field which can be expanded in outgoing partial vector waves , see eq .
( [ eq : e - of - mps - orig ] ) . here
the situation is slightly different .
we seek to determine the incident field that generates a given set of multipole moments inside the object .
in other words , for a given scattered field , which according to eq .
( [ eq : e - of - mps - orig ] ) is given by the multipole moments , we would like to obtain the corresponding incident field .
we expand the incident field as @xmath263 \ , .\ ] ] the relation between the multipole moments and the amplitudes of the incident field is determined by the t - matrix @xmath264 , where @xmath265 is the scattering matrix of the object .
when we solve this relation for the incident field amplitudes @xmath266 , we get @xmath267^{-1}_{lml'm ' } { \mathbf q}^\alpha_{l'm ' } \ , , \ ] ] where @xmath268 is a @xmath269 matrix acting on magnetic and electric multipoles . this relation together with eq .
( [ eq : e - inci - expanded ] ) yields the incident field that generates the multipoles @xmath270 . with this result
we can perform the integration in eq .
( [ eq : em - mac - action - final ] ) which yields @xmath236 in terms of the multipoles , @xmath271^{-1}_{lml'm ' } +    \left [ [ \tilde { { \mathcal t}}^\alpha ] ^{-1}\right]^\dagger_{l'm'lm }   \right\ }   { \mathbf q}^{\alpha}_{l'm'}\ ] ] with the modified t - matrix @xmath272 -t^{\alpha,\textsc{em}}_{lml'm'}&t^{\alpha,\textsc{ee}}_{lml'm ' } \end{pmatrix } \ , , \ ] ] where the @xmath273 are the elements of the matrix @xmath268 with @xmath274 , @xmath275 , @xmath276 labeling magnetic and electric elements , respectively . due to the symmetry @xmath277 t^{\alpha,\textsc{em}}_{lml'm'}&t^{\alpha,\textsc{ee}}_{lml'm ' } \end{pmatrix}^ * \!\!=(-1)^{l+l ' } \begin{pmatrix } t^{\alpha,\textsc{mm}}_{l'm'lm}&-t^{\alpha,\textsc{em}}_{l'm'lm}\\[1em ] -t^{\alpha,\textsc{me}}_{l'm'lm}&t^{\alpha,\textsc{ee}}_{l'm'lm } \end{pmatrix}\ ] ] of the t - matrix eq .
( [ eq : em - diag - action - of - mps ] ) can be simplified to @xmath278^{-1}_{lml'm ' }      { \mathbf q}^{\alpha}_{l'm ' }   \ , .\ ] ]    the total action @xmath279 $ ] for the multipole moments is given by the sum of the actions of eqs .
( [ eq : em - off - dia - mps ] ) , ( [ eq : em - off - diag - terms - mirror - mps ] ) and ( [ eq : em - diag - action - of - mps - simp ] ) . in compact matrix notation
it can be written as @xmath280 = \sum_{\alpha\beta } { \mathbf q}^*_\alpha { \mathbb m}^{\alpha\beta } { \mathbf q}_\beta \ , , \ ] ] with the matrix @xmath281^{-1 }   \delta_{\alpha\beta } - \tilde { \mathbb u}^{\alpha\beta } ( 1-\delta_{\alpha\beta } ) + \tilde { \mathbb u}^{r,\alpha\beta } \right\ } \ , , \ ] ] where @xmath282 , @xmath283 and @xmath284 stand for the matrices with matrix elements given by eqs .
( [ eq : em - mod - t ] ) , ( [ eq : em - u - mat - mod ] ) and ( [ eq : em - u - mat - mod - mirror ] ) , respectively . in analogy to the scalar case , eq .
( [ eq : z - og - kappa - final ] ) , the partition function @xmath170 is obtained by integrating over all multipoles .
notice that by construction of the diagonal parts of the action of eq .
( [ eq : em - total_action - final ] ) the proper weights are assigned to the multipole configurations .
the gaussian integral over multipoles and eq .
( [ eq : casimir - energy ] ) lead to the final result @xmath285 for the electromagnetic casimir energy where the determinant is taken with respect to the partial wave indices @xmath148 , @xmath149 , polarization indices @xmath286 , @xmath276 and the object indices @xmath75 , @xmath10 . the matrix @xmath156 is given by eq .
( [ eq : em - def - matrix - m ] ) with the translation matrices @xmath283 and @xmath284 set to zero .    for one object in front of a perfectly reflecting mirror plane , the casimir energy can be written as @xmath287 the matrix @xmath288 that describes the interaction between fluctuating currents on the object and its mirror image is universal , and the shape and material composition of the object enters through its t - matrix @xmath162 .
in this section we consider a geometry that is most relevant to a large number of experimental studies of casimir interactions carried out in the last decade .
this geometry , consisting of a plane mirror and a sphere , is experimentally favorable since it avoids the problem of parallelism for plane surfaces facing each other . despite its experimental importance
, this geometry lacks a theoretical description of the electromagnetic casimir interaction . for a scalar field with dirichlet boundary conditions at the sphere and the plane , the interaction over a wide range of distances has been obtained recently @xcite , including an analytic expression for the lowest order correction to the proximity force approximation @xcite .
we consider a sphere of radius @xmath289 in front of a plane mirror with distance @xmath290 between the center of the sphere and the mirror , see fig .
[ fig : plane - sphere ] .
for the scalar field we study dirichlet or neumann boundary conditions at the mirror and the sphere .
for the electromagnetic field , the sphere consists of a material with arbitrary dielectric function @xmath291 and permeability @xmath292 while the mirror is assumed to be perfectly reflecting ( implying vanishing of the parallel components of the electric field and the normal component of the magnetic field ) .
we start by studying the casimir interaction at large separations . for separations @xmath290 that
are large compared to the size @xmath289 of the object , the casimir energy can be expressed as an asymptotic series in @xmath293 . using that @xmath294 in eqs .
( [ eq : casimir - energy-1-surface ] ) and ( [ eq : em - casimir - energy-1-surface ] ) , we get by expanding the logarithm @xmath295 where we have defined the operator @xmath296 with @xmath52(@xmath53 ) for a scalar field in the presence of a dirichlet ( neumann ) mirror and @xmath297 for the electromagnetic field and a perfectly reflecting mirror . the operator @xmath298 describes a wave that travels from the mirror to the object and back , involving one scattering at the object .
all information about the shape of the surface and the boundary conditions is provided by the @xmath132-matrix . for a spherically symmetric surface
the matrix is diagonal and completely specified by scattering phase shifts @xmath299 that do not depend on @xmath149 , @xmath300 where @xmath301 is the real frequency .
the scattering phase shifts for a sphere of radius @xmath289 with dirichlet and neumann boundary conditions are @xmath302 where @xmath303 and @xmath98 ( @xmath304 ) are spherical bessel functions of first ( second ) kind . to make use of the general result of eq .
, we have to evaluate the matrix elements for imaginary frequencies @xmath305 .
this can be done using @xmath306 , @xmath307 , where @xmath308 and @xmath309 are modified bessel functions of first and second kind , respectively .
we obtain @xmath310 for dirichlet and neumann boundary conditions where @xmath311 .    the distance dependence of the casimir energy of eq .
is given by the translation matrix that describes the interaction between the original and mirror multipole moments .
since the origins ( centers ) of the original and mirror sphere are related by a translation along the @xmath312-axis , the matrix elements of the translation matrix @xmath313 are given by eqs .   and .
after a rotation of frequency to the imaginary axis we obtain    @xmath314    with the elements of the transition and translation matrices we can compute the casimir interaction from eq .  .
the scaling of the @xmath315-matrix for small @xmath32 shows that partial waves of order @xmath148 start to contribute to the energy at order @xmath316 if the @xmath132-matrix is diagonal in @xmath148 .
hence , we can truncate the infinite matrix @xmath298 in eq .   at finite multipole order to obtain the large distance expansion of the energy at a given order in @xmath293 .
also , we can truncate the series over @xmath317 in eq .   since the @xmath318 power of @xmath298 becomes important only at order @xmath319 .
the first @xmath317 powers of @xmath298 describe @xmath317 scattering off the sphere .
we have used mathematica to perform the matrix operations and to expand the integrand of eq .   in @xmath293 . from this
we find that the casimir energy can be written as @xmath320 where @xmath321 is the coefficient of the term @xmath322 . in the following
we give the results for a mirror plane with dirichlet or neumann boundary conditions separately .
_ dirichlet mirror _  in this case we have @xmath323 .
if the field obeys dirichlet boundary conditions at the sphere , we obtain the for leading coefficients @xmath324 indicating an attractive force .
this result agrees with recent findings of wirzba @xcite .
if we impose neumann boundary conditions at the sphere , we find @xmath325 this result corresponds to a repulsive force which is expected for unlike boundary conditions . to leading order the energy scales as @xmath326 since the first two coefficients vanish .
this behavior can be understood from the absence of low - frequency @xmath327-wave scattering for a sphere with neumann boundary conditions .    _
neumann mirror _  in this case we have @xmath328 .
for a dirichlet sphere we obtain the coefficients @xmath329 and hence the force is again repulsive .
the modulus of the coefficients are smaller than those in eq .   for a dirichlet mirror with the exception of the leading one , @xmath330 and the one with inverted sign , @xmath331 .
notice that the sum of the casimir energies for the dirichlet and neumann mirror , both opposite to a dirichlet sphere , is identical to the casimir energy between two dirichlet spheres at a center - to - center distance @xmath332 @xcite .
the reason for this is that any field configuration can be decomposed into symmetric and antisymmetric modes with respect to the mirror plane which obey neumann and dirichlet boundary conditions on that plane , respectively .
since the energy between two spheres scales as @xmath333 for asymptotically large separations , the modulus of the coefficients @xmath330 in eqs .
( [ eq : coeff_d - d ] ) and ( [ eq : coeff_n - d ] ) must be identical so that the contribution @xmath334 cancels .
it is easily seen that the higher order coefficients of eqs .
( [ eq : coeff_d - d ] ) and ( [ eq : coeff_n - d ] ) combine to the correct coefficients of the large distance expansion for two dirichlet spheres @xcite .    for a neumann sphere
we get @xmath335 for the same reason as before , the first two coefficients vanish .
notice that all given coefficients , with the exception of @xmath331 , are equal to minus the corresponding coefficients for a neumann sphere and a dirichlet mirror in eq .
( [ eq : coeff_d - n ] ) .
this result is consistent with the observation that the casimir energy for two neumann spheres at a center - to - center distance @xmath332 is given by the sum of the energies of a neumann sphere at distance @xmath290 from a dirichlet mirror and a neumann mirror , respectively .
since two neumann spheres interact at large distance with an energy @xmath336 with the next - to - leading term @xmath337 , the sum of the corresponding coefficients , excepting @xmath331 , must vanish .
the coefficients @xmath331 in eqs .
( [ eq : coeff_d - n ] ) , ( [ eq : coeff_n - n ] ) combine to the correct value for two neumann spheres that was obtained in ref .
@xcite .
we consider a dielectric sphere at a center - to - surface distance @xmath290 from a perfectly conducting mirror . due to
spherical symmetry , the electric and magnetic multipoles for all @xmath148 , @xmath149 are decoupled so that the t - matrix is diagonal @xcite ,    @xmath338 - n i_{l+{1\over 2}}(nz ) \left[i_{l+{1\over 2}}(z)+2z i'_{l+{1\over 2}}(z)\right ] } { \eta k_{l+{1\over 2}}(z ) \left[i_{l+{1\over 2}}(nz)+2nzi'_{l+{1\over 2}}(nz)\right ] - n i_{l+{1\over 2}}(nz ) \left[k_{l+{1\over 2}}(z)+2z k'_{l+{1\over 2}}(z)\right ] } \ , , \ ] ]    where the sphere radius is @xmath289 , @xmath311 , @xmath339 , @xmath340 .
@xmath341 is obtained from eq .
( [ eq : t - matrix - elem - sphere ] ) by interchanging @xmath342 and @xmath343 .
the limit of a perfectly conducting sphere is obtained by taking @xmath344 at an arbitrarily fixed @xmath345 which can also vanish .
then the matrix elements become independent of @xmath343 , @xmath346 for all partial waves , the _ leading _ low frequency contribution is determined by the _
static _ electric multipole polarizability , @xmath347r^{2l+1}$ ] , and the corresponding magnetic polarizability , @xmath348r^{2l+1}$ ] . including the next to leading terms
, the t - matrix has the structure @xmath349 \ , , \nonumber\ ] ] and @xmath341 is obtained by @xmath350 , @xmath351 .
the first terms are @xmath352/[5(\mu+2)^2]r^5 $ ] , @xmath353 ^ 2r^6 $ ] , and @xmath354 , @xmath355 are obtained again by interchanging @xmath343 and @xmath342 .
higher order terms can be easily obtained by expanding eq .
( [ eq : t - matrix - elem - sphere ] ) for small @xmath32 .
the translation matrix elements @xmath356 are obtained from eqs .
( [ eq : def - em - u - matrix - mirror ] ) , ( [ eq : em - trans - mat - elements - z - b ] ) and ( [ eq : em - trans - mat - elements - z - c ] ) .
using @xmath121 we get    @xmath357 \sqrt{\frac{(2l+1)(2l'+1)}{l(l+1)l'(l'+1 ) } } ( 2l''+1)\nonumber\\ & & \times\begin{pmatrix}l&l'&l''\\0&0&0\end{pmatrix } \begin{pmatrix}l&l'&l''\\m&-m&0\end{pmatrix }   \frac{k_{l''+1/2}(2\kappa l)}{\sqrt{\pi \kappa l } }   \ , , \\
\label{eq : em - trans - mat - elements - z - c - wk } c_{l'm'lm}(-2l\hat{{\bf z } } ) & = & \delta_{m'm }   2 \kappa l m
( -1)^m   i^{l - l ' } \sum_{l '' } ( -1)^{l '' }   \sqrt{\frac{(2l+1)(2l'+1)}{l(l+1)l'(l'+1 ) } } ( 2l''+1)\nonumber\\ & & \times\begin{pmatrix}l&l'&l''\\0&0&0\end{pmatrix } \begin{pmatrix}l&l'&l''\\m&-m&0\end{pmatrix } \frac{k_{l''+1/2}(2\kappa l)}{\sqrt{\pi
\kappa l } } \ , .\end{aligned}\ ] ]    now we can employ the series representation of the casimir energy in eq .  ( [ eq : energy - scatt - exp ] ) with @xmath358 . since the t - matrix is diagonal in @xmath148 , partial waves of order @xmath148 start to contribute to the energy at order @xmath316 .
also , the @xmath318 power of @xmath298 becomes important only at order @xmath359 .
notice the stronger increase of the exponent with @xmath317 compared to the scalar case where the exponent is @xmath360 .
this can be understood from the absence of @xmath327-waves for the electromagnetic field so that each reflection contributes a factor @xmath361 due to @xmath317-waves .
matrix operations and the expansion in @xmath293 are performed with mathematica . from this
we obtain for a dielectric sphere in front of perfectly conducting mirror plane the casimir energy    @xmath362 \frac{1}{l^7 } + \frac{7}{7200 } \left [ 572 ( \alpha_3^\textsc{e}-\alpha_3^\textsc{m } ) + 675 \left ( 9 (   \gamma_{15}^\textsc{e } -   \gamma_{15}^\textsc{m } ) -55    (   \gamma_{23}^\textsc{e } -   \gamma_{23}^\textsc{m } ) \right)\right ] \frac{1}{l^8 } + \dots \right\ } \ , .\end{aligned}\ ] ]    the electric contribution to the leading term , @xmath363 , was obtained by casimir and polder for the interaction of an atom with static polarizability @xmath364 and a metallic surface @xcite .
later boyer has generalized the leading order result to include magnetic effects described by @xmath365 @xcite .
the higher order terms are new .
they show how higher order polarizabilities and frequency corrections to the static parameters influence the interaction .
there is no @xmath366 term .
notice also that the first three terms of the contribution at order @xmath367 have precisely the structure of the casimir - polder interaction between two atoms with static dipole polarizabilities @xmath365 and @xmath364 but it is reduced by a factor of @xmath368 .
this factor and the distance dependence @xmath369 of this term suggests that it arises from the interaction of the dipole fluctuations inside the sphere with those inside its image at a distance @xmath332 .
the additional coefficient of @xmath370 in the reduction factor @xmath371 can be traced back to the fact that the forces involved in bringing the dipole in from infinity act only on the dipole and not on its image @xcite .    for a perfectly conducting sphere the coefficients of the expansion of the casimir energy in @xmath293 become universal numbers .
these coefficients can be either obtained from eq .
( [ eq : energy - eps - mu - sphere ] ) by using the appropriate values of the parameters for a prefect metal or it can be computed directly from the t - matrix elements given in eqs .
( [ eq : t - matrix - elem - cond - sphere - m ] ) , ( [ eq : t - matrix - elem - cond - sphere - e ] ) . following the latter route
, we get the series @xmath372 where the coefficients up to order @xmath373 are @xmath374 this and the corresponding results for a scalar field appear to be asymptotic series . therefore , the series can not be summed to obtain the interaction at small separations . in the next section
we use a numerical implementation of eq .
( [ eq : em - casimir - energy-1-surface ] ) to compute the interaction at all separations .      to obtain the casimir interaction over a broad range of distances , eq .
( [ eq : em - casimir - energy-1-surface ] ) has to be evaluated numerically . as we employ a partial wave expansion , the computational work increases with decreasing separation between the objects .
however , we shall see below that even at small separations our method yields sufficient precision to obtain the leading corrections to the proximity force approximation ( pfa ) .
the numerical approach is based on the technique presented in refs .
@xcite . using the analytic expressions for the matrix elements of the translation and transition matrices , we compute the determinant and the integral over frequency @xmath32 in eq .
( [ eq : em - casimir - energy-1-surface ] ) numerically .
the matrices are truncated at a finite multipole order @xmath148 which yields a series of estimates @xmath375 for the casimir energy .
the exact casimir energy is then obtained by extrapolating the series to @xmath376 . for the geometry considered here
, we observe an exponentially fast convergence that allows us to obtain accurate results even at small separations from a moderate multipole order .
we compare our results for the interaction energy to the estimate that follows from pfa .
this is important since pfa is often used over a range of separations although its accuracy is unknown even at small distances for most geometries , including the electromagnetic casimir interaction between a plane and a sphere .
the pfa estimate for the latter geometry is given by @xmath377 where @xmath378 is the surface - to - surface distance .
the amplitudes follow from the result for two parallel plates and are given by @xmath379 for a scalar field with like ( @xmath52 ) and unlike ( @xmath53 ) boundary conditions and the electromagnetic field , respectively . in the following we present our results for the plane - sphere interaction over a wide range of separations for scalar fields with dirichlet and neumann boundary conditions and for the electromagnetic field .      ) .
the dashed curve represents the large distance expansion of eq .
( [ eq : spheres - energy - large - l ] ) with the coefficients of eq .
( [ eq : coeff_d - d ] ) .
inset : corrections to the pfa at small separations as function of the surface - to - surface distance @xmath378 .
the dashed curve corresponds to the lowest order correction to the pfa obtained by bordag @xcite . ]     for neumann boundary conditions at the sphere and the plane .
the dashed curve represents the large distance expansion of eq .
( [ eq : spheres - energy - large - l ] ) with the coefficients of eq .
( [ eq : coeff_n - n ] ) . ]     for neumann boundary conditions at the sphere and dirichlet boundary conditions at the plane .
the dashed curve represents the large distance expansion of eq .
( [ eq : spheres - energy - large - l ] ) with the coefficients of eq .
( [ eq : coeff_d - n ] ) . ]     for dirichlet boundary conditions at the sphere and neumann boundary conditions at the plane .
the dashed curve represents the large distance expansion of eq .
( [ eq : spheres - energy - large - l ] ) with the coefficients of eq .
( [ eq : coeff_n - d ] ) . ]
we consider the four cases that correspond to dirichlet or neumann boundary conditions at the sphere and the mirror plane .
the results for the casimir energy are shown in figs .
[ fig : dd]-[fig : dn ] . at large separations
the results of the numerical evaluation of the determinant agree nicely with the large distance expansions of the previous section . for a sphere with dirichlet boundary conditions ,
the energy scales both at large and small separations as @xmath380 so that the curves for @xmath381 in figs .
[ fig : dd ] and [ fig : dn ] tend also to a constant for @xmath382 .
the pfa underestimates the actual interaction energy at all separations . for a neumann sphere
the interaction scales at large separations as @xmath326 so that the curves for @xmath381 in figs .
[ fig : nn ] and [ fig : nd ] tend to zero for @xmath382 . here
the pfa overestimates the actual interaction energy at all separations .
the non - perturbative approach allows us to study also the case of small separations . in that limit
our results can be fitted to a power law of the form @xmath383 \ , .\ ] ] the coefficients @xmath384 obtained from a fit of the function of eq .
( [ eq : e - pfa+corrections ] ) to the data points for the four smallest studied separations are summarized in tab .
[ tab : thetas ] .
the fitted curves are shown as insets in figs .
[ fig : dd]-[fig : dn ] .
the result for @xmath385 for the case of dirichlet boundary conditions at the sphere and the mirror agrees with the analytical result @xmath386 presented in ref .
again for dirichlet conditions , for the second order coefficient a considerably larger numerical estimate of @xmath387 has been obtained from world - line monte carlo sampling @xcite . for the other combinations of boundary conditions
our findings represent the first results for the corrections to the pfa .
.coefficients describing corrections to pfa , see eq .
( [ eq : e - pfa+corrections ] ) .
d = dirichlet , n = neumann ( scalar field ) and em = perfectly conducting boundary conditions ( electromagnetic field ) . [ cols="^,^,^ " , ]      we focus on a perfectly conducting sphere and mirror plane . the casimir energy resulting from a numerical computation of the determinant of eq .
( [ eq : em - casimir - energy-1-surface ] ) is shown for a wide range of separations in fig .
[ fig : em ] . at large separations
the interaction is described by the asymptotic result of eq .
( [ eq : em - energy - series ] ) . at small distances
the casimir energy approaches to pfa estimate and corrections to pfa can be described again by eq .
( [ eq : e - pfa+corrections ] ) .
a corresponding fit to the data points for the four smallest separations is shown as inset in fig .
[ fig : em ] .
the corresponding amplitudes of the correction terms are listed in tab .
[ tab : thetas ] .    ) .
the asymptotic expansion of eq .
( [ eq : em - energy - series ] ) to order @xmath388 is shown as dashed line .
inset : corrections to the pfa at small distances as function of @xmath378 . ]
this result is important for a number of recent measurements of the casimir force between almost flat surfaces .
these experiments have been performed for a plane mirror and a sphere with a radius that is much larger than the distance between the surfaces in order to avoid difficulties from parallelism control .
our results at short distances indicate that for a sphere of radius @xmath389 m the corrections to pfa are below @xmath390 only for surface - to - surface distances @xmath391 that are smaller than @xmath392 nm . for a sphere of that size , pfa fails already by @xmath393 at a distance of @xmath394 m . from this
we conclude that the currently achieved experimental accuracy for measurements of casimir forces in the sphere - plate geometry with a sphere of @xmath389 m at distances below @xmath395 m is within the range or slightly less than the corrections to pfa .
however , deviations from pfa become severe when smaller objects interact with surfaces .
hence it is important to understand the crossover of the casimir interaction between macroscopic objects and the eventual casimir - polder interaction between single atoms and a surface .
a description of this crossover is provided by our results in fig .
[ fig : em ] .
corresponding results can be also obtained for a dielectric sphere or less symmetric objects by the methods presented here .
finally , we note that it should be also possible to compute correction amplitudes like the @xmath384 in eq .
analytically by applying methods similar to those used in ref .
however , the validity range of the corresponding series would be limited to small separations and for an overall description of the interaction one should resort to the device of a numerical evaluation along the lines presented here .
this work emerged from a collaboration with r.  l.  jaffe , m.  kardar , and n.  graham on related problems .
support by the heisenberg program of the deutsche forschungsgemeinschaft is acknowledged .
a. rodriguez , m. ibanescu , d. iannuzzi , f. capasso , j. d. joannopoulos and s. g. johnson , phys .
lett . * 99 * , 080401 ( 2007 ) ; a. rodriguez , m. ibanescu , d. iannuzzi , j. d. joannopoulos , and s. g. johnson , phys .
a * 76 * , 032106 ( 2007 ) .