Arxiv:0902.0560V2 [Gr-Qc] 6 Apr 2009 Tween (See, Gravitation of Theory [2])

Home , Principle of locality

arXiv:0902.0560v2 [gr-qc] 6 Apr 2009 tween (see, gravitation of theory [2]). classical advanced instance, other nonlocal for been the a have of On favor arguments in gravitation. heuristic nonlocal various at rel- hand, special arrive nonlocal to with ativity together employed be cannot it e,cnb nesoda ag hoyo the of theory gauge a as understood sys- be Dirac gauged can a tem, of framework the in electrodynamics, Einstein’s of is GR form in mathematical equivalence precise of for the principle reason that gener- underlying is nonlocal The ob- this direct theory. No Einstein’s a metric? of for the alization itself of suggests terms to method in connection vious expression the nonlocal require the- a example, for Einstein’s be one, in Should implemented ory? be should generalization how But a such process. be a generalization cannot this to theory from generalized relativity exempt general is then relativity theory, special nonlocal equiv- rela- if the special Therefore, of from means tivity. deduced by heuristically is, the- principle, Einstein’s (GR) alence relativity gravitation. general of theory of the ory is view of refer- point the and [1] to refer experi- therein. we possible cited theory, of ences the discussion and of a tests electrodynamics For mental to equation. applied Dirac been the has relativity cial a orcintrshv enpooe htaeinduced are high that sufficiently nonlo- proposed and by been analysis have critical terms detailed correction a cal of subject the been ∗ † lcrncaddress: Electronic lcrncaddress: Electronic omlfaeokfrannoa eeaiaino Einstei of generalization nonlocal a for framework formal A h dabhn hspprepot h nlg be- analogy the exploits paper this behind idea The h etsbett eadesdfo hsnonlocal this from addressed be to subject next The nseilrltvt hoy h rnil flclt has locality of principle the theory, relativity special In ag theories. gauge ewrs olclgaiy eea eaiiy telepara relativity, general gravity, nonlocal 11.10.Lm Keywords: 04.20.Cv, 04.50.Kd, numbers: Tohlin PACS the for mo basis matter. pre a our dark provide in the for can gravity evidence of approach of this nature character that The nonlocal result strated the the matter”. of that “dark generalizat show manifestation nonlocal with we a relativity approximation, for general linear ansatz the an in develop Working to paper this in h nlg ewe lcrdnmc n h translational the and electrodynamics between analogy The eateto hsc n srnm,Uiest fMissour of University Astronomy, and Physics of Department eateto hsc n srnm,Uiest fMissour of University Astronomy, and Physics of Department .INTRODUCTION I. [email protected] [email protected] nttt o hoeia hsc,Uiest fCologne, of University Physics, Theoretical for Institute accelerations nteoehn eko that know we hand one the On strictly 1.Ti olclspe- nonlocal This [1]. Dtd 6Arl2009, April 06 (Dated: oa;therefore, local; armMashhoon Bahram rerc .Hehl W. Friedrich U (1) l NonlocalGrav40.tex file lls gravity llelism oooi oriaeidcs hl re etr indi- indices. letters (tetrad) Greek Lorentz-frame anholonomic while cate indices, coordinate holonomic ydiag(1 by otaain esrdniy n the and density, F tensor contravariant a pctm nie u rm0t n nt r chosen are units and paper, 3 this to In 0 that from such appendices. run the indices of to spacetime mathematical discussion relegated Various A VI. are Sec. details in out. contained is pointed results briefly our are matter” “dark is it hand the other of theory the gauge On a that [3]. known O’Raifeartaigh see group; awl equations Maxwell ino hs olcleutosi em ftestan- the pres- the of interpreta- in terms an but of in theory to ence relativity equations devoted general nonlocal linearized nonlocal is in dard these V principles of Section connec- variational tion in of discussed theories. use are field the subtleties some with Sec. the IV tion in out Sec. equations work In field and III. gravitational approximation matters, linear nonlocal simplify To the linearized nonlocal II. in for Sec. work ansatz in we resulting explained the is and gravitation section this of rest extension. nonlocal a to proceed rela- general start of (GR equivalent to tivity teleparallel helpful so-called be 5]. of this may from method [4, it direct al. theory, a Einstein’s et recognize generalizing Nitsch cannot to instance, we equivalent if for is Consequently, see, Lagrangian, chosen theory; suitably Einstein’s gravity a of for theory teleparallelism that, a yields matter, spinless ij h electromagnetic The h nlg ewe hs rusi one u nthe in out pointed is groups these between analogy The ( = ∗ † o fEnti’ hoyo gravitation. of theory Einstein’s of ion dr matter”. “dark -unteteto h astrophysical the of treatment e-Kuhn ,B E, e,i reydsusd ti demon- is It discussed. briefly is del, , c n olclter seuvln to equivalent is theory nonlocal ing − ag hoyo rvt semployed is gravity of theory gauge || itd“akmte” hc sthe is which matter”, “dark dicted nta.Ti ilide rvd a to way a provide indeed will This instead. ) 02 ¨l,Grayand Köln, Germany 50923 .TeMnosimti esri given is tensor metric Minkowski The 1. = 1 ,Clmi,M 51,USA 65211, MO Columbia, i, ,Clmi,M 51,USA 65211, MO Columbia, i, = ) , ∂ − j H 1 .Electrodynamics A. , ij − − ) = ) sidcs ai etr indicate letters Latin indices, As 1). F ’ hoyo gravitation of theory n’s ji J excitation oain esr ufilthe fulfill tensor, covariant a , i h rpriso h resulting the of properties The ∂ , [ i F H translation jk ij ] ( = 0 = ,H D, , edstrength field ru,for group, = ) −H (1) ji , 2 with i as electric current density. Here i, j, ... =0, 1, 2, 3 Expression for Electrodynamics Transl. Gauge Theory J i i are (holonomic) coordinate indices. The homogeneous conservation law ∂iJ = 0 ∂iTα ∼ 0 equations (1)2 can be solved by the potential ansatz ij i ij i inh. field eqs. ∂j H = J ∂j H α ∼ Tα j α j Fij =2∂[iAj]. (2) force density fi = Fij J fi ∼ Cij Tα α hom. field eqs. ∂[iFjk] = 0 ∂[iCjk] = 0 The two Maxwell equations have to be supplemented by α α a constitutive law that relates the excitation to the field potentials Fij = 2∂[iAj] Cij = 2∂[iej] strength. In conventional vacuum electrodynamics, one TABLE I: Analogy between electrodynamics and gravitody- takes the local, linear, and isotropic ‘constitutive law’ namics. The approximation sign employed in the last column indicates validity only in the linear regime. ij k[i j]l ij = √ gg g Fkl = √ g F . (3) H D − E − This law is equally valid for the flat Minkowski spacetime theory formulated in terms of the calculus of exterior dif- of special relativity as well as for the curved Riemannian ferential forms. Recent works on teleparallelism can be spacetime of general relativity. In the former case, in found, amongst others, in the articles of Aldrovandi et Cartesian coordinates, the metric gij is constant, in the al. [12], Itin [13, 14], Maluf et al. [15, 16], Mielke [17], latter case it becomes a field governed by the Einstein Obukhov et al. [18–20], So et al. [21], and Tung et al. equations. If we substitute Eq. (3) into Eq. (1)1, then [22]; moreover, the two monographs of Blagojević[23] the Maxwell equations are expressed in terms of the field and Ort´ın[24] should be consulted as well as the refer- strength Fij alone. ences cited therein. The inhomogeneous and the homoge- In a nonlocal theory, on the other hand, equation (3) neous gravitational field equations are respectively given is generalized such that the nonlocal constitutive law ex- by presses ij as a nonlocal spacetime relation in terms of H ∂ ij i =∗ i , ∂ C α =0∗ . (7) Fkl. In fact, it is possible to give an effective interpreta- j H α −Eα Tα [i jk] tion of Mashhoon’s nonlocal electrodynamics in this way The star over the equal sign means that the correspond- [6]. That is, one refers the excitation αβ = e αe β ij , i j ing equation is only valid in a suitable frame; this will the field strength F = ei ej F H, and the metricH αβ α β ij be discussed in Sec. II. The nonlinear correction terms in g = ei ej g to non-inertial frames e = ei ∂ that αβ α β ij α α i Eq. (7) represent the energy-momentum tensor density characterize the tetrads of accelerated observers. With 1 of the gravitational gauge field, respect to such frames, the Maxwell equations read [6, 7] i 1 i β jk β ik 1 1 α := e α(Cjk β)+ Cαk β . (8) ∂ αβ C α γβ + C γ αβ = α , (4) E −4 H H βH − 2 βγ H 2 βγ H J As source, we have on the right-hand side of the inho- ∂ F C δF = 0 . (5) [α βγ] − [αβ γ]δ mogeneous field equation the energy-momentum tensor γ i j γ density of matter. In analogy with Eq. (2), the field Here Cαβ := 2e αe β∂[iej] is the object of an- α α holonomity (“structure constants”) measuring the non- strength Cij is defined in terms of the potential ei by inertial nature of the frame. Thus in this approach to α α Cij =2∂[iej] , (9) nonlocal electrodynamics, the nonlocal expression that relates the excitation to the field strength is given by [6] so that the homogeneous gravitational field equations (7)2 are thus identically satisfied. αβ αγ βδ (x) = √ gg g Fγδ(x) Clearly, there is a similarity between the electrody- H − αβγδ ′ ′ 4 ′ namic and gravitational cases, as indicated in Table I. + χ (x, x )Fγδ(x )d x . (6) Z In electrodynamics we have one potential Ai as covec- tor, corresponding to the one-parameter U(1) group; in Eq. (6) represents a new axiom that, together with the the gravitational case, because we have four linearly in- Maxwell equations (1), determines a nonlocal theory of dependent translations in spacetime, we have four covec- α electrodynamics. tors as potentials ei , the tetrad components with α = 0ˆ, 1ˆ, 2ˆ, 3.ˆ Correspondingly, whereas we have six components of the electromagnetic field strength Fij = Fij , we B. Translational gauge theory of gravity have 4 6 components of the gravitational field− strength α × α Cij = Cji . If one compares Eq. (1) with Eq. (7), Let us now turn to the gauge theory of the transla- one recognizes− the quasi-Maxwellian nature of the grav- tion group, which includes GR||. We take the field equa- itational case, only that we have four times more field tions from [8, 9]; however, we change some conventions equations in Eq. (7) and that the gravitational energy- i in order to conform with [7, 10]. For the tensor calculus momentum tensor density α emerges because all phys- background one should compare Schouten [11]. In Ap- ical fields including gravityE are gravitationally ‘charged’, pendix A we provide a brief sketch of translational gauge that is, all fields carry energy-momentum. 3

The analog of Eq. (3), the local, linear, and isotropic our work considerably. In particular, the torsion of the constitutive law in GR|| is [5, 8] Weitzenb¨ock spacetime becomes

α α α α β ∗ α α ij ∗ e 1 ij [ij] [i j]γ T := Dϑ = dϑ +Γβ ϑ = dϑ = C , (14) α = C α Cα +2e αC γ , (10) ∧ H κ 2 − or, equivalently, α where e := det(ei ). Here κ =8πG, where G is Newton’s α α α k α constant of gravitation. If Eq. (10) is substituted into the Tij = 2D[iej] =2∂[iej] 2Γ[ij] ek − inhomogeneous field equations (7)1, the resulting equa- =∗ 2∂ e α = C α . (15) tions have been shown to be equivalent to the Einstein [i j] ij equations [8]. Taking the electromagnetic case (6) as a Accordingly, the gravitational field strength, in the spe- prototype, we tentatively expect that a suitable nonlo- cial “gauge” (13), is represented by the object of an- cal generalization of the constitutive relation (10) would holonomity Cα. This concludes the geometrical set-up lead to a nonlocal generalization of Einstein’s theory of of a Weitzenböck spacetime and we will henceforth drop gravitation. It is interesting to contemplate the nature of the star over the equal sign. this analogy, which is the basis of the present paper. The We now turn to physics in GR||, which is determined linear and possibly nonlocal electrodynamic constitutive via the variation of the action δS = 0, where law is generally valid for sufficiently weak electromagnetic fields, since the constitutive tensor is assumed to 4 S = ( g + m) d x . (16) be independent of the field strength. Therefore, we ex- Z L L pect that the same holds in the linear approximation for the nonlocal theory of gravitation that we develop below Here g and m are respectively the gravitational and on the basis of this analogy. matterL LagrangianL densities. Out of the field strength α Finally, we must mention that this approach to nonlo- C , we can construct a gravitational Lagrangian g that cal gravitation is essentially different from other nonlocal is—in analogy with electrodynamics—quadraticLin the modifications of general relativity; see, for example, [25], field strength. The explicit form of g can be left open for [26], and the references cited therein. the moment. However, we can introduceL quite generally ij the components of the excitation α that are related to those of the field strength by H II. ANSATZ FOR A NONLOCAL THEORY OF GRAVITY ij ∂ g α := 2 L α . (17) H − ∂Cij The arena for teleparallelism is the Weitzenböck space- 1 i time in which we have the orthonormal tetrad (coframe) In terms of differential forms we have Hα = 2 Hijαdx α α i j ij 1 ijkl ijkl ∧ ϑ = ei dx . Then the local (anholonomic) metric is dx , with α = 2 ǫ Hklα, where ǫ is the totally H 0123 ηαβ = diag(1, 1, 1, 1) and we can determine the antisymmetric Levi-Civita symbol with ǫ = 1. Be- − − − (holonomic) coordinate components of the metric by cause of the presumed quadratic nature of g and the Euler theorem on homogeneous functions, weL find α β gij = ηαβei ej . (11) 1 α ∂ g 1 ij α α = C L = C , (18) Simple algebra yields e := det(ei )= √ g. g ij α α ij L 2 ∂Cij −4H If a coframe is a coordinate frame, then− the object of anholonomity which expresses the gravitational Lagrangian in terms of ij α α α the excitation α and the field strength Cij such that C := dϑ (12) ij is linear inH C α. H α ij α The field equations of GR|| given in Eq. (7) follow from vanishes; in general, C is a 2-form which decomposes the variational principle of stationary action (16) with according to Cα = 1 C αdxi dxj such that Eq. (12) in 2 ij i := δ /δe α as the source. We note that, in lin- components results in Eq. (9).∧ Tα Lm i αβ βα αβ i ear approximation, the gravitational energy-momentum The connection 1-form Γ = Γ = Γi dx in a complex (8) will vanish and only the first term of the Weitzenböck spacetime is teleparallel.− That is, we can left-hand side of Eq. (7)1 will survive. introduce a suitable global Cartesian tetrad frame with Having thus established the geometry and the respect to which the connection vanishes Lagrange-Noether machinery for the teleparallel theory αβ ∗ of gravitation, we turn next to the explicit choice of the Γ =0 . (13) gravitational Lagrangian. General relativity is recovered via equation (10), which may be written as Throughout the rest of our paper we will work in a global frame that obeys Eq. (13). Then all covariant deriva- ij e k[i j]l C β tives reduce to partial derivatives and we can simplify α = g g ηαβ kl , (19) H Dκ E 4

′ with the modified field strength Differentiating this equation, we find that Ωai(x, x ) = ′ Ωia(x, x ) are smooth dimensionless two-point tensors α 1 α α α γ such that Cij := Cij C [ij] +2e[i Cj]γ . (20) 2 − ′ lim Ωai(x, x )= gai(x) . (28) x′→x − In particular, along the same line of thought as in [6], we wish to consider a nonlocal constitutive law for Thus a possible nonlocal generalization of Eq. (19) is gravitation. To prepare the ground, let us introduce the given by invariant proper infinitesimal distance ds between two 1 neighboring events in Weitzenböck spacetime ab (x) = g(x) U(x, x′)ΩaiΩbj Ω χ(x, x′) c κ ck H − p− Z 2 i j k ′ ′ 4 ′ ds = g dx dx (21) Cij (x ) g(x ) d x , (29) ij ⊗ × − ′ p ′ and define a geodesic between two fixed events P ′ and where U(x, x ) is unity except when x is in the future of x, in which case U vanishes; in Minkowski spacetime, P to be the path that is an extremum of the spacetime ′0 0 distance between P ′ and P , this means that events with x > x are excluded from the domain of integration in Eq. (29). Moreover, χ(x, x′) P is a scalar given, for instance, by δ ds =0 . (22) ′ ′ ′ Z ′ χ(x, x )= δ(x x )+ Kˆ (x, x ) , (30) P − where the Dirac delta-function δ(x x′) is defined via This path is given by the geodesic equations − d2xi i dxj dxk δ(x x′)ϕ(x′) g(x′) d4x′ = ϕ(x) (31) + =0 , (23) Z − − ds2 jk ds ds p for any smooth function ϕ(x). The scalar kernel Kˆ (x, x′) where in Eq. (30) denotes the nonlocal deviation from a local constitutive law. The main nonlocal constitutive relation i 1 = gil (g + g g ) (24) (29) can therefore be expressed as jk 2 lj,k lk,j − jk,l 1 ab (x)= g(x)[Cab (x) U(x, x′)ΩaiΩbj Ω are the Christoffel symbols. c κ c ck H p− − Z We assume that two causally separated events are con- Kˆ (x, x′)C k(x′) g(x′) d4x′] , (32) unique ij nected by a timelike or null geodesic; more gener- × p− ally, in the spacetime region under consideration, there where the nonlocal deviation from general relativity is exists a unique geodesic joining every pair of events. made explicit. More complicated nonlocal constitutive It then proves useful to employ the world-function Ω, relations are certainly possible; however, Eq. (32) is the which denotes half the square of the proper distance from simplest one involving an unknown scalar kernel Kˆ . ′ ′ P : x = ξ(ζ0) to P : x = ξ(ζ1) along the geodesic path We emphasize that the domain of applicability of xi = ξi(ζ). That is, we define [27] Eq. (32) is not physically restricted; in particular, a nonlocal gravity theory of this type is not directly re- 1 ζ1 dξi dξj lated to nonlocal special relativity [1]. That is, if the Ω(x, x′)= (ζ ζ ) g dζ . (25) 1 0 ij whole heuristic machinery of Einstein’s principle of equiv- 2 − Zζ0 dζ dζ alence could be employed in this case, one would have It turns out that Ω is independent of the affine parameter to conclude—on the basis of nonlocal special relativity— ζ; moreover, the integrand in Eq. (25) is constant by that gravity should be nonlocal for sufficiently high gravi- Eq. (23). The main properties of Ω(x, x′) are summarized tational “accelerations”. However, Einstein’s principle of in Appendix B. equivalence cannot be used in this way, as its precise for- To distinguish coordinate indices that refer to x from mulation in GR is the very embodiment of locality. Thus those that refer to x′, we henceforth use indices a, b, c, ... nonlocal gravity is in this sense decoupled from nonlocal at x and i, j, k, ... at x′. Thus we define special relativity. It is important to note that our ansatz (32) is non- ∂Ω ∂Ω linear as well as nonlocal, since, among other things, the Ω (x, x′)= , Ω (x, x′)= , (26) a ∂xa i ∂x′i metric tensor (11) is quadratic in the gravitational poten- α ′ tials ei . The scalar kernel Kˆ (x, x ) could be given, for ′ and note that covariant derivatives at x and x com- ˆ −1 2 4 instance, by K =Ω +L0, where L0 is a constant char- mute for any bitensor. It follows from the results of Ap- acteristic length. It could also involve the Weitzenböck pendix B that invariants [8] ab ij ijk ijk j ik 2Ω = g ΩaΩb = g ΩiΩj . (27) CijkC , CkjiC , Cij C k (33) 5 at x and x′ as well as scalars formed from the covariant theories is contained in Edelen [28]. The expressions for ′ k ij i derivatives of Ω(x, x ). To illustrate the latter, consider, Cij and C k in terms of ψ j can be easily determined. a for example, Ωa and Ωi in Eq. (26). The vector Ω is We find from Eq. (36) that ′ tangent at x to the geodesic connecting x to x and the α α α length of Ωa is equal to that of this geodesic by Eq. (27); ei = δi + ψ i , (38) in Minkowski spacetime, Ωa = xa x′a. Similarly, the so that from (9), vector Ωi has the same length, is tangent− to the geodesic ′ k k at x , and is directed away from x; in Minkowski space- Cij =2ψ [j,i] . (39) time, Ωi = x′i xi. Hence one can form coordinate a α − i α ′ scalars Ω ea (x) and Ω ei (x ), which turn out to be of Moreover, particular significance in Sec. V. This brief discussion indicates that many options indeed exist for generating a gij = ηij + hij , hij =2ψ(ij) . (40) ˆ ′ scalar kernel K(x, x ) in our ansatz (32). It follows from (20) that In the following section, the nonlocal theory of gravity is treated in the linear approximation; therefore, Cij 1 i j j i [ij] k = hk , hk , + ψ ,k Eq. (32) will be employed in linearized form. In par- −2 − ticular, U(x, x′)Kˆ (x, x′) will be evaluated in Minkowski i j j l j i i l + δ ψ, ψ δ ψ, ψ , (41) spacetime, since Weitzenböck spacetime in the special k − l , − k − l , gauge (13) reduces to Minkowski spacetime when the where ψ = η ψij . α α limit ij gravitational potentials ei reduce to δi . In this , Let us first examine the local part of the Lagrangian we define (37). It follows from ′ ′ ′ K(x, x )= U(x, x )Kˆ (x, x ) . (34) 1 1 δ (local) = Cij δC k δCij C k (42) Lg −4κ k ij − 4κ k ij III. GRAVITATIONAL FIELD EQUATIONS IN that the first term in Eq. (42) is in effect given by LINEAR APPROXIMATION 1 Cij δψk (43) −2κ k,j i The gravitational field equations are obtained from the variation of the action, d4x, where is the La- via Eq. (39) and integration by parts, assuming that δψk L L i grangian density given by theR sum of the correspond- vanishes on the boundary, and hence neglecting such ing quantities for gravitation g and matter m, namely boundary terms in Eq. (43). Applying this same pro- L L = g + m, such that cedure to the second term in Eq. (42), we find after a L L L detailed calculation that the result is again given by δ m = √ gT i . (35) L α α 1 δei − Cij δψk , (44) −2κ k,j i i Here Tα is the energy-momentum tensor of matter, so i i so that in the absence of the nonlocal term (K = 0) that √ gTα = α given in Eq. (7). The gravitational − T the variation of the action can be expressed in the linear part g is given by Eq. (18) in analogy with electrody- namicsL [7]. For the nonlocal theory under consideration, approximation as ij α is now given by Eq. (32); however, to simplify mat- H 1 ij i k 4 ters, we will work in the linear approximation. That is, C k,j + Tk δψ i d x =0 . (45) Z −κ we assume that

i i i Thus for K = 0, the gravitational field equations are e α = δα ψ α , (36) Cij i − ∂j k = κTk ; moreover, we prove in Appendix C that i these are precisely Einstein’s field equations according where the nonzero components of ψ α are such that i to the standard general linear approximation scheme. A ψ α 1. In this approximation, the holonomic co- | | ≪ general proof, valid in the non-linear regime, is already ordinate indices and the anholonomic tetrad indices are i indistinguishable. It then follows from Eq. (32) and Ap- contained in [8]. We note that ∂iTk = 0 is implied by pendix B that the field equations. We now turn to the nonlocal part of the Lagrangian. 1 In this case, the corresponding two terms in the variation = Cij (x)C k(x) Lg −4κ k ij of the gravitational action are 1 k ij 4 Cij (x) K(x, y)C k(y)d y , (37) 1 k ij 4 −4κ Z δ g(nonlocal) = δCij (x) K(x, y) C k(y)d y L −4κ Z where K is the scalar kernel in the weak-field approxima- 1 k ij 4 Cij (x) K(x, y) δC k(y) d y . (46) tion. A discussion of action principles in nonlocal field −4κ Z 6

As before, the ﬁrst term in Eq. (46) is easily shown to matter” and hence moving it to the right-hand side of lead to Eq. (54). In this way, one can recover Einstein’s theory but with “dark matter”. We explore this possibility in 1 k ∂K(x, y) ij 4 δψ i(x) C k(y)d y . (47) the rest of this paper. −2κ Z ∂xj Consider, for instance, the possibility that K is an even function of x y; that is, For the second term in Eq. (46), we can write its contri- − bution to the variation of the total action as K(x, y)= F (x y) , (55) − 1 k ij 4 4 j j K(x, y)Cij (x) δC k(y) d y d x . (48) where F (z) = F ( z). Then, ∂K/∂x = ∂K/∂y and −4κ Z Z − − Eq. (54) can be written as The domains of integration in the double integral are the Cij Cij ∂ k(y) 4 i same; therefore, it is possible to switch x and y in Eq. (48) ∂j k(x)+ K(x, y) j d y = κTk , (56) to get Z ∂y

which implies, via Appendix C, that Tki = Tik. Here we 1 k ij 4 4 k K(y, x)Cij (y) δC k(x) d y d x . (49) have assumed that the derivatives of ψ i vanish on the −4κ Z Z boundary of the spacetime region under consideration. Hence the second term in the nonlocal part of the La- grangian is equivalent to IV. NONLOCAL FIELD EQUATIONS 1 ij k 4 δC k(x) K(y, x)Cij (y)d y . (50) −4κ Z The treatment of the previous section has been based Now applying to Eq. (50) the same procedure we used on the assumption that nonlocal gravitational field equa- in the derivation of Eq. (44), we find after a detailed tions should be obtained from an action principle that calculation that the result is incorporates the nonlocal constitutive relation (32). Al- ternatively, this relation could be employed in the field 1 k ∂K(y, x) ij 4 δψ i(x) C k(y)d y . (51) equations of teleparallelism. The purpose of this section −2κ Z ∂xj is to address certain subtleties involved in these possibil- Let us define the symmetric kernel K by ities. 1 K(x, y)= [K(x, y)+ K(y, x)] . (52) 2 A. The “general” field equations of teleparallelism

Then, combining equations (47) and (51), we find In electrodynamics the Maxwell equations in (1) can 1 be derived from electric charge and magnetic flux con- δ (nonlocal) = δψk (x) servation, see [7]. No action principle is necessary. If Lg −κ i one considers, for example, local and linear magnetoelec- ∂K(x, y) ij 4 ij 1 ijkl C (y)d y . (53) tric matter, the constitutive law is = 2 χ Fkl, with j k H × Z ∂x a constitutive tensor density χijkl = χjikl = χijlk − − This means that the field equations of the nonlocal theory that has 36 independent components. This model can in the linear approximation are given by be derived from a Lagrangian, provided the constitutive tensor obeys additionally the Onsager-type relation χijkl = χklij ; that is, if 15 of its 36 components vanish. ij ∂K(x, y) ij 4 i ∂j C k + C k(y)d y = κTk . (54) Nevertheless, even if this condition is not fulfilled, one Z ∂xj can still find reasonable applications in physics, though i We note that ∂iTk = 0 is still identically satisfied in in this case irreversible processes have to be taken into the nonlocal theory. For a general symmetric kernel K, account. the energy-momentum tensor Tij is not symmetric in gen- The situation is similar for the energy-momentum ten- eral, since the nonlocal part in Eq. (54) is not in general sor density of the electromagnetic field. Starting from the symmetric. This poses no basic difficulty as Tij , given Lorentz force density, substitution of the inhomogeneous in general by Eq. (35), is not symmetric in general. On Maxwell equations and partial integration result in the the other hand, the requirement that Tij be symmetric (canonical) Minkowski energy-momentum tensor density (as in GR, for instance) would impose restrictions on the 1 kernel K. j = δj F kl + F jk . (57) i 4 i kl ik A natural interpretation of the basic nonlocal gravita- T − H H tional field equations (54) can be obtained by consider- Apparently no Lagrangian is necessary for the derivation ing the nonlocal term to be an effective source for “dark of this expression [7, 29]. 7

The gravitational field equations of the translational Alternatively, however, we can equally well substitute the gauge theory, as formulated in Eq. (7), should also have constitutive relation (59) in equation (7)1, which reduces general validity within this framework, independently in the linear regime to the new nonlocal field equations of the existence of the Lagrangian. For the homoge- ∂K(x, y) neous equations (7)2 this is apparent since they represent ij ij 4 i ∂j C k + C k(y)d y = κTk . (60) just the first Bianchi identities in a Weitzenböck space- Z ∂xj time. For the inhomogeneous equations (7)1 this can be seen as follows: We start with the energy-momentum The two nonlocal field equations only differ in their ker- conservation law for matter in the linear approximation nels: Eq. (54) involves K, while Eq. (60) involves K. The question of whether such nonlocal equations as (60) ∂ i 0. We can “solve” it by i = ∂ ij i and iTα ∼ Tα j H α − Eα can be derived from a variational principle is beyond the ij = ji , with some unknown correction term i. H α −H α e e Eα scope of our investigation; for a related study in connec- These are already the inhomogeneous field equations. We tion with acceleration-induced nonlocality see [30]. e e i e just have to determine α explicitly. Clearly it adds to In a similar way as in Sec. III, one can recover a direct i ij E i i α , namely ∂j α = α + α ; then, it is natural to nonlocal generalization of Einstein’s theory with a sym- T H Te E i metric energy-momentum tensor of matter by assuming interpret α ase the energy-momentume tensor density of the gravitationalE field. In analogy with Eq. (57) we ex- that K(x, y) is a function of x y. In this case, Eq. (60) − pect it to havee the form of Eq. (8); we now drop the tildes can be written as i ij from α and α. Therefore the inhomogeneous field Cij E H Cij ∂ k(y) 4 i equations extracted from the energy-momentum conser- ∂j k(x)+ K(x, y) j d y = κTk , (61) vatione law of mattere in this heuristic manner have exactly Z ∂y the same form as (7)1 with (8). Thus we may postulate which should be compared and contrasted with Eq. (56). (7) as the general field equations valid independently of These inequivalent nonlocal field equations both admit the existence of a Lagrangian. a natural interpretation in terms of dark matter as By contrast, if we have a Lagrangian, the field equa- described in detail in the following section; however, tions turn out to be [8] Eq. (61) has an advantage over Eq. (56) in that it involves a causal kernel K, while in the symmetric kernel ∂ ij ei + C β ik = i . (58) j H α − αLg αk H β Tα K of Eq. (56) past and future are treated in the same manner. If = 1 C α ij is substituted into (58), we recover Lg − 4 ij H α (7)1. Consequently, we have shown that the “general” field equations (7)1 with (8) are correct if a Lagrangian V. “DARK MATTER” exists. Accordingly, they are certainly one consistent and reasonable generalization of the Lagrangian-based equa- Let us now consider an important consequence of our tions (58). nonlocal equations for the gravitational field. To keep our discussion general, we start with an equation that B. Ambiguity in the field equations has the same form as Eqs. (54) and (60), except for a kernel , which we take to be equal to K or K, respectively, dependingK on whether one adopts the Lagrangian- We now wish to address a certain ambiguity that is based approach of Sec. III or the more direct approach encountered in implementing a nonlocal constitutive law of Sec. IV. We henceforth assume, for simplicity, that in as in our present work. A local constitutive law relating the former case K(x, y) is an even function of x y and the excitation ij to the gravitational field stregth Cij − H k k in the latter case K(x, y) is simply a function of x y. can be used in the derivation of the field equations of the Then, the arguments of the previous sections imply− that theory in either of two equivalent ways: it could be em- the nonlocal modification of Einstein’s gravitational field ployed in the Lagrangian (18) that then results, via the equations for a symmetric energy-momentum tensor Tij variational principle of stationary action, in the desired may be expressed as field equations or, alternatively, it could be directly substituted in the field equations (7)1. These two methods 4 Gij (x)+ (x, y)Gij (y)d y = κTij (x) , (62) in general produce different results for a nonlocal con- Z K stitutive law, however. It is interesting to illustrate this point for the case at hand in the linear approximation. where Gij is given by In this regime, we have the linear constitutive law i ik G j = ∂kC j (63) ij ij ij 4 κ α(x)= C α(x)+ K(x, y)C α(y)d y , (59) H Z and represents Einstein’s tensor in the linear approximation. which, when inserted in Eq. (18) results, as described The nonlocal term in equation (62) can be interpreted in detail in Sec. III, in the nonlocal field equations (54). in terms of an effective energy-momentum tensor for 8

“dark matter” by writing Eq. (62) as where Q is reciprocal to P , T ∞ Gij (x)= κ [Tij (x)+ ij (x)] , (64) Q(x, y)= P (x, y) . (73) − n nX=1 where the “dark matter” component is given by the symmetric energy-momentum tensor Moreover, Pn and Q are functions of x y, since the integration in Eq. (71) extends over the− entire three- 1 4 dimensional Euclidean space. Thus “dark matter” is Tij (x)= (x, y)Gij (y)d y . (65) −κ Z K the convolution of matter and the reciprocal kernel in this case. The mathematical implications of this fact are Equation (65) may be written in a more transpar- treated in Appendix E. ent form by making use of the methods of Appendix D. It is interesting to consider the confrontation of our That is, suppose that the integral equation (62) can be nonlocal theory in the linear approximation, represented solved by the method of successive substitutions (cf. Ap- by equations (62)-(67), with experimental data. Let us pendix D) and that the resulting infinite series is uni- first determine the Newtonian limit of our nonlocal grav- formly convergent. Then, as shown in Appendix D, it ity. This follows directly from Eq. (64) in the standard is possible to introduce a reciprocal scalar kernel such manner and we find the Poisson equation that R 2Φ=4πG(ρ + ρ ) , (74) ∇ D 4 Gij (x)= κTij(x)+ κ (x, y)Tij (y)d y . (66) Z R where Φ is the Newtonian potential. Here ρ is the density of matter and ρD, Thus the “dark matter” component in equation (64) is 3 given by ρD(t, x)= Q(x, y)ρ(t, y)d y , (75) Z 4 is the corresponding density of “dark matter”. We note Tij (x)= (x, y)Tij (y)d y , (67) Z R that in the general linear approximation of our theory the kernel (x, y) and its reciprocal (x, y) are given func- which is the integral transform of matter Tij by the ker- tions inK Minkowski spacetime andR are thus independent nel (x, y). It is clear from Eq. (67) that in our model, R of any particular physical system. The solar-system tests “dark matter” should be quite similar in its character- of general relativity imply that our “dark matter” must istics to actual matter; for instance, for dust, the corre- be a rather small fraction of the actual matter in the so- sponding “dark matter” is pressure-free, while for radia- k lar system. This can be simply arranged with a suitable tion with Tk = 0, the “dark” energy-momentum tensor choice of the universal reciprocal kernel Q. T ij is traceless as well. It is interesting to illustrate these considerations in the The reciprocal kernel is given as an infinite series in case of the problem of dark matter in spiral galaxies [31– terms of iterated kernelsR constructed from (x, y); that K 37]. Imagine, for instance, the circular motion of stars is, in the disk of a spiral galaxy. Beyond the galactic bulge, ∞ the Newtonian acceleration of gravity for each star at 2 (x, y)= (x, y) . (68) radius x is given by v0 / x toward the center of the −R Kn | | | | nX=1 galaxy. Here v0 is the constant “asymptotic” speed of stars in accordance with the observed rotation curves of To ensure causality, it is useful to assume further that spiral galaxies. It follows from Poisson’s equation that the effective density of dark matter is essentially given (x, y)= δ(x0 y0)P (x, y) , (69) by v 2/(4πG x 2). Comparing this result with equation K − 0 (75) and setting| | where P is a function of x y. Then, all iterated kernels − v 2 as well as the reciprocal kernel are also of this general ρ(t, y)= Mδ(y) , ρ (t, x)= 0 , (76) form; indeed, D 4πG x 2 | | 0 0 where M is the effective mass of the galaxy and the di- n(x, y)= δ(x y )Pn(x, y) , (70) K − mensions of the galactic bulge have been ignored, we find where P1(x, y)= P (x, y) and v 2 Q(x, 0)= 0 . (77) 4πGM x 2 3 Pn+1(x, y)= P (x, z)Pn(z, y)d z . (71) | | − Z The reciprocal kernel Q is a function of x y; therefore, it follows from Eq. (77) that − Therefore, 1 1 Q(x, y)= (78) (x, y)= δ(x0 y0)Q(x, y) , (72) 4πλ x y 2 R − | − | 9

2 α j with a universal length parameter λ = GM/v0 . Thus ei and its dual e β, a metric gij , and a flat connec- αβ βα taking due account of the observed rotation curves of tion Γi = Γi that is chosen to vanish globally αβ − spiral galaxies, Eqs. (74) and (75) with kernel (78) imply (Γi = 0). The gravitational field strength is given by α α Cij and the modified field strength by Cij . In this 1 ρ(t, y)d3y framework, which is capable of nonlocal generalization, 2Φ=4πG ρ(t, x)+ . (79) ∇ 4πλ Z x y 2 GR|| corresponds to a specific gravitational Lagrangian. | − | Working with a nonlocal “constitutive” kernel K(x, x′) in It is important to point out in passing a defect of the spe- the linear approximation, we construct an explicit nonlo- cific form of Eq. (78): the total amount of “dark matter” cal generalization of Einstein’s theory of gravitation that associated with a nonzero matter density ρ is infinite (see is consistent with causality. This theory can be refor- Appendix E); therefore, kernel (78) is too simple to be mulated as linearized general relativity but with “dark quite adequate for the task at hand. matter”, which mimics the contribution of nonlocal grav- For a point mass m with ρ(t, x)= mδ(x), the Newto- ity. We find that the effective energy-momentum tensor nian potential given by Eq. (79) can be expressed as of “dark matter” is simply the integral transform of the energy-momentum tensor of matter. Gm Gm x The application of our nonlocal model in the linear ap- Φ(t, x)= + ln | | . (80) − x λ λ proximation to the dark-matter problem in spiral galax- | | ies is in conflict with the empirical Tully-Fisher relation. The observational data for spiral galaxies indicate that It is possible that the situation can be significantly im- λ is of the order of a kpc, so that the logarithmic term proved with a suitable choice for the kernel in the non- in Eq. (80) is essentially negligible in the solar system. linear regime. However, a more basic theory is needed Moreover, the universality of the kernel in the linear ap- to determine the nonlocal constitutive kernel from first 2 proximation implies, via λ = GM/v0 , that for spiral principles. This is a task for the future. 2 galaxies M v0 , since λ should be independent of any particular physical∝ system. Acknowledgments. We are grateful to Jacob Beken- It is remarkable that in the simple considerations re- stein (Jerusalem), Yakov Itin (Jerusalem), JoséMaluf garding Eqs. (79) and (80)—based on our linear approx- (Brasilia), Eckehard Mielke (Mexico City), and Yuri imation scheme—we have recovered the significant pro- Obukhov (Moscow) for helpful comments. posal put forward by Tohline and further developed by Kuhn et al. [38, 39] to solve the dark matter problem by a natural modification of the Newtonian law of grav- APPENDIX A: TRANSLATIONAL GAUGE itation. An interesting discussion of the Tohline-Kuhn THEORY IN EXTERIOR CALCULUS scheme is contained in the review paper of Bekenstein α α i [39]. Despite various successes, the main drawback of The coframe 1-form ϑ = ei dx represents the gravi- this approach appears to be the violation of the Tully- tational potential; here i, j, ... =0, 1, 2, 3 are (holonomic) 4 Fisher relation, which implies that M v0 [38, 39]. To coordinate indices and α, β, ... = 0ˆ, 1ˆ, 2ˆ, 3ˆ are (anholo- ∝ j agree with the empirical Tully-Fisher law, it seems nec- nomic) frame indices; the frame eβ = e β∂j is dual to essary to go beyond the linear approximation scheme. α i α α j j the coframe, that is, ei e β = δβ and ei e α = δi . The Within the general framework of this work, but going field strength of gravitation is the object of anholonomity beyond the linear approximation as well as the specific 2-form Cα := dϑα. constitutive model employed thus far, it may be possi- Spacetime is described by a Weitzenböck geometry β β i ble in principle to have a relation of the general form with a teleparallel connection 1-form Γα =Γiα dx and of Eq. (67) with a kernel that is highly dependent upon with the curvature 2-form the particular system under consideration. This could provide a natural way to interpret the observational ev- β β γ β 1 β i j Rα := dΓα Γα Γγ = Rijα dx dx (A1) idence for dark matter as the nonlocal manifestation of − ∧ 2 ∧ the gravitational interaction. It remains to elucidate the that vanishes physical origin of the constitutive kernel that has been β the starting point of our investigation. Rα =0 . (A2)

We decompose the connection into a Riemannian (Levi- Civita) part and the contortion, VI. CONCLUSION Γ β = Γ β K β , (A3) α α − α To develop a nonlocal generalization of GR, it proves with the following definitionse of the torsion 2-form and useful to approach Einstein’s theory via its equivalent the contortion 1-form, respectively, within teleparallelism gravity, namely, GR||. Therefore, we work in Weitzenböck spacetime with a tetrad field T α := Dϑα = dϑα +Γ α ϑβ ,T α =: Kα ϑβ . (A4) β ∧ β ∧ 10

As in a Euclidean space, one can likewise choose in a On the other hand, Weitzenb¨ock space suitable frames such that the connec- ∂Ω ∂Ω tion vanishes everywhere: δΩ= δxa + δx′i , (B2) ∂xa ∂x′i Γ β =Γ βdxi =0∗ . (A5) α iα so that Eq. (A5) should be understood as a choice of a certain ∂Ω dxb gauge, which is possible since the curvature R β van- α a = (ζ1 ζ0)gab(x) , ishes everywhere. Many calculations are simpliﬁed in this ∂x − dζ gauge. ∂Ω dx′j = (ζ ζ )g (x′) . (B3) The Lagrangian 4-form is the sum of a gravitational ∂x′i − 1 − 0 ij dζ part and a matter part It is possible to see from Eq. (23) that the integrand in α a a L = Lg(ϑ , C )+ Lm(ψ, dψ, ϑ ) . (A6) Eq. (25) is indeed constant; therefore,

The field equation reads a b ′ 1 2 dx dx Ω(x, x ) = (ζ1 ζ0) gab(x) δLm 2 − dζ dζ dHα Eα = =Σα . (A7) − δϑα 1 dx′i dx′j = (ζ ζ )2g (x′) . (B4) α 1 i 2 1 − 0 ij dζ dζ Here we have the 2-form Hα := ∂Lg/∂C = 2 Hijαdx j − ∧ dx as the excitation and the 3-form Σα as the energy- It follows from Eqs. (B3)-(B4) that Eq. (27) is satisfied; momentum density of matter, whereas the 3-form Eα moreover, Ω = 0 for a null geodesic, Ω = 1 τ 2 for a time- represents the energy-momentum density of the gravita- 2 like geodesic of length τ and Ω = 1 σ2 for a spacelike tional field [10] 2 geodesic of length σ. − E := e L + e Cβ H . (A8) In Minkowski spacetime Ω is given by α α⌋ g α⌋ ∧ β Equations (A7) with (A8) represent the Lagrangian- 1 Ω(x, x′)= η (x′i xi)(x′j xj ) . (B5) based form of the field equations. If we express the 2 ij − − gravitational energy-momentum 3-form E in terms of α According to our convention, η = diag(1, 1, 1, 1); a Minkowski type current [or, equivalently, if we substi- ij − − − tute L = 1 Cα into Eq. (A8)], we find the general hence, g − 2 ∧Hα field equations 1 Ω(x, x′)= (t′ t)2 (x′ x)2 . (B6) 1 2 − − − dH H (e Cβ) Cβ (e H ) =Σ∗ . (A9) α − 2 β ∧ α⌋ − ∧ α⌋ β α In this case, we find that This is the exterior-form version of the field equations ∂2Ω(x, x′) (7)1 with Eq. (8). Ωai = a ′i = ηai , (B7) Usually Hα is, like in electrodynamics, linear in the ∂x ∂x − field strength while gαβ ⋆ (1) β (2) β (3) β 2 2 Hα = (a1 C + a2 C + a3 C ) . (A10) ∂ Ω ∂ Ω κ Ωab = = ηab , Ωij = = ηij . (B8) ∂xa∂xb ∂x′i∂x′j Here (I)Cβ are the different irreducible pieces of Cβ and ⋆ represents the Hodge star. The “Einsteinian” choice (I) β APPENDIX C: SYMMETRY OF THE TENSOR for C , I = 1, 2, 3, turns out to be C j ∂j i k 1 a1 = 1 , a2 =2 , a3 = . (A11) − 2 C j The purpose of this appendix is to show that ∂j i k is In this case, Eq. (A10) is distinguished from all linear a symmetric tensor, so that in the field equations locally Lagrangians in that it becomes Lorentz covariant. ij i ∂j C k = κTk (C1)

APPENDIX B: PROPERTIES OF THE the energy-momentum tensor of the source is symmetric WORLD-FUNCTION Ω (Tik = Tki). In fact Eq. (C1) is identical to Einstein’s ﬁeld equations in the linear approximation. It follows from Eq. (41) that Consider a variation of Eq. (25) that changes the end- points, then ij 1 i ij 1 ij ∂j C k = hk hkj, + h ,kj j ζ1 −2 − 2 ′ dξ i δΩ(x, x ) = (ζ1 ζ0) gij δξ . (B1) 1 i jl i i − dζ δkh ,jl + δk ψ ψ, k , (C2) ζ0 −2 − 11

jl jl where we have used the relation 2ψ ,jl = h ,jl. Define Then, Eq. (D2) can be written as the trace-reversed gravitational potentials h as ik b 1 f(x)= φ(x)+ r(x, y)f(y)dy . (D5) h = h η h , (C3) Za ik ik − 2 ik It is clear from equations (D1) and (D5) that the kernels ij where h = ηij h =2ψ. Then, replacing hij by hij +ηij ψ k and r are reciprocal of each other. everywhere in Eq. (C2) results in Let us note here some properties of these kernels. It follows from Eq. (D3) that j 1 1 j 1 j 1 jl ∂ C = h + h + h η h , (C4) b j i k −2 ki 2 k ,ij 2 i ,kj − 2 ki ,jl kn+p(x, y)= kn(x, z)kp(z,y)dz , (D6) Za which is manifestly symmetric in i and k. From Eqs. (C1) and (C4), we get where p = 1, 2,... . Using Eqs. (D3) and (D4) we find that j j jl hik + hi ,jk + hk ,ji ηikh ,jl =2κTik . (C5) b − − k(x, y)+ r(x, y) = k(x, z)r(z,y)dz These are exactly the same as Einstein’s field equations Za in the linear approximation. b = r(x, z)k(z,y)dz . (D7) Za

APPENDIX D: LIOUVILLE-NEUMANN Moreover, if k is symmetric, k(x, y) = k(y, x), then all METHOD OF SUCCESSIVE SUBSTITUTIONS iterated kernels as well as r(x, y) are likewise symmetric.

Consider a linear integral equation of the second kind APPENDIX E: CONVOLUTION KERNELS given by

b Suppose that the spatial kernels in Sec. V are such that φ(x)= f(x)+ k(x, y)φ(y)dy , (D1) Za P (x, y)= p(x y) ,Q(x, y)= q(x y) . (E1) − − where a and b are constants. We seek a solution of this The purpose of this appendix is to point out the conse- Fredholm equation by the method of successive substitu- quences of the convolution theorem for a system of den- tions. That is, we replace φ in the integrand by its value sity ρ such that the density of “dark matter” ρD is a given by equation (D1). Repeating this process eventu- convolution of ρ and q. In general, the spatial kernels ally leads to an inﬁnite series of the type p and q could depend on the characteristics of the par-

b ticular system under consideration. The main results of φ(x) = f(x)+ k(x, y)f(y)dy interest here are Eqs. (69) to (75). Let us note that the Za particular solution of Poisson’s equation (74) is given by b b + k(x, z)k(z,y)f(y)dy dz + . (D2) ρ(t, y)+ ρD(t, y) 3 Z Z ··· Φ(t, x)= G d y , (E2) a a − Z x y | − | If this series is uniformly convergent, we have a solution provided Φ vanishes sufficiently fast at spatial infinity. It of the integral equation (D1). This solution is unique follows from Eq. (75) that this Newtonian potential can in the space of real continuous functions on the interval also be expressed as [a,b]; a generalization of this result to the space of square- integrable functions is contained in Tricomi [40]. Φ(t, x)= G S(x y)ρ(t, y)d3y , (E3) Let us define the iterated kernels kn, n =1, 2,..., by − Z −

k1(x, y) = k(x, y) , where b 1 q(z)d3z kn+1(x, y) = k(x, z)kn(z,y)dz . (D3) S(r)= + . (E4) Z r Z r z a | | | − | These functions occur in the inﬁnite series of Eq. (D2). Let us assume that the functions of interest here can We deﬁne the reciprocal kernel r(x, y) such that be expressed as Fourier integrals; that is, for a function f, ∞ r(x, y)= kn(x, y) . (D4) ˆ ik·r 3 − νf(r)= f(k)e d k , (E5) nX=1 Z 12 where fˆ(k) is given by It follows that 1 1 ˆ −ik·r 3 qˆ(k)= . (E11) νf(k)= f(r)e d r . (E6) 4(2π)1/2λ k Z | | For the sake of simplicity, we have introduced here a con- Thus for any M > 0, MD = . ∞ stant parameter ν, Once q is determined from Eq. (E8), it is possible to work out its reciprocal kernel. In fact, it follows from 3/2 ν = (2π) . (E7) Eqs. (71) and (73) thatp ˆ1(k)=ˆp(k),

We are interested in Fourier integral transforms of real pˆn+1(k)= νpˆ(k)ˆpn(k) , (E12) functions; therefore, fˆ∗(k)= fˆ( k). − It follows from the convolution− theorem for Fourier in- and tegrals that Eq. (75) implies ∞ qˆ(k)= pˆn(k) . (E13) ρˆD(k)= νρˆ(k)ˆq(k) . (E8) − nX=1 Hereq ˆ(k) is a dimensionless function such that It is then straightforward to show that

MD pˆ(k) νqˆ(0)= , (E9) qˆ(k)= . (E14) M − 1+ νpˆ(k) where M = νρˆ(0) is the mass of the system under con- Therefore, sideration and MD = νρˆD(0) is the corresponding “dark” mass. Consider, for instance, kernel (78) that is associ- qˆ(k) ated with the discussion of the rotation curves of spiral pˆ(k)= . (E15) − 1+ νqˆ(k) galaxies in Sec. V; that is, 1 1 These results are consistent with the fact that p and q q(r)= . (E10) are reciprocal of each other. 4πλ r 2 | |

[1] B. Mashhoon, Ann. Phys. (Berlin) 17, 705 (2008) reprinted (Dover, Mineola, NY, 1989). [arXiv:0805.2926]. [12] R. Aldrovandi, J. G. Pereira, and K. H. Vu, Gen. Relativ. [2] B. Mashhoon, Ann. Phys. (Berlin) 16, 57 (2007) Gravit. 36, 101 (2004) [arXiv:gr-qc/0304106]. [arXiv:hep-th/0608010]. [13] Y. Itin, Gen. Relativ. Gravit. 34, 1819 (2002) [arXiv:gr- [3] L. O’Raifeartaigh, The Dawning of Gauge Theory qc/0111087]. (Princeton University Press, Princeton, NJ, 1997). [14] Y. Itin, in Classical and Quantum Gravity Research [4] J. Nitsch, in Proceedings of the 6th Course of the School Progress, edited by M. N. Christiansen and T. K. Ras- of Cosmology and Gravitation on ‘Spin, Torsion, Ro- mussen (Nova Science Publishers, Hauppauge, NY, 2008) tation, and Supergravity’, Erice, Italy, 1979, edited by [arXiv:0711.4209]. P. G. Bergmann and V. de Sabbata (Plenum, New York, [15] J. W. Maluf, M. V. O. Veiga, and J. F. da Rocha- 1980) p. 63. Neto, Gen. Relativ. Gravit. 39, 227 (2007) [arXiv:gr- [5] J. Nitsch and F. W. Hehl, Phys. Lett. B 90, 98 (1980). qc/0507122]. [6] U. Muench, F. W. Hehl, and B. Mashhoon, Phys. Lett. [16] J. W. Maluf and S. C. Ulhoa, Phys. Rev. D 78, 047502 A 271, 8 (2000) [arXiv:gr-qc/0003093]. (2008) [arXiv:0807.0255]. [7] F. W. Hehl and Yu. N. Obukhov, Foundations of [17] E. W. Mielke, Ann. Phys. (NY) 219, 78 (1992). Classical Electrodynamics: Charge, Flux, and Metric [18] Yu. N. Obukhov and J. G. Pereira, Phys. Rev. D 67, (Birkhäuser, Boston, MA, 2003). 044016 (2003) [arXiv:gr-qc/0212080]. [8] F. W. Hehl, J. Nitsch, and P. Von der Heyde, in General [19] Y. N. Obukhov and G. F. Rubilar, Phys. Rev. D 73, Relativity and Gravitation, Volume 1, edited by A. Held 124017 (2006) [arXiv:gr-qc/0605045]. (Plenum, New York, 1980) p. 329; see, in particular, the [20] Y. N. Obukhov, G. F. Rubilar, and J. G. Pereira, Phys. Appendix. Rev. D 74, 104007 (2006) [arXiv:gr-qc/0610092]. [9] F. W. Hehl, in Proc. of the 6th Course of the School [21] L. L. So and J. M. Nester, in Proceedings of the 10th Mar- of Cosmology and Gravitation on ‘Spin, Torsion, Ro- cel Grossmann Meeting, Rio de Janeiro, Brazil, 2003, tation, and Supergravity’, Erice, Italy, 1979, edited by edited by M. Novello et al. (World Scientific, Hackensack, P. G. Bergmann and V. de Sabbata (Plenum, New York, NJ, 2005), Part B, p. 1498 [arXiv:gr-qc/0612062]. 1980) p. 5. [22] R. S. Tung and J. M. Nester, Phys. Rev. D 60, 021501 [10] F. W. Hehl, J. D. McCrea, E. W. Mielke, and (1999) [arXiv:gr-qc/9809030]. Y. Ne’eman, Phys. Rep. 258, 1 (1995). [23] M. Blagojević, Gravitation and Gauge Symmetries (IoP, [11] J. A. Schouten, Tensor Analysis for Physicists, 2nd ed. Bristol, UK, 2002). 13

[24] T. Ort´ın, Gravity and Strings (Cambridge University F. J. Lockman, Astrophys. J. Lett. 641, L109 (2006); Press, Cambridge, UK, 2004). P. Salucci et al., Mon. Not. R. Astron. Soc. 378, 41 [25] M. E. Soussa and R. P. Woodard, Class. Quant. Grav. (2007). 20, 2737 (2003) [arXiv:astro-ph/0302030]. [35] J. D. Bekenstein, Phys. Rev. D 70, 083509 (2004); er- [26] A. O. Barvinsky, Phys. Lett. B 572, 109 (2003) ratum, Phys. Rev. D 71, 069901(E) (2005) [arXiv:astro- [arXiv:hep-th/0304229]. ph/0403694]. [27] J. L. Synge, Relativity: The General Theory (North- [36] J. D. Bekenstein, Contemp. Phys. 47, 387 (2006) Holland, Amsterdam, 1971). [arXiv:astro-ph/0701848]. [28] D. G. B. Edelen, Nonlocal Variations and Local Invari- [37] J. P. Bruneton, S. Liberati, L. Sindoni, and B. Famaey, ance of Fields (American Elsevier, New York, 1969). [arXiv:0811.3143]. [29] H. Minkowski, Königliche Gesellschaft der Wissenschaf- [38] J. E. Tohline, in IAU Symposium 100, Internal Kine- ten zu Göttingen. Mathematisch-Physikalische Klasse. matics and Dynamics of Galaxies, edited by E. Athanas- Nachrichten, p. 53 (1908). soula (Reidel, Dordrecht, 1983), p. 205. J. R. Kuhn and [30] C. Chicone and B. Mashhoon, Ann. Phys. (Berlin) 16, L. Kruglyak, Astrophys. J. 313, 1 (1987). 811 (2007) [arXiv:0708.2744]. [39] J. D. Bekenstein, in Second Canadian Conference on [31] V. C. Rubin and W.K. Ford, Astrophys. J. 159, 379 General Relativity and Relativistic Astrophysics, A. Co- (1970). ley, C. Dyer, and T. Tupper, eds. (World Scientific, Sin- [32] M. S. Roberts and R. N. Whitehurst, Astrophys. J. 201, gapore, 1988), p. 68. 327 (1975). [40] F. G. Tricomi, Integral Equations (Interscience, New [33] Y. Sofue and V. Rubin, Annu. Rev. Astron. Astrophys. York, 1957). 39, 137 (2001). [34] C. Carignan, L. Chemin, W. K. Huchtmeier, and