We model analytically a relativistic problem consisting of two quasi-particles each with mass m in close orbit around a static Schwarzschild black hole with mass M = 1 situated at the center of mass of the system. The angular momentum l of the system is taken to be 2. We model the mass density of the orbiting particles as a δ-function and we assume that there are no deformations. To model the system, we apply the second-order differential equation obtained elsewhere for a dynamic thin matter shell on a Schwarzschild background. As it is the case in this paper, the framework on which the equation was obtained is Bodi-Sachs. The only change in the equation is that now the quasi-normal mode parameter represents the particle’s orbital frequency from which we are able to analytically compute the gravitational radiation emitted by the system at null infinity. We note that in a real astrophysical scenario the dynamics of the particles paths will be very dynamic and complicated and that the analytical methods used here will have to be developed further to accommodate that.
Until recently, all sorts of relativistic binary systems have been studied only theoretically and on the 14 September 2015 a team of LIGO and Virgo collaborators announced their first detection of a gravitational wave signal from a binary black hole system of about 36 and 29 solar masses. This announcement reaffirmed the predictions of the existance of gravitational waves as predicted by GR and most importantly the affirmation that indeed binary relativistic systems do exist in nature. LIGO, Virgo and all other promising gravitational wave detectors will thus provide with the means to be able to detect all kinds of relativistic binary systems with all sorts of physical properties [
More so, the works in this direction have been in the realm of numerical relativity with a special focus of relativistic two body problems i.e. Black hole-black hole binary, black hole-neutron star binary, and neutron star-neutron star binary. The case of a three-body problem as it is the case here has never been studied before either numerically or analytically even though there is a much compelling possibility that in very dense cluster of galaxies these kinds of systems could in fact be found in the near future. As it was the case at the beginning of the research work in relativistic binary systems in the past decades, it is also very likely that there will be arising technical and computational issues for the gravitating three body problem in full numerical relativity. The analytic method used in this paper has been used by the author in [
In the setting of this paper, we use the analytical results by Bishop [
The Bondi-Sachs formalism uses coordinates x i = ( u , r , x A ) based upon a family of outgoing null hypersurfaces. We label these hypersurfaces by u = const . , null rays by x A ( A = 2 , 3 , x 2 = θ , x 3 = ϕ ) , and the surface area coordinate by r. In this coordinates system the Bondi-Sachs metric [
d s 2 = − [ e 2 β ( 1 + W r ) − r 2 h A B U A U B ] d u 2 − 2 e 2 β d u d r − 2 r 2 h A B U B d u d x A + r 2 h A B d x A d x B , (1)
where h A B h B C = δ B A and d e t ( h A B ) = d e t ( q A B ) , with q A B being a unit sphere metric, U is the spin-weighted field given by U = U A q A . For a Schwarzschild space-time, W = − 2 M . We define the complex quantity J by
J = q A q B h A B / 2 . (2)
For the Schwarzschild space-time, we have J and U being zero and thus they can be regarded as a measure of the deviation from spherical symmetry, and in addition, they contain all the dynamic content of the gravitational field in the linearized regime [
For spherical harmonics we use Z s l m rather than s Y l m as basis functions as follows [
Z s l m = 1 2 [ Y s l m + ( − 1 ) m Y s l − m ] for m > 0 Z s l m = i 2 [ ( − 1 ) m Y s l m − Y s l − m ] for m > 0 Z s l 0 = Y s l 0 , (3)
The s = 0 will be omitted in the case s = 0 , i.e. Z l m = Z 0 l m . The s Z l m are orthonormal and real. We assume the following ansatz
J = R e ( J 0 ( r ) e i σ u ) ð 2 Z l m , U = R e ( U 0 ( r ) e i σ u ) ð Z l m , β = R e ( β 0 ( r ) e i σ u ) Z l m , ω = R e ( ω 0 ( r ) e i σ u ) Z l m , (4)
where r0 is the position of the matter shell, and σ the complex frequency mode which is physical damped and which further means that I m ( σ ) > 0 . In the Bondi frame, the field equations splits into;
• the hypersurface equations and the evolution equations given by
R r r : 4 r β , r = 8 π T r r (5)
q A R r A : 1 2 r ( 4 ð β − 2 r ð β , r + r ð ¯ J , r + r 3 U , r r + 4 r 2 U , r ) = 8 π q A T r A (6)
h A B R A B : ( 4 − 2 ð ð ¯ ) β + 1 2 ( ð ¯ 2 J + ð 2 J ¯ ) + 1 2 r 2 ( r 4 ð U ¯ + r 4 ð ¯ U ) , r − 2 ω , r = 8 π ( h A B T A B − r 2 T ) (7)
q A q B R A B : − 2 ð 2 β + ( r 2 ð U ) , r − 2 ( r − M ) J , r − ( 1 − 2 M r ) r 2 J , r r + 2 r ( r J ) , u r = 8 π q A q B T A B , (8)
• and the constraint equations for off the matter shell in the case of vacuum given by
R u u : 1 2 r 3 ( r ( r − 2 M ) ω , r r + ð ð ¯ ω + 2 ( r − 2 M ) ð ð ¯ β − M r ( ð U ¯ + ð ¯ U ) − r 3 ( ð U ¯ + ð ¯ U ) , u + 2 r ω , u ) = 0, (9)
R u r : 1 4 r 2 ( 2 r ω , r r + 4 ð ð ¯ β − ( r 2 ð U ¯ + r 2 ð ¯ U ) , r ) = 0, (10)
q A R u A : 1 4 r 2 ( 2 r ð ω , r − 2 ð ω + 2 r 2 ( r − 2 M ) ( 4 U , r + r U , r r ) + 4 r 2 U + r 2 ( ð ð ¯ U − ð 2 U ¯ ) + 2 r 2 ð ¯ J , u − 2 r 4 U , u r − 4 r 2 ð β , u ) = 0, (11)
Ref. [
x 3 ( 1 − 2 x M ) d 2 J 2 d x 2 + 2 d J 2 d x ( 2 x 2 + i σ x − 7 x 3 M ) − 2 ( x ( l 2 + l − 2 ) / 2 + 8 M x 2 + i σ ) J 2 = 0 (12)
where J 2 ( x ) ≡ d 2 J 0 + / d x 2 and x = 1 / r , x is the compactification factor in this language. Bishop et al. [
We shall use the following algorithm to calculate the gravitational radiation from the system.
• First we use Equation (12) and the constraints Equations (9)-(11) to get the junction conditions for the Bondi-Sachs matric variables U, ω and J at the boundary i.e. shell,
• Second we test if J , J , r , U , U , r , and ω are smooth across the boundary and if this is true, we then
• Calculate the News function at I + .
We consider a system consisting of two point-particles with equal mass m in quasi-orbit around a stationary Schwarzschild black hole with mass M situated at the center of mass r of the particles when l is 2. We take the orbital radius to be at r0 which means that the distance between the particles is 2r0. We take the initial position of particle 1 to be at r0 with θ and ϕ given by π/2 and νu respectively, ν is the orbital frequency and u the orbital period of the particles. We also take the initial position of a particle 2 to be at r0 with θ and ϕ given by π/2 and νu + π respectively. This imply that the rotation in the following figure is in the yz plane. The initial positions of the objects on the figure should not be confused with the actual initial positions just outline which in actual sense should be along the y axis with the particle 1 on the right and the particle 2 on the left.
The dynamics of this problem is governed by Equation (12) and for our numerical calculation purposes we shall use its Ricatti form [
d v d x = 1 + 2 v x 2 ( 1 − 2 x ) ( ( x − v ) ( 2 + i ν x ) − x ( 7 x + 8 v ) ) (13)
where v is the orbital period of the system.
We model the problem as follows, we start by applying Equation (5) with T r r given by
ρ ( 1 − 2 M r ) − 1 , (14)
where the matter density ρ in the background space-time is given by
ρ = ( 1 − 2 M r ) − 1 r 0 2 δ ( r − r 0 ) ( θ − π 2 ) [ δ ( ϕ − ν u ) + δ ( ϕ − ν u − π ) ] . (15)
Inside the particles orbital radius r < r 0 we set
β = 0 , (16)
and outside the particles orbital radius r > r 0 we set
β = Σ l m β l m Z l m . (17)
Now integrating with respect to r we get
Σ l m β l m Z l m = 2π r 0 ( 1 − 2 M r ) − 2 δ ( θ − π 2 ) [ δ ( ϕ − ν u ) + δ ( ϕ − ν u − π ) ] (18)
By multiplying Equation (18) with Z l ′ m ′ we get
Z l ′ m ′ Σ l m β l m Z l m = 2π r 0 ( 1 − 2 M r ) − 2 δ ( θ − π 2 ) [ Z l ′ m ′ δ ( ϕ − ν u ) + Z l ′ m ′ δ ( ϕ − ν u − π ) ] (19)
and integrating over the sphere it simplifies to
β l ′ m ′ = 2 π r 0 ( 1 − 2 M r ) − 2 [ Z l ′ m ′ ( π 2 , ν u ) + Z l ′ m ′ ( π 2 , ν u + π ) ] . (20)
From Equation (20), for m ′ ≠ 0 we the gravitational radiation otherwise we don’t, and that β l ′ m ′ are generally non-zero for even l and m ′ . We now consider the case l ′ = 2 and we note that
β 21 = 0 , (21)
β 2 , − 1 = 0 , (22)
and that
β 20 ≠ 0. (23)
We note that β 20 mode does not vary in time and hence it does not contain the emitted gravitational radiation. Thus we are only interested in β 22 and β 2, − 2 modes. We use the following normalized spherical harmonics
Y 22 = 1 4 15 2π sin 2 θ e 2 i ϕ , (24)
Y 2 , − 2 = 1 4 15 2π sin 2 θ e − 2 i ϕ , (25)
and the fact that
Z 22 = 1 2 ( Y 22 + Y 2 , − 2 ) , (26)
Z 2 , − 2 = i 2 ( Y 2 , − 2 − Y 22 ) , (27)
to get
Z 22 = 2 4 15 2π sin 2 θ cos 2 ϕ (28)
and
Z 2 , − 2 = 2 4 15 2π sin 2 θ sin 2 ϕ . (29)
Thus from Equation (20)
β 22 = 2π r 0 ( 1 − 2 M r 0 ) − 2 2 4 15 2π cos ( 2 ν u ) 2 = 15π r 0 ( 1 − 2 M r 0 ) − 2 cos ( 2 ν u ) , (30)
and similarly
β 2 , − 2 = 15π r 0 ( 1 − 2 M r 0 ) − 2 sin ( 2 ν u ) (31)
and then finally we write
β = 15π r 0 ( 1 − 2 M r 0 ) − 2 ( R e { e 2 i ν u } Z 22 + R e { − i e 2 i ν u } Z 2, − 2 ) (32)
Now taking M = 1 , Equation (32) then becomes
β = 15π r 0 ( 1 − 2 r 0 ) − 2 ( R e { e 2 i ν u } Z 22 + R e { − i e 2 i ν u } Z 2, − 2 ) (33)
We assume that the orbit is at the innermost stable circular orbit (ISCO), so that r = r 0 = 6 . We then found the change in the Schwarzschild coordinate time T for one complete revolution of 92.3436 from which we found the orbital frequency ν of 0.0680.
To now find the numerical solutions to continue Equation (13) we make the spatial coordinate transformation of x = 1 / r which then imply that the ISCO is now at x m n = 1 / 6 . The numerical computations are done in the domains
D + = { 0 < x < x m n } and D − = { x m n < x < 0.5 } , (34)
with numerical solutions v + ( x ) and v − ( x ) respectively. We start the calculation with the transformed Equation (12) given by
U ( x ) = 2 β 0 x − 1 2 x 4 ( 1 − 2 x M ) d 3 d x 3 J ( x ) − x 3 ( x − 2 x 2 M + i ν ) d 2 d x 2 J ( x ) + x ( 2 x + 2 x 2 M + i ν ) d d x J ( x ) − i ν J ( x ) (35)
where U + ( x ) , U − ( x ) are the Bondi metric functions, and β 0 + , β 0 − are the values of the expansion of the light rays β given by Equation (32) in the exterior and interior domains respectively. Bishop [
We define the general solutions for J 2 ( x ) at x m n outside and inside the orbital radius respectively as
J + ( x ) = c 4 + c 1 x + c 2 J 0 + ( x ) , (36)
J − ( x ) = c 9 + c 6 x + c 7 J 0 − ( x ) , (37)
where c4, c1, c2, c9, c6 and c7 are constants to be determined numerically. The functions J 0 + ( x ) and J 0 − ( x ) are analytic near x m n and therefore can be Taylor expand as
J + ( x ) = J 0 + ( x m n ) + ( x − x m n ) d d x J 0 + ( x ) + ( x − x m n ) 2 2 d 2 d x J 0 + ( x ) + ( x − x m n ) 3 6 d 3 d x J 0 + ( x ) , (38)
J − ( x ) = J 0 − ( x m n ) + ( x − x m n ) d d x J 0 − ( x ) + ( x − x m n ) 2 2 d 2 d x J 0 − ( x ) + ( x − x m n ) 3 6 d 3 d x J 0 − ( x ) , (39)
which then results in Equations (36) and (39) being analytic near x m n . We used Matlab ode45 solver to find numerical solutions of the above derivatives in Equations (38) and (39). We used stringent numerical conditions to get the results to about seven significant figures with RelTol of 10−12, AbsTol of 10−12, and the MaxStep of 0.2 × 10 − 5 and the results we found to be
d d x J 0 + ( x ) = 29144 − 2.280672 × 10 5 i , (40)
d 2 d x J 0 + ( x ) = 2.865551 × 10 6 − 1.52335130 × 10 7 i , (41)
d 3 d x J 0 + ( x ) = 4.8870 × 10 7 − 1.8591431 × 10 9 i , (42)
and
d d x J 0 − ( x ) = 13.04337 − 1.31529 i , (43)
d 2 d x J 0 − ( x ) = 1.54689 × 10 2 − 3.19980 × 10 1 i , (44)
d 3 d x J 0 − ( x ) = − 1.12428 × 10 3 − 1.25311 × 10 3 i . (45)
We have tested for the consistency of the above results by using other Matlab solvers; ode23 and ode15s (which uses the Gears method i.e. backward differentiation formulas) and also observed the accuracy of about 15 significant figures. We went further with the test using ode23t which uses the trapezoidal rule, ode23s which is a modified Rosenbrock formula of order 2, and ode23tb which is an implicit Runge Kutta as opposed to ode45 and ode23 and found the consistency of about 8 significant figures and as opposed to 15 significant figures which is also accurate enough. This illustrate how accurate and valid the results are. These results are very crucial in obtaining the emitted gravitational radiation and hence determining the extent of their convergence is of most paramount importance.
From the hypersurface equation Equation (7) rewritten as
− 2 x 2 ω , x = 2 ( 2 − L 2 ) β 0 + L 2 ( L 2 + 2 ) J − x 4 ( x − 4 L 2 U ) , x (46)
we are able to the Bondi metric function ω + ( r ) and ω − ( r ) . But to find the solution the integration should be done analytically where possible. We only need a solution which is valid in a neighborhood of x = x 0 . Henceforth, it is convenient to make the coordinate transformation x → r = 1 / x . Equation (46) can further be rewritten as
2 ( 2 − L 2 ) β 0 + L 2 ( L 2 + 2 ) J + 1 r 2 ( r 4 L 2 U ) , r = 2 ω , r , (47)
where for l = 2 we have L 2 = − 6 . The constraints equations Equations (9), (10), and (11) now simplifies to
R u u : 1 2 r 3 ( ( r 2 − 2 M r ) ω , r r − 6 ω − 12 ( r − 2 M ) β 0 + 12 M r U − 4 r ( r − 2 M ) i ν β 0 + 12 r 3 i ν U + 2 r i ν ω ) = 0, (48)
q A R u A : 1 2 r 2 ( r ω , r − ω + 4 r 3 U , r + r 4 U , r r + 2 r 2 U − 2 M r 3 U , r r − 8 M r 2 U , r − r 2 i ν J − r 4 i ν U , r − 2 r 2 i ν β 0 ) = 0. (49)
which we then apply in the domains D + and D − . Since these constraints are not completely analytic, this means that we should only evaluate them at the ISCO. We use them among others to eliminate the constants c1, c2, c6, and c7. We now assume that we end up with the solutions
ω + ( x ) = c 5 + ω 0 + ( x ) , ω − ( x ) = c 10 + ω 0 − ( x ) , with ω 0 + ( x 0 ) = ω 0 − ( x 0 ) = 0. (50)
Thus, from the constraints R u u − ( r 0 ) , R u u + ( r 0 ) , R u r − ( r 0 ) , R u r + ( r 0 ) , q A R u A − ( r 0 ) , q A R u A + ( r 0 ) , and the hypersurface Equation (47), we found the metric variables U + ( r 0 ) , U − ( r 0 ) , ω + ( r 0 ) , and ω − ( r 0 ) . From which the expressions of the constants c9, c7, c5, and c10, were found.
We now impose the Bondi gauge conditions:
β 0 + = 0 , c 4 = 0 , (51)
which means that for large r, β 0 + = 0 at I + imply that the coordinate time is the same as proper time and that the regularity at I + require c 4 = 0 . We also impose the following junction conditions at r 0 :
J + ( r 0 ) = J − ( r 0 ) , (52)
2 U + ( r 0 ) = U − ( r 0 ) , (53)
β 0 − = − 2 π r 0 ρ ( 1 − 2 M r 0 ) − 1 (54)
ω + ( r 0 ) − ω − ( r 0 ) = − 4 π r 2 ρ . (55)
From the junction conditions, we were able to find the exact numerical values of the constants c1, c2, and c6 at r 0 = 6 . The exact numerical values of the constants c9, c7, c5, and c10 were then found by substituting the values of c1, c2, and c6 back into their expressions. From here we were then able to plot the graphs of the Bondi metric functions J − ( r 0 ) , J + ( r 0 ) , U − ( r 0 ) , U + ( r 0 ) , ω − ( r 0 ) , and ω + ( r 0 ) as observed in the following graphs.
Physically the metric functions J and U have the smooth asymptotic expansion characteristic through out the entire computational domain and this property is confirmed in
implemented correctly and that our numerical methods and the analytical algorithms we implemented to calculating the gravitational radiation worked properly as intended.
Then finally, since we are in the Bondi gauge, we found the gravitational news to be
N + = 1 2 R e ( c 1 i ν e x p ( i ν u ) ) ( − ( l − 1 ) L 2 ( l + 2 ) ) 2 Z l m , (56)
which then further simplify to
N + = R e ( − 0.3778509291 m − 0.5950899448 i m ) , (57)
with the Bondi mass loss of 0.0114 m2. The author in [
I would like to thank Professor Nigel Bishop for suggestions and validation of the numerical methods and results and on improving the manuscript. I would also like to thank the National Research Foundation of South Africa under GUN 2053724 for financial support.
Kubeka, A.S. (2018) The Gravitational Radiation Emitted by Two Quasi-Particles around a Schwarzschild Black Hole. Journal of Modern Physics, 9, 1195-1214. https://doi.org/10.4236/jmp.2018.96072
1) The constraints computed at r0
R u u − ( r 0 ) = 1.000000000 × 10 − 9 ( − 200 r 0 9 c 7 + 6.005358575 × 10 19 i c 7 r 0 3 + 6.8000000 × 10 7 i r 0 7 c 10 + 1.292409115 × 10 16 i ln ( r 0 ) c 7 r 0 4 − 3.000000000 × 10 9 c 10 r 0 6 − 6.692070654 × 10 17 c 7 r 0 5 − 1.440000000 × 10 11 r 0 5 c 9 − 1000 i c 7 r 0 8 − 1.224000000 × 10 9 i c 9 r 0 7 + 2.72000000 × 10 8 i r 0 5 c 10
+ 4.896000000 × 10 9 i r 0 4 c 6 − 3.227234492 × 10 19 i c 7 r 0 2 + 2.084590843 × 10 16 c 7 r 0 7 − 4.896000000 × 10 9 i c 9 r 0 6 + 7.215316909 × 10 17 ln ( r 0 ) c 7 r 0 5 + 1.044710385 × 10 19 i r 0 4 c 7 − 1000 r 0 8 c 7 + 1.884955592 × 10 11 r 0 8 ρ − 3.600000000 × 10 10 r 0 7 c 9 + 8.16000000 × 10 8 i r 0 7 c 6
− 3.418052808 × 10 9 i r 0 7 ρ − 6.836105613 × 10 9 i r 0 8 ρ + 4.488000000 × 10 9 i r 0 6 c 6 + 3.015928948 × 10 11 r 0 8 ρ + 7.179615196 × 10 17 c 7 r 0 4 + 1.098242223 × 10 19 c 7 r 0 3 + 20000 i r 0 9 c 7 + 2.72000000 × 10 8 i r 0 6 c 10 + 1.806026323 × 10 17 l n ( r 0 ) c 7 r 0 6 − 2.163898247 × 10 17 i c 7 r 0 6
+ 8.160000000 × 10 9 i r 0 5 c 6 + 7.212387448 × 10 17 ln ( r 0 ) c 7 r 0 4 + 8.346481884 × 10 16 c 7 r 0 6 − 1.507964474 × 10 11 r 0 6 ρ + 2.111184000 × 10 19 c 7 − 1.200000000 × 10 10 c 10 r 0 5 + 1.382918067 × 10 18 i c 7 r 0 5 − 1.311705542 × 10 16 i r 0 6 ln ( r 0 ) c 7 − 4.896000000 × 10 9 i c 9 r 0 5 − 4.087019553 × 10 15 i r 0 7 ln ( r 0 ) c 7
+ 1.741939934 × 10 19 c 7 r 0 + 7.323651653 × 10 13 r 0 7 l n ( r 0 ) c 7 − 6.491754249 × 10 20 i c 7 r 0 − 1.440000000 × 10 11 r 0 6 c 9 − 4.113828618 × 10 16 i c 7 r 0 7 − 2.563539606 × 10 9 i r 0 9 p + 1.650271349 × 10 19 c 7 r 0 2 − 3.42398706 × 10 15 i r 0 5 l n ( r 0 ) c 7 − 8.031498194 × 10 20 i c 7 − 1.200000000 × 10 10 c 10 r 0 4 ) / ( r 0 7 ( r 0 + 2 ) 2 ) , (58)
R u u + ( r 0 ) = 1 / 2 ( ( r 0 2 − 2 r 0 ) ( − 915.9586340 i c 2 / r 0 3 + 0.8160000000 i c 4 − 0.6052309472 c 2 + 6745.674492 c 2 / r 0 7 + 80948.09432 c 2 / r 0 5 + 7518.667272 i c 2 / r 0 4 + 90224.00728 i c 2 / r 0 5 + 1.503112547 i c 2 − 69.78185262 i c 2 / r 0 2 − 2733.489212 c 2 / r 0 3 + 12 c 1 / r 0 3 + 36.07819002 c 2 / r 0 2 )
− 7518.667272 i c 2 / r 0 2 − 53.40273063 i c 2 r 0 + 0.1360 i r 0 ( 1253.111212 i c 2 / r 0 2 − 457.9793170 i c 2 / r 0 + 7518.667272 i c 2 / r 0 3 − 0.3026154736 c 2 r 0 2 + 22.10459629 c 2 r 0 + 1124.279082 c 2 / r 0 2 + 6745.674528 c 2 / r 0 3 + 12 r 0 c 4 + 8.900455105 i c 2 r 0
+ 69.78185262 i ln ( r 0 ) c 2 + 0.7515562734 i c 2 r 0 2 + 0.4080000000 i r 0 2 c 4 − 1366.744606 c 2 / r 0 + 6 c 1 / r 0 − 36.07819002 ln ( r 0 ) c 2 + c 5 ) + 1.815692842 c 2 r 0 2 − 132.6275777 c 2 r 0 − 6745.674492 c 2 / r 0 2 − 40474.04717 c 2 / r 0 3 − 72 r 0 c 4 + 2747.875902 i c 2 / r 0
− 45112.00363 i c 2 / r 0 3 − 418.6911157 i ln ( r 0 ) c 2 − 4.509337640 i c 2 r 0 2 + 8200.467636 c 2 / r 0 − 36 c 1 / r 0 + 216.4691401 ln ( r 0 ) c 2 − 6 c 5 + 12 r 0 ( 0.5043591226 × 10 − 1 c 2 − 2.568292000 × 10 − 11 c 2 / r 0 − 9.800840000 × 10 − 11 i c 2 / r 0 − 50.69289058 c 2 / r 0 2
− 0.680 × 10 − 1 i c 4 + 201.0479721 i c 2 / r 0 3 − 0.1252593789 i c 2 − 6745.674528 c 2 / r 0 5 + 228.5542745 c 2 / r 0 3 + 1687.678839 i c 2 / r 0 7 − 38.40352765 i c 2 / r 0 2 − 7518.667272 i c 2 / r 0 5 + 2 c 1 / r 0 2 + 2 c 1 / r 0 3 + 2614.546974 c 2 / r 0 4 )
+ 0.8160 i r 0 3 ( 0.5043591226 × 10 − 1 c 2 − 2.568292000 × 10 − 11 c 2 / r 0 − 9.800840000 × 10 − 11 i c 2 / r 0 − 50.69289058 c 2 / r 0 2 − 0.680 × 10 − 1 i c 4 + 201.0479721 i c 2 / r 0 3 − 0.1252593789 i c 2 − 6745.674528 c 2 / r 0 5 + 228.5542745 c 2 / r 0 3 + 1687.678839 i c 2 / r 0 4
− 38.40352765 i c 2 / r 0 2 − 7518.667272 i c 2 / r 0 5 + 2 c 1 / r 0 2 + 2 c 1 / r 0 3 + 2614.546974 c 2 / r 0 4 − 2.448000000 i r 0 2 c 4 ) / r 0 3 , (59)
R u r − ( r 0 ) = 1.111111111 × 10 − 9 ( − 1.272774256 × 10 19 i c 7 r 0 2 + 1.314989772 × 10 20 i c 7 r 0 − 1.589998513 × 10 18 i c 7 r 0 3 + 76 i c 7 r 0 6 + 3.011811822 × 10 20 i c 7 + 4.846534181 × 10 14 i c 7 r 0 4 + 40952 i r 0 5 c 7 − 7.916940003 × 10 18 c 7 + 2.280721975 × 10 17 c 7 r 0 2
− 3.616361705 × 10 18 c 7 r 0 + 8.260193057 × 10 16 c 7 r 0 3 − 3.02160 × 10 5 c 7 r 0 5 − 1.51080 × 10 5 c 7 r 0 6 + 3.392920066 × 10 10 r 0 6 ρ + 2.704645291 × 10 16 c 7 r 0 4 ) / ( r 0 6 ( r 0 + 2 ) ) , (60)
R u r + ( r 0 ) = 1 / 2 ( r 0 ( − 915.9586340 i c 2 / r 0 3 + 0.8160000000 i c 4 − 0.6052309472 c 2 + 6745.674492 c 2 / r 0 4 + 80948.09432 c 2 / r 0 5 + 7518.667272 i c 2 / r 0 4 + 90224.00728 i c 2 / r 0 5 + 1.503112547 i c 2 − 69.78185262 i c 2 / r 0 2 − 2733.489212 c 2 / r 0 3 + 12 c 1 / r 0 3
+ 36.07819002 c 2 / r 0 2 ) + 6 r 0 2 ( 2.568292000 × 10 − 11 c 2 / r 0 2 + 9.800840000 × 10 − 11 i c 2 / r 0 2 + 101.3857812 c 2 / r 0 3 − 603.1439163 i c 2 / r 0 4 + 33728.37264 c 2 / r 0 6 − 685.6628235 c 2 / r 0 4 − 6750.715356 i c 2 / r 0 5 + 76.80705530 i c 2 / r 0 3 + 37593.33636 i c 2 / r 0 6
− 4 c 1 / r 0 3 − 6 c 1 / r 0 4 − 10458.18790 c 2 / r 0 5 ) + 12 r 0 ( 0.5043591226 × 10 − 1 c 2 − 2.568292000 × 10 − 11 c 2 / r 0 − 9.800840000 × 10 − 11 i c 2 / r 0 − 50.69289058 c 2 / r 0 2 − 0.680 × 10 − 1 i c 4 + 201.0479721 i c 2 / r 0 3 − 0.1252593789 i c 2 − 6745.674528 c 2 / r 0 5
+ 228.5542745 c 2 / r 0 3 + 1687.678839 i c 2 / r 0 4 − 38.40352765 i c 2 / r 0 2 − 7518.667272 i c 2 / r 0 5 + 2 c 1 / r 0 2 + 2 c 1 / r 0 3 + 2614.546974 c 2 / r 0 4 ) ) / r 0 2 , (61)
q A R u A − ( r 0 ) = 1.000000000 × 10 − 10 ( 16000 i c 7 r 0 8 − 1.000000 × 10 6 r 0 8 c 7 + 6.95123928 × 10 18 c 7 r 0 4 + 4.60818516 × 10 20 i c 7 r 0 + 6.811699261 × 10 15 i c 7 r 0 7 + 1.202064574 × 10 18 ln ( r 0 ) c 7 r 0 4 + 1.140199529 × 10 19 i c 7 r 0 5 − 7.391764651 × 10 20 i c 7 r 0 3 − 1.220608500 × 10 14 c 7 r 0 7 + 4.461943440 × 10 21 i c 7
+ 300 r 0 7 p + 5.385037980 × 10 15 i ln ( r 0 ) c 7 r 0 6 − 2.000000000 × 10 10 c 10 r 0 4 + 2.154015192 × 10 16 i ln ( r 0 ) c 7 r 0 5 − 2.671516126 × 10 21 i c 7 r 0 2 + 2.040000000 × 10 9 i r 0 6 c 6 + 3.005161436 × 10 17 ln ( r 0 ) c 7 r 0 6 + 8.160000000 × 10 9 i r 0 4 c 6 − 2.000000000 × 10 10 c 10 r 0 5 − 5.000000000 × 10 9 c 10 r 0 6
+ 1.202064574 × 10 18 l n ( r 0 ) c 7 r 0 5 − 3.817318694 × 10 18 c 7 r 0 5 − 6.543021412 × 10 16 i r 0 6 c 7 + 8.809405089 × 10 19 c 7 r 0 3 + 1.455381063 × 10 20 c 7 r 0 2 + 2.154015192 × 10 16 i l n ( r 0 ) c 7 r 0 4 − 1.44797289 × 10 19 c 7 r 0 + 8.160000000 × 10 9 i r 0 5 c 6 + 4.272566009 × 10 9 i r 0 9 ρ − 3.509551026 × 10 17 c 7 r 0 6
+ 1.001418528 × 10 20 i r 0 4 c 7 − 1.172880000 × 10 20 c 7 + 8.545132018 × 10 9 i r 0 8 ρ ) / ( r 0 6 ( r 0 + 2 ) 2 ) , (62)
q A R u A + ( r 0 ) = 1 / 2 ( − 1253.111212 i c 2 / r 0 2 − c 5 − 1124.279082 c 2 / r 0 2 + 0.3026154736 c 2 r 0 2 − 22.10459629 c 2 r 0 − 12 r 0 c 4 + 36.07819002 l n ( r 0 ) c 2 − 8 r 0 2 ( 2.568292000 × 10 − 11 c 2 / r 0 2 + 9.800840000 × 10 − 11 i c 2 / r 0 4 + 101.3857812 c 2 / r 0 3 − 603.1439163 i c 2 / r 0 4 + 33728.37264 c 2 / r 0 6
− 685.6628235 c 2 / r 0 4 − 6750.715356 i c 2 / r 0 5 + 76.80705530 i c 2 / r 0 3 + 37593.33636 i c 2 / r 0 6 − 4 c 1 / r 0 3 − 6 c 1 / r 0 4 − 10458.18790 c 2 / r 0 5 ) + r 0 ( 457.9793170 i c 2 / r 0 2 + 0.8160000000 i r 0 c 4 + 8.900455105 i c 2 − 0.6052309472 c 2 r 0 + 22.10459629 c 2
− 2248.558164 c 2 / r 0 3 − 20237.02358 c 2 / r 0 4 + 12 c 4 − 2506.222424 i c 2 / r 0 3 − 22556.00182 i c 2 / r 0 4 + 1.503112547 i c 2 r 0 + 69.78185262 i c 2 / r 0 + 1366.744606 c 2 / r 0 2 − 6 c 1 / r 0 2 − 36.07819002 c 2 / r 0 ) + 4 r 0 3 ( 2.568292000 × 10 − 11 c 2 / r 0 2
+ 9.800840000 × 10 − 11 i c 2 / r 0 2 + 101.3857812 c 2 / r 0 3 − 603.1439163 i c 2 / r 0 4 + 33728.37264 c 2 / r 0 6 − 685.6628235 c 2 / r 0 4 − 6750.715356 i c 2 / r 0 5 + 76.80705530 i c 2 / r 0 3 + 37593.33636 i c 2 / r 0 6 − 4 c 1 / r 0 3 − 6 c 1 / r 0 4 − 10458.18790 c 2 / r 0 5 )
+ r 0 4 ( − 5.136584000 × 10 − 11 c 2 / r 0 3 − 1.960168000 × 10 − 10 i c 2 / r 0 3 − 304.1573436 c 2 / r 0 4 + 2412.575665 i c 2 / r 0 5 − 2.023702358 × 10 5 c 2 / r 0 7 + 2742.651294 c 2 / r 0 5 + 33753.57678 i c 2 / r 0 6 − 230.4211659 i c 2 / r 0 4 − 2.255600182 × 10 5 i c 2 / r 0 7
+ 12 c 1 / r 0 7 + 24 c 1 / r 0 5 + 52290.93950 c 2 / r 0 6 ) + 2 r 2 2 ( 0.5043591226 × 10 − 1 c 2 − 2.568292000 × 10 − 11 c 2 / r 0 − 9.800840000 × 10 − 11 i c 2 / r 0 − 50.69289058 c 2 / r 0 2 − 0.680 × 10 − 1 i c 4 + 201.0479721 i c 2 / r 0 3 − 0.1252593789 i c 2 − 6745.674528 c 2 / r 0 5
+ 228.5542745 c 2 / r 0 3 + 1687.678839 i c 2 / r 0 4 − 38.40352765 i c 2 / r 0 2 − 7518.667272 i c 2 / r 0 5 + 2 c 1 / r 0 2 + 2 c 1 / r 0 3 + 2614.546974 c 2 / r 0 4 ) − 8.900455105 i c 2 r 0 − 0.2720 i r 0 2 ( c 4 + c 1 / r 0 + c 2 ( 1 + ( 13.04336905144130 − 1.31528646137769 i ) ( 1 / r 0 − 1 / 6 )
+ ( 77.34402850 − 15.99899824 i ) ( 1 / r 0 − 1 / 6 ) 2 + ( − 187.3798480 − 208.8518687 i ) ( 1 / r 0 − 1 / 6 ) 3 ) ) + 457.9793170 i c 2 / r 0 − 0.680 × 10 − 1 i r 0 4 ( 2.568292000 × 10 − 11 c 2 / r 0 2 + 9.800840000 × 10 − 11 i c 2 / r 0 2 + 101.3857812 c 2 / r 0 2
− 603.1439163 i c 2 / r 0 4 + 33728.37264 c 2 / r 0 6 − 685.6628235 c 2 / r 0 4 − 6750.715356 i c 2 / r 0 5 + 76.80705530 i c 2 / r 0 3 + 37593.33636 i c 2 / r 0 6 − 4 c 1 / r 0 3 − 6 c 1 / r 0 4 − 10458.18790 c 2 / r 0 5 ) − 7518.667272 i c 2 / r 0 3 − 69.78185262 i ln ( r 0 ) c 2
− 0.7515562734 i c 2 r 0 2 − 2 r 0 3 ( − 5.136584000 × 10 − 11 c 2 / r 0 3 − 1.960168000 × 10 − 10 i c 2 / r 0 3 − 304.1573436 c 2 / r 0 4 + 2412.575665 i c 2 / r 0 5 − 2.023702358 × 10 5 c 2 / r 0 7 + 2742.651294 c 2 / r 0 5 + 33753.57678 i c 2 / r 0 6 − 230.4211659 i c 2 / r 0 4 − 2.255600182 × 10 5 i c 2 / r 0 7
+ 12 c 1 / r 0 4 + 24 c 1 / r 0 5 + 52290.93950 c 2 / r 0 6 ) − 0.4080000000 i r 0 2 c 4 − 6745.674528 c 2 / r 0 3 + 1366.744606 c 2 / r 0 − 6 c 1 / r 0 ) / r 0 2 , (63)
2) The Bondi metric variables computed at r0
U + ( r 0 ) = 0.5043591226 × 10 − 1 c 2 − 2.568292000 × 10 − 11 c 2 / r 0 − 9.800840000 × 10 − 11 i c 2 / r 0 − 50.69289058 c 2 / r 0 2 − 0.680 × 10 − 1 i c 4 + 201.0479721 i c 2 / r 0 3 − 0.1252593789 i c 2 − 6745.674528 c 2 / r 0 5
+ 228.5542745 c 2 / r 0 3 + 1687.678839 i c 2 / r 0 4 − 38.40352765 i c 2 / r 0 2 − 7518.667272 i c 2 / r 0 5 + 2 c 1 / r 0 2 + 2 c 1 / r 0 3 + 2614.546974 c 2 / r 0 4 , (64)
U − ( r 0 ) = ( 4.887000002 × 10 7 − 1.859143100 × 10 9 i ) c 7 / r 0 5 + 4 c 7 ( 2.865551000 × 10 6 − 1.523351300 × 10 7 i + ( 4.887000002 × 10 7 − 1.859143100 × 10 9 i ) ( 1 / r 0 − 1 / 6 ) / r 0 4 + ( − 2.443500001 × 10 7 + 9.295715500 × 10 8 i ) c 7 / r 0 4 + ( 2 ( c 6 + c 7 ( 29144 − 2.280672000 × 10 5 i
+ ( 2.865551000 × 10 6 − 1.523351300 × 10 7 i ) ( 1 / r 0 − 1 / 6 ) + ( 2.443500001 × 10 7 − 9.295715502 × 10 8 i ) ( 1 / r 0 − 1 / 6 ) 2 ) ) ) / r 0 3 − c 7 ( 2.865551000 × 10 6 − 1.523351300 × 10 7 i + ( 4.887000002 × 10 7 − 1.859143100 × 10 9 i ) ( 1 / r 0 − 1 / 6 ) ) / r 0 3 − ( 0.680 × 10 − 1 i ) c 7 ( 2.865551000 × 10 6 − 1.523351300 × 10 7 i
+ ( 4.887000002 × 10 7 − 1.859143100 × 10 9 i ) ( 1 / r 0 − 1 / 6 ) / r 0 2 + ( 2 ( c 6 + c 7 ( 29144 − 2.280672000 × 10 5 i
+ ( 2.865551000 × 10 6 − 1.523351300 × 10 7 i ) ( 1 / r 0 − 1 / 6 ) + ( 2.443500001 × 10 7 − 9.295715502 × 10 8 i ) ( 1 / r 0 − 1 / 6 ) 2 ) ) ) / r 0 2 + ( 0.680 × 10 − 1 i ) ( c 6 + c 7 ( 29144 − 2.280672000 × 10 5 i
+ ( 2.865551000 × 10 6 − 1.523351300 × 10 7 i ) ( 1 / r 0 − 1 / 6 ) + ( 2.443500001 × 10 7 − 9.295715502 × 10 8 i ) ( 1 / r 0 − 1 / 6 ) 2 ) ) / r 0 − ( 0.680 × 10 − 1 i ) ( c 9 + c 6 / r 0 + c 7 ( − 2518.141205 + 699.71985 i
+ ( 29144 − 2.280672000 × 10 5 i ) ( 1 / r 0 − 1 / 6 ) + ( 1.432775500 × 10 6 − 7.616756500 × 10 6 i ) ( 1 / r 0 − 1 / 6 ) 2 + ( 8.145000002 × 10 6 − 3.098571834 × 10 8 i ) ( 1 / r 0 − 1 / 6 ) 3 ) ) , (65)
ω + ( r 0 ) = ( 1253.111212 i ) c 2 / r 0 2 − ( 457.9793170 i ) c 2 / r 0 + ( 7518.667272 i ) c 2 / r 0 3 − 0.3026154736 c 2 r 0 2 + 22.10459629 c 2 r 0 + 1124.279082 c 2 / r 0 2 + 6745.674528 c 2 / r 0 3 + 12 r 0 c 4 + ( 8.900455105 i ) c 2 r 0
+ ( 69.78185262 i ) ) ln ( r 0 ) c 2 + ( 0.7515562734 i ) c 2 r 0 2 + ( 0.4080000000 i ) r 0 2 c 4 − 1366.744606 c 2 / r 0 + 6 c 1 / r 0 − 36.07819002 ln ( r 0 ) c 2 + c 5 , (66)
and
ω − ( r 0 ) = − ( 2156.071282 i ) c 7 r 0 2 − 50.26548247 r 0 2 ρ / ( 1 + 2 / r 0 ) + 6 c 6 / r 0 − 5.147566058 × 10 5 c 7 r 0 2 − 63413.86124 c 7 r 0 − 2.932200001 × 10 8 c 7 / r 0 3 − 4.886999995 × 10 7 c 7 / r 0 2 − 2.356233015 × 10 8 c 7 / r 0 − 6.010322872 × 10 7 ln ( r 0 ) c 7 + 12 r 0 c 9 + ( 1.513990017 × 10 7 i ) c 7 r 0
− ( 1.031580813 × 10 9 i ) c 7 / r 0 + ( 1.859143099 × 10 9 i ) c 7 / r 0 2 + ( 1.115485860 × 10 10 i ) c 7 / r 0 3 + ( 0.4080000000 i ) r 0 2 c 9 − ( 1.077007596 × 10 6 i ) ln ( r 0 ) c 7 + c 10 , (67)
The computed constants
c 5 = − 5.925925926 × 10 − 21 ( 2.257876807 × 10 25 i c 2 r 0 + 1.177568763 × 10 22 i ln ( r 0 ) c 2 r 0 4 − 2.537550206 × 10 25 i c 2 − 1.179321011 × 10 22 i c 2 r 0 4 − 6.885000000 × 10 19 i r 0 4 c 1 − 2.276665152 × 10 25 c 2 − 1.379990759 × 10 20 i c 2 r 0 5 + 8.836025110 × 10 9 c 2 r 0 6 − 6.088194566 × 10 21 ln ( r 0 ) c 2 r 0 4
− 3.251077700 × 10 23 i c 2 r 0 3 − 6.873541974 × 10 23 i c 2 r 0 2 − 2.669155863 × 10 20 c 2 r 0 5 + 8.950474631 × 10 21 c 2 r 0 4 − 4.171449115 × 10 23 c 2 r 0 3 + 3.018766507 × 10 9 i c 2 r 0 6 − 2.509446061 × 10 24 c 2 r 0 2 + 2.174295963 × 10 25 c 2 r 0 ) / r 0 4 , (68)
c 7 = ( − 1.683000000 × 10 14 + 1.416666667 × 10 25 i ) c 6 / ( − 3.518394608 × 10 32 + 9.145010221 × 10 31 i − 1.463828647 × 10 23 ln ( r 0 ) + ( 2.060201741 × 10 22 i ) ln ( r 0 ) ) , (69)
c 9 = − 1.111111111 × 10 − 11 ( 1.157988867 × 10 17 r 0 8 c 7 − 3.769911183 × 10 14 r 0 7 ρ + 5.22000 × 10 5 r 0 9 c 7 − 8.250022346 × 10 21 c 7 r 0 7 + 1.884955592 × 10 14 r 0 6 ρ − 2.356194490 × 10 14 r 0 8 ρ − 5.423600000 × 10 11 ln ( r 0 ) c 7 r 0 5 + 9.635910263 × 10 22 c 7 r 0 2 − 1.737134306 × 10 20 c 7 r 0 6
+ 5.0000000 × 10 7 i r 0 7 ln ( r 0 ) c 7 + 6.926736238 × 10 21 c 7 r 0 3 − 8.000000 × 10 6 i r 0 9 c 7 − 3.240686912 × 10 20 c 7 r 0 5 − 1.963050508 × 10 24 i c 7 r 0 2 − 2.726124482 × 10 19 c 7 r 0 7 − 2.500000 × 10 6 r 0 7 ln ( r 0 ) c 7 + 1.159077064 × 10 24 i c 7 r 0 − 1.143558000 × 10 23 c 7 + 2.863090381 × 10 20 i r 0 6 c 7
+ 2.075034463 × 10 15 i c 7 r 0 8 + 1.452672443 × 10 11 r 0 9 ρ + 2.000000000 × 10 9 i ln ( r 0 ) c 7 r 0 6 − 6.319234788 × 10 23 i c 7 r 0 3 + 1.495398102 × 10 13 i r 0 8 ρ + 6.249786891 × 10 19 i c 7 r 0 7 + 4.350394854 × 10 24 i c 7 − 5.000000000 × 10 11 ln ( r 0 ) c 7 r 0 4 + 4.321899263 × 10 22 c 7 r 0 + 6.408849014 × 10 12 i r 0 9 ρ
− 9.236000000 × 10 10 ln ( r 0 ) c 7 r 0 6 + 6.704677816 × 10 21 i c 7 r 0 5 − 4.080000000 × 10 12 i r 0 6 c 6 − 1.020000000 × 10 12 i r 0 7 c 6 − 4.080000000 × 10 12 i r 0 5 c 6 + 6.054991092 × 10 22 i r 0 4 c 7 + 1.387200000 × 10 11 r 0 5 c 6 + 3.468000000 × 10 10 r 0 7 c 6 + 7.000000000 × 10 9 i ln ( r 0 ) c 7 r 0 5 + 1.387200000 × 10 11 r 0 6 c 6
+ 2.500000000 × 10 9 i ln ( r 0 ) c 7 r 0 4 + 4.272566010 × 10 12 i r 0 7 ρ + 7.263362215 × 10 10 r 0 10 ρ ) / ( ( 17 i r 0 2 + 500 r 0 2 + 68 i r 0 + 2000 r 0 + 2000 + 68 i ) r 0 5 ) , (70)
and
c 10 = 2.000000000 × 10 − 10 ( − 3.509551026 × 10 17 c 7 r 0 6 + 1.202064574 × 10 18 ln ( r 0 ) c 7 r 0 5 + 8.809405089 × 10 19 c 7 r 0 3 − 1.447972890 × 10 19 c 7 r 0 − 3.817318694 × 10 18 c 7 r 0 5 + 1.455381063 × 10 20 c 7 r 0 2 − 1.220608500 × 10 14 c 7 r 0 7 + 1.140199529 × 10 19 i c 7 r 0 5 − 1.172880000 × 10 20 c 7
+ 2.040000000 × 10 9 i r 0 6 c 6 + 6.951239280 × 10 18 c 7 r 0 4 + 3.005161436 × 10 17 ln ( r 0 ) c 7 r 0 6 + 4.608185160 × 10 20 i c 7 r 0 + 4.461943440 × 10 21 i c 7 + 2.154015192 × 10 16 i ln ( r 0 ) c 7 r 0 5 + 8.545132018 × 10 9 i r 0 8 ρ + 1.202064574 × 10 18 ln ( r 0 ) c 7 r 0 4 + 1.001418528 × 10 20 i r 0 4 c 7 + 16000 i c 7 r 0 8
+ 5.385037980 × 10 15 i ln ( r 0 ) c 7 r 0 6 + 8.160000000 × 10 9 i r 0 4 c 6 + 6.811699261 × 10 15 i c 7 r 0 7 − 7.391764651 × 10 20 i c 7 r 0 3 − 2.671516126 × 10 21 i c 7 r 0 2 + 2.154015192 × 10 16 i ln ( r 0 ) c 7 r 0 4
− 1.000000 × 10 6 r 0 8 c 7 + 4.272566009 × 10 9 i r 0 9 ρ − 6.543021412 × 10 16 i r 0 6 c 7 + 8.160000000 × 10 9 i r 0 5 c 6 + 300 r 0 7 ρ ) / ( r 0 4 ( 4 + 4 r 0 + r 0 2 ) ) , (71)
c 1 = ( − 3.303262047 − 3.089304784 i ) m , (72)
c 2 = ( − 393.9477195 − 46.65966681 i ) m , (73)
c 6 = ( − 394.7245520 − 8.295859050 i ) m , (74)