General Relativistic Orbital Effects in Compact Binary Stars ( Solution by the Method of Celestial Mechanics )

Perturbation methods are employed to calculate time variation in the orbital elements of a compact binary system. It turns out that the semi-major axis and eccentricity exhibit only periodic variations. The longitude of periastron and mean longitude of epoch exhibit both secular and periodic variation. In addition, the relativistic effects on the time of periastron passage of binary stars are also given. Four compact binary systems (PSRJ0737-3039, PSR1913+16, PSR1543+12 and M33X-7) are considered. Numerical results for both secular and periodic effects are presented, and the possibility of observing them is discussed.


Introduction
In the wake of unceasing development in the post-Newtonian celestial mechanics, at present, the research on the post-Newtonian effects has been exhibited gradually due to the fact that the degree of accuracy of astronomical observation improves unceasingly.Hence several authors devoted their research to this subject and the scopes (Brumberg, (1972(Brumberg, ( , 1985) ) [1] [2], Rubincam (1977) [3], Soffel (1987Soffel ( , 1989) ) [4] [5], Iorio (2005) [6]).These authors research mainly the post-Newtonian effect of the orbital elements of planets and artificial satellites in the solar system.The largest post-Newtonian effects have been exhibited in the orbits of binary systems, especially in the compact binary systems.So, the research in the post-Newtonian orbital effect of the compact binary stars is important and Hence some authors research this subject from theory and observation.Such as, Will (1981Will ( , 2006) ) [7] [8], Damour & Deruelle (1985, 1986) [9] [10], Schāfer & Wex (1993) [11], Wex (1995) [12].Calura (1997) [13], Iorio (2007) [14] studied this subject from theory.On the other hand, Burgay et al. ( 2003) [15], Konacki (2003) [16], Kramer et al. (2005) [17], and Weisberg & Taylor (2005) [18] researched the post-Newtonian effect for the periastron shift of compact binary stars (PSRJ0737-3039, PSR1913+16 PSR1543+12) from observation.These researches are very interested in the theoretical and observing aspects.This paper presents the post-Newtonian effect on all orbital elements of the compact binary stars on the theoretical aspect.The research method of this paper is different from previous studies that this paper uses the method of perturbation theory in the celestial mechanics.

R, S and W Components for the General Relativistic Accelerations in the Two-Body Problem
The relative acceleration of two-body with the Post-Newtonian Parameters is given by Will (1981) [7] ( ) here f denotes the true anomaly.n is a unit vector in the radial direction and λ are unit vectors in the orbital plane.n is directed along the radial direction, and λ is perpendicular to n .In the equation m denotes Gm and the right side should be multiplied by c −2 .G is gravitational constant and c is the speed of light.
In this paper we research the general relativistic effect.In the general relativistic case the Post-Newtonian parameters (Will 1981)  [7].The Equation (1) can be written Using the relative expressions the Equation (3) may be written here boldface denotes vector.
We resolve the acceleration a into a radial component Rn , a component Sλ , normal to Rn and a component W normal to the orbital plane.
i.e, ( ) (the unit vector normal to the orbital plane).
On comparison with the expression (4), we get three scalar accelerative components R, S and W ( ) Substituting the following formulas of the problem of two body into the above formula (Smart, 1953) We obtain 1988) [21] by using the second formula of the Equations ( 6), we get the perturbation equations with an independent variable true anomaly f ( ) where ω is the longitude of periastron, E is the eccentric anomaly and 0 ε is the mean longitude at epoch.

Integration for the Perturbation Equations and Its Perturbation Solutions
Substituting R, S and W for expressions (7) into the perturbation Equation ( 8) by using Kepler's third law ( ) Integrating the above Equation ( 9), we obtain the perturbation secular and periodic solutions where λ denotes the mean longitude of periastron, E denotes the eccentric anomaly.In the last integral expression, we have used the next integral already:

The Secular and Periodic Variation of the Orbital Elements
1) The secular variation per cycle (revolution) By letting 0 0, f = 2π f = , 0 0 E = , 2π E = , the periodic terms are disap- peared and one obtains the secular variables per cycle (revolution): The author Li (2010) [24] obtained the formulas for the post-Newtonian effect on the time variation of periastron passage of binary stars.In that paper we change the symbol of the first expression of (35) (in the case of relativity) as the symbol of the present paper, which is ( ) 2) The secular variable rate: where the period P is denoted in unit of day 3) The periodic variation of amplitudes: In the expressions (10) all terms are the periodic variable terms except for all secular terms.Here we list the maximal and minimal amplitudes of the periodic terms for semi-major axis, a and eccentricity, e from the expressions (10).

Numerical Results for Four Compact Binary Systems
We use the formulae ( 11) -( 13) to calculate the secular of the general relativistic secular effect on the orbital elements of four compact binary systems.It is convenient to reduce the formulas ( 11) -( 13) to practical units m 1 , m 2 and a are denoted by the unit in solar mass here This paper chooses four compact binary systems: PSR1913+16, PSR1543+12, PSRJ0737-3039 and a black hole M33 X-7 as an example.For these compact binary stars, their data for P, a, e, M and m are retrieved from Burgay et al.  2007) [23].Their data are listed in Table 1.
Substituting these data in Table 1 into formulas ( 16) and ( 13) and ( 14) -( 15), we obtain the numerical results for the periodic and secular variation of the orbital elements of four compact binary stars in Tables 2-4.
Table 1.The data of four compact binary stars.It can be seen from Table 3 that the maximum secular variable of longitude per period is PSRJO737-3039.The longest time of the periastron passage is M33-7.
It can be seen from Table 4 that the maximum secular variable rates of longitude is PSRJ0737-3039.The longest time of the periastron passage is M33X-7.

Discussion and Conclusions
1) The comparison of the theoretical results with the observable results.
The theoretical results in this paper as compared with the observable results given by several authors for three compact binary stars are listed in the Table 5.
It can be seen from the above Table 5 the theoretical results are very close to the observed results.2) The possibility of observing effects.
In the solar system the advance of perihelion of Mercury may be observed by the recent instrument.We can see from Table 4 that maximal value of the advance of PSRJ0737-3039 is 61016" per year which correspond to 145000" time value of advance of perihelion of Mercury in solar system.The value of the advance of periastron of compact binary star is largest than that of Mercury in solar system.Therefore the effects of the advance of compact binary stars can be observed too.
We have four conclusions: a) The compact binary stars are the best objects for studying the post-Newtonian effects on the orbits b) Although there are no secular variation for the semi-major axis and eccentricity, there is maximal amplitude of the periodic terms for semi-major axis, such as max 1661.80 km A = .
c) The longitudes of periastron and the mean longitude at epoch exist both secular and periodic variable terms, and the maximal values for ω and 0 ε ar- rive at 16.95/yr and −29.68/yr for PSRJ0737-3039 respectively.
d) The longest time of periastron passage is 19.47 minute per year for black hole M33 X-7.This corresponds to over 3 seconds per a day.
How to cite this paper: Li, L.S. (2017) General Relativistic Orbital Effects in Compact Binary Stars (Solution by the Method L. S. Li DOI: 10.4236/wjm.2017.712027361 World Journal of Mechanics meaningful.

Table 2 .
Periodic variation of the amplitudes of the orbits of semi-major axis and eccentricity.It can be seen from Table2that the maximum amplitude of semi-major axis is the black hole binary system M33 X-7 and the maximum eccentricity is PSR 1913 + 16.

Table 3 .
Secular variation of the orbital elements of four compact binary stars per cycle (Revolution).

Table 4 .
Secular variable rates of the orbital elements of four compact binary stars per year.

Table 5 .
Comparison of the theoretical results with the observable results.