1. Introduction
The purpose of this paper is to introduce a model for the action of the gravitons in a binary star system described in [1], a paper which has a section on binary pulsars based on a comparison to the equation for the decay of a binary star system developed in [2], but we did not give a complete model for the action of the gravitons. This new model describes the rotational phase of gravitons which interact between the masses, where the graviton phase is a locus of points on an ellipse around the central mass, with an eccentricity and period identical to the binary star system ellipse. In this model the gravitons interact directly between the masses in orbit, losing energy in the system by gravitational redshift. We will reanalyze the set of binary pulsars from the previous paper as well as make comparisons to the General Relativity (GR) model. We will also give predictions for two millisecond binary pulsars. We make a detailed comparison of the graviton model with the GR gravitational wave model and show how our model mathematically corresponds to it, although some physical aspects are disparate.
2. Rotational Energy of Gravitons in a Binary Star System
Consider a binary star system composed of a companion star of mass m in orbit around a primary star of mass M. By Kepler’s second law, the companion star in an elliptical orbit will sweep out an equal area in equal time around the primary star, expressed by
(1)
where r is the instantaneous distance between the stars,
is the angle r makes with the line thru the stars at perigee, A and B are the semi-major and semi-minor axes, respectively, and T is the orbital period. We assume that gravitons are responsible for the force between the stars. We define the graviton origin ellipse which encloses the center of mass of the system and has the same eccentricity and the same angular speed as the binary star ellipse. The origin ellipse describes the motion of the line of gravitons interacting between the two stars in the system, having the same period T as the orbit ellipse and expressed by
(2)
where
is the distance from the center of mass to the origin ellipse,
is the small angle swept out in time
,
and
are the semi-major and semi-minor axes of the origin ellipse, respectively, and T is the orbital period of the binary star system. Multiplying Equation (2) by the angular frequency
and simplifying we get,
(3)
where the graviton specific angular momentum
and
, where
is the semi-latus rectum of the graviton origin ellipse and
is the ellipse eccentricity, equal to the binary star system eccentricity. The radial distance
is given by
(4)
The graviton origin ellipse semi-latus rectum
is defined in terms of the Schwarzschild radius
of the combined masses,
(5)
where the correction function
is defined,
(6)
where the parameters
,
, p and q are determined by experiment and c is the speed of light in vacuum. Figure 1 shows the relationship of the binary star ellipse and the graviton origin ellipse.
Figure 1. Binary star system ellipse (outer) and graviton origin ellipse (inner). M is the primary mass, m is the companion mass, r is the distance between the masses,
is the distance from M to the graviton rotational energy point of interaction between the masses and f is the rotation angle of the ellipses. The ellipses are not drawn to scale.
We define the graviton rotational energy
in the center of mass system of the masses M and m,
(7)
where
(8)
is the graviton rotational relativistic mass (not a rest mass since gravitons travel at speed c) and where
is the reduced mass of the system. Since the gravitons having rotational energy
are in free fall in the gravitational field of the binary star system, over a small time
under the field acceleration
, the graviton net rotational energy will change due to gravitational redshift by the amount,
(9)
where we substituted from (7) for
and the change in the free fall velocity
where the minus sign implies a reduced (redshifted) graviton energy1, since the gravitons are moving in the same direction as the velocity
, and where the distance r is given by,
(10)
3. Rate of Change of the Orbital Period
To obtain the time rate of change of the orbital period T, we use Kepler’s third law,
(11)
which upon differentiating with respect to the time t and simplifying yields the rate of change of the orbital period,
(12)
where the total orbital energy
. Substituting from Equations (3) and (5) into Equation (7) we get the graviton rotational energy,
(13)
Dividing Equation (9) by
and substituting for
from (13) gives,
(14)
Then, substituting
from (14) into Equation (12) gives,
(15)
In making fits to binary star data we will take an average of the radial distance r where
for
. Thus, the radial distance
in (15). Finally, using Kepler’s third law, (11), we substitute for A in terms of T into Equation (15) and simplify to obtain,
(16)
Equation (16) is equivalent in form to the derivation given by [2], which was derived from General Relativity with gravitational wave (GW) emission for energy decay. We describe this GW relation in the next section and compare the two methods in a subsequent section.
4. The Energy Loss Due to Gravitational Wave Emission
We give a brief derivation of binary star orbital decay due to GW emission based on [3]. Gravitational waves are emitted by a binary star system due to the time rate of change of the quadrupole moment for the binary, where the quadrupole moment for the simple case of a circular orbit is expressed by,
(17)
where
is the reduced mass, r is the separation of the masses and
is a
traceless matrix, given by,
(18)
The strain due to the wave caused length change
at a distance
from the system [4] is given by,
(19)
At a sufficient distance
from the source we can use a linearized approximation of Einstein’s equations. The rate of change of the strain is given by
(20)
where we have used the dot notation for the time derivative. Then, the power (luminosity)
of the wave is proportional to the square of
integrated over the surface of the volume of space, expressed by
(21)
where
(22)
and the integration over the sphere gives an area of
, cancelling that value in the numerator of (21). From (17) and (18) the summation in (21) becomes
(23)
and by substituting this result into (21) and simplifying yields,
(24)
Using the notation
for the orbital period and r in place of A for the circular orbit in (12) and substituting for
with the negation of (24) because the GW emission causes an energy loss to the orbit, we get,
(25)
where we used Kepler’s law to transform
. Putting in the dependence on the orbit eccentricity
from [2] for an elliptical orbit, we put (25) into the familiar form,
(26)
Application to PSR B1913+16 and Other Binaries
We look at the report on B1913+16, the Hulse-Taylor binary pulsar [5]. This astronomical endeavor spanned 30 years of approximately yearly observations of the binary system. The data for the system is as follows: primary mass
, companion mass
,
,
,
, where M is the primary mass and m is the companion mass,
is the solar mass,
is the binary orbital period,
is the orbital period change and
is the orbital eccentricity. The initial orbital period for B1913+16 is
. For the gravitational constant we use the value
and
for the speeds of gravity (graviton) and light (photon) in vacuum. To make a good fit of (16) to the binary pulsar systems we are examining, we determined best fit parameter values of
,
,
and
, which makes the leading factor
, where the previous study [1] had a leading constant factor of
. Figure 2 shows the correction function
with the specified parameters
and the GR correction function from (26),
(27)
Notice that the correction function
, always being above the horizontal value 1 line in Figure 2. From this fact it can be shown that the semi-major axis
and the semi-minor axis
of the graviton origin ellipse are always greater than the Schwarzschild radius
of the combined masses
.
Substituting these experimental parameters and astronomical values for B1913+16 into (16) yields,
(28)
which is a good match to the experimental corrected value [5] of
(29)
and is also close to the GR theoretical value,
(30)
In Table 1 and Table 2 we list nine PSR’s, [5] - [13] and [14]. In Table 3 we show the error in each prediction, the error computed by
Figure 2. Correction function (solid line) with best fit parameter values of
,
,
and
. GR correction function (dashed line). Horizontal line drawn at 1 for reference.
Table 1. PSR binary systems studied. Column description: 1: PSR name; 2: Primary (pulsating) star mass; 3: Companion star mass; 4: Period (day) Pb or T; 5: Orbit eccentricity.
Table 2. PSR binary systems studied (continued). Column description: 1: PSR name; 2: dPb/dt (observed intrinsic value); 3: dPb/dt (general relativity computation); 4: dT/dt (due to graviton gravitational redshift).
Table 3. Prediction errors of binary systems studied. Column description: 1: PSR name; 2: dPb/dt (general relativity computation error); 3: dT/dt (due to graviton gravitational redshift error). The errors are computed by Equation (31).
(31)
where
is the predicted period decay and
is the measured intrinsic period decay. Excluding PSR J0621-1002 which has a poor intrinsic
measurement, for the PSR prediction errors given in Table 3, the mean error and unbiased standard deviation of the mean error between the observed intrinsic
values and this paper’s predicted
values is
. For a comparison with the standard GR GW emission theory, the mean error and unbiased standard deviation of the mean error is
.
Table 4. PSR binary systems with new predictions. Column description: 1: PSR name; 2: Primary (pulsating) star mass; 3: Companion star mass; 4: Orbit eccentricity.
Table 5. PSR binary systems with new predictions. Column description: 1: PSR name; 3: dPb/dt (general relativity computation); 4: dT/dt (due to graviton gravitational redshift). dPb/dt (observed intrinsic value) not yet determined.
In Table 4 and Table 5 we present predictions for two binary star systems [15] for which the intrinsic
are not yet determined: PSR J1949+3160, a millisecond pulsar and a white dwarf companion, and PSR J1950+2414, also a millisecond pulsar with a possible low mass white dwarf companion. The estimated deviations of the predicted values are determined using the standard deviation of the mean error found in the preceding analyses.
5. Comparison of the Methods
This paper’s approach attempts to use the gravitational redshift mechanism as the cause of the orbital decay found in binary star systems, a mechanism which has no GW emission. Our model emulates the traditional (GW) equations. The goal was to have the graviton energy redshift during free fall in the gravitational field of the binary star account for the observed orbital decay of the binary system. We struck upon the idea of a circulating equation
for the gravitons which surround the nucleus of the binary system, where the position of the orbiting body is tracked in phase by the graviton rotational energy of relativistic mass
. This graviton energy travels at velocity c from the central mass toward the orbiting body and with respect to the frame which is in free fall in the field between the masses, the graviton energy will be redshifted, thus reduced in energy. This mechanism for energy loss is a relativistic effect without emission, just as the case where light loses energy by gravitational redshift of its frequency when traveling away from the surface of a star. This is contrary to the mechanism of GW emission due to rotating binary stars where the wave carries away energy from the system.
Although the addition of a correctional function
was found necessary to make a better equality of the graviton model equation to the astrophysical data, this is a brute force approach even though this model strives to be compliant with the theory of General Relativity as formulated by Einstein, more precisely to its linearized approximation. At this stage, having compared our phenomenological approach to the standard GR method, we can go a step further in our graviton theory by actually equating it to the GR wave emission equation, giving to that equation a new interpretation of a gravitational redshift energy loss which emits no radiation.
Equate (9) and (21) in the form
, expressed by,
(32)
where
is the graviton rotational relativistic mass defined in (8). Moving the acceleration rate from the left side to the right side of (32), the graviton energy
surrounding the nuclear mass
takes the form,
(33)
From (3) and (23) with the dependence on the orbit eccentricity
from [2] for an elliptical orbit, which also converts radial distance r to semi-major axis
, substituting all this into (33) we get,
(34)
Now, using Kepler’s third law for
in (34), with some simplification, yields,
(35)
Then, substituting for
from (5) into (35) and simplifying we get,
(36)
Finally, from (36) we are left with the relationship of the correction function
to the GR correction function
,
(37)
where
is given by (27). Substituting (37) for
into (16) we obtain an orbital period decay equation from our graviton model identical to the standard GW model (26).
6. Conclusions
We presented a model to describe the rotational energy of gravitons in a binary star system. We have based it on an analogy to the General Relativity theory GW energy loss in a binary star system. Our approach enabled an adaptive correction to the graviton rotational energy which minimized the error in the prediction of the decay of the orbital period against the experimental value. Although the gain in accuracy is just a 5% reduction in the prediction error compared to GR, it does suggest that a graviton theory is a viable approach.
There are ongoing projects in the search for continuous low frequency GW’s from known binary pulsars that are far from merging, where the orbital periods are of order > 0.1 day, and the GW frequencies will be
. The LIGO/Virgo detectors can only go down to about 20 Hz [16]. Analysing day long or week long recorded data streams lowers the detection frequency to the range where these continuous GWs could be detected [17], but no detections have been made thus far. Research into continuous GW emission from binary pulsars is ongoing and the LISA telescope will look for these types of signals directly in the 0.1 mHz to 1 Hz range [18]. In a related field, regarding continuous GW emission from isolated neutron stars (pulsars), the latest research in the data of the LIGO/Virgo third observing run (O3) have made no detections [19]. Although there have been GW detections made of relativistically high energy binary blackhole mergers and binary neutron star mergers, these are not the type of low speed binary star events addressed in this discourse. As there has as yet been no detection of gravitational waves from these low speed sources, we deem this as ample justification to consider the thesis we have put forth.
Thanks
We thank the reviewer who has made challenging suggestions for the improvement of this work.
NOTES
1We remark that in the original definitions given in (7) and (9), the graviton rotational mass was defined as the reduced mass
and the gravitational acceleration field was defined as
, but these have now been defined in the physically correct form.