Spin contribution to the perihelion precession in binary systems like OJ287: higher order corrections. Carlos Mar´ın∗ Department of Physics, Universidad San Francisco de Quito, Diego de Robles y V´ıa Interoce´anica Jorge Poveda† Department of Physics, Universidad San Francisco de Quito, Diego de Robles y V´ıa Interoce´anica (Dated: December 15, 2020) Abstract Higher order corrections are obtained for the perihelion precession in binary systems like OJ287, Sagittarius A*-S2 and H1821+643 using both the Schwarzschild metric and the Kerr metric to take into account the spin effect. The corrections are performed considering the third root of the motion ′ equation and developing the expansion in terms of ǫ r / a(1 e2) and ǫ∗ 1 2αE ǫ.The ≡ s − ≡ − cJ results are compared with those obtained in a previous paper . arXiv:2012.06665v1 [gr-qc] 11 Dec 2020 1 I. INTRODUCTION In a previous paper1, higher order corrections (up to n-th order) were obtained for the perihelion precession (see figure (1)) in binary systems like OJ287, Sagittarius A*-S2 and H1821+643 using the Schwarzschild metric and complex integration. The corrections were performed considering the third root of the motion equation and developing the expansion in terms of ǫ r / (a(1 e2)), where r is the Schwarzschild radius.The results were compared ≡ s − s with other expansions that appear in the literature giving corrections to second and third order2–8. In this paper we will consider both Schwarzschild and Kerr metrics to take into 2 account the spin effect and developing the expansion in terms of both ǫ rs/ (a(1 e )) ′ ≡ − and ǫ∗ 1 2αE ǫ, where E′ and J are the energy and angular momentum per unit mass ≡ − cJ and α is a factor proportional to the black hole spin. Kerr’s black holes are very interesting because most of the black holes in the universe probably have a rotational movement (spin) and one of the most important consequences of that spin is that spacetime is dragged around the rotating black hole leading to an effect that is known as ”frame dragging”. In the case of a Schwarzschild black hole, a precession occurs which is the rotation of the elliptical orbit in the fixed plane of that orbit. For Kerr’s solution, the drag of the reference frame introduces an additional precession of the plane of the orbit around the axis of rotation of the black hole in the same direction of rotation of said black hole. To facilitate the calculations we will assume that the axis of rotation of the Kerr black hole is perpendicular to the plane of the orbit. II. SPIN CONTRIBUTION TO THE PERIHELION PRECESSION The Schwarzschild metric describes a body with spherical symmetry, but without electric charge and rotational movement. If we take into account that many black holes that have been found have rotation on their own axis, as for example the binary system OJ287, one might think that this rotation should influence the advance of the perihelion of the elliptical orbits. The rotation around its own axis is given by the angular momentum of spin Sz of the massive body M. To include the spin in the calculation, you can use the Kerr metric. The Kerr metric is a solution to the field equations in vacuum for a body of mass M that rotates on its own axis 2 FIG. 1. Perihelion precession. ω is the initial inclination of the orbit and δω is the angle of precession. 9–15 with an angular momentum Sz. The Kerr metric is : r2 α2 r α2 2r α (ds)2 = c2(dτ)2 = γc2(dt)2 (dr)2 r2 1+ + s (dφ)2 + s cdtdφ, (1) − ∆ − r2 r3 r where α = Sz ,∆= r2 r r + α2, with coordinates x0 = ct, x1 = r, x2 = θ and x3 = φ. Mc − s 2GM rs = c2 is the Schwarzschild radius. We have taken θ = π/2 (equatorial plane). In this paper we only intend to observe how the spin could contribute to the perihelion precession α and if it is relevant to the calculation, for which only first order terms in r will be taken into account. By doing this, the metric is reduced to: 1 2r α (ds)2 = c2(dτ)2 = γc2(dt)2 (dr)2 r2(dφ)2 + s cdtdφ. (2) − γ − r 2 µ ν The arc length ds satisfies the relation ds = gµνdx dx , then, the covariant metric tensor is: 3 rsα γ 0 0 r 0 γ−1 0 0 gµν = − (3) 0 000 rsα 0 0 r2. r − The geodesic equation can be written in an alternative form using the Lagrangian dxµ dxα dxβ L xµ, = g (xµ) (4) dσ − αβ dσ dσ where σ is a parameter of the trajectory of the particle, which is usually taken to be the proper time, τ for a massive particle. Using the Euler-Lagrange equations: ∂L d ∂L µ dxµ =0 (5) ∂x − dσ ∂ dσ ! we get the geodesic equation for the particle: du 1 µ = (∂ g ) uαuβ (6) dσ 2 µ αβ dxµ where uµ = dσ . With σ = τ, for the coordinates ct (µ = 0) and φ (µ = 3), the geodesic equation (6) give us, respectively: d dt r αc dφ γc2 + s =0, (7) dτ dτ r dτ d dφ r αc dt r2 s =0 (8) dτ dτ − r dτ Both of these equations define the following constants along the trajectory of the particle around the massive object: 4 dt r αc dφ γc2 + s = E′, (9) dτ r dτ dφ r αc dt r2 s = J (10) dτ − r dτ where E′ has units of energy per unit mass and J of angular momentum per unit mass. From these equations we get: dφ 1 r α = s E′ + γcJ (11) dτ γcr2 r dt 1 r2 r α = E′ s J (12) dτ γcr2 c − r Equation (2) can be written as: dt 2 1 dr 2 dφ 2 2r α dt dφ c2 = γc2 r2 + s c , (13) dτ − γ dτ − dτ r dτ dτ dt dφ and replacing the values of dτ and dτ : r2 r α 2 dr 2 r α 2 γ2r4c4 = γc2 E′ s J γc2r4 r2 s E′ + γcJ c − r − dτ − r 2r α r α r α + s r2E′ s cJ s E′ + γcJ . r − r r α Simplifying and taking only first order terms in r : r4 2rr α dr 2 r α 2r α γr4c4 = c2 E′2 s JE′ c2r4 r2 2c s JE′ + γc2J 2 + s r2cJE′ . c2 − c − dτ − r r Introducing the value of γ in the left side of the last equation we have: r c2 1 2r α dr 2 J 2 c2 s = E′2 s JE′ γ − r c2 − cr3 − dτ − r2 Finally we obtain an energy conservation equation, similar to the one obtained in the case of the Schwarzschild metric, but with an additional crossed term proportional to JE′ : 5 E′2 dr 2 J 2 r c2 2r α c2 = + γ s + s JE′ (14) c2 − dτ r2 − r cr3 This alllows us to define an effective potential: J 2 r c2 2r α V = γ s + s JE′ (15) r2 − r cr3 From equation (14), we can gete the expression for the radial kinetic energy per unit mass: dr 2 r c2 J 2 J 2r 2r α = A + s + s s JE′, (16) dτ r − r2 r3 − cr3 ′2 E 2 where A = 2 c . c − Recall now that since the orbit is an ellipse, there are two points in which the temporal derivative becomes zero, and they are aphelion and perihelion. For these two points we can write: 2 2 2 rsc J J rs 2rsα ′ A + 2 + 3 3 JE =0 (17) Ra − Ra Ra − cRa 2 2 2 rsc J J rs 2rsα ′ A + 2 + 3 3 JE =0, (18) Rp − Rp Rp − cRp where R = a (1 + e) and R = a (1 e). a p − For an ellipse, the equation of motion have three real and positive roots. Two of the ′ ′ roots are Ra and Rp and the other we will call Ro. To obtain this root (Ro) we use equation (16) : dr 2 r3 A r c2 2r αE′ = r3 + s r2 r + r s =0, dτ J 2 J 2 J 2 − s − cJ and write: A r c2 2r αE′ A r3 + s r2 r + r s = (r R )(r R )(r R′ ). J 2 J 2 − s − cJ J 2 − a − p − o 6 The last equation can be written as: r c2 2r α A s r2 r + r s E′ = | | (R′ + R + R ) r2 J 2 − s − cJ J 2 o a p A A | | (R R′ + R R + R R′ ) r + | |R R R′ . − J 2 p o a p a o J 2 a p o It is important to recall that A is negative for elliptic orbits, so it can be written as A = A . −| | Comparing the coefficients of r0, r and r2 of both sides of the last equation, and replacing the values of Ra and Rp we obtain the following relations: r c2 A s = | | (R′ +2a) , (19) J 2 J 2 o A 1= | | 2aR′ + a2 1 e2 (20) J 2 o − J 2r 2α R′ = s 1 E′ (21) o A (1 e2) a2 − cJ | | − ′ 1 The value of R0 given in (21) is the same result that we obtained in a previous paper , but with an extra factor 1 2α E′. Using equations (20) and (21) we can finally write: − cJ ′ 2 2αE ′ a (1 e ) rs 1 cJ R = − − ′ , (22) 0 a (1 e2) 2r 1 2αE − − s − cJ and 3 J a 2 (1 e2) 1 = 1 (23) − ′ A 2 a (1 e2) 2r 1 2αE 2 | | − − s − cJ dr dr dφ dφ Going back to equation (16), and because dτ = dφ dτ , replacing the expression of dτ given by (11), we get: 7 dr 2 1 r α 2 r c2 J 2 J 2r 2r α s E′ + γcJ = A + s + s s JE′ (24) dφ γ2c2r4 r r − r2 r3 − cr3 α Leaving only first order terms in r we arrive to the following equation: dr 2 2r αE′ A r c2 2r α 1+ s = r4 + s r3 r2 + r s E′ r (25) dφ γrcJ J 2 J 2 − s − cJ that also can be written as: dr 2 2r αE′ A 1+ s = | | (R r)(r R )(r R′ ) r.