Nonlocal Gravitomagnetism

TL;DR

This paper explores nonlocal gravity's theoretical framework, focusing on gravitomagnetic fields of rotating sources in the weak-field limit, and compares these results with standard general relativity's gravitoelectromagnetism.

Contribution

It provides a linear approximation analysis of gravitomagnetic fields in nonlocal gravity and compares them with classical GEM predictions, advancing understanding of nonlocal gravitational effects.

Findings

Nonlocal gravity predicts modifications to gravitomagnetic fields.

Results show differences from standard GEM in weak-field approximation.

Energy-momentum content of the GEM field is analyzed in the nonlocal context.

Abstract

We briefly review the current status of nonlocal gravity (NLG), which is a classical nonlocal generalization of Einstein's theory of gravitation based on a certain analogy with the nonlocal electrodynamics of media. Nonlocal gravity thus involves integro-differential field equations and a causal constitutive kernel that should ultimately be determined from observational data. We consider the stationary gravitational field of an isolated rotating astronomical source in the linear approximation of nonlocal gravity. In this weak-field and slow-motion approximation of NLG, we describe the gravitomagnetic field associated with the rotating source and compare our results with gravitoelectromagnetism (GEM) of the standard general relativity theory. Moreover, we briefly study the energy-momentum content of the GEM field in nonlocal gravity.

Equations

$${}^0{\Gamma }_{\alpha \beta }^{\mu }=\frac{1}{2}{g}^{\mu \nu }(g_{\nu \alpha ,\beta }+g_{\nu \beta ,\alpha }-g_{\alpha \beta ,\nu }).$$

$${}^0{R}^{\alpha }{}_{\mu \beta \nu }={\partial }_{\beta }{}^0{\Gamma }_{\nu \mu }^{\alpha }-{\partial }_{\nu }{}^0{\Gamma }_{\beta \mu }^{\alpha }+{}^0{\Gamma }_{\beta \gamma }^{\alpha }{}^0{\Gamma }_{\nu \mu }^{\gamma }-{}^0{\Gamma }_{\nu \gamma }^{\alpha }{}^0{\Gamma }_{\beta \mu }^{\gamma }.$$

$${}^0{G}_{\mu \nu }+\Lambda g_{\mu \nu }=\kappa T_{\mu \nu },{}^0{G}_{\mu \nu }:={}^0{R}_{\mu \nu }-\frac{1}{2}{g}_{\mu \nu }{}^{0}R,$$

$$g_{\mu \nu }={\eta }_{\hat{\alpha }\hat{\beta }}{e}_{\mu }{}^{\hat{\alpha }}(x)e_{\nu }{}^{\hat{\beta }}(x).$$

$$d{s}^2=g_{\mu \nu }d{x}^{\mu }d{x}^{\nu }={\eta }_{\hat{\alpha }\hat{\beta }}d{x}^{\hat{\alpha }}d{x}^{\hat{\beta }},d{x}^{\mu }=e^{\mu }{}_{\hat{\alpha }}d{x}^{\hat{\alpha }},d{x}^{\hat{\alpha }}=e_{\mu }{}^{\hat{\alpha }}d{x}^{\mu },$$

$${\Gamma }_{\alpha \beta }^{\mu }=e^{\mu }{}_{\hat{\rho }}{\partial }_{\alpha }{e}_{\beta }{}^{\hat{\rho }}.$$

$${\nabla }_{\nu }{e}_{\mu }{}^{\hat{\alpha }}=0,$$

$$C_{\mu\nu}{}^{\alpha}=\Gamma^{\alpha}_{\mu\nu}-\Gamma^{\alpha}_{\nu\mu}=e^{\alpha}{}_{\hat{\beta}}\Big{(}\partial_{\mu}e_{\nu}{}^{\hat{\beta}}-\partial_{\nu}e_{\mu}{}^{\hat{\beta}}\Big{)}\,,$$

$$K_{\mu \nu }{}^{\alpha }={}^0{\Gamma }_{\mu \nu }^{\alpha }-{\Gamma }_{\mu \nu }^{\alpha },$$

$$g_{\alpha \beta ,\gamma }={\Gamma }_{\gamma \alpha }^{\mu }{g}_{\mu \beta }+{\Gamma }_{\gamma \beta }^{\mu }{g}_{\mu \alpha },$$

$$K_{\mu \nu }{}^{\alpha }=\frac{1}{2}{g}^{\alpha \beta }(C_{\mu \beta \nu }+C_{\nu \beta \mu }-C_{\mu \nu \beta }).$$

$$\displaystyle{{}^{0}}G_{\mu\nu}=\frac{\kappa}{\sqrt{-g}}\Big{[}e_{\mu}{}^{\hat{\gamma}}\,g_{\nu\alpha}\,\frac{\partial}{\partial x^{\beta}}\,{\cal H}^{\alpha\beta}{}_{\hat{\gamma}}-\Big{(}C_{\mu}{}^{\rho\sigma}\,{\cal H}_{\nu\rho\sigma}-\frac{1}{4}\,g_{\mu\nu}\,C^{\alpha\beta\gamma}\,{\cal H}_{\alpha\beta\gamma}\Big{)}\Big{]}\,,$$

$${\mathcal{H}}_{\mu \nu \rho }:=\frac{\sqrt{-g}}{\kappa }{\mathfrak{C}}_{\mu \nu \rho },{\mathfrak{C}}_{\alpha \beta \gamma }:=C_{\alpha }{g}_{\beta \gamma }-C_{\beta }{g}_{\alpha \gamma }+K_{\gamma \alpha \beta }.$$

$$\frac{\partial }{\partial x^{\nu }}{\mathcal{H}}^{\mu \nu }{}_{\hat{\alpha }}+\frac{\sqrt{-g}}{\kappa }\Lambda e^{\mu }{}_{\hat{\alpha }}=\sqrt{-g}(T_{\hat{\alpha }}{}^{\mu }+{\mathbb{T}}_{\hat{\alpha }}{}^{\mu }).$$

$$\kappa {\mathbb{T}}_{\mu \nu }:=C_{\mu \rho \sigma }{\mathfrak{C}}_{\nu }{}^{\rho \sigma }-\frac{1}{4}{g}_{\mu \nu }{C}_{\rho \sigma \delta }{\mathfrak{C}}^{\rho \sigma \delta }.$$

$$\frac{\partial}{\partial x^{\mu}}\,\Big{[}\sqrt{-g}\,(T_{\hat{\alpha}}{}^{\mu}+\mathbb{T}_{\hat{\alpha}}{}^{\mu}-\frac{\Lambda}{\kappa}\,e^{\mu}{}_{\hat{\alpha}})\Big{]}=0\,,$$

$$L_g=\frac{c^3}{16\pi G}({}^{0}R-2\Lambda ).$$

$${{}^{0}}R=-\frac{1}{2}\,C^{\alpha\beta\gamma}\,\mathfrak{C}_{\alpha\beta\gamma}+\frac{2}{\sqrt{-g}}\,\frac{\partial}{\partial x^{\delta}}\,\Big{(}\sqrt{-g}\,C^{\delta}\Big{)}\,,$$

$${\mathbb{L}}_g=-\frac{c^3}{32\pi G}(C^{\alpha \beta \gamma }{\mathfrak{C}}_{\alpha \beta \gamma }+4\Lambda ).$$

$$C^{\alpha \beta \gamma }{\mathfrak{C}}_{\alpha \beta \gamma }=\frac{1}{2}{C}_{\alpha \beta \gamma }{C}^{\alpha \beta \gamma }+C_{\alpha \beta \gamma }{C}^{\gamma \beta \alpha }-2C_{\alpha }{C}^{\alpha }.$$

$$R^{\hat{\alpha }}{}_{\hat{\beta }\mu \nu }:=2\left({\partial }_{[\mu }{\Gamma }_{\nu ]\hat{\beta }}^{\hat{\alpha }}+{\Gamma }_{[\mu |\hat{\gamma }}^{\hat{\alpha }}{\Gamma }_{|\nu ]\hat{\beta }}^{\hat{\gamma }}\right)=0.$$

$${\Gamma }_{\mu \hat{\alpha }}^{\hat{\beta }}\stackrel{\ast }{=}0\text{(in a suitable frame)}.$$

$$T_{\mu \nu }{}^{\hat{\alpha }}:=2{\nabla }_{[\mu }{e}_{\nu ]}{}^{\hat{\alpha }}=2\left({\partial }_{[\mu }{e}_{\nu ]}{}^{\hat{\alpha }}+{\Gamma }_{[\mu |\hat{\beta }|}^{\hat{\alpha }}{e}_{\nu ]}{}^{\hat{\beta }}\right).$$

$${}^{(1)}{T}_{\mu \nu \rho }:=T_{\mu \nu \rho }-{}^{(2)}{T}_{\mu \nu \rho }-{}^{(3)}{T}_{\mu \nu \rho },{}^{(2)}{T}_{\mu \nu \rho }:=-\frac{2}{3}{T}_{[\mu }{g}_{\nu ]\rho },{}^{(3)}{T}_{\mu \nu \rho }:=T_{[\mu \nu \rho ]}.$$

$${\mathcal{H}}^{\mu \nu }{}_{\hat{\alpha }}:=-\frac{\partial {\mathcal{L}}_{\mathrm{TG}}}{\partial T_{\mu \nu }{}^{\hat{\alpha }}}.$$

$${\mathcal{L}}_{\text{TG}}\sim \frac{1}{\kappa }\left(a_1{}^{(1)}{T}^{\mu \nu }{}_{\alpha }{}^{(1)}{T}_{\mu \nu }{}^{\alpha }+a_2{}^{(2)}{T}^{\mu \nu }{}_{\alpha }{}^{(2)}{T}_{\mu \nu }{}^{\alpha }+a_3{}^{(3)}{T}^{\mu \nu }{}_{\alpha }{}^{(3)}{T}_{\mu \nu }{}^{\alpha }\right).$$

$${\mathcal{H}}^{\mu \nu }{}_{\alpha }=-\frac{\sqrt{-g}}{\kappa }\left(a_1{}^{(1)}{T}^{\mu \nu }{}_{\alpha }+a_2{}^{(2)}{T}^{\mu \nu }{}_{\alpha }+a_3{}^{(3)}{T}^{\mu \nu }{}_{\alpha }\right).$$

$$a_1=-1,a_2=2,a_3=\frac{1}{2}.$$

$${\mathcal{H}}_{\mu \nu \rho }=\frac{\sqrt{-g}}{\kappa }({\mathfrak{C}}_{\mu \nu \rho }+N_{\mu \nu \rho }),$$

$$N_{\mu \nu \rho }=-\int {\Omega }_{\mu {\mu }^{\prime }}{\Omega }_{\nu {\nu }^{\prime }}{\Omega }_{\rho {\rho }^{\prime }}\mathcal{K}(x,x^{\prime })X^{{\mu }^{\prime }{\nu }^{\prime }{\rho }^{\prime }}(x^{\prime })\sqrt{-g(x^{\prime })}{d}^4{x}^{\prime },$$

Full text

Nonlocal gravity (NLG), Gravitoelectromagnetism (GEM)

I Introduction

The standard formulation of general relativity (GR) involves the extension of classical physics expressed in Minkowski spacetime, with metric $dS^{2}=\eta_{\mu\nu}\,dX^{\mu}\,dX^{\nu}$, first to arbitrary curvilinear (“accelerated”) coordinates via the locality postulate and then to curved spacetime, with metric $ds^{2}=g_{\mu\nu}\,dx^{\mu}\,dx^{\nu}$, by means of Einstein’s principle of equivalence Einstein ; Mash1 ; Mash2 . Here, $\eta_{\alpha\beta}$ is the Minkowski metric tensor given by diag$(-1,1,1,1)$, latin indices run from 1 to 3, while greek indices run from 0 to 3. The theory is thus based on the Levi-Civita connection

[TABLE]

This symmetric connection is torsion free, but has Riemannian curvature

[TABLE]

A left superscript “0” will be employed throughout to designate all geometric quantities that are related to the Levi-Civita connection.

In the curved spacetime of general relativity, free test particles and light rays follow timelike and null geodesics, respectively. The correspondence with Newtonian gravitation is established via Einstein’s field equations Einstein

[TABLE]

where ${{}^{0}}G_{\mu\nu}$ is the Einstein tensor, ${{}^{0}}R_{\mu\nu}={{}^{0}}R^{\alpha}{}_{\mu\alpha\nu}$ is the Ricci tensor and ${{}^{0}}R=g^{\mu\nu}\,{{}^{0}}R_{\mu\nu}$ is the scalar curvature. Moreover, $T_{\mu\nu}$ is the symmetric energy-momentum tensor of matter (and nongravitational fields), $\Lambda$ is the cosmological constant and $\kappa:=8\pi G/c^{4}$. In GR, the gravitational field is identified with the Riemannian curvature of spacetime; therefore, spacetime is flat when gravity is turned off and we then work within the framework of the special theory of relativity.

Einstein’s general relativity has significant observational support. Indeed, GR is at present in good agreement with solar system data as well as data from astronomical binary systems. The recent detection of gravitation radiation due to binary mergers lends further support to Einstein’s theory of gravitation. On the other hand, in the current standard model of cosmology, which assumes the spatial homogeneity and isotropy of the universe, the energy content of the universe consists of about 70% dark energy, about 25% dark matter and about 5% visible matter. Dark energy is a kind of repulsive energy that permeates the universe and not only counteracts the attraction of matter, but causes accelerated expansion of the universe. The nature and origin of dark energy are unknown, but it should have positive energy density and negative pressure. It is uniformly distributed throughout space and though it exists everywhere, it is extremely difficult to detect locally. A possible candidate for dark energy is provided by the cosmological constant $\Lambda$. The existence of dark energy and dark matter indicates that we are almost completely ignorant about our universe. Most of the matter in the universe is currently thought to be in the form of certain elusive particles of dark matter that, despite much effort, have not been directly detected. The existence and properties of this dark matter have thus far been deduced only through its gravity. In modern astronomy, dark matter is needed to explain dynamics of galaxies, clusters of galaxies and structure formation in cosmology. However, it is possible that there is no dark matter at all and the theory of gravitation needs to be modified on the scale of galaxies and beyond in order to take due account of what appears as dark matter in astronomy and cosmology. A suitably extended theory of gravitation could then account for the observational data without any need for dark matter. The present paper is about an attempt in this direction; that is, the nonlocal aspect of gravity in NLG simulates dark matter. The main purpose of this paper is to briefly present the main features of NLG theory and develop a useful linear perturbation scheme involving nonlocal gravitoelectromagnetism.

Einstein’s theory of gravitation can be alternatively formulated within the framework of teleparallelism. In this approach to gravitation, the fundamental fields are the sixteen components of an arbitrary smooth orthonormal tetrad frame $e^{\mu}{}_{\hat{\alpha}}(x)$. The spacetime metric is then defined via the orthonormality condition

[TABLE]

Here, the hatted indices (e.g., $\hat{\alpha}$) refer to anholonomic tetrad—that is, local Lorentz—indices, while ordinary indices (e.g., $\alpha$) refer to holonomic spacetime indices. For instance, in

[TABLE]

the tetrad connects (holonomic) spacetime quantities to (anholonomic) local Lorentz quantities. A coordinate basis is holonomic, while a noncoordinate basis is anholonomic. For instance, given a coordinate system $x^{\mu}$, four coordinate lines pass through each event and for each $\mu=0,1,2,3$, the 1-form $dx^{\mu}$ is exact and hence integrable. On the other hand, for each $\hat{\alpha}=\hat{0},\hat{1},\hat{2},\hat{3}$, the 1-form $dx^{\hat{\alpha}}$ in Eq. (5) is in general not exact and hence nonintegrable. Holonomic systems are integrable, while anholonomic systems are nonintegrable. Holonomic and anholonomic indices are raised and lowered by means of $g_{\mu\nu}(x)$ and $\eta_{\hat{\alpha}\hat{\beta}}$, respectively. To change an anholonomic index of a tensor into a holonomic index or vice versa, we simply project the tensor onto the corresponding tetrad frame. We use units such that $c=1$, unless specified otherwise.

The chosen tetrad frame is employed to define the Weitzenböck connection We

[TABLE]

This nonsymmetric connection is curvature free, but has torsion. It follows from definition (6) that the tetrad frame is covariantly constant

[TABLE]

where $\nabla$ refers to covariant differentiation with respect to the Weitzenböck connection. Equation (7) implies that each leg of the tetrad field is parallel to itself throughout the manifold. That is, for each $\hat{\alpha}$, Eq. (7) is an expression of the parallel transport of the corresponding vector with respect to connection (6). Thus in this theory observers throughout spacetime have access to a global set of parallel vector fields that constitute the components of the tetrad frame field. This circumstance is the essence of teleparallelism; for example, two distant vectors can be considered parallel to each other if they have the same components with respect to the local tetrad frames.

It follows from Eqs. (4) and (7) that $\nabla_{\gamma}\,g_{\alpha\beta}=0$, so that the Weitzenböck connection is compatible with the metric. Thus in the framework under consideration here, we have one spacetime metric and two metric-compatible connections. It is therefore possible to introduce the torsion tensor

[TABLE]

and the contorsion tensor

[TABLE]

which are linearly related. To see this, we note that $\nabla_{\gamma}\,g_{\alpha\beta}=0$ implies

[TABLE]

which, via the Levi-Civita connection (1), leads to

[TABLE]

The torsion tensor is antisymmetric in its first two indices by definition; however, the contorsion tensor turns out to be antisymmetric in its last two indices. The torsion of the Weitzenböck connection and the curvature of the Levi-Civita connection are complementary aspects of the gravitational field within the framework of teleparallelism. Thus it is natural to express Einstein’s field equations in terms of the torsion tensor. The result is the teleparallel equivalent of general relativity, GR*||*, to which we now turn.

I.1 GR*||*

It follows from Eqs. (9) and (11) that one can write Einstein’s field equations in terms of the torsion tensor. To this end, one can prove after much algebra that the Einstein tensor is given by

[TABLE]

where we have introduced auxiliary torsion fields ${\cal H}_{\mu\nu\rho}$ and $\mathfrak{C}_{\alpha\beta\gamma}$,

[TABLE]

Here, $g:=\det(g_{\mu\nu})$, $\sqrt{-g}=\det(e_{\mu}{}^{\hat{\alpha}})$ and $C_{\mu}$ is the torsion vector $C_{\mu}:=C^{\alpha}{}_{\mu\alpha}=-C_{\mu}{}^{\alpha}{}_{\alpha}$. The Einstein field equations can thus be written within the framework of teleparallelism as

[TABLE]

Here, $\mathbb{T}_{\mu\nu}$ is the trace-free energy-momentum tensor of the gravitational field and is given by

[TABLE]

The antisymmetry of ${\cal H}^{\mu\nu}{}_{\hat{\alpha}}$ in its first two indices can be used to show that the law of conservation of total energy-momentum tensor in GR*||*, namely,

[TABLE]

follows from the gravitational field equations.

Let us recall here that GR field equations can be derived from an action principle involving a gravitational Lagrangian given by

[TABLE]

On the other hand, we find

[TABLE]

so that the corresponding Lagrangian for GR*||* is given by

[TABLE]

The special torsion invariant in Eqs. (18)–(19) can be expressed as a linear combination of the three independent algebraic invariants of the torsion tensor, namely,

[TABLE]

I.2 GR*||* as the Gauge Theory of the Translations Group

Fundamentally, teleparallelism and GR*||* can only be understood in the framework of a gauge theory of gravitation BlHe . Nowadays the strong and the electroweak interactions are described by means of gauge theories. For gravity this framework can be utilized as well.

Consider first matter in a Minkowski space. The source of gravity in Newton’s theory is the mass density; within special relativity it should be the energy-momentum tensor $T_{\mu\nu}$ instead. For an isolated material system, energy-momentum is conserved. This is the result of the rigid (often called “global”) translation invariance of the action function of the material system under consideration.

A rigid invariance is in contrast to the idea of field theory. Thus, in adopting the gauge doctrine, we postulate for the action function the invariance under local translations. This forces us to introduce 1+3 translational gauge potentials (nonholonomic frames) $e_{\mu}{}^{\hat{\alpha}}$ thereby deforming the Minkowski space $M_{4}$ to a Weitzenböck space $W_{4}$. Details of this procedure may be found in Ref. Hehl:1979 .

In $W_{4}$, the Lorentz rotations are not gauged, that is, the action is still invariant under rigid Lorentz rotations, exactly like in $M_{4}$. Accordingly, the $W_{4}$ connection $\Gamma_{\mu\hat{\alpha}}^{\hat{\beta}}$ is still flat:

[TABLE]

This guarantees that in a $W_{4}$ the parallel transport is still integrable. Accordingly, like in $M_{4}$, we can choose all over $W_{4}$ a suitable frame such that the connection vanishes everywhere:

[TABLE]

Instead of the curvature, $W_{4}$ carries a translational field strength torsion which, in analogy to electrodynamics, is represented by the curl of the translational potential $e_{\mu}{}^{\hat{\alpha}}$:

[TABLE]

In the teleparallel frame of Eq. (22), we have for the torsion $T_{\mu\nu}{}^{\hat{\alpha}}\stackrel{{\scriptstyle*}}{{=}}C_{\mu\nu}{}^{\hat{\alpha}}\,,$ see Eq. (8), where $C_{\mu\nu}{}^{\hat{\alpha}}$ is the object of anholonomity of Schouten Schouten:1954 . The torsion has three irreducible pieces ${}^{(I)}{}T_{\mu\nu}{}^{\alpha}$, for $I=1,2,3$. With the torsion vector $T_{\mu}:=-T_{\mu\nu}{}^{\nu}$, we have

[TABLE]

So far we reminded ourselves of the kinematics of a translational gauge theory (TG). With the gauge Lagrangian ${\cal L}_{\rm TG}={\cal L}_{\rm TG}(\partial e,e,\Gamma,g)$, we can address the dynamics by defining the gravitational translational field momentum (or translation excitation)

[TABLE]

Should we investigate a physical system which has no Lagrangian—in the case of irreversibility, e.g.—the excitation ${\cal H}^{\mu\nu}{}_{\hat{\alpha}}$ still makes physical sense, as we know, e.g., from electrodynamics and the inhomogeneous Maxwell equation.

The general quadratic TG Lagrangian carries three independent pieces:

[TABLE]

To the Lagrangian (26) we can add a Lagrange multiplier term for enforcing the teleparallel constraint, see Ref. Hehl:1979 . It turns out that we cannot allow spinning matter (other than as test particles) in such a teleparallel space. Accordingly, we have to decree, see page 52 of Ref. Hehl:1979 , that only scalar and electromagnetic matter be allowed in TG, since they do not carry dynamical spin and have, as a consequence, symmetric energy-momentum tensors.

The translational excitation of Lagrangian (26) reads,

[TABLE]

In a teleparallelism theory the three-parametric rigidly Lorentz invariant Lagrangian is a totally acceptable choice. It corresponds to a gauge theory of the translation group. However, as it so happens, amongst these three-parameter Lagrangians, up to an overall constant, there is only one Lagrangian that is locally Lorentz invariant, see Cho Cho . This theory, which we abbreviate by GR*||, is, for scalar and electromagnetic matter, equivalent to GR. The local Lorentz invariance is imposed from the outside, it is not necessary in a translational gauge theory. But it shows that GR can be really understood as a specific translational gauge theory. A Hilbert-Einstein Lagrangian is equivalent to a definite torsion square Lagrangian in the teleparallel limit. This is a big step forward in understanding GR. The constants for GR||* are found to be, see Ref. Muench:1998 :

[TABLE]

This set of constants is called the Einstein choice. Lagrangian (26), together with Eq. (28), and the attached field momentum (27) were the starting point for a classical nonlocal theory of gravity.

I.3 Nonlocal Gravity

A locality assumption runs through the standard theories of special and general relativity Mash1 ; Mash2 . For instance, to render an accelerated system in Minkowski spacetime relativistic, Lorentz transformations are applied in a pointwise manner all along the world line of the accelerated system. An accelerated observer is thus assumed to be physically identical with a hypothetical inertial observer that shares the same state, namely, position and velocity. The locality hypothesis originates from the Newtonian mechanics of classical point particles and its domain of validity is determined by the extent to which physical phenomena could be reduced to pointlike coincidences. However, wave phenomena are generally nonlocal by the Huygens principle. Moreover, Bohr and Rosenfeld have shown that the electromagnetic field measurement requires a certain average over a region of spacetime BR1 ; BR2 . To go beyond the locality assumption, one must include an average over the past world line of the accelerated observer. In this way, a nonlocal special relativity theory has been developed Mash3a ; Mash3b .

Can nonlocal special relativity be extended to include the gravitational interaction by means of Einstein’s principle of equivalence? Einstein’s principle is extremely local, however, and this approach encounters severe difficulties and has been abandoned. Instead, we use Einstein’s fundamental insight regarding the connection between inertia and gravitation as a guiding principle and develop nonlocal general relativity patterned after the nonlocal electrodynamics of media. To this end, we exploit the formal analogy between GR*||* and electrodynamics and introduce an average of the gravitational field into the field equations via a causal constitutive kernel HM1 ; HM2 ; Mash3c . In nonlocal gravity, the gravitational field is local, but satisfies partial integro-differential field equations.

In nonlocal gravity, as in the electrodynamics of media, we retain the gravitational field equations (14), but change the local constitutive relation (13) to

[TABLE]

where the new tensor $N_{\mu\nu\rho}$ involves a linear average of the torsion tensor over past events. More specifically, we assume that

[TABLE]

where $\Omega(x,x^{\prime})$ is Synge’s world function Sy , ${\cal K}$ is the scalar causal kernel of the nonlocal theory and $X_{\mu\nu\rho}(x)$ is a tensor that is antisymmetric in its first two indices and is given by

[TABLE]

Here, $p\neq 0$ is a constant dimensionless parameter and $\check{C}^{\mu}$ is the torsion pseudovector defined via the Levi-Civita tensor $E_{\alpha\beta\gamma\delta}$ by

[TABLE]

Finally, the gravitational field equation of nonlocal gravity (NLG) is given by

[TABLE]

where the energy-momentum tensor of the gravitational field, $\mathcal{T}_{\mu\nu}$, is now given by

[TABLE]

The total energy-momentum conservation law then takes the form

[TABLE]

No exact nontrivial solution of the nonlocal field equation (33) is known. In this connection, the main source of difficulty appears to be the complicated relation that introduces nonlocality into the theory, namely, Eq. (30). In a recent paper Puetzfeld:2019wwo , a simpler form of Eq. (30) has been suggested, where the bitensor $\Omega_{\mu\mu^{\prime}}$ is replaced by the parallel propagator $-g_{\mu\mu^{\prime}}$. It remains to determine whether this simplification could help in generating exact nontrivial solutions of NLG.

The arbitrary tetrad frame we adopted to develop GR*||* could be any smooth tetrad frame field in spacetime. At each event, any two tetrad frame fields are related by an element of the local Lorentz group. This circumstance is in agreement with the invariance of Einstein’s GR under the local Lorentz group, since Einstein’s theory ultimately depends only upon the metric tensor $g_{\mu\nu}$. The introduction of nonlocality into the theory may remove this pointwise 6-fold degeneracy of GR*||*. However, as expected, NLG remains invariant under the global Lorentz group.

II Nonlocal GEM

We are interested in the stationary gravitational field of a rotating astronomical body, which is assumed to be confined to a compact region of space. We work in the linear approximation of nonlocal gravity, since the gravitational field is assumed to be weak. In this regime, a certain analogy with classical electrodynamics HeOb turns out to be fruitful. Indeed, in linearized GR, the framework of gravitoelectromagnetism (GEM) has proved rather useful in describing and interpreting the gravitational effects of rotating masses. It is therefore interesting to develop this method in NLG. A preliminary account is already contained in Ref. Mash3c and will be further developed in this paper.

II.1 Linearized NLG

In the weak-field regime, we can write the chosen tetrad frame field in the form

[TABLE]

where $\psi_{\mu\nu}(x)$ is the first-order perturbation away from a background global inertial reference frame in Minkowski spacetime such that

[TABLE]

In this linear approximation scheme, the distinction between spacetime and tetrad indices disappears and it follows from Eq. (4) that $g_{\mu\nu}=\eta_{\mu\nu}+\psi_{\mu\nu}+\psi_{\nu\mu}$. Therefore, we can write

[TABLE]

The gravitational perturbation $\psi_{\mu\nu}$ is thus comprised of a symmetric metric part $\frac{1}{2}\,h_{\mu\nu}$ and an antisymmetric tetrad part $\frac{1}{2}\,\phi_{\mu\nu}$. In connection with the metric part, it is useful to introduce, as in GR, the trace-reversed potentials

[TABLE]

In the teleparallel approach to gravity, the sixteen gravitational potentials consist of ten metric potentials familiar from GR and six local Lorentz potentials connected with the local choice of the tetrad system involving three rotations and three boosts. Gravitational potentials are gauge dependent. Under an infinitesimal coordinate transformation, $x^{\mu}\mapsto x^{\prime\mu}=x^{\mu}-\epsilon^{\mu}(x)$, the potentials change to linear order in accordance with $\psi_{\mu\nu}\mapsto\psi^{\prime}_{\mu\nu}=\psi_{\mu\nu}+\epsilon_{\mu,\nu}$. Hence,

[TABLE]

As expected, the linearized gravitational field as well as the corresponding field equations remains invariant under gauge transformations. For instance, it is straightforward to show that the torsion tensor

[TABLE]

and the auxiliary torsion tensor

[TABLE]

do not change under a gauge transformation. To obtain the field equations of linearized NLG, we set $\Lambda=0$ and note that Eq. (33) reduces in the linear regime to

[TABLE]

Here, $T_{\mu\nu}$, $T^{\mu\nu}{}_{,\mu}=0$, is the conserved symmetric energy-momentum tensor of matter. We can write

[TABLE]

where $\mathchoice{\vbox{\hrule height=0.4pt\hbox{\vrule width=0.4pt,height=6.0pt\kern 6.0pt\vrule width=0.4pt}\hrule height=0.4pt}}{\vbox{\hrule height=0.4pt\hbox{\vrule width=0.4pt,height=6.0pt\kern 6.0pt\vrule width=0.4pt}\hrule height=0.4pt}}{\vbox{\hrule height=0.3pt\hbox{\vrule width=0.3pt,height=4.2pt\kern 4.2pt\vrule width=0.3pt}\hrule height=0.3pt}}{\vbox{\hrule height=0.3pt\hbox{\vrule width=0.3pt,height=3.0pt\kern 3.0pt\vrule width=0.3pt}\hrule height=0.3pt}}:=\eta^{\alpha\beta}\partial_{\alpha}\partial_{\beta}$. Furthermore,

[TABLE]

where $X_{\mu\nu\rho}$ is given by Eq. (31) and $\mathcal{K}(x,y)$ reduces to a universal convolution kernel $K(x-y)$ in the linear approximation Mash3c . The linearized field equations of NLG thus take the GR form

[TABLE]

For further discussion, see Refs. Blome:2010xn ; Mash4a ; Mash4b and the references cited therein.

II.2 Kernel of Linearized NLG

A detailed discussion of the causal universal convolution kernel $K(x-y)$ and its reciprocal $R(x-y)$ is contained in Ref. Mash3c . Here, we simply summarize their main properties.

We assume that the convolution kernels $K$ and $R$ are $L^{1}$ and $L^{2}$ functions on spacetime. They are reciprocal of each other and satisfy the reciprocity integral equation

[TABLE]

If $\hat{K}(\xi)$ and $\hat{R}(\xi)$ are Fourier integral transforms of $K(x)$ and $R(x)$, respectively, then, $(1+\hat{K})(1+\hat{R})=1$. Thus, given $\hat{R}(\xi)$, one can in principle determine $K$ from

[TABLE]

provided $1+\hat{R}(\xi)\neq 0$.

We assume that the reciprocal kernel is given by

[TABLE]

where $\nu^{-1}$ is a galactic length scale, $\Theta$ is the Heaviside unit step function such that $\Theta(t)=0$ for $t<0$ and $\Theta(t)=1$ for $t\geq 0$. Moreover, the Newtonian reciprocal kernel $q(\mathbf{x}-\mathbf{y})$ has been determined on the basis of the observational data regarding the rotation curves of spiral galaxies. Two possible forms for $q$ are given by

[TABLE]

and

[TABLE]

where $r=|\mathbf{x}-\mathbf{y}|$ and $\lambda_{0}$, $\mu_{0}$ and $a_{0}$ are constant parameters such that $\lambda_{0}$, the fundamental length scale of NLG, is expected to be $\approx 3$ kpc and $\mu_{0}^{-1}\approx 17$ kpc RaMa . The short-distance nonlocality parameter $a_{0}$ is expected to be much smaller than $\lambda_{0}$. From the solar system data for the orbit of Saturn, one expects approximately that $a_{0}$ is greater than or about the size of the solar system Ior ; DX ; ChMa .

Let us now return to Eq. (49) for the reciprocal kernel and note that

[TABLE]

This has an important implication for gravitational fields that are independent of time. For instance, let $Z(y)$ be a smooth function that is independent of $y^{0}$, then,

[TABLE]

It follows from the reciprocity relation that

[TABLE]

where $\chi$ here is the Newtonian kernel reciprocal to $q$; that is,

[TABLE]

II.3 GEM in Linearized NLG

Consider the stationary gravitational field of a rotating astronomical source in the linear approximation. Equation (45) then takes the form

[TABLE]

Next, from $\partial\chi/\partial x^{i}=-\partial\chi/\partial y^{i}$ and Gauss’s theorem, we find

[TABLE]

It proves convenient to impose the gauge conditions

[TABLE]

which correspond to the transverse gauge condition in the metric part and the vanishing of the local Lorentz potentials. The gauge is not quite fixed, however, since a gauge transformation with $\epsilon_{\mu}=\partial_{\mu}\zeta$ and $\mathchoice{\vbox{\hrule height=0.4pt\hbox{\vrule width=0.4pt,height=6.0pt\kern 6.0pt\vrule width=0.4pt}\hrule height=0.4pt}}{\vbox{\hrule height=0.4pt\hbox{\vrule width=0.4pt,height=6.0pt\kern 6.0pt\vrule width=0.4pt}\hrule height=0.4pt}}{\vbox{\hrule height=0.3pt\hbox{\vrule width=0.3pt,height=4.2pt\kern 4.2pt\vrule width=0.3pt}\hrule height=0.3pt}}{\vbox{\hrule height=0.3pt\hbox{\vrule width=0.3pt,height=3.0pt\kern 3.0pt\vrule width=0.3pt}\hrule height=0.3pt}}\,\zeta=0$ is still possible.

With the imposition of gauge condition (58), the torsion pseudovector vanishes, $\check{C}^{\mu}=0$, by Eqs. (32) and (41); then, the gravitational field equations of linearized NLG for a stationary source take the form

[TABLE]

since $X_{\mu}{}^{\sigma}{}_{\nu}=\mathfrak{C}_{\mu}{}^{\sigma}{}_{\nu}$ and $\partial_{\sigma}\,\mathfrak{C}_{\mu}{}^{\sigma}{}_{\nu}=\,^{0}G_{\mu\nu}$. Equation (44) implies that with the transverse gauge condition, we have $\mathchoice{\vbox{\hrule height=0.4pt\hbox{\vrule width=0.4pt,height=6.0pt\kern 6.0pt\vrule width=0.4pt}\hrule height=0.4pt}}{\vbox{\hrule height=0.4pt\hbox{\vrule width=0.4pt,height=6.0pt\kern 6.0pt\vrule width=0.4pt}\hrule height=0.4pt}}{\vbox{\hrule height=0.3pt\hbox{\vrule width=0.3pt,height=4.2pt\kern 4.2pt\vrule width=0.3pt}\hrule height=0.3pt}}{\vbox{\hrule height=0.3pt\hbox{\vrule width=0.3pt,height=3.0pt\kern 3.0pt\vrule width=0.3pt}\hrule height=0.3pt}}\,\overline{h}_{\mu\nu}=-2\,^{0}G_{\mu\nu}$. It then follows from the temporal independence of gravitational potentials and the reciprocity relation that

[TABLE]

We assume the source consists of slowly moving matter ($|\mathbf{v}|\ll c$) of density $\rho$, pressure $P$ and matter current ${\bf j}=\rho\,{\bf v}$. The matter energy-momentum tensor can thus be written as $T_{00}=\rho c^{2}$, $T_{0i}=-c\,j_{i}$ and $T_{ij}\sim\rho v_{i}v_{j}+P\delta_{ij}$. The corresponding gravitational potentials are $\overline{h}_{00}=-4\Phi/c^{2}$, $\overline{h}_{0i}=-2A_{i}/c^{2}$ and $\overline{h}_{ij}=O(c^{-4})$. In the gravitational potentials, we neglect all terms of $O(c^{-4})$. The static gravitoelectric and gravitomagnetic potentials are thus given by

[TABLE]

and

[TABLE]

respectively. Here, $\rho_{D}(\mathbf{x})$ and $\mathbf{j}_{D}(\mathbf{x})$ are the effective dark matter density and current, respectively. The transverse gauge condition requires that $\bm{\nabla}\cdot\mathbf{A}=0$. On the other hand, $\bm{\nabla}\cdot\mathbf{j}=0$ follows from the energy-momentum conservation law. It follows that the dark matter current is conserved as well; that is, $\bm{\nabla}\cdot\mathbf{j}_{D}=0$.

Let us next introduce the GEM fields $\mathbf{E}_{g}=\bm{\nabla}\Phi$ and $\mathbf{B}_{g}=\bm{\nabla}\times\mathbf{A}$ such that

[TABLE]

[TABLE]

We remark in passing that our GEM conventions are in conformity with the gravitational Larmor theorem Mashhoon:2003ax .

The corresponding GEM metric takes the form

[TABLE]

whose geodesics can be employed to investigate the motion of test particles and null rays in nonlocal GEM. It is possible to show, for instance, the existence of the gravitational analog of the Lorentz force law Mashhoon:2003ax . We note that $\Phi(\mathbf{x})$ is the gravitoelectric potential of nonlocal gravity in the Newtonian regime and has been investigated in some detail Mash3c ; therefore, we concentrate here first on the gravitomagnetic vector potential $\mathbf{A}(\mathbf{x})=O(c^{-1})$.

It is interesting to compare nonlocal GEM with the standard GR treatment Mashhoon:2000he ; Mashhoon:2003ax . In NLG, the steady-state assumption is rather necessary and leads to great simplification. Thus topics such as time-varying gravitomagnetism or gravitational induction that are standard in the GR treatment are beyond the reach of nonlocal GEM. Furthermore, the steady-state requirement limits any further gauge freedom; however, it is possible to shift the magnitude of the the gravitoelectric potential by a constant in metric (65). To this end, let us recall that the remaining gauge freedom is in the form $\epsilon_{\mu}=\partial_{\mu}\zeta$, where $\mathchoice{\vbox{\hrule height=0.4pt\hbox{\vrule width=0.4pt,height=6.0pt\kern 6.0pt\vrule width=0.4pt}\hrule height=0.4pt}}{\vbox{\hrule height=0.4pt\hbox{\vrule width=0.4pt,height=6.0pt\kern 6.0pt\vrule width=0.4pt}\hrule height=0.4pt}}{\vbox{\hrule height=0.3pt\hbox{\vrule width=0.3pt,height=4.2pt\kern 4.2pt\vrule width=0.3pt}\hrule height=0.3pt}}{\vbox{\hrule height=0.3pt\hbox{\vrule width=0.3pt,height=3.0pt\kern 3.0pt\vrule width=0.3pt}\hrule height=0.3pt}}\,\zeta=0$. With $\zeta=-\beta\,(3t^{2}+|\mathbf{x}|^{2})/6$, metric (65) changes to

[TABLE]

where $\beta$ is a constant parameter.

III Gravitomagnetism in Nonlocal GEM

In Eq. (62), the gravitomagnetic vector potential depends on the choice of the reciprocal kernel $q$. To indicate which Newtonian reciprocal kernel is under consideration, we introduce a parameter $\delta$ such that $\delta=1$ for $q_{1}$ and $\delta=2$ for $q_{2}$. Let us write the solution of Eq. (62) in the form

[TABLE]

Using Eq. (62) and the explicit form of the Newtonian reciprocal kernels (50)–(51), we find

[TABLE]

where $r=|\mathbf{r}|$, $\alpha_{0}:=2/(\lambda_{0}\,\mu_{0})\approx 11$ and ${}^{\delta}\mathbb{N}$ is given by

[TABLE]

Here, $E_{1}$ is the exponential integral function given by

[TABLE]

so that for $x:0\to\infty$, $E_{1}(x)$ is a positive monotonically decreasing function that diverges as $-\ln x$ near $x=0$ and falls off exponentially as $x\to\infty$. Moreover, we have introduced a dimensionless quantity $\varsigma$ such that $\varsigma:=a_{0}\,\mu_{0}<1$. For the exterior of the Earth, we may assume $r$ is small compared to an astronomical unit and therefore $r\ll a_{0}$, as $a_{0}$ is about the size of the solar system. For $r\ll a_{0}$, we have $\mu_{0}\,r\ll\varsigma$; then, the Taylor expansion of $E_{1}(\varsigma+\mu_{0}\,r)$ about $E_{1}(\varsigma)$ and repeated differentiation of Eq. (70) result in

[TABLE]

Putting these results together and neglecting terms of $O\left(r^{3}/a_{0}^{3}\right)$, we find

[TABLE]

Assuming $|\mathbf{x}|>|\mathbf{y}|$, which is appropriate for the exterior of the source, and expanding $|\mathbf{x}-\mathbf{y}|$ to first order in $|\mathbf{y}|/|\mathbf{x}|$, we get

[TABLE]

and

[TABLE]

Let the compact gravitational source reside in the interior of a finite closed spatial domain $\mathcal{D}$ that completely surrounds the source. This means that $\mathbf{j}$ vanishes on the surface of $\mathcal{D}$ and beyond. Then, the conservation of matter current implies

[TABLE]

where $f(\mathbf{y})$ is a smooth function. Applying Gauss’s theorem and setting the integral on $\partial\,\mathcal{D}$ equal to zero, we get

[TABLE]

For $f(\mathbf{y})=1$ and $f(\mathbf{y})=y^{i}$, we find the following relations

[TABLE]

respectively. Let

[TABLE]

be the total proper angular momentum of the gravitational source. Then, it is straightforward to show using Eq. (77) that

[TABLE]

It then follows from these results that the gravitomagnetic vector potential is given by

[TABLE]

where the relevant nonlocality length scale $L_{N}$ is given by

[TABLE]

The nonlocal contribution to $\mathbf{A}$ at the level of approximation under consideration is nonzero for $q_{1}$ but vanishes for $q_{2}$. The length scale $L_{N}\gtrsim 1$ pc, so that the nonlocal contribution to $\mathbf{A}$ in the exterior of the Earth is relatively quite small and less than about $10^{-10}$ of the standard GR value.

Finally, the gravitomagnetic field can be calculated from Eq. (80) and the result is

[TABLE]

The gravitomagnetic field of the Earth has been directly measured via the GP-B experiment and the GR prediction has been verified to about 19% Francis . The nonlocal contribution to the gravitomagnetic field of the Earth is at most ten orders of magnitude smaller than the GR value and is thus beyond current measurement capabilities for the foreseeable future. A similar estimate holds for nonlocal gravitomagnetic effects in the motion of the Moon. In connection with the lunar laser ranging experiment, we note that the main relativistic effects in the motion of the Moon are due to the gravitational field of the Sun and have been calculated in Refs. MaTh1 ; MaTh2 . The Earth-Moon system with its orbital angular momentum acts as an extended gyroscope in the gravitomagnetic field of the Sun. The nonlocal modification of this field is given by Eq. (82) and the corresponding nonlocal gravitomagnetic effects in the motion of the Moon are then about ten orders of magnitude smaller than the GR predictions as well. Another consequence of the existence of the gravitomagnetic field is the Lense-Thirring effect, see Refs. MaHeT ; ILRC ; Renz and the references cited therein.

III.1 Nonlocal Contributions to the Metric

The spacetime metric (66) in our nonlocal GEM contains gravitoelectric and gravitomagnetic potentials. The latter is given by Eq. (80). It is therefore necessary to find the corresponding gravitoelectric potential ${{}^{\delta}}\Phi(\mathbf{x})$, which is given by

[TABLE]

To simplify matters, we assume that the gravitational source has a spherically symmetric matter distribution. This means that $\rho(\mathbf{y})=\rho(|\mathbf{y}|)$; then, we go through essentially the same steps as in Eqs. (68)–(74), except that

[TABLE]

as a consequence of spherical symmetry for the matter distribution. Therefore, Eq. (72) impies

[TABLE]

and

[TABLE]

Here, $M$ is the mass of the spherical rotating source in our linear approximation scheme, namely,

[TABLE]

We note that Eq. (85) here is consistent with Eqs. (8.39) and (8.40) of Ref. Mash3c . With gravitoelectric potential (85) and gravitomagnetic potential (80), the GEM metric (65) can now be used consistently in investigating nonlocal effects in GEM.

III.2 Gravitomagnetic Clock Effect in NLG

There is a special temporal structure around a rotating mass that is best expressed via the gravitomagnetic clock effect CoMa ; MGT ; MGL ; MIL . To illustrate this effect in NLG, let us assume that the gravitational source rotates about the $z$ axis, $\mathbf{J}=J\,\hat{\mathbf{z}}$, and write the GEM metric in the corresponding spherical polar coordinates. Under the transformation $x^{\mu}\mapsto(t,r,\theta,\phi)$, metric (65) takes the form

[TABLE]

where $g_{tt}:=-1-2\,\Phi$, $g_{t\phi}:=-2\,r\,\sin\theta\,A_{\phi}$ and we have neglected in our GEM approach the contribution of $\Phi$ to the spatial part of the metric. In the local theory (GR), we have $\Phi=-GM/r$ and $A_{\phi}=(GJ/r^{2})\sin\theta$. These potentials change in our nonlocal approach as follows

[TABLE]

where

[TABLE]

We are interested in the nonlocal modification of Keplerian periods of the equatorial circular orbits in this spacetime.

The geodesic equation for the radial coordinate takes the form

[TABLE]

where $\tau$ is the proper time. This equation can be solved for $r=$ constant and $\theta=\pi/2$. The solution in the linear approximation under consideration is given by

[TABLE]

Indeed, for $\theta=\pi/2$, we have

[TABLE]

It follows from a detailed analysis that, as expected, deviations exist from the standard GR results for $\delta=1,2$. For an equatorial circular orbit with Keplerian frequency $\omega_{K}=(GM/r^{3})^{1/2}$ and Keplerian period $T_{K}=2\pi/\omega_{K}$, we find for the periods of co-rotating (+) and counter-rotating (-) orbits in terms of coordinate time

[TABLE]

and in terms of proper time

[TABLE]

Here, we work to linear order in perturbation quantities and the nonlocal contributions are given by terms proportional to $r^{2}/L_{N}^{2}$, $\Delta_{M}$ and $\Delta_{J}$, where

[TABLE]

It is interesting to note that the prograde period is longer than the retrograde period, namely,

[TABLE]

In GR, the gravitomagnetic clock effect for circular equatorial orbits around the Earth is given by $\tau_{+}-\tau_{-}\approx 2\times 10^{-7}$ sec. This prediction of GR has not yet been verified by observation. The GR effect is indeed rather difficult to measure since the Keplerian period of a near-Earth orbit increases by about $2\times 10^{-7}$ sec when the orbital radius is increased by $0.015$ cm, see Refs. Tar1 ; Tar2 ; LGM ; Ior1 ; Ior2 ; ILM ; IoLi ; LiIoM ; Hackmann:2014aga and the references cited therein. The magnitude of the nonlocal contribution to the gravitomagnetic clock effect for the Earth is smaller than about $10^{-10}$ of the GR value.

III.3 Gravitational Larmor Theorem in NLG

In classical electrodynamics, Larmor’s theorem establishes a local relation between the motion of a charged test particle in an electromagnetic field and its motion in the absence of the field, but in an accelerated system of reference. The gravitational version of this theorem is essentially Einstein’s principle of equivalence expressed within the GEM framework Mashhoon:2003ax ; Mashhoon:2000he . It is useful to point out that the theorem extends to nonlocal GEM as well.

Let us imagine an accelerated observer following a world line $\bar{X}^{\mu}(\tau)$ in Minkowski spacetime. Here, $\tau$ is the observer’s proper time. The observer carries an orthonormal tetrad frame $\lambda^{\mu}{}_{\hat{\alpha}}(\tau)$ along its path such that

[TABLE]

where $\Psi_{\hat{\alpha}\hat{\beta}}=-\Psi_{\hat{\beta}\hat{\alpha}}$ is the antisymmetric acceleration tensor. In analogy with the electromagnetic field tensor, we can decompose $\Psi_{\hat{\alpha}\hat{\beta}}$ into its “electric” and “magnetic” parts, namely, $\Psi_{\hat{0}\hat{i}}=\gamma_{\hat{i}}$ and $\Psi_{\hat{i}\hat{j}}=\epsilon_{\hat{i}\hat{j}\hat{k}}\,\omega^{\hat{k}}$. Here, $\bm{\gamma}$ and $\bm{\omega}$ represent the invariant translational and rotational accelerations of the observer, respectively. Let us now introduce a geodesic coordinate system in the neighborhood of the accelerated observer. At a given proper time $\tau$, the straight spacelike geodesics normal to $\bar{X}^{\mu}(\tau)$ form a Euclidean hyperplane. An event on this hyperplane with inertial coordinates $X^{\mu}$ will be assigned geodesic (Fermi) coordinates $x^{\hat{\mu}}$ such that $x^{\hat{0}}=\tau$ and $X^{\mu}-\bar{X}^{\mu}(\tau)=x^{\hat{i}}\,\lambda^{\mu}{}_{\hat{i}}(\tau)$. With these transformations, $dS^{2}=\eta_{\mu\nu}\,dX^{\mu}\,dX^{\nu}$ becomes $dS^{2}=g_{\hat{\mu}\hat{\nu}}\,dx^{\hat{\mu}}\,dx^{\hat{\nu}}$, where

[TABLE]

A detailed discussion of these local coordinates and their admissibility is contained in Ref. Mash3c . In general, $\bm{\gamma}$ and $\bm{\omega}$ are functions of proper time $\tau$. However, in the present context of steady-state GEM, we assume that these accelerations are constants and do not vary with proper time.

A comparison of this flat metric at the linear order with metric (66) once we neglect its spatial curvature reveals that an accelerated observer in Minkowski spacetime is locally equivalent to an observer in a GEM field provided $\Phi+\beta=\bm{\gamma}\cdot\mathbf{x}$ for a suitable choice of the constant $\beta$ and $-2\,\mathbf{A}=\bm{\omega}\times\mathbf{x}$, which means that $\mathbf{E}_{g}=\bm{\gamma}$ and $\mathbf{B}_{g}=-\bm{\omega}$, respectively. These GEM fields contain nonlocal effects; in this way, the gravitational Larmor theorem has been extended to the nonlocal regime.

An interesting application of the gravitational Larmor theorem involves the interaction of spin with the gravitational field. The coupling of intrinsic spin with the gravitomagnetic field has been discussed extensively and a brief review of the subject is contained in Ref. Mashhoon:2003ax . The effect is related to spin-rotation coupling via the gravitational Larmor theorem. The spin-rotation coupling for neutrons has recently been measured via neutron interferometry DDSH ; DDKWLSH . The extension of the gravitational Larmor theorem to the nonlocal regime means that spin-gravity coupling can likewise be extended to the nonlocal regime.

IV Gravitational Energy-Momentum Tensor

The traceless gravitational energy-momentum tensor $\mathcal{T}_{\mu\nu}$ of NLG is given by Eq. (34). Let us first compute the local part of this tensor $\mathbb{T}_{\mu\nu}$, which is traceless as well, for the GEM case. To this end, we write Eq. (15) in the form

[TABLE]

and express the components of the torsion tensor in terms of the GEM potentials. That is,

[TABLE]

and

[TABLE]

It follows that

[TABLE]

where we have used the relation

[TABLE]

It is now possible to compute the components of the traceless energy-momentum tensor, which are

[TABLE]

[TABLE]

and

[TABLE]

where

[TABLE]

These local results must be supplemented with nonlocal terms. That is, we must go back to Eq. (34) and compute $\mathcal{T}_{\mu\nu}$, which contains nonlocal terms of the form

[TABLE]

where $\chi$ is the kernel of NLG theory in the Newtonian regime Mash3c . The explicit calculation of this kernel is rather complicated and is beyond the scope of this paper.

It is interesting to compare and contrast the local Eqs. (106)–(109) with those obtained via the Landau-Lifshitz pseudotensor $t_{\mu\nu}$ of GR L+L within the standard GEM framework Mashhoon:2003ax ; Mashhoon:2000he . To this end, it is necessary to assume a steady-state GR configuration (i.e., $\partial\Phi/\partial t=0$ and $\partial\mathbf{A}/\partial t=0$). Then,

[TABLE]

The similarity between these different gravitational results and the corresponding electromagnetic ones is noteworthy. In particular, imagine a steady-state configuration involving a rotating astronomical source with mass $M$ and angular momentum $\mathbf{J}=J\,\hat{\mathbf{z}}$. Then, it follows from the gravitational Poynting vector that there is a steady circulation of gravitational energy in the same sense as the rotation of the source with an azimuthal flow speed given in spherical polar coordinates by

[TABLE]

The proportionality constant depends on the underlying theory of gravitation Mashhoon:2003ax ; Mashhoon:2000he .

V Discussion

We have developed gravitoelectromagnetism (GEM) within the framework of nonlocal gravity (NLG). Except for the trivial solution of field equations involving flat spacetime, NLG has no other known exact solution at present. We must therefore resort to the linearized theory, where GEM is possible for steady-state configurations. We have examined the nonlocal GEM corrections to the stationary gravitational field of an isolated rotating mass in the weak-field and slow-motion approximations. Due to the existence of galactic length scales in NLG, the nonlocal GEM effects around the Earth or the Sun turn out to be at most about ten orders of magnitude smaller than the corresponding GR effects.