a detailed knowledge about the dynamics of spinning compact binaries is of interest to both astrophysicists and theoretical physicists .
the astrophysical interest arises from the highly accurate radio observations of relativistic binary pulsars , which can be used to test general relativity in strong field regimes @xcite . moreover , inspiralling spinning compact binaries are promising sources of gravitational radiation for the ground based and proposed space based laser interferometric gravitational - wave detectors @xcite .
the theoretical interest is due to the possibility of constructing numerically spinning compact binary configurations in full general relativity [ see ref .
@xcite and references therein ] .
further , it is of some interest to explore if the dynamics of spinning compact binaries will be deterministic or not @xcite .    in many of the above mentioned scenarios ,
the dynamics of spinning compact binaries can be described in great detail by the post - newtonian ( pn ) approximation to general relativity .
the pn approximation allows one to express the equations of motion of a compact binary as corrections to newtonian equations of motion in powers of @xmath1 , where @xmath2 , @xmath3 , and @xmath4 are the characteristic orbital velocity , the total mass and the typical orbital separation of the binary . for spinning compact binaries , the dynamics is determined , not only by the orbital equations of motion for these objects , but also by the precessional equations for the orbital plane and the spin vectors themselves @xcite .
note that for spinning compact binaries , spin effects enter the dynamics for the first time through spin - orbit interaction .
therefore , the hamiltonian describing the leading order spin - orbit interaction can be added to the pn accurate conservative hamiltonian which gives the dynamics of nonspinning compact binaries .
it is highly desirable to obtain keplerian - type parametric solution to the conservative dynamics of spinning compact binaries .
this is so as 1pn accurate quasi - keplerian parametrization , obtained by damour and deruelle @xcite , is heavily employed to analyze radio observations of relativistic binary pulsars and to test general relativity in strong field regimes @xcite .
moreover , 2pn accurate keplerian - type parametric solution for nonspinning compact binaries , obtained in refs .
@xcite , is heavily employed to construct highly efficient `` ready to use '' search templates for nonspinning compact binaries moving in inspiralling eccentric orbits @xcite .
recently , 3pn accurate generalized quasi - keplerian parametrization for the conservative dynamics of spinning compact binaries in eccentric orbits was obtained in refs .
more precisely , the pn accurate parametric solutions are available for the following two distinct cases : ( i ) the binary is composed of two compact objects of equal mass and two arbitrary spins , and ( ii ) the two objects have unequal masses but only one object has a spin . the pn accurate parametric solution may be used to describe the following two astrophysical situations . as the measured neutron star ( ns )
masses in ns
ns binaries usually cluster around @xmath5 , case ( i ) may be used to describe such ns
ns binaries .
the case ( ii ) may be applicable in the instance of a neutron star orbiting a much larger , rapidly rotating black hole . in that case
, the black - hole spin will dominate the orbital precession , and the neutron star spin may be neglected .
the above parametrization , for the two cases considered , shows that the dynamics is integrable and can not be chaotic @xcite .
the parametric solution will be employed to construct `` ready to use '' search templates of spinning compact binaries in inspiralling orbits @xcite .
further , it will also be useful to construct a fully 2pn accurate `` timing formula '' , required to analyze binary pulsar observations by the square kilometre array ( ska ) @xcite .    in this paper , we compute , using our parametrization , the observable dimensionless relativistic periastron advance parameter @xmath0 , relevant for binary pulsars .
the expression for @xmath0 that we derive will be applicable , like the parametrization computed in ref .
@xcite , only for the two distinct cases mentioned above .
our expression for @xmath0 is in agreement with the one obtained by damour and schfer @xcite .
it should be noted that damour and schfer employed the evolution equation for the laplace - runge - lenz vector , computed from a hamiltonian , to obtain the spin - orbit contributions to the parameter @xmath0 . in this paper
, we employ our generalized quasi - keplerian parametrization , derived in ref .
@xcite , to compute the parameter @xmath0 .
there are two reasons that motivated us to do the current computation .
the 2pn accurate expression for @xmath0 , derived by damour and schfer , valid for arbitrary masses and spins , is extremely important due to the discovery of the double pulsar psr j07373039 @xcite .
this is because the accurate timing of the double pulsar over a period of few years should lead to the measurement of the moment of inertia of one of the pulsars @xcite .
it is from the measurement of @xmath0 to 2pn order , which includes leading order spin - orbit contributions , that one will be able to measure the moment of inertia of the fast rotating pulsar in this double pulsar @xcite .
the measurement of the moment of inertia of a neutron star for the first time will lead to strong constrains on its equation of state @xcite .
we compute @xmath0 , as mentioned earlier , using keplerian - type parametrization with the aim of presenting an alternative way of obtaining this observable quantity .
a second reason for computing @xmath0 is the existence of certain intriguing features appearing in our parametrization .
we observed that the factor that modulates the true anomaly , due to spin - orbit interaction , contains a contribution , independent of the spins themselves , and of 1pn order [ see eq .
( 5.6k ) in ref .
we want to make sure that this feature will not appear in the observable parameter @xmath0 , computed using our generalized quasi - keplerian parametrization .
the paper is structured in the following manner . in sec .
[ sec2 ] , we briefly describe the 3pn accurate keplerian - type parametrization for the conservative dynamics of spinning compact binaries , moving in eccentric orbits , when spin effects are restricted to the leading order spin - orbit interaction , derived in ref .
@xcite , and point out its salient features .
[ sec3 ] and sec .
[ sec4 ] deal with the issues related to the derivation of the parameter @xmath0 .
we show for the first time how to derive @xmath0 using the parametrization , described in sec .
[ sec2 ] , after describing how damour and schfer obtained @xmath0 .
we summarize our results in sec .
[ secconclusions ] .
appendices [ appendix_binary_geometry ] and [ appendix_lrl_vector ] provide details of our computation and consistency checks .
recently , 3pn accurate generalized quasi - keplerian parametrization for compact binaries in eccentric orbits with the first - order spin - orbit interaction was obtained in refs .
@xcite for the following two distinct cases : ( i ) the binary is composed of two compact objects of equal mass and two arbitrary spins , and ( ii ) the two objects have unequal masses but only one object has a spin .
the underlying dynamics is fully specified by the following pn accurate ( reduced ) hamiltonian @xmath6 , written symbolically as @xmath7 where @xmath8 , @xmath9 , @xmath10 , and @xmath11 are , respectively , the newtonian , first , second and third pn contributions to the conservative dynamics of compact binaries , when the spin effects are neglected .
the leading order spin - orbit coupling to the binary dynamics is given by @xmath12 . in the above equation , @xmath13 , @xmath14 , and @xmath15 , where @xmath16 and @xmath17 are the relative separation vector and its conjugate momentum vector , respectively .
the familiar symbols @xmath3 and @xmath18 have the usual meaning , namely the total mass and the reduced mass .
the explicit expressions for pn corrections , associated with the motion of nonspinning compact binaries , were obtained in refs .
the spin - orbit contribution , available in refs .
@xcite , reads @xmath19 where the effective spin @xmath20 is defined by @xmath21 with    @xmath22    the finite mass ratio @xmath23 is given by @xmath24 .
the ( reduced ) spin vectors @xmath25 and @xmath26 are related to the individual spins @xmath27 and @xmath28 by @xmath29 and @xmath30 , respectively .
we recall that @xmath16 , @xmath17 , @xmath27 , and @xmath28 are canonical variables , such that the orbital variables commute with the spin variables , e.g. , see refs .
@xcite .    in ref .
@xcite , we obtained the keplerian - type parametrization in an orbital triad @xmath31 . in this triad ,
the unit vector @xmath32 is along the line of nodes , associated with the intersection of the orbital plane with the invariable plane @xmath33 , which is the plane perpendicular to the total ( reduced ) angular momentum @xmath34 , the only constant vector in our analysis .
the unit vector @xmath35 is defined by @xmath36 , where @xmath37 is the ( reduced ) orbital angular momentum vector , and hence @xmath35 is perpendicular to the orbital plane ( see fig .
[ fig:3fullplanes ] ) .
we display below , in its entirety , the parametric solution to the conservative 3pn accurate dynamics of spinning compact binaries moving in eccentric orbits , in adm - type coordinates , with spin effects restricted to the leading order spin - orbit interactions .
the dynamical vectors of the problem are given by @xmath38 we note that in the equal - mass case ( i ) the ( reduced ) total spin @xmath39 is the sum of two arbitrary spins @xmath25 and @xmath26 , and in the single - spin case ( ii ) @xmath39 is either @xmath25 _ or _ @xmath26 . in our restricted cases , namely ( i ) and ( ii ) , the ( reduced ) orbital angular momentum @xmath40 and the ( reduced ) total spin @xmath39 precess about the conserved ( reduced ) total angular momentum vector @xmath41 at the same rate @xcite .
this configuration is called _ simple precession _ in the literature , e.g. , see ref .
@xcite .
the time dependent basic vectors @xmath31 are explicitly given by    [ eq : ijk_in_the_combined_sol ] @xmath42    the angle @xmath43 denotes the constant precession cone angle of @xmath44 around @xmath34 and is given by    [ eq : theta_in_terms_of_slj ] @xmath45    where the magnitude of the total angular momentum is given by @xmath46 and @xmath47 is the constant angle between @xmath40 and @xmath39 .
the time evolution for @xmath48 , @xmath49 , and @xmath50 is given in a parametric and pn accurate form , which reads    [ eq : finalparamadm_3pn_plus_so ] @xmath51 \,.\end{aligned}\ ] ]    the pn accurate expressions for the orbital elements @xmath52 , @xmath53 , @xmath54 , @xmath55 , @xmath56 , and @xmath57 and the pn orbital functions @xmath58 , @xmath59 , @xmath60 , @xmath61 , @xmath62 , @xmath63 , @xmath64 , @xmath65 , @xmath66 , @xmath67 , @xmath68 , and @xmath69 , in terms of @xmath70 , @xmath71 , @xmath72 , @xmath23 , and @xmath47 , are available in ref .  @xcite .
the mass dependent coupling constant @xmath73 connects the effective spin with the total ( reduced ) spin via @xmath74 and is given by @xmath75 we display below the 3pn accurate expression for the _ internal _
parameter @xmath56 , appearing in eq .
[ this parameter was denoted by @xmath0 in ref .
@xcite , see eqs .
( 5.5c ) and ( 5.6k ) there ] .
it consists of two parts : the first one due to the 3pn accurate conservative dynamics of nonspinning point masses and the second one due to the leading order spin - orbit interaction @xmath76 where    [ eq : parameter_k_3pn_and_k_so ] @xmath77 \biggr\ } \ , , \\ \label{eq : just_parameter_k_so } \tilde{k}_\text{so } & = \frac{3}{c^2 l^2 } \biggl\ { - \frac { \chi_\text{so } } { 3 } - \chi_\text{so } \frac{s}{l } \cos\alpha \biggr\ } \,.\end{aligned}\ ] ]    it should be noted that @xmath78 gives the observable advance of periastron in the case of nonspinning compact binaries .
this results from the construction of @xmath56 [ see sec
.  v in ref .
@xcite for details and ref .
however , when spin - orbit interactions are added to the dynamics of compact binaries , @xmath56 does not directly lead to the measurable advance of periastron , relevant for binary pulsars .
this parameter @xmath56 is a measure for the advance of periastron _ only inside _ our parametrization and has to be related to the observational one .
further , note the presence of the @xmath73 term ( which does not depend on the spins ) , appearing in @xmath79 at the 1pn order .
it is due to these points , namely the @xmath73 term that appears at the 1pn order and @xmath56 does not give the observable periastron advance for spinning compact binaries , that we compute the dimensionless observable periastron advance parameter @xmath0 for spinning compact binaries in this paper .
note , that we have corrected the sign in front of @xmath80 in @xmath79 , important for the proper nonspinning limit and the following calculation of the observable periastron advance parameter @xmath0 in sec .
[ sec4 ] .
before we go onto the details of computing @xmath0 , let us list a few salient featurs of the parametrization .
the parameter @xmath81 , appearing in eq .  ,
is given by the newtonian contribution to say @xmath82 .
note that at newtonian order all three eccentricities @xmath82 , @xmath83 , and @xmath84 , are identical to each other and at pn orders they are related by pn accurate relations [ see eqs .
( 5.7 ) in ref .
therefore , it is possible to describe the binary dynamics in terms of one of the eccentricities .
we emphasize that , while adapting the parametrization to describe the binary dynamics , care should be taken to restrict orbital elements to the required pn order .
we are aware that the above parametrization does not lead to the classic keplerian parametrization simply by putting @xmath85 .
the detailed explanation of the nonspinning limit and the apparent lack of a _ simple _ nonspinning limit is given in sec .
iv  b in ref .
@xcite .
the above parametric solution is usually refered to as the genaralized quasi - keplerian parametrization for the pn accurate orbital dynamics of spinning compact binaries .
it extends the computations of damour , deruelle , schfer , and wex @xcite .
finally , we want to emphasize that the parametric solution , given by eqs .
, describes the entire conservative pn accurate dynamics of a spinning compact binary in an eccentric orbit , when the spin effects are restricted to the leading order spin - orbit interaction .
this means the parametrization consistently describes not only the precessional motion of the orbit inside the orbital plane , but also the precessional motions of the orbital plane and the spins themselves .
further , the deterministic nature of the dynamics described by the hamiltonian , given by eq .   and treated as a self - contained dynamical system , naturally follows from the above parametric solution .    in the next section ,
we develop the basics required to compute @xmath0 from our parametrization .
in this section , we will introduce two new triads , in terms of which it will be easy to obtain @xmath0 .
we will also discuss in detail how to connect these triads and hence the parametrization expressed in terms of these triads .
as we are interested in computing @xmath0 using our parametrization , it is natural that we introduce a periastron - based triad @xmath86 .
later , we will express the relative separation vector @xmath87 , given by eq .
, in terms of this new triad .
let us introduce this periastron - based triad in a pedagogical way .
we start by writing the angle @xmath88 , given by eq .
, in the following way @xmath89 remember , @xmath88 is the polar angle inside the orbital plane measured from the line of nodes @xmath32 , associated with the intersection of the orbital plane with the invariable plane @xmath33 , and @xmath2 is the `` true anomaly '' , i.e. , the angle between the ( instantaneous ) periastron and the position of the relative separation vector @xmath87 .
the angle @xmath90 may be considered as the `` argument of the periastron '' , i.e. , the angle between the line of nodes @xmath32 and the ( instantaneous ) periastron . by comparing eqs .   and
, we can easily infer the parametric evolution for @xmath90 .
therefore , it is obvious that the angle @xmath90 contains via @xmath56 the secular contributions to the periastron advance .    inserting eq .   into eq .
and making use of familiar trigonometric relations leads to @xmath91 where the new basic vectors are given by    @xmath92    now , it is obvious that @xmath93 is the unit vector pointing towards the ( instantaneous ) periastron , as @xmath87 is periodic in the angular variable @xmath2 with period @xmath94 in these new basic vectors @xmath93 and @xmath95 .
hence it is straightforward to define the periastron - based orthonormal triad @xmath96 , by rotating the @xmath31 frame around the orbital angular momentum unit vector @xmath35 , by an angle @xmath90 .
this gives us @xmath97 where @xmath90 is the angle between @xmath32 and the ( instantaneous ) periastron unit vector @xmath93 , as explained above .
the periastron - based vectors @xmath98 are co - rotating with the precessing ( instantaneous ) newtonian ellipse .
this is similar to the behaviour of a body fixed reference system usually encountered in classical mechanics @xcite .
hence , the set of angles @xmath99 is equivalent to the eulerian angles defined there .    in the next subsection
, we will introduce an observer - based triad and investigate how to connect the observer and periastron - based triads .
the second new triad that we introduce is a fixed one associated with the observer on earth and denoted by @xmath100 . in this fixed triad @xmath101 and @xmath102
are two orthogonal unit vectors in the plane of the sky , i.e. , the plane transverse to the radial direction @xmath103 , pointing from the source ( compact binary ) to the observer ( see fig .  [ fig:3fullplanes ] ) .    in order to compute @xmath0 using our parametrization ,
it is natural to expect that we need to connect the periastron - based triad @xmath86 to the observer triad @xmath100 .
it is possible to connect these two frames with the help of rotational matrices quite similar to the way in which the @xmath31 frame was connected to the @xmath100 frame in ref .
@xcite [ see sec .
vi in ref .
@xcite for details ] . in this paper
, we come from the observer triad @xmath100 to the periastron - based triad @xmath86 with the help of _ four _ successive rotations : first , a rotation around @xmath101 via the _ constant _ angle @xmath104 to come to the reference triad @xmath105 , second , a rotation around @xmath106 by an angle @xmath107 , third , a rotation around the direction of the line of nodes @xmath32 [ associated with the intersection of the orbital plane with the invariable plane @xmath33 ] by an angle @xmath43 , and , finally fourth , a rotation around the orbital angular momentum direction @xmath35 by an angle @xmath90 . in terms of rotational matrices  denoted by @xmath109 $ ]  we have @xmath110        \ ;   { \cal r}[\theta , { \boldsymbol{i}}(t ) ]        \ ;   { \cal r}[\upsilon(t ) , { \boldsymbol{e}}_{z } ]    \nonumber \\ & \quad \times { \cal r}[i_{0 } , { \boldsymbol{p } } ] \begin{pmatrix } { \boldsymbol{p } } \\ { \boldsymbol{q } } \\ { \boldsymbol{n } } \end{pmatrix } \,,\end{aligned}\ ] ] where the first rotation appears at the right - hand side , the second one appears left from the one before , and so on .
we note , that the approach to generate the @xmath86 frame via these rotations may be also adapted for scenarios , where the angle @xmath43 depends on time .
the constant behaviour of @xmath43 is a direct consequence of choosing two distinct cases ( i ) and ( ii ) , where our parametrization is valid and hence we will use the fact that @xmath43 is constant only at the end of our analysis .    it should be noted that there are other possible ways to connect the orthonormal triads @xmath100 and @xmath86 . a second way to generate the @xmath86 frame
is described in great detail by damour and schfer in sec .  5 of ref .
@xcite and involves _ three _ rotational matrices .
this way is closely linked to the observation as it involves rotational matrices defined with the help of the observational longitude of periastron @xmath111 and the orbital inclination angle @xmath112 .
damour and schfer connected the observer triad @xmath100 to the periastron - based triad @xmath86 via the following _ three _ successive rotations : first , a rotation around @xmath103 by an angle @xmath113 , second , a rotation around the direction of the line of nodes @xmath114 [ associated with the intersection of the orbital plane with the plane of the sky @xmath115 by an angle @xmath112 , and , finally third , a rotation around the orbital angular momentum direction @xmath35 by an angle @xmath111 .
the corresponding rotational matrices read @xmath116   \ ;   { \cal r}[i(t ) , { \boldsymbol{i}}_{i}(t ) ]    \ ; { \cal r}[\omega(t ) , { \boldsymbol{n } } ] \begin{pmatrix } { \boldsymbol{p } } \\
{ \boldsymbol{q } } \\ { \boldsymbol{n } } \end{pmatrix } \,.\end{aligned}\ ] ]    we note , that for the unit vector @xmath114 along the line of nodes , associated with the intersection of the orbital plane with the plane of the sky , the subscript @xmath112 is employed as @xmath114 is related to the orbital inclination angle @xmath112 and to avoid any confusion with the unit vector @xmath32 .    in the next section ,
we obtain @xmath0 with the help of the parametrization and the ideas discussed in this section .
in this section , we will compute @xmath0 , the observable dimensionless periastron advance parameter due to @xmath56 . following damour and schfer ,
the pn accurate total observable dimensionless periastron advance parameter @xmath0 can be written as @xmath117 where @xmath118 is the usual pn accurate periastron advance parameter up to 3pn order , corresponding to the nonspinning case [ see refs .
@xcite for the 2pn and 3pn contributions , respectively ] , and @xmath119 denotes the contribution due to spin - orbit interaction . it should be kept in mind that the 3pn contribution to @xmath118 will not be required while probing the dynamics of relativistic binary pulsars even with ska .
let us first briefly summarize how damour and schfer obtained @xmath0 , especially @xmath119 , in the next subsection .
timing of binary pulsars allows via the observed `` roemer time delay '' the direct measurement of the argument of periastron @xmath111 @xcite .
therefore , damour and schfer defined the relevant _ secular _ periastron advance as @xmath120 , or equivalently the timing parameter ( not to be confused with the orbital angular - momentum unit vector @xmath35 ) @xmath121 where @xmath54 is the mean orbital motion that appears in the pn accurate kepler equation , given by eq .  , and @xmath122 denotes the orbital average .
damour and schfer computed the contribution to @xmath123 due to spin - orbit interaction in the following way .
first , they computed the precessional equations for the orbital angular momentum vector @xmath124 and the laplace - runge - lenz vector @xmath125 , which tracks the periastron advance , with the help of the 2pn accurate version of eq .  .
these precessional equations where orbital averaged to obtain the rotational velocity vector , denoted by @xmath126 in ref .
the expression for @xmath123 was extracted from @xmath126 by decomposing @xmath126 along three _ non - orthogonal _ directions [ see eq .
( 5.11 ) in ref .
@xcite ] . in this way
a parametric solution to the dynamics of spinning compact binaries was not required to compute the 2pn accurate measurable dimensionless periastron advance parameter @xmath0 in ref .
@xcite , and hence the formula of damour and schfer for @xmath0 is more general .    in the next subsection
, we will compute @xmath123 using our keplerian - type parametric solution and the two different ways of connecting the @xmath100 frame to the @xmath86 frame .
this subsection deals with the computation of @xmath123 , arising due to spin - orbit interaction , using our keplerian - type parametrization .
following damour and schfer , we introduce the rotation velocity vector @xmath126 , containing all involved angular velocities , and extract @xmath123 from it in the following way .
we have already noted that there are two different ways to connect the observer - based triad to the pariastron - based one [ see eqs .   and and related discussions there ] .
the first way , given by eq .  , involves the angular variables @xmath127 ,
@xmath43 , @xmath50 , and @xmath104 , appearing in our parametrization , and the second way , given by eq .  , involves the angles @xmath128 , @xmath129 , and @xmath130 , where two of them are relevant for the observations , namely @xmath128 and @xmath129 .
therefore , we can decompose the rotation velocity vector @xmath126 in the following two ways :    [ eq : vec_omega_in_both_ways ] @xmath131    note , that the above two different decompositions are along four / three _ non - orthogonal _ directions , respectively .
the elementary rotation velocities , appearing in eq .  , can be defined with the help of our parametrization , whereas in eq .
we have @xmath123 : the rate of change of the argument of the periastron .
further note , that eq .
is identical to eq .
( 5.11 ) of ref .
@xcite .    in order to extract @xmath123 and to relate it to the angular variables associated with our parametrization we proceed as follows . since all vectors , appearing in eqs .
, can be expressed in terms of @xmath86 , it is clear that the rotation velocity vector @xmath126 may also be written as the sum of three elementary rotation velocities , along three _ orthogonal _ directions @xmath93 , @xmath95 , and @xmath35 , namely @xmath132 it is then straightforward to compute the components @xmath133 , @xmath134 , and @xmath135 in terms of the elementary rotation velocities and angles appearing in eqs .  . the components of @xmath126 in the co - rotating periastron - based @xmath86 frame  using the parametrization - based variables and the observation - based variables , respectively  are given by    [ eq : components_of_omega_in_abk ] @xmath136    these _ general _ set of equations provides the transformation from the parametrized angles to the angles of interest for the observation .
we note , that the last component , @xmath135 , gives the so - called _ intrinsic _ contribution to the periastron advance @xmath137 , defined by @xmath138 .
this is the component of @xmath126 acting _ within _ the orbital plane .
however , in the timing of binary pulsars the parameter @xmath0 enters in the ( relativistic ) argument of the periastron @xmath111 , and so we have to extract @xmath123 from the above set of equations .
a close inspection of fig .  [ fig:3fullplanes ] indicates that @xmath104 is constant , since it connects two invariant reference frames .
for the two distinct cases ( i ) and
( ii ) , where our parametrization is valid , we also have @xmath43 as a constant angle .
this allows us to come from eqs .
to the _ reduced _ set of equations for the components @xmath133 , @xmath134 , and @xmath135 , which reads    [ eq : reduced_components_of_omega_in_abk ] @xmath139    solving the last equation , namely eq .  , for @xmath123 leads to @xmath140 the first and the second contribution to @xmath123 are directly provided by our parametrization .
we also figured out that they can be obtained via the time evolution of the laplace - runge - lenz vector without any parametrization ( see appendix  [ appendix_lrl_vector ] for details ) .
the third contribution to @xmath123 in eq .
, namely @xmath141 , needs a closer inspection . to obtain @xmath142 ,
we multiply eq .   with @xmath143 and
add it to the product of eq .
times @xmath144 , which results in @xmath145 using now the explicit expressions for @xmath143 and @xmath144 , given by eqs .   in appendix
[ appendix_binary_geometry ] , in the above equation , we get @xmath146 where @xmath147 , @xmath148 , @xmath149 , and @xmath150 are shorthand notations for the constant trigonometric functions @xmath151 , @xmath152 , @xmath153 , and @xmath154 , respectively . the time evolution for the orbital inclination @xmath112 , appearing in eqs .
 , is simply given by eqs .   in appendix  [ appendix_binary_geometry ] .
let us now explicitly compute @xmath155 , @xmath156 and hence @xmath142 with the help of our parametric solution . using eqs .
, , and for @xmath90 , and eq .   for @xmath107 , we display below the relevant parametric expressions symbolically    [ eq : display_the_separation ] @xmath157    the structure of eqs .
indicates that we can handle the pn and the spin - orbit induced contributions to @xmath123 separately .    to obtain the spin - orbit contributions to @xmath158 , we compute with the help of eqs .   and the orbital averaged spin - orbit ( so ) contributions to the above required time derivatives    [ eq : so_contr_to_pi_upsilon_omega_reduced ] @xmath159    the quantity @xmath160 gives the orbital average of @xmath161 , @xmath162 being the scaled radial separation , to the newtonian order [ see , e.g. , ref .
@xcite ] @xmath163    with the help of eqs .  , , and , it is straightforward to obtain the spin - orbit induced contribution to @xmath0 , defined by eqs .   and , and it reads @xmath164 \,,\end{aligned}\ ] ] where @xmath165 is given in terms of our _ reduced _ variables by @xmath166 the proper nonspinning limit of eqs .   and , can be easily deduced as when @xmath167 , @xmath34 is identical to @xmath40 , which implies @xmath168 and @xmath169 , @xmath112 being the angle between @xmath103 and @xmath40 , and therefore @xmath119 vanishes . the disappearance of @xmath119 as @xmath167 indicates that the presence of the 1pn accurate @xmath73 term in @xmath79 , given by eq .  , does not have any observational consequence .
our result for @xmath119 is in agreement with eqs .
( 5.17 ) and ( 5.17a ) of ref .
@xcite , when we specialize these equations to the _ simple precession _
( for the equal - mass case and for the single - spin case ) @xmath170 \,,\end{aligned}\ ] ] where @xmath171 , @xmath172 and the coefficient @xmath165 is given here in terms of the _ non - reduced _ variables , used by damour and schfer in ref .
@xcite and in appendix  [ appendix_lrl_vector ] , by @xmath173 the dimensionless mass dependent coupling constant @xmath174 is related to @xmath73 by @xmath175 ( see appendix  [ appendix_lrl_vector ] for details ) . to check the above agreement between our result and the one of damour and schfer , we have used eqs .
and eq .  .
the two expressions for the coefficient @xmath165 , namely eqs .   and , are identical .
let us now briefly consider the nonspinning 3pn accurate contributions to @xmath0 . these contributions , enter only via @xmath90 , orbital averaged analogue to eq .
, our computation , and result in @xmath176 .
however , it should be noted that the expression for @xmath118 , as given by eq .  , is not yet ready for application .
this is so as this expression for @xmath118 depends on the conserved _ theoretical _ quantities @xmath70 and @xmath71 rather than on _ observable _ quantities like mean motion @xmath54 and proper - time eccentricity @xmath177 .
however , we have checked that up to 2pn order our expression for @xmath118 is in agreement with the one of damour and schfer , obtained in ref .
@xcite .
therefore , the 3pn accurate observable periastron advance parameter @xmath0 becomes @xmath178 where @xmath118 is the _
3pn contribution , corresponding to the nonspinning case , given by the 3pn accurate contribution to eq .
, namely eq .
. the spin - orbit contribution @xmath119 is given by eqs .   and or equivalently by eqs .   and
we recapitulate that @xmath176 and @xmath179 , where @xmath78 and @xmath79 are the corresponding contributions to @xmath56 , as given by eqs .   and .
this results from the construction of @xmath56 in ref .
@xcite by combining the parameter @xmath180 of the nonspinning case with the internal parameter @xmath79 , given by eq .
( 4.45e ) in ref .
@xcite , of the spin - orbit parametrization [ see sec .
v in ref .
@xcite for details ] .    finally , let us present the following geometric explanation why @xmath181 .
the difference between @xmath56 and @xmath0 results from the fact , that the angle @xmath90  associated with the parameter @xmath56  is measured from @xmath32 and in contrast to that , the angle @xmath111  associated with the observable parameter @xmath0  is measured with respect to @xmath114 .
the difference between @xmath111 and @xmath90 is denoted by @xmath182 in fig .
[ fig:3fullplanes ] ( see appendix  [ appendix_binary_geometry ] for details ) .
in addition , we want to point out , that the last rotation @xmath183 $ ] in eq
.   can be decomposed into two parts @xmath184 & = { \cal r}[\pi(t),{\boldsymbol{k}}(t ) ] \ ; { \cal r}[\delta(t),{\boldsymbol{k}}(t ) ]   \,,\end{aligned}\ ] ] and so it is obvious that @xmath185 and correspondingly we have @xmath186 .
the contribution of @xmath187 causes the difference between @xmath119 and @xmath79 .
we have presented a new method to compute the pn accurate observable dimensionless relativistic periastron advance parameter @xmath0 of a spinning compact binary using our keplerian - type parametric solution to the pn accurate conservative dynamics of spinning compact binaries moving in eccentric orbits @xcite .
the expression for @xmath0 is applicable , like our parametric solution , only for the following two distinct cases , namely , case ( i ) , the equal - mass case and case ( ii ) , the single - spin case .
our result is consistent with a more general formula for @xmath0 , derived by damour and schfer , using a different procedure @xcite .
the computation also helped us to clarify certain intriguing features that we observed in our parametrization .
it is our pleasure to thank gerhard schfer for helpful discussions and encouragements .
we gratefully acknowledge the financial support of the deutsche forschungsgemeinschaft ( dfg ) through sfb / tr7 `` gravitationswellenastronomie '' .    the algebraic computations , appearing in this paper ,
were performed using maple and mathematica .
in this section , we make contact with the variables and angles , commonly used in the literature @xcite .
it is obvious that we have three planes in the binary geometry , namely the orbital plane , the invariable plane and the plane of the sky .
the associated orthonormal triads are @xmath31 , @xmath105 , and @xmath100 .
therefore , we may define three lines of nodes and associated inclination angles in the following way .
we infer from fig .
[ fig:3fullplanes ] that the plane of the sky meets the invariable plane along a line of nodes , characterized by @xmath188 , as we choose the line - of - sight unit vector @xmath103 to lie in the @xmath189-@xmath190 plane .
the associated constant angle @xmath104 between @xmath103 and @xmath191 may be treated as the inclination angle , related to the above line of nodes . in the second case
, we consider the line of nodes @xmath32 , created by the intersection of the orbital plane with the invariable plane .
the corresponding inclination angle , namely the angle between the orbital angular momentum unit vector @xmath192 and @xmath191 , is denoted by @xmath43 . finally , we consider the intersection of the orbital plane with the plane of the sky .
we associate @xmath114 with this line of nodes and the related inclination angle @xmath112 , the angle between @xmath35 and the line of sight unit vector @xmath103 , is the usually defined orbital inclination @xmath112 and hence the use of the subscript @xmath112 in @xmath114 .
it is easy to figure out that the unit vectors , associated with the last two lines of nodes , @xmath32 and @xmath114 are defined by    @xmath193    using @xmath35 , expressed in terms of @xmath105 and @xmath100 , given by @xmath194 we compute the following expressions for @xmath32 and @xmath114 as    @xmath195 \,,\end{aligned}\ ] ]    where @xmath196 . the time evolution for the ( observational ) orbital inclination @xmath112 is simply given by    [ eq : time_evolution_of_i ] @xmath197^{1/2 }   \nonumber \\ & = \left [ 1 - ( c_{i_{0 } } c_\theta - s_{i_{0 } } s_\theta \cos\upsilon ) ^2 \right]^{1/2 }   \,.\end{aligned}\ ] ]    the equation for @xmath198 is in agreement with eq .
( 6 ) in ref .
we may also obtain a differential equation for the time evolution of @xmath112 , as given in ref .
@xcite . to get that , we differentiate @xmath199 with respect to time and make use of @xmath200 with @xmath201 , @xmath196 and eq .  .
in this way we deduce the following differential equation for the orbital inclination @xmath112 : @xmath202 which agrees with eq .
( 5.15 ) of ref .
@xcite , whereas we note that the quantity characterizing the _ ascending node _ in ref .
@xcite , namely @xmath32 , is identical to our @xmath114 .
following ref .
@xcite , let @xmath113 be the longitude of the line of nodes @xmath114 , as measured from @xmath101 in the @xmath101-@xmath102 plane ( plane of the sky ) , which implies @xmath203 the time evolution for @xmath113 is defined by either of the following two equations    @xmath204    finally , we note that it is also possible to obtain the time evolution for the longitude of periastron @xmath111 , measured from @xmath114 in the orbital plane , as defined in ref .
@xcite . to do that
, it is advantageous to define a periastron - based orthonormal triad @xmath205 .
this is easily obtained by a rotation of the @xmath31 frame around the orbital angular momentum unit vector @xmath35 , by an angle @xmath90 .
this gives us @xmath206 where @xmath90 is the angle between @xmath32 and the ( instantaneous ) periastron unit vector @xmath93 .
the time evolution of @xmath111 is given by either of the following relations    [ eq : cos_omega_sin_omega ] @xmath207 \ , , \\
\label{eq : sin_omega_in_appendix } \sin\omega & = |{\boldsymbol{i}}_{i } \times { \boldsymbol{a}}| = \frac{1}{\sin i } \left [ ( c_{i_{0 } } s_\theta + s_{i_{0 } } c_\theta   \cos\upsilon ) \sin\pi \right . \nonumber \\ & \quad\left.+ s_{i_{0 } } \sin\upsilon \cos\pi \right ] \,.\end{aligned}\ ] ]    the parametric solution to @xmath90 , which also gives its time evolution , is given by our parametrization via eq .  .
in addition , the parametric solution to @xmath90 can be also obtained by noting that @xmath208 , where @xmath125 is the laplace - runge - lenz vector given by @xmath209 \,.\end{aligned}\ ] ] further details are described in appendix  [ appendix_lrl_vector ] .
eq .   for @xmath143
is in agreement with eq .
( 6 ) in ref .
the additional check @xmath210 leads to @xmath211^{1/2 } \,,\end{aligned}\ ] ] which is equivalent to eq .
( using trigonometric relations ) .    since @xmath114 , @xmath32 , and @xmath93
lie in the orbital plane , we may define @xmath182 such that it is the angle between @xmath114 and @xmath32 .
this implies that @xmath212 ( see fig .  [ fig:3fullplanes ] ) and the time evolution for @xmath182 is given by    @xmath213    finally , let @xmath214 be the polar angle of @xmath87 , measured from @xmath114 , as defined by damour and schfer in ref .  @xcite .
as we define @xmath88 from our @xmath32 , this leads to @xmath215 . using these inputs , it is possible to write the relative separation vector @xmath87 in @xmath98 in terms of @xmath214 and @xmath111 as @xmath216 this is consistent with eqs .
( 5.5b ) , ( 5.6 ) , and ( 5.10 ) in ref .  @xcite . in our angular variables ,
.   reads @xmath217
in this section , we present an elegant way , starting from the famous laplace - runge - lenz vector , to obtain the time evolution of the apsidal motion without any parametrization .
the quantities in this section are non - reduced for a simple comparism with the work of damour and schfer  @xcite , and the work of wex and kopeikin  @xcite .
the first equation gives the orbital angular momentum vector @xmath124 and the second one the apsidal vector @xmath125 , which is the famous supplementary first integral of the kepler problem .
the vector @xmath125 is usually named after runge and lenz or laplace , though he was first discovered by lagrange .        where @xmath222 denotes the poisson brackets and @xmath223 is the full hamiltonian for the relative motion of our binary system . as @xmath124 and @xmath125 are first integrals of @xmath224 , only @xmath225 contribute to the right - hand sides of eqs .  .
the effects of @xmath226 are already studied in the literature , and so we just want to consider the contribution of @xmath227 where the effective spin @xmath228 is given by @xmath229 we recall that @xmath16 , @xmath17 , @xmath27 , and @xmath28 are canonical variables , such that the orbital variables commute with the spin variables , e.g. , see refs .
@xcite . using the basic poisson brackets @xmath230 , @xmath231 , @xmath232 , @xmath233 ,
we find      after orbital averaging these expressions , we can write down the equations for the precessional motion , induced by the spin - orbit interaction , in the vectorial form with the help of the rotational velocity vector @xmath235        in case of the so - called _ simple precession _ , where the orbital angular momentum @xmath124 and the total spin @xmath239 precess about the conserved total angular momentum vector @xmath240 at the same rate @xcite , the precession of the orbital plane and the longitude of periastron are best described by the angles @xmath107 and @xmath90 , defined with respect to the invariable plane , i.e. , the plane perpendicular to @xmath241 ( see fig .
[ fig:3fullplanes ] for details ) :      where @xmath243 denotes the unit vector in the direction of the orbital angular momentum vector @xmath124 .
the conditions for the _ simple precession _ are fulfilled in our restricted cases , namely ( i ) , the equal - mass case and ( ii ) , the single - spin case . in these cases the effective spin @xmath228 and the total spin @xmath239 are related by @xmath244 , where the mass dependent coupling constant is given by @xmath245 note , that @xmath175 .
using @xmath244 in eq .   and comparing eqs .   with an averaged version of eqs .
, we obtain the precession in terms of the angles @xmath107 and @xmath90              following urls provide a wealth of information about the ground and space based gravitational wave interferometers : http://www.ligo.caltech.edu , http://www.virgo.infn.it , http://www.geo600.uni-hannover.de , http://tamago.mtk.nao.ac.jp , and http://lisa.nasa.gov .