astrometric space missions , especially the esa ( european space agency ) cornerstone mission gaia , see e.g. _ perryman , et al . _
@xcite , are in preparation to attain microarcsecond ( @xmath2 ) level of accuracy in absolute positional measurements of stars and other celestial objects .
this unprecedented accuracy of astrometric observations makes it necessary to account for many subtle effects which were totally negligible before .
a practical model for astrometric observations with an accuracy of @xmath3 has been formulated by _
@xcite , where the influence of the gravitational fields inside the solar system were taken into account .
furthermore , a number of additional effects potentially observable at this level of accuracy due to various gravitational fields generated outside of the solar system were also briefly discussed in that investigation .
one of them is the gravitational light deflection of one companion of a binary system in the gravitational field of the other companion ( without loss of generality , throughout the investigation the light deflection of component b at component a is considered , hence , component a is considered to be the massive body , while component b is the light - source ) .
this effect would change the apparent position of one component of a binary system . while it is clear that this effect is relatively small and even at the level of @xmath3 observable only for edge - on binary systems ,
the huge amount of binaries generate some hope that there are relevant systems where this light deflection effect becomes detectable .
for instance , gaia will observe @xmath4 stars brighter than @xmath5 apparent magnitude .
detailed numerical simulations predict the detection of about @xmath6 ( resolved , astrometric , eclipsing , spectroscopic ) binary systems by gaia mission , see _ zwitter & munari _
@xcite , which is a considerable increase compared to the @xmath7 binary systems known so far , see `` _ washington double star catalog _ ''
@xcite .
let us consider this argument a bit more quantitatively .
for a simple estimate of the expected order of magnitude in light deflection , the classical lens equation ( in the form given by eq.(67 ) in _ fritelli , et al . _
@xcite , eq .
( 24 ) in _ bozza _ @xcite or eq .
( 23 ) in _ zschocke _
@xcite ) is applied . in terms of orbital elements of a binary system
it can be written as follows : @xmath8    here , @xmath9 is the light deflectin angle , @xmath10 is the inclination , @xmath11 is the semi - major axis and @xmath12 is the distance of center - of - mass of the binary system from the observer , and the schwarzschild radius is @xmath13 , where @xmath14 is the gravitational constant and @xmath15 is the speed of light , and @xmath16 is the stellar mass of component a of the binary system . from ( [ classical_lens_1 ] )
one obtains the maximal possible light deflection for edge - on binaries , i.e. attained for the case when the inclination is exactly @xmath17 : @xmath18    here , @xmath19 is the astronomical unit , @xmath20 stands for parallax of one arcsecond , and @xmath21 is the solar mass . according to this formula
, the choice of moderate values like @xmath22 and @xmath23 would result into significant light deflection effect on microarcsecond level even at large distances of about @xmath24 .
a meaningful value for the density of stars in the solar neighborhood , @xmath25 , implies already @xmath7 binaries inside a sphere of @xmath26 .
thus , one might conclude among them there are a few relevant edge - on binary systems .
however , in order to estimate quantitatively the number of relevant systems , one needs to know the probability for the occurrence of such edge - on binaries which have a given light deflection depending on their orbital parameters like inclination , mass , distance and semi - major axis .
such a relation between a given light deflection and orbital parameters is given by a so - called inclination formula .
since the inclinations of binary systems are of course randomly distributed , it is meaningful to resolve such an inclination formula in terms of inclination . as it has been shown by _
klioner , et al .
_ @xcite ( see appendix [ appendix_kms ] for some basic steps ) , such an inclination formula can be obtained by means of the analytical solution of light deflection in standard post - newtonian approach , and is given as follows : @xmath27    for a better illustration of the inclination formula , relation ( [ kms_d ] ) is rewritten in terms of angular degrees instead of radians : @xmath28    where also @xmath29 has been used . according to this relation , the inclination @xmath10 of a binary system with stellar mass @xmath16 and at distance @xmath12 must not deviate from the edge - on value @xmath17 too much in order to have a given light deflection @xmath9 .
for example , for a hypothetical binary star with @xmath22 situated at a distance of @xmath30 the light deflection effect attains @xmath3 only if @xmath31 , which means that the probability to observe this binary at a favorable inclination is only about @xmath32 .    but even this estimation is still much too optimistic . as a concrete example of todays high - precision astrometry ,
let us consider one important parameter about the astrometric accuracy of the gaia mission : the accuracy of one individual positional measurement in the most ideal case ( bright star , i.e. @xmath33 magnitude ) amounts to be @xmath34 , which implies @xmath35 .
furthermore , inside a sphere of @xmath36 around the sun almost every star and binary system is known already by the data of _ research consortium on nearby stars ( recons ) _ @xcite . since inside that sphere there is no binary system having a light deflection on microarcsecond level
, one has to take at least @xmath37 . by taking into account these both remarks
, one obtains @xmath38 in relation ( [ condition_0 ] ) .
therefore , even in the best case one has to conclude @xmath39 , that means ( for @xmath22 ) , the probability to observe such a binary at a favorable inclination is practically only about @xmath40 . and the binaries must be , in fact , almost edge - on in order to have a light deflection which can be detected by todays astrometry .
accordingly , while relation ( [ classical_lens_2 ] ) triggers the hope about the existence of many relevant binary systems , from relation ( [ condition_0 ] ) one concludes that the number of relevant binary systems is considerably reduced .
an estimation of the total amount of binaries depends on many different parameters , like mass , semi - major axis , inclination and distances of the binaries .
therefore , a simple estimation is not so straightforward as one might believe .
moreover , the relation ( [ classical_lens_2 ] ) has been obtained with the aid of classical lens equation , while relation ( [ condition_0 ] ) has been obtain by means of standard post - newtonian approach .
these both approaches have different regions of validity .
however , a rigorous treatment of the problem of light deflection in binary systems implies the need of an analytical formula which is valid for such kind of extreme astrometric configurations like the binary systems are .
recently , a generalized lens equation has been derived by _
zschocke _ @xcite , which allows to determine the light deflection of binary systems on microarcsecond level .
one aim of this study is , therefore , to reobtain the criteria ( [ classical_lens_2 ] ) and ( [ condition_0 ] ) as stringent conditions from one and the same approach , i.e. with the aid of generalized lens equation .
this is possible , because in the corresponding limits the generalized lens equation agrees with the classical lens equation and the standard post - newtonian solution .
another aim is , to derive an inclination formula like ( [ kms_d ] ) for binary systems which allows to determine the needed inclination for a given light deflection angle and as a function of the orbital parameters . for that one
has to take into account the distribution of stellar masses and the distribution of semi - major axes in binary systems .
finally , the aim of this study is , to determine the total number of relevant binaries having a light deflection on microarcsecond level , and to investigate the possibility to detect this effect of light deflection by todays high precision astrometry .
the article is organized as follows : in section [ section1 ] some basics about orbital elements of binary systems are given . the generalized lens equation and the inclination formula
are presented in section [ section2 ] . in section [ condition1 ] two stringent conditions on the orbital elements of binary systems ( astrometric , spectroscopic , eclipsing and resolved binaries )
are presented , which allow to determine whether or not the binary system will have a light deflection of a given magnitude .
the total number of binaries which have a given light deflection for an infinite time of observation is estimated in section [ section3a ] , while the more practical case of a finite time of observation is considered in section [ section3b ] .
the special case of resolved binaries is considered in section [ condition2 ] .
for that , the specific instrumentation of gaia mission is considered in some detail as the most modern astrometric mission with todays highest possible accuracy .
a summary is given in section [ summary ] .
consider a binary system , component @xmath11 with mass @xmath41 at coordinate @xmath42 and component @xmath43 with mass @xmath44 at coordinate @xmath45 . in order to express the light deflection effect in terms of orbital elements , spherical coordinates
are introduced , illustrated by fig .
[ fig : orbital_elements ] .    , semi - major axis @xmath11 , inclination @xmath10 , eccentricity @xmath46 , eccentric anomaly @xmath47 , periapsis @xmath48 , and mass ration @xmath49 .
the orbit of the binary system spans the @xmath50-plane and the @xmath51-axis is perpendicular to the orbital plane .
the @xmath52-axis is oriented along the semi - major axis of the orbit of the binary system , while the @xmath53-axis is perpendicular to the @xmath52-axis .
the vector @xmath54 is directed from the center - of - mass ( cms ) of the binary system , see eq .
( [ inclination_3 ] ) , to the observer .
the center of spherical coordinate system is located at the cms of binary system , i.e. @xmath55 . the inclination @xmath56 is the angle between @xmath54 and @xmath51-axis ; @xmath57 is called edge - on and @xmath58 corresponds to retrograde orbit .
the dotted line indicates the projection of @xmath54 onto orbital @xmath50-plane , i.e. z - component of @xmath54 equals zero .
the angle between this projection and @xmath52-axis is called argument of periapsis @xmath59 . the orbital elements
semi - major axis @xmath11 , eccentricity @xmath60 and mass ratio @xmath49 govern uniquely the geometric shape of both ellipses .
the eccentric anomaly @xmath61 ( not plotted here ) , is defined in eq .
( [ eccentric_anomaly_5 ] ) of appendix [ appendixb ] and determines the actual position of the bodies a and b on their orbit.,width=453 ]    the center of coordinate system is located at the cms , i.e. @xmath62    thus , the vector @xmath54 , which points from cms to the observer , is given by @xmath63{c } \displaystyle r\,\cos \omega\;\sin i \\
\nonumber\\ \displaystyle r\,\sin \omega\;\sin i \\
\nonumber\\ \displaystyle r\,\cos i \end{array}\right)\ , , \label{inclination_5}\ ] ]    where @xmath64 , the argument of periapsis is denoted by @xmath48 , and @xmath10 is the inclination , see fig .
[ fig : orbital_elements ] .
the solution of equation of motion yields for vectors @xmath42 and @xmath45 the expression given by eqs .
( [ appendixb_56 ] ) - ( [ appendixb_62 ] ) .
the vector @xmath65 points from the mass center of massive body to the observer , and vector @xmath66 points from the mass center of massive body to the source , see also fig .
[ fig : orbital_elements ] .
the coordinates of these vectors can be expressed by the orbital elements of the binary star as follows : @xmath67{c } \displaystyle r\,\cos \omega\;\sin i - \frac{a\ , \left(\cos e - { \rm e } \right)}{1 + \frac{\displaystyle m_a}{\displaystyle m_b}}\\ \nonumber\\ \displaystyle r\,\sin \omega\;\sin i - \frac{a\,\sqrt{1 - { \rm e}^2}\sin e}{1 + \frac{\displaystyle m_a}{\displaystyle m_b}}\\ \nonumber\\ \displaystyle r\,\cos i \end{array}\right)\ , , \label{inclination_10}\ ] ]    @xmath68{c } a \ ; \left(\cos e - { \rm e } \right )   \\
\nonumber\\ \displaystyle a\ ; \sqrt{1 - { \rm e}^2 } \sin e \\
\nonumber\\ \displaystyle 0 \end{array}\right)\ ,
. \label{inclination_15}\ ] ]    here , @xmath11 is the semi - major axis , @xmath46 is the eccentricity , and @xmath47 is the eccentric anomaly , see appendix [ appendixb ] .
the vectors ( [ inclination_10 ] ) and ( [ inclination_15 ] ) will be used to express the light deflection in terms of orbital elements of the binary system .
a scheme of light propagation of a signal emitted at component b in the gravitational field of component a is shown in fig .
[ fig : distance1 ] .
the vector @xmath65 points from the mass center of massive body to the observer , and vector @xmath66 points from the mass center of massive body to the source , and we define @xmath69 , the absolute value @xmath70 and unit vector by @xmath71 . furthermore , the impact vector @xmath72 is defined , and its absolute value is denoted by @xmath73 .
the schwarzschild radius of massive body , i.e. of component a of the binary system is denoted by @xmath74 .        for determining the light deflection in weak gravitational fields
there are two essential approaches : standard post - newtonian approach , e.g. _ brumberg _
@xcite , and classical lens equation e.g. eq.(67 ) in _ fritelli , et al . _
@xcite , eq .
( 24 ) in _ bozza _
@xcite , or eq .  ( 23 ) in _ zschocke _ @xcite . while the first approach is restricted by the condition @xmath75
, the second approach is only valid for the case that source and observer are far from the massive body , especially for @xmath76 and @xmath77 ( for geometrical illustration of @xmath78 , @xmath79 and @xmath80 see fig .
[ fig : distance1 ] ) .
however , in binary systems extreme configurations are possible like @xmath81 or @xmath82 .
therefore , in order to investigate the light deflection in binary systems one needs a generalized lens equation which is valid in such extreme configurations where the standard post - newtonian approach as well as the classical lens equation can not be applied .
recently , _ zschocke _ @xcite has derived a generalized lens equation , which allows to determine the light deflection in such extreme astrometric configurations like the binary systems are : @xmath83    actually , the lens equation has two solutions , but here only one solution is considered , while the second solution represents just the second image of one and the same source which is not relevant in our investigation .
the generalized lens equation is valid up to terms of the order @xmath84 , and the absolute value of their total sum can be shown to be smaller or equal to @xmath85 . here
, @xmath86 is chandrasekhar s impact parameter , see _ chandrasekhar _
@xcite , where @xmath87 being the orbital momentum and @xmath47 is the energy of the photon in the gravitational field of massive body .
basically , the light - ray of component b can not be observed if @xmath88 is smaller than the radius of massive body a. for stars , the radius is much larger than schwarzschild radius @xmath89 , hence @xmath90 .
furthermore , the generalized lens equation ( [ generalized_lens_1 ] ) is finite for @xmath91 and @xmath92 , both of which are possible astrometric configurations in binary systems .
furthermore , in _ zschocke _ @xcite it has been shown that the generalized lens equation ( [ generalized_lens_1 ] ) yields in the appropriate limits the correct standard post - newtonian solution and the classical lens equation , hence provides a bridge between these essential approaches .    ) or ( [ generalized_lens_3 ] ) , respectively .
the parameters chosen are : distance @xmath93 , semi - major axis @xmath94 inclination @xmath95 , mass @xmath96 , mass ratio @xmath97 , eccentricity @xmath98 , argument of periapsis @xmath99.,width=302 ]    in the following , eq .
( [ generalized_lens_1 ] ) is applied in order to determine the light deflection in binary systems . for that , the coordinates @xmath66 and @xmath65 are used in the form as given by eqs .
( [ inclination_10 ] ) and ( [ inclination_15 ] ) , respectively . a typical light - curve of a binary system , calculated by means of generalized lens equation ( [ generalized_lens_1 ] )
is shown in fig .
[ fig : lightcurve ] .
now , an inclination formula is derived from the generalized lens equation ( [ generalized_lens_1 ] ) .
the impact of eccentricity is neglected , circular orbits @xmath100 are considered , implying @xmath101 .
thus , the coordinates @xmath66 and @xmath65 are simplified to the expressions ( [ appendixc_5 ] ) and ( [ appendixc_10 ] ) given in appendix [ appendixc ] .
furthermore , the maximal value of light deflection of a binary system is of interest , i.e. the astrometric configuration @xmath102 is considered here .
then , by inserting these coordinates in the generalized lens equation ( [ generalized_lens_1 ] ) one obtains ( see eq .
( [ appendixc_24_c ] ) in appendix [ appendixc ] ) up to terms of the order @xmath103 : @xmath104    the minimal value @xmath105 and maximal value @xmath106 of light deflection for the astrometric position @xmath102 follow from eq .
( [ inclination_22 ] ) : @xmath107    where in ( [ inclination_24 ] ) terms of the order @xmath108 have been neglected .
the expression ( [ inclination_22 ] ) can be reconverted in terms of inclination ( see appendix [ appendixd ] ) : @xmath109    where @xmath110    the inclination formula ( [ inclination_30_a ] ) yields the upper limit for @xmath111 of a binary system in order to have a given value of light deflection @xmath9 . note ,
that the values of @xmath9 can not be chosen arbitrarily , but they are restricted by @xmath112 and @xmath113 given by eqs .
( [ inclination_24 ] ) and ( [ inclination_25 ] ) , respectively .    the inclination formula ( [ inclination_30_a ] ) can considerably be simplified . from ( [ inclination_35 ] ) and ( [ inclination_40 ] ) one obtains by series expansion @xmath114    by means of ( [ inclination_25 ] ) , the last term in ( [ comparison_15 ] ) can be estimated to be smaller than @xmath115 . here
, it should be underlined that @xmath116 even at large distances @xmath117 and small values for semi - major axis @xmath118 .
thus , one obtains @xmath119    up to terms of the order @xmath120 and @xmath121 ; here the relation @xmath122 for @xmath123 has been used .
it should be underlined , that the applicability of ( [ inclination_49 ] ) is restricted by the condition ( [ appendix_kms_e ] ) given in appendix [ appendix_kms ] and by @xmath124 . here , it should be noticed that @xmath125 , even in such an extreme case like @xmath93 , @xmath126 and @xmath127 one obtains a small number @xmath128 ; for an analytical proof use the exact expression for @xmath112 .
due to @xmath116 , the impact of semi - major axis is of lower order and can be neglected in the inclination formula . the inclination formula ( [ inclination_49 ] )
can be expressed in terms of dimensionless quantities as follows : @xmath129    note , that expression ( [ inclination_50 ] ) agrees with an inclination formula ( [ kms_d ] ) derived at the first time by _ klioner , et al . _ @xcite ; the arguments of their work are represented in appendix [ appendix_kms ] .
the simplified inclination formula ( [ inclination_50 ] ) is not only useful for straightforward estimations about the order of magnitude , but it s simple structure provides also an obvious comprehension about the interplay of the individual terms .
in this section , two stringent conditions on the orbital elements are highlighted for binaries having a given light deflection .
these strict conditions are valid for any binary system : astrometric , spectroscopic , eclipsing and resolved binaries .
the first stringent condition follows from the maximal light deflection angle ( [ inclination_25 ] ) , given by @xmath130    it represents a strict criterion for the maximal light deflection of a binary system with given schwarzschild radius of component a , given semi - major axis @xmath11 , and given distance @xmath12 between binary system and the observer .
the second stringent condition on the orbital elements follows from the inclination formula in the simplified form as given by eq .
( [ inclination_50 ] ) . for a better illustration
, this condition is given in terms of angular degrees instead of radians .
using @xmath29 one obtains @xmath131    according to this strict condition , the inclination @xmath10 of a binary system with mass @xmath16 ( recall @xmath16 is the stellar mass of component a ) and at distance @xmath12 must not exceed the given value in order to have a light deflection @xmath9 .
both these stringent conditions eqs .
( [ condition_1 ] ) and ( [ condition_2 ] ) were already stated in the introductory section by eqs .
( [ classical_lens_2 ] ) and ( [ condition_0 ] ) , respectively .
however , is should be underlined here , that both ( [ condition_1 ] ) and ( [ condition_2 ] ) were obtained with the aid of one and the same approach , namely the generalized lens equation , while ( [ classical_lens_2 ] ) and ( [ condition_0 ] ) were obtained by means of the classical lens equation and the post - newtonian solution , that means by two different approaches .
the observability of light deflection effect in binaries implies the realization of both these stringent conditions ( [ condition_1 ] ) and ( [ condition_2 ] ) simultaneously for a given binary system .
but even if a given binary system satisfies these both conditions , the observability of light deflection effect is not guaranteed , because the astrometric position @xmath102 has to be reached during the time of observation .
nonetheless , as soon as the orbital elements @xmath12 , @xmath11 , @xmath89 and @xmath10 of the binaries are known , these both stringent conditions ( [ condition_1 ] ) and ( [ condition_2 ] ) allow to find a possible candidate for being a relevant binary system for a given light deflection @xmath9 depending on the instrumentation of the observer . however , as it will be shown in section [ section3a ] and section [ section3b ] , the existence of such systems is highly improbable .
in the previous section , the conditions on orbital parameters for a binary system have been determined in order to have a given magnitude of light deflection @xmath9 . in this section , the total number of such relevant binaries is determined . in order to estimate the total number of binaries having a given light deflection @xmath9 , the following formula is applied : @xmath132    here , @xmath133 is the density of binaries , @xmath134 is the semi - major axis distribution of binary systems , @xmath135 is the mass distribution of stars where @xmath136 is the mass - ratio of massive body ( component a ) and solar mass . the probability distribution @xmath137 to find
a binary system with given inclination @xmath56 , is a function of distance @xmath12 , semi - major axis @xmath11 , mass ratio @xmath138 , and the given light deflection angle @xmath9 . according to ( [ inclination_30_a ] ) ,
the probability distribution @xmath137 is given by ( the inclination of binary systems is of course a random distribution ) @xmath139    where @xmath140 and @xmath141 are given by eqs .
( [ inclination_35 ] ) and ( [ inclination_40 ] ) .
for the minimal distance of a binary system from the sun one can safely assume @xmath142
. from ( [ inclination_25 ] ) follows that beyond a sphere with radius @xmath143 only a very few exceptional binary systems might have a light deflection of more than @xmath3 .
the sun is located in the _ orion arm _ which has about @xmath144 across and approximately @xmath145 in length . for the estimate according to ( [ numerical_results_5 ] ) ,
it is meaningful to assume that the stars are homogeneously distributed inside the _ orion arm_. for the uniform star density @xmath146 , a value which is in line with the data of _ recons _ @xcite .
furthermore , a common presumption is , that about @xmath147 percent of all stars are components of a binary or multiple system , see _ duquennoy & mayor _ @xcite and _ halbwachs , et al .
then one obtains for the density of binaries @xmath148    let us now consider the distribution of semi - major axis @xmath11 in binary systems .
statistical investigations show that the distribution of binary semi - major axis is flat in a logarithmic scale over the range of six orders of magnitude , that means assumed to be valid in the large range @xmath149 , see _ kouwenhoven & de grijs _ @xcite .
the lower limit @xmath150 is determined by the semi - major axis at which roche lobe overflow occurs , while the upper limit @xmath151 depends on how large the averaged star density is . the logarithmic distribution is known as `` _ pik s law _ '' ( 1924 ) named after its discoverer , and is given by @xmath152 , a law which has also been confirmed by recent investigations , see _ poveda , et al .
this distribution is a consequence of the process of star formation as well as of the dynamical history of the binary system , and one can take this law as a given fact for numerical studies . accordingly , see appendix [ probability_distribution ] : @xmath153    furthermore , for the mass distribution @xmath154 let us recall the initial mass function ( imf ) which is the probablity that a star is newly formed with a stellar mass @xmath16 and is frequently assumed to be a power law @xmath155 .
originally , the imf has been introduced by _ salpeter _
@xcite for solar neighborhood region , who assumed the value @xmath156 and a validity region for stars with masses between @xmath157 and @xmath158 . during the past
decades the imf has been refined by several investigations .
especially , the numerical values of slope parameter @xmath159 and regions of validity have been proposed in subsequent investigations , e.g. _ scalo _
@xcite , _ robin , et al .
_ ninkovic & trajkovska _ @xcite , _ ninkovic _
@xcite , _ kroupa _ @xcite ; for a review see _ kroupa _ @xcite .
moreover , the imf does not necessarily coincide with the real mass distribution of stars , because imf describes mass distribution of a star formation , while the solar neighborhood mainly consists of evolved stars of main sequence . here ,
for simplicity this distribution is used as a given fact with @xmath160 , and the proposed region of validity @xmath161 , and @xmath162 is assumed . according to imf , one finds for @xmath163 ( see appendix [ probability_distribution ] ) @xmath164    in order to motivate that distribution further , one can compare ( [ mass_distribution_7 ] ) with the _ recons _ data @xcite where one finds a fair agreement . using ( [ inclination_30_b ] ) - ( [ mass_distribution_7 ] ) ,
the results of the estimate ( [ numerical_results_5 ] ) are shown in fig .
[ fig : number_binaries_2 ] ; recall , when evaluating this integral one has to take into account the boundary conditions given by eqs .
( [ inclination_24 ] ) and ( [ inclination_25 ] ) .    according to fig .
[ fig : number_binaries_2 ] , in total there are about @xmath165 binaries having a light deflection of at least @xmath166 . here , the observation time is assumed to be infinite ,
so that it is guaranteed that the optimal configuration @xmath102 is reached .    ) , having parameters such that the light deflection of component b at component a is larger than a given value for @xmath9 .
note , that for one individual observation the astrometric accuracy of gaia is about @xmath34 and the end - of - mission accuracy is about @xmath167 in the most ideal case ( bright star ) .
the observation time is assumed to be infinite , so that it is guaranteed that the optimal configuration @xmath102 is reached.,width=302 ]
in eq .  ( [ numerical_results_5 ] ) the number of binaries with a given maximal possible light deflection @xmath9 has been determined , just by taking for eccentric anomaly the value @xmath102 , that means the ideal configuration where the light deflection takes its maximal value ( note , the eccentricity @xmath168 ) .
it is , however , obvious that during the most part of the orbital motion one will have @xmath169 and the light deflection will be much smaller than the maximal possible light deflection angle @xmath9 . on the other side , the orbital period @xmath170 of relevant binaries , given by eq .
( [ appendixb_55 ] ) , will easily exceed the time of observation ; for instance , the gaia mission time is about @xmath171 .
therefore , it will be not very probable , that the component b will be just at the relevant position near the value @xmath102 , where the light deflection becomes observable on microarcsecond level . in order to determine that number of observable relevant binaries
, one has to extend eq .
( [ numerical_results_5 ] ) as follows , @xmath172    here , @xmath173 is the probability for the binary system to be in the region @xmath47 , where the light deflection is larger than a given value for @xmath9 .    ) where the binary system reaches the optimal configuration @xmath102 during an assumed observation time of @xmath174 , and having orbital parameters such that the light deflection of component b at component a is larger than a given value for @xmath9 .
note , that for one individual observation the astrometric accuracy of gaia is about @xmath34 and the end - of - mission accuracy is about @xmath167 in the most ideal case ( bright star).,width=302 ]    in the very same way , as applied for the derivation of the inclination formula ( [ inclination_30_a ] ) , one can reconvert ( [ appendixc_24_d ] ) in terms of eccentric anomaly @xmath47 and one finds the eccentric anomaly formula : @xmath175    where @xmath140 and @xmath141 are given by eqs .  ( [ inclination_35 ] ) and ( [ inclination_40 ] ) . for a given value of light deflection @xmath9 , the formula ( [ observable_30 ] ) yields the value of eccentric anomaly @xmath47 of a binary system characterized by semi - major axis @xmath11 and schwarzschild radius ( or mass ) @xmath89 at a distance @xmath12 .
however , the values of @xmath9 can not be chosen arbitrarily , instead they are restricted by @xmath176 and @xmath177 given by ( of course , only astrometric positions with @xmath178 are taken into account , because for the area @xmath179 the light deflection is negligible ) : @xmath180    where in ( [ observable_20 ] ) terms of the order @xmath108 have been neglected .
these expressions resemble the corresponding expressions in eqs .
( [ inclination_24 ] ) and ( [ inclination_25 ] ) . according to eq .
( [ observable_30 ] ) , the region where the binary system has a light deflection larger or equal @xmath9 is given by @xmath181 .
one has also to take into account , that during gaia mission time @xmath182 the component b moves along the orbit and could move into the region @xmath181 .
therefore , the probability @xmath173 that the binary system is during gaia mission time at least ones inside the relevant astrometric position with the value @xmath47 in ( [ observable_30 ] ) , is given by @xmath183    where the operator @xmath184 is defined by @xmath185{c } \displaystyle x \quad { \rm if } \quad x < 1\,,\\ \displaystyle 1 \quad { \rm if } \quad x \ge 1\ , .
\end{array } \label{observable_37}\end{aligned}\ ] ]    the probability distribution ( [ observable_35 ] ) has to be implemented in eq .
( [ observable_5 ] ) in order to determine the number of binary systems having a given light deflection @xmath9 and to be observable during gaia mission time @xmath182 .
the results of eq .
( [ observable_5 ] ) are shown in fig .
[ fig : number_binaries_3 ] ; recall , when evaluating this integral one has to take into account the boundary conditions given by eqs .
( [ inclination_24 ] ) and ( [ inclination_25 ] ) .
accordingly , there are only a very few binaries @xmath186 having a light deflection of at least @xmath166 , if the astrometric position @xmath102 is reached during an assumed observation time of @xmath174 .
in this section , the special case of a resolved binary system is considered . in order to investigate the detectability of the light deflection effect in binaries , todays most modern astrometric mission , the esa cornerstone mission gaia , and its instrumentation which provides the highest possible astrometric accuracy at the moment ,
will be considered here as a concrete example .
the core of gaia optical instrumentation consists of two identical mirror telescopes , astro-1 and astro-2 , with a rectangular pupil whose dimensions are @xmath187 , @xmath188 , and @xmath189 is the effective focal length .
the intensity is given by , see _ lattanzi , et al .
_ @xcite and _ lindgren _ @xcite : @xmath190    where @xmath191 , @xmath192 , @xmath193 and @xmath194 the width and length of rectangular mirror , @xmath195 is the wavelength of incident light - ray , and @xmath196 are the angle of observation , i.e. the angle between the axis of the rectangular aperture and the line between aperture center and observation point
. the intensity of incident light - ray at @xmath197 is denoted by @xmath198 .
the function @xmath199 in eq .
( [ rectangular_5 ] ) is the ( by @xmath200 normalized ) point spread function ( psf ) for monochromatic incident light with wavelength @xmath195 for a rectangular aperture .
the optical spectrum of stars is @xmath201 . in fig .
[ fig : resolving_power1 ] the psf for an incident monochromatic light - ray with @xmath202 is represented for gaia telescopes .    ) .
the incident monochromatic light - ray has a wavelength of @xmath203 .
the parameters of the rectangular telescope are : @xmath204 , @xmath205.,width=453 ]    most of the light is concentrated in the central bright rectangular shaped pattern .
the length @xmath206 and width @xmath207 of this rectangle is determined by the first zero - roots of ( [ rectangular_5 ] ) at @xmath208 and @xmath209 , respectively . from
that it follows , that @xmath210 , and @xmath211 . furthermore ,
if the diffraction pattern is shown on a screen at a distance @xmath212 , then the length and width is given by , see _ hog , et al .
_ @xcite : @xmath213    where @xmath212 is the focal length of the optic , i.e. of the rectangular gaia mirror .
accordingly , the given numerical values @xmath187 , @xmath188 , and @xmath214 result in @xmath215 and @xmath216 for the length and width of `` _ airy rectangle _ '' of gaia optics . note , that the `` _ airy rectangle _ '' has the same order of magnitude than the pixel size ( @xmath217 ) of the @xmath218 ccds ( charge - coupled device ) sensors of astrometric field part of the focal plane . in order to separate two pointlike sources , the distance between their centers of the rectangle has to be larger than either @xmath219 or @xmath220 . since @xmath221
, in this study the better resolution value @xmath220 is used , which corresponds to a resolution angle of @xmath222    the resolving power is the minimal angular distance between two objects to get separable by gaia instrumentation . with the parameters given above one
obtains the resolving power @xmath223 of gaia optics : @xmath224    in what follows this parameter is of fundamental importance in order to determine the ability of gaia to determine the light deflection in binary systems .      in this section , the question is addressed which and how many binary systems can be separated by gaia instrumentation among all those relevant binaries found in the previous section ; see fig .  [
fig : number_binaries_3 ] . in average
, gaia will observe each object @xmath225 times , but will not constantly observe these objects during mission time .
however , for simplicity the scanning law of gaia is approximated by assuming a permanent observation of all objects during the whole mission time .
furthermore , visual binaries are considered in this section , i.e. binaries which are separable by gaia telescopes .
the both largest telescopes of gaia have a resolution angle @xmath223 discussed in the previous section , see eqs .
( [ rectangula_30 ] ) and ( [ minimal_alpha_10 ] ) . for binary systems ,
this resolution angle @xmath223 corresponds to a minimal distance between the components a and b to get separable within gaia optics . using ( [ expression_d ] ) and ( [ rectangula_30 ] ) one obtains the condition @xmath226    where @xmath12 is the distance between the remote objects and gaia observer .
this condition is by far much more important than taking into account the effect of finite radius of the stars , which would imply @xmath227 , where @xmath228 being the radius of component @xmath11 . by inserting the extreme case @xmath229 into ( [ inclination_22 ] ) , one obtains @xmath230    where @xmath231 has been used , and terms of higher order @xmath232 are neglected .
relation ( [ visual_5 ] ) is an expression for the maximal light deflection angle of a binary system when taking into account the resolving power of gaia .
( [ visual_5 ] ) is a much stricter restriction than the generalized lens equation ( [ inclination_22 ] ) , because ( [ visual_5 ] ) determines the light deflection angle only of those binary systems having a resolution angle @xmath223 of gaia optical instrumentation , while ( [ inclination_22 ] ) determines the light deflection angle of any possible binary system .
note , that from eq .
( [ visual_5 ] ) follows the maximal possible distance of visual binaries for a given light deflection : @xmath233    where in the last expression the optimal values @xmath234 and @xmath235 has been used .
the condition ( [ condition_3 ] ) can also be written by @xmath236    these both conditions ( [ condition_3 ] ) and ( [ condition_4 ] ) imply rather extreme orbital parameters on visual binaries .
for instance , condition ( [ condition_3 ] ) implies a maximal distance of @xmath237 for solar - mass type binaries even with a huge semi - major axis of @xmath238 , while condition ( [ condition_4 ] ) implies a large semi - major axis for solar - mass type binaries outside a sphere of @xmath239 .
it is almost certain , that such extreme parameters will not be realized in nature .
in this study , the light deflection in binary systems has been considered . while there is absolutely not any doubt about the existence of this relativistic effect ,
is has not been observed so far . to investigate this effect of light deflection ,
an inclination formula ( [ inclination_30_a ] ) has been derived by means of generalized lens equation ( [ generalized_lens_1 ] ) obtained recently by _ zschocke _
@xcite , and these both equations are the theoretical basis for investigating the light deflection effect in binary systems . a simplified inclination formula has been presented by eq .
( [ inclination_50 ] ) and its validity has been discussed in some detail .
this simplified inclination formula has also been obtained by _ klioner , et al . _
@xcite by an independent approach .
furthermore , two stringent conditions on the orbital parameters have been given by eqs .
( [ condition_1 ] ) and ( [ condition_2 ] ) .
these both stringent conditions are relations between the orbital elements of a ( resolved , astrometric , eclipsing , spectroscopic ) binary system for a given magnitude of light deflection , and allow to find a relevant binary system in a straightforward way .    in section [ section3a ] , the total number of binaries with a given light deflection has been determined by means of the semi - major axis distribution according to `` _ pik s law _ '' and the mass distribution according to `` _
salpeter s mass distribution _ '' .
since the inclinations are randomly distributed , the inclination formula allows to estimate the total number of relevant binaries with the aid of eq .
( [ numerical_results_5 ] ) .
it turns out , that in total there exist about @xmath0 binaries having orbital parameters such that the light deflection amounts to be at least @xmath3 ; see fig .
[ fig : number_binaries_2 ] .    in section [ section3b ] a finite time of observation of @xmath240 years ( gaia mission time ) has been considered , which considerably reduces the total number of relevant binaries . clearly , this case is of more practical importance , since a restricted time window of observation is in more agreement with reality than the first scenario . by taking into account the probability to find the system in the ideal astrometric position @xmath102 where the light deflection becomes maximal , it has been found by evaluating the corresponding integral ( [ observable_5 ] ) that there is not any of the relevant binary systems in the ideal position @xmath102 during gaia mission time ;
see fig .
[ fig : number_binaries_3 ] .
thus , while in principle a few binaries will have a significant light deflection , the effect could not be detected due to the restricted finite time window of observation .
furthermore , the special case of resolved binaries has been considered in section [ condition2 ] .
the astrometric instrumentation of the esa cornerstone mission gaia , see e.g. _ perryman , et al . _
@xcite , has been considered in some detail in order to decide whether or not this subtle effect of light deflection can be observed .
two conditions for resolved binaries were presented in eqs .
( [ condition_3 ] ) and ( [ condition_4 ] ) for such special kind of binary systems .
it has been shown , however , that even for the gaia mission , which is an outstanding milestone of progress in astrometry , such binary systems must have rather extreme orbital parameters in order to reach todays level of detectability , i.e. on microarcsecond level .
the existence of such exotic binaries is , however , highly improbable .    in summary ,
the main results are presented by the inclination formulae ( [ inclination_30_a ] ) and its simplified version ( [ inclination_50 ] ) , the stringent conditions ( [ condition_1 ] ) and ( [ condition_2 ] ) and by the diagrams fig .  [
fig : number_binaries_2 ] and fig .
[ fig : number_binaries_3 ] .
accordingly , one comes to the conclusion that the detectability of light deflection in binary systems reaches the technical limit of todays high precision astrometry and might be detected only in case of a very few and highly exotic binary systems .
it is , however , very unlikely that such extreme binaries might exist .
it seems that the detection of the light deflection effect in binary systems needs an astrometric accuracy of better than about @xmath241 .
thus , only astrometric missions of the next generation can accept the challenge to detect this subtle effect of relativity .
this work was partially supported by the bmwi grants 50qg0601 and 50qg0901 awarded by the deutsche zentrum fr luft- und raumfahrt e.v .
( dlr ) . enlighting discussions with prof .
sergei a. klioner , prof .
michael h. soffel and prof .
francois mignard are gratefully acknowledged .
99    bozza , v. , 2008 , phys . rev .
d , 78 , 103005    brumberg , v. a. 1991 , _ essential relativistic celestial mechanics _ , ( bristol : adam hilder )    chandrasekhar , s. , 1983 , _ the mathematical theory of black holes _
( oxford : clarendon )    duquennoy , a. & mayor , m. , 1991 , a&a , 248 , 485    fritelli , s. , kling , t.p . &
newmann , e.t . , 2000 , phys .
d , 64 , 064021    halbwachs , j.l . ,
mayer , m. , udry , s. & arenou , f. , 2003 , a&a , 397 , 159    hog , e. , bastian , u. & seifert , w. , 1997 , _ optical design of gaia _ , www.en.scientificcommons.org    klioner , s.a . , 2003 , aj 125 , 1580    klioner , s.a . ,
mignard , f. & soffel , m. , 2003 , _
astrometric signature of gravitational microlensing on the components of edge - on binary systems _ , unpublished    klioner , s.a .
& zschocke , s. , 2010 class .
quantum grav . , 27 , 075015    kouwenhoven , m.b.n . & de grijs , r. , 2008 , a&a , 480 , 103    kroupa , p. , 2001 ,
mnras , 322 , 231    kroupa , p. , 2002 , science , 295 , 82 .
landau , l.d . & lifshitz , e.m .
, 1976 , _ mechanics _ , vol . 1 , ( oxford pergamon press )    lattanzi , m. , carollo , d. & gai , m. , 1998 , _ aberrated point spread function s for the mms configuration and astrometric error _ , technical report , sag - ml-014 ; available at the gaia document archive http://www.rssd.esa.int/llink/livelink
.    lindgren , l. , 1998 , _ point spread functions for gaia including aberrations _ , technical report from lund observatory , sag - ll-025 ; available at the gaia document archive http://www.rssd.esa.int/llink/livelink    ninkovic , s. , 1995 , bull .
belgrade , 151 , 1    ninkovic , s. & trajkovska , v. , 2006 , serb .
aj , 172 , 17    pik , e.j . , 1924 , tartu obs .
publ . , 25    perryman , m.a.c .
, de boer , k. , gilmore , g. et al . , a&a , 2001 , 369 , 339    poveda , a. , allen , c. & hernandez - alcantara , a. , 2007 , in proc .
240 ; arxiv : astro - ph/0705.2021    recons , 2012 , http://www.recons.org/    robin , a.c . ,
reyle , c. , derriere s. & picaud , s. , 2003 , a&a , 409 , 523    salpeter , e.e . , 1955 , aj , 121 , 161    scalo , j.m . , 1986 , fundam
, 11 , 1    wds , 2011 , http://ad.usno.navy.mil/wds/wds.html    zschocke , s. 2011 , _ a generalized lens equation for light deflection in weak gravitational fields _ , class .
quantum grav .
28 , 125016    zwitter , t. & munari , u. , 2004 , _ the environment and evolution of double and multiple stars _ , revmexaa ( sc ) , 21 , 25
the calculations in this appendix follow mainly _ landau & lifshitz _ @xcite . consider two massive bodies , one component having a mass @xmath41 and spatial coordinate @xmath42 , and second component with a mass @xmath44 and spatial coordinate @xmath45 , respectively .
they orbit around their common center of mass @xmath242 , @xmath243    the lagrangian @xmath244 of the two - body problem is given by @xmath245    where @xmath246 is the potential . with
the aid of relative coordinate @xmath247 and reduced mass @xmath248 , the two - body problem can be transformed into into one - body problem , @xmath249    polar coordinates @xmath250 yield @xmath251    the orbital angular momentum @xmath87 is conserved @xmath252    by means of which one obtains for the total energy of the two - body system the expression @xmath253    from eq .
( [ appendixb_25 ] ) one deduces @xmath254 - \frac{l^2}{\overline{m}^2\;r_{ab}^2}\right)\ ; , \label{appendixb_30}\end{aligned}\ ] ]    and from eq .
( [ appendixb_30 ] ) one obtains up to an integration constant @xmath255 - \frac{l^2}{\overline{m}^2\;r^2_{ab}}\right)^{-1/2}\ , , \label{appendixb_35 } \\
\nonumber\\ \phi & = & \int d r_{ab}\ ; \frac{\overline{m}}{r^2_{ab}}\ ; \left ( 2 \,\overline{m } \left [ e - u \left(r_{ab}\right ) \right ] - \frac{l^2}{r_{ab}^2}\right)^{-1/2}\ , , \label{appendixb_36}\end{aligned}\ ] ]    where in the second relation eq .  ( [ appendixb_20 ] ) has been used ; note that ( [ appendixb_36 ] ) is the relation between @xmath256 and @xmath257 and is called orbital equation .
the eqs .
( [ appendixb_35 ] ) and ( [ appendixb_36 ] ) are the general integral solutions of a two - body problem . in order to integrate these eqs .
( [ appendixb_35 ] ) and ( [ appendixb_36 ] ) one has to specify the potential @xmath246 . in case of kepler problem
one has @xmath258
( [ appendixb_36 ] ) can be integrated and yields @xmath259    where the axes are chosen such that the mentioned integration constant vanishes .
furthermore , by introducing the eccentricity @xmath260 ( possible values of eccentricity are between @xmath261 ; @xmath262 corresponds to a circular orbit ) , @xmath263    the solution ( [ appendixb_41 ] ) can be written as @xmath264    note the expressions of semi - major axis @xmath11 and semi - minor axis @xmath43 , @xmath265    to solve the integral ( [ appendixb_35 ] ) , one substitutes @xmath266    where @xmath47 is called eccentric anomaly .
then , one obtains for the integral in eq .
( [ appendixb_35 ] ) the expression @xmath267    and the solution is given by @xmath268    where the integration constant vanishes , i.e. the particle at @xmath269 is in periastron .
the eqs .
( [ appendixb_46 ] ) and ( [ appendixb_50 ] ) are the general solutions of two - body problem .
they can be rewritten as @xmath270    in case of ellipse , @xmath271 in periastron , @xmath272 in apastron , and for a complete orbit @xmath47 runs from @xmath271 to @xmath273 .
thus , one obtains for the orbital period the expression @xmath274    note , the solution @xmath54 in cartesian coordinates , @xmath275 and @xmath276 : @xmath277{l } x \\ \displaystyle y \end{array}\right ) , \label{appendixb_56 } \\ \nonumber\\ x & = & a \left ( \cos e - { \rm e } \right ) \ , , \label{appendixb_57 } \\
\nonumber\\ y & = & a \left ( 1 - { \rm e}^2 \right)^{1/2 } \;\sin e \ , .
\label{appendixb_58}\end{aligned}\ ] ]    the coordinates of the bodies a and b , i.e. their orbits , are given by @xmath278    accordingly , the geometry of the orbit is determined by two orbital parameters : semi - major axis a and eccentricity e. in order to know the position of one celestial body , either component a or component b , two additional orbital parameters are needed , namely orbital period @xmath170 and true anomaly @xmath279 . a geometrical representation of the coordinates of the components of a binary star is given in fig .  [ fig : ellipse ] for the case of @xmath280 .     and
@xmath281 , respectively .
the semi - major axis of the binary system is chosen @xmath282 and eccentricity is taken @xmath283 .
the coordinates of mass center are @xmath55 .
the massive bodies a and b are always in opposition to each other.,width=226 ]
for the inclination formula the impact of eccentricity on light deflection is neglected , thus @xmath284 , implying that @xmath285 is taken .
then , for the vectors from massive body to observer @xmath65 and from massive body to source @xmath66 , one has @xmath286{c } \displaystyle \sin i - \epsilon_1\,\cos e\\ \nonumber\\ \displaystyle - \epsilon_1\,\sin e
\\ \nonumber\\ \displaystyle \cos i \end{array}\right ) , \label{appendixc_5}\ ] ]    @xmath287{c } \displaystyle \cos e\\ \nonumber\\ \displaystyle \sin e \\
\nonumber\\ \displaystyle 0 \end{array}\right ) , \label{appendixc_10}\ ] ]    where the small parameter @xmath288    has been introduced . from eqs .
( [ appendixc_5 ] ) and ( [ appendixc_10 ] ) one obtains for vector @xmath71 , where @xmath69 , the expression @xmath289{c } \displaystyle \sin i + \epsilon_2\,\cos e\\ \nonumber\\ \displaystyle \epsilon_2\,\sin e \\
\nonumber\\ \displaystyle \cos i \end{array}\right ) , \label{appendixc_15}\ ] ]    where the small parameter @xmath290    has been introduced .
note that ( [ appendixc_10 ] ) and ( [ appendixc_15 ] ) yields @xmath291    using ( [ appendixc_5 ] ) - ( [ epsilon_2 ] ) , the generalized lens equation ( [ generalized_lens_1 ] ) reads @xmath292    where @xmath293 and @xmath294    by series expansion one obtains up to terms of order @xmath295 : @xmath296    where @xmath297 with @xmath298 . the minimal and maximal light deflection angle
are @xmath299    in this study , the maximal possible light deflection effect is of interest .
accordingly , two configurations are relevant : @xmath300    up to order @xmath301 which is just eq .
( [ inclination_22 ] ) , and @xmath302    up to order @xmath295 .
furthermore , it is useful to take into account only astrometric positions with @xmath178 , because otherwise the light deflection is for sure negligible .
from ( [ appendixc_24_c ] ) one obtains @xmath303    from ( [ appendixd_10 ] ) one obtains @xmath304    eq .
( [ appendixd_15 ] ) represents an quadratic equation for the expression @xmath305 , which has the following both solutions for the inclination @xmath10 : @xmath306    where @xmath307    eq .  ( [ appendixd_21 ] ) represents two solutions , however only the one with the plus - sign is valid .
this can be shown as follows . for the value @xmath308
the light deflection has to be @xmath309 , according to ( [ inclination_25 ] ) . by inserting @xmath113 in eqs .
( [ appendixd_30 ] ) and ( [ appendixd_35 ] ) one obtains @xmath310 and @xmath311 .
if one inserts @xmath308 for @xmath140 and @xmath141 into eq .
( [ appendixd_21 ] ) one obtains the relation @xmath312    obviously , relation ( [ appendixd_40 ] ) is only correct for the upper sign . a very similar proof can also be done using @xmath112 which also yields that the upper sign is the correct and only solution .
therefore , the inclination formula is given by ( note , that in the region under consideration @xmath313 ) @xmath314{c } \displaystyle \arcsin \left ( - \frac{p}{2 } + \sqrt{\frac{p^2}{4 } - q } \right )   \quad { \rm for } \quad 0 \le i \le \frac{\pi}{2 } \ , , \\ \\
\displaystyle \pi - \arcsin \left ( - \frac{p}{2 } + \sqrt{\frac{p^2}{4 } - q } \right ) \quad { \rm for } \quad \frac{\pi}{2 }
< i \le \pi\ , . \end{array } \label{appendixd_45}\end{aligned}\ ] ]    for the complete region @xmath56 one obtains for the inclination formula the following expression : @xmath315    where @xmath140 and @xmath141 are given by eqs .
( [ appendixd_30 ] ) and ( [ appendixd_35 ] ) .
in this appendix , some basic steps of the calculations of _ klioner , et al . _ @xcite are presented .
a scheme of the light propagation in a binary system is presented by fig .
[ fig : distance1 ] . the light signal from component b , considered to be the light source ,
is deflected by component a , considered to be the massive body which curves the space - time .
the vector @xmath65 points from the mass center of massive body to the observer , and vector @xmath66 points from the mass center of massive body to the source ; @xmath69 and its absolute value @xmath70 , and @xmath74 is the schwarzschild radius of massive body ; the explicit label a is omitted .
furthermore we define the impact vector @xmath72 is defined , and its absolute value @xmath73 , see also fig .
[ fig : distance1 ] .    according to _ klioner & zschocke _
@xcite , the transformation of @xmath316 to the unit tangent vector @xmath317 of light - trajectory at observer is in standard post - newtonian approach given by @xmath318    note , the ppn parameter of parameterized post - newtonian formalism , which characterizes a possible deviation of the physical reality from general theory of gravity , is set equal to @xmath1 here for simplicity .
this expression is valid as long as @xmath124 , but diverges for @xmath91 .
thus , it is not valid for all possible binary configurations , instead one has take care to consider only such astrometric configurations with @xmath124 . by means of ( [ appendix_kms_a ] )
one obtains for the light deflection angle @xmath9 , i.e. for the angle between @xmath317 and @xmath316 , the expression @xmath319    where @xmath320 , @xmath321 has been used , and @xmath322 is the angle between @xmath66 and @xmath65 . the expression ( [ appendix_kms_b ] ) diverges for @xmath323 , which corresponds with the mentioned divergence of ( [ appendix_kms_a ] ) for @xmath91 .
obviously , @xmath324 ( from eq .  ( [ observable_30 ] ) it is clear that eccentric anomaly @xmath47 of binary system should be very close to zero for the light deflection effect to be observable at the level of microarcsecond , i.e. one actually could even assume @xmath325 ) and one obtains @xmath326    where @xmath327 has been used , @xmath328 , and the asymmetry of function @xmath329 . from ( [ appendix_kms_c ] )
one obtains @xmath330    where @xmath331 has been used ( recall @xmath16 is the mass of component a and @xmath21 is the solar mass ) , and the numerical values @xmath332 , @xmath333 and @xmath334 have been inserted , so that @xmath335 .
the simplified inclination formula ( [ appendix_kms_d ] ) has first been obtained by _ klioner , et al .
_ @xcite ; here the fact is noticed , that due to the divergence of the post - newtonian solution ( [ appendix_kms_a ] ) for @xmath336 , which corresponds to @xmath323 , the applicability of ( [ appendix_kms_d ] ) is restricted by the condition @xmath124 . using @xmath337
one obtains the validity condition for the applicability of ( [ appendix_kms_d ] ) : @xmath338    equation ( [ appendix_kms_d ] ) provides a relation between inclination @xmath10 and light deflection @xmath9 for a binary system characerized by it s distance @xmath12 from the observer and it s stellar mass @xmath16 ; note that ( [ appendix_kms_d ] ) agrees with the simplified inclination formula given by equation ( [ inclination_50 ] ) .
assume a probability distribution of a quantity @xmath52 is given by @xmath339 .
the probability @xmath340 , to find a system in the interval @xmath341 is given by @xmath342    where the region of validity of probability distribution @xmath339 is given by @xmath343 and @xmath344 .
in the infinitesimal limit @xmath345 , one obtains by series expansion the following explicit form for the here used probability distributions : for a power law @xmath346 with @xmath163 one finds @xmath347