Radio Interferometric Planet Search II: Constraints on sub-Jupiter-Mass Companions to GJ 896A

Geoffrey C. Bower11affiliation: Astronomy Department & Radio Astronomy Laboratory,University of California, Berkeley, CA 94720 , Alberto Bolatto11affiliation: Astronomy Department & Radio Astronomy Laboratory,University of California, Berkeley, CA 94720 22affiliation: Department of Astronomy, University of Maryland, College Park, MD, 20742-2421 , Eric B. Ford11affiliation: Astronomy Department & Radio Astronomy Laboratory,University of California, Berkeley, CA 94720 33affiliation: Astronomy Department, University of Florida, Gainesville, FL 32611-2055 , Adam Fries11affiliation: Astronomy Department & Radio Astronomy Laboratory,University of California, Berkeley, CA 94720 , Paul Kalas11affiliation: Astronomy Department & Radio Astronomy Laboratory,University of California, Berkeley, CA 94720 , Karol Sanchez11affiliation: Astronomy Department & Radio Astronomy Laboratory,University of California, Berkeley, CA 94720 , Phoebe Sanderbeck44affiliation: Department of Physics, Brandeis University, Abelson-Bass-Yalem 107, MS 057, 415 South Street, Waltham, MA 02453 , Vincent Viscomi11affiliation: Astronomy Department & Radio Astronomy Laboratory,University of California, Berkeley, CA 94720

We present results from the Radio Interferometric Planet (RIPL) search for companions to the nearby star GJ 896A. We present 11 observations over 4.9 years. Fitting astrometric parameters to the data reveals a residual with peak-to-peak amplitude of mas in right ascension. This residual is well-fit by an acceleration term of mas/y. The parallax is fit to an accuracy of 0.2 mas and the proper motion terms are fit to accuracies of 0.01 mas/y. After fitting astrometric and acceleration terms residuals are 0.26 mas in each coordinate, demonstrating that stellar jitter does not limit the ability to carry out radio astrometric planet detection and characterization. The acceleration term originates in part from the companion GJ 896B but the amplitude of the acceleration in declination is not accurately predicted by the orbital model. The acceleration sets a mass upper limit of 0.15 at a semi-major axis of 2 AU for a planetary companion to GJ 896A. For semi-major axes between 0.3 and 2 AU upper limits are determined by the maximum angular separation; the upper limits scale from the minimum value in proportion to the inverse of the radius. Upper limits at larger radii are set by the acceleration and scale as the radius squared. An improved solution for the stellar binary system could improve the exoplanet mass sensitivity by an order of magnitude.

astrometry, stars: activity, stars: early-type, stars: planetary systems, radio continuum: stars

1 Introduction

The study of extrasolar planets provides an important link between the study of star formation, the study of our own solar system, and the search for extraterrestrial life.  Radial velocity searches (Butler et al., 2006), transit searches with Kepler (Borucki et al., 2011), and microlensing searches (Gaudi et al., 2008) have surveyed significant components of the exoplanet mass-separation parameter space exposing a wide variety of phenomenology.

Despite their abundance, low mass stars have been less well-characterized for a number of reasons. In particular, low mass stars are often too faint or too strongly variable for high precision radial velocity and transit searches. Microlensing searches have identified companions to low mass stars (Bennett et al., 2008; Beaulieu et al., 2006; Gould et al., 2006), which suggest that such companions may be very common. On the other hand, radial velocity searches suggest very low planet occurrence rates for close-in companions (Endl et al., 2006; Johnson et al., 2007). The frequency and distribution of companions to low mass stars can be an important diagnostic of methods for planet formation (Laughlin et al., 2004; Ida & Lin, 2005; Boss, 2006; Kennedy et al., 2007).

Astrometry can be a powerful tool for detection and characterization of nearby, low-mass planetary systems. Muterspaugh et al. (2010) demonstrate that optical astrometric measurements with the Palomar Testbed Interferometer are capable of achieving sub-milliarcsecond accuracy in the separation between close binaries. They summarize the known properties of sub-stellar companions to low mass stars in binary systems and conclude that giant planet companions in such systems are common. The Carnegie Astrometric Planet Search Program (Boss et al., 2009) has also obtained sub-milliarcsecond astrometry on the low mass star NLTT 48256 over 2 years of observations; this program will ultimately survey 100 nearby M, L, and T dwarfs.

The Very Long Baseline Array (VLBA) can routinely achieve an astrometric accuracy of as in a single epoch at frequencies GHz and as at higher frequencies. The VLBA is capable of accuracies as high as 8 as under favorable circumstances (Fomalont & Kopeikin, 2003).   To use this technique the target source must have a sufficiently high brightness temperature to be detected by a high resolution radio interferometer.  For the VLBA, brightness temperatures must be K, requiring nonthermal emission.  Thus, the systems which can be studied are limited to the most active stars. VLBA astrometry has been important for a wide range of stellar measurements including parallax (e.g., Sandstrom et al., 2007; Menten et al., 2007; Loinard et al., 2007; Torres et al., 2009), galactic structure (Xu et al., 2006, 2011; Brunthaler et al., 2011), binary orbit characterization (Guirado et al., 2006), and the radio/optical reference-frame link (Lestrade et al., 1999). Recently, Forbrich & Berger (2009) detected a brown dwarf with the VLBA and explored the possibilities of companion detection. Bower et al. (2007) reviewed the prospects for future radio astrometric surveys with respect to other methods.

RIPL, the Radio Interferometric Planet Search, is a 4-year program using the VLBA and the 100m Green Bank Telescope (GBT) to search for the astrometic signatures of massive planets around nearby ( pc), low-mass (M dwarf) stars. We are observing a sample of 30 stars a total of 12 times over a 4-year period with the expectation of achieving milliarcsecond astrometric accuracy per epoch, and ultimately having the sensitivity to detect Jupiter-mass planets at 1 AU. This astrometric search is an important complement to radial velocity, transit, microlensing, and optical astrometric searches for planets on account of greater sensitivity to 1) low mass stars, 2) long-period planets, 3) active stars, and, 4) planets that may be imaged with extreme adaptive optics. Many RIPL stars are in binary systems, which may have higher planetary fractions.

In Paper I (Bower et al., 2009), we used Very Large Array (VLA) observations to identify the sample of stars appropriate for RIPL and we used the VLBA to characterize the detectability and short-term ( days) stability of the astrometric images of these stars. The short-term stability addresses a significant concern of radio astrometry that stellar activity will lead to large astrometric errors. The absence of a significant difference in the historical optical proper motions and the radio proper motion from this RIPL precursor survey permitted us to exclude brown dwarf companions to the four stars observed. RIPL observing began in October 2007 and is anticipated to end in late 2011.

In this paper, we present initial results from RIPL for the star GJ 896A. This star was also observed as part of the precursor survey in paper I. Together these observations demonstrate astrometric stability over a 4.9 year span. GJ 896A (EQ Peg A) is a variable star in a multiple system with GJ 896 B, C, and D. GJ 896B is its closest companion; the orbit has been characterized, but it is not well constrained (Heintz, 1984). GJ 896 A and B have spectral types M3.5 and M4.5. These two stars have been detected previously in the radio and seen to have point-like structure and strong radio variability (Pallavicini et al., 1985; Jackson et al., 1989; Benz et al., 1995; Gagne et al., 1998). In §2, we describe our observations, analysis, and results. In §3, we present our astrometric analysis of the data. In §4, we examine the constraints that we can place on companions to GJ 896A. In §5, we summarize our results.

2 Observations, Data Analysis, and Results

Observations of GJ 896A were obtained as part of the RIPL survey. RIPL observations spanned October 2007 through June 2011, and are ongoing. In each epoch, we observed two stars over a span of 8 hours. Fast switching with a phase reference calibrator and secondary (and sometimes tertiary) astrometric calibrators reduced total observing time on each star to approximately two hours. The total span from start to end of observation for a particular star within an epoch varied, depending on the relative coordinates of the stars and the allocated time slot. In some instances, the stars were observed sequentially in four hour blocks. In other instances, the stars were interleaved and the span was over the full eight hours.

Observations typically included all ten antennas of the Very Long Baseline Array as well as the 100 meter Green Bank telescope. Data were recorded at a rate of 512 Mbps in dual circular polarization mode with observing bands at 8.4 GHz using the standard VLBA mode v4cm-512-8-2. More precisely, this mode corresponds to four dual-polarization frequency bands centered at 8384.49, 8400.49, 8416.49, and 8432.49 MHz, each sampled at 32 Msamples/second with 2 bit sampling. Occasional observations of fringe calibrators were interspersed throughout each epoch. The switching cycle included repeated observations of 1 minute on the phase reference calibrator followed by 3 minutes on the star. Every 30 minutes short duration observations of the secondary and tertiary calibrators phase referenced to the primary calibrator were obtained.

Data were reduced with an AIPS (Greisen, 2003) pipeline that performed standard calibrations of amplitude and phase. Corrections for the atmospheric total electron content (TEC) and post-correlation updates to the Earth orientation parameters (EOP) were also included. All sources were imaged and deconcolved using the CLEAN algorithm. The typical synthesized beam for GJ 896A was 3.0 1.1 mas, oriented North-South. We fit Gaussian models to all sources using the AIPS task JMFIT.

In the case of GJ 896A, the primary phase reference calibrator is J2328+1929, which has an ICRF2 position that is accurate to 1.84 mas 1.58 mas 111 The assumed position for J2328+1929 is (23 28 24.874773, 19 29′58.03010″); these updated values are slightly offset from the assumed position of paper I by (0.254, 0.059) mas in right ascension and declination, respectively. All positions for other sources are relative to this position. The secondary calibrator is J2334+2010. No tertiary calibrator was used. We apply a correction for tropospheric phase gradients on the astrometric position of the star by removing variations in the position of the secondary calibrator. This correction is applied on an epoch-by-epoch basis. A more sophisticated approach using the AIPS task ATMCA can make corrections on shorter timescales. The magnitude of corrections to GJ 896A are mas. This approach may overestimate trophospheric errors slightly given the greater distance and higher declination of J2334+2010 from the primary calibrator. The close spacing and nearly linear arrangement of the star and calibrators, however, make this a good calibration set (Fig. 1). Typically, 3C 454.3 was observed as the fringe calibrator. As of June 2011, we have 8 good epochs for GJ 896A.

Both calibrators are dominated by a strong point source and show only weak evidence for non-point-like structure. Images of the calibrators are included in Paper I. The primary calibrator, J2328+1929, has a typical flux density of 60 mJy; the secondary calibrator, J2334+2010, has a typical flux density of 20 mJy. In Table 1, we summarize flux densities and positions for the secondary calibrator as a function of epoch; we also include the mean value from the three epochs of paper I. The positions for J2334+2010 have an rms variation of (0.11, 0.09) mas and a mean offset (0.06, 0.22) mas from the astrometric position obtained in paper I. That mean offset is within the error bars of the original measurement. These measures for the secondary calibrator provide an estimate of the astrometric accuracy of the experiment.

We summarize positions and flux densities for GJ 896A in Table 2 as a function of epoch. We also include measurements from Paper I. GJ 896A is strongly variable: flux densities vary from 0.140 to 9.512 mJy during the course of the experiment. The images of GJ 896A are predominantly point-like, as found in Paper I. The smaller errors on the flux density and the positions are reflective of the increased sensitivity and coverage of RIPL over the precursor experiment due to the doubling of the observing time, doubling of the bit rate, and addition of the GBT to the array. The statistical errors are significantly less than the systematic errors associated with the astrometry. For the analysis discussed below, we add in quadrature an error of 0.2 mas to each coordinate as an estimate of the systematic error.

3 Astrometric Fitting

In addition to the eight epochs from RIPL, we obtained three epochs in the RIPL precursor survey over three days in March 2006. We consider the data from all 11 epochs in our astrometric fitting here. The data span 4.9 years.

GJ 896A was not originally part of the Hipparcos catalog but it is included in the reanalysis of the Hipparcos data (van Leeuwen, 2007). We summarize the known astrometric parameters in Table 3. GJ 896A is also known to be in a double system with GJ 896B; parameters are determined from years of optical observations (Table 4; Heintz, 1984; Mason et al., 2001). Since, the baseline for this orbit determination is times the estimated period, the orbital parameters are poorly determined. This is reflected in the binary orbit grade of 5 (“indeterminate”) on a scale of 1 to 5.

In Table 3, we show the nominal astrometric data available from optical sources as well as best-fit solutions based on RIPL data. We performed two kinds of fits: one with just the astrometric parameters, and the other with astrometric parameters and an acceleration term. Figures 23, and 4 show the astrometric data, the models, and the residuals for the cases of the initial model, the astrometric model, and the astrometric plus acceleration model, respectively. The astrometric plus acceleration model that we fit to our data at epoch is


where is the equinox J2000 position at time , the epoch 2000.0, is the proper motion, are the projections of the parallax onto the two coordinates, (, ) is the local acceleration, and and are zero-points for the acceleration in each coordinate. The initial model clearly shows a substantial error in and a mas error in . Fitting the astrometric parameters (, , , , ) substantially improves the quality of the fit but leaves a parabolic residual in . An acceleration term reduces the rms in from to 0.26 mas. The residual in declination is improved from 0.33 to 0.26 mas through the addition of an acceleration term. The astrometric fit alone with 22 measurements and 5 parameters producing 17 degrees of freedom results in . The acceleration term significantly reduces to 29.7 with 9 parameters producing 13 degrees of freedom. The reduced in this final fit.

We estimate errors in individual parameters by mapping the surface as a function of each parameter. For the purely astrometric case, the curves are parabolic. We plot the curves for the astrometric and acceleration terms in Figure 5. These curves are also parabolic and symmetric about their best-fit values. We estimate the errors from these curves. The astrometric parameter estimates are significantly more accurate when the acceleration terms are included. We do not estimate the full covariance matrix of the parameter estimates due to the computational cost of exploring a 9-dimensional space. We did explore the covariance between the parallax and the right ascension acceleration term, two parameters that one might expect to be strongly linked. We find that the surface in this two-dimensional domain is well-characterized by a paraboloidal surface with axes that are misaligned from the principle axes by less than 30 degrees; that is, the parameters are essentially independent.

The estimated values for the astrometric parameters are determined one to two orders of magnitude more accurately than the Hipparcos values. The values for , , , and are in agreement with the Hipparcos values. The proper motion in right ascension, , however, differs by mas/y. Over the time baseline between the Hipparcos epoch (1991.25) and the mean RIPL epoch (2008.75), this implies an acceleration of mas/y. This has the same sign and a comparable magnitude to the acceleration term measured from the RIPL data alone. The agreement in implies an upper limit of 0.06 mas/y over this interval.

Motion of the star during the observing epoch can be an important effect, especially if the flux density is variable on short timescales. For GJ 896A, the median absolute value of motions from astrometric parameters are (0.10, 0.048) mas/hour in the two coordinates, respectively, and are always lower than (0.14, 0.065) mas/hour. This can lead to smearing out of the source and reduction of the peak flux density. If the source flux density is strongly variable during the observation, then the mean time of the epoch may not correspond to the time of the astrometric detection, leading to errors be as large as (0.6,0.3) mas. The small residuals we obtain from fitting indicate that this effect is not likely to be important in most epochs. Future reductions of this data will attempt to determine any light curve variability and its effect on the astrometry and/or include corrections to the visibilities for the assumed astrometric motion of the star. In high signal to noise ratio cases, we may be able to simultaneously fit a position and velocity within a given epoch.

4 Constraints on Companions to GJ 896A

The best fit accelerations for GJ 896A are (0.458, 0.087) mas/y. We compare this to the known parameters of the binary system. We plot relative separation, velocity, and acceleration determined from the binary parameters (Figure 6). Note that the definition of the orbit gives the position of B relative to A. Thus, the sign of velocity and acceleration terms is inverted relative to our measurement. If the stars have equal mass, then the motion relative to the center of mass will be half of the value plotted. The mean calculated accelerations for A relative to B over our observing epochs is (, ) mas/y. The calculated acceleration term in agrees reasonably well with the measured value and is comparable to the mean acceleration over the entire orbit mas/y. The calculated acceleration term does not agree at all with the measured value. The RIPL measurements occur near the predicted periastron, and, therefore, a maximum in the value of . Our results clearly indicate a significant error in the model for this binary system and appear to indicate that the period may be longer than 359 years, or that the eccentricity or time of periastron have not been properly estimated.

The acceleration allows us to a place a limit on the stellar companion mass, , using


where is the current separation projected onto the sky and is the mass of Jupiter. If we translate the measured acceleration to physical units, we find AU/y and AU/y, or the total acceleration AU/y. Applied to the stellar companion at a projected distance of 5.5″(the optically measured separation in 2007), we estimate a mass for GJ 896B of . This mass estimate is in rough agreement with expectations for a star of spectral type M4.5 (Kraus et al., 2011).

We can also apply the acceleration as a limit to estimate planetary mass companions. If we take the acceleration as an upper limit then we can place a limit on a planetary mass, . This acceleration limit applies to companions with a period y, or at AU for a stellar mass of 0.2 . For closer orbits, we can set a constraint based on the maximal amplitude of angular displacement. We set this limit by computing the periodogram for the residuals in each coordinate after subtraction of the astrometric plus acceleration model (Figure 7). The amplitude of the periodogram falls well below the 95%-confidence threshold; we find no evidence for a short-period companion. The maximal power in the periodogram corresponds to an upper limit of 0.4 mas in angular displacement. We are sensitive to orbital periods as short as the characteristic separation between epochs days; this limit translates to an inner radius cutoff of AU. In Figure 8, we combine the acceleration limits and the angular displacement limits. The maximal sensitivity is at 2 AU with an upper limit of 0.15 . Thus, we are sensitive to a planet with approximately 3 Neptune masses at large radii from GJ 896A. Note that this limit on the mass applies strictly to face-on systems; orbits at all inclinations will have a maximum acceleration equal to the face-on case but will have smaller instantaneous values depending on the phase of the orbit.

A caveat to this calculation is that the acceleration we measure may not be a upper limit to the effect of a planetary companion if the companion and GJ 896B produce accelerations that almost cancel each other. The agreement in between measured and calculated values suggests that this is not likely; the disagreement in could be due to a planetary companion acceleration that wholly or partially cancels the binary acceleration. A Jupiter mass companion at 2 AU would produce the magnitude of acceleration necessary for the cancellation. This cancellation, however, requires a coincidence that is unlikely to persist over a period of observation of this duration. The agreement in between the Hipparcos and RIPL measurements indicates that such a cancellation must persist over a period of y, which is improbable. Improved optical astrometry of the binary system can determine the contribution of GJ 896B to the motion of GJ 896A.

The rms residuals indicate that the centroid of the stellar emission is confined to a physical region . This is comparable to or less than the photospheric radius of a late-type star (e.g., Kraus et al., 2011). These residuals occur over a range of source flux densities that vary by a factor of . We conclude that the magnetic activity associated with GJ 896A is symmetric about the stellar surface and/or confined primarily to the disk of the star.

5 Conclusions

We have presented here an analysis of RIPL observations of the star GJ 896A. We measure the parallax to an accuracy of 0.2 mas and the proper motion to an accuracy of 0.01 mas/y. We measure an acceleration term of mas/y. This term is significantly less than expected based on the model for the orbit of the GJ 896AB system, suggesting an error in the optically-determined orbital parameters. Attributing this acceleration to a planetary mass companion, we demonstrate sensitivity to a planetary companion at 2 AU of mass 0.15 .

The residuals that we obtain from our fitting of 0.25 mas in each coordinate over 4.9 y indicate the feasibility of planet searching and characterization with radio astrometry of nearby low mass stars. These residuals can potentially be improved for this star and others through technical improvements to the data analysis including accounting for the motion of the star during the observing epoch and more exact removal of tropospheric fluctuations. The residuals indicate that astrometric jitter from stellar activity on scales of a stellar radius or larger is not a limiting factor in planet astrometry. Improvements in our analysis of systematic errors from stellar motion and tropospheric phase errors will permit us to explore sub-stellar-radius astrometric accuracy.

Future RIPL papers will explore analysis of additional data from this star, improved analysis techniques, and analysis of other stars in the sample. Upgrades to the recording bandwidth of the VLBA (and GBT) will produce significant improvements in the capability to carry out observing programs of this kind. The improved sensitivity will allow more rapid detection of bright stars such as GJ 896A, thereby reducing any smearing effect from stellar motion, allow the use of closer calibrators for improved tropospheric corrections and reduced time on calibrator, provide higher signal to noise ratio detections of fainter stars, and enable us to detect a larger sample of stars.

The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. This research was supported by NASA Origins grant NNX07AP19G.


  • Beaulieu et al. (2006) Beaulieu, J.-P., Bennett, D. P., Fouqué, P., et al. 2006, Nature, 439, 437
  • Bennett et al. (2008) Bennett, D. P., Bond, I. A., Udalski, A., et al. 2008, ApJ, 684, 663
  • Benz et al. (1995) Benz, A. O., Alef, W., & Guedel, M. 1995, A&A, 298, 187
  • Borucki et al. (2011) Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 728, 117
  • Boss (2006) Boss, A. P. 2006, ApJ, 643, 501
  • Boss et al. (2009) Boss, A. P., Weinberger, A. J., Anglada-Escudé, G., et al. 2009, PASP, 121, 1218
  • Bower et al. (2007) Bower, G. C., Bolatto, A., Ford, E., Kalas, P., & Ulvestad, J. 2007, ArXiv e-prints
  • Bower et al. (2009) Bower, G. C., Bolatto, A., Ford, E. B., & Kalas, P. 2009, ApJ, 701, 1922
  • Brunthaler et al. (2011) Brunthaler, A., Reid, M. J., Menten, K. M., et al. 2011, Astronomische Nachrichten, 332, 461
  • Butler et al. (2006) Butler, R. P., Wright, J. T., Marcy, G. W., et al. 2006, ApJ, 646, 505
  • Endl et al. (2006) Endl, M., Cochran, W. D., Kürster, M., et al. 2006, ApJ, 649, 436
  • Fomalont & Kopeikin (2003) Fomalont, E. B., & Kopeikin, S. M. 2003, ApJ, 598, 704
  • Forbrich & Berger (2009) Forbrich, J., & Berger, E. 2009, ApJ, 706, L205
  • Gagne et al. (1998) Gagne, M., Valenti, J., Johns-Krull, C., et al. 1998, in Astronomical Society of the Pacific Conference Series, Vol. 154, Cool Stars, Stellar Systems, and the Sun, ed. R. A. Donahue & J. A. Bookbinder, 1484–+
  • Gaudi et al. (2008) Gaudi, B. S., Bennett, D. P., Udalski, A., et al. 2008, Science, 319, 927
  • Gould et al. (2006) Gould, A., Udalski, A., An, D., et al. 2006, ApJ, 644, L37
  • Greisen (2003) Greisen, E. W. 2003, in Astrophysics and Space Science Library, Vol. 285, Astrophysics and Space Science Library, ed. A. Heck, 109
  • Guirado et al. (2006) Guirado, J. C., Martí-Vidal, I., Marcaide, J. M., et al. 2006, A&A, 446, 733
  • Heintz (1984) Heintz, W. D. 1984, AJ, 89, 1063
  • Ida & Lin (2005) Ida, S., & Lin, D. N. C. 2005, ApJ, 626, 1045
  • Jackson et al. (1989) Jackson, P. D., Kundu, M. R., & White, S. M. 1989, A&A, 210, 284
  • Johnson et al. (2007) Johnson, J. A., Butler, R. P., Marcy, G. W., et al. 2007, ApJ, 670, 833
  • Kennedy et al. (2007) Kennedy, G. M., Kenyon, S. J., & Bromley, B. C. 2007, Ap&SS, 311, 9
  • Kraus et al. (2011) Kraus, A. L., Tucker, R. A., Thompson, M. I., Craine, E. R., & Hillenbrand, L. A. 2011, ApJ, 728, 48
  • Laughlin et al. (2004) Laughlin, G., Bodenheimer, P., & Adams, F. C. 2004, ApJ, 612, L73
  • Lestrade et al. (1999) Lestrade, J.-F., Preston, R. A., Jones, D. L., et al. 1999, A&A, 344, 1014
  • Loinard et al. (2007) Loinard, L., Torres, R. M., Mioduszewski, A. J., et al. 2007, ApJ, 671, 546
  • Mason et al. (2001) Mason, B. D., Wycoff, G. L., Hartkopf, W. I., Douglass, G. G., & Worley, C. E. 2001, AJ, 122, 3466
  • Menten et al. (2007) Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515
  • Muterspaugh et al. (2010) Muterspaugh, M. W., Lane, B. F., Kulkarni, S. R., et al. 2010, AJ, 140, 1657
  • Pallavicini et al. (1985) Pallavicini, R., Willson, R. F., & Lang, K. R. 1985, A&A, 149, 95
  • Sandstrom et al. (2007) Sandstrom, K. M., Peek, J. E. G., Bower, G. C., Bolatto, A. D., & Plambeck, R. L. 2007, ApJ, 667, 1161
  • Torres et al. (2009) Torres, R. M., Loinard, L., Mioduszewski, A. J., & Rodríguez, L. F. 2009, ApJ, 698, 242
  • van Leeuwen (2007) van Leeuwen, F. 2007, A&A, 474, 653
  • Xu et al. (2011) Xu, Y., Moscadelli, L., Reid, M. J., et al. 2011, ApJ, 733, 25
  • Xu et al. (2006) Xu, Y., Reid, M. J., Zheng, X. W., & Menten, K. M. 2006, Science, 311, 54
Figure 1: Positions for GJ 896A and primary and secondary calibrators. The distance from the primary calibrator, J2328+1929, to GJ 896A is 0.9 deg. The distance from the secondary calibrator, J2334+2010, to GJ 896A is 0.6 deg. The nearly linear arrangement and small separations make this an excellent set of calibrators for this star.
Figure 2: Initial model and astrometric data for GJ 896A. Top row shows astrometric data (black crosses) with the initial (optically-determined) astrometric model (green circles at the epoch of observation and red lines for the continuous calculation). Positions in the top row are given in arcsec relative to a fiducial position of 2331, 1956’. Bottom row shows residuals after subtraction of the model. The symbols in the top row are much larger than the errors; in the bottom row, errors include statistical and systematic contributions.
Figure 3: Fitted astrometric model and astrometric data for GJ 896A. Symbols are the same as in Figure 2.
Figure 4: Fitted astrometric plus acceleration model and astrometric data for GJ 896A. Symbols are the same as in Figure 2.
Figure 5: as a function of each parameter. The x-axis values are offsets from the best-fit value. The parameters explored are (1) right ascension (hours), (2) declination (degrees), (3) proper motion in right ascension (mas/y), (4) proper motion in declination (mas/y), (5) acceleration in right ascension (mas/y), (6) acceleration in declination (mas/y), (7) minimum time for the right ascension acceleration (days), (8) maximum time for the declination acceleration (days), and (9) parallax (mas).
Figure 6: Orbital motion for GJ 896B relative to GJ 896A based on optical astrometric data. We plot the relative separation, velocity, and acceleration in the upper left, upper right, and lower left, respectively. We plot in the lower right the historical data (black crosses) from which this orbit derived (Heintz, 1984) and the predictions of the orbital model (green circles at the time of measurement, red line for continuous time). We include an additional datum from 2007 that was not part of the fit (Mason et al., 2001). The plotted model curve extends from 1932 to 2017.
Figure 7: Periodogram power as a function of frequency for residuals after removal of astrometric and acceleration terms for the two coordinates. We also plot the 95% confidence threshold as the blue dashed line. The absence of a significant peak in the periodogram places limits on short-period companions.
Figure 8: Limits on companion mass as a function of semi-major axis for GJ 896A. Limits are derived from acceleration and angular separation terms. The space below the solid curve is allowed.
JD Flux
(mJy) (J2000) (mas) (J2000) (mas)
2453820 23 34 14.1564960 0.135 20 10 28.882640 0.279
2454482 23 34 14.1564897 0.007 20 10 28.882867 0.013
2454585 23 34 14.1564998 0.013 20 10 28.882864 0.022
2454680 23 34 14.1564977 0.018 20 10 28.882824 0.033
2454793 23 34 14.1565001 0.015 20 10 28.882793 0.025
2454978 23 34 14.1564972 0.010 20 10 28.882838 0.017
2455101 23 34 14.1564926 0.007 20 10 28.883053 0.015
2455226 23 34 14.1565001 0.004 20 10 28.882733 0.006
2455467 23 34 14.1565155 0.006 20 10 28.882853 0.010
2455591 23 34 14.1565109 0.006 20 10 28.882888 0.010

Note. – The first line is the mean position reported in Paper I. Errors in the position are errors from fitting and do not include systematic effects.

Table 1: Positions for Secondary Calibrator J2334+2010
JD Flux
(Jy) (J2000) (mas) (J2000) (mas)
2453818.310 23 31 52.4336850 0.113 19 56 13.712351 0.171
2453820.310 23 31 52.4342630 0.183 19 56 13.714736 0.202
2453821.310 23 31 52.4345400 0.028 19 56 13.715419 0.045
2454482.443 23 31 52.4967356 0.031 19 56 13.570094 0.081
2454679.922 23 31 52.5343729 0.062 19 56 13.704882 0.128
2454792.651 23 31 52.5301598 0.015 19 56 13.570517 0.031
2454978.109 23 31 52.5703944 0.033 19 56 13.606525 0.060
2455100.776 23 31 52.5719845 0.018 19 56 13.599548 0.051
2455226.443 23 31 52.5809398 0.059 19 56 13.443975 0.101
2455466.776 23 31 52.6126099 0.021 19 56 13.537846 0.045
2455591.443 23 31 52.6215506 0.018 19 56 13.383628 0.044

Note. – The first three measurements are the results from Paper I.

Table 2: Positions for GJ 896A
Parameter Unit Optical RIPL Astrometric RIPL Astrometric + Accel
23 31 52.17898 23 31 52.179555 23 31 52.179565
(mas) 12.31 0.409 0.123
19 56 14.1505 19 56 14.141024 19 56 14.139842
(mas) 9.51 0.434 0.136
Table 3: Astrometric Parameters for GJ 896A
Period Semi-Major Axis Inclination Periastron Eccentricity
(y) (arcsec) (deg) (deg) (y) (deg)
359 6.87 123.5 82.1 2008 0.20 354.0
Table 4: Orbital Parameters for the GJ 896 AB System From Optical Astrometry
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description