Arxiv:1603.08180V3 [Gr-Qc]

Home , Gravitoelectromagnetism

arXiv:1603.08180v3 [gr-qc] 1 Jun 2017 a.Bsdo rvoswrs(e 3 n references and [3] (see works Car- by previous way developed on general forms Based most of the calculus tan. in exterior the expressed using be can magnetism its and time (see and properties [2]). and space also in fields physical of questions with nature fundamental relation delves profound most the quest the namely a of physics, Such one into [1]. physics devel- directly important theoretical driving in uni- occurred, opments have a towards theory field attempts not fied is historical it many interactions, that Given these surprising between environments. analogies particular, and strong in bodies the astrophysical important, many electro- extremely in scales are these fields at magnetic even other but the dynamics, and On at astrophysical cosmological the governs also scales. substantially gravity relevant subatomic hand, extremely and although atomic is ex- scales, molecular, field physical macroscopic electromagnetic our at pervading the directly ones more the perience These are programs. fields rich research classical incredibly experimental an drive and to continue theoretical still to They first the studied. were be forces gravitational and electromagnetic ∗ † [email protected] [email protected] ti elkonta h hoydsrbn electro- describing theory the that known well is It the Nature, in interactions physical known the From lcrdnmc n pctm emty srpyia ap Astrophysical geometry: spacetime and Electrodynamics a neetcniee eecrepn setal othe to curvatur essentially spacetime correspond the here by considered induced interest are cal I phenomena spacetime. related curved in and electrodynamics of consequences cal ewrs lcrdnmc,(suo-imn emty a geometry, (pseudo)-Riemann electrodynamics, 04.40.-b Keywords: 04.30.-w, 04.20.Cv, numbers: PACS variability gravitational (radial) strong spatial for a objects of astrophysical obje possibility spherical rotating the around gra and of spacetime fields presence the the of in gravitomagnetism fields implicat electromagnetic the physical of the decay assumpti the discuss the with as: and both accordance solutions relax in simple we coordinate explore geometries radial symmetric the axially on pe the electric depend the vacuum fields where in Maxwell-Ampère the laws and conside between we Gauss relations alized geometries constitutive spherical the For in of tenso assumption priori. permeability properties a and magnetic isotropic permittivity and and the electric Therefore, Post-Newton the Parameterised tries. that the mas use hypothesis spherical we the latter, rotating the a In around regime. spacetime the and spacetime nttt eAto´sc ieca oEspa¸co, Faculdad Ciências do Astrof´ısica e de Instituto fe re eiwo h onain f(r-erc elec (pre-metric) of foundations the of review brief a After .INTRODUCTION I. dfı i 8 ap rne -7906Lso,Portugal. Lisbon, P-1749-016 Grande, Campo Edif´ıcio C8, rnic Cabral Francisco Dtd ue2 2017) 2, June (Dated: ∗ n rnic .N Lobo N. S. Francisco and ossol o ecniee ooeeu n isotropic and homogeneous ( considered ten- be permeability not Likewise, should magnetic sors and situations. permittivity physical electric specific the of geometry local led ttefudtoa ee hr sa niaelink intimate an is there level foundational the Indeed, at gravity. already cou- and the fields of electromagnetic context between the pling in experimentally and tested can be theories should electromagnetic extended or alternative smtisaenot are isometries spacetime (locally). the isometries follow assuming which elec- symmetries of (or have spacetime consequence trovacuum) of a properties electromagnetic as the isotropy, understood that and of be homogeneity idea can of the which requirement mentioned the briefly we go [10], letting example. work previous for a [9], In theory non-local [8] Mashhoon-Volterra are electrodynamics or changes non-linear Different in such Heisenberg-Euler freedom of from three. those greater arriving other a say, the models could than electromagnetic one postulate is, this it- There changing spacetime 7]. the [6, ex- for understood self by relations be constitutive supported can expressing vacuum well in as so relations not These is homogeneous periment. relations, local, isotropic linear, and considers postu- usually last the which while basis late, experimental good a have three we h ed ( fields be- relations the con- constitutive tween the flux and magnetic force ax- conservation, Lorentz or servation, charge postulates are main for- These four this on ioms. that established well clarified is [3–7] malism Obukhov and Hehl therein), priori a ned ntesii fgnrlrltvt G) such (GR), relativity general of spirit the in Indeed, tohsclapplications strophysical fields. eCˆnisd nvriaed Lisboa, de Universidade Ciências da de e ftevlct flgti aumaround vacuum in light of velocity the of ssol o ecniee homogeneous considered be not should rs osrltdt ieetpeoeasuch phenomena different to related ions eea,nweetoantccouplings electromagnetic new general, n olwn emtis h Schwarzschild the geometries: following iy antctrsi as a u to due law Gauss in terms magnetic vity, h ffc frlxn h homogeneity the relaxing of effect the r n xiain.Ti ffcstegener- the affects This excitations. and t,afaedagn ffc nelectric on effect frame-dragging a cts, ,rgrls fsaeielclsmere.Such symmetries. local spacetime of regardless ), rmgeim eepoesm physi- some explore we tromagnetism, a PN omls.W loexplore also We formalism. (PPN) ian n fhmgniyadiorp.We isotropy. and homogeneity of ons mtiiyadmgei permeability magnetic and rmittivity aumrflc h pctm isome- spacetime the reflect vacuum ntewa edadso rotation slow and field weak the in s .Teapiain fastrophysi- of applications The e. h oa smtiso pc.For space. of isometries local the E priori a , B † n xiain ( excitations and ) ie,isedte eetthe reflect they instead given, plications D , H .Tefirst The ). 2 between spacetime geometry, in particular, the conformal The applications in the following sections can have structure, and electrodynamics via the constitutive rela- some astrophysical relevance but the main purpose here tions [6, 7]. In this work, we start by considering linear, is to illustrate the effect of spacetime geometry in electric local, homogeneous and isotropic constitutive relations and magnetic fields motivated by relatively realistic as- and study the effect of spacetime geometry on electro- trophysical scenarios. An appropriate treatment of mag- magnetic fields for specific physical situations. With this netic fields in relativistic astrophysical situations (with restriction on the constitutive relations we are then re- strong gravity) would require the equations governing stricting the study to Maxwell electrodynamics in curved fluids and fields of magnetohydrodynamics in the back- spacetime around astrophysical objects. We will also ground of an appropriate spacetime metric. Furthermore, briefly study the effects of relaxing the condition of homo- for cases with very strong magnetic fields, which can con- geneity for spherical (isotropic) geometries and relaxing tribute to the gravitational field, the coupled Einstein- both homogeneity and isotropy for the axial symmetric Maxwell equations (or its generalizations) need also to case. The inhomogeneous equations, the Coulomb-Gauss be included in the analysis. Such procedures require ap- and Oersted-Ampère laws and the wave equations, con- propriate numerical methods and simulations. stitute the basis for the physical applications we explore The structure of this work is as follows. In Section in the following sections. The equations will be shown ex- II, we briefly review the foundations of electrodynamics plicitly and some physical consequences explored for two and its deep relation to spacetime geometry [10], pre- different geometrical backgrounds: 1) the Schwarzschild senting the basic equations to be explored in the applica- metric and 2) the geometry outside a stationary spher- tions. In Section III, we explore some physical implica- ical gravitational mass with slow rotation in the weak tions for two important geometrical backgrounds, namely field limit. (We considered the background of gravita- the Schwarzschild metric and the exterior geometry of a tional waves in a separate work [11]). The second case stationary spherical gravitational mass with slow rota- will be analysed using the Parameterized Post-Newtonian tion in the weak field limit. In Section IV, we discuss our (PPN) formalism, first developed for solar system tests, results and conclude. which allows to parameterize different theories of gravity (see [12–22]). As previously mentioned, both gravitational and elec- II. GENERAL FORMALISM tromagnetic fields are fundamental at astrophysical scales governing the dynamics and driving the complex pro- A. Pre-metric formalism plus spacetime cesses of structure formation (stars, galaxies, clusters, constitutive relations etc) as well as creating the conditions for extreme ther- modynamical states of matter and nuclear reactions. Before considering Maxwell fields in curved spacetime Many highly energetic phenomena in astrophysics involve let us recall the more general (premetric) formalism at very strong gravity and electromagnetic fields interacting the foundations of electrodynamics. In the pre-metric with very hot relativistic plasmas (see [23]). Therefore, formalism of electrodynamics, charge conservation gives the study of the coupling between gravity and electro- the inhomogeneous field equations while magnetic flux magnetic fields is fundamental to a deeper understand- conservation is expressed in the homogeneous equations ing of many observed phenomena in high energy astro- [3–7] physics. It is also relevant for the processes behind the formation and evolution of Active Galactic Nuclei (AGN) dJ =0 dG = J, ⇒ and to the physics of compact objects such as pulsars and dF =0 dA = F. (2.1) black holes. For example, it is well know that gravito- ⇒ magnetic fields play an important role in models for the These are fully general, coordinate free, covariant equa- collimation of astrophysical jets along a well established tions without the need of an affine or metric structure axis [24, 25], but such models should include the coupling of spacetime. Here J is the charge current density 3- of electromagnetic fields with gravity. Such couplings, in form, G is the 2-form representing the electromagnetic particular, with the gravitomagnetic term, might be use- excitation, F is the Faraday 2-form and A is the elec- ful to deepen the understanding of astrophysical jets in tromagnetic potential 1-form. The operator d stands for radio (active) galaxies. Magnetic fields are particularly exterior derivative. The equations above were presented relevant, since they pervade the physical Cosmos at plan- in a symmetric way to emphasize that, technically, the etary, stellar, galactic and extragalactic scales with differ- field G is a potential for the charge J in complete analogy ent magnitudes (see [26] and references therein). These with A, which is a potential for the field strength F. fields also play a vital role in the complex interaction What provides the metric structure are the constitu- between protostars and the environment of the hosting tive relations between G and F. We will assume linear, molecular giant clouds, which drives the processes of star local, homogeneous and isotropic relations given by the formation and constitute a fundamental ingredient for expression the understanding of (phenomenological) stellar formation rates in galaxies [27]. G = λ0 ⋆ F, (2.2) 3 where dimensional analysis allows us to write λ0 = 0 to 3 while Latin indices are space indices from 1 to ε0/µ0 with ε0 and µ0 being the electric permittivity 3. To explore electromagnetism in this geometry we will and the magnetic permeability of vacuum, respectively. use the tensor formalism in order to better identify spe- Thep Hodge star operator ⋆ introduces the metric and cific effects due to the components of the metric, for ex- maps k-forms to (n k)-forms, with n representing the ample, gravitomagnetic terms arising from off-diagonal dimension of the spacetime− manifold. With this assump- time-space components which are a characteristic of ax- tion, the inhomogeneous equations can then be written ially symmetric spacetimes. by Integration theory and geometrical methods in 3- −1 dimensional space can be used to derive the following d(⋆F)= λ0 J, (2.3) electromagnetic equations which emerge from the prin- or in terms of the potential by ciples of charge conservation, magnetic flux conservation and Lorentz force [3–7] d ⋆ dA = λ−1J. (2.4) 0 j ijk i i ∂j D = ̺, ǫ ∂j Hk ∂tD =  , In component form the above constitutive relations (2.2) − ∂ Bj =0, ǫijk∂ E + ∂ Bi =0, (2.7) are given by j j k t i i j k λ where ̺ √ gρ and  √ gj . Here D and B are = 0 √ ggαλgβγǫ F , (2.5) ≡ − ≡ − Gµν 2 − µνλγ αβ the components of vector densities (related to D and B k 1 kij k 1 kij βγ through the expressions D = ǫ Dij , B = ǫ Bij ). where g is the inverse metric and ǫµνλγ is the 4- 2 2 dimensional Levi-Civitta pseudo-tensor. The factor These equations also follow from (2.1) using (2.6). To αλ βγ √ gg g ǫµνλγ is conformally invariant which means solve these equations one needs to consider the constitu- that− at the foundational level electromagnetic fields are tive relations. It is easy to see [10] that the linear rela- coupled to the conformal part of the metric. As a re- tions in (2.5) or (2.2) imply that, in general, for the quan- sult, the light cone can be derived from linear and local tities in Eq. (2.7) we get the expressions D = D [E,B] electrodynamics [6]. The 2-form fields G and F can be and H = H [E,B], which include a mixing of electric and separated into their electric and magnetic parts accord- magnetic components and the influence of the spacetime ing to metric. As a consequence, magnetic terms appear in the Gauss law for geometries with non-vanishing time-space F = E dσ + B, G = H dσ + D, (2.6) metric components. This is easily seen by considering the ∧ − ∧ more traditional 4-dimensional tensor formalism. We will where E and H are the electric field and magnetic ex- consider this below. citation 1-forms, while B and D are the magnetic field The constitutive relations in (2.5) are implicitly as- and electric excitation 2-forms, respectively. Having per- sumed in the usual electromagnetic action of the more formed a spacetime foliation, these fields are defined on traditional variational or gauge approach. The set the 3-dimensional hypersurfaces orthogonal to a time di- of Maxwell’s equations in curved spacetime (pseudo- rection. Riemann geometry) are the well known expressions The field equations in (2.1) have a well established empirical foundation (charge and magnetic flux conser- µν µν 1 µν ν vations) while the constitutive relations are essentialy a µF = ∂µF + ∂µ(√ g)F = µ0j , ∇ √ g − postulate. Hehl and Obukhov have studied in detail more − general constitutive relations beyond the usual local and ∂[αFβγ] =0. (2.8) linear relations [7]. In [10], we suggested that the permit- µν αµ βν tivity and permeability tensors for vacuum should reflect where F = g g Fαβ and we have used in the inho- locally the spacetime isometries and should not be con- mogeneous equations the general expression for the di- sidered a priori homogeneous and isotropic. This hypoth- vergence of anti-symmetric tensors in pseudo-Riemann 1 esis changes the constitutive relations and therefore the geometry, Θµν = ∂ (√ gΘµν ). ∇µ √ g µ − dynamical equations for the electromagnetic fields. In − this work, we will also address these changes in specific With the definitions physical applications. With the constitutive relations the Ek i main axiomatic (foundational) formalism of electrody- F0k , Fjk Bjk = ǫijkB , (2.9) ≡ c ≡− − namics is complete. the homogeneous equations are the usual Faraday and magnetic Gauss laws B. Electrodynamics in curved spacetime i ijk j ∂tB + ǫ ∂j Ek =0, ∂j B =0, (2.10) We will consider a 4-dimensional spacetime manifold with pseudo-Riemannian geometry with a (+ ) sig- while the inhomogeneous equations can be separated into nature. Greek letters are spacetime indices ranging−−− from the generalized Coulomb-Gauss and Oersted-Ampère 4 laws. These are, respectively New electromagnetic effects induced by spacetime geometry include an inevitable spatial variability (non- αkj ∂ E + γj E cgmµgn0ǫ ∂ Bk k j j − kmn µ uniformity) of electric fields whenever we have non- mn0 k ρ vanishing geometric functions γk, electromagnetic os- cσ ǫkmnB = , (2.11) − ε0 cillations (therefore waves) induced by gravitational and radiation and also additional electric contributions to Maxwell’s displacement current in the generalized 1 µji 1 ji k mµ ni Maxwell-Ampère law. Notice that the functions ξii van- α ∂µEj + Ej ξ ǫkmn∂µB g g c c − ish for stationary spacetimes but might have an impor- k mni i B ǫkmnσ = µ0j , (2.12) tant contribution for strongly varying gravitational waves − (high frequencies), since they depend on the time deriva- with the following geometric coefficients tives of the metric. Moreover, besides these predictions, as previously said, for axially symmetric spacetimes grav- αkj g0kgj0 gjkg00 , αµji g0µgji gjµg0i , ≡ − ≡ − itomagnetic effects induce magnetic contributions to the Gauss law, with even static magnetic fields as possible sources of electric fields. These are physical, observable γj ∂ g0kgj0 gjkg00 ≡ k − effects of spacetime geometry in electromagnetic fields h 1 expressed in terms of the extended Gauss and Maxwell- + ∂ (√ g) g0kgj0 gjkg00 , √ g k − − Ampère laws which help the comparison with the usual − i inhomogeneous equations in Minkowski spacetime, making clearer the physical interpretations of such effects. mnβ mµ nβ 1 mµ nβ Finally, we review the field equations in terms of the σ ∂µ(g g )+ ∂µ(√ g)(g g ) , ≡ √ g − electromagnetic 4-potential which in vacuum are also use- − ful for the issue of electromagnetic waves. From (2.4) or

ji 0µ ji jµ 0i (2.8), we get ξ ∂µ g g g g ≡ − µ ν λν ε ν µ ν h 1 µ A g RελA ( µA )= µ0j , (2.17) + ∂ (√ g) g0µgji gjµg0i . ∇ ∇ − − ∇ ∇ √ g µ − − − where Rελ is the Ricci tensor. Using the expression i We will refer to these as simply the Gauss and for the (generalized) Laplacian in pseudo-Riemann mani- 1 Maxwell-Ampère laws. One sees clearly, that new elec- folds, µψ = ∂ √ ggµλ∂ ψ , and considering tromagnetic phenomena is expected due to the presence ∇µ∇ √ g µ − λ − of extra electromagnetic couplings induced by spacetime vacuum we arrive at curvature. In particular, the magnetic terms in the Gauss µ ν 1 µλ ν law are only present for non-vanishing off-diagonal time- ∂µ∂ A + ∂µ √ gg ∂λA √ g − space components g0j , which in linearized gravity cor- − λν ε ν µ respond to the gravitomagnetic potentials. These terms g RελA ( µA )=0, (2.18) − − ∇ ∇ are typical of axially symmetric geometries (see [28, 29]). For diagonal metrics, the inhomogeneous equations, which is a generalized Proca-like equation with vari- i.e., the Gauss and Maxwell-Ampère laws, can be recast able (spacetime dependent) effective mass induced by the into the familiar forms curved geometry. The second term in Eq. (2.18) can also be written in terms of the Levi-Civita connection, ̺ ∂ E˜k = , (2.13) αβ λ 1 αλ k through the formula g Γ = ∂α √ gg , ε0 αβ −√ g − − valid in pseudo-Riemann geometry. In usual Proca-like wave equations there is no such a term dependent on ǫ ∂ B¯iijjk = µ (i + i ), (2.14) ijk j 0 D the first derivative of the (massive) vector field. Simi- where lar terms appear for wave phenomena with longitudinal modes. For a diagonal metric in vacuum, we get j jj 00 E˜ g g √ gEj , ̺ √ gρ, (2.15) ≡− − ≡ − µ ν 1 µµ ν ∂µ∂ A + ∂µ √ gg ∂µA and √ g − νν − ε ν µ g Rεν A ( µA)=0, (2.19) B¯iijjk giigjj √ gBk, i √ gji, − − ∇ ∇ ≡ − ≡ − i 00 ii 2 ii with no contraction assumed in ν. In general, and con- D ε0√ g g g ∂tEi + c Eiξ . (2.16) ≡− − trary to electromagnetism in Minkowski spacetime, the (These equations are special cases of the more general equations for the components of the electromagnetic 4- expression which can be obtained from (2.8), namely potential are coupled even in the (generalized) Lorenz ∂ F˜µν = ν , where F˜µν √ gF µν and ν √ gjν .) gauge ( Aµ = 0). On the other hand, for Ricci-flat µ ≡ − ≡ − ∇µ 5 spacetimes, the term containing the Ricci tensor van- where ρm is the mass density. To solve this equation ishes. Naturally, the vacuum solutions of GR are ex- one needs the equation of state describing the thermody- amples of such cases. New electromagnetic phenomena namical properties of the interior of the star, as well as are expected to be measurable, for gravitational fields the other two remaining differential equations describing where the geometry dependent terms in Eq. (2.18) are the inner structure of spherical relativistic stars, which significant. are derived from Einstein’s equations. These correspond The inhomogeneous equations (Gauss and Maxwell- to the equations for ρm(r) and m(r). The coordinate Ampère laws) constitute the basis for the physical ap- pathology happening for r = 2GM/c2 in the black hole plications we explore in the following section. The equa- case, has no similar problem here because a careful anal- tions will be shown explicitly and some physical conse- ysis of the interior solution for ordinary stars shows that quences explored for different geometrical backgrounds: r> 2Gm(r)/c2 always [15]. A simple pedagogical model 1) Schwarzschild metric and 2) the geometry outside a which is not realistic (in fact it predicts an infinite speed stationary spherical gravitational mass with slow rotation of sound inside the star) is that of a star with constant in the weak field limit. The second case will be analysed density. Other more realistic well known exact solution using the PPN formalism (first developed for solar sys- is that of Buchdahl (1981) [15], which assumes an equa- 2 tem tests) which allows to parameterize different theories tion of state of the form ρm/c = 12√pp¯ 5p, wherep ¯ is of gravity. The case corresponding to the background of an arbitrary constant, which can be made− causal (with gravitational waves was presented in a separate work [11]. a speed of sound lower than that of light) and which for low pressures reduces to the equation of state of a n =1 polytrope in Newtonian theory of stellar structure [15]. III. APPLICATIONS Beyond ordinary stars, some old neutron stars can have a negligible rotation and have an approximately spheri- A. Electromagnetic fields in spherically symmetric cal gravitational metric field both inside and outside the spacetime matter. Nevertheless, the theoretical research on the equation of state for neutron stars is more delicate and 1. The spherically symmetric spacetime geometry it still has many open questions due to our lack of understanding of the properties of matter at supranuclear In this section, we consider the spherically symmet- densities characteristic of the inner cores of such compact ric spacetime around a spherical mass M, given by the objects. We hope to clarify many of the physical issues Schwarzschild solution of Einstein’s equations involved with the advent of precision gravitational wave astronomy, in particular with its applications to astero- ds2 = [1 Ψ(r)] c2dt2 [1 Ψ(r)]−1 dr2 − − − seismology. For simplicity, we will not consider interior r2 dθ2 + sin θ2dϕ2 , (3.1) solutions. − We now proceed with the analysis of electrodynamics where Ψ(r) 2GM/c2r. This is the metric represent- ≡ in the background of the spherical geometry in Eq. (3.1). ing the geometry outside a (non-rotating) spherical mass Let us start by observing that in general, the first term due to a star for example, or outside a spherical (non- in Maxwell’s inhomogeneous equations in Eq. (2.8) for rotating) blackhole. When we consider the above metric 2 any fixed ν, corresponds to the usual divergence, while in the next applications we assume that r > 2GM/c . the factors in the second term which are being contracted We recall that the coordinate r in general has no direct with (the contravariant components of) the Faraday ten- relation to the physical (proper) distance. Rather, it was 2 ′ 2 sor correspond to the components of a gradient. These chosen such that given the line element dl = f(t, r )dΩ terms must be computed using the appropriate expres- (where dΩ2 = dθ2 + sin θ2dϕ2) for the 2-spheres of con- ′ sions for a given system of coordinates. From these constant t and r , under the appropriate coordinate change siderations in spherical coordinates, with the metric in r′ r (such that f(t, r)= r2), the perimeter and area −→ Eq. (3.1), the following expressions are obtained for the of the 2-spheres with constant t and r, are given by the geometric functions that enter in the Gauss law (2.11) expressions 2πr and πr2, respectively [15]. Inside stars with radius R∗, imposing the continuity 4 2cot θ γr(x)= , γθ(x)= , γϕ(x)=0, for the metric functions at r = R∗, we get r r3 [1 Ψ(r)] − ds2 = e2Φ(r)c2dt2 [1 ψ(r)]−1 dr2 r2 dθ2 + sin θ2dϕ2 , (3.4) − − − (3.2) and in the Maxwell-Ampère law (2.12) 2 where ψ(r) 2Gm(r)/c r, m(R∗) M and the surface rϕϕ 1 rθθ 2 Ψ(r) of the star is≡ defined such that p(R≡)=0,(p is the pres- σ (x)= [2 Ψ(r)] , σ (x)= − , ∗ r3 sin2 θ − r3 sure inside the star). The potential Φ(r) is obtained from θrr 2cot θ θϕϕ the r component of the energy-momentum conservation σ (x)= 3 [1 Ψ(r)] , σ (x)=0, µν r − laws µT = 0, which gives [13, 15] ∇ σϕθθ(x)=0, σϕrr(x)=0, ξii(x)=0, dΦ dp (ρ c2 + p) = , (3.3) m dr − dr respectively. 6

2. Electrostatics in the Schwarzschild geometry Electric field due to charged plates in equatorial orbit in the (gravitational) weak field limit: The Gauss law (2.11) in the Schwarzschild geometry Previously, we assumed that the electric field had spher- (3.1), for the case of a static radial field with spherical ical symmetry just like the metric. Now we consider the symmetry becomes simplest case in which the electric field has a single component Ex (in a certain fixed system of coordinates). In dEr 4 ρ this case, outside the charge distribution, the Gauss law + Er = . (3.5) dr r ε0 (2.11) provides The charge distribution can be that of a spherical charge γx ∂ E = E . (3.9) or a spherical shell. In any case, a characteristic radius x x gxxg00 x can be defined, which we simply represent by R. In vacuum, outside the gravitational source, be it a star or a Let us consider the weak field limit of the Schwarzschild black hole (in which case we consider r > 2GM/c2) and solution (expanding the metric in powers of 2GM/c2r outside the charge distribution, the solution to Gauss’ up to first order). In isotropic Cartesian coordinates the law is metric is then given by C 2GM 2GM E = , (3.6) ds2 = c2 1 dt2 1+ r r4 − c2r˜ − c2r˜ × 2 2 2 −2 dx + dy + dz ,(3.10) contrasting with the simple r dependence in usual elec- × trostatics. Through dimensional analysis we can write where,r ˜ x2 + y2 + z2. We then have the following C = αkQR2, where, as mentioned before, R can be solution ≡ viewed as the radius of the charged body, k 1/4πε0 p ≡ and α is an arbitrary (dimensionless) constant which can x x ´x γ gxxg00dx r/R˜ sch 1 be fixed by comparing with the inner solution at r = R Ex = E0xe 0 = E0x − . (3.11) and imposing continuity. Near the surface of the mas- sr/R˜ sch +1 sive body the deviation of our solution with respect to In flat spacetime, the electric field is uniform with mag- the case without gravity is largest. In particular, such nitude given by E . Sufficiently far away from the grav- effects can be considered for charged stars or charged 0x itational source (star, planet or black hole) the electric black holes, by identifying R with R in the case of stars ∗ field is uniform. The effect of gravity on the electric field or with 2GM/c2 for the black hole case. can then be characterized by the dimensionless quantity One can also assume, as an interesting exercise, the case in which the charge distribution is due to a spherical r/R˜ sch 1 shell (with outer radius R) immersed in the Schwarzschild Ex E0x /E0x =1 − , (3.12) geometry (3.1). In such a case, inside the spherical | − | − sr/R˜ sch +1 charged body but still outside the black hole or outside which vanishes whenr ˜ . Recall that electric fields the star, by assuming uniform charge density, the solu- → ∞ tion of Gauss’ law is given by with quasi parallel field lines can result from charged plates. In particular, this result means that in the vicin- C˜ ρ ity of a spherical mass, the electric field created by a Er = 4 + r, (3.7) single charged plate or by a plane capacitor (oriented r 5ε0 along the x axis in this case) is no longer uniform. The 2 where 2GM/c Quasar, in the x body and increases more slowly inside, in comparison axis (of the Cartesian system centred in the gravitational with the same physical situation in flat spacetime. In source). Once the system is in orbit, that direction will particular in vacuum, we get be aligned with the radial direction relative to the central mass, twice per cycle. Alternatively this direction will be flat curved 4 Er (r) Er (r) R perpendicular to the radial line, also twice a cycle (when flat− =1 α 4 , r R. (3.8) Er (r) − r ≥ the system is crossing the y axis). The electric field, 7 which is always aligned with the x axis, has a magnitude be linearly polarized and is completely induced by the which varies with the distance to the center and therefore, effect of gravity in the electric field inside the capacitor twice a cycle the change in magnitude is either along the and the orbital motion of the spacecraft. line connecting the two plates or along the perpendicular to that direction. The maximum change in magnitude Electric field in the gravitational field of a mas- for both cases is equal if D = L. On the other hand, if sive spherical object – Case of non-homogeneous, D L and L is sufficiently small, then the field is ap- ≫ isotropic constitutive relations: In [10] we consid- proximately uniform whenever the system is crossing the ered the possibility of changing the constitutive relations y axis and non-uniform in the rest of the orbit. If, when- between the electromagnetic fields and excitations by re- ever the capacitor is crossing the x axis the maximum laxing the assumptions of homogeneity and isotropy. The change in magnitude is ∆Ex , then the dimensionless idea behind this suggestion comes from the very deep re- quantity of experimental| relevance| is then ∆E /E , | x| 0x lation between spacetime geometry and electrodynamics, which measures the strength of the effect of gravity in already present at the foundations of electromagnetism, 2 the weak field regime. The value E0x = Q/ε0L is the and so well explained by Hehl and Obukhov. But, it also usual value of the uniform field inside the capacitor in comes from the notion that the physical properties of the absence of gravity. vacuum (or electrovacuum) should not be a priori given. The strength of this effect when the capacitor is cross- This follows the spirit of GR which is a background inde- ing the x axis, with one plate at position x and the other pendent theory and therefore the local geometry of space- at x + D, is given by time is not a priori given, rather it has to be considered for each physical system as a solution of the dynamical Ex(x) Ex(x + D) 2 X˜ 1 equations. Likewise, we postulate that the electric per- | − | = 1 − , mittivity and magnetic permeability tensors for vacuum E0x s − 1+ X˜ + D˜ − sX˜ +1 (3.13) should reflect the local symmetries of spacetime geometry. where the distances are in units of R , X˜ r/R˜ , sch ≡ sch Let us consider spherical symmetry. Relaxing the as- D˜ D/R . This dimensionless quantity measures the ≡ sch sumption of homogeneity around spherical bodies, the maximum change of the magnitude of the electric field in- linear, local, isotropic and non-homogeneous constitutive side the capacitor due to the gravitational field of the as- relation between the electric field and electric excitation trophysical spherical mass. In principle, this effect could is [10] be used to test GR and modified theories of gravity, pro- viding another complementary (weak field) test in the i 00 ij D = g g √ gε0(r)Ej . (3.15) Solar System. − − Since the values of the Schwarzschild radius for the Sun The Gauss law can be written as and Jupiter are approximately 2.95 103m and 2.2m, × i respectively, the effect should be very small unless one ∂iD = ρ√ g. (3.16) gets extremely close to their surfaces. For example, for a − laboratory in orbit around Jupiter at a distance approxi- We will use the following Ansatz mately equal to three times Jupiter’s radius (X˜ 9 107) ∼ × 2GM we have the following values ε (r)= ε 1+¯γ . (3.17) 0 0 c2r ∆E /E 10−14 (D˜ = 100) , | x| 0x ∼ −13 which can be seen as the linear approximation of a Taylor ∆E /E 10 (D˜ = 1000). (3.14) 2 | x| 0x ∼ expansion in powers of 2GM/(c r). Hereγ ¯ is a dimen- Naturally, these tiny values represent an enormous ex- sionless parameter. Then, considering the Schwarzschild perimental challenge in terms of the sensitivities and geometry, and choosingγ ¯ = 1 we find the following solu- noise control requirements. The fact that we are not us- tion in vacuum ing test masses, but electric fields instead to study grav- C 1 E = . (3.18) ity, complicates further the experiment due to various r r4 (1 + R /r) possible environmental effects related to space weather, sch in particular, solar and planetary magnetospheres. In This result includes the effect of a radial dependence of any case, it is always better to measure voltage differ- the electric permittivity of vacuum for spherical gravita- ences than to measure the electric field itself, since better tional fields. The strength of the effect is given by sensitivities can be obtained. The Voltage drop between the two plates will be less than in the absence of gravity. ε0(r) ε0 ε0 Rsch/r Er Er /Er = , (3.19) In principle, since different points inside the capacitor | − | 1+ Rsch/r will be at different radial coordinates which change with ε0 time, we also expect the generation of electromagnetic ra- where Er corresponds to the solution in (3.6). diation with a frequency related to the orbital frequency According to the hypothesis we explore here in the of the spacecraft. This electromagnetic radiation should curved geometry around a massive object, the electric 8 permittivity of vacuum is not homogeneous, but has a ra- Magnetic field created by a ring of circulating dial dependence instead. For black holes, this effect influ- plasma around a spherical mass: Consider a spher- ences electric fields more strongly near the horizon. The ical astrophysical mass, which could be a black hole, with strength of the effect is maximum at the horizon where a stationary current loop around it on the equatorial the magnitude of the electric field is 0.5 times weaker plane (θ = π/2). This current distribution, which could than the case with a constant permittivity tensor. be due to a disk of very hot plasma, would create an If we now repeat the analysis for the case of the electric axially symmetric magnetic field of the general form field inside a plane capacitor in orbit around a spherical r θ mass, we arrive at the result B = B (r, θ)er + B (r, θ)eθ , (3.26) 1 R2 /r˜2 which implies E = E − sch . (3.20) r˜ 0x (1 R /r˜)2 p sch r θ − H = Hr(r, θ)e + Hθ(r, θ)e . (3.27) The strength of the effect can be quantified again by comparison with the case where the electric permittivity is Then, the relevant Maxwell equation is constant, where we get 1 1 ϕ ∂r(rHθ) ∂θHr = √ gj , (3.28) ε0(˜r) ε0 ε0 r/R˜ sch r − r − Er˜ Er˜ /Er˜ =1 . (3.21) | − | − 1+˜r/Rsch which can be solved outside the current distribution by In principle, this effect could be tested using electric fields setting the RHS to zero. We will consider this equation inside capacitors and these expressions can be generalized in the equatorial plane and outside the disk (or ring) of to include the PPN approach. currents, i.e., for r > R where R is a mean represen- tative of the radius of the circular current distribution. Setting θ = π/2 the magnetic field will only have the θ 3. Magnetostatics in the Schwarzschild geometry – component Astrophysical applications B = Bθ(r,π/2)e , (3.29) ± θ Suppose that in some reference frame we have a static where the refers to the cases in which the circulating magnetic field due to a (stationary) current and no elec- ± tric field. Then, the generalized Maxwell-Ampère law in current is moving in the direction of positive or nega- Eq. (2.7) is given by the following equations tive ϕ, respectively. We then get Hθ(r,π/2) = f(π/2)/r, which implies 1 1 r ∂θ (sin θHϕ) ∂ϕHθ = √ gj , (3.22) r sin θ − r sin θ − θ µ0 1 B (r,π/2) = Hθ 5 . (3.30) √ ggθθ ∝ r 1 1 − ∂ H ∂ (rH )= √ gjθ, (3.23) r sin θ ϕ r − r r ϕ − We conclude that the magnetic field (on the equatorial plane) due to the circular current distribution decays faster with the radial distance, in the curved spacetime 1 1 ϕ ∂r(rHθ) ∂θHr = √ gj . (3.24) of Schwarzschild geometry, than in the flat (Minkowski) r − r − case (although, strictly speaking the radial coordinate We can solve this for the homogeneous and isotropic here does not correspond to a physical distance). constitutive relations [10], that can be derived from Eq. We recall that at this stage we are neglecting rota- (2.5): tion, therefore these calculations can be viewed as having an approximate validity around quasi-static spherical H = µ−1√ gg Bj . (3.25) k 0 − jk masses with an electric current due to a highly ionized gas The following applications can have some astrophys- in the orbiting accretion disk. Since the electromagnetic ical relevance although, as previously mentioned, the equations were considered and solved outside the cur- main purpose here is to illustrate the effect of space- rent distribution creating the magnetic field and in the time geometry in magnetostatics motivated from mini- exterior (Schwarzschild) spacetime, another astrophysi- mally realistic astrophysical scenarios. An appropriate cal scenario compatible with the approach taken here is treatment of magnetic fields in relativistic astrophysical that of a (quasi) spherical mass with negligible rotation situations (with strong gravity) would require the equa- with a magnetic field generated by electric currents in tions governing fluids and fields of magnetohydrodynam- its interior, as long as the magnetic field energy density ics in the background of an appropriate spacetime metric. has a negligible effect in the gravitational field (no back Furthermore, for cases with very strong magnetic fields, reaction). which can contribute to the gravitational field, the coupled Einstein-Maxwell (or its generalizations) also need Magnetic field created by an astrophysical jet: to be included in the analysis. Consider now another idealized astrophysical scenario in 9 which the spherical body emits a stationary jet of charged the constitutive relations and the physical dimensions of particles vertically, defining an axis. We set this to be the the magnetic field. Indeed, the magnetic permeability z axis. In this case, a magnetic field will arise with axial appears in the denominator in (3.34) while the electric symmetry and along the ϕ direction with symmetrical permittivity appears in the numerator in (3.15). Like in configurations above or below the θ = π/2 plane the case of the electric field, sufficiently far way from the source, the effect becomes neglibible. B = Bϕ(r, θ)e . (3.31) ± ϕ The relevant Maxwell-Ampère equations outside the cur- B. Maxwell fields around slowly rotating objects in rent distribution are the weak field limit within the PPN formalism 1 1 ∂ (sin θH )=0, ∂ (rH )=0, (3.32) r sin θ θ ϕ − r r ϕ The appropriate metric describing the local geometry outside a stationary rotating mass is the Kerr metric [30]. and therefore H (r, θ) (r sin θ)−1, which implies ϕ ∝ Expanding this geometry to first order in the angular µ 1 momentum J, we get the geometry outside the source in Bθ(r, θ)= 0 H . (3.33) ϕ 5 4 the limit of slow rotation [13] √ ggϕϕ ∝ r sin θ − Here we considered a constant magnetic permeability 2 2 4GJ 2 dsKerr dsSchawrzschild + sin θdϕdt. (3.37) tensor for vacuum, but following our hypothesis that the ≈ c2r properties of this tensor should reflect the local space- Recall that the deformations of the object’s spherical time isometries it should be interesting to compute the symmetry depend quadratically with angular momentum changes to these results if the magnetic permeability has J, while the metric already changes at the linear level. a radial dependence. This is what we do next. Therefore, this geometry is a good approximation to that of a slowly rotating (quasi) spherical gravitational mass. Magnetic field around a spherical gravitational In the weak field limit (non-relativistic weak gravita- field, with non-homogeneous constitutive rela- tional sources) the above metric is given by the following tions: Let us consider the effect of relaxing the con- expressions in spherical and (isotropic) Cartesian coordi- dition of homogeneity in the constitutive relations, as- nates, respectively [13] suming a radial dependence of the permeability tensor in vacuum, i.e., 2GM 2GM ds2 = c2 1 dt2 1+ dr2 − c2r − c2r H = µ−1(r)√ gg Bj . (3.34) k 0 − jk 4GJ r2 dθ2 + sin θ2dϕ2 + sin2 θdϕdt,(3.38) We introduce the ansatz − c2r r ¯ Schw µ0(r)= µ0 1+ γ¯ , (3.35) 2GM 2GM r ds2 = c2 1 dt2 1+ − c2r˜ − c2r˜ × and choose γ¯ = 1. The above results are therefore gen- 4GJ eralized into dx2 + dy2 + dz2 + dt(xdy ydx),(3.39) × c2r˜3 − θ µ0 rSchw B (r,π/2) 1+ , ∝ r5 r wherer ˜ x2 + y2 + z2, as before. The first three µ r θ 0 Schw terms correspond≡ to the Schwarzschild geometry in the B (r, θ) 4 1+ , ∝ r5 sin θ r first order approximation.p In fact, these expressions are for the magnetic field due to a ring of current (in the a particular case of the general line-element for non- equatorial plane) and to the astrophysical jet, respec- relativistic stationary sources, understood as a linear per- tively. The strength of this effect in comparison with the turbation of Minkowski background spacetime [12] case with a homogeneous magnetic permeability tensor, 2Φ(g) 2Φ(g) is stronger near the horizon ds2 = c2 1 dt2 1+ − c2 − c2 × δB rSchw | | = . (3.36) (dx1)2 + (dx2)2 + (dx3)2 +2A(g)dxidt,(3.40) B r × i i Contrary to what happened with the electric field case, with gαβ = ηαβ + hαβ, where Φ(g) and A(g) are the grav- the radial dependence of the permeability tensor en- itoelectric and gravitomagnetic potentials respectively, hances the magnetic field in comparison with the case which can be defined through the expressions below [12] of magnetic homogeneity of vacuum. In the limit, at the horizon of black holes, the magnetic field is stronger by 4Φ Ai h¯00 (g) , h¯0i (g) , (3.41) a factor of 2. This fact comes from the very nature of ≡ c2 ≡ c 10 and so-called Post-Newtonian corrections while the expansion up to first order represent the weak field limit and A dxi = η Aidxj = A~ dx,~ i ij − · is sometimes referred to as the Newtonian limit. How- 1 ever, strictly speaking there can be corrections to Newton h¯ h h h, h = ηµν h . (3.42) µν ≡ µν − 2 µν µν gravity even in the first order (weak field) limit, as in the For a brief review on gravitolelectromagnetism in the per- case of slow rotation where the gravitomagnetic potential turbative as well as geometric approaches, see [31]. In and corresponding vector field allows non-Newtonian pre- this perspective, we can clearly see that the influence dictions, such as frame dragging or Lens-Thirring effects of the object’s angular momentum in the surrounding [12]. In GR, we have β = γ = 1. Many of the so called spacetime express the gravitomagnetic effect (measured gravitational classical tests can be expressed in terms of around Earth in the Gravity Probe B experiment [32]). these parameters [12, 13]. For example, the deflection of The PPN (Parameterized Post-Newtonian) formal- light due to the Sun’s gravitational field, the precession ism allows one to parameterize geometrical gravitational of the perihelion of planetary orbits and the Shapiro time theories within the Solar System, or more generally delay of light signals in the Sun’s field (see [13], for exam- around spherically symmetric stationary (possibly rotat- ple). Will [16, 17], Ni [18, 19] and Misner et al [20] used ing) gravitational sources. Many dimensionless parame- ten parameters in the Beta-Delta notation and later, a ters appearing in Taylor expansions of the metric are thus different set of 10 parameters was used in the Alpha-Zeta defined describing deviations from GR. For example, two notation (see for example [21, 22]), but γ and β coincide of the most important ones, β and γ, can arise from an in both notations. In the first notation ∆1 and ∆2 are expansion in powers of GM/c2r of the most general static intimately related to the off-diagonal elements character- spherically symmetric spacetime [13] istic of a Kerr-like metric, since ∆1 measures how much dragging of inertial frames is produced by unit linear mo- 2 2 2 2 2 2 2 2 ds = D(r)c dt B(r)dr r dθ + sin θ dϕ mentum and ∆2 measures the difference between radial − − = weak field (Newtonian) + post Newt., (3.43) and transverse momentum on dragging of inertial frames. − All parameters are potentially useful to constrain al- obtaining ternative or extended geometric theories of gravity. For the slowly rotating object of our interest, when the ex- D(r)=1 Ψ(r) + (β γ)Ψ(r)2 + ..., 2 − − pansions on the dimensionless quantities GM/c r and B(r)=1+ γΨ(r)+ ..., (3.44) GJ/c3r2 are taken up to first order (in the most general axially symmetric metric), one arrives at the following where γ is a measure of the spatial curvature produced expressions in spherical and (isotropic) Cartesian coordi- by a unit rest mass, while β is a measure of how much nates non-linearity is present in the superposition law for gravity. The higher order terms in the expansion give the

2GM 2GM α 2GJ ds2 = c2 1 dt2 1+ γ dr2 r2 dθ2 + sin2 θdϕ2 + 1+ γ + 1 sin2 θdϕdt, (3.45) − c2r − c2r − 4 c2r 2GM 2GM α 2GJ ds2 = c2 1 dt2 1+ γ dx2 + dy2 + dz2 + 1+ γ + 1 dt(xdy ydx), (3.46) − c2r˜ − c2r˜ 4 c2r˜3 −

respectively, where α1 7∆1 + ∆2 4γ 4, measures 1. Gravitomagnetic coupling between electric and magnetic the extent of preferred frame≡ effects and− is− equal to zero fields in GR. The following bounds taken from Will [22] due to local (Solar System) tests should be taken under consid- We now consider Maxwell’s inhomogeneous equations eration by any gravitational theory in the background spacetime around a spherical mass with slow rotation in the weak field (first order) limit −5 −3 −4 (using Eq (3.45)). New gravitomagnetic terms appear γ 1 =2.3 10 , β 1 =3 10 , α1 < 10 , | − | × | − | × in the Gauss and Maxwell-Ampère laws due to the off- diagonal time-space metric component. The generalized where the first result was obtained through light deflec- Gauss law, Eq. (2.11), has now a mixture of electric and tion and time delay, the second is due to perihelion shift magnetic components. The coupling to spacetime geome- and the last from orbit polarization in Lunar Laser Rang- try and the gravitomagnetic effect provides magnetic cor- ing (LLR). rections, which vanish in the non-rotating regime. This 11 puts forward the interesting possibility of having even (t, R,φ,z˜ ), the metric becomes static magnetic fields as sources of electric fields. A possible application concerns the magnetic field 2GM around rotating neutron stars. This field feels the pres- ds2 = c2 1 dt2 ence of very strong gravity where the curvature of space-  − 1/2  c2 R˜2 + z2 time should not be neglected. In such astrophysical con-   ditions the theory here exposed suggests an induced elec-   tric field component due to the coupling between the 2GM 1+ γ dR˜2 + R˜2dφ2 + dz2 magnetic field and the gravitomagnetic character of grav-  1/2  − 2 ˜2 2 c R + z ity.   Another application comes from supermassive rotat-   α1 2GJ ˜2 ing black holes in the center of disk-shaped galaxies. It + 1+ γ + 3/2 R dtdφ, (3.47) 4 2 ˜2 2 is well know that the frame-dragging character of the c R + z gravitomagnetic field around the rotating mass can be understood as a differential rotation of the curved space- where R˜ is a radial coordinate related to the physical time around the rotation axis of the object (as seen from distance to the rotation axis. The Gauss law, in vacuum, an observer far away from the source). Such a rotation is for this case is analogous to what happens in tornadoes where, in contrast to rigid bodies, the angular velocity is higher to- 1 ˜ R˜ z 1 φ wards the center and decays with radial distance. Ac- ∂R˜ (RD )+ ∂zD + ∂φD =0. (3.48) R˜ R˜ cordingly, from the coupling between gravitomagnetism and electromagnetic fields within Maxwell’s equations, it We will assume for simplicity no magnetic fields. Then is natural to expect a frame-dragging effect on these fields we have the following linear (homogeneous and isotropic) (see also [33, 34]). Therefore, an electric field around a constitutive relations [10] super-massive rotating black hole would feel the differential rotation of spacetime resulting in a spiral pattern for k 0i 0k 00 ik D = ε0√ g(g g g g )Ei . (3.49) the electric field lines along the galactic equator. If this Let us consider an electric− field− with axial symmetry [E = relativistic (non-Newtonian) effect might provide some ˜ R˜ ˜ φ light into the understanding of the formation processes ER˜ (R, z)e + Eφ(R,z)e ], where we have assumed that of spiral structures in galaxies, remains up to this mo- the Ez component is negligible near the equatorial plane ment a challenging and open question. (z = 0). Therefore, for this approximation, we have [D = R˜ ˜ φ ˜ Changing the Cartesian (isotropic) coordinates in D (R,z)eR˜ +D (R,z)eφ]. Thus, the Gauss law provides Eq. (3.46) to (axi-symmetric) cylindrical coordinates the solution DR˜(R,z˜ )= f(z)/R˜.

In the equatorial plane, setting z = 0, the radial component of the electric ﬁeld in the spacetime around the rotating massive object is given by

1/2 r2 J 2(4 + α +4γ)2 + 16c2m2(1 r /R˜)(1 + γr /R˜)R˜2 f(0) Sch − Sch Sch E ˜ = , (3.50) R 2 4ε0mcR˜ h (1 + γrSch/R˜) i

and for J =0 becomes clear that electrostatics and magnetostatics are no longer separated, but instead become coupled in the C (1 rSch/R˜) presence of gravitomagnetism. We assumed for simplicity ER˜ = − , (3.51) a vanishing magnetic field, but as previously mentioned, R˜ s(1 + γrSch/R˜) in general, even a static magnetic field will contribute where C is a constant. On the other hand, the Faraday to the electric field via this gravitomagnetic coupling induced by astrophysical sources with rotation. law gives ∂zEφ = ∂R˜ Eφ = 0, therefore Eφ is an arbitrary constant that can be set to zero without a significant loss of generality. We see that the gravitomagnetic term affects the gravitationally induced decay of the radial component of the electric field. Case of non-homogeneous and anisotropic elec- These astrophysical scenarios intend to illustrate pos- tric permittivity: Following our hypothesis that the sible gravitomagnetic effects, affecting electromagnetism electromagnetic properties of vacuum should reflect directly through the very nature of the field equations the local spacetime isometries, the constitutive relation in curved spacetime. In fact, in the most general case it (3.49) in this case (neglecting magnetic fields for simplic- 12 ity) is generalized into and isotropy of the electromagnetic properties of vacuum (the electromagnetic oscillations induced by gravitational R˜ R˜ 00 R˜R˜ D = ε (R,z˜ )√ gg g E ˜ , (3.52) waves were presented in another work [11]). In particular, − 0 − R φ φ 0φ 2 00 φφ the physical (possibly measurable) effects of spacetime D = ε (R,z˜ )√ g (g ) g g Eφ. (3.53) 0 − − geometry in electromagnetic fields, expressed in terms of The important physical idea here is that the electric per- the extended Gauss and Maxwell-Ampère laws, helps the mitivity should have axial symmetry and depend on the comparison with the usual results obtained from electro- ij R˜ ˜ φ ˜ z ˜ direction, i.e., (ε0) = diag[ε0 (R,z),ε0 (R,z),ε0(R,z)]. magnetism in Minkowski spacetime, making clearer the Let us assume as an Ansatz, the following expression physical interpretations of such effects. In the following, we briefly summarize the topics explored in this work: ˜ In the spacetime around spherical sources, the results R˜ 2GM 2GJR ε = ε0 1+¯γ + ∆¯ ,(3.54) confirm that gravity induces a (geometric) contribution 0  c2R˜ 3/2  c3 R˜2 + z2 to the decay of electric and magnetic fields along any   radial direction. This can be seen in Eqs. (3.6), (3.8),   whereγ ¯ and ∆¯ are dimensionless parameters. This cor- (3.30) and (3.33). In principle, even electric fields due to responds to the linear approximation of a Taylor expan- plane charged plates, which are uniform in the absence of sion in powers of the relevant dimensionless quantities gravity will manifest a spatial variability as is clear from related to mass and angular momentum. By considering Eq. (3.9). This effect could be tested in principle, un- for simplicity the equatorial plane, i.e z = 0, then the der appropriate experimental conditions similar to those used in the GP-B experiment, namely with recourse to result in Eq. (3.50) is generalized by replacing ε0 with ˜ drag-free motion of satellites in polar or equatorial orbits εR(R,z˜ = 0). As a consequence, in this case the electric 0 around a spherical mass, housing the probe (in this case field has a radial dependence which includes both the con- a capacitor) in vacuum under extreme low temperatures tribution from spacetime curvature (gravity) and a vari- achieved by cryogenics. able electric permittivity of vacuum. This is analogous to According to the hypothesis we explore here, in the what we had in the spherical symmetric case generalized curved geometry around a massive object the electric for the axial symmetric geometries. By turning off the permittivity of vacuum is not homogeneous, but has a angular momentum, we recover the spherical case. Fur- radial dependence instead. For black holes, this effect ther research is required to understand how the effects of influences electric fields more strongly near the horizon. having non-homogeneous (and anisotropic) permittivity The strength of the effect is maximum at the horizon and permeability tensors can have an impact on phys- where the magnitude of the electric field is 0.5 times ical observables. Electric and magnetic fields interact weaker than the case with a constant permittivity ten- with astrophysical plasmas in black hole accretion disks sor. In principle, this effect could be tested using electric and around neutron stars. These interactions depend on fields inside capacitors. The expressions (3.19) and (3.21) the coupling constants which are basically the electric we obtained can be generalized to include the PPN ap- and magnetic properties of the medium. Such interac- proach. Similar considerations apply to magnetic fields tions need to be carefully taken care of, for instance, us- around massive spherical (non-rotating) objects. ing magnetohydrodynamical computations. On the other The hypothesis for non-homogeneous permittivity hand, by changing the permittivity and permeability ten- and permeability (electromagnetic) properties of vaccum sors, the propagation properties of electromagnetic waves leads to the result that the speed of electromagnetic is affected. We briefly discuss these issues in the conclu- waves is not the same in every point around a massive sions. However, testing these hypotheses via observations object. Instead, it must have a radial dependence. Us- need further analysis using the wave equations and also ing the Ansatz considered in this work, Eqs. (3.17) and taking into account environmental effects, although this (3.35), in the first order approximation (with respect to is not considered in the present work. Taylor expansions in powers of rSch/r), we get

c0 1 IV. SUMMARY AND DISCUSSION c(r)= , c0 . (1 + rSch/r) ≡ ε0µ0 Building on previous work [10], in this paper, we ex- As a consequence, local observers could still agree about plored the physical applications of electrodynamics in the the velocity of light and the local conformal structure background of a (pseudo) Riemann spacetime manifold. of spacetime, i.e., the local light-cone, but the change The main electromagnetic effects induced by spacetime in the (local) light cone structure from one point to an- curvature addressed here include the following: gravita- other now has both the influence of the spacetime curva- tional contributions for the decay of electric and magnetic ture and the fact that the permittivity and permeability fields in spherically symmetric spacetime, magnetic con- tensors change. This prediction for a non-homogeneous tributions to the Gauss law due to the gravitomagnetic (but isotropic) speed of light in vacuum, in the spherical character of the spacetime around rotating objects and gravitational fields around massive objects, should have the effects of relaxing the assumptions of homogeneity observable consequences that need to be tested experi- 13 mentally. effect might provide some light into the understanding of In the axisymmetric spacetime around rotating sources, the formation processes of spiral structures in galaxies, electrostatics and magnetostatics are no longer sepa- remains a challenging and open question. rated, but instead become coupled due to the presence In any case, we found the radial component of an ax- of off-diagonal time-space metric components. We con- ially symmetric electric field, solution to the Gauss law, sidered the metric in Eq. (3.45) which corresponds to in the geometry given by Eq. (3.47). This solution again the spacetime around a rotating spherical mass in the confirms that gravitational fields can decrease the mag- weak field and slow rotation regime. This metric has nitude of electromagnetic fields in vacuum and the grav- the off diagonal component proportional to the angular itomagnetic term (proportional to the source’s angular momentum of the source. Such metric components, cor- momentum) also contributes to this effect. It is also ped- respond in linearised gravity to the components of the agogical to illustrate the role of the coupling between gravitomagnetic potential characteristic of frame drag- gravity and electromagnetism for testing different theo- ging (Lens-Thirring) effects around axisymmetric astro- ries of gravity. We also generalized the constitutive re- physical systems in rotation. The coupling to space- lations in this case, to introduce non-homogeneity and time geometry, in particular to gravitomagnetism, in- anisotropy in the permittivity tensor, corresponding to a duces magnetic corrections to the Gauss law, i.e., Eq. spacetime with axial symmetry. (2.11). This coupling suggests that even static magnetic In this work we briefly explored the idea that the elec- fields can act as sources of electric fields via the gravito- tromagnetic properties of physical vacuum should follow magnetic frame dragging character of spacetime around the spacetime isometries. This idea reinforces the debate rotating objects. In fact, the magnetic field around rotat- and research about the deep relation between gravity ing neutron stars feels the presence of very strong grav- and electromagnetic fields, spacetime and the nature of ity and therefore spacetime curvature should not be ne- vacuum. In particular, the conformal (casual) structure glected. In such astrophysical conditions, the theory here of spacetime is fundamentally connected to electromag- exposed, suggests an induced electric field component, netism and the constitutive relations. The electromag- generated by the coupling between the magnetic field and netic properties of vacuum and the conformal part of the the geometrodynamical character of gravity. Some work metric appear in the constitutive relations between the has been done in the past related to these issues (see fields and excitations, therefore, in principle the electro- [33, 34]) but much more can be investigated. This is an magnetic properties of vacuum should not be considered illustration of a gravitomagnetic effect affecting electro- to be independent from the spacetime isometries or more magnetism directly through the very nature of the field specifically to the local conformal structure. Following equations in curved spacetime. It opens our perspectives this reasoning one can even reduce the number of con- in the way we imagine the astrophysical environment of cepts in an attempt to unify different approaches to the such compact objects and also other sources of strong same problem. In this sense one can speak about the astromagnetic fields electromagnetic properties of the “spacetime medium”. As mentioned in Section III B, another possible appli- Following these ideas, the (local) conformal structure cation comes from supermassive rotating black holes in of spacetime becomes a more fundamental concept than the center of disk-shaped galaxies. It is believed that the the metric structure. This can suggest the change of interaction of these supermassive black holes with the paradigm in spacetime physics according to which space- surrounding galactic environment is an important ingre- time is not absolute but rather it is the local causal struc- dient in the formation and evolution of the whole Galaxy ture that is invariant. Therefore all observers agree on and in the formation and evolution of AGNs and stel- the (local) light cone related to some spacetime point but lar formation bursts. On the other hand, it is also well spacetime distances between any two events can be dif- known that the frame-dragging character of the gravito- ferent for different observers. This change of paradigm magnetic field around a rotating mass can be understood naturally implies theories of gravity in which the (lo- as a differential rotation of curved spacetime around the cal) gauge symmetries are not Poincarésymmetries but rotation axis of the object. Such a rotation is analogous rather those of the 15 parametric conformal group which to what happens in tornadoes where, in contrast to rigid includes the Weyl group and the proper conformal trans- bodies, the angular velocity is higher towards the center formations. Dilatations (which belong to the Weyl group, and decays with radial distance. Accordingly, from the together with the Poincarétransformations) and proper coupling between gravitomagnetism and electromagnetic conformal transformations change the spacetime line ele- fields within Maxwell’s equations, it is natural to expect a ment and obervers related by proper conformal transfor- frame-dragging effect on these fields. Therefore, an elec- mations cannot be inertial observers. This wider gauge tric field produced by a supermassive charged, rotating group naturally invites non-riemann geometries, in par- black hole, would feel the differential rotation of space- ticular, general affine geometries with curvature, torsion time resulting in a spiral pattern for the electric field lines and non-metricity tensors. along the galactic equator. This could induce currents Regarding electromagnetic waves in the presence of and stationary charge density spiral patterns on the sur- gravity, one can show that extra terms appear in the gen- rounding ionized gas. If this relativistic (non-Newtonian) eralized wave equations which deserves further research. 14

Indeed, going beyond the geometrical optics limit, light show that the gauge invariant wave equation for the deflection (null geodesics) and gravitational redshift are coupled electric and magnetic fields also includes similar not the only effects arising from the coupling between gravitomagnetic terms. These terms not only contribute light and gravity. More generally, all electromagnetic to the decay of the wave amplitude, but will also provide waves can experience gravitational effects on the ampli- a geometrically induced coupling between the dynamics tudes, frequencies and polarizations [11] (see also [35]). of the various electric and magnetic components which Besides, electric and magnetic wave dynamics become will most probably affect the polarization. Gravitomag- coupled in general, even in the Lorenz gauge. The cou- netic effects on electromagnetic waves deserve further pling between the dynamics of different components sug- research with potential applications for relativistic gests polarization effects. Moreover, there are terms in astrophysics related to Pulsars, AGNs and for the study the wave equations depending on the first derivatives of the electromagnetic counterpart of gravitational wave of the electromagnetic fields and similar terms for the sources. wave equations written in terms of the potentials (2.19). These terms might be responsible for a gravitational contribution to the decay of the oscillations, but formally, these are also compatible with the existence of longitudinal modes induced by spacetime curvature, since similar ACKNOWLEDGMENTS terms appear in wave equations for vector fields with longitudinal modes. The fundamental reason behind this is FC acknowledges financial support of the Funda¸cão the fact that in the presence of gravity the electric field para a Ciência e Tecnologia (FCT) through the grant in vacuum is no longer divergent-free (in the sense that PD/BD/128017/2016 and Programa de Doutoramento k ∂ Ek = 0). FCT, PhD::SPACE Doctoral Network for Space Sciences One6 can easily show that in the spacetime around (PD/00040/2012). FSNL acknowledges financial support a rotating astrophysical object, the equation for the of an Investigador FCT Research contract, with refer- potential includes the influence of the gravitomagnetic ence IF/00859/2012, funded by FCT/MCTES (Portu- term and therefore, one expects that the electromagnetic gal). This article is based upon work from COST Action field will experience a frame-dragging (Lens-Thirring) CA15117, supported by COST (European Cooperation effect due to the gravitomagnetism. In fact, one can in Science and Technology).

[1] Hubert F. M. Goenner, “On the history of unified field [10] F. Cabral and F. S. N. Lobo, “Electrodynamics and theories”, Living Rev. Relativity, 7, (2004), 2. Spacetime Geometry: Foundations”, Found. Phys. 47, [2] M. Blagojevic and F. W. Hehl, “Gauge Theories of Grav- no. 2, 208 (2017) [arXiv:1602.01492 [gr-qc]]. itation,” arXiv:1210.3775 [gr-qc]. [11] F. Cabral and F. S. N. Lobo, “Gravitational waves and [3] F. W. Hehl and Y. N. Obukhov, “Foundations of electrodynamics: New perspectives”, Eur. Phys. J. C 77, Classical Electrodynamics: Charge, Flux, and Metric”, no. 4, 237 (2017) [arXiv:1603.08157 [gr-qc]]. Birkhuser, Boston, 2003. [12] Hobson, Efstathiou, Lasenby, “General Relativity an in- [4] F. Gronwald, F. W. Hehl and J. Nitsch, “Axiomatics of troduction for physicists”, Cambridge University Press, classical electrodynamics and its relation to gauge field 2006. theory”, physics/0506219. [13] James B. Hartle, “Gravity. An Introduction to Einsteins [5] F. W. Hehl and Y. N. Obukhov, “A gentle introduction to General Relativity”, Pearson Addison Wesley, 2003. the foundations of classical electrodynamics: The mean- [14] L. D. Landau and E. M. Lifshitz, “The Classical Theory ing of the excitations (D,H) and the field strengths (E, of Fields”, Elsevier, 1975. B)”, physics/0005084. [15] B. F. Schutz, “A first course on general relativity”, Cam- [6] F. W. Hehl and Y. N. Obukhov, “Spacetime met- bridge University Press, 2003. ric from local and linear electrodynamics: A New ax- [16] C. M. Will, “Theoretical frameworks for testing relativis- iomatic scheme”, Lect. Notes Phys. 702 (2006) 163 tic gravity: The parametrized post-Newtonian formal- [gr-qc/0508024]. ism”, Dissertation (Ph.D.), California Institute of Tech- [7] “F. W. Hehl and Y. N. Obukhov, How does the electro- nology (1971). magnetic field couple to gravity, in particular to metric, [17] C. M. Will and K. Nordtvedt, Jr., “Conservation nonmetricity, torsion, and curvature?”, Lect. Notes Phys. Laws and Preferred Frames in Relativistic Gravity. I. 562 (2001) 479, [gr-qc/0001010]. Preferred-Frame Theories and an Extended PPN Formal- [8] W. Heisenberg, H. Euler, “Consequences of Diracs theory ism”, Astrophys. J. 177 (1972) 757. of positrons”, Z. Phys. 98 (1936) 714-732. [18] K. S. Thorne, C. M. Will and W. T. Ni, “Theoreti- [9] B. Mashhoon, “Nonlocal electrodynamics, in: Cosmol- cal Frameworks For Testing Relativistic Gravity: A Re- ogy and Gravitation”, Proc. VII Brasilian School of Cos- view”, OAP-247. mology and Gravitation, Rio de Janeiro, August 1993, [19] W. T. Ni, “Theoretical Frameworks For Testing Rela- M. Novello, editor. Editions Frontieres, Gif-sur-Yvette tivistic Gravity: 4. A Compendium Of Metric Theories (1994) pp. 245-295. Of Gravity And Their Post-newtonian Limits”, OAP- 15

268. netized neutron star”, Mon. Not. Roy. Astron. Soc. 322 [20] C. W. Misner, K. S. Thorne and J. A. Wheeler, “Gravi- (2001) 723 [astro-ph/0011316]. tation”, San Francisco 1973. [29] D. Bini, C. Germani and R. T. Jantzen, “Gravitoelec- [21] C. M. Will, “Theory And Experiment In Gravitational tromagnetism and the integral formulation of Maxwell’s Physics”, Cambridge, Uk: Univ. Pr. (1981). equations”, Int. J. Mod. Phys. D 10 (2001) 633 [22] C. M. Will, “The Confrontation between general rel- [gr-qc/0012068]. ativity and experiment”, Living Rev. Rel. 9 (2006) 3 [30] M. Visser, “The Kerr spacetime: A Brief introduction”, [gr-qc/0510072]. arXiv:0706.0622 [gr-qc]. [23] G. Brodin, M. Marklund and P. K. S. Dunsby, “Nonlin- [31] B. Mashhoon, “Gravitoelectromagnetism: A Brief re- ear gravitational wave interactions with plasmas”, Phys. view”, gr-qc/0311030. Rev. D 62 (2000) 104008 [gr-qc/0006030]. [32] C. W. F. Everitt et al., “Gravity Probe B: Final Results [24] C. Chicone and B. Mashhoon, “Gravitomagnetic Jets”, of a Space Experiment to Test General Relativity”, Phys. Phys. Rev. D 83 (2011) 064013 [arXiv:1005.1420 [gr-qc]]. Rev. Lett. 106 (2011) 221101 [arXiv:1105.3456 [gr-qc]] [25] A. Merloni, M. Vietri, L. Stella and D. Bini, “On grav- [33] J. Bik, V. Karas and T. Ledvinka, “Black holes itomagnetic precession around black holes”, Mon. Not. and magnetic fields”, IAU Symp. 238 (2007) 139 Roy. Astron. Soc. 304 (1999) 155 [astro-ph/9811198]. [astro-ph/0610841]. [26] E. Barausse, V. Cardoso and P. Pani, “Can envi- [34] V. Karas, O. Kopacek and D. Kunneriath, “Influence ronmental effects spoil precision gravitational-wave as- of frame-dragging on magnetic null points near rotat- trophysics?”, Phys. Rev. D 89 (2014) no.10, 104059 ing black hole”, Class. Quant. Grav. 29 (2012) 035010 [arXiv:1404.7149 [gr-qc]]. [arXiv:1201.0009 [astro-ph.HE]]. [27] B. W. Carrol,D. A. Ostlie, “An Introduction to Modern [35] M. Marklund, G. Brodin and P. K. S. Dunsby, “Radio Astrophysics”, Addison-Wesley, 1996. wave emissions due to gravitational radiation”, Astro- [28] L. Rezzolla, B. J. Ahmedov and J. C. Miller, “General rel- phys. J. 536 (2000) 875 [astro-ph/9907350]. ativistic electromagnetic fields of a slowly rotating mag-