1. Introduction
A binary gravitational rotator, also called the two-body problem, is a pair of masses m1, m2 moving around their center-of-mass (com) in their own gravitational field.
The problem can be formulated as a single rotator under certain conditions.
In Newtonian gravitation, the two-body problem can be described by a single reduced mass (single rotator)
orbiting around the total mass
situated in com in the distance
, which is the distance between the two original masses.
In Newtonian case, the com has zero acceleration, i.e. it is a well-defined Lorentzian frame.
In section 2, we discuss the two-body problem in Newtonian, Schwarzschild and Kerr spacetime context, and show that the (single) rotator is well-defined also in Schwarzschild spacetime, and in Kerr spacetime it is well-defined, if the rotational parameters of the two spinning Kerr-masses are equal.
We describe the general GR-ansatz for the two-body problem in section 3, and formulate the corresponding Kerr orbit equations in section 4.
In section 5, we adapt the Kerr rotational parameter to the Newtonian correction of the rotator-potential.
In section 6, we present the exact GR-solution (Manko-Ruiz) of the two-body problem, which is a generalized Kerr spacetime with four parameters, and adapt it to the Newtonian correction.
In section 7, we discuss extensions: the Einstein radiation power formula and complex rotation parameters in the generalized Kerr spacetime.
In subsection 8.1, we carry out numerical orbit calculations for observer orbit in rotator potential for all involved models and compare them.
In subsection 8.2, we derive the formula for the angular momentum of the Manko-Ruiz rotator in its own field, and calculate and discuss the exact rotator orbits in two examples.
The formulas for metrics, orbit equations, and Christoffel symbols are calculated symbolically from fundamental equations in [1] and [2], and are therefore errorless. In the text, they are typed by hand where needed, but are inserted additionally as images in the appendix to serve as protection against typos.
2. The Rotator in Newtonian, Schwarzschild and Kerr Spacetime
2.1. Schwarzschild and Kerr Spacetime
We start with exact solutions of Einstein equations in spherical coordinates for the non-rotating (Schwarzschild) and rotating (Kerr) black-hole.
The Kerr line element reads [3] [4] (with metric signature
)
(1)
where
is the Schwarzschild radius, and
is the angular momentum radius (amr), α has the dimension of a distance:
, and J is the angular momentum.
In the limit α→ 0 the Kerr line element becomes the standard Schwarzschild line element
(2)
The total energy for a mass m in Newtonian gravitation field of a mass M is:
(3)
where Et is the total energy and
the relative total energy. We use in the following the terminology of [4] for the GR energy and radial orbit equation:
, where
is the (dimensionless) relativistic velocity factor, and
the (here negative) relative total energy,
the absolute relative total energy.
Because of conservation of angular momentum L is
, l = reduced angular momentum is a constant. Using this relation, (3) becomes the Newtonian orbit differential equation for the orbit radius r, with the parameters l and
to be determined from the initial condition.
From the first Schwarzschild orbit equation (see below) we get
[4], where F is the above relativistic velocity factor and has the dimensionality of velocity
.
In the general relativistic Schwarzschild case the Newtonian approximation (3) becomes the exact relativistic energy equation (radial orbit equation [5] [6] [7]):
(4)
or with Schwarzschild radius
.
, the relative total energy is negative for a bound state:
, we use in the following also the absolute relative energy
.
We can write the energy equation using the effective potential [5]
(4a)
or still simpler using the general effective potential
(4b)
The second form is convenient for expansion of the energy equation in section 6.
The equation is solvable in integral form
For the Newtonian case, the
term is missing, and the integral can be calculated in closed form (see 5.1).
2.2. Newtonian, Schwarzschild, Kerr Rotator
Newtonian rotator
In Newtonian gravitation, the movement of two masses with location vectors
and
in their own gravitational potential takes a very simple form, when formulated in their center-of-mass (com) reference frame [8] [9].
The location of com is
and we get for its acceleration from the third Newtonian law
, so
, com has constant velocity and is a well-defined reference frame, in which the com relation
is valid.
For the displacement vector
we have then for a central force
then
, or
and
where
is the reduced mass: i.e. the movement of m1 and m2 is described by the movement of mass
with distance
in com reference frame, this orbit is the single rotator.
Schwarzschild rotator
We consider now the case of the Schwarzschild spacetime gravitational rotator, also called the Schwarzschild two-body problem: two point masses m1 and m2 with
rotating around their center-of-mass (com), in the com reference-frame with orbit radii r1 and r2 respectively.
We make the approximation, that the Schwarzschild potentials of the two masses add up to the total spacetime, which is approximate, since the Einstein equations are non-linear. Under this assumption, every mass “sees” the Schwarzschild potential of the other mass, like in Newtonian case.
The original Schwarzschild energy equation reads [5] [6]
with
(5)
and for the two masses
we use the two Schwarzschild orbit invariants l and F and get
Because of the com-condition
, r1 and r2 can be calculated from the distance r0between m1 and m2
and
where
is the reduced mass, and
is the total mass, furthermore we have
for the rotation angles, and the reduced angular momenta
,
, where
, so
is the reduced angular momentum of the rotator, with the relation
.
From the above follows
,
, so
,
and the com acceleration vanishes:
, so the com frame has constant velocity relative to the observer, the Schwarzschild single rotator is well-defined.
We add the two orbit equations
and from
we get the well-known single rotator orbit equation
(6)
The single rotator has the orbit of the reduced mass
in the potential
with energy
and reduced angular momentum
.
In order to calculate the individual relative energies
, we consider
, since
for
. From
follows
,
Kerr rotator
In case of two rotating relativistic masses, the underlying individual spacetime is the Kerr spacetime.
We consider the 1/r-expanded Kerr energy equation to third order and get (dimensionless
,
)
(7)
Like in the Schwarzschild case, we get (here with error
) and setting
,
(8)
,
with rotation parameters
and
.
Now the condition for a well-defined com reference frame is
, i.e.
.
Under this condition, the movement of masses m1 and m2is described by the binary Kerr rotator orbit
with the radial potential
, with reduced angular momentum
, rotation parameter
, and relative energy
.
2.3. The GR Field and Orbit Equations
The Einstein field equations are [5] [6] [8]:
(9)
where
is the Ricci tensor, R0 the Ricci curvature,
,
is the energy-momentum tensor,
is the cosmological constant. In the following
is neglected, i.e. set
, because it is important only on cosmological scale, and here we consider a distance scale of star binaries, i.e. d = 10, …, 300 AU.
Further on, we use the Christoffel symbols (second kind)
(10)
and the Ricci tensor
(11)
The geodesic orbit equations O1…O4 in vacuum (
) are:
,
(12)
where τ = proper time.
The four equations O1, O2, O3, O4 are the relativistic orbit equations for
i.e. time t and the three spatial coordinates x, y, z in dependence on the relativistic proper time τ For τ we get for the line-element ds = dτ and therefore trivially:
(13)
This relation yields for the Kerr- and Schwarzschild-spacetimes the GR energy relation, we choose the denomination E1s for it, the orbit equations for the Kerr spacetime [1] [3] O1s, O2s, O3s, O4s are given in the appendix [9].
Setting [τ] = π/2 solves O3s = 0, and simplifies considerably the other equations.
The simplified expressions E1s1, O1s1, O2s1, O4s1 are given in the appendix [9].
In the following, we use the expression for the τ-derivative with dot or with prime:
.
3. The General Ansatz for the GR Two-Body Problem
The General Formulation of GR Two-Body Problem
The GR GR two-body problem is the GR-gravitationally bound binary system of two masses m1 and m2 moving around their center-of-mass (com) with respective distances (r1,r2)from com [10] [11].
The GR two-body problem has 4 independent parametersthe total mass
(resp. corresponding Schwarzschild radius
), the mass ratio
, and two orbit parameters: the reduced angular momentum
and the relative total energy
(resp. its absolute value
).
Equivalently we can use the two masses m1 and m2,
= mean distance, and l, which defines the eccentricity by the formula
.
The corresponding spacetime has the form
In comparison, Schwarzschild spacetime has only one parameter
.
(see matrix form in the appendix) and Kerr spacetime has 2 parameters:
and
the angular momentum radius
(see matrix form in the appendix)
The general formulation of the GR two-body-problem is as follows [11] [12].
We formulate the problem in coordinates relative to center-of-mass,
and
for the two masses, and
the metric of the GR two-body-system, with the com-equations
,
, and
proper-time of the metric
. We can eliminate
using the com-equations and are left with
and
as the 10 + 4 unknown variables of the problem.
Equivalently, we can express
and
by the rotator 4-vector
in center-of-mass (com) spatial spherical coordinates, with r being the distance between m1 and m2,
and
where
is the reduced mass, and
is the total mass. In Newtonian limit, this is a correct formulation of the two-body problem (see above), so it is also valid in the GR case.
We have then 10 Einstein-equations for the metric
,
where the energy-mass tensor is [4]:
, and where for black-holes
(or dimensionless
), and the 4 geodesic orbit equation for
So in principle, we can solve the 14 coupled partial deq’s for the 14 unknowns imposing appropriate boundary conditions.
We show in section 6.3 that the Manko-Ruiz spacetime
is the exact GR-solution of the relativistic two-body problem. However, the Manko-Ruiz solution is very complicated, since it has a complexity of
terms.
If the direct evaluation of the exact solution is too demanding, one is forced to use some iterative scheme.
The common approach is the so-called post-Newtonian approximation [10].
Here we make a series ansatz for the GR-Hamiltonian
in the reduced location 3-vector
and reduced momentum 3-vector
, where
(from com-equations), unit vector
, reduced mass ratio
.
Note that here the reduced momentum has the dimensionality
of velocity, and the reduced location has the dimensionality
of 1/velocity2, so the dimensionality of the Hamiltonian
is energy/mass.
The Newtonian Hamiltonian is then
and the first post-Newtonian Hamiltonian is
Now we can apply the Hamiltonian equations to get the approximate GR equation-of-movement
.
4. The Kerr Orbit Equations
Of special importance is the solution of O1, which gives the derivative
.
In the Newtonian approximation, is of course
and
.
In the Schwarzschild spacetime, O1 can be solved analytically, and the well-known solution is [5] [6]
, where
the velocity factor, and reduced energy
are constants of the orbit and result from initial conditions.
Furthermore in Schwarzschild spacetime, angular momentum L is conserved:
, where l is the reduced angular momentum, and the azimuthal angle is constant, we can choose
.
In the Kerr spacetime, the solution cannot be given in analytical form, but it can be expressed as a series in r and α, it seems that it is derived here for the first time.
First, we bring O1 into a new form using
,
and
, thus eliminating
and
:
(14)
This has the general form
and after multiplication with a function
it can be made a total differential
with
And with this condition the formal solution can be derived immediately:
and
with an integration constant F.
In the Schwarzschild case with α = 0 and
this results immediately in
and
.
In the Kerr case,
,
and after turning the integral in the numerator of
into a series in α and 1/r,
becomes to first order in α
which for α = 0 results again in
.
So the first-order Kerr-correction to
is of the order
from the F-term (total energy) and of the order
from the l-term (rotational energy).
Now we can eliminate
,
,
, using the above relations in the Kerr energy equation E1d, and we get the radial equation for
in the form
, where
is called the effective potential.
Inserting the Schwarzschild (dimensionless) invariants:
and
yields the simplified Kerr energy equation in α, l, and positive relative total energy
(see also appendix [9]), with
(15)
and for pure Schwarzschild α = 0 and we recover the Schwarzschild energy equation
E1 with Kerr-replacement (full Kerr)
(16)
The Kerr radial equation E1dS0 can be expanded in 1/r powers:
(17)
5. Adaptation of the Kerr-Parameter α to the Newtonian Rotator Correction
5.1. The Newtonian and Schwarzschild Orbit
The Newtonian Kepler orbit is a fundamental formula of classical celestial mechanics. We shall use it in the following sub-chapter and describe it here as a reminder.
The Kepler orbit results from the Newtonian gravitational potential
where
is the Schwarzschild radius of the central mass M.
The Kepler orbit is an ellipse with eccentricity e described by [9].
, where the reduced angular momentum is
, the eccentricity is
, and
is the relative total energy from the radial orbit equation,
for bound states with elliptical orbits,
is the absolute value of the (negative) relative total energy.
The orbit rotation angle
is time-dependent and obeys the relation
.
The Newtonian orbit Equation (4) is solvable in integral form
(18)
The resulting formula is [13] (dimensionless, time in Schwarzschild-units
, radius in Schwarzschild-units
, c = 1).
(19)
The zeros of the root in the integral are the minimal and maximal radius, i.e. the semi-minor and the semi-major axis of the orbit ellipse [13].
,
dimensionless
,
full dimensional, mean distance is
.
For circular orbit,
,
,
.
The minimal and maximal radii
are zeros of the Newtonian effective potential (4b) (dimensionless, i.e.
).
, which yields a quadratic equation for the two radii.
In case of the Schwarzschild metric with the corresponding effective potential
the minimal and maximal radii
are roots of a cubic equation in r.
The first root
is the turning point of a parabolic orbit, the other two are
of a closed Schwarzschild orbit [2], see appendix .
When
become complex, there is no closed Schwarzschild orbit, this happens for small distances from the center. The mean radius
has the form ( [2])
for small energies
, it becomes
.
5.2. The Orbit Equations for the Remote Observer in Kerr Spacetime and in the Newtonian Limit
First, we calculate the two-body correction
to the Newtonian gravitational potential from the Newtonian Kepler orbits of the rotatoras shown in Figure 1.
Second, we take into account the Kerr correction
in the Kerr energy E1dS0 and compare both expressions, from this we get a formula for
.
The z-axis of the observer is the z-axis of the rotator: perpendicular to the rotator plane, the origin is at rotator-com.
We make the assumptions:
- Observer distance from com is large compared to the rotator’s orbit diameter, so
= rotation angle of the rotator, changes much faster than
= rotation angle of the observer orbit.
- The distance vector between observer and center-of-mass:
, i.e. observer orbit movement is much slower than rotator movement
We have the following denominations:
= vector(observer, com rotator)
Cartesian observer coordinates
spherical coordinates relative to com
,
distance between masses
,
,
masses:
,
,
,
, mass ratio
rotator orbit angle
= angle(rotator-axis, observer x-projection) observer mass-distance projections:
,
distances between observer and masses
,
are
We calculate now the correction of the Newtonian potential due to the two-body rotation of the rotator.
The Newtonian potential of the rotator mass
is
,
Figure 1. The Newtonian correction to the observer orbit in rotator potential.
We develop the distances
and
Now, we have to average over a rotation period over the angle
. Theoretically, we have to insert
from the orbit equation
and invert
to get
, and then integrate over a rotation period
.
It is hopeless to get an analytical expression for
, but we can assume that
is approximately linear
, which is exact for a circular orbit:
. Then the term
cancels by averaging over a rotation period, and we get the expression
, where the average
,
is the mean diameter of the rotator orbit ellipse,
,
With this approximation, we get the correction of the Newtonian potential
or dimensionless
(20)
The term
in E1dS0 is with
dimensionless
Now, for circular orbit dimensionless
, so
, and we get for
, and
and
, as it should be, because for
the rotator becomes a single mass, and the cylindrical symmetry becomes spherical, the Kerr-spacetime becomes Schwarzschild-spacetime with
.
6. The Manko-Ruiz Solution as a Generalized Kerr Spacetime
6.1. Tomimatsu-Sato Solutions as Generalizations of Kerr Spacetime
The Kerr solution belongs to the class of Tomimatsu–Sato solutions describing exterior gravitational fields of stationary rotating axisymmetric sources introduced by Tomimatsu and Sato (1972, 1973).
In Weyl’s canonical (cylindrical) coordinates (ρ, z), the axisymmetric Papapetrou space-time metric ( [14] 19.17) is given by
(21)
with functions U, k, A.
We get the Einstein field equations in the form
where
,
For
we obtain in Weyl coordinates the Ernst equation ( [14] 19.41)
(22)
from which we calculate the metric functions
,
, where
The Ernst equation with Γ can be reformulated by introducing the new potential
by
We go over to the prolate spheroidal coordinates
connected with Weyl’s canonical coordinates
by the relations
,
,
In these coordinates, the Ernst equation becomes ( [14] 20.37)
Based on this ansatz, the Tomimatsu-Sato solution with integer δ can be obtained via a limiting process from the non-linear superposition of δ Kerr solutions with common symmetry axis (Tomimatsu and Sato, 1981).
The potential ξ of these solutions is a quotient
, α and β being polynomials in the coordinates x and y, with parameters
where
. Here, the constant σ in the relation between the Weyl coordinates
and the prolate spheroidal coordinates
, the angular momentum J and the quadrupole moment Q are given by
,
,
, where δ is an integer
For δ = 1 we get the Kerr solution in prolate coordinates
,
,
The full Kerr metric can be obtained from ξ as ( [14] 20.36)
,
,
where mq = a, mp = σ
,
,
,
,
For δ = 2 we get the Kinnersley and Chitre (1978) solution
,
,
6.2. Metric for Kerr-Kerr Binary Configuration
We consider the spacetime of two Kerr black-holes with masses m1 and m2, rotation parameters a1 and a2, in fixed distance R, and we are looking for the corresponding vacuum solution of Einstein equations, which is asymptotically flat, i.e. it becomes Minkowski metric at infinity (Manko & Ruiz [15]), it is also called the Kinnersley-Chitre five-parameter vacuum solution, Manko & Ruiz [16].
This spacetime generalizes the Tomimatsu-Sato solution for δ = 2 [16].
Ernst complex potential Γ of the exact solution for two aligned Kerr sources is defined on the symmetry axis (ρ = 0) by the expression (Stephani [14] (34.94), Manko & Ruiz [15])
where m is the total mass of the binary system, a is the rotational parameter, and the real quantities s, μ, τ, δ are related to the physical characteristics:
,
,
,
We have the 5 real physical parameters {m1, m2, a1, a2, R}, with distance R, total mass m and total angular momentum J
,
, rotational parameters
,
, total rotational parameter a depends on J and is obtained from the cubic equation in
and
are the half-horizon radii
,
where
Ernst complex potential Γ in the whole
space is then extended to
(23)
with the radial functions
,
,
which satisfies the Ernst equation (Cabrera [13])
(24)
The line element becomes now [13] in the Weyl-form from 6.1 with new denominations
,
,
(24a)
with the functions f, γ, ω of the Weyl coordinates ρ and z
,
,
, (24b)
where
6.3. Transformation to Spherical Coordinates
It is advisable to transform the metric from the Weyl cylindrical coordinates
to spherical coordinates
.
The transformation reads:
,
, where
is the total mass [14].
The differentials in the line element also have to be transformed correspondingly:
and
. The transformed metric doubles in complexity [1] [2],
,
, but handling the Christoffel symbols and the Einstein equations in familiar spherical coordinates is easier.
The resulting Christoffel symbols have the complexity
terms, and the Einstein equations have the complexity
terms, so symbolic verification of the Einstein equations is out of the question, but numeric verification by stochastic insertion works (although it is time-consuming), and the result is indeed zero within the limit of precision for this complexity (10−4) [2].
In the following and in the program code, we use exclusively this transformed Manko-Ruiz metric.
6.4. Adaptation of the Manko-Ruiz Kerr-Kerrpotential to the Newtonian Correction for Observer Orbit
The Manko-Ruiz Kerr-Kerr-solution
has originally 5 parameters: two rotator masses m1, m2, the distance between the masses R, two orbit rotation parametera1, a2. We interpret R as the mean rotator distance rc, and a1, a2 as the two individual reduced angular momenta of the rotator
masses. In this case, the rotation parameters have the ratio
, where
(see 2.1), so we are left with 4 independent parameters:
, where
,
.
The 4 parameters have to be adapted to the Newtonian rotator correction: in order
:
, and in order
:
, preserving the mass ratio
; in order
and
the Manko-Ruiz metric is Schwarzschild.
The general effective potential of the Manko-Ruiz metric in the energy (or radial) equation is
, it depends on the 4 metric parameters
and the two orbit parameters
of the observer orbit.
The two free parameters
in
are in fact functions of the rotator orbit parameter
, i.e.
,
. The coefficient of
in
is
.
It is a quadratic polynomial in
, as evaluation shows [2]:
, the same is true for
,
.
Adaptation means: in the effective orbit potential
the
independent coefficient satisfies
(25)
and the coefficient
.
This results in two equations for the two parameters
, which are solvable, usually in real numbers, as the numerical example shows.
The adaptation cannot be done symbolically, because the radial equation potential
for Manko-Ruiz has a complexity of
terms, and to calculate the expansion in
-powers takes several weeks on a
desktop. In comparison, for Kerr-spacetime
with the expansion computing time
. The assessment for the computing time is
, so we get
.
A numerical alternative is to fit the expansion in
-powers of the radial equation
with inserted parameters
, for the two expansion coefficients: for
coefficient of
, and for
coefficient of
.
This is done in [2] in section “Interpolation Manko-Ruiz coeff-expansion (1/r2), coeff-expansion (1/r3), in R, a” as a fit in
with fixed
.
We find the solution of the relation (25) to be [2]:
We use the approximate values
for the calculation of Manko-Ruiz observer orbits with
in subsection 8.1.
6.5. Manko-Ruiz Rotator Orbits in Its Own Potential
The Manko-Ruiz Kerr-Kerr-solution
has originally 5 parameters, we replace the two masses
by the total mass
, and the ratio
. The rotational parameters
are proportional to the reduced orbit angular momentum of the two masses:
and
, as in the case of the Kerr metric, where
. We have shown in 2.2 that the angular momentums of the rotator have the ratio
, so we can set
and
, where a is the rotational parameter of the rotator.
Now the Manko-Ruiz metric becomes
as a function in spherical coordinates
with 3 parameters
and the mass m. We make the metric dimensionless by setting, as usual,
and
, which is equivalent to calculation in Schwarzschild units
and
.
In order to calculate Manko-Ruiz rotator orbits in its own potential, we have to derive a formula for the orbit parameters
in dependence of the three Manko-Ruiz parameters
:
and
.
Relative energy
Here
is the mean radius of the orbit, and we know from 5.1 that the Kepler formula for
is
.
As for
, in the case of Manko-Ruiz rotator, we calculate
at fixed
by iteration with the ansatz
: we start with
and
, calculate the orbit from the Manko-Ruiz orbit equations, calculate
, and adapt the factor
in such a way to make
, The iteration normally converges after 2 - 3 steps with accuracy
, as demonstrated in [2].
In principle, we can calculate symbolically the roots
of the Manko-Ruiz effective potential
, calculate
, and solve the equation
for
.
In reality, the equation is of course hopelessly complicated.
However, we can solve the equation numerically for given discrete parameter values
on lattice in a certain parameter region, and then fit the values
in three dimensions. This has been done in [2] in subsection “fit epsc = c1(rc)/(4rc)” and the result is
, where the coefficient
is fitted by a rational quadratic function
(26)
which depends only on R, and the fit-data plot is shown below Figure 2:
Figure 2. The ε-coefficient
as a function of R in Schwarzschild units.
For R > 20,
, and we have the Newtonian relation (19).
Reduced angular momentum
.
In order to derive the formula
, we start with the off-diagonal metric component
of the Kerr metric. In the orbit plane
and we get
, i.e.
is the coefficient of
in the expansion of
in
-powers.
We apply this to the Manko-Ruiz metric. Symbolic calculation is out of the question, but we can fit the expansion of
numerically in the parameters
, as we did in 6.4. in parameters
.
This is done in [1] in subsections “interpolation M14” and “fit M14”.
The result is the formula
(27)
which we use for the calculation of rotator orbits in 8.2.
7. Extended Interpretations of Rotator Energy Equation
7.1. The Gravitational Radiation Power
The Einstein-formula gravitational radiation power of a two-body rotator is as presented in [5]:
(28)
where
is the mean distance of the rotator. or expressed by the Schwarzschild radius
, which shows manifestly dimensionality energy/time, or dimensionless
We compare it to the Kerr energy term
in E1dS0 with
, where we have to set
, since we consider now the rotator orbit in his own metric,
and calculate the corresponding power = energy/period dimensionless setting
and have
, or dimensionless
Now, for circular Keplerian orbits,
(dimensionless), therefore
, so the term is approximately
We see, that none of the term has the correct order in
, and, since the term
in E1dS0 contains no
-term, there is no appropriate term in Kerr potential, the Kerr approximation cannot explain radiation loss.
In fact, there is another important precondition for radiation lass, this energy term is extracted from the gravitational system, therefore it must be imaginary in the potential
. In optics and quantum mechanics, imaginary terms in energy are related to dissipation and wave damping.
What about the Manko-Ruiz metric?
One can see from numeric
-expansion of the Manko-Ruiz potential
, that there is a
-term in the
coefficient, so the Manko-Ruiz potential yields an appropriate term for radiation energy. However, in order to be imaginary, the Manko-Ruiz metric itself has to be made complex and still satisfy the Einstein equations.
In fact, there is a natural way to do it.
In the basic equations of the Manko-Ruiz metric, in (24b) the line element coefficient of
is defined as
In fact, this equation is a little artificial, it is made so, to ensure that the metric is real.
The setting
is natural, and of course it satisfies the Einstein equations as well.
The Manko-Ruiz metric thus becomes manifestly complex, and the orbits are calculated from complex orbit equations. With this ansatz, the solution
becomes complex, and we simply take the absolute of the variables
, as all of them must be positive, as below in subsection 7.2.
So in Manko-Ruiz metric, the radiation power can be derived from the potential as the imaginary part of the term
, but we will deal with it in a separate paper.
7.2. Taking into Account Self-Rotation of the Participating Masses m1 and m2
Self-rotation with the Kerr ansatz
The Kerr ansatz is valid as long as there is a
-symmetry, i.e. the system is independent of
.
If there is self-rotation (around z-axis) for the masses m1 and m2, the θ-symmetry is not disturbed and self-rotation (spin-) angular momentums L1 and L2 add up and contribute to the Kerr parameter α according to the formula
. For a rotating black hole
, where
is the Schwarzschild radius, κ is the inertia-factor (κ = 2/3 for a spherical shell) and
the angular frequency.
Now we can make the ansatz that in the case of a binary rotator with self-rotation of the black-holes with masses m1 and m2, this contribution adds a part to the (real)
of the binary rotator, where from above
.
For the masses (m1, m2) (dimensionless, i.e. m = 1) rotating with angular frequencies (ω1, ω2) we get then
For
the rotator becomes a single rotating Kerr- blackhole with
, as it should be. Also, the contribution
from gravitational rotation becomes zero for
, and the spacetime becomes the normal real Kerr spacetime of a rotating black-hole.
We can generalize the self-rotation to arbitrary rotation axis: the z-component of the axis is yields the (real) Kerr parameter
as above, the perpendicular component, i.e. the projection on the rotator plane yields the imaginary Kerr parameter
, so that the total Kerr parameter becomes
With this ansatz, the solution
becomes complex, and we simply take the absolute of the variables
, as all of them must be positive. The numerical solution of the radial equation with complex coefficient, as well as that of the full orbit equations with complex coefficient presents no additional difficulties, so the ansatz is perfectly feasible.
The program code for the calculated examples can be found in [2].
Self-rotation with Manko-Ruiz rotator
In case of the Manko-Ruiz rotator, the metric parameter a is determined from the orbit reduced angular momentum
according to the formula
(see 6.5 (27)).
If we take into account individual rotation of the two rotator masses, then the additional Kerr parameters
add up to the relative angular momentum
directly:
. If the individual rotation is inclined under angle
to the rotator (i.e. orbit) z-axis, then the corresponding perpendicular components enter the formula as imaginary numbers, as described above:
We get then, as above, orbit equations with complex coefficients, solve them in complex orbit radius
, and the physical radius is then taken as the absolute value
.
8. Numerical Examples
In the following subsections 8.1 and 8.2, we present orbits around the binary rotator calculated numerically from different models.
In subsection 8.3, we present results of calculation of rotator orbits in its own metric for mean radius
, and
.
The details of the calculation, the results and graphics are accessible in the Mathematica source code [2].
8.1. Binary Rotator with Mass Ratio μ = 2, Observer Orbit
We consider an example of a rotator consisting of a binary black-hole with mass ratio
, in dimensionless spacetime, i.e. r in Schwarzschild-units rs, and t in Schwarzschild-units
.
We consider an observer orbit around the rotator, with orbit parameters
, or equivalently
.
For the Kerr radial equation, we have to specify the initial value for
, the orbit parameters
,
, and the rotator parameters
and
. Here we use the Kerr parameter value derived for the Newtonian
correction in 5.2:
, where
is the mean distance of the
binary rotator. We have two scales here:
determines the maximum orbit radius (orbit scale),
is the mean rotator distance (rotator scale).
We use the following values:
,
,
,
,
,
, relative absolute energy
, eccentricity
,
pure Newtonian correction, resp. complex
with additional self-rotation
for both black-holes 45˚ to rotator axis
.
The reduced angular momentum
is set on purpose relatively low in order to make eccentricity large, to test the validity of the formula for
for non-circular orbits: the results show that the approximation of the Newtonian correction is still good in this case.
In the case of the Manko-Ruiz spacetime, we have made the adaptation of the orbit parameters by semi-manual gradient procedure.
In this way [2], we got the Manko-Ruiz parameter values with fixed mass ratio
,
,
.
The result values are:
,
,
,
,
.
We use the following calculation models without self-rotation: pure Newtonian N1, Newtonian with rotator correction N1s, Schwarzschild S1, Kerr approximation (E1dS0 Schwarzschild-Kerr) K1, Kerr approximation (E1dA0 full Kerr-approx.) L1, exact Kerr (full Kerr orbit equations) A1, the Manko-Ruiz 4-parameter Kerr-extension M1.
We use the following calculation models with self-rotation: Kerr approximation (E1dS0 Schwarzschild-Kerr) K1s, Kerr approximation (E1dA0 full Kerr) L1s, exact Kerr (full Kerr orbit equations) A1s.
8.2. Numerical Results
Orbit parameter table
The following table (Table 1) contains the principal orbit values for the models without self-rotation: minimal and maximal radius rmin, rmax, radius drift Δr, first rotation period T1, period drift ΔT, all in Schwarzschild units.
Table 1. Orbit parameters of observer orbits in the rotator field in different models.
Orbit radius plots r(τ)
Kerr exact, Kerr apS Schwarz-app, Schwarzschild, Newton
Here we present models, which approximate well the uncorrected original pure Newton orbit with values: rmin: ~22, period T1 ~4400. The main deviation in period is in the model Kerr apSas as shown in Figure 3.
Kerr apF full-Kerr-approx., Schwarzschild, Newton corrected, Manko-Ruiz
Here we present models as shown in Figure 4, which approximate well the rotator-corrected Newton orbit with values: rmin: ~34, period T1 ~4470, and for comparison Schwarzschild orbit with rmin = 21.5. The main deviation in period is in the Schwarzschild model.
Kerr exact, Kerr exact with spin
Here we show the influence of 2x spin (self-rotation)
added to the orbit angular momentumas shown in Figure 5, where the spin is inclined 45˚ to rotator axis. The exact Kerr solution shows only a small deviation
and
.
Kerr apS Schwarzschild-approx, Kerr apS with spin
Here we show the influence of self-rotation on the Kerr solution with Schwarzschild-approximationas shown in Figure 6.
The Kerr solution shows a larger deviation
and
.
Kerr apF full Kerr-approx, Kerr apF with spin
Here we show the influence of self-rotation on the Kerr solution with full Kerr-approximation.
The Kerr solution shows a larger deviation
and
.
8.3. Manko-Ruiz Rotator Orbits in Its Own Field (Metric)
We consider the Manko-Ruiz rotator in its own field, with parameters
, relative orbit energy
for
, and reduced angular orbit momentum
, where the two latter relations have been derived in subsection 6.5, we set the total mass
and calculate in Schwarzschild units, i.e. dimensionless.
Figure 3. Orbit radius r(τ): Kerr exact, Kerr apS = Schwarz-app, Schwarzschild, Newton.
Figure 4. Orbit radius r(τ): Kerr apF full-Kerr-approx., Schwarzschild, Newton corrected, Manko-Ruiz.
Figure 5. Orbit radius r(τ) without and with spin: Kerr exact.
A rotator orbit for a rotator with given mass ratio
is determined by the mean radius
(which determines the relative orbit energy according to
) and the reduced angular momentum
(which determines the eccentricity of the orbit). We extract the metric parameters from it:
, and a from the above formula
.
We calculate for the rotator with
two orbits with mean radii
resp.
and
.
We calculate also the corresponding Newton and Schwarzschild orbits with the same orbit parameters
,
.
The result are presented in Schwarzschild units in the following Table 2 and Figure 7, resp. Table 3 and Figure 8.
Comparison of the two tables reveals strong differences in behavior of the Manko-Ruiz orbits and Kepler orbits as shown in Figure 9.
Figure 6. Orbit radius r(τ) without and with spin: Kerr apS Schwarzschild-approx.
Figure 7. Orbit radius r(τ) without and with spin: Kerr apF full Kerr-approx.
Figure 8. The Manko-Ruiz
orbit and corresponding Newton orbit.
Figure 9. The Manko-Ruiz
orbit and corresponding Newton and Schwarzschild orbit.
Table 2. Orbit parameters of observer orbits in its own field,
.
Table 3. Orbit parameters of observer orbits in its own field,
.
The Kepler law
does not apply anymore, we have
.
Also, the formula for the eccentricity is different from the Keplerian
: e is larger for Manko-Ruiz, and the energy-radius relation
deviates from the Newtonian case.
However, we have to keep in mind that in the limit
(i.e. here
) Manko-Ruiz metric becomes the first Schwarzschild and then flat Minkowski metric and the orbit become the first Schwarzschild and then Newtonian orbits.
9. Conclusion
The central theme here is the gravitational two-body problem.
New results are:
- Precise formulation and calculation
The two-body problem is precisely formulated and calculated as (single) rotator with a single mass based on the Newtonian, Schwarzschild, and Kerr effective potential
in the energy (radial) equation.
- Adaptation of the Kerr and Manko-Ruiz spacetime for observer orbit
Kerr and Manko-Ruiz spacetime are adapted to the Newtonian rotator correction for the observer orbit, i.e. parameters are fitted resp. calculated in order to describe the Newtonian correction in order
.
In the case of the Manko-Ruiz spacetime, the Newtonian-limit-principle in GR guarantees that with adaptation to the Newtonian correction, the Manko-Ruiz spacetime
for observer orbit is the exact GR-solution for an orbit in the rotator potential, since it is a solution of the Einstein equation with exactly 4 independent parameters matching the 4 parameters
of the two-body problem, and becomes the Newtonian solution in the limit
.
- Calculation of rotator orbit in its own field
The verification of Manko-Ruiz Einstein equations, the symbolic calculation of the Manko-Ruiz general effective potential
, its development in
-powers are new results.
Furthermore, we derive a formula for the reduced angular momentum
of a Manko-Ruiz rotator in its own field from the rotational parameter a, the orbit radius R, and the mass ratio μ, and for the relative total energy
. These two parameters determine the rotator orbit, and we calculate this orbit in two examples and present the results in table and graphical form.
- Numerical calculation of Manko-Ruiz orbits and of the effective potential
The presented numerical calculation of observer orbits in the Newton-adapted Manko-Ruiz potential is also a new result.
Appendix A1
Schwarzschild spacetime in matrix form (dimensionless, r in rs units)
Kerr spacetime in matrix form (Boyer-Lindquist coordinates), Kerr parameter
Christoffel symbols
(Schwarzschild) have the values