We present the visual orbits of two long period spectroscopic binary stars, HD 8374 and HD 24546, using interferometric observations acquired with the CHARA Array and the Palomar Testbed Interferometer. We also obtained new radial velocities from echelle spectra using the APO 3.5 m and Fairborn 2.0 m telescopes. By combining the visual and spectroscopic observations, we solve for the full, three-dimensional orbits and determine the stellar masses and distances to within 3% uncertainty. We then estimate the effective temperature and radius of each component star through Doppler tomography and spectral energy distribution analyses, in order to compare the observed stellar parameters to the predictions of stellar evolution models. For HD 8374, we find masses of M1 = 1.636 +/- 0.050 Msun and M2 = 1.587 +/- 0.049 Msun, radii of R1 = 1.84 +/- 0.05 Rsun and R2 = 1.66 +/- 0.12 Rsun, temperatures of Teff1 = 7280 +/- 110 K and Teff2 = 7280 +/- 120 K, and an estimated age of 1.0 Gyr. For HD 24546, we find masses of M1 = 1.434 +/- 0.014 Msun and M2 = 1.409 +/- 0.014 Msun, radii of R1 = 1.67 +/- 0.06 Rsun and R2 = 1.60 +/- 0.10 Rsun, temperatures of Teff1 = 6790 +/- 120 K and Teff2 = 6770 +/- 90 K, and an estimated age of 1.4 Gyr. HD 24546 is therefore too old to be a member of the Hyades cluster, despite its physical proximity to the group.

We report 370 measures of 170 components of binary and multiple star systems, obtained from speckle imaging observations made with the Differential Speckle Survey Instrument at Lowell Observatory's Discovery Channel Telescope in 2015 through 2017. Of the systems studied, 147 are binary stars, 10 are seen as triple systems, and 1 quadruple system is measured. Seventy-six high-quality non-detections and fifteen newly resolved components are presented in our observations. The uncertainty in relative astrometry appears to be similar to our previous work at Lowell, namely linear measurement uncertainties of approximately 2 mas, and the relative photometry appears to be uncertain at the 0.1 to 0.15 magnitude level. Using these measures and those in the literature, we calculate six new visual orbits, including one for the Be star 66 Oph, and two combined spectroscopic-visual orbits. The latter two orbits, which are for HD 22451 (YSC 127) and HD 185501 (YSC 135), yield individual masses of the components at the level of 2 percent or better, and independent distance measures that in one case agrees with the value found in the Gaia DR2, and in the other disagrees at the 2-$\sigma$ level. We find that HD 22451 consists of an F6V+F7V pair with orbital period of $2401.1 \pm 3.2$ days and masses of $1.342 \pm 0.029$ and $1.236 \pm 0.026$ $ M_{\odot}$. For HD 185501, both stars are G5 dwarfs that orbit one another with a period of $433.94 \pm 0.15$ days, and the masses are $0.898 \pm 0.012$ and $0.876 \pm 0.012$ $ M_{\odot}$. We discuss the details of both the new discoveries and the orbit objects.

We report extensive differential V-band photometry and high-resolution spectroscopy for the 1.14 day, detached, double-lined eclipsing binary BT Vul (F0+F7). Our radial-velocity monitoring and light curve analysis lead to absolute masses and radii of M1 = 1.5439 +/- 0.0098 MSun and R1 = 1.536 +/- 0.018 RSun for the primary, and M2 = 1.2196 +/- 0.0080 MSun and R2 = 1.151 +/- 0.029 RSun for the secondary. The effective temperatures are 7270 +/- 150 K and 6260 +/- 180 K, respectively. Both stars are rapid rotators, and the orbit is circular. A comparison with stellar evolution models from the MIST series shows excellent agreement with these determinations, for a composition of [Fe/H] = +0.08 and an age of 350 Myr. The two components of BT Vul are very near the zero-age main sequence.

We report extensive differential V-band photometry and high-resolution spectroscopic observations of the early F-type, 1.06-day detached eclipsing binary V506 Oph. The observations along with times of minimum light from the literature are used to derive a very precise ephemeris and the physical properties for the components, with the absolute masses and radii being determined to 0.7% or better. The masses are 1.4153 +/- 0.0100 M(Sun) and 1.4023 +/- 0.0094 M(sun) for the primary and secondary, the radii are 1.725 +/- 0.010 R(Sun) and 1.692 +/- 0.012 R(Sun), and the effective temperatures 6840 +/- 150 K and 6780 +/- 110 K, respectively. The orbit is circular and the stars are rotating synchronously. The accuracy of the radii and temperatures is supported by the resulting distance estimate of 564 +/- 30 pc, in excellent agreement with the value implied by the trigonometric parallax listed in the Gaia/DR2 catalog. Current stellar evolution models from the MIST series for a composition of [Fe/H] = -0.04 match the properties of both stars in V506 Oph very well at an age of 1.83 Gyr, and indicate they are halfway through their core hydrogen-burning phase.

The Dharma Planet Survey (DPS) aims to monitor about 150 nearby very bright FGKM dwarfs (within 50 pc) during 2016$-$2020 for low-mass planet detection and characterization using the TOU very high resolution optical spectrograph (R$\approx$100,000, 380-900nm). TOU was initially mounted to the 2-m Automatic Spectroscopic Telescope at Fairborn Observatory in 2013-2015 to conduct a pilot survey, then moved to the dedicated 50-inch automatic telescope on Mt. Lemmon in 2016 to launch the survey. Here we report the first planet detection from DPS, a super-Earth candidate orbiting a bright K dwarf star, HD 26965. It is the second brightest star ($V=4.4$ mag) on the sky with a super-Earth candidate. The planet candidate has a mass of 8.47$\pm0.47M_{\rm Earth}$, period of $42.38\pm0.01$ d, and eccentricity of $0.04^{+0.05}_{-0.03}$. This RV signal was independently detected by Diaz et al. (2018), but they could not confirm if the signal is from a planet or from stellar activity. The orbital period of the planet is close to the rotation period of the star (39$-$44.5 d) measured from stellar activity indicators. Our high precision photometric campaign and line bisector analysis of this star do not find any significant variations at the orbital period. Stellar RV jitters modeled from star spots and convection inhibition are also not strong enough to explain the RV signal detected. After further comparing RV data from the star's active magnetic phase and quiet magnetic phase, we conclude that the RV signal is due to planetary-reflex motion and not stellar activity.

Because of the recent technological advances, the key technologies needed for precision space optical astrometry are now in hand. The Microarcsecond Astrometry Probe (MAP) mission concept is designed to find 1 Earth mass planets at 1AU orbit (scaled to solar luminosity) around the nearest ~90 FGK stars. The MAP payload includes i) a single three-mirror anastigmatic telescope with a 1-m primary mirror and metrology subsystems, and ii) a camera. The camera focal plane consists of 42 detectors, providing a Nyquist sampled FOV of 0.4-deg. Its metrology subsystems ensure that MAP can achieve the 0.8 uas astrometric precision in 1 hr, which is required to detect Earth-like exoplanets in our stellar neighborhood. MAP mission could provide ~10 specific targets for a much larger coronagraphic mission that would measure its spectra. We argue for the development of the space astrometric missions capable of finding Earth-2.0. Given the current technology readiness such missions relying on precision astrometry could be flown in the next decade, perhaps in collaboration with other national space agencies.

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process.

We report spectroscopic and differential photometric observations of the A-type system V482 Per that reveal it to be a rare hierarchical quadruple system containing two eclipsing binaries. One has the previously known orbital period of 2.4 days and a circular orbit, and the other a period of 6 days, a slightly eccentric orbit (e = 0.11), and shallow eclipses only 2.3% deep. The two binaries revolve around their common center of mass in a highly elongated orbit (e = 0.85) with a period of 16.67 yr. Radial velocities are measured for all components from our quadruple-lined spectra, and are combined with the light curves and with measurements of times of minimum light for the 2.4 day binary to solve for the elements of the inner and outer orbits simultaneously. The line-of-sight inclination angles of the three orbits are similar, suggesting they may be close to coplanar. The available observations appear to indicate that the 6 day binary experiences significant retrograde apsidal motion in the amount of about 60 degrees per century. We derive absolute masses for the four stars good to better than 1.5%, along with radii with formal errors of 1.1% and 3.5% for the 2.4 day binary and about 9% for the 6 day binary. A comparison of these and other physical properties with current stellar evolution models gives excellent agreement for a metallicity of [Fe/H] = -0.15 and an age of 360 Myr.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We report the detections of a giant planet (MARVELS-7b) and a brown dwarf candidate (MARVELS-7c) around the primary star in the close binary system, HD 87646. It is the first close binary system with more than one substellar circum-primary companion discovered to the best of our knowledge. The detection of this giant planet was accomplished using the first multi-object Doppler instrument (KeckET) at the Sloan Digital Sky Survey (SDSS) telescope. Subsequent radial velocity observations using ET at Kitt Peak National Observatory, HRS at HET, the "Classic" spectrograph at the Automatic Spectroscopic Telescope at Fairborn Observatory, and MARVELS from SDSS-III confirmed this giant planet discovery and revealed the existence of a long-period brown dwarf in this binary. HD 87646 is a close binary with a separation of $\sim22$ AU between the two stars, estimated using the Hipparcos catalogue and our newly acquired AO image from PALAO on the 200-inch Hale Telescope at Palomar. The primary star in the binary, HD 87646A, has Teff = 5770$\pm$80K, log(g)=4.1$\pm$0.1 and [Fe/H] = $-0.17\pm0.08$. The derived minimum masses of the two substellar companions of HD 87646A are 12.4$\pm$0.7M$_{\rm Jup}$ and 57.0$\pm3.7$M$_{\rm Jup}$. The periods are 13.481$\pm$0.001 days and 674$\pm$4 days and the measured eccentricities are 0.05$\pm$0.02 and 0.50$\pm$0.02 respectively. Our dynamical simulations show the system is stable if the binary orbit has a large semi-major axis and a low eccentricity, which can be verified with future astrometry observations.

Visible-light long baseline interferometry holds the promise of advancing a number of important applications in fundamental astronomy, including the direct measurement of the angular diameters and oblateness of stars, and the direct measurement of the orbits of binary and multiple star systems. To advance, the field of visible-light interferometry requires development of instruments capable of combining light from 15 baselines (6 telescopes) simultaneously. The Visible Imaging System for Interferometric Observations at NPOI (VISION) is a new visible light beam combiner for the Navy Precision Optical Interferometer (NPOI) that uses single-mode fibers to coherently combine light from up to six telescopes simultaneously with an image-plane combination scheme. It features a photometric camera for calibrations and spatial filtering from single-mode fibers with two Andor Ixon electron multiplying CCDs. This paper presents the VISION system, results of laboratory tests, and results of commissioning on-sky observations. A new set of corrections have been determined for the power spectrum and bispectrum by taking into account non-Gaussian statistics and read noise present in electron-multipying CCDs to enable measurement of visibilities and closure phases in the VISION post-processing pipeline. The post-processing pipeline has been verified via new on-sky observations of the O-type supergiant binary $\zeta$ Orionis A, obtaining a flux ratio of $2.18\pm0.13$ mag with a position angle of $223.9\pm1.0^{\circ}$ and separation $40.6\pm1.8$ mas over 570-750 nm, in good agreement with expectations from the previously published orbit.

We report extensive high-resolution spectroscopic observations and V-band differential photometry of the slightly eccentric 7.02-day detached eclipsing binary V501 Mon (A6m+F0), which we use to determine its absolute dimensions to high precision (0.3% for the masses and 1.8% for the radii, or better). The absolute masses, radii, and temperatures are M(A) = 1.6455 +/- 0.0043 M(Sun), R(A) = 1.888 +/- 0.029 R(Sun), and T(A) = 7510 +/- 100 K for the primary, and M(B) = 1.4588 +/- 0.0025 M(Sun), R(B) = 1.592 +/- 0.028 R(Sun), and T(B) = 7000 +/- 90 K for the secondary. Apsidal motion has been detected, to which General Relativity contributes approximately 70%. The primary star is found to be a metallic-line A star. A detailed chemical analysis of the disentangled spectra yields abundances for more than a dozen elements in each star. Based on the secondary, the system metallicity is near solar: [Fe/H] = +0.01 +/- 0.06. Lithium is detected in the secondary but not in the primary. A comparison with current stellar evolution models shows a good match to the measured properties at an age of about 1.1 Gyr.

To measure the stellar and orbital properties of the metal-poor RS CVn binary o Draconis (o Dra), we directly detect the companion using interferometric observations obtained with the Michigan InfraRed Combiner at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array. The H-band flux ratio between the primary and secondary stars is the highest confirmed flux ratio (370 +/- 40) observed with long-baseline optical interferometry. These detections are combined with radial velocity data of both the primary and secondary stars, including new data obtained with the Tillinghast Reflector Echelle Spectrograph on the Tillinghast Reflector at the Fred Lawrence Whipple Observatory and the 2-m Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We determine an orbit from which we find model-independent masses and ages of the components (M_A = 1.35 +\- 0.05 M_Sun, M_B = 0.99 +\- 0.02 M_Sun, system age = 3.0 -\+ 0.5 Gyr). An average of a 23-year light curve of o Dra from the Tennessee State University Automated Photometric Telescope folded over the orbital period newly reveals eclipses and the quasi-sinusoidal signature of ellipsoidal variations. The modeled light curve for our system's stellar and orbital parameters confirm these ellipsoidal variations due to the primary star partially filling its Roche lobe potential, suggesting most of the photometric variations are not due to stellar activity (starspots). Measuring gravity darkening from the average light curve gives a best-fit of beta = 0.07 +\- 0.03, a value consistent with conventional theory for convective envelope stars. The primary star also exhibits an anomalously short rotation period, which, when taken with other system parameters, suggests the star likely engulfed a low-mass companion that had recently spun-up the star.

The dwarf stars in the 26 year period binary alpha Com were predicted to eclipse each other in early 2015. That prediction was based on an orbit model made with over 600 astrometric observations using micrometers, speckle interferometry, and long baseline optical interferometry. Unfortunately, it has been realized recently that the position angle measurements for three of the observations from ~100 years ago were in error by 180 degrees, which skewed the orbital fit. The eclipse was likely 2 months earlier than predicted, at which point the system was low on the horizon at sunrise.

The dwarf stars in the 26 year period binary Alpha Com are predicted to eclipse each other in early 2015. Simulations of randomly selected orbital parameters based on the most recently published orbital elements and uncertainties show the eclipse will happen within 3 days of MJD 57047 (January 25, 2015; 1 sigma; all simulations fall within 2 weeks of this) and will have a depth of several tenths of a magnitude. The eclipse will last between 28 and 45 hours, depending on the exact orbital parameters. Observers are encouraged to monitor this unique event.

We acquired Johnson BV photometry of the binary Be disk system delta~Scorpii during its 2009, 2010, 2011, and 2012 observing seasons and used it to probe the innermost regions of the disk. We found that several disk building events have occurred during this time, resulting in an overall brightening in the V-band and reddening of the system. In addition to these long-term trends, we found cyclical variability in each observing season on timescales between 60 and 100 days. We were able to reproduce the changes in the magnitude and colour of delta Sco using our theoretical models and found that variable mass-loss rates in the range 2.5-7.0x10^{-9} M_{sun}/yr over ~35 days can reproduce the observed increase in brightness.

AP And is a well-detached F5 eclipsing binary star for which only a very limited amount of information was available before this publication. We have obtained very extensive measurements of the light curve (19097 differential V magnitude observations) and a radial velocity curve (83 spectroscopic observations) which allow us to fit orbits and determine the absolute properties of the components very accurately: masses of 1.277 +/- 0.004 and 1.251 +/- 0.004 solar masses, radii of 1.233 +/- 0.006 and 1.1953 +/- 0.005 solar radii, and temperatures of 6565 +/- 150 K and 6495 +/- 150 K. The distance to the system is about 400 +/- 30 pc. Comparison with the theoretical properties of the stellar evolutionary models of the Yonsei-Yale series of Yi et al. shows good agreement between the observations and the theory at an age of about 500 Myr and a slightly sub-solar metallicity.

We have determined spectroscopic orbits for five single-lined spectroscopic binaries, HD 100167, HD 135991, HD 140667, HD 158222, HD 217924. Their periods range from 60.6 to 2403 days and the eccentricities, from 0.20 to 0.84. Our spectral classes for the stars confirm that they are of solar type, F9 to G5, and all are dwarfs. Their [Fe/H] abundances, determined spectroscopically, are close to the solar value and on average are 0.12 greater than abundances from a photometric calibration. Four of the five stars are rotating faster than their predicted pseudosynchronous rotational velocities.

We used the Palomar Testbed Interferometer (PTI) to resolve 2.2 $\mu$m emission from the classical nova V458 Vul 2007 over the course of several days following its discovery on 2007 August 8.54 UT. We also obtained K-band photometric data and spectra of the nova during the early days of the outburst. We also used photometric measurements from the AAVSO database. This is a unique data set offering a 3-technique approach: high-resolution imaging, spectroscopy and photometry. Our analysis shows that the nova ejecta can be modeled as an inclined disk at low inclination i.e. low ellipticity which is consistent with the nova being in the fireball phase at which the outflowing gas is optically thick, confirmed by the presence of strong P-Cygni Balmer lines in the spectra. The expansion velocity is $\approx$1700 $\rm km\ s^{-1}$, derived from the H$\alpha$ line. By combining the nova's angular expansion rate measured by PTI with the expansion rate measured from spectroscopy, the inferred distance to the nova is 9.9-11.4 kpc. We also used the K-band fluxes and the derived size of the emission to estimate the total mass ejected from the nova $\approx 4\times 10^{-4} M_{\odot}$. The quick transition of the nova from Fe II to He/N class makes V458 Vul 2007 a hybrid nova.

We present the orbital and physical parameters of a newly discovered low-mass detached eclipsing binary from the All-Sky Automated Survey (ASAS) database: ASAS J011328-3821.1 A - a member of a visual binary system with the secondary component separated by about 1.4 seconds of arc. The radial velocities were calculated from the high-resolution spectra obtained with the 1.9-m Radcliffe/GIRAFFE, 3.9-m AAT/UCLES and 3.0-m Shane/HamSpec telescopes/spectrographs on the basis of the TODCOR technique and positions of H_alpha emission lines. For the analysis we used V and I band photometry obtained with the 1.0-m Elizabeth and robotic 0.41-m PROMPT telescopes, supplemented with the publicly available ASAS light curve of the system. We found that ASAS J011328-3821.1 A is composed of two late-type dwarfs having masses of M_1 = 0.612 +/- 0.030 M_sun, M_2 = 0.445 +/- 0.019 M_sun and radii of R_1 = 0.596 +/- 0.020 R_sun, R_2 = 0.445 +/- 0.024 R_sun, both show a substantial level of activity, which manifests in strong H_alpha and H_beta emission and the presence of cool spots. The influence of the third light on the eclipsing pair properties was also evaluated and the photometric properties of the component B were derived. Comparison with several popular stellar evolution models shows that the system is on its main sequence evolution stage and probably is more metal rich than the Sun. We also found several clues which suggest that the component B itself is a binary composed of two nearly identical ~0.5 M_sun stars.

We present the most precise to date orbital and physical parameters of the well known short period (P=5.975 d), eccentric (e=0.3) double-lined spectroscopic binary BY Draconis, a prototype of a class of late-type, active, spotted flare stars. We calculate the full spectroscopic/astrometric orbital solution by combining our precise radial velocities (RVs) and the archival astrometric measurements from the Palomar Testbed Interferometer (PTI). The RVs were derived based on the high resolution echelle spectra taken between 2004 and 2008 with the Keck I/HIRES, Shane/CAT/HamSpec and TNG/SARG telescopes/spectrographs using our novel iodine-cell technique for double-lined binary stars. The RVs and available PTI astrometric data spanning over 8 years allow us to reach 0.2-0.5% level of precision in M sin3(i) and the parallax but the geometry of the orbit (i=154 deg) hampers the absolute mass precision to 3.3%, which is still an order of magnitude better than for previous studies. We compare our results with a set of Yonsei-Yale theoretical stellar isochrones and conclude that BY Dra is probably a main sequence system more metal-rich than the Sun. Using the orbital inclination and the available rotational velocities of the components, we also conclude that the rotational axes of the components are likely misaligned with the orbital angular momentum. Given BY Dra's main sequence status, late spectral type and the relatively short orbital period, its high orbital eccentricity and probable spin-orbit misalignment are not in agreement with the tidal theory. This disagreement may possibly be explained by smaller rotational velocities of the components and the presence of a substellar mass companion to BY Dra AB.

Precise near-infrared radial velocimetry enables efficient detection and transit verification of low-mass extrasolar planets orbiting M dwarf hosts, which are faint for visible-wavelength radial velocity surveys. The TripleSpec Exoplanet Discovery Instrument, or TEDI, is the combination of a variable-delay Michelson interferometer and a medium-resolution (R=2700) near-infrared spectrograph on the Palomar 200" Hale Telescope. We used TEDI to monitor GJ 699, a nearby mid-M dwarf, over 11 nights spread across 3 months. Analysis of 106 independent observations reveals a root-mean-square precision of less than 37 m/s for 5 minutes of integration time. This performance is within a factor of 2 of our expected photon-limited precision. We further decompose the residuals into a 33 m/s white noise component, and a 15 m/s systematic noise component, which we identify as likely due to contamination by telluric absorption lines. With further development this technique holds promise for broad implementation on medium-resolution near-infrared spectrographs to search for low-mass exoplanets orbiting M dwarfs, and to verify low-mass transit candidates.

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.

Determining the type of matter that is inside a neutron star (NS) has been a long-standing goal of astrophysics. Despite this, most of the NS equations of state (EOS) that predict maximum masses in the range 1.4-2.8 solar masses are still viable. Most of the precise NS mass measurements that have been made to date show values close to 1.4 solar masses, but a reliable measurement of an over-massive NS would constrain the EOS possibilities. Here, we investigate how optical astrometry at the microarcsecond level can be used to map out the orbits of High-Mass X-ray Binaries (HMXBs), leading to tight constraints on NS masses. While previous studies by Unwin and co-workers and Tomsick and co-workers discuss the fact that the future Space Interferometry Mission should be capable of making such measurements, the current work describes detailed simulations for 6 HMXB systems, including predicted constraints on all orbital parameters. We find that the direct NS masses can be measured to an accuracy of 2.5% (1-sigma) in the best case (X Per), to 6.5% for Vela X-1, and to 10% for two other HMXBs.

The astrometry technique is an important tool for detecting and characterizing exoplanets of different type. In this review, the different projects which are either operating, in construction or in discussion are presented and their performance discussed in the framework of the Blue Dots study. We investigate the sensitivity of astrometry to different sources of noise and we show that astrometry is a key technique in the path of discovering and characterizing new types of planets including the very challenging category of Earth-like planets orbiting the habitable zone of solar-type stars.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

We present high precision radial velocities (RVs) of double-lined spectroscopic binary stars HD78418, HD123999, HD160922, HD200077 and HD210027. They were obtained based on the high resolution echelle spectra collected with the Keck I/Hires, Shane/CAT/Hamspec and TNG/Sarge telescopes/spectrographs over the years 2003-2008 as a part of TATOOINE search for circumbinary planets. The RVs were computed using our novel iodine cell technique for double-line binary stars. The precision of the RVs is of the order of 1-10 m/s. Our RVs combined with the archival visibility measurements from the Palomar Testbed Interferometer allow us to derive very precise spectroscopic/astrometric orbital and physical parameters of the binaries. In particular, we derive the masses, the absolute K and H band magnitudes and the parallaxes. The masses together with the absolute magnitudes in the K and H bands enable us to estimate the ages of the binaries. These RVs allow us to obtain some of the most accurate mass determinations of binary stars. The fractional accuracy in m*sin(i) only and hence based on the RVs alone ranges from 0.02% to 0.42%. When combined with the PTI astrometry, the fractional accuracy in the masses ranges in the three best cases from 0.06% to 0.5%. Among them, the masses of HD210027 components rival in precision the mass determination of the components of the relativistic double pulsar system PSRJ0737-3039. In the near future, for double-lined eclipsing binary stars we expect to derive masses with a fractional accuracy of the order of up to ~0.001% with our technique. This level of precision is an order of magnitude higher than of the most accurate mass determination for a body outside the Solar System - the double neutron star system PSRB1913+16.

We present preliminary results of the first and on-going radial velocity survey for circumbinary planets. With a novel radial velocity technique employing an iodine absorption cell we achieve an unprecedented RV precision of up to 2 m/s for double-lined binary stars. The high resolution spectra collected with the Keck I/Hires, TNG/Sarg and Shane/CAT/Hamspec telescopes/spectrographs over the years 2003-2008 allow us to derive RVs and compute planet detection limits for ten double-lined binary stars. For this initial sample of targets, we can rule out planets on dynamically stable orbits with masses as small as ~0.3 to 3 MJup for the orbital periods of up to ~5.3 years. Even though the presented sample of stars is too small to make any strong conclusions, it is clear that the search for circumbinary planets is now technique-wise possible and eventually will provide new constraints for the planet formation theories.

We derive the absolute physical and orbital parameters for a sample of 18 detached eclipsing binaries from the \emph{All Sky Automated Survey} (ASAS) database based on the available photometry and our own radial velocity measurements. The radial velocities (RVs) are computed using spectra we collected with the 3.9-m Anglo-Australian Telescope and its \emph{University College London Echelle Spectrograph} and the 1.9-m SAAO Radcliffe telescope and its \emph{Grating Instrument for Radiation Analysis with a Fibre Fed Echelle}. In order to obtain as precise RVs as possible, most of the systems were observed with an iodine cell available at the AAT/UCLES and/or analyzed using the two-dimensional cross-correlation technique (TODCOR). The RVs were measured with TODCOR using synthetic template spectra as references. However, for two objects we used our own approach to the tomographic disentangling of the binary spectra to provide observed template spectra for the RV measurements and to improve the RV precision even more. For one of these binaries, AI Phe, we were able to the obtain an orbital solution with an RV $rms$ of 62 and 24 m s$^{-1}$ for the primary and secondary respectively. For this system, the precision in $M \sin^3{i}$ is 0.08%. For the analysis, we used the photometry available in the ASAS database. We combined the RV and light curves using PHOEBE and JKTEBOP codes to obtain the absolute physical parameters of the systems. Having precise RVs we were able to reach $\sim$0.2 % precision (or better) in masses in several cases but in radii, due to the limited precision of the ASAS photometry, we were able to reach a precision of only 1% in one case and 3-5 % in a few more cases. For the majority of our objects, the orbital and physical analysis is presented for the first time.

Astrometry can detect rocky planets in a broad range of masses and orbital distances and measure their masses and three-dimensional orbital parameters, including eccentricity and inclination, to provide the properties of terrestrial planets. The masses of both the new planets and the known gas giants can be measured unambiguously, allowing a direct calculation of the gravitational interactions, both past and future. Such dynamical interactions inform theories of the formation and evolution of planetary systems, including Earth-like planets. Astrometry is the only technique technologically ready to detect planets of Earth mass in the habitable zone (HZ) around solar-type stars within 20 pc. These Earth analogs are close enough for follow-up observations to characterize the planets by infrared imaging and spectroscopy with planned future missions such as the James Webb Space Telescope (JWST) and the Terrestrial Planet Finder/Darwin. Employing a demonstrated astrometric precision of 1 microarcsecond and a noise floor under 0.1 micro-arcseconds, SIM Lite can make multiple astrometric measurements of the nearest 60 F-, G-, and K-type stars during a five-year mission. SIM Lite directly tests theories of rocky planet formation and evolution around Sun-like stars and identifies the nearest potentially habitable planets for later spaceborne imaging, e.g., with Terrestrial Planet Finder and Darwin. SIM was endorsed by the two recent Decadal Surveys and it meets the highest-priority goal of the 2008 AAAC Exoplanet Task Force.

We ask if Earth-like planets (terrestrial mass and habitable-zone orbit) can be detected in multi-planet systems, using astrometric and radial velocity observations. We report here the preliminary results of double-blind calculations designed to answer this question.

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 TEDI (TripleSpec - Exoplanet Discovery Instrument) will be the first instrument fielded specifically for finding low-mass stellar companions. The instrument is a near infra-red interferometric spectrometer used as a radial velocimeter. TEDI joins Externally Dispersed Interferometery (EDI) with an efficient, medium-resolution, near IR (0.9 - 2.4 micron) echelle spectrometer, TripleSpec, at the Palomar 200" telescope. We describe the instrument and its radial velocimetry demonstration program to observe cool stars.

Externally Dispersed Interferometry (EDI) is the series combination of a fixed-delay field-widened Michelson interferometer with a dispersive spectrograph. This combination boosts the spectrograph performance for both Doppler velocimetry and high resolution spectroscopy. The interferometer creates a periodic spectral comb that multiplies against the input spectrum to create moire fringes, which are recorded in combination with the regular spectrum. The moire pattern shifts in phase in response to a Doppler shift. Moire patterns are broader than the underlying spectral features and more easily survive spectrograph blurring and common distortions. Thus, the EDI technique allows lower resolution spectrographs having relaxed optical tolerances (and therefore higher throughput) to return high precision velocity measurements, which otherwise would be imprecise for the spectrograph alone.

We have resolved the classical nova V1663 Aql using long-baseline near-IR interferometry covering the period from 5--18 days after peak brightness. We directly measure the shape and size of the fireball, which we find to be asymmetric. In addition we measure an apparent expansion rate of 0.21 +/- 0.03 mas/day. Assuming a linear expansion model we infer a time of initial outburst approximately 4 days prior to peak brightness. When combined with published spectroscopic expansion velocities our angular expansion rate implies a distance of 8.9 +/- 3.6 kpc. This distance measurement is independent of, but consistent with, determinations made using widely available photometric relations for novae.

Searches for planets in close binary systems explore the degree to which stellar multiplicity inhibits or promotes planet formation. There is a degeneracy between planet formation models when only systems with single stars are studied--several mechanisms appear to be able to produce such a final result. This degeneracy is lifted by searching for planets in binary systems; the resulting detections (or evidence of non-existence) of planets in binaries isolates which models may contribute to how planets form in nature. In this chapter, we consider observational efforts to detect planetary companions to binary stars in two types of hierarchical planet-binary configurations: first ``S-type'' planets which orbit just one of the stars, with the binary period being much longer than the planet's; second, ``P-type'' or circumbinary planets, where the planet simultaneously orbits both stars, and the planetary orbital period is much longer than that of the binary. The S-type planet finding techniques are different for binaries that can or cannot be spatially resolved. For wider systems, techniques reviewed include dualstar interferometric differential astrometry and precision radial velocities. Alternatively, unresolved binaries can be studied using modified dualstar "PHASES-style" differential astrometry or a modification of the radial velocity technique for composite spectra. Should a fortunately aligned--but still long period--binary be found, eclipse timing can also reveal the presence of S-type planets. Methods for detecting P-type planets include the composite-spectra variant of the radial velocity technique and eclipse timing.

We report observations of the recurrent nova RS Oph using long-baseline near-IR interferometry. We are able to resolve emission from the nova for several weeks after the February 2006 outburst. The near-IR source initially expands to a size of approximately 5 milli-arcseconds. However, beginning around day 10 the IR source appears to begin to shrink, reaching approximately 2 milli-arcseconds by day 100. We combine our measured angular diameters with previously available interferometric and photometric data to derive an emission measure for the source, and hence are able to determine the mass-loss rate of the nova in the days following the outburst.

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.

Hostile tidal forces may inhibit the formation of Jovian planets in binaries with semimajor axes of $\la$$50\au$, binaries that might be called ``close'' in this context. As an alternative to in situ planet formation, a binary can acquire a giant planet when one of its original members is replaced in a dynamical interaction with another star that hosts a planet. Simple scaling relations for the structure and evolution of star clusters, coupled with analytic arguments regarding binary-single and binary-binary scattering, indicate that dynamical processes can deposit Jovian planets in $<$1% of close binaries. If ongoing and future exoplanet surveys measure a much larger fraction, it may be that giant planets do somehow form frequently in such systems.

Following the recent outburst of the recurrent nova RS Oph on 2006 Feb 12, we measured its near-infrared size using the IOTA, Keck, and PTI Interferometers at multiple epochs. The characteristic size of ~3 milliarcseconds hardly changed over the first 60 days of the outburst, ruling out currently-popular models whereby the near-infrared emission arises from hot gas in the expanding shock. The emission was also found to be significantly asymmetric, evidenced by non-zero closure phases detected by IOTA. The physical interpretation of these data depend strongly on the adopted distance to RS Oph. Our data can be interpreted as the first direct detection of the underlying RS Oph binary, lending support to the recent ``reborn red giant'' models of Hachisu & Kato. However, this result hinges on an RS Oph distance of ~< 540 pc, in strong disagreement with the widely-adopted distance of ~1.6 kpc. At the farther distance, our observations imply instead the existence of a non-expanding, dense and ionized circumbinary gaseous disk or reservoir responsible for the bulk of the near-infrared emission. Longer-baseline infrared interferometry is uniquely suited to distinguish between these models and to ultimately determine the distance, binary orbit, and component masses for RS Oph, one of the closest-known (candidate) SNIa progenitor systems.

During the last two years we have used the Palomar Testbed Interferometer to observe several explosive variable stars, including V838 Monocerotis, V1663 Aquilae and recently RS Ophiuchi. We observed V838 Monocerotis approximately 34 months after its eruption, and were able to resolve the ejecta. Observations of V1663 Aql were obtained starting 9 days after peak brightness and continued for 10 days. We were able to resolve the milliarcsecond-scale emission and follow the expansion of the nova photosphere. When combined with radial-velocity information, these observations can be used to infer the distance to the nova. Finally we have resolved the recurrent nova RS Oph and can draw some preliminary conclusions regarding the emission morphology.

One of the biggest challenges associated with a nulling interferometer-based approach to detecting extra-solar Earth-like planets comes from the extremely stringent requirements of pathlength, polarization and amplitude matching in the interferometer. To the extent that the light from multiple apertures are not matched in these properties, light will leak through the nuller and confuse the search for a planetary signal. Here we explore the possibility of using the coherence properties of the starlight to separate contributions from the planet and nuller leakage. We find that straightforward modifications to the optical layout of a nulling interferometer will allow one to measure and correct for the leakage to a high degree of precision. This nulling calibration relaxes the field matching requirements substantially, and should consequently simplify the instrument design.

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.

Atmospheric turbulence is a serious problem for ground-based interferometers. It places tight limits on both sensitivity and measurement precision. Phase referencing is a method to overcome these limitations via the use of a bright reference star. The Palomar Testbed Interferometer was designed to use phase referencing and so can provide a pair of phase-stabilized starlight beams to a second (science) beam combiner. We have used this capability for several interesting studies, including very narrow angle astrometry. For close (1-arcsecond) pairs of stars we are able to achieve a differential astrometric precision in the range 20-30 micro-arcseconds.

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.

We have used the Palomar Testbed Interferometer to perform very high precision differential astrometry on the 0.25 arcsecond separation binary star HD 171779. In 70 minutes of observation we achieve a measurement uncertainty of approximately 9 micro-arcseconds in one axis, consistent with theoretical expectations. Night-to-night repeatability over four nights is at the level of 16 micro-arcseconds. This method of very-narrow-angle astrometry may be extremely useful for searching for planets with masses as small as 0.5 Jupiter Masses around a previously neglected class of stars -- so-called ``speckle binaries.'' It will also provide measurements of stellar parameters such as masses and distances, useful for constraining stellar models at the 10^-3 level.