M. M. Colavita

EXOZODIACAL DUST LEVELS FOR NEARBY MAIN-SEQUENCE STARS: A SURVEY WITH THE KECK INTERFEROMETER NULLER

The Keck Interferometer Nuller (KIN) was used to survey 25 nearby main-sequence stars in the mid-infrared, in order to assess the prevalence of warm circumstellar (exozodiacal) dust around nearby solar-type stars. The KIN measures circumstellar emission by spatially blocking the star but transmitting the circumstellar flux in a region typically 0.1-4 AU from the star. We find one significant detection (η Crv), two marginal detections (γ Oph and α Aql), and 22 clear non-detections. Using a model of our own solar systems zodiacal cloud, scaled to the luminosity of each target star, we estimate the equivalent number of target zodis needed to match our observations. Our three zodi detections are η Crv (1250 ± 260), γ Oph (200 ± 80), and α Aql (600 ± 200), where the uncertainties are 1σ. The 22 non-detected targets have an ensemble weighted average consistent with zero, with an average individual uncertainty of 160 zodis (1σ). These measurements represent the best limits to date on exozodi levels for a sample of nearby main-sequence stars. A statistical analysis of the population of 23 stars not previously known to contain circumstellar dust (excluding η Crv and γ Oph) suggests that, if the measurement errors are uncorrelated (for which we provide evidence) and if these 23 stars are representative of a single class with respect to the level of exozodi brightness, the mean exozodi level for the class is <150 zodis (3σ upper limit, corresponding to 99% confidence under the additional assumption that the measurement errors are Gaussian). We also demonstrate that this conclusion is largely independent of the shape and mean level of the (unknown) true underlying exozodi distribution.

The Visual Orbit of Pegasi

We have determined the visual orbit for the spectroscopic binary ι Pegasi with interferometric visibility data obtained by the Palomar Testbed Interferometer in 1997. ι Peg is a double-lined binary system whose minimum masses and spectral typing suggests the possibility of eclipses. Our orbital and component diameter determinations do not favor the eclipse hypothesis: the limb-to-limb separation of the two components is 0.151±0.069 mas at conjunction. Our conclusion that the ι Peg system does not eclipse is supported by high-precision photometric observations. The physical parameters implied by our visual orbit and the spectroscopic orbit of Fekel & Tomkin are in good agreement with those inferred by other means. In particular, the orbital parallax of the system is determined to be 86.9±1.0 mas, and masses of the two components are determined to be 1.326±0.016 and 0.819±0.009 M, respectively.

CONSTRAINING THE EXOZODIACAL LUMINOSITY FUNCTION OF MAIN-SEQUENCE STARS: COMPLETE RESULTS FROM THE KECK NULLER MID-INFRARED SURVEYS

Forty-seven nearby main-sequence stars were surveyed with the Keck Interferometer mid-infrared Nulling instrument (KIN) between 2008 and 2011, searching for faint resolved emission from exozodiacal dust. Observations of a subset of the sample have already been reported, focusing essentially on stars with no previously known dust. Here we extend this previous analysis to the whole KIN sample, including 22 more stars with known near-and/or far-infrared excesses. In addition to an analysis similar to that of the first paper of this series, which was restricted to the 8-9 µm spectral region, we present measurements obtained in all 10 spectral channels covering the 8-13 µm instrumental bandwidth. Based on the 8-9 µm data alone, which provide the highest signal-to-noise measurements, only one star shows a large excess imputable to dust emission (η Crv), while four more show a significant (> 3σ) excess: β Leo, β UMa, ζ Lep, and y Oph. Overall, excesses detected by KIN are more frequent around A-type stars than later spectral types. Astatistical analysis of the measurements further indicates that stars with known far-infrared (y ≥ 70 µm) excesses have higher exozodiacal emission levels than stars with no previous indication of a cold outer disk. This statistical trend is observed regardless of spectral type and points to a dynamical connection between the inner (zodi-like) and outer (Kuiper-Belt-like) dust populations. The measured levels for such stars are clustering close to the KIN detection limit of a few hundred zodis and are indeed consistent with those expected from a population of dust that migrated in from the outer belt by Poynting-Robertson drag. Conversely, no significant mid-ilinfrared excess is found around sources with previously reported near-infrared resolved excesses, which typically have levels of the order of 1% over the photospheric flux. If dust emission is really at play in these near-infrared detections, the absence of a strong mid-infrared counterpart points to populations of very hot and small (submicron) grains piling up very close to the sublimation radius. For solar-type stars with no known infrared excess, likely to be the most relevant targets for a future exo-Earth direct imaging mission, we find that their median zodi level is 12±24 zodis and lower than 60 (90) zodis with 95% (99%) confidence, if a lognormal zodi luminosity distribution is assumed.

Keck Interferometer Nuller Data Reduction and On-Sky Performance

We describe the Keck Interferometer nuller theory of operation, data reduction, and on-sky performance, particularly as it applies to the nuller exozodiacal dust key science program that was carried out between 2008 February and 2009 January. We review the nuller implementation, including the detailed phasor processing involved in implementing the null-peak mode used for science data and the sequencing used for science observing. We then describe the Level 1 reduction to convert the instrument telemetry streams to raw null leakages, and the Level 2 reduction to provide calibrated null leakages. The Level 1 reduction uses conservative, primarily linear processing, implemented consistently for science and calibrator stars. The Level 2 processing is more flexible, and uses diameters for the calibrator stars measured contemporaneously with the interferometer’s K-band cophasing system in order to provide the requisite accuracy. Using the key science data set of 462 total scans, we assess the instrument performance for sensitivity and systematic error. At 2.0 Jy we achieve a photometrically-limited null leakage uncertainty of 0.25% rms per 10 minutes of integration time in our broadband channel. From analysis of the Level 2 reductions, we estimate a systematic noise floor for bright stars of ~0.2% rms null leakage uncertainty per observing cluster in the broadband channel. A similar analysis is performed for the narrowband channels. We also provide additional information needed for science reduction, including details on the instrument beam pattern and the basic astrophysical response of the system, and references to the data reduction and modeling tools.

AN INTERFEROMETRIC STUDY OF THE FOMALHAUT INNER DEBRIS DISK. II. KECK NULLER MID-INFRARED OBSERVATIONS

We report on high-contrast mid-infrared observations of Fomalhaut obtained with the Keck Interferometer Nuller (KIN) showing a small resolved excess over the level expected from the stellar photosphere. The measured null excess has a mean value of 0.35% ± 0.10% between 8 and 11 μm and increases from 8 to 13 μm. Given the small field of view of the instrument, the source of this marginal excess must be contained within 2 AU of Fomalhaut. This result is reminiscent of previous VLTI K-band (≃2μm) observations, which implied the presence of a ~0.88% excess, and argued that thermal emission from hot dusty grains located within 6 AU from Fomalhaut was the most plausible explanation. Using a parametric two-dimensional radiative transfer code and a Bayesian analysis, we examine different dust disk structures to reproduce both the near- and mid-infrared data simultaneously. While not a definitive explanation of the hot excess of Fomalhaut, our model suggests that the most likely inner few AU disk geometry consists of a two-component structure, with two different and spatially distinct grain populations. The 2-11 μm data are consistent with an inner hot ring of very small (≃10-300 nm) carbon-rich grains concentrating around 0.1 AU. The second dust population—inferred from the KIN data at longer mid-infrared wavelengths—consists of larger grains (size of a few microns to a few tens of microns) located further out in a colder region where regular astronomical silicates could survive, with an inner edge around 0.4 AU-1 AU. From a dynamical point of view, the presence of the inner concentration of submicron-sized grains is surprising, as such grains should be expelled from the inner planetary system by radiation pressure within only a few years. This could either point to some inordinate replenishment rates (e.g., many grazing comets coming from an outer reservoir) or to the existence of some braking mechanism preventing the grains from moving out.

MASSES, LUMINOSITIES, AND ORBITAL COPLANARITIES OF THE μ ORIONIS QUADRUPLE-STAR SYSTEM FROM PHASES DIFFERENTIAL ASTROMETRY

μ Orionis was identified by spectroscopic studies as a quadruple-star system. Seventeen high-precision differential astrometry measurements of μ 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 (μas). Updated orbits for δ Equulei, κ Pegasi, and V819 Herculis are also presented.

PHASES differential astrometry and the mutual inclination of the V819 Herculis triple star system

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.

THE PHASES DIFFERENTIAL ASTROMETRY DATA ARCHIVE. II. UPDATED BINARY STAR ORBITS AND A LONG PERIOD ECLIPSING BINARY

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—α 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.

PHASES Differential Astrometry and Iodine Cell Radial Velocities of the κ Pegasi Triple Star System

κ Pegasi is a well-known, nearby triple star system. It consists of a wide pair with semimajor axis = 235 mas, one component of which is a single-line spectroscopic binary (semimajor axis = 2.5 mas). Using high-precision differential astrometry and radial velocity observations, the masses for all three components are determined and the relative inclination between the wide and narrow pairs orbits is found to be 438 ± 30, just over the threshold for the three-body Kozai resonance. The system distance is determined to be 34.60 ± 0.21 pc and is consistent with trigonometric parallax measurements.

PHASES High-Precision Differential Astrometry of δ Equulei

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 δ Equ 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 μas. 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 one 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 to improve this deficiency and discuss the outlook for further study of this binary.