Arxiv:1804.11106V2 [Gr-Qc] 2 Jun 2018 Oe Yscinii Nwihw Ou Nteproperties the on Focus Fol- We Is Which This in III, Equation

Home , Geodesic deviation

arXiv:1804.11106v2 [gr-qc] 2 Jun 2018 oe yscinII nwihw ou nteproperties the on focus fol- we is which This in III, equation. section deviation by lowed generalized exact derivation the the in needed of concepts the introduce briefly we well [36]. as theory Poincaré [35] gauge theory to par- Einstein-Cartan as in to background; apply valid Riemann-Cartan they are a ticular, results in Our theory any on lines. events for world links (adjacent) which general vector two the connecting space- describes the equation of for This dynamics allowing torsion. with background, endowed times Riemann-Cartan a in to seems devi- geometries overdue. the be post-Riemannian of to extension means further equation to by a ation coupling theories bodies, gravitational test natural of structured a of tests is possible in there of degrees Therefore, view – intrinsic spacetime. of spin with features post-Riemannian like bodies the – test freedom that of [32–34] of known description well the already is in it particular, frame- In geometrical work. Riemannian the reach situa- beyond [29–31] physical significantly theories background. many gravitational Riemannian modern in results several justified a tions, the are in theories However, such theories While to lines. limited world deviation the were general of applicability to of findings range equation the Our al- extended for [5] generalizations. [6–28] in also and see derivations equation, ternative deviation the of versions be test can neutral structureless and of description assumptions bodies. the equa- certain for the This only under of used means holds [1–4]. by only equation described tion (Jacobi) is deviation bodies test geodesic of motion tive ∗ † [email protected] ikpezedzr.n-rmnd;http://puetzfeld. [email protected]; h tutr fteppri sflos nscinII, section In follows: as is paper the of structure The equation deviation generalized a derive we work this In generalized out worked have we [5], work previous a In rela- the Relativity, General of theory the Within nvriyo rmn etro ple pc ehooyan Technology Space Applied of Center Bremen, of University ewrs eito qain prxmto ehd;Equ methods; Approximation equation; Deviation 04.25.-g Keywords: 04.20.Cv; ge 04.50.-h; the numbers: beyond backgrounPACS go Riemann-Cartan which a theories in t gravitational on of theory events class any links large for which valid vector is connecting result the of dynamics the edrv eeaie eito qaini Riemann-Cart in equation deviation generalized a derive We .INTRODUCTION I. eito qaini imn-atnspacetime Riemann-Cartan in equation Deviation usa cdm fSine,BTlky 2 111Moscow, 115191 52, B.Tulskaya Sciences, of Academy Russian hoeia hsc aoaoy ula aeyInstitute Safety Nuclear Laboratory, Physics Theoretical org Dtd ue5 2018) 5, June (Dated: uiN Obukhov N. Yuri ikPuetzfeld Dirk eut fordrvto r ie napni B. appendix intermediate in and given are details derivation Some our I. of table results in applica- and summarized possible are results A conventions their the appendix and of of notations discussion outlook Our a an tions. with with VI and section obtained in conclude We paper equation. our deviation the the un- and of discuss interpretation value the the operational we on affect V theory, depending ver- section gravity choices, derlying approximate In coordinate expanded equation. different an in deviation how the at applied of arrive then sion are to results IV Riemann- These a section in autoparallels background. on based Cartan function world a of otegnrlzddvaineuto sdrvdi [5]. curves in general derived two as compare to equation want deviation We generalized the to na rirr pctm aiod Here manifold. spacetime arbitrary an in enwcnettopoints case, Riemannian two the connect to time contrast now In proper we curves. the given necessarily the not on i.e. parameters, general aallt tefwt epc otecneto nthe lines. on geodesic autoparallel connection with space the coincide Riemannian curves a to In transported respect manifold. is with vector spacetime itself velocity the to which parallel autoparallel along An curve a unique). (we is is points autoparallel two this the that joining assume autoparallel the by curves two h ol ucin eoe ttepoint of the derivative at covariant denoted the function, is world points the two vec- the connecting the between to tor object closest the conceptually and be .Floigorcnetos h eeec uv ilbe will curve reference the Y conventions, our Following 1. ( e sbifl eaiuaeterlvn tp htlead that steps relevant the recapitulate briefly us Let ln h uoaallw aetewrdfunction world the have we autoparallel the Along t I OL UCINADDEVIATION AND FUNCTION WORLD II. irgaiy(AM,239Bee,Germany Bremen, 28359 (ZARM), Microgravity d ,adw en h eeaie oncigvco to vector connecting generalized the define we and ), toso oin imn-atnspacetime Riemann-Cartan motion; of ations ∗ † ea eaiitcframework. relativistic neral nsaeie h qaindescribes equation The spacetime. an ognrlajcn ol ie.Our lines. world adjacent general wo ,i atclr ti plcbet a to applicable is it particular, in d, EQUATION η , y Russia := − x σ ∈ y X and y Y y by ( ∈ t t and ) σ and Y y f fig. cf. , nthe on t ˜ X are (1) ( σ t ˜ ) , 2

and inserted into (3)

2 D2 dt˜ ηy1 = −σy1 uy2 uy3 − σy1 ay2 − σy1 a˜x2 dt2 y2y3 y2 x2 dt Dσy4 −2σy1 uy2 Kx3 uy4 − Hx3 y2x3 y4 y4 dt Dσy4 −σy1 Kx2 uy4 − Hx2 x2x3 y4 y4 dt Dσy5 × Kx3 uy5 − Hx3 y5 y5 dt 2˜ y3 y1 dt d t x2 y3 x2 Dσ −σ x K y u − H y . (5) 2 dt˜dt2 3 3 dt Note that we may determine the factor dt/dt˜ by requiring that the velocity along the curve X is normalized, i.e. x u˜ u˜x = 1, in which case we have

FIG. 1. Sketch of the two arbitrarily parametrized world lines dt˜ Dσy2 =u ˜ Kx1 uy2 − u˜ Hx1 . (6) Y (t) and X(t˜) and the (dashed) autoparallel connecting two dt x1 y2 x1 y2 dt points on these world lines. The generalized deviation vector along the reference world line Y is denoted by ηy. Equation (5) is the exact generalized deviation equation, it is completely general and can be viewed as the extension of the standard geodesic deviation (Jacobi) equation to any order. In particular, it allows for a comparison of Taking its covariant total derivative, we have two general, i.e. not necessarily geodetic or autoparallel, world lines in spacetime. ˜ D y1 y1 y2 y1 x2 dt η = −σ y2 u − σ x2 u˜ , (2) dt dt III. WORLD FUNCTION IN RIEMANN-CARTAN SPACETIME where the velocities along the two curves Y and X are deﬁned as uy := dY y/dt, andu ˜x := dXx/dt˜. Denoting In this section we work out the basic properties of a the accelerations by ay := Duy/dt, anda ˜x := Du˜x/dt˜, world function σ based on autoparallels in a Riemann- the second derivative becomes Cartan spacetime, which, in contrast to a Riemannian spacetime, is endowed with an asymmetric connection Γ c. Relevant references that contain some results in a D2 dt˜ ab y1 − y1 y2 y3 − y1 y2 x3 Riemann-Cartan context are [37–44]. 2 η = σ y2y3 u u 2σ y2x3 u u˜ dt dt For a world function σ based on autoparallels, we have 2 dt˜ the following basic relations in the case of spacetimes −σy1 ay2 − σy1 u˜x2 u˜x3 y2 x2x3 dt with asymmetric connections: 2 dt˜ d2t˜ x y − y1 x2 − y1 x2 σ σx = σ σy =2σ, (7) σ x2 a˜ σ x2 u˜ 2 . (3) dt dt x2 x1 x1 σ σx2 = σ , (8) x3 σx1x2 − σx2x1 = Tx1x2 ∂x3 σ. (9) The generalized deviation equation is obtained from (3) Note, in particular, the change in (9) due to the presence by expressing all quantities therein along the reference x3 −1 of the spacetime torsion Tx1x2 , which leads to σx1x2 =6 σ word line Y . With the help of the inverse – obtained σx x , in contrast to the symmetric Riemannian case, in −1 −1 2 1 y1 x y1 x1 y x1 s 1 from σ xσ y2 = δ y2 and σ yσ x2 = δ x2 – and by which σx1x2 = σx2x1 holds. −1 x1 x1 In many calculations the limiting behavior of a biten- deﬁning the Jacobi propagators H y2 := − σ y2 and −1 sor B...(x, y) as x approaches the references point y is x1 x1 y1 K y2 := − σ y1 σ y2 , the velocity along X may be expressed as

y1 1 We use “s” to indicate relations which only hold for symmetric x3 x3 y2 dt x3 Dσ dt u˜ = K y u − H y , (4) connections and denote Riemannian objects by the overbar. 2 dt˜ 1 dt dt˜ 3 required. This so-called coincidence limit of a bitensor (third) order the presence of the torsion leads to B (x, y) is a tensor ... 1 [σ ]= (T + T + T )= K , x1x2x3 2 y1y3y2 y2y3y1 y1y2y3 y2y1y3 [B...] = lim B...(x, y), (10) x→y (16) where in the last line we made use of the contortion K, at y and will be denoted by square brackets. In partic- cf. also appendix A for an overview of the geometrical ular, for a bitensor B with arbitrary indices at diﬀerent quantities.2 With the help of (11), we obtain for the points (here just denoted by dots), we have the rule [45] other combinations with three indices

[σx1x2y3 ]= −[σx1x2x3 ] = [σy1x2y3 ] = [σx2x3x1 ]= Ky2y3y1 , [B...];y = [B...;y] + [B...;x] . (11) [σx1y2y3 ]= −[σx2x3x1 ] = [σx1x3x2 ] = [σy1x2x3 ]= Ky3y1y2 ,

We collect the following useful identities for the world [σy1y2y3 ]= −[σx3x2x1 ]= Ky2y1y3 , function σ: [σx1y2x3 ]= −[σx1x3x2 ]= Ky3y2y1 ,

[σy1y2x3 ] = [σx3x2x1 ]= Ky2y3y1 . (17) [σ] = [σ ] = [σ ]=0, (12) x y The non-vanishing of these limits leads to added com- [σx1x2 ] = [σy1y2 ]= gy1y2 , (13) plexity in subsequent calculations compared to the Rie-

[σx1y2 ] = [σy1x2 ]= −gy1y2 , (14) mannian case. At the fourth order we have [σx3x1x2 ] + [σx2x1x3 ]=0. (15) y y y Ky1 y2 Ky3yy4 + Ky1 y3 Ky2yy4 + Ky1 y4 Ky2yy4 Note that up to the second covariant derivative the coin- +[σx x x x ] + [σx x x x ] + [σx x x x ]=0, (18) cidence limits of the world function match those in space- 4 1 2 3 3 1 2 4 2 1 3 4 times with symmetric connections. However, at the next and in particular

1 1 1 [σ ] = ∇ (K + K )+ ∇ (3K − K )+ ∇ (3K − K )+ π , x1x2x3x4 3 y1 y3y2y4 y4y2y3 3 y3 y2y1y4 y1y2y4 3 y4 y2y1y3 y1y2y3 y1y2y3y4 (19) 1 1 1 [σ ] = − ∇ (K + K ) − ∇ (3K − K )+ ∇ K − π , (20) x1x2x3y4 3 y1 y3y2y4 y4y2y3 3 y3 y2y1y4 y1y2y4 3 y4 y1y2y3 y1y2y3y4 1 1 1 [σ ] = ∇ (K + K ) − ∇ K − ∇ K + π , (21) x1x2y3y4 3 y1 y4y2y3 y3y2y4 3 y4 y1y2y3 3 y3 y1y2y4 y1y2y4y3 1 1 1 [σ ] = − ∇ (K + K )+ ∇ K + ∇ K + ∇ K − π , (22) x1y2y3y4 3 y1 y3y4y2 y2y4y3 3 y3 y1y4y2 3 y2 y1y4y3 y4 y3y1y2 y1y4y3y2 1 1 1 [σ ] = ∇ (−2K + K ) − ∇ K − ∇ K − ∇ K + π , (23) y1y2y3y4 3 y4 y2y3y1 y1y3y2 3 y2 y4y3y1 3 y1 y4y3y2 y3 y2y4y1 y4y3y2y1 1 π := K y (K + K ) − K y (K + K ) − K y (K + K ) y1y2y3y4 3 y1y2 y3y4y y4y3y y1y3 y4y2y yy2y4 y1y4 y3y2y yy2y3 h y y y −3Ky2y1 Ky3y4y + Ky3y1 Kyy2y4 + Ky4y1 Kyy2y3 + Ry1y3y2y4 + Ry1y4y2y3 . (24) i

y Again, we note the added complexity compared to terms of a parallely propagated tetrad e(a), which in s y the Riemannian case, in which we have [σx1x2x3x4 ] = turn allows for the transport of objects, i.e. V = 1 y x y1y2 y1 y2 x1x2 3 Ry2y4y1y3 + Ry3y2y1y4 at the fourth order. In par- g xV , V = g x1 g x2 V , etc., along an au- ticular, we observe the occurrence of derivatives of the toparallel: contortion in (19)-(23). Finally, let us collect the basic properties of the so- y1 x y1 x1 y x1 g xg y2 = δ y2 , g yg x2 = δ x2 , (25) y y (a) called parallel propagator g x := e ex , deﬁned in x x1 y x1 (a) σ ∇xg y1 = σ ∇yg y1 =0, x y1 y y1 σ ∇xg x1 = σ ∇yg x1 =0, (26) y x σx = −g xσy, σy = −g yσx. (27) 2 K The contortion y2y1y3 should not be confused with the Jacobi Kx propagator y. Note, in particular, the coincidence limits of its deriva- 4 tives y3 1 σy1y2x3 = gx3 Ky2y3y1 − ∇y3 (Ky2y4y1 + Ky1y4y2 ) ( 3" x0 y0 [g y1 ]= δ y1 , (28) − ∇y2 Ky3y4y1 − ∇y1 Ky3y4y2 x0 x0 [g y1;x2 ] = [g y1;y2 ]=0, (29) x0 x0 x0 y4 2 [g y ;x x ]= − [g y ;x y ] = [g y ;x x ] +3πy3y4y2y1 σ + O σ , (38) 1 2 3 1 2 3 1 2 3 # ) 1 x0 y0 = − [g y1;y2y3 ]= R y1y2y3 . (30) 2 y2 y3 σy1x2x3 = gx2 gx3 Ky3y1y2 ( In the next section we will derive an expanded approx- 1 imate version of the deviation equation. For this we ﬁrst + ∇y2 (Ky4y3y1 + Ky1y3y4 ) work out the expanded version of quantities around the 3" reference world line Y . In particular, we make use of the − ∇y4 (Ky2y3y1 +3Ky3y1y2 ) covariant expansion technique [45, 46] on the basis of the autoparallel world function. y4 2 − ∇y1 Ky2y3y4 +3πy2y3y4y1 σ + O σ .(39) # ) The Jacobi propagators are approximated as

IV. EXPANDED RIEMANN-CARTAN x1 x1 x1 y3 2 H y2 = g y2 + Ky3y2 σ + O σ , (40) DEVIATION EQUATION x1 x1 x1 x1 y3 2 K y2 = g y2 + (Ky2 y3 + Ky3y2 )σ + O σ , (41) For a general bitensor B... with a given index structure, we have the following general expansion, up to the third which in turn allows for an expansion of the recurring order (in powers of σy): term entering (5):

′ y2 ′ y yn+1 Dσ Dσ x1 y2 x1 x1 ′ y By1...yn = Ay1...yn + Ay1...yn+1 σ K y2 u − H y2 = g y u − dt " dt 1 yn+1 yn+2 3 + Ay1...yn+1yn+2 σ σ + O σ , (31) 2 ′ ′ y y y2 y3 2 + Ky y + Ky y u σ + O σ . (42) Ay ...y := [By ...y ] , (32) 2 3 3 2 1 n 1 n # Ay1...yn+1 := By1...yn;yn+1 − Ay1...yn;yn+1 , (33)

y0 Ay1...yn+2 := By1...yn;yn+1yn+2 − Ay1...yny0 σ yn+1yn+2 A. Synchronous parametrization −A − 2A .(34) y1...yn;yn+1yn+2 y1...yn(yn+1;yn+2) Before writing down the expanded version of the gen- With the help of (31) we are able to iteratively expand eralized deviation equation, we will simplify the latter any bitensor to any order, provided the coincidence lim- by choosing a proper parametrization of the neighboring its entering the expansion coeﬃcients can be calculated. curves. The factors with the derivatives of the parame- The expansion for bitensors with mixed index structure ters t and t˜ appear in (5) due to the non-synchronous can be obtained from transporting the indices in (31) by parametrization of the two curves. It is possible to means of the parallel propagator. make things simpler by introducing the synchronization To develop an approximate form of the generalized de- of parametrization. Namely, we start by rewriting the viation equation (5) up to the second order, we need the velocity as following expansions – note that we give some explicit dY y dt˜dY y intermediate results in appendix B – of the derivatives of uy = = . (43) the world function: dt dt dt˜ That is, we now parametrize the position on the ﬁrst y3 2 σy1y2 = gy1y2 + Ky2y1y3 σ + O σ , (35) curve by the same variable t˜ that is used on the second y y3 2 curve. Accordingly, we denote σy1x2 = −gy1x2 + gx2 Ky3yy1 σ + O σ , (36) dY y 1 uy = . (44) σy1y2y3 = Ky2y1y3 + ∇y4 (Ky2y3y1 + Ky1y3y2 ) e ˜ 3" dt

− ∇y2 Ky4y3y1 − ∇y1 Ky4y3y2 − 3∇y3 Ky2y4y1 By diﬀerentiation, we then derive

2 2 y4 2 y d t˜ y dt˜ y +3πy4y3y2y1 σ + O σ , (37) a = u + a , (45) # dt2 e dt e 5 where into a simpler form

2 D y y y y x y y y 2 y 1 1 2 1 2 1 2 3 D D Y η = −σ y2 a − σ x2 a˜ − σ y2y3 u u ay = uy = . (46) dt˜2 e e e 2 e dt˜e dt˜ y4 y1 y2 x3 y4 x3 Dσ −2σ y x u K y u − H y 2 3 e 4 e 4 dt˜ Analogously, we derive for the derivative of the deviation y4 y1 x2 y4 x2 Dσ vector −σ x2x3 K y4 u − H y4 e dt˜ y5 x3 y5 x3 Dσ 2 y 2 y 2 2 y × K y5 u − H y5 . (48) D η d t˜Dη dt˜ D η e dt˜ = + . (47) dt2 dt2 dt˜ dt dt˜2 B. Explicit expansion of the deviation equation Now everything is synchronous in the sense that both ˜ curves are parametrized by t. Substituting the expansions (35)-(42) into (48), we ob- As a result, the exact deviation equation (5) is recast tain the ﬁnal result

2 y3 D y1 y1 y1 y1 y2 Dη y1 y2 y1 y2 y1 y2 y3 y4 2 η =˜a − a + Ty y u − Ky y a − Ky y a˜ + ∆ y y y u u η + O σ , (49) dt˜2 e 2 3 e dt˜ 2 4 e 4 2 2 3 4 e e where we introduced the abbreviation ′ ′ y ∆y1y2y3y4 := 2πy3y4y2y1 − πy4y3y2y1 − πy2y3y4y1 + Ty y2y1 Ty4y3 − 2∇y2 K(y1y3)y4 + ∇y1 Ky2y3y4 − ∇y4 Ky2y3y1 . (50)

It should be understood that the last expression is con- curring covariant total derivatives need to be rewritten tracted with uy2 uy3 and hence the symmetrization is nat- and an appropriate coordinate choice needs to be made. e e The lhs of the deviation equation takes the form: urally imposed on the indices (y2y3). Equation (49) allows for the comparison of two general 2 D ηa b ◦◦ b d ◦ b c d world lines in Riemann-Cartan spacetime, which are not =u ˙ ∇bηa+ η a −2u Γba ηd −u u Γcb ∂dηa dt2 necessarily geodetic or autoparallel. It therefore repre- − b c e − d e − d e sents the generalization of the deviation equation derived u u ηe ∂cΓba Γcb Γda Γca Γbd . in [5, (35)]. (53) ◦ Here we used η a := dηa/dt for the standard total deriva- C. Riemannian case tive. Observe that the ﬁrst term on the rhs vanishes in the A great simpliﬁcation is achieved in a Riemannian case of autoparallel curves (u ˙ a := Dua/dt = 0). Also background, when note the symmetrization of the connection imposed by the velocities in some terms. ∆ =2π − π − π y1y2y3y4 y3y4y2y1 y4y3y2y1 y2y3y4y1 Rewriting the connection in terms of the contortion

= Ry1y3y2y4 , (51) and switching to normal coordinates [49–55] along the and (49) is reduced to world line, which we assume to be an autoparallel, yields 2 y 2 ◦◦ ◦ D y1 s y1 y1 1 y2 y3 y4 2 D ηa |Y b d b c d η =ã − a − R y y y u u η + O σ . = η a +2u Kba ηd +u u Kcb ∂dηa dt˜2 e 2 3 4 e e dt2 (52) 2 +ubucη ∂ K e − R e + K dK e − K dK e . e c ba 3 c(ba) cb da ca bd Along geodesic curves, this equation is further reduced to the well known geodesic deviation (Jacobi) equation. (54) Note the appearance of a term containing the partial (not V. CHOICE OF COORDINATES ordinary total) derivative of the deviation vector, in contrast to the Riemannian case. The first term in the second line may be rewritten as an To utilize the deviation equation for measurements or ◦ b c e b e in a gravitational compass setup [5, 45, 47, 48], the oc- ordinary total derivative, i.e. u u ηe∂cKba = u ηeKba , 6 but this is still inconvenient when recalling the compass consequently, the correct type of test body. Taking the equation, which will contain terms with covariant deriva- example of a theory with spin-curvature coupling, like tives of the contortion. Einstein-Cartan theory, this would eventually lead to the usage of test bodies with vanishing spin – since those still move on standard geodesics, and therefore lead to A. Operational interpretation an identical procedure as in the general relativistic case, i.e. one adopts coordinates in which the motion of those At this point, some thoughts about the operational test bodies becomes rectilinear. interpretation of the coordinate choice are in order. In particular, it should be stressed that so far we did not specify any physical theory in which the deviation equation (48) should be applied. Stated the other way round, the derived deviation equation is of completely geomet- VI. CONCLUSIONS & OUTLOOK rical nature, i.e. it describes the change of the deviation vector between points on two general curves in Riemann- Cartan spacetime. In this work we investigated the generalization of the From the mathematical perspective, the choice of co- deviation equation in a Riemann-Cartan geometry. As a ordinates should be solely guided by the simplicity of the novel technical result, we have developed Synge’s world resulting equation. In this sense, our previous choice of function approach in the non-Riemannian spacetime with normal coordinates appears to be appropriate. But what curvature and torsion. Our expanded version of the de- about the physical interpretation or, better, the opera- viation equation (48) can be directly compared to re- tional realization of such coordinates? sult in the Riemannian context [5]. The generalization Let us recall the coordinate choice in General Relativ- should serve as a foundation for the test of gravitational ity in a Riemannian background. In this case, normal theories that make use of post-Riemannian geometrical coordinates also have a clear operational meaning, which structures. is related to the motion of structureless test bodies in As we have discussed in detail, the operational us- General Relativity. As is well known, such test bodies ability of the Riemann-Cartan deviation equation differs move along the geodesic equation. In other words, we from the one in a general relativistic context, which was could – at least in principle – identify a normal coordi- also noticed quite early in [29]. In particular, it remains nate system by the local observation of test bodies. If to be shown which additional concepts and assumptions other external forces are absent, normal coordinates will 3 are needed in order to fully realize a gravitational com- locally lead to straight line motion of test bodies. In pass [5, 45, 47, 48] in a Riemann-Cartan background. In this sense, there is a clear operational procedure for the contrast to the Riemannian case, an algebraic realization realization of normal coordinates. of a gravitational compass on the basis of the deviation However, now we are in a more general situation, since equation is out of the question due to the appearance of we have not yet specified which gravitational theory we derivatives of the torsion even at the lowest orders. are considering in the geometrical Riemann-Cartan background. The physical choice of a gravity theory will be An interesting question for future works is the possi- crucial for the operational realization of the coordinates. ble application of (48) to the analysis of motion of (mi- Recall the form of the equations of motion for a very cro)structured test bodies. In particular, it seems worth- large class [33, 34] of gravitational theories, which also while to search for new ways to map the gravitational allow for additional internal degrees of freedom, in par- field with the help of such a deviation equation and work ticular for spin. In this case, the equations of motion out its implications for various applications, aiming for are no longer given by the geodesic equation or, as it novel tests of relativistic gravity theories. is sometimes erroneously postulated in the literature, by the autoparallel equation. In such theories, test bodies exhibit an additional spin-curvature coupling, which leads to non-geodesic motion, even locally. How does this impact the operational realization of ACKNOWLEDGEMENTS normal coordinates in such theories? Mainly, one just has to be aware of the fact that for the experimental realization of the normal coordinates, one now has to This work was supported by the Deutsche Forschungs- make sure to use the correct equation of motion and, gemeinschaft (DFG) through the grant PU 461/1-1 (D.P.). The work of Y.N.O. was partially supported by PIER (“Partnership for Innovation, Education and Re- search” between DESY and Universität Hamburg) and 3 Here, “locally” refers to the observers laboratory on the reference by the Russian Foundation for Basic Research (Grant world line. No. 16-02-00844-A). 7

Appendix A: Notation & Conventions TABLE I. Directory of symbols. Symbol Explanation The curvature and the torsion are deﬁned w.r.t. the c general connection Γab as follows:

Geometrical quantities d d d d n d n Rabc := ∂aΓbc − ∂bΓac +Γan Γbc − Γbn Γac , g Metric ab (A1) δa b Kronecker symbol c c c xa,ya Coordinates Tab := Γab − Γba . (A2) c i Γab Connection The symmetric Levi-Civita connection Γkj , as well as c Γab Levi-Civita connection all other Riemannian quantities, are denoted by an addi- d Rabc Curvature tional overline. For a general tensor A of rank (n,l) the c Tab Torsion commutator of the covariant derivative thus takes the c Kab Contortion form:

σ World function c1...cn e c1...ck (∇a∇b − ∇b∇a) A d ...d = −Tab ∇eA d ...d y 1 l 1 l η Deviation vector k l y0 g x0 Parallel propagator ci c1...e...ck e c1...ck + Rabe A d1...dl − Rabdj A d1...e...dl . i=1 j=1 X X Misc (A3) ˜ i Y (t), X(t) (Reference) world line In addition to the torsion, we deﬁne the contortion Kkj ua Velocity with the following properties b a Acceleration i i i x x Kkj := Γkj − Γkj , (A4) K y,H y Jacobi propagators 1 Ay1...yn Expansion coeﬃcient K = − (T + T + T ) , (A5) kji 2 kji ikj ijk πy1y2y3y4 Auxiliary quantities i i Tkj = −2K[kj] . (A6) Operators The signature of the spacetime metric is assumed to be ∂i, “ ,” Partial derivative (+1, −1, −1, −1). y y ∇i, “ ;” Covariant derivative As usual, σ x1...y2... := ∇x1 ... ∇y2 ... (σ ) denote the D higher-order covariant derivatives of the world function. dt =“˙” Total cov. derivative d ◦ dt =“ ” Total derivative “[. . . ]” Coincidence limit Appendix B: Intermediate results “ ” Riemannian object Here we give some intermediate results of the derivation of the expansions the world function derivatives around the reference world line Y :

1 1 σ = g + [σ ]σy3 + [σ ] − [σ ] − [σ ][σy5 ] σy3 σy4 + O σ3 , (B1) y1y2 y1y2 y1y2y3 2 y1y2y3y4 y1y2y3 ;y4 2 y1y2y5 y3y4 1 σ = [σ ]+ [σ ] − [σ ] σy4 + [σ ] − [σ ] − 2[σ ] y1y2y3 y1y2y3 y1y2y3y4 y1y2y3 ;y4 2 y1y2y3y4y5 y1y2y3 ;y4y5 y1y2y3y4 ;y5

y0 y0 y4 y5 3 −[σy1y2y3y0 ][σ y4y5 ] + [σy1y2y3 ];y0 [σ y4y5 ] σ σ + O σ , (B2) 1 σ = −g + g y[σ ]σy3 + g y[σ ] − g y2 g [g x ] y1x2 y1x2 x2 y1xy3 2 x2 y1xy3y4 x2 y1y y2 ;y3y4

y y y5 y3 y4 3 −2gx2 [σy1x(y3 ];y4) − gx2 [σy1xy5 ][σ y3y4 ] σ σ + O σ . (B3) 8

[1] T. Levi-Civita. Sur l’écart géodésique. Math. Ann., 97: tions. Class. Quantum Grav., 18:4725, 2001. 291, 1927. [24] S. Manoff. Deviation operator and deviation equations [2] J. L. Synge. The first and second variations of the length over spaces with affine connections and metrics. J. Geom. integral in Riemannian space. Proc. Lond. Math. Soc., Phys., 39:337, 2001. 25:247, 1926. [25] J. W. van Holten. World-line deviations and epicycles. [3] J. L. Synge. On the geometry of dynamics. Phil. Trans. Int. J. Mod. Phys. A, 17:2764, 2002. R. Soc. Lond. A, 226:31, 1927. [26] C. Chicone and B. Mashhoon. The generalized Jacobi [4] F. A. E. Pirani. On the physical significance of the Rie- equation. Class. Quant. Grav., 19:4231, 2002. mann tensor. Acta Phys. Pol., 15:389, 1956. [27] V. Perlick. On the generalized Jacobi equation. Gen. [5] D. Puetzfeld and Y. N. Obukhov. Generalized deviation Rel. Grav, 40:1029, 2008. equation and determination of the curvature in General [28] J. Vines. Geodesic deviation at higher orders via covari- Relativity. Phys. Rev. D, 93:044073, 2016. ant bitensors. Gen. Rel. Grav, 47:59, 2015. [6] J. Plebański. Conformal geodesic deviations. Acta Phys. [29] F. W. Hehl, P. von der Heyde, G. D. Kerlick, and J. M. Pol., 28:141, 1965. Nester. General relativity with spin and torsion: Foun- [7] D. E. Hodgkinson. A modified theory of geodesic devia- dations and prospects. Rev. Mod. Phys., 48:393, 1976. tion. Gen. Rel. Grav., 3:351, 1972. [30] M. Blagojevićand F. W. Hehl. Gauge theories of grav- [8] S. L. Ba˙zański. The relative energy of test particles in itation: A reader with commentaries. Imperial College General Relativity. Nova Acta Leopoldina, 39:215, 1974. Press, London, 2013. [9] S. L. Ba˙zański. A geometric formulation of the Taylor [31] V. N. Ponomarev, A. O. Barvinsky, and Y. N. Obukhov. theorem for curves on affine manifolds. J. Math. Phys. Gauge approach and quantization methods in gravity the- (N.Y.), 17:217, 1976. ory. Nauka, Moscow, 2017. [10] S. L. Ba˙zański. Kinematics of relative motion of test [32] F. W. Hehl, Yu. N. Obukhov, and D. Puetzfeld. On particles in general relativity. Ann. H. Poin. A, 27:115, Poincarégauge theory of gravity, its equations of motion, 1977. and Gravity Probe B. Phys. Lett. A, 377:1775, 2013. [11] S. L. Ba˙zański. Dynamics of relative motion of test parti- [33] D. Puetzfeld and Yu. N. Obukhov. Equations of motion cles in general relativity. Ann. H. Poin. A, 27:145, 1977. in metric-affine gravity: a covariant unified framework. [12] M. Novello, I. Damião Soares, and J. M. Salim. On Jacobi Phys. Rev. D, 90:084034, 2014. fields. Gen. Rel. Grav, 8:95, 1977. [34] Yu. N. Obukhov and D. Puetzfeld. Multipolar test [13] A. N. Aleksandrov and K. A. Piragas. Geodesic struc- body equations of motion in generalized gravity theo- ture: I. Relative dynamics of geodesics. Theoretical and ries. ”Equations of Motion in Relativistic Gravity”, D. Mathematical Physics, 38:48, 1978. Puetzfeld et. al. (eds.), Fundamental theories of Physics, [14] S. Manoff. Lie derivatives and deviation equations in Springer, 179:67, 2015. Riemannian spaces. Gen. Rel. Grav, 11:189, 1979. [35] A. Trautman. Einstein-Cartan theory. Encyclopedia of [15] R. Schattner and M. Trümper. World vectors, Jacobi Mathematical Physics, edited by J.-P. Francoise, G.L. vectors and Jacobi one-forms on a manifold with a linear Naber and S.T. Tsou, Oxford: Elsevier, 2:189, 2006. symmetric connection. J. Phys. A: Math. Gen., 14:2345, [36] Yu. N. Obukhov. Poincarégauge gravity: selected topics. 1981. Int. J. Geom. Meth. Mod. Phys., 03:95, 2006. [16] N. S. Swaminarayan and J. L. Safko. A coordinate-free [37] W. H. Goldthorpe. Spectral geometry and SO(4) gravity derivation of a generalized geodesic deviation equation. in a Riemann-Cartan spacetime. Nucl. Phys. B, 170:307, J. Math. Phys. (N.Y.), 24:883, 1983. 1980. [17] B. Schutz. On generalized equations of geodesic devi- [38] H. T. Nieh and M. L. Yan. Quantized Dirac field in ation. In: Galaxies, Axisymmetric Systems, and Rel- curved Riemann-Cartan background: I. Symmetry prop- ativity, Ed. M.A.H. MacCallum, Cambridge University erties, Green’s function. Annals of Physics (N.Y.), 138: Press, Cambridge, page 237, 1985. 237, 1982. [18] I. Ciufolini. Generalized geodesic deviation equation. [39] N. H. Barth. Heat kernel expansion coefficient: I. An Phys. Rev. D, 34:1014, 1986. extension. J. Phys. A: Math. Gen., 20:857, 1987. [19] S. L. Ba˙zański and N. N. Kostyukovich. Kinematics of [40] S. Yajima. Evaluation of the heat kernel in Riemann- relative motion of charged test particles in general rela- Cartan space. Class. Quant. Grav., 13:2423, 1996. tivity. I. The first electromagnetic deviation. Acta Phys. [41] S. S. Manoff. Auto-parallel equation as Euler-Lagrange’s Pol. B, 18:601, 1987. equation in spaces with affine connections and metrics. [20] S. L. Ba˙zański and N. N. Kostyukovich. Kinematics of Gen. Rel. Grav., 32:1559, 2000. URL gr-qc/0010048. relative motion of charged test particles in general rel- [42] S. S. Manoff. Deviation equations of Synge and Schild ativity. II. The second electromagnetic deviation. Acta over spaces with affine connections and metrics. Int. J. Phys. Pol. B, 18:621, 1987. Mod. Phys., A16:1109, 2001. URL gr-qc/0012015. [21] S. L. Ba˙zański. Hamilton-Jacobi formalism for geodesics [43] B. Z. Iliev. Deviation equations in spaces with a transport and geodesic deviations. J. Math. Phys. (N.Y.), 30:1018, along paths. JINR communication, E2-94-40, Dubna, 1989. 1994, 2003. URL math-ph/0303002. [22] L. Vanzo. A generalization of the equation of geodesic [44] R. J. van den Hoogen. Towards a covariant smoothing deviation. Nuovo Cim. B, 107:771, 1992. procedure for gravitational theories. J. Math. Phys., 58: [23] R. Kerner, J. W. van Holten, and R. Colistete Jr. Rel- 122501, 2017. ativistic epicycles: another approach to geodesic devia- [45] J. L. Synge. Relativity: The general theory. North- 9

Holland, Amsterdam, 1960. spaces. Cambridge: University Press, 1934. [46] E. Poisson, A. Pound, and I. Vega. The motion of point [52] J. A. Schouten. Ricci-Calculus. An introduction to tensor particles in curved spacetime. Living Reviews in Relativ- analysis and its geometric applications. Springer, Berlin, ity, 14(7), 2011. 2nd edition, 1954. [47] P. Szekeres. The gravitational compass. J. Math. Phys., [53] I. G. Avramidi. A covariant technique for the calculation 6:1387, 1965. of the one-loop effective action. Nucl. Phys. B, 355:712, [48] I. Ciufolini and M. Demianski. How to measure the cur- 1991. vature of space-time. Phys. Rev. D, 34:1018, 1986. [54] I. G. Avramidi. Covariant methods for the calculation of [49] O. Veblen. Normal coordinates for the geometry of paths. the effective action in quantum field theory and investi- Proc. Nat. Acad. Sci. (USA), 8:192, 1922. gation of higher-derivative quantum gravity. 1995. URL [50] O. Veblen and T. Y. Thomas. The geometry of paths. hep-th/9510140. Trans. Amer. Math. Soc., 25:551, 1923. [55] A. Z. Petrov. Einstein spaces. Pergamon Press: Oxford, [51] T. Y. Thomas. The differential invariants of generalized New York, 1969.