The Palomar High-precision Astrometric Search for Exoplanet Systems monitored 51 subarcsecond binary systems to evaluate whether tertiary companions as small as Jovian planets orbited either the primary or secondary stars, perturbing their otherwise smooth Keplerian motions. Six binaries are presented that show evidence of substellar companions orbiting either the primary or secondary star. Of these six systems, the likelihoods of two of the detected perturbations to represent real objects are considered to be "high confidence", while the remaining four systems are less certain and will require continued observations for confirmation.

The Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) monitored 51 sub-arcsecond binary systems to determine precision binary orbits, study the geometries of triple and quadruple star systems, and discover previously unknown faint astrometric companions as small as giant planets. PHASES measurements made with the Palomar Testbed Interferometer (PTI) from 2002 until PTI ceased normal operations in late 2008 are presented. Infrared differential photometry of several PHASES targets were measured with Keck Adaptive Optics and are presented.

Differential astrometry measurements from the Palomar High-precision Astrometric Search for Exoplanet Systems have been combined with lower precision single-aperture measurements covering a much longer timespan (from eyepiece measurements, speckle interferometry, and adaptive optics) to determine improved visual orbits for 20 binary stars. In some cases, radial velocity observations exist to constrain the full three-dimensional orbit and determine component masses. The visual orbit of one of these binaries---alpha Com (HD 114378)---shows that the system is likely to have eclipses, despite its very long period of 26 years. The next eclipse is predicted to be within a week of 2015 January 24.

The Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) monitored 51 subarcsecond binary systems to evaluate whether tertiary companions as small as Jovian planets orbited either the primary or secondary stars, perturbing their otherwise smooth Keplerian motions. Twenty-one of those systems were observed 10 or more times and show no evidence of additional companions. A new algorithm is presented for identifying astrometric companions and establishing the (companion mass)-(orbital period) combinations that can be excluded from existence with high confidence based on the PHASES observations, and the regions of mass-period phase space being excluded are presented for 21 PHASES binaries.

(Abridged) Differential astrometry measurements from the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) are used to constrain the astrometric orbit of the previously known lesssim 2 day subsystem in the triple system 63 Gem A and have detected a previously unknown 2 year Keplerian wobble superimposed on the visual orbit of the much longer period (213 years) binary system HR 2896. The very small astrometric perturbation caused by the inner pair in 63 Gem A stretches the limits of current astrometric capabilities, but PHASES observations are able to constrain the orientation of the orbit. The two bright stars comprising the HR 2896 long period (213 year) system have a combined spectral type of K0III and the newly detected object's mass estimate places it in the regime of being a M dwarf. The motion of the stars are slow enough that their spectral features are always blended, preventing Doppler studies. The PHASES measurements and radial velocities (when available) have been combined with lower precision single-aperture measurements covering a much longer timeframe to improve the characterization of the long period orbits in both binaries. The visual orbits of the short and long period systems are presented for both systems, and used to calculate two possible values of the mutual inclinations between inner and outer orbits of 152 pm 12 degrees or a less likely value of 31 pm 11 degrees for 63 Gem A and 10.2 pm 2.4 degrees or 171.2 pm 2.8 degrees for HR 2896. The first is not coplanar, whereas the second is either nearly coplanar or anti-coplanar.

We have used high precision differential astrometry from the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) project and radial velocity measurements covering a time-span of 20 years to determine the orbital parameters of the 88 Tau A system. 88 Tau is a complex hierarchical multiple system comprising a total of six stars; we have studied the brightest 4, consisting of two short-period pairs orbiting each other with an 18-year period. We present the first orbital solution for one of the short-period pairs, and determine the masses of the components and distance to the system to the level of a few percent. In addition, our astrometric measurements allow us to make the first determination of the mutual inclinations of the orbits. We find that the sub-systems are not coplanar.

mu Orionis was identified by spectroscopic studies as a quadruple star system. Seventeen high precision differential astrometry measurements of mu Ori have been collected by the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES). These show both the motion of the long period binary orbit and short period perturbations superimposed on that caused by each of the components in the long period system being themselves binaries. The new measurements enable the orientations of the long period binary and short period subsystems to be determined. Recent theoretical work predicts the distribution of relative inclinations between inner and outer orbits of hierarchical systems to peak near 40 and 140 degrees. The degree of coplanarity of this complex system is determined, and the angle between the planes of the A-B and Aa-Ab orbits is found to be 136.7 +/- 8.3 degrees, near the predicted distribution peak at 140 degrees; this result is discussed in the context of the handful of systems with established mutual inclinations. The system distance and masses for each component are obtained from a combined fit of the PHASES astrometry and archival radial velocity observations. The component masses have relative precisions of 5% (component Aa), 15% (Ab), and 1.4% (each of Ba and Bb). The median size of the minor axes of the uncertainty ellipses for the new measurements is 20 micro-arcseconds. Updated orbits for delta Equulei, kappa Pegasi, and V819 Herculis are also presented.

The Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) has monitored 37 sub-arcsecond binary systems to determine whether their Keplerian orbits are perturbed by faint astrometric companions to either star. Software has been developed to evaluate the regions in a companion mass-period phase space in which the PHASES observations can exclude the possibility of face-on orbit perturbations. We present results for 8 systems for which astrometric companions with masses as small as those of giant planets can be excluded.

A new observing mode for the Palomar Testbed Interferometer was developed in2002-2003 which enables differential astrometry at the level of 20 micro-arcseconds for binary systems with separations of several hundred milli-arcseconds (mas). This phase-referenced mode is the basis of the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES), a search for giant planets orbiting either the primary or secondary star in fifty binary systems. We present the first science results from the PHASES search. The properties of the stars comprising binary systems are determined to high precision. The mutual inclinations of several hierarchical triple star systems have been determined. We will present upper limits constraining the the existence of giant planets in a few of the target systems.

V819 Herculis is a well-studied triple star system consisting of a ``wide'' pair with 5.5 year period, one component of which is a 2.2-day period eclipsing single-line spectroscopic binary. Differential astrometry measurements from the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) are presented and used to determine a relative inclination between the short- and long-period orbits of 23.6 +- 4.9 degrees. This represents only the sixth unambiguous determination of the mutual inclination of orbits in a hierarchical triple system. This result is combined with those for the other five systems for analysis of the observed distribution of mutual inclinations in nearby triple systems. It is found that this distribution is different than that which one would expect from random orientations with statistical significance at the 94% level; implications for studying the spatial distribution of angular momentum in star forming regions is discussed.

kappa Pegasi is a well-known, nearby triple star system. It consists of a ``wide'' pair with semi-major axis 235 milli-arcseconds, one component of which is a single-line spectroscopic binary (semi-major axis 2.5 milli-arcseconds). Using high-precision differential astrometry and radial velocity observations, the masses for all three components are determined and the relative inclinations between the wide and narrow pairs' orbits is found to be 43.8 +/- 3.0 degrees, just over the threshold for the three body Kozai resonance. The system distance is determined to 34.60 +/- 0.21 parsec, and is consistent with trigonometric parallax measurements.

delta Equulei is among the most well-studied nearby binary star systems. Results of its observation have been applied to a wide range of fundamental studies of binary systems and stellar astrophysics. It is widely used to calibrate and constrain theoretical models of the physics of stars. We report 27 high precision differential astrometry measurements of delta Equulei from the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES). The median size of the minor axes of the uncertainty ellipses for these measurements is 26 micro-arcseconds. These data are combined with previously published radial velocity data and other previously published differential astrometry measurements using other techniques to produce a combined model for the system orbit. The distance to the system is determined to within a twentieth of a parsec and the component masses are determined at the level of a percent. The constraints on masses and distance are limited by the precisions of the radial velocity data; we outline plans improve this deficiency and discuss the outlook for further study of this binary.

The Palomar Testbed Interferometer has observed several binary star systems whose separations fall between the interferometric coherence length (a few hundredths of an arcsecond) and the typical atmospheric seeing limit of one arcsecond. Using phase-referencing techniques we measure the relative separations of the systems to precisions of a few tens of micro-arcseconds. We present the first scientific results of these observations, including the astrometric detection of the faint third stellar component of the kappa Pegasi system.

The structure of the radio source RC B1146+052 of the ``Cold'' catalogue is investigated by data of the MIT-GB-VLA survey at 4850 MHz. This source belongs to the steep spectrum radio sources subsample of the RC catalogue. Its spectral index is $\alpha$ = -1.04. The optical image of this source obtained with 6m telescope is analysed. The radio source center is situated in a group of 8 galaxies of about 24$^m$ in the R-filter. The possible explanations of the complex structure of radio components are considered.

The gravitational lens 0957+561 was monitored with the Very Large Array from 1979 to 1997. The 6 cm light curve data from 1995-1997 and the 4 cm data from 1990-1997 are reported here. At 4 cm, the intrinsic source variations occur earlier and are twice as large as the corresponding variations at 6 cm. The VLBI core and jet components have different magnification factors, leading to different flux ratios for the varying and non-varying portions of the VLA light curves. Using both the PRHQ and Dispersion statistical techniques, we determined the time delay, core flux ratio, and excess non-varying B image flux density. The fits were performed for the 4 cm and 6 cm light curves, both individually and jointly, and we used Gaussian Monte Carlo data to estimate 68% statistical confidence levels. The delay estimates from each individual wavelength were inconsistent given the formal uncertainties, suggesting that there are unmodeled systematic errors in the analysis. We roughly estimate the systematic uncertainty in the joint result from the difference between the 6 cm and 4 cm results, giving 409+-30 days for the PRHQ statistic and 397+-20 days for the Dispersion statistic. These results are consistent with the current optical time delay of 417+-3 days, reconciling the long-standing difference between the optical and radio light curves and between different statistical analyses. The unmodeled systematic effects may also corrupt light curves for other lenses, and we caution that multiple events at multiple wavelengths may be necessary to determine an accurate delay in any lens system. Now that consensus has been reached regarding the time delay in the 0957+561 system, the most pressing issue remaining for determining H_0 is a full understanding of the mass distribution in the lens.

We report the discovery of a binary quasar, MGC2214+3550 A,B, whose components have similar optical spectra at a redshift z=0.88. The quasars are separated on the sky by 3.0", and have a magnitude difference of Delta(m_I)=0.5 mag. The VLA radio map at 3.6 cm shows a single 47mJy radio source with a core-jet morphology that is coincident with the brighter optical quasar A. Gravitational lensing is ruled out by the lack of radio emission from quasar B, and the lack of any visible galaxies to act as the lens. We conclude that MGC 2214+3550 A and B are physically associated. With a projected separation of 12.7h^{-1} kpc (Omega_0=1), MGC 2214+3550 A,B is one of the smallest z>0.5 binary quasars.

We present radio and optical observations of MG0248+0641, which contains a kiloparsec-scale ring-like structure in one of its radio lobes. The radio observations show a typical core-double morphology: a central core between two lobes, each of which has a hotspot. The western radio lobe appears as a nearly continuous ring, with linear polarization electric field vectors which are oriented in a radial direction from the ring center. We consider several different interpretations for the nature of this ring, including gravitational lensing of a normal jet by a foreground galaxy. Even though simple lensing models can describe the ring morphology reasonably well, the high linear polarization seen around the ring cannot be easily explained. The chance interposition of a galactic supernova remnant, nova, planetary nebula, or H II region, has been ruled out. The highly polarized ring of MG0248+0641 is much like the prominent ring seen in 3C219, and the multiple ones in 3C310 and Hercules A, suggesting that similar physical processes are producing shell structures in these radio galaxies. The ring in MG0248+0641 may be caused by the formation of ``bubbles'', as a result of instabilities in the energy flow down the western radio jet. It may also be possible that the required instabilities are triggered by the infall of gas, via tidal interaction of the central source with a nearby galaxy. This scenario may be indicated by our marginal detection of an optical source close to the western hotspot.

We have produced deep radio maps of the double quasar 0957+561 from multiple-epoch VLA observations. To achieve high sensitivity to extended structure we have re-reduced the best available 1.6 GHz observations and have combined 5 GHz data from multiple array configurations. Regions of faint emission approximately 15 arcsec north and south of the radio source G are probably lobes associated with the lensing galaxy. An arc 5 arcsec to the east of G may be a stretched image of emission in the background quasar's environment. 1.4 arcsec southwest of G we detect a source that we interpret as an image of emission from the quasar's western lobe, which could provide a constraint on the slope of the gravitational potential in the central region of the lens. We explore the consequences of these new constraints with simple lens models of the system.

We have designed and built a new type of CCD that we call an orthogonal transfer CCD (OTCCD), which permits parallel clocking horizontally as well as vertically. The device has been used successfully to remove image motion caused by atmospheric turbulence at rates up to 100 Hz, and promises to be a better, cheaper way to carry out image motion correction for imaging than by using fast tip/tilt mirrors. We report on the device characteristics, and find that the large number of transfers needed to track image motion does not significantly degrade the image either because of charge transfer inefficiency or because of charge traps. For example, after 100 sec of tracking at 100 Hz approximately 3% of the charge would diffuse into a skirt around the point spread function. Four nights of data at the Michigan-Dartmouth-MIT (MDM) 2.4-m telescope also indicate that the atmosphere is surprisingly benign, in terms of both the speed and coherence angle of image motion. Image motion compensation improved image sharpness by about 0.5 arcsec in quadrature with no degradation over a field of at least 3 arcminutes.

We report the discovery of a new gravitationally lensed radio source. Radio maps of MG0751+2716 show four lensed images, which, at higher resolution, are resolved into long arcs of emission. A group of galaxies is present in optical images, including the principal lensing galaxy, with a much brighter galaxy just a few arcseconds away. We have measured the redshift of this brighter galaxy. No optical counterpart to the background source has been detected. Lens models that can readily reproduce the lensed image positions all require a substantial shear component. However, neither the very elongated lens nor the bright nearby galaxy are correctly positioned to explain the shear. Lens models which associate the mass with the light of galaxies in the group can produce an acceptable fit, but only with an extreme mass-to-light ratio in one of the minor group members.

We report on 15 years of VLA monitoring of the gravitational lens B0957+561 at 6 cm. Since our last report in 1992, there have been 32 additional observations, in which both images have returned to their quiescent flux density levels and the A image has brightened again. We estimate the time delay from the light curves using three different techniques: the chi-squared analysis of Press, Rybicki, & Hewitt (1992a,b), the dispersion analysis of Pelt et al. (1994, 1996), and the locally normalized discrete correlation function of Leh\'ar et al. (1992). Confidence intervals for these time delay estimates are found using Monte Carlo techniques. With the addition of the new observations, it has become obvious that five observations from Spring 1990 are not consistent with the statistical properties of the rest of the light curves, so we analyze the light curves with those points removed, as well as the complete light curves. The three statistical techniques applied to the two versions of the data set result in time delay values in the range 398 to 461 days (or 1.09 to 1.26 years, A leading B), each with ~5% formal uncertainty. The corresponding flux ratios (B/A) are in the range 0.698 to 0.704. Thus, the new features in the light curve show that the time delay is less than 500 days, in contrast with analysis of earlier versions of the radio light curves. The large range in the time delay estimates is primarily due to unfortunate coincidences of observing gaps with flux variations.