Binary Star System Decay by Graviton Interaction

The action of gravitons in a binary star system is modelled as the locus of points on an ellipse synchronous to the elliptic orbit of the binary star. In their interaction between the masses in the system the rotational energy of the gravitons is reduced by gravitational redshift, which accounts for the decay of the binary star orbital period. This model is able to fit a broad range of eccentricities of binary pulsar orbits and orbital period decay comparable to the General Relativistic gravitational wave model. the graviton model with the GR gravitational wave model and show how our model mathematically corresponds to it, although some physical aspects are disparate.


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 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 where r is the distance from the center of mass to the origin ellipse, dφ is the small angle swept out in time dt , a and b 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) where the graviton specific angular momentum 2 d dt φ = h r and ( ) 2 2 3 2 1 ε − = l ab , where l 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 r is given by The graviton origin ellipse semi-latus rectum l is defined in terms of the Schwarzschild radius s R of the combined masses, where the correction function ( ) σ ε is defined, 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, r 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, 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 t δ under the field acce- , the graviton net rotational energy will change due to gravitational redshift by the amount, , where we substituted from (7) for Ξ and the change in the free fall velocity where the minus sign implies a reduced (redshifted) graviton energy 1 , since the gravitons are moving in the same direction as the velocity v δ , and where the distance r is given by, 1 We remark that in the original definitions given in (7) and (9), the graviton rotational mass was defined as the reduced mass GMm M m r + , but these have now been defined in the physically correct form.

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, which upon differentiating with respect to the time t and simplifying yields the rate of change of the orbital period, where the total orbital energy ( ) . Substituting from Equations (3) and (5) into Equation (7) we get the graviton rotational energy, Dividing Equation (9) by t δ and substituting for Ξ from (13) gives, Then, substituting d d .
In making fits to binary star data we will take an average of the radial distance r where (15). Finally, using Kepler's third law, (11), we substitute for A in terms of T into Equation (15) and simplify to obtain, 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.

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, is the reduced mass, r is the separation of the masses and ij I is a 3 3 × traceless matrix, given by, The strain due to the wave caused length change ij d δ at a distance L d from the system [4] is given by, At a sufficient distance L d from the source we can use a linearized approximation of Einstein's equations. The rate of change of the strain is given by where we have used the dot notation for the time derivative.
and the integration over the sphere gives an area of 2 4 L d π , cancelling that value in the numerator of (21). From (17) and (18) and by substituting this result into (21) and simplifying yields, Using the notation 2 b P T ω = = π for the orbital period and r in place of A for the circular orbit in (12) and substituting for d d E t with the negation of (24) because the GW emission causes an energy loss to the orbit, we get, 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,

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 Notice that the correction function ( ) 1 σ ε > , always being above the horizontal value 1 line in Figure 2. From this fact it can be shown that the semi-major axis a and the semi-minor axis b of the graviton origin ellipse are always greater than the Schwarzschild radius s R of the combined masses M m + .
Substituting these experimental parameters and astronomical values for B1913+16 into (16) yields, 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    where p P is the predicted period decay and i P is the measured intrinsic period decay. Excluding PSR J0621-1002 which has a poor intrinsic b P 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 b P values and this paper's predicted d d T t values is d d 0.0608 0.0084 T t MeanErr = ± . For a comparison with the standard GR GW emission theory, the mean error and unbiased standard deviation of the mean error is 0.1097 0.0231 bGR P MeanErr = ± . In Table 4 and Table 5 we present predictions for two binary star systems [15] for which the intrinsic b P 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.

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 ω Ξ =Mh 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 M . 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 phenome-nological 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) where M 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 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
Then, substituting for l from (5)

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 ( ) 4 2 1 0.1 86400 1.2 10 Hz f − < × ≈ × . 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.

Conflicts of Interest
The author declares no conflicts of interest regarding the publication of this paper.