einstein s theory of general relativity ( gr ) predicts that space - time is curved and warped by the mass and the angular momentum of a gravitating body .
the kinematics of objects in the vicinity of such a body are expected to differ from those based on newton s theory of gravity .
the nasa / stanford spaceborne mission gravity probe  b ( ) was designed to measure these differences with four gyroscopes essentially freely falling around earth , in a polar orbit at an altitude of about 640 km . according to gr
, the gyroscopes spin axes should precess by @xmath46  yr@xmath5 in a north - south direction due to the curving of space - time by the earth s mass , and by @xmath439  in an eastward direction due to the warping of space - time by the earth s angular momentum .
these precessions are dubbed the
geodetic " effect and the  frame - dragging , " ( `` lense - thirring , '' or `` gravitomagnetic '' ) effect , respectively , and can be considered as rotations of any near - earth inertial frame relative to the distant universe .    the gyroscopes were set spinning within a quartz block with respect to which , through superconducting quantum interference devices ( squids ) ,
each of the spin axes precessions could be measured with an accuracy of 100  within @xmath440 days .
the quartz block was bonded to a @xmath415  cm diameter quartz optical telescope .
all of these devices  the quartz block , the squids , and the quartz telescope ( together termed the probe )  were placed in a dewar filled with liquid helium and cooled to a temperature of 1.8 k. ideally , the orientation of the probe could be best provided by pointing the telescope to one suitable quasar , one of the most distant compact types of single object in the universe .
however , quasars were too dim for the telescope to track .
only stars in our galaxy were bright enough .
the challenge was therefore to find a
guide star , " bright enough for the on - board telescope to track , sufficiently isolated to limit light contamination , and suitably located to not unnecessarily decrease the accuracy of the measurement of the precessions of the gyroscopes .
then the guide star s motion on the sky needed to be measured with respect to quasars or distant galaxies so that the precession of the gyroscopes could be determined with respect to the distant universe .
we described the astronomical effort undertaken with the guide star for the support of the  mission ,
vlbi for gravity probe b , " in seven papers :  overview " ( * ? ? ?
* paper i ) ;  monitoring of the structure of the reference sources 3c 454.3 , b2250 + 194 , and b2252 + 172 " ( * ? ? ?
* paper ii ) ;  a limit on the proper motion of the ` core ' of the quasar 3c 454.3 " ( * ? ? ? * paper iii ) ;  a new astrometric analysis technique and a comparison with results from other techniques " ( * ? ? ?
* paper iv ) ;  proper motion and parallax of the guide star , i m pegasi " ( * ? ? ? * paper v ) ;  the orbit of i m pegasi and the location of the source of radio emission " ( * ? ? ? * paper vi ) ; and  the evolution of the radio structure of i m pegasi " ( * ? ? ? * paper vii ) . here
we give a synopsis of these papers .
a suitable guide star for  had to simultaneously meet a number of requirements :    * the motion of the star needed to be known or measurable with an uncertainty small enough not to significantly increase the prelaunch anticipated standard error of  of @xmath60.5  in each of the two precessions .
the specific requirement was that the standard error be @xmath70.14  for measurements of the star s motion in each coordinate . in 1990 , when we entered the first stage of our search for a suitable guide star , no star s ( or any other celestial object s ) motion on the sky
was known with such accuracy in the optical .
such accuracy had previously been achieved only in the radio .
pioneered by @xcite , astrometric measurements with the radio technique of very long baseline interferometry ( vlbi ) of a quasar relative to another quasar nearby on the sky yielded a proper - motion estimate with a standard error less than half of the requirement for  @xcite .
therefore vlbi appeared to be at that time the only option to determine the motion of a suitable guide star with the desired accuracy . *
the star needed to be sufficiently bright in the radio , preferably with a flux density at least of order 1  mjy , for vlbi astrometry measurements to be successful . * the star needed to be sufficiently bright in the optical , at least 7th magnitude , for the on - board telescope to detect . *
the star s angular distance , @xmath8 , from the ecliptic needed to be between 20@xmath9 and 40@xmath9 .
the lower limit was set by the requirement that , with a suitable sun shield at the telescope , sunlight not enter the dewar and boil off the liquid helium .
the upper limit was set by the power requirements of the spacecraft , since the solar panels were fixed to the spacecraft , which would continually point to the guide star , making the power received by the solar panels dependent on @xmath8 . *
the star needed to be located as close as possible to the celestial equator , so as to maximize the sensitivity of the measurement of the frame - dragging precession , which essentially decreases with the cosine of the star s declination . * the star needed to be sufficiently isolated on the sky .
neighbouring stars or reflection nebulae could possibly intolerably decrease the accuracy of the pointing of the telescope to the star s center of brightness .    to find the most suitable guide star for , we conducted a vla survey at 8.4 ghz from 1990 to 1992 of @xmath41200 stars with v magnitude @xmath76.0 and a declination between @xmath10 and @xmath11 .
no single stars were found in the survey with a flux density above @xmath41  mjy .
in fact , we detected only previously known radio stars , all of them binaries of the rs canum venaticorum type ( see also , for a later survey , * ? ? ?
we identified four potentially suitable guide stars : @xmath12  andromeda  ( hr  8961 , @xmath13 declination ) , hr  1099 ( @xmath14 ) , hr  5110  ( @xmath15 ) , and ( hr  8703 , @xmath16 ) . after investigating each of the four candidates with vlbi , studying the sky fields around them in the optical , and checking on possible vlbi reference sources for each of them , we selected  as the most suitable for
.      is a binary with a giant primary and sun - type secondary .
we summarize its optical characteristics and previously known properties in table [ topt ] .
the radius of the primary is @xmath413.3 @xmath3 0.6 , which translates at a distance of  of @xmath4100 pc , to an angular radius of @xmath40.64 @xmath3 0.03 mas .
this radius is comparable to the full - width at half - maximum ( fwhm ) of the synthesized beam of a global vlbi array operating at 8.4 ghz , and therefore in principle allows for an accurate determination of the location of the radio emission with respect to the center of the primary .
a prerequisite for successful guide star tracking was that the onboard telescope be able to lock sufficiently accurately on the center of the optical disk of the star ( or a point with an offset from this center constant during the mission ) .
one concern here was the appearance of variable dark spots on the primary , shown by doppler imaging and photometry to cover @xmath415% of the star s surface .
if the variations in the size or location of these spots resulted in a linear trend over the course of the  mission in the centroid of optical brightness , relative to the center of the disk of the i m peg primary , then that trend would have the same effect on the  data reduction as an equal error in our vlbi measurement of the proper motion of i m peg .
such a trend is not necessarily implausible , given the year - to - year variations in i m peg s optical brightness and the astrophysical plausibility of variations in typical spot latitudes analogous to the systematic variations of sunspot locations on the sun .
however , @xcite showed , through extensive spot mapping spanning the duration of the  mission , that any error due to variability in the spot pattern would yield no more than a  0.04  drift in the optical center .
a second concern was the surroundings of i m peg on the sky .
any star or nebula within the field of view of the telescope when locked on , possibly combined with s photometric variability , could cause systematic astrometric errors .
however , for , no other star within , e.g. @xmath17 , is brighter than v magnitude 10 , which can be compared with the corresponding relative large brightness of , that varied during the mission only between v magnitude 5.7 and 6.0 . all in all , through extensive optical and even millimetric observations of the co(j=1@xmath180 ) line ( the latter to search for any molecular cloud that could be associated with a reflection nebula ) , we constrained any astrometric systematic error to a negligible value : the combined error from dark spots , neighbouring stars , and a nebula , combined with the photometric variability of , we estimated to be smaller than 0.05  .
l c c@  c c c    _ hipparcos _ , @xmath19 , epoch 1991.25 & & & & 1 + _ hipparcos _ , @xmath20 , epoch 1991.25 & & & & 1 + _ hipparcos _ , parallax ( mas ) & & & & 1 + _ hipparcos _ , distance ( pc ) & & & & + & & & & & +   + & & & & & + v magnitude range during mission & & 5.7 to 6.0 & & & 2 + mass ( @xmath21 ) & & @xmath22 & @xmath23 & & 3,3 + spectral type & & k2  iii & g  v & & 4 + @xmath24 ( k ) & & @xmath25 & @xmath26 & & 4,3 + radius ( @xmath27 ) & & @xmath28 & @xmath29 & & 4,3 + radius ( mas ) & & @xmath30 & @xmath31 & & 4,3 + & & & & & +   + & & & & & + @xmath32 ( @xmath27 ) & & @xmath33 & @xmath34 & & 3,3 + @xmath32 ( mas ) & & @xmath35 & @xmath36 & & + @xmath37 ( days ) & & & & 3 + @xmath38 ( @xmath9 ) & & & & 4,5 + @xmath39 & & & & 4 + @xmath40 ( hjd ) & & & & 3 + [ topt ]
l l l l l l l 3c  454.3 & quasar & & & 7  10 & 0.859 & 1610 + b2250 + 194 & galaxy & @xmath41 & 3.6 & 0.35  0.45 & 0.28 & 880 + b2252 + 172 & unidentified & 0.4 & 1.4 & 0.017 & & +   + i m peg & rs cvn & @xmath42 & 0.7 & 0.0002  0.08 & 0.0 & 0.0 + [ t2refsources ]
our goal was to determine the motion on the sky of  relative to the distant universe , with a standard error not exceeding 0.14  in either coordinate . to achieve this goal , we selected radio reference sources with respect to which the motion of  could be measured .
these sources needed to be at cosmological distances , be sufficiently strong and compact for vlbi observations , and located nearby to  on the sky so as to minimize systematic astrometric errors .
we selected three reference sources , with the strong quasar , , as the main one .
the other two were the active nucleus of the radio galaxy , , and the unidentified source ( most likely also a quasar or a radio galaxy ) , .
we display their positions on the sky , together with that of , in figure  [ f1skypos ] , and list their characteristics in table  [ t2refsources ] .
( east - west ) and @xmath43 ( north - south ) directions on the plot are on the same scale .
figure taken from . ]    in figure  [ f2images ] we present typical vlbi images at 8.4 ghz of  and each of the three references sources .
contour ( fwhm ) in the inset .
north is up and east to the left . for more information ,
see  and .
figure taken from , but rearranged.,title="fig : " ]   contour ( fwhm ) in the inset .
north is up and east to the left . for more information ,
see  and .
figure taken from , but rearranged.,title="fig : " ]   contour ( fwhm ) in the inset .
north is up and east to the left . for more information ,
see  and .
figure taken from , but rearranged.,title="fig : " ]   contour ( fwhm ) in the inset .
north is up and east to the left . for more information ,
see  and .
figure taken from , but rearranged.,title="fig : " ]    all the sources are approximately located along the same axis ( north - south ) , allowing perhaps for a more accurate relative correction of tropospheric and ionospheric effects . however , each of the three sources has advantages and disadvantages as reference sources for .
the quasar 3c  454.3 is both the strongest and the closest on the sky to .
it was used as a reference source for  from 1991 to 1994 @xcite , and could therefore extend the time baseline of our vlbi observations .
the disadvantage is the complexity of 3c  454.3 s structure and its superluminal nature . in figure  [ f2images ]
we identified six components of which the easternmost component , c1 , is still compact at the highest angular resolution , has a flat or inverted high - frequency spectrum @xcite , and is located almost exactly at the same position ( ) as the core of the quasar seen at 43 ghz @xcite .
the component c1 was therefore considered to be the core at 8.4 ghz , the part most closely related to the supermassive black hole and to the center of mass of the quasar ; it was thus the reference point in the quasar s structure for our astrometric measurements .
the other components are moving away from c1 , with apparent speeds of up to @xmath45  c ( ) .
the remaining sources were both rather compact , so that the brightness peak could simply be taken as the reference point in the structure ; their disadvantages were that b2250 + 194 was located relatively far away from  and b2252 + 172 was rather faint ( with an unknown redshift ) .
in addition , we used geodetic and astrometric vlbi observations regularly conducted over many years to determine the limit on the proper motion of two of the three reference sources , 3c  454.3 and b2250 + 194 .
these determinations were made without consideration of a particular reference point in the structure of the two sources but were made with respect to a celestial reference frame , crf .
this frame is defined by a large number of extragalactic radio sources distributed over the sky and is for our purposes essentially the same as the international celestial reference frame , icrf2 @xcite .
all in all , we took c1 in 3c  454.3 as the primary reference point with respect to which the motion of  was measured .
the other measurements served to determine the stability of c1 on the sky so that finally the motion of  could be determined with respect to the distant universe .
to determine the motion of  relative to c1 and test c1 s stability relative to the brightness peaks of the other two reference sources , we took into account :    *  is catalogued as a binary .
however , we could not rule out that it is perhaps part of an extended triple or even larger multiple system . therefore , to ensure our obtaining a sufficiently accurate value of the proper motion of  for the @xmath415 months  data - collection interval in 2004 and 2005
, we spread a sufficiently large number of vlbi measurements over eight and a half years to be able to determine any nearly constant change of the proper motion with time (  proper acceleration " ) , or place a limit on it . * for the best determination of s parallax , we spread the vlbi measurement epochs appropriately to cover the different seasons each year .
similarly , to allow us to model any orbital component of the motion of the radio emission , we carefully distributed the epochs over the phase of the binary orbit , which could be computed from earlier spectroscopic results . * to measure the time - variable brightness distributions of  and the three reference sources , and in particular to identify c1 in 3c 454.3 , we selected each vlbi observation session to be long enough to ensure that the  coverage was sufficiently dense for us to obtain excellent imaging of the sources . * since we expected  to fluctuate strongly in radio brightness , and since we knew that radio emission from b2252 + 174 was relatively weak , we needed our vlbi array to be sufficiently sensitive to nearly guarantee detection of the sources in each observing session . *
lastly , we needed an observing frequency for which receivers were available at each of the telescopes of the vlbi array , and which guaranteed sufficiently high angular resolution and low system temperatures .
all said , 35 sessions of astrometric vlbi observations were made between 1997 january 16 and 2005 july 16 with about four sessions every year , each about 11 to 15 h long .
the observations were made with a global array of 12 to 16 radio telescopes comprised of most or all of : mpifr s 100 m telescope at effelsberg , germany ; nasa / caltech / jpl s 70 m dsn telescopes at robledo , spain , goldstone , ca , and tidbinbilla , australia ; nrao s ten 25 m telescopes of the vlba , across the u.s.a ; nrao s phased vla , equivalent to a 130 m telescope , near socorro , nm ; and , at early times , nrcan s 46 m algonquin radio telescope near pembroke , ont , canada , and nrao s 43  m telescope in green bank , wv .
the observing frequency was 8.4 ghz .
all of these vlbi observations were made to extract interferometric phase delays related to the difference of arrival times of a radio wave from a celestial source at each pair of antennas of a vlbi array .
such observations allow the most accurate astrometric measurements @xcite and can yield relative positions and proper motions referenced to particular points in the structure of the sources with uncertainties as low as @xmath410   and @xmath410  , respectively ( for the earliest such measurements , see , e.g. , marcaide & shapiro 1983 ; bartel et al . 1986 ; and for a recent review , see reid & honma 2014 )
.    the phase delay is given by @xmath44 where @xmath45 is the interferometric phase or  phase " at observing angular frequency , @xmath46 , and time , @xmath47 .
the phase can be expressed as @xmath48   + 2\pi n(\omega , t ) , \label{phi}\ ] ] where @xmath49 is the `` geometric delay , '' the difference in the arrival times of the radio wave in vacuum at the two antennas ; @xmath50 describes the difference in the instrumental delays ( including clock behavior ) at the two antenna sites ; @xmath51 gives the difference in radio wave propagation times to the two antennas due to all tropospheric and ionospheric effects ; @xmath52 is the delay contribution from source structure to account for any non - pointlike brightness distribution of the celestial source ; @xmath53 is the ( thermal ) noise contribution to the phase measurement ; and @xmath54 describes the integer number of 2@xmath55 ambiguities , or `` phase wraps , '' in the measurement .
the astrometric information of a vlbi measurement is given by the geometric delay which , apart from relativistic contributions and those due to earth s motion relative to the solar system barycenter , is given by @xmath56 , \label{geom}\ ] ] where @xmath57 is the speed of light in vacuum , @xmath58 is the 3-dimensional vector between two antennas of a baseline of a vlbi array , and @xmath59 is the unit vector in the direction of the observed source .    for phase - delay vlbi observations to succeed , to be able to remove the 2@xmath55 phase ambiguities , and to reduce many sources of astrometric errors , we needed to switch the antennas rapidly between sources .
for the first 23 sessions , we used a typical sequence of 3c  454.3 ( 80  s ) - i m  peg ( 170  s ) - b2250 + 194 ( 80  s ) . for the remaining 12 sessions , b2252 + 172 was included in the observations and the sequence was altered to 3c  454.3 ( 80  s ) - i m   peg ( 125  s ) - b2250 + 194 ( 80  s ) - 3c  454.3 ( 80  s ) - i m  peg ( 125  s ) - b2250 + 194 ( 80  s ) - b2252 + 172 ( 90  s ) .
in addition to the astrometric 8.4 ghz observations , we also observed once at 5.0  ghz and once at 15.0  ghz , to allow investigation of the compactness and spectral properties of the components of the sources .
almost all observations were recorded in both right and left circular polarization , and all were processed with the vlba correlator at socorro , nm .
we used a custom software package ( ) that allowed us to analyze the data efficiently .
first , we connected the phases and removed the 2@xmath55 phase ambiguities to generate the correct phase delays .
then we corrected the phase delays for the effects of the troposphere and the ionosphere with _ a priori _ estimates .
we further processed them with a kalman filter to model the residuals of the troposphere and the ionosphere , and the clock offset at each telescope from the one chosen as the reference clock for the vlbi array .
then , for the more extended source 3c  454.3 , the structure effects were removed from the phase delays based on the reference point , c1 , defined with images made with nrao s aips software package .
the other two reference sources were sufficiently compact that no correction for their source structure was needed .
the relative weakness of s radio emission , and sometimes also that of b2252 + 172 , made it necessary for us to develop a special technique for the analysis of  vlbi data .
it combines the advantages of parametric model - fitting via weighted least - squares applicable for the relatively strong sources 3c 454.3 and b2250 + 194 with the sensitivity of phase - referenced mapping , the latter required for  and b2252 + 172 .
the merged analysis technique yielded superior results to the phase - referenced mapping technique ( ) .
the radio emission from  varied from session to session , reflecting the complex astrophysical nature of the environment of the star . in most cases
the images of  showed only one clearly defined component . in these cases
the coordinates were determined for the maximum of a two - dimensional gaussian fit to that component . in the remaining cases
, there was more than one local maximum of approximately equal brightness in the emission region : in eight cases there were two such maxima , and in one case there were three , all of them significant .
in such cases , we took the position of i m peg to be the mean position of the different maxima .
the final result of the analysis of the  vlbi observations was a set of coordinates with statistical standard errors for each source and each observing epoch .
these sets were the basis for the study of the motion , or upper bounds on it , of the sources with respect to each other .
in addition , we analyzed geodetic astrometric vlbi observations and made use of interferometric group delays .
group delays are given by @xmath60 these are observations of thousands of sources from all over the sky .
these observations are done on a routine basis , independent of the  program . for 3c  454.3
we used all available data from a total of 1,119 observing sessions , from 1980 to 2008 . in support of , the second reference source , b2250 + 194 ,
was included in 38 sessions of routine geodetic astrometric group - delay observations between 1997 and 2008 .
group - delay observations have the advantage that 2@xmath55 ambiguities are not inherent in the data and that the positions are determined relative to a celestial reference frame ( crf ) .
the disadvantage , however , is that the parameter , @xmath52 , was not determined and that therefore the exact reference point in the brightness distribution of the sources for the position measurement remained undefined .
all available geodetic astrometric data with a total of 6.5 million group - delay determinations were processed via least - squares and yielded estimates of coordinates with statistical standard errors for 3c 454.3 and b2250 + 194 for each observing session , keeping the coordinates of the other sources constant .
this set of solutions forms our crf , which for s purposes is virtually identical ( see , e.g. , * ? ? ?
* ) to the icrf2 @xcite , the most fundamental such frame presently in use .
the phase - delay determinations of the position of c1 , the core of 3c 454.3 , relative to those of the brightness peaks of b2250 + 194 and b2252 + 172 , are plotted in figure  [ f3posref ] .
the group - delay determinations of the position of b2250 + 194 relative to the crf are also plotted in figure  [ f3posref ] .
see  for a similar plot of the group - delay determinations of the position of 3c 454.3 relative to the crf .     are also given .
bottom panel : the relative coordinates of b2250 + 194 as determined from routine geodetic group - delay vlbi observations of up to @xmath44000 extragalactic sources scattered over the sky .
figure taken from . , title="fig:",scaledwidth=80.0% ]   are also given .
bottom panel : the relative coordinates of b2250 + 194 as determined from routine geodetic group - delay vlbi observations of up to @xmath44000 extragalactic sources scattered over the sky .
figure taken from .
, title="fig:",scaledwidth=72.0% ]    weighted least - squares straight line fits to each set of position determinations yielded the position at epoch and proper motion for each coordinate separately , along with their statistical standard errors and correlations . in each case
the fits to the data gave @xmath61 ( @xmath62 per degree of freedom ) larger than unity , indicating systematic errors in the position determinations .
these errors can for instance be due to uncorrected effects from source structure and/or atmospheric and ionospheric variations , and are difficult to quantify .
therefore the total standard errors of the position determinations were determined partly empirically by including a particular constant for each fit .
this constant was added in quadrature to the statistical standard errors for each set of position determinations and for each coordinate separately so as to obtain @xmath63 for each fit .
the resultant proper - motion values with their standard errors are listed for each set of position determinations in the panels of figure  [ f3posref ] .
the upper panel of figure  [ f3posref ] indicates that there appears to be some motion of c1 relative to b2250 + 194 at the 3.5@xmath64 significance level .
the source 3c 454.3 has a long jet , indicating strong activity .
some motion of a component like c1 even if located close to the black hole could therefore be expected .
no significant motion is discernible for b2250 + 194 relative to b2252 + 172 .
no matter whether there is some motion due to activity near the black hole or not , the 1@xmath64 upper limits of the proper - motion values for c1 relative to b2250 + 194 and b2250 + 194 relative to the crf are each smaller than 40 .
combining the phase - delay proper - motion determination of c1 relative to b2250 + 194 with the group - delay proper - motion determination of b2250 + 194 relative to the crf gives for the time period from 1998 to 2005 a 1@xmath64 standard error for the mean position of c1 in the crf of 45 and 68 .
the proper motion of c1 relative to the crf for that same period is @xmath65  and @xmath66  in  and , respectively , which we do not consider significant .
we therefore give a 1@xmath64 upper limit on the magnitude of the proper - motion components of 46 and 56  , respectively .
the latter pair of determinations is our limit on the level of stationarity of the reference point c1 in 3c 454.3 relative to the distant universe ( ) .
with the stationarity of the reference point in 3c 454.3 determined , the motion of  relative to the distant universe could be analyzed .
five kinds of motion needed to be considered :    * proper motion mainly due to the motion of  relative to the motion of the sun in the galaxy , * annual parallax due to the orbital motion of earth around the sun , together with the much smaller contribution of the motion of the sun with respect to the solar system barycenter , * change of proper motion over time due to a hypothetical third member of the  binary system , * orbital motion around the common center of mass of the binary , and * erratic motion of the center of radio emission within the binary system .      to determine the proper motion , parallax , proper acceleration , and orbital motion , we used weighted least - squares to fit a linearized model to the 35 determinations of the positions of the reference point in  relative to the position of c1 in 3c 454.3 .
the uncertainties of these position determinations are hard to compute theoretically .
we therefore took them from an empirical analysis based on the computation of the root - mean - square ( rms ) scatter separately in @xmath67 and @xmath43 , of the postfit residuals obtained for the 12 positions of  relative to c1 and the 35 positions of  also relative to c1 .
the resulting uncertainty was 0.06 mas in each coordinate .
in addition there is the previously discussed uncertainty of the position of c1 and the upper limit of the proper motion of c1 in the crf that we considered for the standard errors of the position determinations of  relative to the crf .
the parameters of the model fit were the position of  at epoch , proper motion , parallax , and four scalar parameters for the description of the projection on the sky of an orbit with zero - eccentricity @xcite and a period known from optical spectroscopic observations @xcite .
the fit resulted in rms values of the residuals somewhat different for each coordinate but much larger than the standard errors of the position determinations of .
the source of this large scatter is the astrophysical nature of the radio emission , characterized by motion related to radio brightness changes @xcite and plausibly in part to stellar magnetic field changes implied by spot maps from optical spectroscopy ( e.g. , * ? ? ?
for the fit we used uniform weighting independently for  and , and allowed their errors to have non - zero correlation .
after iteration , we obtained the final set of parameter estimates with statistical standard errors and @xmath63 .
this set is presented in table  [ tfinal ] .    for the determination of the proper acceleration
, we enlarged the set of positions of 35 epochs by four additional positions of   obtained at epochs between 1991 and 1994 , also at 8.4 ghz , by @xcite in support of the _ hipparcos _ mission .
these observations considerably extended our time baseline ; however , they were only of limited use for the estimate of parameters other than the proper acceleration because of inferior  coverage , lower angular resolution , and the position determinations not being referred to c1 .
we found no significant proper acceleration ; i.e. , none larger than one statistical standard error ; and so we include no such proper acceleration in the fit used to obtain our final results in table  [ tfinal ] .
figure  [ f4proppar ] displays the 39 position determinations and the fit to the 35 position determinations , with the four earliest positions not used in the fit .
the effect of the proper motion and the parallax can be clearly seen in the data as well as in the fit .
figure  [ f5paraellipse ] displays the effect of the parallax still more clearly .
it shows the 35 position determinations with the estimated position at epoch , proper motion , and the primary s orbital motion subtracted .
it also shows the estimated parallax ellipse .
figure  [ f6morbit ] displays the orbit clearly .
again we plot the 35 position determinations with the estimated position at epoch , proper motion , but now with parallax rather than orbital motion subtracted . also shown
are the estimated orbit and the corresponding positions predicted from the orbit model .
figure  [ f7resvt ] shows the 39 position determinations and the fit to the above 35 positions , but now with all model contributions , namely the estimated position at epoch , proper motion , parallax , and the primary s orbital motion subtracted .
we compare them to the corresponding flux densities of .
the scatter of the position residuals appears to be random with certainly no correlation to the measured flux densities .
the rms of the scatter is @xmath40.4 mas in each coordinate , much larger than the standard errors of the position determinations and most likely dominated by the fluctuation of the positions of the stellar radio emission relative to the center of the primary of .
figure  [ f8imschemradio ] gives an artist s three - dimensional rendition of  with the primary as a giant with dark spots and the secondary as a sun - like star , each in its estimated orbit .
errors ( here and hereafter in the figure captions in the sense of the full length of @xmath68 error bars ) are much smaller than the symbols .
the curve shows the fit of the solution , with orbital terms omitted for clarity .
in addition , the four positions from 1991 to 1994 are shown ( upper left ) ; however , they were not used in the fit .
the arrow indicates the direction and its length the magnitude of the proper motion in one year .
figure taken from . ]
errors ( of position determinations relative to c1 ) are about a third of the size of the circles .
the relatively large scatter is most likely dominated by the fluctuations of the positions of the stellar radio emission relative to the center of the disk of the  primary .
the model positions ( not shown ) are almost exactly on the plotted parallax ellipse and close to the observed positions .
the parallax ellipse , however , is slightly time - dependent due to the influence of planets . over the course of our observations of several years , the parallax ellipse in the figure shifts by about the width of the line .
, scaledwidth=70.0% ]     errors ( for the position determinations relative to c1 ) are about twice the size of the circles .
a solid line connects each observed position with the corresponding fitted position indicated by a dot on the estimated orbit .
the large scatter of the observed positions about the model orbit ellipse is again most likely dominated by the fluctuations of the positions of the stellar radio emission relative to the center of the disk of the  primary .
figure adapted from .
, scaledwidth=70.0% ]    .
the secondary is the smaller yellow star .
the projected orbit of the primary is the same as the inferred radio source orbit shown in figure  [ f6morbit ] .
the projected orbit of the secondary and the diameters of the two stars correspond to the nominal values given in table  [ topt ] .
the system is shown with the primary at its ascending node .
note that the primary s visible pole is near its southeast side .
the size of the figure is approximately 3.8  mas by 3.4  mas . , scaledwidth=95.0% ]      a detailed analysis of systematic errors affecting each of the nine fit parameters of s motion relative to the distant universe led us to the following conclusions .
for the position at the mean epoch of the vlbi observations of 2001.29 , the systematic error in each coordinate was computed from 0.5 times the angular radius of the primary along the sky projection of the normal to the binary orbital plane ( p.a .
= 130.5   @xmath3 8.6 ) . for the later epoch of 2005.08 in table  [ tfinal ]
, the value for each component was added in quadrature with a possible maximum rms drift during that time interval of the radio brightness reference point of .
the rms drift rate was assumed to be not larger than one stellar radius over the vlbi observation time span of 8.5 years or 0.06  in @xmath67  and 0.05  in @xmath43 .
the systematic error of the proper motion estimate was found to be mostly comprised of the upper limit of the proper motion of c1 in the crf and of the aforementioned possible maximum rms drift of the radio brightness reference point of .
the systematic errors of the parallax and orbit estimates were found to be most clearly defined by the position values of b2250 + 194 relative to those of c1 in 3c 454.3 , since the true parallax and orbit parameter values should be zero .
the same parameter fit used for  resulted in a parallax estimate for b2250 + 194 relative to c1 of @xmath69 @xmath3 0.074 mas , with the statistical standard error decreasing to 0.026 mas for @xmath63 .
for the five - fold smaller separation of  from 3c  454.3 the systematic uncertainty was estimated to be @xmath700.015  mas in each coordinate .
the same parameter fit also resulted in orbit parameter estimates for b2250 + 194 relative to c1 much smaller than the estimates for  itself .
for each of the parameters in table  [ tfinal ] the total error is also listed . in general
it is computed as the root - sum - square of the statistical standard error and the estimated systematic error . for the position and proper motion values , however , the statistical errors are doubled before computing the root - sum - square to allow for correlated noise in the vlbi positions . for the parallax and orbit parameter estimates ,
the systematic errors were too small to contribute appreciably to the total standard errors .
l c c c c @xmath67 at epoch 2005.08 ( errors in mas ) & & 0.12 & 0.33 & 0.40 + @xmath43 at epoch 2005.08 ( errors in mas ) & & 0.13 & 0.29 & 0.39 + @xmath71 ( ) & @xmath72 & 0.026 & 0.073 & 0.090 + @xmath73 ( ) & @xmath74 & 0.030 & 0.074 & 0.095 + parallax ( mas ) & @xmath75 & 0.074 & @xmath700.015 & 0.074 +  ( mas ) & @xmath76 & 0.10 & @xmath770.1 & 0.10 +  ( mas ) & @xmath78 & 0.11 & @xmath770.1 & 0.11 +  ( mas ) & @xmath79 & 0.09 & @xmath770.1 & 0.09 +  ( mas ) & @xmath80 & 0.11 & @xmath770.1 & 0.11 + semimajor axis ( mas ) & 0.89 & 0.09 & @xmath770.1 & 0.09 + axial ratio & 0.30 & 0.13 & @xmath770.1 & 0.13 + p.a .  of ascending node ( deg ) & 40.5 & 8.6 & @xmath778 & 8.6 + @xmath40 ( heliocentric jd ) & 2450342.56 & 0.44 & @xmath770.4 & 0.44 +      comparing our results with the most precise previous measurements of proper motion and parallax
, we find that our estimates agree with those in the _ hipparcos _ catalogue @xcite and those of @xcite within their respective larger standard errors .
however , we found that our estimates slightly disagree with the values in the _ hipparcos _ re - reduction @xcite by 1.6 and 2.4 times the combined standard error in  and in parallax , respectively .
this discrepancy in  is , however , more than 10 times smaller than the standard error of either the geodetic or the frame - dragging effect @xcite , and therefore of no consequence for the  results .    comparing the orbital parameters
, we find that our estimate of the time of conjunction is consistent with that of @xcite and the combination of our estimates of orbital inclination and semimajor axis are consistent with their value of @xmath32 .
this agreement , in turn , is consistent with the reasonable expectation that the radio emission is on average centered on the primary of the  system , and that it orbits with the same inclination .
any offset in phase of the radio orbit from that of the primary corresponds to less than one - forth of the radius of the primary .      with position at epoch , proper motion , parallax , orbital motion , and average location relative to the center of the primary determined , how can the radio emission be described and where does it appear relative to the primary of  from epoch to epoch ?
figure  [ f9image ] displays a set of vlbi images of  relative to the center of the primary in its orbit about the binary s barycenter .
the rms background ; and contours at 50% and above are drawn in black .
the center of each panel ( red dot ) is the fit position of the star s center , as derived from the astrometric results in . in other words ,
the coordinate origin in our radio image should approximately represent the center of the disk of the primary star .
the red circle indicates the angular size of the primary star ( radius of @xmath28  ; * ? ? ?
the cyan dotted ellipse shows the orbit of the primary as in  and .
figure taken from , where also the complete set of our vlbi images is given .
, scaledwidth=80.0% ]    the analysis of the images showed that the radio emission is highly variable and frequently partly circularly polarized .
the morphology is also variable with the average size of the radio emission slightly larger than the disk of the primary .
the positions of the peaks of the emission regions are scattered over an area on the sky slightly larger than the disk of the primary .
the scatter has a tendency to be distributed along the orbit normal and expected primary s spin axis .
comparison with simulations suggest that the brightness peaks preferentially occur near the polar regions similar to the dark spots in the optical @xcite .
the height of the emission is likely close to the surface of the primary with 2/3 of the emission peaks located within 0.25 times the stellar radius above the surface .
the radio emission may be due to flares linked to a possibly dipolar magnetic field , whose axis is normal to the orbit plane , as could be expected giving the tidally locked rotation of the primary . for a movie of the star based on these vlbi images , see the website , www.yorku.ca/bartel/impeg.mpg , or the online version of .
our series of 35 vlbi sessions between 1997 and 2005 resulted in the most comprehensive and detailed radio investigations ever made of a star .    *
the proper motion of  relative to the distant universe is @xmath0 @xmath3 0.09  in  and @xmath2 @xmath3 0.09  in , where the errors are intended to represent one - standard deviation errors and include estimates of systematic errors . *
these results met the pre - launch requirements of the  mission of standard errors no larger than 0.14  in each coordinate so as not to discernibly degrade the estimates of the geodetic and frame - dragging effects . *
the parallax of  is 10.37 @xmath3 0.07  mas , corresponding to a distance of 96.4 @xmath3 0.7  pc . *
the proper - motion and parallax estimates agree with previous estimates , including those in the _ hipparcos _ catalogue , within the combined standard errors , with the exception of the revised _ hipparcos _ results of @xcite where the agreement is within 1.6 and 2.4 times the respective combined standard errors .
any discrepancy at these levels is however of no consequence for the  results . *
the radio emission of  is highly variable . *
the vlbi images of  show emission regions slightly larger than the disk of the primary and their peaks scattered about a similarly sized area , consistent with being located preferentially near the polar regions and within @xmath40.25 stellar radii above the surface . *
the vlbi images of  were assembled for a movie of the star ( see the website of the first author , www.yorku.ca/bartel/impeg.mpg , or the online version of ) .    acknowledgements .
we thank each of the following people , listed alphabetically , for their contributions of observations and data reductions to our investigations of  and its sky neighborhood , for providing us with unpublished updated catalog data , and/or for enlightening discussions on various aspects of this research program : s.  berdyugina , a.  buffington , d.  buzasi , n.  caldwell , r.  campbell , j.  chandler , t.  dame .
r.  donahue , g.  fazio , d.  fischer , o.  franz , g.  gatewood , d.  gordon , e.  guinan , p.  hemenway , a.  henden , m.  holman , the late j.  huchra , p.  kalas , g.  keiser , j.  kolodziejczak , r.  kurucz , b.  lange , d.  latham , j.  lederman , j .- f .
lestrade , j.  lipa , p.  luca , s.  marsden , h.  mcalister , d.  monet , n.  nunes , l.  petrov , r.  schild , g.  torres , w.  traub , d.  trilling , n.  turner , the late h.  wendker , and the late c.  worley , this research was primarily supported by nasa , through a contract from stanford university to sao , with a major subcontract from sao to york university .
we obtained observations made with the 100-m telescope of the max - planck - institut fr radioastronomie at effelsberg , with antennas of the national radio astronomy observatory ( nrao ) , which is a facility of the national science foundation operated under cooperative agreement by associated universities , inc .
, with antennas of the dsn , operated by jpl / caltech , under contract with nasa , and with the nasa / esa hubble space telescope , which is operated by the association of universities for research in astronomy , inc .
, under nasa contract nas 5 - 26555 .
we made use of the sao / nasa astrophysics data system abstract service , initiated , developed , and maintained at sao .
, a.  l. , gordon , d. , & jacobs , c.  s. , ( eds . ) 2009 , iers technical note no .
35 , http://www.iers.org/iers/en/publications/technicalnotes/tn35.html ( frankfurt : verlag des bundesamts fr kartographie und geodsie )                        , a. , krichbaum , t.  p. , graham , d.  a. , et  al .
2004 , in european vlbi network on new developments in vlbi science and technology ( trieste : proceedings of science ) , ed .
r.  bachiller , f.  colomer , j .- f .
desmurs , & p.  de vicente , 710