Maxwell equations in curved spacetime

In curved spacetime, Maxwell’s equations can be expressed in forms valid in Minkowski background, with the effect of the metric (gravity) appearing as effective polarizations and magnetizations. The electric and magnetic (EM) fields depend on the observer’s frame four-vector. We derive Maxwell’s equations valid in general curved spacetime using the fields defined in the normal frame, the coordinate frame, and two other non-covariant methods used in the literature. By analyzing the case in the generic frame we show that the EM fields, as well as the charge and current densities, defined in non-covariant ways do not correspond to physical ones measured by an observer. We show that modification of the homogeneous part is inevitable to any observer, and such a modification is difficult to interpret as the effective medium property. The normal frame is the relevant one to use as it gives the EM fields measured by an Eulerian observer.


I. INTRODUCTION
The special theory of relativity originated from Maxwell's electromagnetism.It showed that the electric and magnetic (EM) fields, E and B, are not independent concepts, and their nature depends on the observer's motion [1].Minkowski, by introducing a field strength tensor F ab , has constructed the spacetime covariant form of Maxwell's equations valid in the flat spacetime of special relativity [2].With the completion of the general theory of relativity, Einstein extended the spacetime covariant form of Maxwell's equations by Minkowski, now valid in curved spacetime [3].In a curved spacetime, besides their dependence on the motion of the observer, the EM fields also depend on the metric, the curved nature of the spacetime, through observer's four-vector.
Using a generic timelike four-vector u a , F ab can be decomposed into E a and B a in a spacetime covariant manner [4][5][6]; the charge and current densities, ̺ and j a , are similarly decomposed from the four-current J a based on the same four-vector.The four-vector u a depends on the observer.Often used ones are the fluid four-vector corresponding to the Lagrangian observer (comoving with the fluid), the normal four-vector n a corresponding to the Eulerian observer, the coordinate four-vector na corresponding to an observer attached with spatial coordinates, the coordinate observer; in the following, the normal four-vector, the normal frame, and the normal (or Eulerian) observer are used with the same meaning indicating n a .
In this work, we will consider the normal observer, the coordinate observer, and the all encompassing generic observer.As further two cases, we consider non-covariant definitions of the EM fields, which are historically interesting, partly coincide with flat spacetime results, and are often used in the literature related to the gravitational wave detection and the gravity-medium analogy.However, these two non-covariant definitions fail to identify the corresponding observer's four-vectors, and thus fail to introduce the external charge and current densi-ties properly [7].We will show this failure by analysing the generic frame; i.e., the non-covariant definitions do not have corresponding frame four-vectors.Without the observer's four-vector identified, the actual existence of such a set of EM fields as measurable quantities by any observer is in doubt.This also implies that the homogeneous part of Maxwell's equations is inevitably modified by gravity to any observer.Thus, the gravity-medium analogy is not valid.
In [7], we derived Maxwell's equations in the above four different definitions of the EM fields, assuming linear order perturbations in the Minkowski background metric, thus assuming weak gravity.Maxwell's equations were presented in the Minkowski background, with the effect of weak gravity appearing as the effective polarizations and magnetizations, Ps and Ms, often appearing in the homogeneous and/or inhomogeneous parts of Maxwell's equations.We showed why the non-covariant definitions of the EM fields are not only arbitrary but also wrong by missing the corresponding charge and current densities and missing the observer who can measure such fields.That is, the variables introduced in that manner are not EM fields in any sense.
Here we extend the case to the general spacetime metric.Equations will be presented using the EM fields and current density associated with the Euclidean threespace metric δ ij mixed with the Arnowitt-Deser-Misner (ADM) metric variables.The ADM metric variables can be expressed in terms of the fully nonlinear and exact (FNLE) perturbation metric variables where vector and tensor metric perturbation variables are associated with the Euclidean three-space metric δ ij and its inverse.We additionally present the case of the generic four-vector and its weak gravity limit.
In Sec.II we present the spacetime covariant formulation and the spatially covariant ADM formulation of Maxwell's equations.Section III presents Maxwell's equations in four different definitions of the EM fields.Section IV presents the equations in the generic frame and proves the inevitable nature of the gravity-modified homogeneous part.In Secs.V and VI we clarify a cou-ple of misconceptions in the literature and errors in the commonly used medium interpretation of gravity in electrodynamics.In Sec.VII, we discuss the result.Appendix A presents the ADM metric variables in terms of variables associated with δ ij and its inverse using the FNLE formulation.Appendix B presents relations to the linear order metric perturbation which will be useful to handle electrodynamics coupled with gravitational waves and the weak gravity.

II. COVARIANT AND ADM FORMULATIONS
In the presence of helical coupling, the electromagnetic part of Lagrangian is where the field strength tensor is φ is a scalar field with f = φ for the conventional axion coupling [8].A system together with Einstein's gravity, scalar field, and a general fluid is studied in [9].
Maxwell's equations in vacuum are The second equation is the same as F * ab ;b = 0. Using it, the first one becomes with this can be written as Using a generic time-like four-vector u a , with u a u a ≡ −1, we define the EM fields associated with the fourvector [4][5][6] We have E a ≡ F ab u b and B a ≡ F * ab u b with E a u a ≡ 0 ≡ B a u a .Thus, we have with The effect of axion can be regarded as the effective polarization and magnetization.The current four-vector is decomposed using the same four-vector as where ̺ and j a are charge and current densities, respectively, measured by an observer with the same fourvelocity u a .Using the decompositions, Eq. ( 3) and the second one in Eq. ( 2), respectively, give where g ≡ det(g ab ).As we have η abcd = ǫ abcd / √ −g with ǫ abcd an anti-symmetric symbol, Maxwell's equations can be written in forms with the effect of gravity (curved spacetime metric) behaving as the effective polarizations and magnetizations appearing in the homogeneous and inhomogeneous parts for a generic four-vector u a .In the following, we will show that modification in the homogeneous part due to gravity is inevitable to any observer even to the linear order perturbation.From now on, for clarity, we often indicate the spacetime covariant quantities using an overtilde.
We often use spatially covariant formulation of the curved spacetime.The ADM metric and its inverse are [10] where h ij is an inverse of the three-space intrinsic metric h ij , thus h ik h jk ≡ δ i j , and the index of N i is raised and lowered by h ij and its inverse; i, j . . .are spatial indices.We have where h ≡ det(h ij ), and indices of η ijk and η ijk are raised and lowered using h ij and δ ij , respectively, and their respective inverses.The normal four-vector is introduced as For the EM fields and current density in the normal frame, we set and similarly for E a , D a H a , and j a with E i ≡ E i etc.; indices of E i , etc., are raised and lowered by h ij and its inverse.Based on the normal frame, from Eq. ( 5) we have and similarly for H ab and H * ab with (D i , H i ) replacing (E i , B i ).
Maxwell's equations in the ADM formulation are presented in textbooks of numerical relativity [11][12][13]; see Eqs. (65)-(73) in [9].Equations can be arranged, including the axion coupling, as Using the vector notation, we have where the dot and cross products are associated with the intrinsic metric h ij , and ∇ is a covariant derivative associated with h ij .

III. MAXWELL'S EQUATIONS
In this section, we present exact forms of Maxwell's equations in the general curved spacetime using the EM field vectors associated with δ ij as the metric.We present four different forms based on four different ways of defining the EM fields.Two covariant methods are based on the normal frame and the coordinate frame as the observer's four-vectors.The other two non-covariant methods are based on the coincidences with the special relativistic equations noticed by Einstein in his way to guess the spacetime covariant form of Maxwell's equation valid in general relativity [3]; these non-covariant methods are also popularly used in the literature.
The two non-covariant definitions are not based on the observer's four-vector: the EM fields are directly read as components of special relativistic forms of F ab and F * ab , both with covariant indices.It turns out that, in this way, one cannot identify the frame four-vectors allowing such definitions of the EM fields, see Sec.IV.The situations are troublesome because without the observer's four-vectors specified, one cannot introduce the charge and current densities associated with the respective EM fields.Without the associated four-vector identified, the EM variables introduced in non-covariant manner are not the EM fields that can be measured by any observer.

A. Normal observer
The normal four-vector n a is normal to the hypersurface and is the four-velocity of an observer instantaneously at rest in the chosen time slice.It can be interpreted as an Eulerian observer as its motion follows the hypersurface independently of the coordinates chosen [12,14,15].
We introduce or directly, where the index of B i is raised and lowered using δ ij and its inverse; similarly we introduce E i , D i , H i , and j i .Equations ( 5) and (8) give and similarly for H ab .Notice that we are using mixed notation associated with δ ij and h ij and their respective inverses: the indices of N i and h ij are associated with the intrinsic metric h ij , whereas indices of η ijk , E i , etc., are raised and lowered using δ ij and its inverse.In Appendix A we will express the ADM metric variables in terms of metric perturbation variables associated with δ ij and its inverse.In this way, all variables in the following can be expressed entirely using variables associated with δ ij .Using Eq. ( 26), Eq. ( 2) gives where ∇ i is an ordinary derivative associated with δ ij .Useful relations needed for the derivation can be found in [9].In terms of the effective polarization and magnetization vectors, these can be arranged as where the effective Ps and Ms caused by the metric are To the normal frame observer, the gravitational field appears as effective Ps and Ms in both the homogeneous and inhomogeneous parts of Maxwell's equations.

B. Coordinate observer
In the literature, we often find a frame choice with u i ≡ 0 [11,15,16].This corresponds to the coordinate frame where the observer is attached with the coordinate, see below Eq. ( 63).In the coordinate frame, we have and the observer is at rest in the spatial coordinate, thus is a coordinate observer.
For the EM fields and current density in the coordinate frame, we set and similarly for Ēa , Da , Ha , and ja with E i ≡ Ēi etc.; indices of Bi , etc., are raised and lowered using δ ij and its inverse.Equations ( 5) and (8) give and similarly for H ab .Maxwell's equations follow from Eq. ( 2) as Similarly as in Eqs. ( 31)-( 34), these equations can also be expressed using effective Ps and Ms appearing in both the homogeneous and inhomogeneous parts of Maxwell's equations.
Although we suuggest the normal frame as the one to use because it corresponds to the frame of an Eulerian observer, as both frames are introduced in covariant manner, any frame is fine mathematically.Relations among parameters (the EM fields and charge and current densities) between the two frames can be easily read by evaluating Eqs. ( 5) and ( 8) in both frames.These are and similarly between (D i , H i ) and ( Di , Hi ).The coordinate frame coincides with the normal frame for N i = 0.
C. Special relativistic F ab or F * ab In this subsection we ignore the axion contribution for clarity of the argument.In his way of arriving at Eq. ( 2) in general relativity, in [3] Einstein noticed a couple of exact coincidences with the special relativistic Maxwell's equations based on two different definitions of the EM fields.These are and where indices of Êi , Bi , Ȇi , and Bi are raised and lowered using δ ij and its inverse.Notice that these definitions of EM fields, directly associating the physical EM fields to components of tensors, are not covariant.The merit of these two definitions noticed by Einstein [3] is that, in these ways, the second part of Eq. (2) using Eq.(44) and Eq. ( 4) using Eq. ( 45), in the absence of the external current four-vector, naturally gives the homogeneous part of Maxwell's equations in exactly special relativistic form as This set is a combination of parts of two sets of Maxwell's equations associated with two different definitions of the EM fields in Eqs. ( 44) and ( 45).The more complicated complementary parts are ignored, see below.
Notice that for the external source terms, we kept the components of the four-current without decomposing it into the charge and current densities.The two noncovariant ways of defining EM fields in Eqs. ( 44) and (45) do not guarantee the presence of the corresponding observer's four-vectors: the normalized four-vector only has three degrees of freedom, whereas in Eq. (44) or Eq. ( 45) six degrees of freedom are used.It turns out that even to the linear order perturbations in the metric we cannot find the four-vectors for the EM fields satisfying these conditions, see [7].Without the observer's four-vector, we cannot introduce the charge and current densities associated with the observer, which is true even for the EM fields as well.For further trouble with these non-covariant ways of identifying physical quantities as tensor components, see [17].
Comparing Eq. ( 26) with Eqs. ( 44) and (45) we have relations to EM fields in the normal frame These can be inverted to give Equations ( 46)-(49) look deceivingly simple: in the absence of the four-current, and if we ignore (we cannot, of course) the difference between ( Êi , Bi ) and ( Ȇi , Bi ), these are exactly the same as Maxwell's equations in special relativity.If we express the two sets using consistent notations, the other parts of equations become quite complicated.From and its inverse relation, we can derive the relation between the two EM fields which Einstein called "pretty complex relationship" [3].Using these relations between the two definitions in Eq. ( 53), two complete sets of Maxwell's equations can be written individually.In terms of Êi and Bi , we have where PE and ME are defined in Eq. ( 53).The effective P and M appear only in the inhomogeneous part of Maxwell's equations, and in such a case the effect of gravity can be interpreted as the effective medium property.
In terms of B and Ȇ, however, Eqs. ( 58)-(57) give where PB and MB are defined in Eq. (53).The effective P and M now appear in the homogeneous part which is opposite to the previous choice.
The above two non-covariant definitions are widely used in the literature of handling the effect of gravitational waves [18][19][20][21], and interpreting the effect of gravity as a medium property [4,16,[22][23][24][25][26].In these noncovariant definitions one cannot identify the four-vectors of corresponding observers who may measure such EM fields.Consequently, one cannot introduce the corresponding charge and current densities.We will elaborate our criticism on previous literature in Secs.V and VI.
Einstein in [3] used these two separate coincidences with special relativity as the guide to suggest the covariant form of Maxwell's equations in (2).However, he suggested non of these two as the definition of the EM fields [3].

IV. GENERIC OBSERVER
In the literature of gravitational wave detection using electromagnetic means, a special choice of the EM fields, assuming F ab (with two covariant indices) independent of the metric perturbations, is popular [18][19][20][21].In this way, the homogeneous part of Maxwell's equations becomes independent of gravity with apparent merit: the perturbed metric (including gravitational waves) appears only as the effective polarization and magnetization in ordinary Maxwell's equations in flat spacetime, similarly as in the axion case.
In the following, using a generic observer we will show that in the curved spacetime, such an ad hoc choice of the EM fields in Eqs.(44) or (45) is not possible for any observer.For the gravitation, it is inevitable to have effective polarizations and magnetizations appearing in both the homogeneous and inhomogeneous parts of Maxwell's equations.Gravitational wave detection should take this complication compared with the axion into account; for correct equations in the normal frame, see [7].
The fluid four-vector of a generic observer is introduced as [9] with the Lorentz factor where indices of V i and V i are raised and lowered by h ij and its inverse; V i is the fluid velocity measured by the Eulerian observer with n a [12,14,15], and V i ≡ u i /u 0 = dx i /dx 0 is the coordinate velocity of the fluid measured by the coordinate observer [11,15].For V i = 0 we have the normal frame of an Eulerian observer, and for V i ≡ 0 we have the coordinate frame na .We introduce B i ≡ B i with the index of B i associated with h ij as the metric.We have and similarly for E i , D i , H i , and j i .
Equations ( 9) and (10), with a = 0 and i, respectively, give The left-hand sides of Eqs. ( 65)-( 68) can be arranged so that the effect of gravity can be interpreted as the effective polarizations and magnetizations appearing in both the homogeneous and inhomogeneous parts as where indices of P i E , etc., are associated with h ij .The effective Ps and Ms caused by the metric and the axion are The effect due to axion is included in D i and H i .
A. The non-covariant definitions do not lead to measurable EM fields We will show that in the curved spacetime, the noncovariant definitions of the EM fields in Eqs.(44) or (45) are not possible for any observer.As a consequence, we will show that gravity inevitably causes modifications in both the homogeneous and inhomogeneous equations.From Eqs. ( 5) and ( 6), we have where the choice of observer's frame corresponds to choosing V i .These conditions relate the six F ab to the six EM fields.We cannot achieve arbitrary relations between them using only three conditions available in V i , which turns out to be the case for achieving F ab (or whatever combinations including F * ab ) independent of the metric perturbations even to the linear order [7].
To the linear order in metric perturbations and V i , we have Thus, even to the linear order in metric perturbations, one cannot remove the metric dependence of F ab by choosing V i .This conclusion is independent of the gauge conditions.For example, for the normal frame with V i ≡ 0, synchronous gauge condition with N = 1 and N i = 0 still leaves dependence on √ h and h ij through η ijk and B i .
Therefore, the non-covariant definitions of the EM fields introduced in [4,20,21,26], naïvely assigning the EM fields to components of the field strength tensor as in Eq. ( 44) cannot be regarded as any measurable EM fields.That is, as these variables cannot be measured as EM fields by any observer with their own four-vector, these are disqualified as the EM field variables.We can similarly show the non-covariant definition in Eq. ( 45) is not allowed.

B. Using flat spacetime variables
In terms of variables associated with δ ij , we introduce B i ≡ B i with the index of B i associated with δ ij , thus B i = h ij B j , and similarly for E i , D i , H i , and j i .We also introduce V i ≡ v i /c with the index of v i associated with δ ij , thus V i = h ij v j /c.Equations ( 69)-( 73) become where the effective Ps and Ms caused by the metric and the axion are ) For v i ≡ 0 we recover Eqs. ( 31)-( 35) in the normal frame, and for V i ≡ N i /N we recover Eqs. ( 39)-(42) in the coordinate frame.

C. Linear order metric perturbations
To the linear order, using Eq.(B2), we have and with From Eq. ( 81) we can see that it is not possible to remove dependence on the metric perturbation (h ab ) by choosing a frame (v i ) [7].For v i = 0 and v i /c = h 0i we recover the equations in the normal frame and the coordinate frame, respectively [7].

V. MISCONCEPTIONS IN THE LITERATURE
There are two wrong arguments [20,21] used to support the absence of metric in the definition of EM fields using F ab with two covariant indices; equivalently, the absence of metric in the homogeneous Maxwell's equations.
Authors of [20] argued for metric independence of the homogeneous Maxwell's equation using the topological nature of the equation in terms of form algebra. On this regards, however, the inhomogeneous equation is also topological, with no relation to metric.According to [27], "Remarkably, neither equation makes any reference whatsoever to metric.[T]he concepts of form and exterior derivative are metric-free.Metric made an appearance only in one place, in the concept of duality ("perpendicularity") that carried attention from F to the dual structure F * ."Thus, the authors of [20] confused the topological nature of the homogeneous Maxwell's equation as the absence of metric in the EM field decomposition.The metric appears when we related the two parts, and also in the decomposition of field strength tensor into the EM fields.
Authors of [21] similarly argued for metric independence of F ab using the absence of the metric in the relation between F ab and the four-potential A a , i.e., This, however, does not imply the absence of metric between F ab and the EM fields.In the normal-frame, Eq. ( 26) can be inverted to gives Using Eq. (87), we have Thus, metric is heavily involved in the relation between the EM fields and four-potential.Only in Minkowski spacetime, we recover the well-known relations E = ∇A 0 − A ,0 and B = ∇ × A. Thus, the authors of [21] confused the absence of metric in its relation to the four-potential as the absence of the metric in the EM field decomposition of F ab .The second one (homogeneous part) in Eq. ( 2) is naturally satisfied by using A a ; i.e., written in terms of the fourpotential, the homogeneous part is identically satisfied without metric influence.However, in terms of the EM fields, the homogeneous part is still coupled with the metric because the EM fields are non-trivially related to the four-potential for any choice of the frame four-vector; see Eq. (89) for the relation to the EM fields in the normal frame.As we show below, in curved spacetime the EM variables defined in the simple relations in Eq. ( 90) are nothing but Êi and Bi which do not correspond to measurable EM fields.
Authors of [19] derived perturbed Maxwell's equations via variation of the four-potential A a .As we mentioned, this corresponds to assuming F ab with two covariant indices independent of the metric which is fine if one uses A a .If one associate the potential to the EM fields naîvely using the special relativistic relations, the Maxwell's equations derived in this way correspond to the ones using Êa and Ba , i.e., Using Eq. ( 89), we can show that the relation between (E i , B i ) and ( Êi , Bi ) is the same as in Eq. ( 51).We have expounded in this work that Êi and Bi do not correspond to EM fields measured by any observer.Thus, the EM variables used in Maxwell's equations derived in such a way should be translated to the proper EM fields based on the observer's four-vector, for example, using Eq. ( 51) or (89) for the Eulerian observer.
The presence of metric in the decomposition of F ab and F * ab into the EM fields is partly related to the presence of metric in the four-vector in gravitational environments.The generic four-vector in Eq. (62) shows that for any choice of V i , it is not easy to have the special relativistic form of the normal four-vector n a = (−1, 0, 0, 0).In the normal frame, by imposing synchronous gauge condition with N ≡ 1 and N i ≡ 0, we can achieve this form with n a = (1, 0, 0, 0).But, even in such a case, Eq. (26) shows that F ij still depends on the metric through h kℓ and h.
The presence of metric in the decomposition of F ab into EM fields directly leads to the presence of metric in the homogeneous part of Maxwell's equations.Therefore, proposals for gravitational wave detection using the axion haloscope in [20,21] are based on wrong Maxwell's equations.Although the effect of axion can be treated as an effective medium, the effect of gravity is more complicated, as explained in this work.Proper study should use Maxwell's equations in the normal frame in Eqs. ( 27)- (30) or Eqs. ( 31)-(35).

VI. ERRORS IN MEDIUM INTERPRETATION
It is well known that the effect of axion coupling to electrodynamics in flat spacetime can be interpreted as an effective medium property in electrodynamics [8,28,29].There are numerous works interpreting the effect of gravity similarly as an effective medium property.To our knowledge, however, all the previous literature on the subject are, in fact, based on the two non-covariant identifications of the EM fields [4,16,[22][23][24][25][26], and thus are in error.
Most of the literature concerning medium interpretation of gravity follow Plebanski [24].Plebanski defined the EM fields in vacuum as where indices of E i , B i , D i , and H i are associated with δ ij and its inverse.Using Eq. ( 12), Eq. ( 92) is the same as Using Eqs. ( 91) and (92), Eq. ( 2), ignoring the axion coupling, gives These are ordinary Maxwell's equations in flat spacetime which Minkowski has derived in [2]; notice that we kept the four-current instead of the charge and current densities, see below.
In order to derive the effective constitutive relations, Plebanski used √ −gF ab = √ −gg ac g bd F cd .
Using Eq. (92), we can show with where ε ij and µ ij are the permittivity and permeability tensors, respectively, and γ ij is the electromagnetic mixing term; indices of ε ij , etc., are associated with δ ij and its inverse.Notice that definitions in Eq. (91) and Eq. ( 93) exactly correspond to the non-covariant definitions in Eqs. ( 44) and (45), respectively.Thus, (E i , B i ) and (D i , H i ) are nothing but our non-covariant definitions ( Êi , Bi ) and ( Di , Hi ), respectively.As explained in previous sections, in these non-covariant ways of defining the EM fields, one cannot identify the four-vector allowing these definitions and consequently cannot introduce the accompanying charge and current densities, even to the linear order in gravity [7]; this is why we kept the four-current in Eqs.(94) and (95).
The non-covariant definitions of EM fields introduced in [4,26] differ from the ones in Plebanski by using the three-space covariant metric γ ij , defined as γ ij ≡ g ij , instead of δ ij used in [24]; see Sec. 115 in [4] and Sec.90 in [26].In this way, we can show that the constitutive relations (thus the medium property) differ from Plebanski's [30].Using the non-covariant definitions of EM fields in gravity systems, the dependence of medium property on the spatial decomposition method and the coordinate used are also noticed in [31].We will analyse the errors in previous literature concerning the medium interpretation in detail in a subsequent work [30].
In Sec.IV we showed that the above non-covariant definitions do not have corresponding observer's four-vectors, thus these are not measurable EM fields by any observer.We also proved that any choice of the observer's fourvector leads to the presence of metric in the relation between the field strength tensor (in any form) and the EM fields.As a direct consequence, gravity causes modification of Maxwell's equations in both the homogeneous and inhomogeneous parts.The presence of metric in the homogeneous Maxwell's equations cannot be interpreted as the medium property in ordinary sense.In this way, the medium interpretation of gravity in previous literature is in error.

VII. DISCUSSION
Our main results are Maxwell's equations in a general curved spacetime presented in forms valid in Minkowski background with effects of gravity appearing as effective Ps and Ms.The ADM metric variables are expressed in terms of the metric perturbation variables associated with δ ij using the FNLE formulation, see Eq. (A2).We consider four different definitions of the EM fields, two based on the covariant decomposition, and two others using non-covariant identifications directly matching the EM fields with tensor components of F ab and F * ab .We also presented the case for the generic observer.
The two non-covariant definitions of EM fields lead to trouble as these fail to identify the observer's fourvectors, who may measure such fields, thus failing to identify the corresponding external charge and current densities, which also depend on the same four-vector (observer).This trouble was noticed even in the weak gravity case [7].Here, we prove that such definitions are not possible for the EM fields using the analysis for a generic observer.As a consequence, to any observer the gravity causes the homogeneous part of Maxwell's equations modified in addition to the inhomogeneous part.We propose the normal frame as the one to use as it corresponds to the Eulerian observer.Relations to the other definitions of the EM fields, including the non-covariant ones, are derived in this work.
In Laboratory situations, like in the Earth's gravita-tional field, in measuring passing gravitational waves using electromagnetic means, in inverse designing medium property for desired optical path using curved geometry, and in given weak gravitational fields, Maxwell's equations are enough with the gravity encoded in a given spacetime metric.The previous literature, however, rely on defining the EM fields in non-covariant manner.Here we showed that the EM fields defined in non-covariant manner are not the EM fields measured by any observer.Consequently, gravity causes modifications in both homogeneous and inhomogeneous parts, and the gravitymedium analogy is not valid.In the presence of gravity, for optical properties we should analyze Maxwell's equations directly.In the normal-frame, these are Eqs.( 27)- (30) or Eqs. ( 31)-(35) in general curved spacetime.
In realistic astrophysical situations, Maxwell's equations should be dynamically combined with Einstein's equation with accompanying fluids and other fields.The complete sets of Einstein's equation with a fluid and a scalar field (with additional helical coupling with the EM fields) in the covariant, the ADM, and the FNLE formulations are presented in [9].Maxwell's equations in the general curved spacetime derived in different forms in this work are applicable in all these situations.The cosmological FNLE formulation of Maxwell's equations presented in [9] considered the normal frame and ignored the tensor-type perturbations.Cosmological extension of the present work using the FNLE formulation is trivial.verse (C i ,i ≡ 0), and C ij is symmetric, transverse, and trace-free (C j i,j ≡ 0 ≡ C i i ).The perturbation variables may have arbitrary amplitudes only subject to the gravitational field equation.Without losing generality and mathematical convenience we can impose γ ≡ 0 ≡ C i as the spatial gauge condition; this is the unique spatial gauge condition which removes the spatial gauge degrees of freedom completely, thus the remaining variables (together with similar temporal gauge condition) are gauge invariant to fully nonlinear order in perturbations [32,34]; in this gauge condition we set B i = χ i .
The ADM metric quantities are derived in terms of the FNLE metric for the above general metric perturbations in [33].Imposing γ ≡ 0 ≡ C i as the spatial gauge conditions, we have For the general case without imposing the spatial gauge condition, see Sec. 2 of [33].Notice that we have not imposed the temporal gauge (hypersurface or slicing) condition; for several gauge conditions valid to fully nonlinear order, see [32][33][34].With these ADM metric quantities determined in terms of the metric variables associated with δ ij , Maxwell's equations in this work can be expressed as the ones in flat spacetime with the effect of gravity (metric) appearing as effective Ps and Ms.

Appendix B: Linear perturbation
To the linear order perturbations in the metric, using the notation in [7], i.e., g ab ≡ η ab + h ab with indices of h ab associated with η ab as the metric, we have These relations are valid even in the presence of γ and C i , thus for general linear metric perturbations without any condition.Using these relations between the FNLE metric and the linear perturbations, all the Maxwell's equations in this work reduce to the ones in [7] which are valid to linear order in metric perturbations.