Journal of Modern Physics
Vol.09 No.06(2018), Article ID:84499,20 pages
10.4236/jmp.2018.96072

The Gravitational Radiation Emitted by Two Quasi-Particles around a Schwarzschild Black Hole

Amos S. Kubeka

Department of Mathematical Sciences, University of South Africa, Pretoria, South Africa

Copyright © 2018 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: July 19, 2017; Accepted: May 12, 2018; Published: May 15, 2018

ABSTRACT

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.

Keywords:

Gravitational Radiation, Schwarzschild Black Hole, Quasi-Particles, Quasi-Normal Mode Null Formalism

1. Introduction

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 [1] - [23] .

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 [24] to validate other analytical results obtained by [25] [26] for a binary system.

In the setting of this paper, we use the analytical results by Bishop [27] to study analytically the relativistic triple system consisting of two point-particles in quasi-orbit around a static Schwarzschild black hole. In practice the particles could be either two black holes, two neutron stars or in another setting a combination of both. Our objective thus is to determine the amount of the emitted gravitational radiation by the system at I + in Bondi-Sachs formalism. This paper is structured as follows: Section 2 gives the background material. Section 3 defines the physical problem to be studied. Section 4 calculates the emitted gravitational radiation at I + .

2. Background

2.1. Formalism

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 [28] [29] takes the form

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 [30] . Usually we can describe this space-time by β = 0 and W = 2 M , or by β = β c ( constant ) and W = ( e 2 β c 1 ) r 2 M .

For spherical harmonics we use Z s l m rather than s Y l m as basis functions as follows [27]

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. [27] got the following second order differential equation when solving the above systems of ordinary differential equations for the Schwarzschild background;

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. [31] solved Equation (12) numerically and obtained interesting quasi-normal modes results of a Schwarzschild white hole. However in this paper, we are going to solved it for a different problem since we can apply the same physical settings in the Bondi-frame to model our problem with σ having a different physical meaning as we shall see later.

2.2. An analytic Algorithm for Calculating the Gravitational News

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 + .

3. The Problem

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 [31]

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.

4. The Emitted Gravitational Radiation

4.1. The Linear Expansion of the Light Rays From the System to I +

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 = 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 = 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 sin 2 θ e 2 i ϕ , (24)

Y 2 , 2 = 1 4 15 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 sin 2 θ cos 2 ϕ (28)

and

Z 2 , 2 = 2 4 15 sin 2 θ sin 2 ϕ . (29)

Thus from Equation (20)

β 22 = r 0 ( 1 2 M r 0 ) 2 2 4 15 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)

4.2. The Gravitational Radiation

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 [27] has indicated that the derivatives of J should not be worked out numerically, but should be worked out analytically in terms of J 1 , J 2 and v from Equation (13) with ν = 0.0680 .

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 1012, 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 Figure 1 and Figure 2. The metric function ω do not have this physical property as can be confirmed in Figure 3 but this function is crucial in the calculation procedure of the gravitation radiation in the entire domain. Physically the function J in the only one that have the time derivative and thus carries the gravitational radiation information to calculated at I + and that all the other Bondi metric functions are intergrated radially from Γ to I + . The above results indicate that the junction conditions at r 0 = 6 where

Figure 1. The graph of I m ( J 2 + ( r 0 ) ) for I m ( J ) in the entire domain, and R e ( J 2 + ) for R e ( J ) also in the entire for the Schwarzschild space-time. ν = 0.07 and l = 2 .

Figure 2. The graph of R e ( U ( r 0 ) ) , I m ( U ( r 0 ) ) and R e ( U + ( r 0 ) ) , I m ( U + ( r 0 ) ) for the Schwarzschild space-time. ν = 0.07 and l = 2 .

Figure 3. The graph of R e ( ω ( r 0 ) ) , I m ( ω ( r 0 ) ) and R e ( ω + ( r 0 ) ) , I m ( ω + ( r 0 ) ) for the Schwarzschild space-time. ν = 0.07 and l = 2 .

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 [24] has done a similar work for a single point particle in close orbit around a Schwarzschild black hole in the Bondi-frame and obtained the Bondi mass loss of 0.00089897 m2. He succeeded in validating the results by comparing it with that of the 5.5 PN formalism by Poisson [25] and Sasaki et al. [26] for the same problem. Thus the methods used in this article open up the possibilities in numerical relativity to be able to study analytically the gravitational radiation emitted by a systems consisting of one black hole and two equal orbiting black holes/neutron stars. With further improvement, the method can be develop to look at two unequal orbiting black holes or neutron stars or a combination of both with efficiency and accuracy as demonstrated in [24] for single orbiting black hole/neutron star.

Acknowledgements

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.

Cite this paper

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

References

  1. 1. Baumgarte, T.W., Skoge, M.L. and Shopiro, S.L. (2004) Physical Review D, 70, Article ID: 064040. https://doi.org/10.1103/PhysRevD.70.064040

  2. 2. Carter, B. and Luminet, J.P. (1983) Astronomy & Astrophysics, 121, 97-113.

  3. 3. Carter, B. and Luminet, J.P. (1985) Astronomy & Astrophysics, 121, 23.

  4. 4. Chandrasekhar, S. (1969) Ellipsoidal Figures of Equilibrium. Yale University Press, New Heaven.

  5. 5. Faber, J.A., Baumgarte, T.W., Shapiro, S.L. and Taniguchi, K. (2006) The Astrophysical Journal, 641, L93. https://doi.org/10.1086/504111

  6. 6. Faber, J.A., Baumgarte, T.W., Shapiro, S.L., Taniguchi, K. and Rasio, F.A. (2006) Physical Review D, 73, Article ID: 024012. https://doi.org/10.1103/PhysRevD.73.024012

  7. 7. Fishbone, L.G. (1973) The Astrophysical Journal, 185, 43-67. https://doi.org/10.1086/152395

  8. 8. Grandclément, P. (2006) Physical Review D, 74, Article ID: 124002. https://doi.org/10.1103/PhysRevD.74.124002

  9. 9. Grandclément, P. (2007) Physical Review D, 75, Article ID: 129903(E).

  10. 10. Ishii, M., Shibata, M. and Mino, Y. (2005) Physical Review D, 71, Article ID: 044017. https://doi.org/10.1103/PhysRevD.71.044017

  11. 11. Kobayashi, S., Laguna, P., Phinney, E.S. and Mészáros, P. (2004) The Astrophysical Journal, 615, 855-865. https://doi.org/10.1086/424684

  12. 12. Lai, D. and Wiseman, A.G. (1996) Physical Review D, 54, 3958-3964. https://doi.org/10.1103/PhysRevD.54.3958

  13. 13. Lee, W.H. (2000) Monthly Notices of the Royal Astronomical Society, 318, 606-624. https://doi.org/10.1046/j.1365-8711.2000.03870.x

  14. 14. Löffler, F., Rezzollas and Ansorg, M. (2006) Physical Review D, 74, Article ID: 104018. https://doi.org/10.1103/PhysRevD.74.104018

  15. 15. Mashhoon, B. (1975) The Astrophysical Journal, 197, 705-716. https://doi.org/10.1086/153560

  16. 16. Miller, M.C. (2005) The Astrophysical Journal, 626, L41-L44. https://doi.org/10.1086/431583

  17. 17. Postnov, K.A. and Yungelson, L.R. (2006) Living Reviews in Relativity, 9, 6. http://www.livingreviews.org/lrr-2006-6

  18. 18. Rosswog, S., Speith, R. and Wynn, G.A. (2004) Monthly Notices of the Royal Astronomical Society, 351, 1121-1133. https://doi.org/10.1111/j.1365-2966.2004.07865.x

  19. 19. Shibata, M. and Uryū, K. (2006) Physical Review D, 74, Article ID: 121503(R).

  20. 20. Shibata, M. and Uryū, K. (2007) Classical and Quantum Gravity, 24, S125-S137. https://doi.org/10.1088/0264-9381/24/12/S09

  21. 21. Sopuerta, C.F., Sperhake, U. and Laguna, P. (2006) Classical and Quantum Gravity, 23, S579-S598. https://doi.org/10.1088/0264-9381/23/16/S15

  22. 22. Taniguchi, K., Baumgarte, T.W., Faber, J.A. and Shapiro, S.L. (2006) Physical Review D, 74, Article ID: 041502(R).

  23. 23. Taniguchi, K., Baumgarte, T.W., Faber, J.A. and Shapiro, S.L. (2007) Physical Review D, 75, Article ID: 084005. https://doi.org/10.1103/PhysRevD.75.084005

  24. 24. Kubeka, A.S. (2013) Applications of Linearized Perturbations of the Schwarzschild Geometry within the Bondi-Sachs Formalism. PHD Thesis.

  25. 25. Poisson, E. (1993) Physical Review D, 47, 1497-1510. https://doi.org/10.1103/PhysRevD.47.1497

  26. 26. Sasaki, M. and Tagoshi, H. (2003) Living Reviews in Relativity, 6, 6. https://doi.org/10.12942/lrr-2003-6

  27. 27. Bishop, N.T. (2005) Classical and Quantum Gravity, 22, 2393-2406.

  28. 28. Bondi, H., van der Burg, M.J.G. and Metzner, A.W.K. (1962) Proceedings of the Royal Society of London A, 269, 21-52. https://doi.org/10.1098/rspa.1962.0161

  29. 29. Sachs, R.K. (1963) Physical Review Letters, 150, 66.

  30. 30. Bishop, N.T., Gómez, R., Lehner, L., Maharaj, M. and Winicour, J. (1997) Physical Review D, 56, 6298-6309. https://doi.org/10.1103/PhysRevD.56.6298

  31. 31. Bishop, N.T. and Kubeka, A.S. (2009) Physical Review D, 80, Article ID: 064011. https://doi.org/10.1103/PhysRevD.80.064011

Appendix

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 i1.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)