Apr 15, 2010 - a guidance for future implementation of metamaterial invisibility cloaks. The geometry of the multilayer cylindrical cloak model is shown in Fig. 1, where a cylin- der of perfect electric conductor (PEC) with radius RM+1 is cloaked by
Apr 15, 2010 - Baile Zhangâ, and Bae-Ian Wu. Research Laboratory of ... practical development of metamaterial invisibility cloaking is still far from real application at this stage, since the rigorous ... by applying the genetic optimization ,
diffraction of noble gas trimers an outstanding enterprise. Firstly, diffraction presently is the only ...... coordinate R (not shown) corresponds to the center of mass position and is, therefore, identical for all three sets. ...... only on âÂ±1 a
Jan 26, 2018 - As a result, this method shares the simplicity of the mode- matching method in the ...... f1,II(q) = f3,II(q) = f5,II(q) = 3. 2 qeâiÏ/6+iÎ¸q , f2,II(q) = f4 ...
Aug 13, 2008 - trance to the gates of the temple of science are written the ..... tromagnetic waves in linear, homogeneous and isotropic media can be defined5.
Aug 9, 2006 - anisotropy Î± = Î¾â¥/Î¾ . Our simulations were typically carried out using relatively large system sizes (L = 512) and averaged over 100 runs. Fig. 1(a) shows typical results for the surface width as a function of film thickness for de
Jan 12, 2016 - using the standard approach. The distances travelled in each layer and the resulting weights are given in Appendix A, Table 2. In such a case, expression (2) becomes. cTI. 1212 = n. â i=1 di c2323i n. â j=1 dj. ,. (3) where di is t
Aug 13, 2008 - approach to electromagnetic scattering theory for cylinders, that is calculated analogue to Mie scattering  and to an exact solution of Maxwell's equations for cylindrical SPP. The validity .... roman era (4th century AD), which sh
Mar 7, 2010 - Bs, 42.25.Fx, 02.70.Bf, 73.20.Mf ... equations, show the numerical flow chart for matching the TF/SF boundary condition, and discuss several ..... 4(b)], the magnetic field pattern clearly reveals the interference between. 11 ...
Analysis of the evolution dynamics â consideration of simple cases . ..... extensively over the past 30 years and have been further advanced in recent years under the heading of complexity research, play an ..... new instruments into scientific pra
Katy BÃ¶rner, Cyberinfrastructure for Network Science Center, School of Library and Information Science, ... psychological perspective on individual scholars and groups; scientometrics, informetrics, and webometrics; the history of ...... http://garf
This paper describes a specific project, the current situation leading to it, its project design and first results. In particular, we ... The promise of the web-based LOD Cloud is to free up data, metadata and information to a ... and another website
Oct 27, 2011 - In this paper, key predictions of DSA at relativistic shocks that are germane to ..... This solution is accomplished outside the simulation program, and the jump conditions are then input manually as initial conditions for the simulati
Jul 23, 2009 - isting data on interferential mirror birefringence. We also report new measurements of two sets of .... The mirror reflectiv- ity can be inferred by the cavity finesse F = 6600 given by the authors at Î» = 633 nm: R ..... trend is cruc
Aug 2, 2011 - index s = 2 (Axford, Leer & Skadron 1977, Krimskii 1977,. Bell 1978 ... 2. A.R. Bell, K.M. Schure and B. Reville equation for the CR distribution function is. âf. ât. +(vx+u). âf. âx. â. âu. âx px. âf. âpx. â u c. â
Advanced Institute of Science and Technology (KAIST), Daejeon 305-701, Korea. 2. Center for Spintronics ... at an oblique angle (Î¸) to the spin direction (the oblique Hanle effect, OHE) in. CoFe/MgO/semiconductor (SC) ...... (Color online) (a)-(e) S
Jun 1, 2015 - discord state corresponds to an orthonormal basis, while a ... properly contained inside the set of zero-oblique-discord states, and the set of ...
Sep 7, 2007 - in-cell simulations in the case of oblique incidence of a laser pulse. It is shown that the ... high-Z atoms and a thin coating of low-Z atoms. Exten- sive computer ... box is removed to reveal the internal structure. (b,d): Dis- tribut
Keywords: Photonic crystal, Zero-n gap, Loss factor, Double-negative material. 1. Introduction. Photonic ... Email address: [email protected] (M. Barati) ..... Opt. Technol. Lett. 34 (2002) 409.  J. Li, L. Zhou, C.T. Chan, P. Sheng, Phys. Rev. Let
for the trace element analysis. The unique features of XRF technique make it a potential candidate for majority of applications in various fields of sciences.
This paper presents one aspect of a work-in-process PhD proposal that currently has the overarching objective of ... model design fault it appears not to have been extensively studied nor considered to be the source of .... single location in the mod
Jan 31, 2017 - their access to TB diagnosis and care is subject to their ability to access health care services and health insurance coverage, provided either by the state or the employer. Illegal migrants ..... Using the above model, we consider dif
May 4, 2009 - These measurements can constrain the star formation rate history (Hamann ..... Pittsburgh, University of Portsmouth, Princeton University, the ...
Scharnhorst effect at oblique incidence Stefano Liberati,1, ∗ Sebastiano Sonego,2, † and Matt Visser3, ‡ 1
arXiv:quant-ph/0010055v3 6 Mar 2001
International School for Advanced Studies, Via Beirut 2-4, 34014 Trieste, Italy 2 Universit` a di Udine, Via delle Scienze 208, 33100 Udine, Italy 3 Physics Department, Washington University, Saint Louis MO 63130-4899, USA (Dated: 16 October 2000; revised 13 December 2000; LATEX-ed June 2, 2018) We consider the Scharnhorst effect (anomalous photon propagation in the Casimir vacuum) at oblique incidence, calculating both photon speed and polarization states as functions of angle. The analysis is performed in the framework of nonlinear electrodynamics and we show that many features of the situation can be extracted solely on the basis of symmetry considerations. Although birefringence is common in nonlinear electrodynamics it is not universal; in particular we verify that the Casimir vacuum is not birefringent at any incidence angle. On the other hand, group velocity is typically not equal to phase velocity, though the distinction vanishes for special directions or if one is only working to second order in the fine structure constant. We obtain an “effective metric” that is subtly different from previous results. The disagreement is due to the way that “polarization sums” are implemented in the extant literature, and we demonstrate that a fully consistent polarization sum must be implemented via a bootstrap procedure using the effective metric one is attempting to define. Furthermore, in the case of birefringence, we show that the polarization sum technique is intrinsically an approximation. PACS numbers: 12.20.Ds, 41.20.Jb; quant-ph/0010055.
In 1990 Scharnhorst demonstrated that the propagation of light in the Casimir vacuum [1, 2, 3] is characterized by an anomalous speed . In fact photons propagating perpendicular to the plates travel at a speed c⊥ which slightly exceeds the usual speed of light c in the Minkowski vacuum [4, 5, 6, 7]. The propagation of photons parallel to the plates instead occurs at the usual speed ck = c. Unfortunately this anomalous propagation is far too small to be experimentally detectable; the relative modifications to the speed of light are of order 10−2 α2 /(me a)4 , where α is the fine structure constant, a is the distance between the plates, and me is the electron mass. [We work in units such that ¯h = c = 1. Greek and Latin indices run from 0 to 3 and from 1 to 3, respectively. ηµν = diag (−1, +1, +1, +1) denotes the Minkowski metric.] It is nevertheless an important point of principle that a quantum vacuum which is polarized by an external constraint can behave as a dispersive medium with refractive index n(ω) which remains less than unity to arbitrarily high frequency. Similar effects have been discovered in the case of quantum vacuum polarization in gravitational fields [8, 9, 10, 11], and several other situations . It is to be emphasized that the Scharnhorst effect, albeit small, is of fundamental theoretical importance.
Though the present calculations are carried out in the “soft photon” approximation (wavelengths much larger than an electron Compton wavelength) there is an argument based on dispersion relations which strongly suggests the actual signal speed is also modified (indeed, increased) [6, 7] — the physics here is very different from that of the resonance-induced “apparent” superluminal velocities currently of experimental interest [13, 14]; here we are dealing with quantum polarization induced “true” superluminal velocities, albeit well outside the realm of present day experimental technique. A common feature of the gravitational analogs of the Scharnhorst effect is the presence of birefringence: Photons with different polarizations propagate at different speeds. Similarly in generic situations nonlinear electrodynamics often leads to birefringence [15, 16, 17]. Nevertheless the occurrence of birefringence in nonlinear electrodynamics is not universal. Indeed, in the Casimir vacuum between parallel conducting plates the possibility of birefringence is tightly constrained: There is a residual (2+1) Lorentz invariance which prevents birefringence for photons that propagate parallel to the plates, and similarly there is an O(2) rotational invariance that prevents birefringence for photons that propagate perpendicular to the plates. It is only for photons that propagate obliquely to the plates that there is even any possibility of birefringence. We shall be interested in a complete analysis of both propagation speed and polarization states as a function of angle, and as a side effect report the (melancholy) conclusion that birefringence is completely absent in the Casimir vacuum. As a consequence, an effective metric description is sufficient for completely specifying photon propagation in the Casimir vacuum—though it will in general not be sufficient for other more general vacuum states. The effective metric we find is subtly at
2 variance with previous results obtained by performing a “polarization sum” [18, 19]. The disagreement is due to the way that “polarization sums” are implemented in the extant literature, where they are defined in terms of the flat Minkowski metric. We show that in the absence of birefringence a fully consistent polarization sum must be implemented in terms of a bootstrap procedure; using as input the effective metric one is attempting to define. Furthermore, in the presence of birefringence, we demonstrate that the polarization sum technique is intrinsically an approximation. Fortunately the differences first show up at order α4 , and in almost all circumstances are completely negligible. Throughout, we emphasize the use of symmetry arguments as a way of extracting general information that is (as much as possible) independent of the particular choice of Lagrangian. Finally, we have a few words to say concerning the utility of effective metric approaches in other contexts well beyond the Casimir vacuum. II.
Anomalous photon propagation can most easily be interpreted in terms of nonlinear electrodynamics. After integrating out virtual electron loops, the Maxwell Lagrangian should be replaced (in the absence of boundaries, at the one-loop level, and provided the distance scales defined by the field gradients and photon wavelength are much larger than the electron Compton wavelength) by Schwinger’s effective Lagrangian  LSchwinger = L (F , G) .
That is, the considerations of the following section are not limited to the Casimir parallel plate geometry but apply whenever the soft photon approximation (and the ancillary linearization procedure and restricted eikonal approximation) make sense. A.
Equations of motion
The complete equations of motion for nonlinear electrodynamics consist of the Bianchi identity, F[µν,λ] = 0, plus the dynamical equation ∂L = 0. ∂ν ∂Fµν
We now adopt a linearization procedure: Split the electromagnetic field into a background plus a propagating photon background photon Fµν = Fµν + fµν .
Then assuming the background satisfies the equations of motion, and retaining only linear terms in the propagating photon, we have photon f[µν,λ] = 0,
! ∂2L photon f = 0. ∂Fµν ∂Fαβ background αβ
The precise functional form of Schwinger’s Lagrangian will not be needed: At the end of the calculation we shall see that it is sufficient to retain only the quartic terms beyond the Maxwell Lagrangian (the Euler–Heisenberg Lagrangian ). Now in the Casimir geometry between parallel plates there is one additional invariant one could consider, namely ~ · ~n)2 − (~n × B) ~ 2, H ≡ nµ Fµσ F σν nν = (E
Leffective = L (Fµν (x), x) .
Here we have adopted the now common variables [18, 19] 1 1 ~ 2 ~ 2 F ≡ Fµν F µν = B −E , 4 2 1 ~ · B. ~ G ≡ Fµν ⋆ F µν = −E 4
of h ¯ /(me ca), but in the interests of generality we retain them for the time being and simply write:
where nµ ≡ (0, 0, 0, 1) is a unit vector orthogonal to the plates. (A similar invariant obtained by interchanging ~ and B ~ is actually a redundant linear combination of E the invariants F and H.) In addition the coefficients in the effective Lagrangian can now explicitly depend on the z coordinate. When it comes to practical calculations all these effects are suppressed by additional factors
∂2L , = ∂Fµν ∂Fαβ background
equation (10) can be rewritten in the somewhat more compact form photon = 0. (12) ∂α Ωµανβ fνβ
Note that the tensor Ωµναβ is symmetric with respect to exchange of the pairs of indices µν and αβ, and antisymmetric with respect to exchange of indices within each pair. We now apply a restricted form of the eikonal approximation by introducing a slowly varying amplitude fµν and a rapidly varying phase φ: photon fµν = fµν eiφ .
3 The wave vector (actually a one-form, but we shall loosely refer to both vectors and one-forms as “vectors” in the text, because of the mapping provided by ηµν and its inverse η µν ) is then defined as kµ = ∂µ φ. This approximation is similar to, but not quite identical with, the usual eikonal approximation. This is because one assumes that φ varies on scales much smaller than those of the background, while, on the other hand, use of the Lagrangian (5) also implies that the components of k are much smaller than the values fixed by the electron mass (soft-photon regime). Under these hypotheses, Ωµανβ kα fνβ = 0.
But the background field is itself subject to quantum fluctuations, and to take this into account the coefficients of this equation are identified with the expectation value of the corresponding quantum operators in the background state |ψi: hψ|Ωµανβ |ψi kα fνβ = 0.
In taking this expectation value we are using the fact that the fluctuations in the background fields are determined by the geometry (in the specific case of the Casimir geometry, by the distance between the plates), so that in the spirit of the restricted eikonal approximation there is a separation of scales between the background fluctuations and the propagating photon. The Bianchi identity (9) constrains fµν to be of the form fµν = kµ aν − kν aµ ,
where we have introduced the gauge potential a for the propagating field. Inserting (16) into (15) we find hψ|Ω
|ψi kα kβ aν = 0.
In general, the tensorial quantity hψ|Ωµανβ |ψi can be decomposed into an isotropic part plus anisotropic contributions, that we group together into a term ∆µανβ with the same symmetries as Ωµανβ : hψ|Ωµανβ |ψi = d1 η µν η αβ − η µβ η αν + ∆µανβ , (18)
where d1 is a function that can in principle be computed directly from the effective Lagrangian. B.
Equation (17) represents a condition for a as a function of k — it constrains a to be an eigenvector, with zero eigenvalue, of the k-dependent matrix µν
A (k) = hψ|Ω
|ψi kα kβ .
Any non-zero solution corresponds to a physically possible field polarization, that can be identified by a unit
polarization vector ǫ (provided a is not a null vector — a possibility that can always be avoided by a suitable gauge choice). A necessary and sufficient condition for the eigenvalue problem Aµν aν = 0 to have non-zero solutions is det (Aµν ) = 0; however, this gives us no information at all. Indeed, any a parallel to k is always a non-zero solution, so the condition det (Aµν ) = 0 is actually an identity. On the other hand, a k k is merely an unphysical gauge mode that corresponds to fµν = 0 by (16), so we need to find other, physically meaningful, solutions of the eigenvalue problem. To this end, we exploit gauge invariance under a → a + λ k and fix a gauge, thus removing the spurious modes. It is particularly convenient to adopt the temporal gauge a0 = 0. Then we can define a polarization vector ǫµ ≡ aµ / (aν aν )1/2 , and the eigenvalue problem Aµν ǫν = 0 splits into the equation A0i ǫi = 0,
plus the reduced eigenvalue problem Aij ǫj = 0 . The latter admits a nontrivial solution only if det Aij = 0.
The condition (22) plays the same role as the Fresnel equation in crystal optics  — it is a scalar equation for k and thus gives the dispersion relation for light propagating in our “medium”. An explicit calculation (see Appendix B) gives (23) det Aij = ω 2 P4 (k),
where ω ≡ −k0 and P4 (k) is a homogeneous fourth-order polynomial in the variables ω and ki . This means that in the most general case there are four dispersion relations, corresponding to the four roots of the equation P4 (k) = 0.
But, if k = (−ω, ~k) is a root then so is −k = (ω, −~k); thus, by CPT invariance, only two of these dispersion relations are physically distinct. Different polarization states are represented by linearly independent solutions of the eigenvalue problem (21), under the condition (22). [Obviously, (20) cannot be independent of (21), since we know that det (Aµν ) ≡ 0.] Thus, the space of polarizations is at most twodimensional. Since (22) gives rise to two dispersion relations, the polarization states actually satisfy two (in general, different) eigenvalue equations, µν
A(r) ǫ (r) ν = 0,
where r = 1, 2 labels the dispersion relations and A(r) is obtained from Aµν by imposing the corresponding condition on k. Thus, in the general case, modes of the field
4 with different polarizations have different dispersion relations, hence propagate in different ways; this leads to the phenomenon of birefringence [15, 16, 17]. In some special cases (in particular, this behaviour is quite common in nonlinear electrodynamics) the polynomial P4 (k) factorizes into two quadratic forms, αβ µν (26) kα kβ , P4 (k) = γ(1) kµ kν γ(2) in which case we obtain two second-order dispersion relations: µν γ(1) kµ kν = 0
µν and γ(2) kµ kν = 0.
two polarization states propagate with the same dispersion relation (28), and there is no birefringence. Substituting (28) back into the propagation equation (31) we find another relationship typical of this case, γ µν kµ aν = 0.
Formally, the above equation looks like a gauge condition. This might seem puzzling, because nowhere in the present subsection have we fixed a gauge. In fact, (32) is a consequence of the dynamical equation Aµν aν = 0, when the “on-shell” condition (28) is satisfied, and it does not imply any gauge fixing. Finally, it is interesting to notice that there is now a self-consistency or “bootstrap” condition,
3 Λ γ µν = hψ|Ωµανβ |ψi gαβ .
We now want to consider the special situation in which µν µν . (That is, not only does (26) hold, but also γ(1) = γ(2) the fourth-order polynomial P4 (k) is a perfect square.) In this case one ends up with a single quadratic dispersion relation of the familiar form γ
kµ kν = 0,
where γ µν is some symmetric tensor. It should be clear from our previous discussion that a necessary condition for this to happen is the absence of birefringence. Remarkably, we see that the wave vector is now null with respect to a (unique) “effective inverse metric” γ µν . Therefore, the propagation of light can be described in terms of an effective geometry, defined by the metric tensor gµν obtained by inverting γ µν , such that γ µν gνρ = δ µ ρ . (Warning: Even when an effective metric is defined, we always raise and lower indices using the flat Minkowski metric ηµν . In particular, note that gµν 6= ηµρ ηνσ γ ρσ ; this justifies our use of different symbols, g and γ, for the effective metric and its inverse.) This situation implies that hψ|Ωµανβ |ψi must be algebraically constructible solely in terms of γ µν . In view of the symmetries of Ωµανβ we know, without need for detailed calculation, that it must be of the form (29) hψ|Ωµανβ |ψi = Λ γ µν γ αβ − γ µβ γ να for some function Λ. Conversely, if Ωµανβ is of the form (29), then the matrix (19) is Aµν = Λ γ µν γ αβ kα kβ − γ µα kα γ νβ kβ , (30)
and the photon propagation equation Aµν aν = 0 becomes γ αβ kα kβ γ µν aν − γ αβ aα kβ γ µν kν = 0. (31) This equation is obviously satisfied by the uninteresting gauge modes a k k, with no constraints on k. Solutions corresponding to a non-vanishing fµν exist only if the coefficient of γ µν aν is zero, i.e., if (28) holds. Thus, the
This can be re-stated directly in terms of the fundamental coefficients as 3 Λ γ µν = d1 η µν η αβ gαβ − d1 η µα gαβ η βν +∆µανβ gαβ .
We stress that these relations depend only on the assumed existence of a single unique effective metric gµν — they do not make any reference to other specifics of the quantum state. III.
Let us now consider a region of empty space delimited by two perfectly conducting parallel plates placed orthogonal to the z axis at positions z = z0 and z = z0 + a, in the quantum state |ψi → |Ci corresponding to the Casimir vacuum [1, 2, 3]. A.
Symmetries and effective metric
In the Casimir vacuum considerable information can be extracted by using only symmetry considerations, similarly to what Bryce DeWitt did for the stress-energymomentum tensor in reference . Always working in the soft photon approximation, the presence of a preferred direction and the symmetry of the configuration allow us to claim that the function d1 can only depend on the z coordinate, and to write ∆µα νβ in the form ∆µα νβ = d2 (z) δ µ ν nα nβ − δ µ β nα nν +δ α β nµ nν − δ α ν nµ nβ ,
where d2 (z) is another function. The calculation of d1 (z) and d2 (z) in this case is straightforward, using identities that follow from (18) and (35): d1 (z) =
1 hC |Ωµν µν − 2 Ωµν αν nµ nα | Ci ; 6
5 1 d2 (z) = − hC |Ωµν µν − 4 Ωµν αν nµ nα | Ci . 6
(21). We find that ~ǫ must be a solution of the problem ij A ǫj = 0, where
On defining the function
d2 (z) , d1 (z)
and the tensor γ µν = η µν + ξ(z) nµ nν , hC|Ωµανβ |Ci takes the particularly simple form hC|Ωµανβ |Ci = d1 (z) γ µν γ αβ − γ µβ γ να .
Thus, from the discussion in subsection II C, we deduce without further argument that there is no birefringence, and that γ µν given in (39) plays the role of an effective inverse metric for the propagation of light in the Casimir vacuum. The corresponding effective metric is gµν = ηµν −
ξ(z) nµ nν . 1 + ξ(z)
Thus the dispersion relation determined from (22) takes in this case the form (28), as expected. We note here that, although an equation of the form (28) has been already obtained by Dittrich and Gies , the expression for the effective metric underlying their dispersion relation is subtly at variance with our results (39) and (41). We shall comment about the reason for this discrepancy in Appendix C.
Equation (42) then becomes
sin θ ǫ1 + (1 + ξ) cos θ ǫ3 = 0,
so two independent polarizations are:
which also follows from (32) when ǫ0 = 0. A tedious but straightforward calculation gives 2 det Aij = − (1 + ξ) d13 ω 2 γ µν kµ kν .
= −γ il kl γ jm km .
This eigenvalue problem, however, is manifestly equivalent to the single equation (42), which is thus the only constraint that the polarization states must satisfy. Therefore, the space of polarization states is twodimensional, as expected. A basis for such a space can be easily constructed by considering generic propagation in the xz plane [not a restrictive hypothesis, because of O(2) invariance with respect to rotations around the z axis], described by the 3-vector (these are taken to be covariant components, index down): ~k = |~k| sin θ, 0, |~k| cos θ . (45)
We stress that for this derivation of the effective metric to make sense it is necessary that the spatial variation in ξ(z) be compatible with the restricted eikonal approximation: that is ξ(z) should vary slowly on the scale of the photon wavelength. This is certainly true for QED at lowest nontrivial order where we shall soon see that ξ is in fact independent of z. The conclusion about the absence of birefringence can also be obtained more explicitly. Choosing the temporal gauge, and following the steps outlined in the general analysis of subsection II B, equation (20) becomes γ ij ki ǫj = 0,
In order to explicitly find the polarization states, let us evaluate (19) using the expression for hC|Ωµανβ |Ci and (28), and then consider the reduced eigenvalue problem
~ǫ (1) = (0, 1, 0); 1 ~ǫ (2) = (1 + ξ) cos θ, 0, − sin θ , N (θ) N (θ) ≡ 1 + 2 ξ cos2 θ + ξ 2 cos2 θ
is a normalization coefficient. Note in particular that ~ǫ (2) is not perpendicular to ~k when viewed in terms of the Minkowski metric η µν , though they are perpendicular when viewed in terms of the effective metric γ µν . Furthermore we have chosen to make ~ǫ (2) a unit vector with respect to the Minkowski metric, not with respect to the effective metric. C.
Phase, signal, and group velocities
If ~k has the form (45), equation (28) becomes (the norms are with respect to the Euclidean spatial metric induced by ηµν ): ω 2 = |~k|2 + ξ |~k|2 cos2 θ.
The phase velocity is given by vphase (θ) =
1/2 ω = 1 + ξ cos2 θ , |~k|
and is independent of the polarization. This is again a µν consequence of the fact that the quantity γ kµ kν apij pears squared in det A , so equation (22) describes a degenerate fourth-degree surface in the space of the vectors ~k/ω. Hence, there is only one dispersion relation,
6 equation (49), independent of the polarization state, and only one phase speed for each value of |~k|. This confirms that birefringence does not take place in the Casimir vacuum. Note that if ξ > 0 and θ 6= π/2, we have vphase (θ) > 1. Since, in the limit ω → +∞, vphase equals the front velocity (signal velocity) , one is tempted to argue that the propagation is superluminal for all values of θ different from π/2. This conclusion can also be inferred directly from (28), which for ξ > 0 implies that k is timelike (with respect to the undisturbed Minkowski metric). Since k can also be interpreted as the four-vector orthogonal to a surface of discontinuity of the field [16, 17], it follows that such a surface is spacelike (with respect to the Minkowski metric), i.e., that electromagnetic signals travel “faster than light”. (This phrasing is standard but unfortunate, and is logically indefensible. It should always be interpreted in the sense “faster than light would have travelled in an undisturbed portion of normal vacuum”.) Equivalently, (41) implies that for ξ > 0 the lightcones of the effective metric gµν are wider than those of the Minkowski metric ηµν . Unfortunately, while it is certainly true that our treatment is valid at high frequencies (with respect to those associated to the background scales), it nevertheless also requires ω ≪ me — the condition under which one can use the Lagrangian (1) —, so we have no direct information about the strict ω → +∞ limit. See, however, reference  for an indirect argument based on the Kramers–Kronig dispersion relation that combined with the present calculation is sufficient to establish superluminal propagation for the signal velocity. In the special case θ = π/2 one has η µν kµ kν = 0: Photons propagating parallel to the plates travel at the standard speed of light. The group velocity is a little tricky, it equals i vgroup ≡
∂ω 1 i k + ξ (k · n) ni , = ∂ki ω
and is not quite the same as the direction of propagation of the wave packet. Fortunately at both normal incidence (θ = 0) and parallel propagation (θ = π/2) the distinction between group and phase velocities disappears. Furthermore, as we shall see below, if we work to second order in the fine structure constant, the difference between group and phase speed is negligible at all angles, although a difference of order α2 still remains between the directions of ~vgroup and ~k. The group velocity (53) is greater than 1 when ξ > 0. Thus, the group speed is larger than 1 whenever the phase speed is, although in general their individual values are different.
FIG. 1: Group and phase speeds as a function of angle, from θ = 0 to θ = π/2 radians. (The group speed is the upper curve.) For clarity we have greatly exaggerated the physically expected value of ξ by setting ξ = 1/2.
so equations (45) and (50) give ~vgroup (θ) =
1 vphase (θ)
sin θ, 0, (1 + ξ) cos θ .
In particular, the group velocity ~vgroup is not parallel to the wave-vector ~k, though it is always orthogonal to the polarization vector. Taking the norm (in the physical Minkowski metric), vgroup (θ) =
The group speed is at all angles slightly greater than (or at worst equal to) the phase speed. (See figure 1.) Indeed 2 2 vgroup (θ) = vphase (θ) +
ξ 2 cos2 θ sin2 θ . 2 vphase (θ)
(See figure 2.) Care should be taken to realize that here θ is the angle between the wave vector and the normal,
SIZE OF THE EFFECT
So far we have seen that, based solely on symmetry considerations, we can eliminate the possibility of birefringence in the Casimir vacuum and place strong constraints on the general features of photon propagation in this vacuum. Our results are in fact generic to any form of nonlinear electrodynamics subject to the boundary conditions appropriate for the Casimir vacuum, and are not specifically restricted to QED. Where QED is important is in determining the specific form of the functions d1 (z) and d2 (z), which determine ξ(z), and hence determine the size (but not the qualitative features) of the Scharnhorst effect. Unless one wants to perform calculations to orders higher than α2 , it is sufficient to consider the Euler– Heisenberg Lagrangian  which, in the F –G formalism
7 c2 G µν ǫ αβ 2 c1 c2 ⋆ µν ⋆ + F µν Fαβ + F Fαβ , 2 2 +
so in the Casimir vacuum
d1 (z) = − 0.03
1 + O(α2 ) 16π
while d2 (z) is of order α2 . Hence, to first order in α2 , we have ξ(z) = −16π d2 (z), and
vphase (θ) = vgroup (θ) = 1 − 8πd2 (z) cos2 θ.
FIG. 2: Difference between the squares of the group and phase speeds as a function of angle, from θ = 0 to θ = π/2 radians. For clarity we have greatly exaggerated the physically expected value of ξ by setting ξ = 1/2.
LEH = −
1 F + c1 F 2 + c2 G 2 , 4π
α , 90π 2 m4e
Though in principle the coefficient d2 (z) could depend on z, the position relative to the two plates, we shall see that in the specific case of the Casimir vacuum it is simply a position-independent number. To establish this, start with (37), and insert the specific form (57). Using the algebraic identity [18, 20] ⋆
F µν ⋆ Fαν = F µν Fαν − 2 F δ µ α ,
7α . 360π 2 m4e
The terms proportional to F 2 and G 2 of this Lagrangian are quartic in the field, and describe the low-energy limit of the box diagram in QED, when four photons couple to a single virtual electron loop. Thus, the Lagrangian (55) is only accurate to order α2 , and it is meaningless to retain higher order terms within this model. In particular, deviations from (3+1) Lorentz symmetry due to the presence of the plates (the plates reduce the symmetry group to that of (2+1) Lorentz invariance) must on physical grounds vanish as the plate separation goes to infinity, where one must recover the full (3+1) Lorentz symmetry. On dimensional grounds this implies that such terms will be suppressed by some function of ¯h/(me ca) and can, to lowest nontrivial order, be neglected in the effective Lagrangian. On the other hand, even for lowest nontrivial order in the effective Lagrangian (that is, for the Euler–Heisenberg Lagrangian), the presence of the plates leads to a nontrivial expectation value for hC|Fµν Fαβ |Ci and will in this way contribute to the effective metric. For the Euler–Heisenberg Lagrangian, the tensor Ωµν αβ is 1 c1 F µν Ω αβ = − δµα δν β − δµ β δν α + 16π 2
one easily gets Ωµν αν =
adopted above, takes the form
1 c1 F − 3 δµα + 16π 2 (c1 − c2 ) F µ +2π (c1 + c2 ) T µ α + δ α , (61) 2
where T µ α is Maxwell’s stress-energy-momentum tensor. Performing the indicated traces, and using the fact that the Maxwell stress-energy tensor is traceless, we find d2 (z) =
4π (c1 + c2 ) hC |Tzz | Ci . 3
The symmetries of the Casimir vacuum stress-energy (as analyzed for instance by DeWitt ), then imply d2 (z) = 4π(c1 + c2 ) hC |T00 | Ci .
Finally the well-known result hC|T00 |Ci = −π 2 /720a4 allows us to write d2 = −
11πα2 . 64800 a4 m4e
And in particular, this implies d2 (z) is position independent. Thus, at first order in α2 , vphase (θ) = vgroup (θ) = 1 +
11π 2 α2 cos2 θ. 8100 a4 m4e
This expression reproduces Scharnhorst’s result in the case θ = 0, generalizing it to an arbitrary direction of propagation. It is interesting to notice that the correction is essentially determined by the expectation value of the energy density, hC|T00 |Ci, in agreement with the general results of [12, 18].
We have outlined a general scheme that allows one to write down the dispersion relation, and to find the polarization states, for electromagnetic radiation propagating in a region where nonlinear effects cannot be ignored. The deviations from the behaviour in Maxwell’s theory are completely described by the tensorial quantity ∆µανβ , that has its origin in the anisotropy of the medium. In particular, we have investigated the case of propagation in the Casimir vacuum, where there is a privileged direction in space, identified by a unit vector n. We have seen that symmetry considerations alone imply that ∆µα νβ must be of the form (35). This, in turn, implies that there is only one value for the speed of light, independent of the polarization. Thus, the possibility of birefringence in the Casimir vacuum is completely ruled out by very general arguments (the abstract form of the Lagrangian for nonlinear electrodynamics, plus the tensor structure of ∆µα νβ dictated by the geometry). Because of its generality, this conclusion applies not just to QED itself but to arbitrary types of nonlinear electrodynamics (such as for instance, Born–Infeld theories). How are we to interpret this result in view of the fact that nonlinear electrodynamics generically does lead to birefringence [15, 16, 17]? The key point is that in those analyses the background field (F µν )background is nonzero. More generally, even after averaging over quantum fluctuations of the background, hψ|F µν |ψi 6= 0 in the quantum state appropriate to those analyses. In the Casimir vacuum on the other hand, the expectation values linear in the field do vanish, hC|F µν |Ci = 0, and it is only the quadratic expectation values (and higher) that are nonzero (hC|Ωµναβ |Ci 6= 0). This key difference makes the analyses of [15, 16, 17] inapplicable to the Casimir vacuum—technically speaking, the polarization basis used in those papers fails to be meaningful for the Casimir vacuum. On the other hand, the discussion of the present paper could be easily adapted to treat light propagation in the presence of an electromagnetic background field. It is sufficient to keep in mind that in this case there is a preferred 2-dimensional plane defined by the 2-form hψ|F µν |ψi. Correspondingly, the form of ∆µανβ will be more complicated [see (A2)]. We reiterate that the absence of birefringence is crucial for the use of the “effective geometry” approach (adopted in this type of context by Latorre et al , and further developed by Novello et al , and by Dittrich and Gies ). Only if the propagation of light does not depend on its polarization and is thus, in a sense, universal, it is meaningful to describe it by a single effective metric. As a consequence of our analysis, photon propagation in the Casimir vacuum can indeed be phrased entirely in terms of the effective metric (39). This observation is potentially important in that “effective metric” approaches similar in spirit to the above are currently attracting attention in fields as diverse as acoustics [24, 25], optics [26, 27], superfluid quasiparticles , and Bose–
Einstein condensates . Finally, we mention related work on the light cone condition for a thermalized QED vacuum due to Gies  (wherein the analysis implicitly relies on taking the quantum expectation value of Ωµναβ in a manner somewhat analogous to the present paper), and intriguing results on photon propagation in rather general linear theories in classical backgrounds due to Obukhov and coworkers [31, 32, 33]. Acknowledgements
It is a pleasure to thank Holger Gies and Klaus Scharnhorst for several extremely useful comments that led to improvements in the presentation. MV was supported by the US Department of Energy. APPENDIX A: CONNECTION WITH OTHER FORMULATIONS
We have tried to set up the formalism in a streamlined and self-contained manner. Nevertheless to aid comparison with other results in the literature it is useful to give the explicit form of the tensor Ωµναβ for generic Schwinger-like Lagrangians of the form (1). For such Lagrangians 1 Ωµναβ = (∂F L) η µα η νβ − η µβ η να 4 1 1 + (∂G L) ǫµναβ + M µναβ , (A1) 4 4 where 2 M µναβ = F µν F αβ ∂F L + ⋆ F µν ⋆ F αβ ∂G2 L + F µν ⋆ F αβ + ⋆ F µν F αβ ∂F G L (A2) has the same symmetries as Ωµναβ . As soon as one inserts this tensor into the photon equation of motion (14), the completely antisymmetric part proportional to the Levi–Civita tensor drops out, because of the Bianchi identity (9). The remaining pieces reproduce the photon equation of motion in the perhaps more usual form considered by Dittrich and Gies , or Novello and co-workers [16, 17]. Also note that, depending on the details of the geometry and the quantum state, M µναβ can contribute to both d1 and to ∆µναβ . APPENDIX B: FRESNEL EQUATION IN NONLINEAR ELECTRODYNAMICS
Our purpose is to prove equation (23). The spatial components of the matrix (19) are Aij = ω 2 hψ|Ωi0j0 |ψi − ω km hψ|Ωi0jm + Ωimj0 |ψi +km kn hψ|Ωimjn |ψi. (B1)
9 Let us define the unit vector kˆ = ~k/|~k|. Then the components of Aij in a basis with one axis directed along kˆ are Aij kˆj = ω 2 hψ|Ωi0j0 |ψi kˆj −ω km hψ|Ωimj0 |ψi kˆj ≡ ω V i . (B2) In particular Aij kˆi kˆj = ω 2 hψ|Ωi0j0 |ψi kˆi kˆj ≡ ω 2 S. Then the matrix Aij has the following structure: 2 ω S ω VJ , ω V I T IJ
ˆ where I and J label the two directions orthogonal to k. Evaluating the determinant by expanding in the first row or column, it is easy to see that every term will contain at least two factors of ω, which establishes equation (23) as desired.
hψ|Ωµανβ |ψi kα kβ e ǫν(r) = 0.
µ ν = 0, d1 k 2 e ǫ(r) + ∆µα νβ kα k β e ǫ(r)
= 0 this can be written
µ where e ǫ(r) = η µν e ǫν . One then takes the scalar product of each equation with the corresponding polarization vector. Finally, the resulting two equations are summed. Furthermore, it was explicitly asserted that one could effectively replace  2 X r=1
µ ν → η µν e ǫ (r) e ǫ(r)
which is equivalent to equation (14) of reference . This has the form of the dispersion relation γ e µν kµ kν = 0, with e µν ∝ 2 d1 η µν + ∆µανβ ηαβ . γ
µ ν = η µν + uµ uν − wµ wν . e ǫ(r) e ǫ(r)
Working in the Lorentz gauge, as in reference , the two-dimensional subspace of Minkowski spacetime, orthogonal to e ǫ(1) and e ǫ(2) and spanned by the unit vectors u and w, contains the wave vector k, so k µ = − (k · u) uµ + (k · w) wµ .
Following the same steps that led to (C4), but using (C7) instead of the replacement (C3), we find, in place of the single term ∆µα να kµ k ν occurring in (C4), the two terms ∆µα να kµ k ν + k 2 ∆µανβ uµ uν wα wβ ,
where we have used (C8) and the symmetry properties of ∆µανβ . Clearly, the second term in the expression (C9) cannot be ignored unless k 2 = 0. We conclude that it is the subtly incorrect replacement (C3) which implies that the derivation of (C4), (C5), and (C6) is limited to first order in α2 . For the Casimir vacuum, ∆µανβ uµ uν wα wβ = d2 (n · u)2 − (n · w)2 # " 2 X 2 . n·e ǫ(r) = d2 −1 + r=1
in the terms containing ∆µα νβ . Under this hypothesis one obtains, following the steps described above, 2d1 k 2 + ∆µα να kµ k ν = 0,
2 d2 nµ nν , 2 d1 + d2
clearly different from (39). Nevertheless, when one considers QED corrections to the Maxwell Lagrangian, d1 is of order α0 , while d2 is of order α2 (see section IV above); therefore, the two metrics (39) and (C6) agree to order α2 . The reason for the discrepancy is that the formal replacement (C3) is limited in a subtle manner. Strictly speaking, it is correct as it stands only if photon propagation is governed by the Minkowski metric ηµν ; which is exactly the situation we are trying to get away from. In order to see this explicitly, let us introduce an orthonormal tetrad u, e ǫ(1) , e ǫ(2) , w , with u timelike and w spacelike, so one can write r=1
In this appendix we discuss the reason for the difference between our expressions (39) and (41) for the effective metric and those implicit in the “light cone condition” in reference . In the spirit of that paper a dispersion relation would be derived as follows (a certain amount of “translation” is required to get that formalism to match with the current notation). First, one writes equation (17) for two polarization states e ǫ (1) and e ǫ (2) orthogonal to each other: Working in Lorentz gauge k · e ǫ as
e µν = η µν + γ
APPENDIX C: POLARIZATION SUM
Specializing to the Casimir vacuum, and normalizing appropriately, one has
Using now the relation k 2 n · e ǫ(r) = 0, which follows from (C2) upon contraction with n in the Lorentz gauge, we see that k 2 ∆µανβ uµ uν wα wβ = −d2 k 2 ,
so (C4), corrected by the additional term (C10), indeed reproduces the metric (39). In general, if the photon dispersion is in fact determined by some (unique) effective metric gµν then one
10 can choose a tetrad like u, e ǫ(1) , e ǫ(2) , w , but which is orthonormal with respect to gµν , so now 2 X r=1
µ ν = γ µν + uµ uν − wµ wν . e ǫ(r) e ǫ(r)
Applying this polarization sum to the photon propagation equation written in the form (17), and using the properties (28) and (32), we see that the strictly correct replacement for the polarization sum is
[see (25)], with the same matrix A (k). However, if e ǫ (1) and e ǫ (2) correspond to two different dispersion relations, their propagation equations will actually differ, µν µν because now A(1) 6= A(2) . Therefore, a procedure of the form described at the beginning of this appendix is fundamentally flawed in the general case, as it fails to take into account the possibility of birefringence, and is in fact (strictly speaking) incompatible with birefringence. Indeed, summing the equations µν
2 X r=1
µ ν e ǫ(r) e ǫ(r) → γ µν .
This modified polarization sum, applied to (17), is consistent with the “bootstrap” condition (33). In the particular case of the Casimir vacuum, it is a simple matter to verify, using our expression for γ µν , that the bootstrap condition is indeed satisfied. It is now also clear why the formalism of  works to first order in α2 : Since the difference between η µν and γ µν is itself of order ξ (that is, of order α2 ) the difference between the (subtly incorrect) formula (C6), derived on the basis of (C3), and the consistency condition (33) is automatically of order ξ 2 (α4 ). This may also be explicitly verified using the approximation ξ ≪ 1 so that gµν ≈ ηµν − ξ nµ nν .
A(1) e ǫµ(1) e ǫν(1) = 0
A(2) e ǫµ(2) e ǫν(2) = 0, (C16)
which do take birefringence into account, appears rather impractical, and it is not obvious whether the structure of the resulting equation will allow the use of a polarization sum at all. Thus, producing “average dispersion relations” by polarization sums is at some deep level internally inconsistent, except in the cases when no polarization sum is needed. It is only if one is limiting interest to lowest-order (α2 ) corrections away from Minkowski space that the polarization sum makes sense for a birefringent system, and only then in the sense that 2 X r=1
µ ν e ǫ(r) e ǫ(r) → η µν + O(α2 ).
Substituting this into (33) self-consistently reproduces γ µν to the required order. Note that in almost all cases of interest one is quite content to work to first order in α2 in which case all of these subtleties are moot. We conclude with some general comments about the derivation of a metric from a polarization sum: If there is only one dispersion relation for all polarizations, as for instance in the Casimir vacuum discussed above, then e ǫ (1) and e ǫ (2) both satisfy the same equation (C15)
That is: Although equation (17) is exact, we could choose to contract it with approximate polarization states appropriate to flat Minkowski space. As long as the deviations from the ordinary speed of light are small, so are the deviations of the true polarization states from the usual ones, and then an approximate polarization sum of the above form is useful. We conclude that the polarization sum technique is most useful if there is no birefringence, and that in the presence of birefringence it is, at best, only useful when strictly limited to lowest-order corrections to Minkowski space propagation.
 H. B. Casimir, “On the attraction between two perfectly conducting plates,” Kon. Ned. Akad. Wetensch. Proc. 51, 793–795 (1948).  V. M. Mostepanenko and N. N. Trunov, The Casimir Effect and its Applications (Oxford Science Publications, Clarendon Press, Oxford, 1997).  G. Plunien, B. Muller and W. Greiner, “The Casimir effect,” Phys. Rep. 134, 87–193 (1986).  K. Scharnhorst, “On propagation of light in the vacuum between plates,” Phys. Lett. B236, 354–359 (1990).  G. Barton, “Faster-than-c light between parallel mirrors. The Scharnhorst effect rederived,” Phys. Lett. B237, 559–562 (1990).  G. Barton and K. Scharnhorst, “QED between parallel mirrors: light signals faster than c, or amplified by the
vacuum,” J. Phys. A 26, 2037–2046 (1993).  K. Scharnhorst, “The velocities of light in modified QED vacua,” Ann. Phys. (Leipzig) 7, 700–709 (1998) [hepth/9810221].  I. T. Drummond and S. J. Hathrell, “QED vacuum polarization in a background gravitational field and its effect on the velocity of photons,” Phys. Rev. D 22, 343–355 (1980).  R. D. Daniels and G. M. Shore, “ ‘Faster than light’ photons and charged black holes,” Nucl. Phys. B425, 634–650 (1994) [hep-th/9310114].  R. D. Daniels and G. M. Shore, “ ‘Faster than light’ photons and rotating black holes,” Phys. Lett. B367, 75–83 (1996) [gr-qc/9508048].  G. M. Shore, “ ‘Faster than light’ photons in gravitational
(k) e ǫ (r) ν = 0
  
   
fields — Causality, anomalies and horizons,” Nucl. Phys. B460, 379–394 (1996) [gr-qc/9504041]. J. I. Latorre, P. Pascual and R. Tarrach, “Speed of light in non-trivial vacua,” Nucl. Phys. B437, 60–82 (1995) [hep-th/9408016]. A. D. Jackson, A. Lande and B. Lautrup, “Apparent superluminal behavior,” physics/0009055. A. M. Steinberg, “No thing goes faster than light,” Phys. World 13 (9), 21 (2000). S. L. Adler, “Photon splitting and photon dispersion in a strong magnetic field,” Ann. Phys. (N.Y.) 67, 599–647 (1971); “Comment on “Photon splitting in strongly magnetized objects revisited”,” astro-ph/9601156. S. L. Adler and C. Schubert, “Photon splitting in a strong magnetic field: Recalculation and comparison with previous calculations,” Phys. Rev. Lett. 77, 1695–1698 (1996) [hep-th/9605035]. M. Novello, V. A. De Lorenci, J. M. Salim and R. Klippert, “Geometrical aspects of light propagation in nonlinear electrodynamics,” Phys. Rev. D 61, 045001 (2000) [gr-qc/9911085]. V. A. De Lorenci, R. Klippert, M. Novello and J. M. Salim, “Light propagation in non-linear electrodynamics,” Phys. Lett. B482, 134–140 (2000) [grqc/0005049]. W. Dittrich and H. Gies, “Light propagation in nontrivial QED vacua,” Phys. Rev. D 58, 025004 (1998) [hepph/9804375]. W. Dittrich and H. Gies, Probing the Quantum Vacuum, Springer Tracts in Modern Physics 166 (Springer, 2000). J. Schwinger, “On gauge invariance and vacuum polarization,” Phys. Rev. 82, 664–679 (1951). W. Heisenberg and H. Euler, “Consequences of Dirac’s theory of positrons,” Z. Phys. 98, 714–732 (1936). L. D. Landau, E. M. Lifshitz and L. P. Pitaevskii, Electrodynamics of Continuous Media, 2nd edition (Oxford, Pergamon Press, 1984). B. S. DeWitt, “Quantum gravity: the new synthesis,” in General Relativity, edited by S. W. Hawking and W. Israel (Cambridge, Cambridge University Press, 1979), pp. 680–745. W. G. Unruh, “Experimental black hole evaporation?” Phys. Rev. Lett. 46, 1351–1353 (1981); “Dumb holes and the effects of high frequencies on black hole evaporation,” Phys. Rev. D 51, 2827–2838 (1995), [gr-qc/9409008]. (Title changed in journal: “Sonic analog of black holes and...”) M. Visser, “Acoustic propagation in fluids: An unexpected example of Lorentzian geometry,” gr-qc/9311028; “Acoustic black holes: Horizons, ergospheres, and Hawking radiation,” Class. Quantum Grav. 15, 1767–1791
(1998) [gr-qc/9712010]; “Hawking radiation without black hole entropy,” Phys. Rev. Lett. 80, 3436–3439 (1998) [gr-qc/9712016]; “Acoustic black holes,” gr-qc/9901047. S. Liberati, S. Sonego and M. Visser, “Unexpectedly large surface gravities for acoustic horizons?” Class. Quantum Grav. 17, 2903–2923 (2000) [gr-qc/0003105]. U. Leonhardt and P. Piwnicki, “Optics of nonuniformly moving media,” Phys. Rev. A 60, 4301–4312 (1999) [physics/9906038]; “Relativistic effects of light in moving media with extremely low group velocity,” Phys. Rev. Lett. 84, 822-825 (2000) [cond-mat/9906332]. M. Visser, “Comment on Relativistic effects of light in moving media with extremely low group velocity”, Phys. Rev. Lett. 85, 5252 (2000) [gr-qc/0002011]. U. Leonhardt and P. Piwnicki, “Reply to the Comment on Relativistic Effects of Light in Moving Media with Extremely Low Group Velocity”, Phys. Rev. Lett. 85, 5253 (2000) [gr-qc/0003016]. G. E. Volovik, “Simulation of quantum field theory and gravity in superfluid 3 He,” Low Temp. Phys. (Kharkov) 24, 127–129 (1998) [cond-mat/9706172]. N. B. Kopnin and G. E. Volovik, “Critical velocity and event horizon in pair-correlated systems with “relativistic” fermionic quasiparticles,” Pisma Zh. Eksp. Teor. Fiz. 67, 124–129 (1998) [condmat/9712187]. G. E. Volovik, “Gravity of monopole and string and gravitational constant in 3 He-A,” Pisma Zh. Eksp. Teor. Fiz. 67, 666–671 (1998); JETP Lett. 67, 698–704 (1998) [cond-mat/9804078]. G. E. Volovik, “Links between gravity and dynamics of quantum liquids,” gr-qc/0004049. L. J. Garay, J. R. Anglin, J. I. Cirac and P. Zoller, “Sonic analog of gravitational black holes in Bose-Einstein condensates,” Phys. Rev. Lett. 85, 4643–4647 (2000) [grqc/0002015]; “Sonic black holes in dilute Bose-Einstein condensates,” Phys. Rev. A 63, 023611 (2001) [gr-qc/0005131]. H. Gies, “Light cone condition for a thermalized QED vacuum,” Phys. Rev. D 60, 105033 (1999) [hepph/9906303]. Y. N. Obukhov and F. W. Hehl, “Spacetime metric from linear electrodynamics,” Phys. Lett. B458, 466–470 (1999) [gr-qc/9904067]. F. W. Hehl, Y. N. Obukhov and G. F. Rubilar, “Spacetime metric from linear electrodynamics. II,” grqc/9911096. Y. N. Obukhov, T. Fukui and G. F. Rubilar, “Wave propagation in linear electrodynamics,” Phys. Rev. D 62, 044050 (2000) [gr-qc/0005018].