the quantum vacuum is a fundamental concept in theoretical physics and its properties has been widely investigated in the literature of quantum gravity and string theory .
the theory of quantum fields in curved spacetime which is believed to be the low energy limit of the ultimate theory of quantum gravity , has predicted famous quantum effects in the presence of gravity . in general , due to the lack of the global symmetries in the spacetime manifold , quantum effect considerations in curved spacetime are mainly limited to the analysis of the local quantities such as the vacuum expectation value of the energy - momentum tensor , i.e. @xmath0 , in some point @xmath1 .
in fact , the most famous results of the semi - classical theory of gravity like the hawking radiation and the particle production in the expanding universe has been achieved from the analysis of the @xmath2 in the related curved backgrounds .
the most famous vacuum state effect is the casimir effect .
an important aspect in the researches around the casimir effect in curved spacetime is that the characteristics of the vacuum state are apparently dependent to the geometry of the background spacetime .
we are also motivated to see explicitly such a dependency in this paper .
furthermore , computation of the energy , i.e. @xmath3 , has also been done for large number of problems in various spacetimes and sometimes @xcite has helped us to confirm the validity of the principle of correspondence in the context of the casimir effect . finding the total gravitational force on a set of two conducting casimir plates @xcite is a typical example .
so , due to the importance of the stress - energy tensor , we will consider the @xmath2 for the plates in a general weak background .
the casimir effect arises when there is a boundary in our problem and it predicts a force between two uncharged conducting metals in the presence of a quantum field .
the effect has been measured to a great accuracy @xcite .
we use the zero - point energy approach here although it is possible to find the casimir force and energy without any reference to the zero - point energy @xcite .
we may also have the casimir effect without having a boundary at all .
in fact , some non - trivial topologies in curved spacetime do the same job as a boundary does@xcite,@xcite .
the casimir energy in curved spacetime has been also analysed by many authors ( @xcite and references there in ) . recently ,
a casimir apparatus consisting of two ideal conducting parallel plates in the weak field limit of the kerr and the horava - lifshitz spacetimes has been studied in @xcite .
a purpose of this paper is to generalize the above analysis for scalar field doing exact solution of the klein - gordon equation in a general weak gravitational field .
also we extend the method for the case the electromagnetic field is present inside the plates .
the organization of the paper is as follows . in section
ii using the fact that the apparatus has composed of tiny pales , we will find the small deviations of the metric from flat spacetime inside the apparatus .
in fact we will expand the metric up to first order in terms of the parameters of the general weak gravitational field .
in section iii the klein - gordon equation will be solved exactly inside the apparatus using the metric obtained in the previous section . in section
iv mode frequencies inside the apparatus will be obtained under the influence of both dirichlet and neumann boundary conditions for the scalar field .
generalization to the case where the electromagnetic field is present will be done using an interesting property of the fermi spacetime in section v. computation of the energy momentum tensor for the scalar field for both neumann and dirichlet boundary conditions are done in section vi .
also the electromagnetic energy density is obtained in this section .
well known weak gravitational fields are examined in section vii .
the electromagnetic energy density for the far field limit of the kerr spacetime and the horava - lifschitz theory of gravity are of special interest .
the final section is devoted to the conclusion .
in the 1 + 3 formalism of general relativity , stationary spacetime metric is defined by @xmath4 where @xmath5 is the so called gravitomagnetic potential and @xmath6 in the weak field slowly rotating limit ( @xmath7 ) , the metric ( [ eq01 ] ) is equivalent to @xmath8 the explicit form of a general weak gravitational field line element is as follows ( see ( 19.13 ) in @xcite ) @xmath9dt^2-[4\epsilon_{ijk}s^j\frac{x^k}{r^3}+o(\frac{1}{r^3})]dx^idt \\ & -\{(1+\frac{2m}{r}+\frac{3m^2}{2r^2})\delta_{ij}+gravitational\;radiation\ ; terms\;that\;die\ ; out\ ; as\ ; o(\frac{1}{r } ) \}dx^i dx^j , \end{split}\end{aligned}\ ] ] comparison between ( [ eq03 ] ) and ( [ eq04 ] ) shows that @xmath10 is the newtonian potential , @xmath11 is the gravitomagnetic potential and @xmath12 .
we need to have an isotropic coordinate representation of the metric in the next sections and equations ( [ eq01]),([eq02 ] ) and ( [ eq03 ] ) shows us the way we can construct an isotropic coordinate representation of any weak gravitational field .    .
the boundaries are at z=0,z = l ]    fig . [
fig : apparatus ] demonstrates the apparatus in a weak gravitational field .
two ideal conducting plates are in a small coordinate distance @xmath13 from each other . inside the plates ,
the metric has a small variation relative to the flat spacetime metric . to find this variation ,
we adapt a rectangular coordinate system having the origin at one of the plates ( the one which is closer to the source and have coordinate distance @xmath14 from it ) and expand the metric ( [ eq04 ] ) inside the apparatus in the neighborhood of the point @xmath15 .
the overall size of the apparatus is so small that we can assume @xmath16 : @xmath17 in which @xmath18 . for the case of a static spacetime , the components of the metric ( [ eq05 ] )
can be written in the following form : @xmath19 where @xmath20 are constants and use is made of @xmath21 . concerning the form of @xmath22 , the above expansion satisfies @xmath23 provided that @xmath24 i.e. in the equatorial plane .
so ( [ eq05]),([eq06 ] ) are also valid for the case of the far field limit of the kerr spacetime and we will back to it in the examples in the section vi .    motivated by the above discussion , in this paper we analyse the general case of the spacetime of the form @xmath25,\end{aligned}\ ] ] in which @xmath26 .    to solve the klein - gordon equation , it is better to recast the metric ( [ eq07 ] ) to the known fermi metric .
we use the linearized weak field regime of general relativity and change the variables with the aid of the following gauge transformation : @xmath27 which we have assumed @xmath28 to force the spatial sector of the metric ( [ eq07 ] ) to be flat .
more explicitly we have : @xmath29 the metric then takes the following form up to first order in the parameters @xmath30 @xmath31 the gauge transformation , however , changes our primary problem as follows .
in fact according to the last equation in ( [ eq08 ] ) the boundaries must be transformed from @xmath32 and @xmath33 in the spacetime ( [ eq07 ] ) to @xmath34 and @xmath35 in the spacetime ( [ eq10 ] ) .
another change that the gauge transformation brings into the problem is the rescaling of time in ( [ eq09 ] ) by the factor @xmath36 .
this corresponds , in turn , to dividing the mode frequencies @xmath37 by @xmath36 because of the presence of the factor @xmath38 in our solution of the klein - gordon equation in the next section .
so time must be re - inverted after solving the klein - gordon equation when we obtain the mode frequencies . _ the net effect of the rescaling of time is that the final mode frequencies must be multiplied by the factor @xmath36_. we drop the dashes @xmath39 on @xmath40 from now on using the new boundary conditions instead .
the massless klein - gordon equation is @xmath41= 0 \;\;\;\;\;\;,\;\;\ ; \mathfrak{g}\equiv { \rm det}g^{\mu\nu}.\end{aligned}\ ] ] since the spacetime is spatially flat we assume the following form for the solution @xmath42 where @xmath43 is normalization constant determined through the commutation relations @xmath44 the scalar product is defined as @xmath45^{\frac{1}{2}}n^{\mu } d\sigma.\end{aligned}\ ] ] in which @xmath46 and @xmath47 spans the space between the plates . under the above assumptions equation ( [ eq11 ] ) reads @xmath48 where @xmath49 denotes derivation with respect to z and @xmath50 .
another variable change @xmath51 yields @xmath52 appearance of factor 3 in front of the second term , introduces a significant simplification when we change the variable to @xmath53 .
it recasts ( [ eq16 ] ) into @xmath54 a simple reparametrization of this last equation via @xmath55 end up with @xmath56 this is the known kummer s differential equation @xmath57 in which @xmath58 .
kummer s differential equation is not suitable for next considerations and we transform it via @xmath59 to another known form called wittaker s differential equation : @xmath60 in which in our case @xmath61 .
equation ( [ eq20 ] ) has two independent set of solutions @xmath62 , @xmath63 and their asymptotic behaviour is as follows @xcite    @xmath64    @xmath65    in general , @xmath63 is the acceptable physical solution which is finite at infinity . in the problem under consideration
we must choose a linear combination of the two , due to the fact that both are finite in between the plates .
the exact mode functions are : @xmath66 e^{-i\omega t - ik_xx -ik_yy}\;\;\;\;\;\;\;\;\\end{aligned}\ ] ] the asymptotic form of ( [ eq23 ] ) for small value of @xmath67 can be written as follows ( see appendix a ) @xmath68 in which @xmath69 . in the next section we use this asymptotic form to extract mode frequencies for both neumann and dirichlet boundary conditions on the plates .
we put the approximation @xmath70 into ( [ eq24 ] ) . from the dirichlet boundary condition @xmath71 we have @xmath72 and from @xmath73 we have : @xmath74 after careful expansion of ( [ eq25 ] ) , the mode frequencies proved to satisfy the following relation : @xmath75\}=\sqrt{k_\perp^2+(\frac{n\pi}{l+\gamma_1 l+ \frac{1}{2 } \lambda_1 l^2})^2}\end{aligned}\ ] ] note that the factor @xmath76 is nothing but the _ proper distance _ between the plates and so we have @xmath77 in which @xmath78 denotes proper ( or the corresponding flat space ) mode frequencies in the local lorentz frame of an observer comoving with the plates .
as stated in the previous section , the final mode frequencies will be obtained by multiplication of the factor @xmath36 due to rescaling of time during the gauge transformation ( [ eq09 ] ) .
so we have the final mode frequencies inside the casimir apparatus for the spacetime ( [ eq07 ] ) : @xmath79 mode frequencies are influenced only by @xmath80 component of the metric from the point of view of a proper observer .
the neumann boundary condition @xmath81 imposed on ( [ eq24 ] ) gives @xmath82 : @xmath83 another neuman boundary condition @xmath84 end up with @xmath85 which in turn result in ( [ eq26 ] ) and ( [ eq27 ] ) again .
the electromagnetic field has two physical degrees of freedom and it is known that in the rindler spacetime the electromagnetic field can be represented in terms of two scalar fields satisfying the klein - gordon equation separately @xcite .
the photon propagator and the energy - momentum tensor of the electromagnetic field in a weak gravitational of the fermi spacetime has also been obtained in @xcite .
however in @xcite , the computations has been done through a lengthy and cumbersome method of the green functions . in a paper by the author @xcite
, it is shown that the energy density ( i.e. the 0 - 0 component of the energy - momentum tensor ) of the electromagnetic field in fermi spacetime is exactly the some of the energy density of the two scalar fields mentioned in @xcite .
the method was done without any address to the green function method frequently used in the literature .
here , we briefly review the relationship between the two scalar fields and the electromagnetic field in fermi spacetime .
the spin one vector field in curved spacetime in the lorentz gauge satisfies @xmath86 it can be shown that in fermi metric , the ricci tensor satisfies @xmath87 and the second term of the wave equation ( [ eq31 ] ) must be ignored . furthermore , because the metric ( [ eq10 ] ) is spatially flat , the lorentz gauge in ( [ eq31 ] ) can be broken into two independent parts @xcite : @xmath88 in which @xmath89 and @xmath90 we know also that both @xmath91 and @xmath92 satisfy the klein - gordon eqaution separately(see the appendix in @xcite ) .
boundary condition for the electric field on the plates is @xmath93 and @xmath94 which in turn can be recast into boundary conditions on @xmath91 and @xmath92 according to equation ( [ eq33 ] ) . in @xcite it has been shown that boundary conditions for the electric and magnetic fields return a dirichlet boundary condition on @xmath92 and a neumann boundary condition on @xmath91 .
we have shown in previous section that both of this conditions will be ended up to a same frequency shift . as a result mode frequencies in ( [ eq28 ] ) are also valid for the electromagnetic field .
this section has three subsections . in the first subsection
, the electromagnetic energy - momentum tensor will be represented in terms of the energy momentum of the scalar fields mentioned in ( [ eq33 ] ) . in the other two sections
the energy momentum tensor of the scalar and vector fields will be calculated .
the vacuum expectation value of the quantum energy - momentum tensor is defined as @xmath95.\end{aligned}\ ] ] the classical energy - momentum tensors for the scalar and vector fields are @xmath96 , \\
t_{\mu\nu}^{gauge } & = -c^*_{;\mu}c_{;\nu}-c^*_{;\nu}c_{;\mu}-g_{\mu\nu}g^{\lambda\theta}c_{;\lambda}c_{;\theta } , \\
t_{\mu\nu}^{maxwell } & = \frac{1}{4}g_{\mu\nu}f^{\lambda\theta}f_{\lambda\theta}-f_{\mu^\theta}f_{\theta\nu}. \end{split}\end{aligned}\ ] ]    as is well - known , in the quantum level , the contributions of ghost and gauge fields in the electromagnetic energy - momentum tensor cancel each other @xcite and we only concern the maxwell sector of the energy - momentum tensor .
expansion of @xmath97 in terms of scalar fields @xmath91 and @xmath92 shows that ( see appendix b ) @xmath98 in which @xmath99 . using @xmath100 and ( [ eq32 ] )
we have : @xmath101 in which we have used the general form of the wave function ( [ eq12 ] ) .
quadratic products of fields @xmath102 and @xmath103 produce terms like @xmath104 which have no contribution when the expectation value is taken because of the fact that @xmath91 and @xmath92 are not correlated and belongs to independent hilbert spaces .
we calculate @xmath105 component of the energy - momentum tensor for the electromagnetic field : @xmath106 @xmath107 \ } \\ & = k_\perp^2[<0|t_{00}^{\phi}|0>+<0|t_{00}^{\psi}|0 > ] \end{split}\end{aligned}\ ] ] after a lengthy but straightforward calculation , other components of the energy - momentum tensor in both sides are related to each other as follows : @xmath108 note that @xmath109 can be absorbed into @xmath43 .
this section is devoted to the calculation of the energy density for the casimir apparatus via the direct method without any reference to the traditional green function method . using the approximations
@xmath110 we expand the wave function ( [ eq24 ] ) and find up to first order in @xmath111 : @xmath112 sin(\sqrt{b}z)+\frac{a}{4\sqrt{b}}z^2 cos(\sqrt{b}z)\ } \;\;\ , dirichlet \\
z(z ) & = z_0 \{[1-(\frac{\lambda}{2}+\frac{a}{4b})z ] cos(\sqrt{b}z)-\frac{a}{4\sqrt{b}}z^2 sin(\sqrt{b}z)\ } \;\;\ , newmann \end{split}\end{aligned}\ ] ] @xmath113 can be absorbed also in @xmath114 .
the energy density is defined by @xmath115 where @xmath116 is the lapse vector normal to the hypersurface @xmath117 i.e. @xmath118 .
the mean energy density so has the following form : @xmath119 calculating the factor @xmath120 for both wave functions in ( [ eq45 ] ) results in @xmath121   \\ & + \frac{a}{4\sqrt{b}}[\frac{z}{2b}sin(2\sqrt{b}z ) + ( \frac{1}{2b}-z^2)\frac{cos(2\sqrt{b}z)}{2\sqrt{b } } ]   \}|_\alpha^\beta , \\
f_2 = & \{\frac{z}{2}-\frac{sin(2\sqrt{b}z)}{4\sqrt{b}}-(\lambda+\frac{a}{2b})[\frac{z^2}{4}-\frac{z sin(2\sqrt{b}z)}{4\sqrt{b}}-\frac { cos(2\sqrt{b}z)}{8b } ] \\ & + \frac{a}{4\sqrt{b}}[\frac{z}{2b}sin(2\sqrt{b}z)+(\frac{1}{2b}-z^2)\frac{cos(2\sqrt{b}z)}{2\sqrt{b } } ]   \}|_\alpha^\beta ,   \\
f_3 = & \ { b[\frac{z}{2}+\frac{sin(2\sqrt{b}z)}{4\sqrt{b}}+(\frac{a}{2b}-\lambda)[\frac{z^2}{4}+\frac{z sin(2\sqrt{b}z)}{4\sqrt{b}}+\frac { cos(2\sqrt{b}z)}{8b } ] ] \\ & -\frac{a\sqrt{b}}{4 } [ \frac{z}{2b}sin(2\sqrt{b}z)+(\frac{1}{2b}-z^2)\frac{cos(2\sqrt{b}z)}{2\sqrt{b } } ] + ( \frac{\lambda}{4}+\frac{a}{8b } ) cos(2\sqrt{b}z)\}|_\alpha^\beta .
\\ \end{split}\end{aligned}\ ] ] @xmath122   \\ & -\frac{a}{4\sqrt{b}}[\frac{z}{2b}sin(2\sqrt{b}z ) + ( \frac{1}{2b}-z^2)\frac{cos(2\sqrt{b}z)}{2\sqrt{b } } ]   \}|_\alpha^\beta , \\ f_2 = & \{\frac{z}{2}+\frac{sin(2\sqrt{b}z)}{4\sqrt{b}}-(\lambda+\frac{a}{2b})[\frac{z^2}{4}+\frac{z sin(2\sqrt{b}z)}{4\sqrt{b}}+\frac { cos(2\sqrt{b}z)}{8b } ] \\ & -\frac{a}{4\sqrt{b}}[\frac{z}{2b}sin(2\sqrt{b}z)+(\frac{1}{2b}-z^2)\frac{cos(2\sqrt{b}z)}{2\sqrt{b } } ]   \}|_\alpha^\beta ,   \\
f_3 = & \ { b[\frac{z}{2}-\frac{sin(2\sqrt{b}z)}{4\sqrt{b}}+(\frac{a}{2b}-\lambda)[\frac{z^2}{4}-\frac{z sin(2\sqrt{b}z)}{4\sqrt{b}}-\frac { cos(2\sqrt{b}z)}{8b } ] ] \\ & + \frac{a\sqrt{b}}{4 } [ \frac{z}{2b}sin(2\sqrt{b}z)+(\frac{1}{2b}-z^2)\frac{cos(2\sqrt{b}z)}{2\sqrt{b } } ] -(\frac{\lambda}{4}+\frac{a}{8b } ) cos(2\sqrt{b}z)\}|_\alpha^\beta . \end{split}\end{aligned}\ ] ]    simplification of the terms like @xmath123 and @xmath124 is possible using the fact that in boundaries the wave function ( [ eq24 ] ) must vanish and : @xmath125 in which @xmath126 takes one of two boundary values @xmath32 and @xmath127 . equation ( [ eq47 ] )
result in @xmath128 and the final result is as follows @xmath129(\beta^2-\alpha^2 ) \}.\end{aligned}\ ] ] exactly the same @xmath130 is obtained for dirichlet boundary condition .
the constant @xmath43 was defined in ( [ eq07 ] ) and can be determined simply as : @xmath131 sin^2(\sqrt{b}z)+\frac{a}{4\sqrt{b}}z^2 sin(2\sqrt{b}z)\}dz\}^{-1 } \;\ , dirichlet \\
c^2 & = \frac{1}{2(2\pi)^2\omega } \{\int_0^l \{[1-(2\lambda+\frac{a}{2b})z ] cos^2(\sqrt{b}z)-\frac{a}{4\sqrt{b}}z^2 sin(2\sqrt{b}z)\}dz\}^{-1 } \;\ , neumann \end{split}\end{aligned}\ ] ] just like @xmath130 , the constant @xmath43 has also the same form for both boundary conditions up to first order in @xmath111 @xmath132 the final result for the energy - density after all is : @xmath133 in which @xmath134 is the corresponding flat spacetime casimir energy density @xmath135 @xcite .    up to now
, we have shown that for neumann and dirichlet boundary conditions there exists the relation @xmath136 between flat and curved energy density contents of the casimir apparatus in the weak spacetime of the metric ( [ eq07 ] ) .      using the first equation of ( [ eq41 ] ) it is evident that the same shift in the energy density obtained in ( [ eq52 ] ) holds also for the electromagnetic field : @xmath137 in the next section we will analyse some well - known weak gravitational fields and find the parameters @xmath138 in each case .
the problem of the divergences near a perfect generic conductor first studied systematically by deutsch and candelas @xcite .
they found that the energy - momentum tensor near the surface behaves like : @xmath139 where @xmath140 are constants and @xmath141 is the distance from the surface of the ideal boundary . in case where
there is a conformal invariance in the action @xmath142 .
the divergences originates from the unphysical nature of classical ideal conductor boundary conditions .
it has been shown @xcite that we can remove the infinities of the total energy ( and not the energy - momentum tensor ) of the plates using the zeta function regularization unless for the case that the zeta function has poles itself .
on the other hand , the cut off regularization method suggests the removal of the divergences via _ ad hoc _ although in this method it still remains a logarithmic ambiguity @xcite in the energy density .
for the imperfect conductors ( more realistic boundary conditions ) we can easily remove the divergences introducing some suitable cut off frequencies although the boundary effect may become quit large ( but finite)@xcite .
now the question is that what happens when we go to curved spacetime ? .
does the surface divergences of the energy - momentum tensor are ignorable in the semi - classical einstein s equations ? .
the answer is negative according to @xcite . in @xcite
however , the authors has found a way to get ride of the surface divergences of the einstein s field equations for the case the boundary is a parallelepiped . they have used a suitable cut off along with the so - called estrada - kanwal distribution theory of asymptotics to regularize / renormalize the infinities and show that the energy - momentum tensor near a plane boundary , as a source , converges to a consistent theory when the cutoff is removed .
remarkably , the process of curing the divergences in curved spacetime has been continued by milton et al @xcite where they have shown that most of the energy between the plates is restored near the surfaces and a part of it resides exactly on the plates .
they finally have shown that this energies responds to gravity just like any other finite energy following the newtonian relation @xmath143 .
this is expectable as this large energies are simply a part of the total energy of the system .
we show here that , in our case , the divergences do not present in the first order of approximation that we have used here and they appear only in the higher orders of approximation . to do so , we write down the explicit structure of possible divergences in the einstein s field equations along the lines depicted in @xcite .
the one - loop effective action @xmath144 for the semi - classical theory of gravity is : @xmath145 in which @xmath146 is the dewitt - schwinger - feynmann s propagator . using the dewitt - schwinger representation of the action ,
the asymptotic expansion of @xmath147 is as follows : @xmath148 in which @xmath149 is a length scale to fix the dimensional issues .
the potentially divergent part of the effective action ( the first three terms ) @xcite reads : @xmath150 \ } [ \frac{1}{n(n-2)}4m^2a_0-\frac{1}{(n-2)}2m^2a_1+a_2 ] \\.
\end{split}\end{aligned}\ ] ] in which @xmath151 is the mass of the scalar field(m=0 in our case ) .
the coefficients @xmath152 are : @xmath153 in which @xmath154 is the riemann s curvature tensor .
the total gravitational lagrangian density can be shown to have the form ( see ( 6.49 ) in @xcite ) : @xmath155 \ } a_2 , \\ & a=\frac{4m^4}{(4\pi)^{\frac{n}{2 } } n(n-2 ) } \{\frac{1}{n-4}+\frac{1}{2}[\gamma+\ln{\frac{m^2}{\mu^2 } } ] \ } \\ & b=\frac{m^2}{(4\pi)^{3\frac{n}{2 } } ( n-2 ) } \{\frac{1}{n-4}+\frac{1}{2}[\gamma+\ln{\frac{m^2}{\mu^2 } } ] \}. \\ \end{split}\end{aligned}\ ] ] however , in the massless case of ours , the only non - vanishing potentially ultraviolent term is the one related to @xmath156 in ( [ eq52 ] ) and @xmath157 in ( [ eq52 ] ) vanish ( see ( 6.101 ) in @xcite ) .
our calculation for the metric ( [ eq07 ] ) shows that @xmath156 is of second order of approximation : @xmath158 so the total lagrangian density ( [ eq59 ] ) reduces to the the standard bare density @xmath159 . in conclusion , the potentially divergent term @xmath156 vanishes within the first order of approximation
in this section number of spacetimes are investigated and the parameters appeared in ( [ eq28]),([eq52]),([eq53 ] ) will be found . to this end , we will try to find their weak field form according to ( [ eq03]),([eq04 ] ) and ( [ eq07 ] ) .
recently bezerra et al @xcite studied the apparatus for the scalar fields in the weak field limit of the kerr spacetime in the equatorial plane . in that work the apparatus co -
rotates with the local angular velocity of the spacetime i.e. the measurements had been assumed to be done in the point of view of a zero angular momentum observer(zamo ) .
they found the metric inside the apparatus through some two stage successive approximation method as follows : @xmath160\;\;\ ; 0-th\ ; order\ ; of\ ; approximation , \\ ds^2 \approx & ( 1 + 2b\lambda z)dt^2 - ( 1 -2\lambda z)[dx^2+dy^2+dz^2]\;\;\ ; first\;order\;of\;approximation . \end{split}\end{aligned}\ ] ] in which @xmath161 , @xmath162 and @xmath163 .
@xmath164 is the angular momentum per mass and @xmath165 the local angular velocity of the kerr spacetime .
the authors set @xmath166 for the first order of approximations .
this does not works because @xmath111 is related to @xmath167 through @xmath168 .
putting one of them equal to zero forces the other one to vanish also .
evidently , in the first order of approximations we must keep both @xmath167 and @xmath169 and the metric must be written as follows instead of ( [ eq60 ] ) : @xmath170\;\;\ ; 0-th\ ; order\ ; of\ ; approximation , \\ ds^2 \approx & [ 1 + 2b(\phi_0+\gamma z)]dt^2 - [ 1 -2(\phi_0+\gamma z)][dx^2+dy^2+dz^2]\;\;\ ; first\ ; order\ ; of\ ; approximation . \end{split}\end{aligned}\ ] ] taking the above comment and the metric ( [ eq61 ] ) into account , we find @xmath171 and arrive at the following relations for the e__lectromagnetic casimir energy density _ _ and the mode frequencies : @xmath172 , \\ < 0|t_{00}^{maxwell}|0 > & = -\frac{\pi^2}{720 l^4}[1+b\phi_0+b\frac{\gamma l_p}{2 } ] .
\end{split}\end{aligned}\ ] ] the scalar field casimir energy density inside the apparatus which has been sketched in eq.(32 ) in @xcite must also be corrected as follows : @xmath173 , \\
< 0|t_{00}^{scalar}|0 > & = -\frac{\pi^2}{1440 l^4}[1+b\phi_0+b\frac{\gamma l_p}{2 } ] .
\end{split}\end{aligned}\ ] ] _ far field limit of the schwarzschild spacetime _ is also covered by setting @xmath174      the fermi spacetime is described by the following metric @xmath175 the importance of this metric is that it is traditionally recognized as the spacetime of a static accelerating observer near the surface of the source of a constant gravitational field @xcite .
comparison of this metric to the general metric in ( [ eq07 ] ) gives @xmath176 : @xmath177 , \\
< 0|t_{00}^{maxwell}|0 > & = -\frac{\pi^2}{720 l^4}[1+\frac{a l_p}{2 } ] , \\
< 0|t_{00}^{scalar}|0 > & = -\frac{\pi^2}{1440 l^4}[1+\frac{a l_p}{2 } ] .
\end{split}\end{aligned}\ ] ] the above energy densities are exactly the results ( 5.2 ) in @xcite , ( 3.4 ) in @xcite and ( 5.4 ) in @xcite .
the hoava - lifshitz ( hl ) gravity is a renormalizable theory of gravity that is invariant under the lifshitz scaling transformation @xmath178 .
this transformations manifestly break the space and time covariance .
the anisotropy between space and time , in turn , may affects the casimir effect as well .
it is interesting to investigate the vacuum characteristics of the theory .
recently , the effect of the hl theory on the casimir energy of the apparatus has been studied in @xcite .
the authors recommended to set a constraint on spacetime anisotropies in such a way that the casimir energy modifications remain within the experimental bounds .
recently in @xcite , @xcite the same problem considered in curved spacetime in the context of a spherical symmetric solution of the hl theory .
finite temperature casimir energy in spacetime ( [ eq07 ] ) has been analysed by the author in @xcite and the following weak field limit for the hl theory has been calculated : @xmath179 + [ -2\frac{\widehat{m}}{r^2}-\frac{3\widehat{m}^2}{2r^3}+\frac{\widehat{m}^2}{2\widehat{\omega } r^5}]z\ } dt^2 \\ & -\{1+[-\frac{\widehat{m}}{r}+\frac{\widehat{m}^2}{4r^2}+\frac{\widehat{m}^2}{4\widehat{\omega}r^4 } ] + [ \frac{\widehat{m}}{r^2}-\frac{\widehat{m}^2}{4r^3}-\frac{\widehat{m}^2}{4\widehat{\omega}r^5}]z \ } ( dx^2+dy^2+dz^2 ) \end{split}\end{aligned}\ ] ] where @xmath180 .
this spacetime is the weak field limit of the park s spherical symmetric solution to the ir limit of the hl theory _ in the presence of a cosmological constant@xcite _ @xcite .
@xmath181 may have the same role as the cosmological constant but not necessarily being a small parameter and @xmath37 is a constant frequently used to regulate the uv limit of the hl theory .
comparison between ( [ eq66 ] ) and the general metric ( [ eq07 ] ) shows : @xmath182 another spacetime , which is a solution to the hl theory _ without cosmological constant _ is the kehagias - sfetsos ( ks ) solution .
this spacetime has been discussed in @xcite to obtain the energy density of the apparatus and is as follows : @xmath183 where @xmath184 .
@xmath185 is a radial coordinate and @xmath37 is the free parameter of the hl theory .
putting @xmath186 in the park s solution recover the ks solution and so the coefficients in equation ( [ eq67 ] ) are also valid for the ks solution .
we analysed the energy density of a casimir apparatus consisting of two nearby conducting parallel plates in a general weak gravitational field .
the metric in the equation ( [ eq07 ] ) denotes the deviation of the weak gravitational field from flat spacetime inside the apparatus .
we transformed the metric ( [ eq07 ] ) through a gauge transformation into the fermi metric and then solved the klein - gordon equation exactly .
the mode frequencies were found for the scalar field inside the apparatus for both neumann and dirichlet boundary conditions in terms of the weak gravitational field parameters @xmath138 .
this result was shown to be valid also for the electromagnetic field in section v. the energy density of the apparatus was found for both scalar and electromagnetic fields in terms of the weak field parameters .
some examples of weak gravitational fields were analysed in section vii .
specially the electromagnetic energy density and mode frequencies in the far field limit of the kerr spacetime in its equatorial plane were obtained .
the weak field limit of the hoava - lifshitz gravity with a cosmological constant was also investigated and the weak field parameters were sketched .
consistency of the results with the literature was checked by considering the fermi metric .
in this section we find an explicit and simple asymptotic form for the wave function .
as is apparent we need to have an asymptotic expansion with both argument and first parameter being large .
whittaker functions have such an expansion in terms of airy functions @xcite : @xmath187\sum_{n=0}^{\infty}\frac{a_n(\zeta)}{(4\kappa)^{2n}}+\frac{ai^{'}[(4\kappa)^{\frac{2}{3}}\zeta]}{(4\kappa)^{\frac{4}{3}}}\sum_{n=0}^{\infty}\frac{b_n(\zeta)}{(4\kappa)^{2n}}\}\end{aligned}\ ] ] where @xmath188 is defined as @xmath189 and in our case @xmath190 and so @xmath191 if @xmath192 then both the argument and the first parameter go to infinity . as @xmath193
the second term in the bracket is of order @xmath194 and must be ignored to stay within the first order of expansion in @xmath111 .
also the summation in the first term reduces to only @xmath195 term and @xmath196 .
@xmath197\end{aligned}\ ] ] the argument of the above airy function can be written as @xmath198\}^{\frac{2}{3}}\equiv - v(x ) \;\;\ , \;\;\
v\rightarrow\infty\end{aligned}\ ] ] the airy function with large argument is as follows ( section 9.7 from @xcite ) @xmath199 as @xmath200 , the @xmath201-part and @xmath202 in the @xmath203-part must be ignored also .
the observation that @xmath204 and @xmath205 moves the situation forward significantly : @xmath206 a same process can be applied to the other whittaker s function @xmath62 except that according to @xcite instead of the a argument @xmath207 in ( [ a5 ] ) we have @xmath208 .
this difference changes the above process as follows : @xmath209 @xmath210 as a result the total wave function in ( [ eq22 ] ) when @xmath211 reads @xmath212 where @xmath213 and @xmath157 came from ( [ eq22 ] ) and the following relation has been used : @xmath214
in this section we find @xmath215 .
other components will be find by similar techniques .
we do the computations in the fermi spacetime i.e. @xmath216 .
greece indices run from 0 to 3 and latins from 1 to 3 . from ( [ eq36 ] )
we have @xmath217 the first term has the following form : @xmath218 where we have used @xmath219 .
the second term in ( [ b1 ] ) simplifies as follows : @xmath220b^bb^n \\ &
\;\;\;\;\;\;\;\;\;\;\;= ... =-g^{ij}b^2+b^ib^j .
\end{split}\end{aligned}\ ] ] based on ( [ b1])-([b3 ] ) we find @xmath221 after lowering the indices we have finally : @xmath222    now we find @xmath223 . from ( [ eq36 ] )
we have @xmath224 in which the first term already has been obtained in ( [ b2 ] ) .
the second term has been also obtained in ( [ b3 ] ) and after lowering the indices again we have @xmath225    now we find @xmath223 . from ( [ eq36 ] )
we have @xmath226 from which , after lowering the indices again , we find @xmath227
the author would like to thank university of tehran for supporting this project under the grants provided by the research council .
30 k. a. milton et al . , j.phys.a 40,10935 - 10943,2007 ; s. a. fulling , k. a. milton , p. parashar , a. romeo , k.v .
shajesh , j. wagner phys .
d 76 , 025004 ( 2007 ) k. a. milton , k. v. shajesh , s. a. fulling , prachi parashar , phys .
d 89 , 064027 ( 2014 ) .
s. lamoreaux , phys .
78 , 5 ( 1996 ) ; g. bressi et al .
88 041804 ( 2002 ) .
a. jaffe , phys .
d 72 021301 ( 2005 ) .
l. h. ford , phys .
d 14 , 3304 ( 1976 )
. b. s. dewitt , phys .
19 , 295 ( 1975 ) .
bezerra , h.f.mota , c.r .
phys.rev.d 89,044015(2014 ) . f. sorge , phys .
d 90 , 084050 ( 2014 ) .
muniz , v.b .
bezerra , m.s .
cunha , phys .
d 88 , 104035(2013 ) .
muniz , v.b .
bezerra , m.s .
cunha , annals of physics 359 ( 2015 ) 55 - 63 .
b. geyer , g. l. klimchitskaya , and v. m. mostepanenko phys .
a 67 , 062102 ( 2003 ) ; m. bordag , g.l .
klimchitskaya , u. mohideen , v.m .
mostepanenko , advances in casimir effect , oxford university press , 2009 ; m. bordag , proceedings of the fourth workshop on quantum field theory under the influence of external conditions : the casimir effect 50 years later , world scientific , 1998 .
g. bimonte et al . , phys .
d , 74 ( 2006 ) 085011;phys .
d 78 ( 2008 ) 024010 ; erratum ibid .
d 75 , 049904(2007 ) ; erratum ibid .
d 75 , 089901 ( 2007 ) ; erratum ibid .
d 77 , 109903 ( 2008 ) .
f. sorge , class .
( 2005 ) 5109 - 5119 .
g. esposito , g. m. napolitano , l rosa , phys.rev.d 77,105011(2008 ) .
g. m. napolitano , g. esposito , l rosa , phys.rev.d 78,107701(2008 ) .
m. nouri - zonoz , b. nazari , phys .
d 82 ( 2010 ) 044047 .
borzoo nazari , m. nouri - zonoz , phys .
d 85 ( 2012 ) 044060 .
j. s. dowker , r. critchley , j.phy.a : math.gen
. 9 535(1976 ) .
d. lynden - bell , m. nouri - zonoz , rev .
phys.70 , 427 ( 1998 ) l. d. landao , e. m. lifshitz , the classical theory of fields , pergamon press , 1987 ( reprint ) n.d .
birrell , p.c.w .
davies , quantum field theory in curved spacetime , cambridge universsity press , 1982 .
p. hoava , phys . rev .
d 79 , 084008 ( 2009 ) .
m. park , jhep 0909 ( 2009 ) 123 .
ferrari et al , mod .
lett . a 29 ( 2013 ) 1350052 .
misner , k. s. thorne , j. h. wheeler , gravitation , w.h .
freeman and company , 1973 .
p. deutsch , d. candelas , proc .
a. 354 , 79 - 99 ( 1977 ) .
p. deutsch , d. candelas , phys .
d 20 , 3063(1979 ) 79 - 99 ( 1977 ) .
r. estrada , s. a fulling , l. kaplan , k. kirsten , z. liu and k. a. milton , j. phys .
a : math . theor .
41 ( 2008 ) 164055 ; r. estrada , s. a fulling , f. d mera , j. phys . a : math .
theor . 45 ( 2012 ) 455402 .
olde daalhuis , nist handbook of mathematical functions , cambridge university press , new york , ny , 2010 , ch.13 i.s .
gradshteyn , i.m .
ryzhik , table of integrals , series and products , 7th edition , elsevier inc .
, 2007 milton abramowitz , irene a. stegun , handbook of mathematical functions , 10th printing , national bureau of standards us , 1972 , ch.13 w. rindler , relativity : special , gerneral and cosmological , cambridge university press , 2006 ( chapter 9 ) f. w. j. olver , asymptotics and special functions , academic press , 1974 .
( chapter 11 , ex .
7.3 ) a. erdlyi , c. a.swanson asymptotic forms of whittaker s confluent hypergeometric functions .
1957 ( 1957 ) , no .
25 , 49 pp .
( equation 10.3 ) borzoo nazari , thermal corrections to the casimir effect in a weak gravitational field , under review