1 Introduction

The Casimir force, which acts on two (uncharged) conducting parallel plates in vacuum is, in general, attractive [1, 2]. Casimir force can be explained as a consequence of quantum vacuum fluctuations, but it can be alternatively understood in terms of Van der Waals interactions [3]. Although this force is attractive in vacuum and depends only in the plates separation, a different behavior can be manifested under other conditions. When a medium made of different materials (or metamaterials), with diverse properties or geometrical shapes, is inserted between the plates, a repulsive Casimir force may appear (see for instance, Refs. [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]). The Casimir force can even produce a torque in the plates in optically anisotropic materials [19].

Recently, Jiang and Wilczek [4] have shown that a tunable repulsive and attractive Casimir force can also be obtained when electromagnetic waves propagate in chiral material media between two conducting parallel plates. In such cases, the behavior of the Casimir force can be controlled by varying the electric permittivity and magnetic permeability of the material. In this chiral material, right– and left–circularly polarized electromagnetic waves propagate with wavevectors \(k_+\) and \(k_-\) respectively. In this material, both wavevectors are described by \(k_\pm =k\pm \delta k\), meaning that the difference between the wavevectors of the two polarized waves is given by \(k_+-k_-=2\delta k\). This difference is responsible for producing tunable repulsive or attractive Casimir forces when the separation l between plates is varied.

In such cases, the Casimir energy is sensitive to this phase difference, and it turns out to be [4]

$$\begin{aligned} E_C= & {} \int _0^\infty \frac{d\varpi }{2\pi }\int _{-\infty }^\infty \frac{d^2 k_\parallel }{4\pi ^2} \nonumber \\&\quad \ln \left[ {1+e^{-4\varrho l}-2e^{-2\varrho l}\cos (2\delta k l)}\right] \, , \end{aligned}$$
(1)

where \(\varpi =i\omega \), \(\varrho =({\varpi ^2+k_\parallel ^2})^{1/2}\), \(k_\parallel =({k_1^2+k_2^2})^{1/2}\), with the wave frequency \(\omega \), and the wavevector components \(k_i\) (\(i=1,2,3\)). From here, the Casimir force per unit area \(\mathcal{{F}}_C=- d E_c/dl\) experienced between two plates is readily found to be

$$\begin{aligned} {\mathcal{{F}}_C}= & {} \frac{1}{2\pi ^3}\int _0^\infty d\varpi \int _{-\infty }^\infty d^2 k_\parallel \nonumber \\&\times \left( \frac{\varrho e^{-4\varrho l}-\varrho e^{-2\varrho l}\cos (2\delta k l)-\delta k e^{-2\varrho l}\sin (2\delta k l)}{1+e^{-4\varrho l}-2e^{-2\varrho l}\cos (2\delta k l)}\right) \, .\nonumber \\ \end{aligned}$$
(2)

This force can be either attractive or repulsive depending on the magnitude of \(2\delta k l\), and it has been shown in Ref. [4] that repulsive Casimir force may appear in Faraday and Optical active materials. On the contrary, in a non-chiral material, as there is no difference between phase velocities of right– and left–circularly polarized waves \(\delta k=0\), the force (2) becomes simply \(-{\pi ^2}/({240 l^4})\), recovering the value for vacuum Casimir force.

In Ref. [4], only material media were considered. The main purpose of this work is to show that gravitational fields described by stationary spacetimes (metrics with rotation, such as Kerr black holes or Gödel universe), or by anisotropic cosmological spacetimes (such as Bianchi–I metric), produce a behavior similar to the one appearing in chiral media when electromagnetic waves propagate on it.

Related gravitational effects due to spacetime curvature on Casimir force have been explored previously [20,21,22,23,24,25,26,27,28]. Also, the effect of spacetime curvature correction due to Kerr metric on Casimir forces has been studied for both scalar [29,30,31,32,33] and electromagnetic fields [34]. Our work explores how the coupling of polarization with the angular momentum or anisotropy of spacetime modifies the Casimir force. We show this effects by describing the analogy between electromagnetic wave propagation in curved spacetimes with the electromagnetic fields in media. We obtain solutions for Kerr, Gödel and Bianchi–I spacetimes, where this coupling creates a difference between right– and left–circularly polarized electromagnetic waves. We finally show how the chiral behavior of light on those spacetimes modifies the Casimir force.

2 Maxwell equations in curved spacetime

We start by describing the dynamics of electromagnetic fields on a gravitational background field (from now on we use natural units \(c=\hbar =1\)). In general, covariant Maxwell equations in curved spacetime may be written as

$$\begin{aligned} \nabla _\alpha F^{\alpha \beta }=0\, ,\quad \nabla _\alpha F^{*\alpha \beta }=0\, , \end{aligned}$$
(3)

in terms of the antisymmetric electromagnetic field tensor \(F^{\alpha \beta }\) and its dual \(F^{*\alpha \beta }\). Here, \(\nabla _\alpha \) is the covariant derivative defined by a metric \(g_{\mu \nu }\).

For the case of electromagnetic waves, several articles [35,36,37,38,39,40,41,42] have shown that light does not propagate along null geodesics in curved spacetimes, giving rise to a dispersion relation that depends on their polarization. This occurs because gravitational fields behave as effective material media with non–trivial effective permeability and effective permittivity, both of them modifying the (vacuum) Maxwell equations. Thus, the electromagnetic field amplitude and polarization (besides its phase), couple to spacetime curveture. This is explicitly shown by defining the corresponding electric \(E_i=F_{i0}\), \(D^i=\sqrt{-g} F^{0i}\), and magnetic \(B^i=\varepsilon ^{0ijk} F_{jk}\), \(\varepsilon ^{0ijk}H_k=\sqrt{-g} F^{ij}\) fields, where \(\varepsilon ^{0ijk}\) is the Levi-Civita symbol and g the metric determinant with latin indices used to denote (three–dimensional) space coordinates. The above vector fields are related by [35,36,37, 39, 40, 43]

$$\begin{aligned} D^i= & {} \epsilon ^{ij}E_j-\varepsilon ^{0ijk}\mu _j H_k\, \nonumber \\ B^i= & {} \epsilon ^{ij}H_j+\varepsilon ^{0ijk}\mu _j E_k\, , \end{aligned}$$
(4)

explicitly showing the analogy with electric and magnetic fields in the presence of a medium. Here

$$\begin{aligned} \epsilon ^{ij}=-\sqrt{-g}\frac{g^{ij}}{g_{00}}\, , \quad {\mu _j}=-\frac{g_{0j}}{g_{00}}\, , \end{aligned}$$
(5)

are the effective permittivity and the effective vector permeability of the curved spacetime, respectively. Both are defined in terms of spacetime metric (and its inverse \(g^{\mu \nu }\)). Using the above electric and magnetic fields, Maxwell equations (3) can now be written in a vectorial fashion as

$$\begin{aligned} \partial _i D^i= & {} 0\, ,\quad \partial _0 D^i=\varepsilon ^{0ijk}\partial _j H_k\, ,\nonumber \\ \partial _i B^i= & {} 0\, , \quad \partial _0 B^i=-\varepsilon ^{0ijk}\partial _j E_k\, , \end{aligned}$$
(6)

where \(\partial _0\) and \(\partial _i\) stands for the time and spatial partial derivatives. Equation (6) are completly equivalent to the covariant Maxwell equations (3).

An alternative and enlightening representation of electromagnetic fields in curved spacetime [35, 41] can be obtained by defining the (generalized) Riemann–Silberstein vectors

$$\begin{aligned} F_{j\pm }= & {} E_j\pm i H_j\, ,\nonumber \\ S_\pm ^j= & {} D^j\pm i B^j\, , \end{aligned}$$
(7)

where the \(+\) (−) symbol represents right (left) polarization of the electromagnetic fields. Using relations (4), we obtain that

$$\begin{aligned} S_{\pm }^j= \epsilon ^{ij} F_{j\pm }\pm i \varepsilon ^{0ijk}\mu _jF_{k\pm }\, , \end{aligned}$$
(8)

and therefore, Maxwell equations (6) reduce to the simplest form

$$\begin{aligned} \partial _j S^j_\pm= & {} 0\, ,\nonumber \\ \pm i\, \partial _0 S_\pm ^j= & {} \varepsilon ^{0jkm}\partial _k F_{m\pm }\, . \end{aligned}$$
(9)

We can now use Riemann–Silberstein polarization vectors, to finally write Maxwell equations curved spacetime as [35]

$$\begin{aligned} \pm i\, \partial _0\left( \epsilon ^{jk} F_{k\pm }\right) =\varepsilon ^{0jkm}\left[ \partial _k F_{m\pm }+\partial _0\left( \mu _k F_{m\pm }\right) \right] \, . \end{aligned}$$
(10)

Equations (10) are completely equivalent to Eq. (3). From here is clear that, in general, both polarization states (±) propagate differently. Spacetime couples with different components of electromagnetic fields (this fact does not occur in vacuum), even producing rotation of the polarization state of light [37, 40, 42]. Isotropic and symmetric spacetimes (with \(\mu _j=0\) and \(\epsilon ^{ij}=\epsilon \delta ^{ij}\)) allow null geodesic light propagation [37, 42], while a general curved spacetime behaves as an effective medium for travelling electromagnetic waves. Nevertheless, when \(\mu _j\ne 0\) spacetime acts as a chiral material. Similarly, when \(\epsilon ^{ij}=\epsilon ^i \delta ^{ij}\), with \(\epsilon ^i\ne \epsilon ^j\) for \(i\ne j\), spacetimes produce birefringence [37, 42].

3 Polarized electromagnetic waves in rotating and anisotropic spacetimes

It can be straightforwardly shown that Kerr and Gödel curved spacetimes, both having \(\mu _j\ne 0\), produce different solutions for electromagnetic waves with right- or left-polarization. On the other hand, cosmological anisotropy with \(\epsilon ^{ij}\approx \epsilon ^i\delta ^{ij}\) and \(\mu _j=0\), produces birefringence that also modifies the propagation of polarized waves.

Thereby, these effects induce changes in the magnitude of Casimir force, which are directly associated to the interaction of light polarization and gravitational fields.

3.1 Electromagnetic waves propagating in a slowly Kerr black hole background

For the sake of simplicity, we focus in the case of an electromagnetic wave travelling in the exterior of the gravitational field of a slowly rotating Kerr black hole with mass M and spin J. This problem was solved by Mashhoon [35], showing that these propagating electromagnetic waves experience chiral effects in such a metric due to coupling with their polarization.

For a slowly rotating black hole, the stationary Kerr metric in isotropic coordinates becomes [35, 44]

$$\begin{aligned} ds^2\approx & {} -\left( 1-\frac{2M}{r}\right) dt^2-\frac{4J}{r^3}\left( x_1 dx_2-x_2dx_1\right) dt\nonumber \\&+\delta _{ij}dx^i dx^j\, , \end{aligned}$$
(11)

for a rotation in the \(x_3\)-direction, and with \(r=|\mathbf{x}|\). Its spin (angular momentum of magnitude J) is \(J^j=\delta ^{j3} J\) in these coordinates. This intrinsic angular momentum is measured far from the black hole [44]. Using metric (11) in (5), we get that

$$\begin{aligned} \epsilon ^{ij}\approx \delta ^{ij}\left( 1+\frac{2 M}{|\mathbf{x}|}\right) \approx \delta ^{ij}\, ,\quad \mu _j\approx -2\varepsilon _{0jkm}\frac{J^k x^m}{|\mathbf{x}|^3}\, , \end{aligned}$$
(12)

are the permittivity and permeability of the effective medium for an electromagnetic wave propagating far from the black hole, with \(|\mathbf{x}|\sim D\gg M\), where D is the given distance where the light propagates. For our slowly rotating metric, \(|\mu |\sim J/D^2\ll 1\). Therefore, Maxwell equations (10) reduce to [35]

$$\begin{aligned}&\pm \omega \, \delta ^{jk}f_{k\pm }=\varepsilon ^{0jkm}\left[ \partial _k -i\, \omega \, \eta \, \left( x_2 \delta _{1k}- x_1\delta _{2k}\right) \right] f_{m\pm } \, ,\nonumber \\ \end{aligned}$$
(13)

where we have used \(F_{j\pm }(t, \mathbf{x})=f_{j\pm }(\mathbf{x}) e^{-i\omega t}\), with wave frequency \(\omega \). Besides, \(\eta =2 J/D^3\). Furthermore, the scale of the electromagnetic wave is much less than the one of the black hole, such that \(\omega M\gg 1\).

Solutions to Eq. (13) represent right– and left–circularly polarized electromagnetic waves propagating in the \(x_2-x_3\) plane, very far from the Kerr black hole (where spacetime is almost flat). Performing the change of variables \(f_{j\pm }(\mathbf{x}) =\exp (i\omega \eta x_1 x_2)\Psi _{j\pm }(\mathbf{x})\), from Eq. (13) we find

$$\begin{aligned} \mp i \, \omega \, \Psi ^j_{\pm }=-i\, \varepsilon ^{0jkm}\partial _j \Psi _{m\pm }+ 2\omega \,\eta \, x_1 \varepsilon ^{0jkm} \delta _{2k} \Psi _{m\pm } \, . \end{aligned}$$
(14)

This equation can be readily solved by using a new change of variables \(\Psi _{j\pm }(\mathbf{x}) =\exp (i k_2 x_2+i k_3^\pm x_3)\phi _{j\pm }(x_1)\), where \(k_{2}\) and \(k_3^\pm \) are the wavevectors in directions \(x_2\) and \(x_3\) respectively, for the corresponding polarizations. Each component of \(\phi _{j\pm }\) satisfies different equations. The component along the \(x_3\)-direction has the solution \(\phi _{3\pm }(\xi )\propto \exp (-\xi ^2/2) H_{n^\pm }(\xi )\) [35], where \(\xi =(2\eta \omega x_1+k_2)/\sqrt{2\eta \omega }\), and \(H_{n^\pm }\) is the Hermite polynomial, with \(n^\pm =0,1,2,\ldots \). For this solution, the dispersion relation

$$\begin{aligned} \omega ^2-\left( k_3^\pm \right) ^2\pm 2\eta \, k_3^\pm =2\eta \,\omega (2n^\pm +1)\, , \end{aligned}$$
(15)

must be satisfied. Solutions for the remaining components \(\phi _{1\pm }\) and \(\phi _{2\pm }\) are also known, and we refer the reader to Ref. [35] in order to find their thorough mathematical study.

Dispersion relation (15) defines the propagation of an electromagnetic wave in Kerr metric (and in a Gödel universe as we will see). As \(k_3^+\ne k_3^-\), this relation shows that right– and left–circularly polarized waves move with different phase velocities when propagating along the direction of the rotation axis of rotation of curved spacetime. In general, we find that the difference between the wavevectors is

$$\begin{aligned} k_3^+-k_3^-= & {} 2\delta k=2\eta +\sqrt{(\omega -\eta )^2-4\eta \omega n^+}\nonumber \\&\quad -\sqrt{(\omega -\eta )^2-4\eta \omega n^-}\, , \end{aligned}$$
(16)

which is a function of the frequency of the wave. However, because of our previous assumption on the nature of the black hole and the scales of electromagnetic waves, we always find that \({\omega }/{\eta }={\omega M}{(M/D)^{-1}(J/D^2)^{-1}}\gg 1\), describing waves in the high-frequency limit. In this case, the physical solutions for the electromagnetic wave components \(\phi _{1\pm }\) and \(\phi _{2\pm }\) imply that \(n^+-1=n^-+1\) [35]. All these allow us to find that the final wavelength difference becomes simply

$$\begin{aligned} 2\delta k= & {} -2\eta \, . \end{aligned}$$
(17)

Difference (17) is a physical consequence of the reaction of wave polarization to the rotation axis of the spacetime. Each spin of the circularly polarized wave couples differently to curvature. It is not an effect of coordinate choice, as the angular momentum of Kerr metric is well-defined as a measure of the frame dragging of this spacetime [44, 45]. The above result shows that these left– and right–circularly polarized electromagnetic waves propagate in an fashion analogous to waves experiencing Faraday rotation effect, which occurs in materials [46, 47] and magnetized plasmas [48].

3.2 Electromagnetic waves propagating in a slowly rotating Gödel universe background

Let us consider a Gödel universe, with the stationary metric in cartesian coordinates [35, 49]

$$\begin{aligned} ds^2= & {} -dt^2+dx_1^2+dx_3^2\nonumber \\&+2\sqrt{2}\left[ 1-e^{\sqrt{2}x_1\Omega }\right] dx_2dt\nonumber \\&+\left[ 4e^{\sqrt{2}x_1\Omega }-e^{2\sqrt{2}x_1\Omega }-2 \right] dx_2^2\, , \end{aligned}$$
(18)

where \(\Omega \) is a constant related to the angular velocity of the rotating universe.

In the case of a slowly rotating universe, we have \(x_1\Omega \ll 1\), and from (5) we simply find

$$\begin{aligned} \epsilon ^{ij}\approx \delta ^{ij}\, , \quad {\mu _j}=-2\, x_1 \Omega \, \delta _{2j}\, , \end{aligned}$$
(19)

Using this in Eq. (10), it is straightforward to show that electromagnetic waves propagates satisfying the following equation

$$\begin{aligned} \mp i \, \omega \, \Psi ^j_{\pm }=-i\, \varepsilon ^{0jkm}\partial _j \Psi _{m\pm }+ 2\omega \,\Omega \, x_1 \varepsilon ^{0jkm} \delta _{2k} \Psi _{m\pm } \, , \end{aligned}$$
(20)

with the transformation \(F_{j\pm }(t, \mathbf{x})=\Psi _{j\pm }(\mathbf{x}) e^{-i\omega t}\) and wave frequency \(\omega \).

Equations (20) is exactly the same as Eq. (14) (replacing \(\eta \) by \(\Omega \)). Hence, we can find the same solutions, satisfying the dispersion relation (15), and wavelength difference (17). Anew, the polarized behavior of light is a consequence of its coupling with the spacetime of Gödel universe, whose rotation is a consequence of the total vorticity of the flow that generates it [49].

3.3 Electromagnetic waves propagating in anisotropic cosmology

Let us consider an anisotropic Bianchi–I cosmological model [50], with an interval given by

$$\begin{aligned} ds^2=-dt^2+a_i^2 dx^i dx^i\, , \end{aligned}$$
(21)

in cartesian coordinates (\(i=1,2,3\)). Here, \(a_i=a_i(t)\) are arbitrary time dependent functions, such as \(a_1\ne a_2\ne a_3 \ne a_1\), in general. With this metric, we find from Eq. (5) that \(\mu _j=0\) and \(\epsilon ^{ij}=\epsilon ^i\delta ^{ij}\), with

$$\begin{aligned} \epsilon ^i=\frac{a_1a_2a_3}{a_i^2}\, . \end{aligned}$$
(22)

This cosmology behaves in analogue fashion to a birefringent medium, as \(\epsilon ^1\ne \epsilon ^2\ne \epsilon ^3\ne \epsilon ^1\) in general [37]. Friedmann–Lemaître–Robertson–Walker cosmology is the particular case when \(a_1=a_2=a_3=\epsilon _1=\epsilon _2=\epsilon _3\).

Using the general metric (21), Eqs. (10) are now written as

$$\begin{aligned} \pm \, \partial _0\left( \epsilon ^i \delta ^{ij} \zeta _{j\pm }\right) =\varepsilon ^{0ijm}k_j \zeta _{m\pm }\, . \end{aligned}$$
(23)

where we have now considered that \(F_{i\pm }(t,\mathbf{x})=\zeta _{j\pm }(t) e^{i\mathbf{k}\cdot \mathbf{x}}\), with constant \(\mathbf{k}=(k_1,k_2,k_3)\). The three coupled Eqs. (23) can be solved for any component. As an example, let us obtain the wave equation for polarization \(\zeta _{3\pm }\). First, manipulating (23), the following condition is obtained

$$\begin{aligned} \zeta _{2\pm }=\frac{-k_2k_3 \epsilon ^3 \zeta _{3\pm }\pm k_1\partial _\xi (\epsilon ^3\zeta _{3\pm })}{\epsilon ^1 k_1^2+\epsilon ^2 k_2^2}\, , \end{aligned}$$
(24)

where the derivative \(\partial _\xi \) is in terms of the effective cosmological time

$$\begin{aligned} \xi =\int \frac{dt}{\epsilon ^1}\, . \end{aligned}$$
(25)

Using again (24) in Eqs. (23), we finally can get

$$\begin{aligned}&\partial _\xi \left[ \frac{\epsilon ^1 k_1^2}{\epsilon ^1 k_1^2+\epsilon ^2 k_2^2}\partial _\xi \left( \epsilon ^3\zeta _{3\pm }\right) \right] \nonumber \\&\qquad -\partial _\xi ^2\left( \epsilon ^3\zeta _{3\pm }\right) \mp \partial _\xi \left( \frac{\epsilon ^1 k_1k_2k_3}{\epsilon ^1 k_1^2+\epsilon ^2 k_2^2} \right) \epsilon ^3\zeta _{3\pm }\nonumber \\&\quad =\epsilon ^1 k_2^2\left( 1+ \frac{\epsilon ^3 k_3^2}{\epsilon ^1 k_1^2+\epsilon ^2 k_2^2}\right) \zeta _{3\pm }\, . \end{aligned}$$
(26)

From this equation can be readily seen that different polarizations propagate at different speed when \(k_1\ne 0\), \(k_2\ne 0\) and \(k_3\ne 0\). Equations similar to (26) can be found for \(\zeta _{1\pm }\) and \(\zeta _{2\pm }\).

In order to obtain a sensible and simple solution, let us consider the case of an almost–isotropic cosmology, such that \(a_2=a\), and \(a_1=a_3=a(1+\vartheta /2)\), with \(\vartheta \ll 1\). The anisotropy is in the direction of \(k_2\). In this case, we get \(\epsilon ^1=\epsilon ^3=a\) and \(\epsilon ^2\approx a(1+\vartheta )\). Furthermore, let us assume for simplicity the particular case when \(k_1=k_2=k_3=k\). Thus, Eq. (26) reduces to

$$\begin{aligned} \partial _\xi ^2\left( a\, \zeta _{3\pm }\right) +\frac{1}{2}\partial _\xi \vartheta \, \partial _\xi \left( a\, \zeta _{3\pm }\right) +\lambda ^\pm a\, \zeta _{3\pm }=0\, , \end{aligned}$$
(27)

where

$$\begin{aligned} \lambda ^\pm =\left( 3-\frac{\vartheta }{2}\right) k^2\mp \frac{k}{2}\partial _\xi \vartheta \, . \end{aligned}$$
(28)

The different behavior of polarizations is due to the difference of the variations of anisotropy parameters only. Eq. (27) can be solved for \(\zeta _{3\pm }\), as [51]

$$\begin{aligned} \zeta _{3\pm }(t)=\frac{\zeta _{3o}}{a}\exp \left( -\frac{\vartheta }{4}\right) \exp \left( -i\int \frac{W_{3}^{\pm }}{a}{dt}\right) \, , \end{aligned}$$
(29)

where \(\zeta _{3o}\) is a constant, and \(W_{3}^{\pm }\) fulfill the equations

$$\begin{aligned} \left( W_{3}^{\pm }\right) ^2= & {} \frac{3}{4}\left( \frac{\partial _\xi W_{3}^{\pm }}{W_{3}^{\pm }}\right) ^2-\frac{\partial _\xi ^2 W_{3}^{\pm }}{2\, W_{3}^{\pm }}\nonumber \\&+\lambda ^\pm -\frac{(\partial _\xi \vartheta )^2}{8} -\frac{\partial _\xi ^2\vartheta }{4}\, . \end{aligned}$$
(30)

Similar solutions can be found for \(\zeta _{1\pm }\) and \(\zeta _{2\pm }\).

The behavior of \(\zeta _{3\pm }\) depends on \(W_{3}^{\pm }\) through its evolution given in Eq. (30). This equation can be considered the dispersion relation for that particular propagation in this anisotropic cosmology. Hence, the wave frequency for each polarization can be defined as \(\omega _{3}^{\pm }= W_{3}^{\pm }/a\). This definition takes into account the cosmological redshift. An approximated solution can be obtained assuming that variations of \(W_{3}^{\pm }\) are negligible compared with variations of \(\vartheta \). In this case, we get

$$\begin{aligned} a\, \omega _{3}^{\pm }\approx \sqrt{\lambda ^\pm }\approx \sqrt{3k^2\mp \frac{k}{2}\partial _\xi \vartheta }\, . \end{aligned}$$
(31)

Wave frequency can be written as \(a\, \omega _{3}^{\pm }=\sqrt{3} k\pm \delta k\), as a deviation from a light–like behavior when \(k\gg \partial _\xi \vartheta \). Thereby

$$\begin{aligned} a\left( \omega _{3+}-\omega _{3-}\right)= & {} 2\delta k\nonumber \\= & {} \sqrt{3k^2- \frac{k}{2}\partial _\xi \vartheta }-\sqrt{3k^2+ \frac{k}{2}\partial _\xi \vartheta }\nonumber \\\approx & {} -\frac{a\, \partial _0\vartheta }{2\sqrt{3}}\, . \end{aligned}$$
(32)

Under the above assumption, it is clear that different polarizations propagate differently only due to their interaction with the anisotropic characteristic of this cosmology.

4 Casimir force due polarized electromagnetic waves

In Ref. [4] it was shown that polarized electromagnetic waves can give origin to repulsive or attractive Casimir forces when the wavevector of their left and right polarizations fulfill \(k^+\ne k^-\) in chiral media. This is exactly the behavior of differences (17) and (32) for polarized electromagnetic waves propagating in a Kerr and Gödel background metrics and in anisotropic cosmologies, respectively.

As a consequence, we can treat these curved spacetimes as analogue media for photon propagation with the sole purpose of calculating the associated Casimir force [33, 34]. In our case, polarized high–frequency electromagnetic waves fulfilling differences (17) and (32), modify the Casimir force decreasing its intensity, so it becomes less attractive.

4.1 Casimir force due to a rotating spacetime

In this case, the wavelenght difference is \(\delta k=-\eta \) for Kerr sacetimes, and \(\delta k=-\Omega \) for Gödel spacetimes. Now, let us consider two parallel plates with a length separation l, at some distance \(D\gg l\) far from the rotating black hole. Plates size are negligible compared to black hole dimensions. The plates are perpendicular to the rotation axis of the rotating spacetime, and in between them, polarized electromagnetic waves propagate. The Casimir force per unit area (2) experienced between two plates is then

$$\begin{aligned} {\mathcal{{F}}_C}= & {} \frac{1}{2\pi ^3}\int _0^\infty d\varpi \int _{-\infty }^\infty d^2 k_\parallel \nonumber \\&\times \left( \frac{\varrho e^{-4\varrho l}-\varrho e^{-2\varrho l}\cos (2\eta l)-\eta e^{-2\varrho l}\sin (2\eta l)}{1+e^{-4\varrho l}-2e^{-2\varrho l}\cos (2\eta l)}\right) \, ,\nonumber \\ \end{aligned}$$
(33)

where we have used \(\eta \) for both Kerr and Gödel cases. Although this result is general, we must notice that for Kerr metric \(\eta l=\left( {2J}/{D^2}\right) \left( {l}/{D}\right) \ll 1\), whereas for a Gödel universe \(\Omega l=(D\Omega )\left( {l}/{D}\right) \ll 1\). Using this, we can explicitly evaluate (2) to obtain

$$\begin{aligned} \mathcal{{F}}_C\approx -\left( 1-\frac{10}{\pi ^2}\eta ^2 l^2\right) \frac{\pi ^2}{240 l^4}\, , \end{aligned}$$
(34)

Casimir force per unit area (34) shows that spacetimes describing rotating bodies (such as a Kerr black hole and Gödel universe) induce a reduction on the magnitude of attractive Casimir forces. This effect is due to the chiral behavior of polarized electromagnetic waves triggered by the coupling of different circularly polarized waves and the rotation of spacetime, in an analogue fashion to what occurs in a chiral medium.

Strictly speaking, in order to obtain the result (34) we have used the flat–spacetime calculation for electromagnetic waves in a medium, developed in Ref. [4]. However, forces (33) and (34) are evaluated on curved spacetimes. This may be done because we have taken into consideration the effect of spacetime curvature into the electromagnetic wave propagation, by defining the effective medium described by the permittivity and permeability (12). Therefore, a medium treatment for Casimir energy becomes completely analogous to waves in curved spacetime. On the other hand, Casimir effect has already been studied in this limit on Kerr spacetimes for scalar fields [29,30,31, 52], and electromagnetic fields [31]. For both kind of fields, the gravitational correction to Casimir energy depends on \(\mathcal {O}(M/D)\), producing a reduction of the attractive force. The above used assumptions in our calculations require that \(M/D\ll 1\), as it can be seen in (12). This implies that we are neglecting the gravitational potential (Schwarzschild–like) correction to the Casimir force in our calculations. However, the net effect of the rotation spacetime does not vanish, as \(\eta \) is independent of the M. Previous results for the Casimir force produced by electromagnetic fields in a weak Kerr gravitational field have shown the M/D dependence without considering the polarization coupling with the rotation of the black hole [31]. As soon as this effect is taken into account, new corrections of order \(\eta ^2l^2\), not previously envisaged, emerge as a source of reduction of the magnitude of this force.

4.2 Casimir force due to anisotropic cosmology

In this case, the difference between left– and right–polarization is produced by anisotropies in spacetime, and is equal to \(\delta k=-a \partial _0\vartheta /(4\sqrt{3})\). This induces a modification on the Casimir force (2) for two parallel plates with a proper length separation \(a\, l\).

Consider a separation smaller than the scale of variation of the anisotropy, \(a\, l\, \partial _\xi \vartheta \ll 1\). In this case, force (2) for the \(\zeta _{3\pm }\) polarizations is simplified to

$$\begin{aligned} \mathcal{{F}}_C\approx -\left( 1-\frac{5 a^4}{24\pi ^2}\left( l\, {\partial _0\vartheta }\right) ^2\right) \frac{\pi ^2}{240 \, a^4 l^4}\, , \end{aligned}$$
(35)

showing how anisotropic cosmology variations reduces the strength of Casimir force due to electromagnetic interactions, as this expanding spacetime behaves as a medium. In the case of an isotropic cosmology with \(\vartheta =0\), Casimir force (35) reduces to its known \(\propto a^{-4}\) behavior [53].

5 Discussion

Results described in this work belong to a family of different behaviors of Casimir forces due to spacetime curvature. These have been studies in de Sitter [20], Schwarzschild [21,22,23], Gödel [54], wormhole [24] spacetimes, in cosmology [53, 55, 56], or in quantum cosmology [25], due to scalar [29, 30, 52] and vectorial fields [31], also showing that Casimir energy complies with the equivalence principle [26,27,28].

However, and differently, in this work we show that polarization of electromagnetic vacuum fluctuating fields interacting with curved spacetime (in a similar fashion to what occurs in the presence of materials) also modify the Casimir force. This effect manifests itself (in the calculations) only when electromagnetic waves are not treated as light rays (in the eikonal limit). It is the coupling between the spacetime curvature with the electromagnetic wave extended properties which modify the strength of Casimir forces, in addition to the curvature effects produced by the gravitational field. In this way, the spacetime background has a direct physical consequence that cannot be replicated by scalar fields or electromagnetic fields in the eikonal limit.

Finally, forces (34) and (35) hint that a more general spectrum of Casimir forces can be found for propagating electromagnetic waves in these kinds of curved backgrounds with inequivalent directions. The results presented in this article can be generalized to any spacetime with \(\mu _j\ne 0\) or non–constant \(\epsilon ^{jk}\). In general, for high–frequency electromagnetic waves in rotating spacetimes, the difference between the propagation of different polarizations should be proportional to \(|\mu |\), or variations of \(|\epsilon |\). Those effects create the small corrections to the Casimir force. However, one should expect to have larger effects if more general solutions of wave equation (10) can be found. Whenever the current assumptions can be relaxed, the effect induced by the different propagation for different polarizations of electromagnetic waves with longer wavelengths could be solved. Then it would be possible to study the repulsive effect on Casimir forces due to curved spacetimes, in an analogue fashion to what occurs in general chiral materials [4]. These endeavors are left as the subject for future research.