We present a flexible code created for imaging from the bispectrum and visibility-squared. By using a simulated annealing method, we limit the probability of converging to local chi-squared minima as can occur when traditional imaging methods are used on data sets with limited phase information. We present the results of our code used on a simulated data set utilizing a number of regularization schemes including maximum entropy. Using the statistical properties from Monte-Carlo Markov chains of images, we show how this code can place statistical limits on image features such as unseen binary companions.

A crucial diagnostic that can tell us about processes involved in the formation and dynamical evolutionof planetary systems is the angle between the rotation axis of a star and a planet's orbital angular momentum vector ("spin-orbit" alignment or "obliquity"). Here we present the first spin-orbit alignment measurement for a wide-separation exoplanetary system, namely on the directly-imaged planet $\beta$ Pictoris b. We use VLTI/GRAVITY spectro-interferometry with an astrometric accuracy of 1 $\mu$as (microarcsecond) in the Br$\gamma$ photospheric absorption line to measure the photocenter displacement associated with the stellar rotation. Taking inclination constraints from astroseismology into account, we constrain the 3-dimensional orientation of the stellar spin axis and find that $\beta$ Pic b orbits its host star on a prograde orbit. The angular momentum vectors of the stellar photosphere, the planet, and the outer debris disk are well-aligned with mutual inclinations $<3\pm5^{\circ}$, which indicates that $\beta$ Pic b formed in a system without significant primordial misalignments. Our results demonstrate the potential of infrared interferometry to measure the spin-orbit alignment for wide-separation planetary systems, probing a highly complementary regime to the parameter space accessible with the Rossiter-McLaughlin effect. If the low obliquity is confirmed by measurements on a larger sample of wide-separation planets, it would lend support to theories that explain the obliquity in Hot Jupiter systems with dynamical scattering and the Kozai-Lidov mechanism.

Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

We present models of the inner region of the circumstellar disk of RY Tau which aim to explain our near-infrared ($K$-band: $2.1\,\mu$m) interferometric observations while remaining consistent with the optical to near-infrared portions of the spectral energy distribution. Our sub-milliarcsecond resolution CHARA Array observations are supplemented with shorter baseline, archival data from PTI, KI and VLTI/GRAVITY and modeled using an axisymmetric Monte Carlo radiative transfer code. The $K$-band visibilities are well-fit by models incorporating a central star illuminating a disk with an inner edge shaped by dust sublimation at $0.210\pm0.005\,$au, assuming a viewing geometry adopted from millimeter interferometry ($65^{\circ}$ inclined with a disk major axis position angle of $23^{\circ}$). This sublimation radius is consistent with that expected of Silicate grains with a maximum size of $0.36-0.40\,\mu$m contributing to the opacity and is an order of magnitude further from the star than the theoretical magnetospheric truncation radius. The visibilities on the longest baselines probed by CHARA indicate that we lack a clear line-of-sight to the stellar photosphere. Instead, our analysis shows that the central star is occulted by the disk surface layers close to the sublimation rim. While we do not see direct evidence of temporal variability in our multi-epoch CHARA observations, we suggest the aperiodic photometric variability of RY~Tau is likely related temporal and/or azimuthal variations in the structure of the disk surface layers.

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.

The innermost astronomical unit in protoplanetary disks is a key region for stellar and planet formation, as exoplanet searches have shown a large occurrence of close-in planets that are located within the first au around their host star. We aim to reveal the morphology of the disk inner rim using near-infrared interferometric observations with milli-arcsecond resolution provided by infrared interferometry. We provide reconstructed images of 15 objects selected from the Herbig AeBe survey carried out with PIONIER at the VLTI, using SPARCO. We find that 40% of the systems are centrosymmetric at the angular resolution of the observations. For the rest of the objects, we find evidence for asymmetric emission due to moderate-to-strong inclination of a disk-like structure for 30% of the objects and noncentrosymmetric morphology due to a nonaxisymmetric and possibly variable environment (30%). Among the systems with a disk-like structure, 20% show a resolved dust-free cavity. The image reconstruction process is a powerful tool to reveal complex disk inner rim morphologies. At the angular resolution reached by near-infrared interferometric observations, most of the images are compatible with a centrally peaked emission (no cavity). For the most resolved targets, image reconstruction reveals morphologies that cannot be reproduced by generic parametric models. Moreover, the nonaxisymmetric disks show that the spatial resolution probed by optical interferometers makes the observations of the near-infrared emission sensitive to temporal evolution with a time-scale down to a few weeks. The evidence of nonaxisymmetric emission that cannot be explained by simple inclination and radiative transfer effects requires alternative explanations, such as a warping of the inner disks. Interferometric observations can, therefore, be used to follow the evolution of the asymmetry of those disks at a sub-au scale.

In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.

We present the visual orbit of the double-lined eclipsing binary, HD 185912, from long baseline interferometry with the CHARA Array. We also obtain echelle spectra from the Apache Point observatory to update the spectroscopic orbital solution and analyze new photometry from Burggraaff et al. to model the eclipses. By combining the spectroscopic and visual orbital solutions, we find component masses of M1 = 1.361 +/- 0.004 Msun and M2 = 1.331 +/- 0.004 Msun, and a distance of d = 40.75 +/- 0.30 pc from orbital parallax. From the light curve solution, we find component radii of R1 = 1.348 +/- 0.016 Rsun and R2 = 1.322 +/- 0.016 Rsun. By comparing these observed parameters to stellar evolution models, we find that HD 185912 is a young system near the zero age main sequence with an estimated age of 500 Myr.

We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.

Exoplanet science is one of the most thriving fields of modern astrophysics. A major goal is the atmospheric characterization of dozens of small, terrestrial exoplanets in order to search for signatures in their atmospheres that indicate biological activity, assess their ability to provide conditions for life as we know it, and investigate their expected atmospheric diversity. None of the currently adopted projects or missions, from ground or in space, can address these goals. In this White Paper we argue that a large space-based mission designed to detect and investigate thermal emission spectra of terrestrial exoplanets in the MIR wavelength range provides unique scientific potential to address these goals and surpasses the capabilities of other approaches. While NASA might be focusing on large missions that aim to detect terrestrial planets in reflected light, ESA has the opportunity to take leadership and spearhead the development of a large MIR exoplanet mission within the scope of the "Voyage 2050" long-term plan establishing Europe at the forefront of exoplanet science for decades to come. Given the ambitious science goals of such a mission, additional international partners might be interested in participating and contributing to a roadmap that, in the long run, leads to a successful implementation. A new, dedicated development program funded by ESA to help reduce development and implementation cost and further push some of the required key technologies would be a first important step in this direction. Ultimately, a large MIR exoplanet imaging mission will be needed to help answer one of mankind's most fundamental questions: "How unique is our Earth?"

Diffraction fundamentally limits our ability to image and characterize exoplanets. Current and planned coronagraphic searches for exoplanets are making incredible strides but are fundamentally limited by the inner working angle of a few lambda/D. Some crucial topics, such as demographics of exoplanets within the first 50 Myr and the infrared characterization of terrestrial planets, are beyond the reach of the single aperture angular resolution for the foreseeable future. Interferometry offers some advantages in exoplanet detection and characterization and we explore in this white paper some of the potential scientific breakthroughs possible. We demonstrate here that investments in 'exoplanet interferometry' could open up new possibilities for speckle suppression through spatial coherence, a giant boost in astrometric precision for determining exoplanet orbits, ability to take a census of young giant exoplanets (clusters <50 Myr age), and an unrivaled potential for infrared nulling from space to detect terrestrial planets and search for atmospheric biomarkers. All signs point to an exciting future for exoplanets and interferometers, albeit a promise that will take decades to fulfill.

New images of young stars are revolutionizing our understanding of planet formation. ALMA detects large grains in planet-forming disks with few AU scale resolution and scattered light imaging with extreme adaptive optics systems reveal small grains suspended on the disks' flared surfaces. Tantalizing evidence for young exoplanets is emerging from line observations of CO and H-alpha. In this white paper, we explore how even higher angular resolution can extend our understanding of the key stages of planet formation, to resolve accreting circumplanetary disks themselves, and to watch planets forming in situ for the nearest star-forming regions. We focus on infrared observations which are sensitive to thermal emission in the terrestrial planet formation zone and allow access to molecular tracers in warm ro-vibrational states. Successful planet formation theories will not only be able to explain the diverse features seen in disks, but will also be consistent with the rich exoplanet demographics from RV and transit surveys. While we are far from exhausting ground-based techniques, the ultimate combination of high angular resolution and high infrared sensitivity can only be achieved through mid-infrared space interferometry.

The ultimate astronomical observatory would be a formation flying interferometer in space, immune to atmospheric turbulence and absorption, free from atmospheric and telescope thermal emission, and reconfigurable to adjust baselines according to the required angular resolution. Imagine the near/mid-infrared sensitivity of the JWST and the far-IR sensitivity of Herschel but with ALMA-level angular resolution, or imagine having the precision control to null host star light across 250m baselines and to detect molecules from the atmospheres of nearby exo-Earths. With no practical engineering limit to the formation's size or number of telescopes in the array, formation flying interferometry will revolutionize astronomy and this White Paper makes the case that it is now time to accelerate investments in this technological area. Here we provide a brief overview of the required technologies needed to allow light to be collected and interfered using separate spacecrafts. We emphasize the emerging role of inexpensive smallSat projects and the excitement for the LISA Gravitational Wave Interferometer to push development of the required engineering building-blocks. We urge the Astro2020 Decadal Survey Committee to highlight the need for a small-scale formation flying space interferometer project to demonstrate end-to-end competency with a timeline for first stellar fringes by the end of the decade.

An international group of scientists has begun planning for the Planet Formation Imager (PFI, www.planetformationimager.org), a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions and taking planetary system 'snapshots' of young systems to understand exoplanet architectures. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band beam combiner will be especially sensitive to thermal emission from young exoplanets (and their circumplanetary disks) with a high spectral resolution mode to probe the kinematics of CO and H2O gas. In this brief White Paper, we summarize the main science goals of PFI, define a baseline PFI architecture that can achieve those goals, and identify key technical challenges that must be overcome before the dreams of PFI can be realized within the typical cost envelope of a major observatory. We also suggest activities over the next decade at the flagship US facilities (CHARA, NPOI, MROI) that will help make the Planet Formation Imager facility a reality. The key takeaway is that infrared interferometry will require new experimental telescope designs that can scale to 8 m-class with the potential to reduce per area costs by x10, a breakthrough that would also drive major advances across astronomy.

T Tauri stars are low-mass young stars whose disks provide the setting for planet formation. Despite this, their structure is poorly understood. We present new infrared interferometric observations of the SU Aurigae circumstellar environment that offer 3 x higher resolution and better baseline position angle coverage over previous observations. We investigate the characteristics of circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. The CHARA array offers opportunities for long baseline observations, with baselines up to 331 m. Using the CLIMB 3-telescope combiner in the K-band allows us to measure visibilities as well as closure phase. We undertook image reconstruction for model-independent analysis, and geometric modeling. Additionally, the fitting of radiative transfer models constrains the physical parameters of the disk. For the first time, a dusty disk wind is introduced to the radiative transfer code TORUS to model protoplanetary disks. Our implementation is motivated by theoretical dusty disk winds, where magnetic field lines drive dust above the disk plane close to the sublimation zone. Image reconstruction reveals an inclined disk with slight asymmetry along its minor-axis, likely due to inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission/scattering of the far-side. Geometric modelling of a skewed ring finds the inner rim at 0.17+/-0.02 au with an inclination of 50.9+/-1.0 degrees and minor axis position angle 60.8+/-1.2 degrees. Radiative transfer modelling shows a flared disk with an inner radius at 0.18 au which implies a grain size of 0.4 um and a scale height of 15.0 au at 100 au. Among the tested radiative transfer models, only the dusty disk wind successfully accounts for the K-band excess by introducing dust above the mid-plane.

This work presents a study of two Herbig Ae transitional discs, Oph IRS 48 and HD 169142; which both have reported rings in their dust density distributions. We use Keck-II/NIRC2 adaptive optics imaging observations in the L' filter (3.8 micron) to probe the regions of these discs inwards of ~20AU from the star. We introduce our method for investigating these transitional discs, which takes a forward modelling approach: making a model of the disc (using the Monte Carlo radiative transfer code RADMC), convolving it with point-spread functions of calibrator stars, and comparing the convolved models with the observational data. The disc surface density parameters are explored with a Monte Carlo Markov Chain technique. Our analysis recovers emission from both of the discs interior to the well known optically thick walls, modelled as a ring of emission at ~15AU in Oph IRS 48, and ~7AU for HD 169142, and identifies asymmetries in both discs. Given the brightness of the near-symmetric rings compared to the reported companion candidates, we suggest that the reported companion candidates can be interpreted as slightly asymmetric disc emission or illumination.

Super-high spatial resolution observations in the infrared are now enabling major advances in our understanding of supermassive black hole systems at the centers of galaxies. Infrared interferometry, reaching resolutions of milliarcseconds to sub-milliarcseconds, is drastically changing our view of the central structure from a static to a very dynamic one by spatially resolving to the pc-scale. We are also starting to measure the dynamical structure of fast moving gas clouds around active supermassive black holes at a scale of less than a light year. With further improvements on resolution and sensitivity, we will be able to directly image the exact site of the black hole's feedback to its host galaxy, and quantify the effect of such interaction processes. Near-future high angular resolution studies will definitely advance our mass determinations for these black holes, and we might even witness the existence of binary black hole systems at the center of galaxies.

We implement an electron avalanche photodiode (e-APD) in the MIRC-X instrument, upgrade of the 6-telescope near-infrared imager MIRC, at the CHARA array. This technology should improve the sensitivity of near-infrared interferometry. We first used the classical Mean-Variance analysis to measure the system gain and the amplification gain. We then developed a physical model of the statistical distribution of the camera output signal. This model is based on multiple convolutions of the Poisson statistic, the intrinsic avalanche gain distribution, and the observed distribution of the background signal. At low flux level, this model constraints independently the incident illumination level, the total gain, and the excess noise factor of the amplification. We measure a total transmission of $48\pm3\%$ including the cold filter and the Quantum Efficiency. We measure a system gain of 0.49 ADU/e, a readout noise of $10$ ADU, and amplification gains as high as 200. These results are consistent between the two methods and therefore validate our modeling approach. The measured excess noise factor based on the modeling is $1.47\pm0.03$, with no obvious dependency with flux level or amplification gain. The presented model allows measuring the characteristics of the e-APD array at low flux level independently of preexisting calibration. With $<0.3$\,electron equivalent readout noise at kilohertz frame rates, we confirm the revolutionary performances of the camera with respect to the PICNIC or HAWAII technologies. However, the measured excess noise factor is significantly higher than the one claimed in the literature ($<$1.25), and explains why counting multiple photons remains challenging with this camera.

Direct determination of fundamental stellar parameters has many profound and wide-ranging impacts throughout astrophysics. These determinations are rooted in high angular resolution observations. In particular, as long-baseline optical interferometry has matured over the past decade, increasingly large survey samples are serving to empirically ground the basic parameters of these building blocks of the universe. True imaging and improved parametric fitting are becoming routinely available, an essential component of fully characterizing stars, stellar environments, and planets these stars may host.

Evolved stars dominate galactic spectra, enrich the galactic medium, expand to change their planetary systems, eject winds of a complex nature, produce spectacular nebulae and illuminate them, and transfer material between binary companions. While doing this, they fill the HR diagram with diagnostic loops that write the story of late stellar evolution. Evolved stars sometimes release unfathomable amounts of energy in neutrinos, light, kinetic flow, and gravitational waves. During these late-life times, stars evolve complexly, with expansion, convection, mixing, pulsation, mass loss. Some processes have virtually no spatial symmetries, and are poorly addressed with low-resolution measurements and analysis. Even a "simple" question as how to model mass loss resists solution. However, new methods offer increasingly diagnostic tools. Astrometry reveals populations and groupings. Pulsations/oscillations support study of stellar interiors. Optical/radio interferometry enable 2-3d imagery of atmospheres and shells. Bright stars with rich molecular spectra and velocity fields are a ripe opportunity for imaging with high spatial and spectral resolution, giving insight into the physics and modeling of later stellar evolution.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We are now in an era where we can image details on the surfaces of stars. When resolving stellar surfaces, we see that every surface is uniquely complicated. Each imaged star provides insight into not only the stellar surface structures, but also the stellar interiors suggesting constraints on evolution and dynamo models. As more resources become operational in the coming years, imaging stellar surfaces should become commonplace for revealing the true nature of stars. Here, we discuss the main types of stars for which imaging surface features is currently useful and what improved observing techniques would provide for imaging stellar surface features.

Cool, evolved stars are the main source of chemical enrichment of the interstellar medium (ISM), and understanding their mass loss and structure offers a unique opportunity to study the cycle of matter in the Universe. Pulsation, convection, and other dynamic processes in cool evolved stars create an atmosphere where molecules and dust can form, including those necessary to the formation of life (e.g.~Carbon-bearing molecules). Understanding the structure and composition of these stars is thus vital to several aspects of stellar astrophysics, ranging from ISM studies to modeling young galaxies and exoplanet research. Recent modeling efforts and increasingly precise observations now reveal that our understanding of cool stars photospheric, chromospheric, and atmospheric structures is limited by inadequate knowledge of the dynamic and chemical processes at work. Here we outline promising scientific opportunities for the next decade. We identify and discuss the following main opportunities: (1) identify and model the physical processes that must be included in current 1D and 3D atmosphere models of cool, evolved stars; (2) refine our understanding of photospheric, chromospheric, and outer atmospheric regions of cool evolved stars, their properties and parameters, through high-resolution spectroscopic observations, and interferometric observations at high angular resolution; (3) include the neglected role of chromospheric activity in the mass loss process of red giant branch and red super giant stars, and understand the role played by their magnetic fields; (4) identify the important shaping mechanisms for planetary nebulae and their relation with the parent asymptotic giant branch stars.

We present the visual orbit of the double-lined spectroscopic binary HD 224355 from interferometric observations with the CHARA Array, as well as an updated spectroscopic analysis using echelle spectra from the Apache Point Observatory 3.5m telescope. By combining the visual and spectroscopic orbital solutions, we find the binary components to have masses of M1 = 1.626 +/- 0.005 Msun and M2 = 1.608 +/- 0.005 Msun, and a distance of d = 63.98 +/- 0.26 pc. Using the distance and the component angular diameters found by fitting spectrophotometry from the literature to spectral energy distribution models, we estimate the stellar radii to be R1 = 2.65 +/- 0.21 Rsun and R2 = 2.47 +/- 0.23 Rsun. We then compare these observed fundamental parameters to the predictions of stellar evolution models, finding that both components are evolved towards the end of the main sequence with an estimated age of 1.9 Gyr.

The disks around some Herbig Be stars have been observed to be more compact than the expected dust sublimation radius for such objects, with highly refractory dust grains and optically thick gas emission having been proposed as possible explanations for this phenomenon. Our aim is to search for direct evidence for the presence of hot gas inside of the expected dust sublimation radius of MWC147. By combining VLTI/AMBER spectro-interferometry (R=12,000) with CRIRES spectroscopy (R=100,000) we can both spectrally and spatially resolve the Br-gamma line-emitting gas around MWC147. Our interferometric visibility modelling of MWC147 indicates the presence of a compact continuum disk with a close to face-on orientation. We model the continuum with an inclined Gaussian, as well as a ring with a radius of 0.60mas (0.39au) which is well within the expected dust sublimation radius of 1.52au. We detect no significant change in the measured visibilities across the Br-gamma line, indicating that the line-emitting gas is located in the same region as the continuum-emitting disk. We fit our AMBER spectro-interferometry data with a kinematic model of a disk in Keplerian rotation, with both the line-emitting and continuum-emitting components of the disk originating from the same compact region close to the central star. The presence of line-emitting gas in the same region as the K-band continuum supports the interpretation that the K-band continuum traces an optically-thick gas disk. Our spatially and spectrally resolved observations of MWC147 reveal that the K-band continuum and Br-gamma emission both originate from a similar compact region, with Br-gamma emitted from the accretion disk or disk wind region and exhibiting a rotational velocity profile. We conclude that we detect the presence of a compact, gaseous accretion disk in Keplerian rotation around MWC147.

We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature.

We aim at detecting and characterizing the main-sequence companions of a sample of known and suspected Galactic binary Cepheids. We used the multi-telescope interferometric combiners MIRC and PIONIER to detect and measure the astrometric positions of the high-contrast companions orbiting 16 bright Galactic Cepheids. We made use of the CANDID algorithm to search for the companions and set detection limits. We also present new high-precision RVs which were used to fit radial pulsation and orbital velocities. We report the detection of the companions orbiting the Cepheids U Aql, BP Cir, and S Mus, possible detections for FF Aql, Y Car, BG Cru, X Sgr, V350 Sgr, and V636 Sco, while no component is detected around U Car, YZ Car, T Mon, R Mus, S Nor, W Sgr and AH Vel. For U Aql and S Mus, we performed a preliminary orbital fit combining astrometry with new high-precision RVs, providing the full set of orbital and pulsation parameters. Assuming the distance, we estimated preliminary masses of M(U Aql) = 4.97+/-0.62Msun and M(S Mus) = 4.63+/-0.99Msun. For YZ Car, W Sgr, V350 Sgr, and V636 Sco, we revised the spectroscopic orbits using new high-precision RVs, while we updated the pulsation parameters for BP Cir, BG Cru, S Nor and AH Vel. Our interferometric observations also provide measurements of the angular diameters, that can be used in a Baade-Wesselink type analysis. We have now several astrometric detections of Cepheid companions. When RVs of the companions will be available, such systems will provide accurate and independent masses and distances. Orbital parallaxes with an accuracy <5% will be particularly useful for a better calibration of the P-L relation. The final Gaia parallaxes will be also particularly helpful for single-line spectroscopic systems, where mass and distance are degenerate. Mass measurements are necessary for a better understanding of the age and evolution of Cepheids

V1247 Orionis harbours a pre-transitional disc with a partially cleared gap. Earlier interferometric and polarimetric observations revealed strong asymmetries both in the gap region and in the outer disc. The presence of a companion was inferred to explain these asymmetric structures and the ongoing disc clearing. Using an extensive set of multi-wavelength and multi-epoch observations we aimed to identify the origin of the previously detected asymmetries. We have observed V1247 Ori at three epochs spanning $\sim678$ days using sparse aperture masking interferometry with Keck/NIRC2 and VLT/NACO. In addition, we search for signs of accretion through VLT/SPHERE-ZIMPOL spectral differential imaging in H$\alpha$ and R-band continuum. Our SMA sub-millimetre interferometry in 880 $\mu$m continuum and in the CO(3-2) line allows us to constrain the orientation and direction of rotation of the outer disc. We find the L'-band emission to be dominated by static features which trace forward-scattered dust emission from the inner edge of the outer disc located to the north-east. In H- and K-band, we see evidence for a companion candidate that moved systematically by 45$^{\circ}$ within the first $\sim$345 days. The separation of the companion candidate is not well constrained, but the observed position angle change is consistent with Keplerian motion of a body located on a 6 au orbit. From the SMA CO moment map, the location of the disc rim, and the detected orbital motion, we deduced the three-dimensional orientation of the disc. We see no indication of accretion in H$\alpha$ and set upper limits for an accreting companion. The measured contrast of the companion candidate in H and K is consistent with an actively accreting protoplanet. Hence, we identify V1247 Ori as a unique laboratory for studying companion-disc interactions and disc clearing.

The physical processes occurring within the inner few astronomical units of proto-planetary disks surrounding Herbig Ae stars are crucial to setting the environment in which the outer planet-forming disk evolves and put critical constraints on the processes of accretion and planet migration. We present the most complete published sample of high angular resolution H- and K-band observations of the stars HD 163296 and HD 190073, including 30 previously unpublished nights of observations of the former and 45 nights of the latter with the CHARA long-baseline interferometer, in addition to archival VLTI data. We confirm previous observations suggesting significant near-infrared emission originates within the putative dust evaporation front of HD 163296 and show this is the case for HD 190073 as well. The H- and K-band sizes are the same within $(3 \pm 3)\%$ for HD 163296 and within $(6 \pm 10)\%$ for HD 190073. The radial surface brightness profiles for both disks are remarkably Gaussian-like with little or no sign of the sharp edge expected for a dust evaporation front. Coupled with spectral energy distribution analysis, our direct measurements of the stellar flux component at H and K bands suggest that HD 190073 is much younger (<400 kyr) and more massive (~5.6 M$_\odot$) than previously thought, mainly as a consequence of the new Gaia distance (891 pc).

We present results for the inflated super-Neptune HATS-8b from MOPSS, The Michigan Optical Planetary Spectra Survey. This program is aimed at creating a database of optical planetary transmission spectra all observed, reduced, and analyzed with a uniform method for the benefit of enabling comparative exoplanet studies. HATS-8b orbits a G dwarf and is a low density super-Neptune, with a radius of 0.873 R$_{Jup}$, a mass of 0.138 M$_{Jup}$, and a density of 0.259 g/cm$^3$. Two transits of HATS-8b were observed in July and August of 2017 with the IMACS instrument on the Magellan Baade 6.5m telescope. We find an enhanced scattering slope that differs between our two nights. These slopes are stronger than one due only to Rayleigh scattering and cannot be fully explained by unocculted star spots. We explore the impact of condensates on the scattering slope and determine that MnS particulates smaller than 10$^{-2}\mu$m can explain up to 80\% of our measured slope if the planet is warmer than equilibrium, or 50\% of the slope at the equilibrium temperature of the planet for a low mean molecular weight atmosphere. The scattering slope that we observe is thus beyond even the most extreme case predicted by theory. We suggest further follow up on this target and host star to determine if the temporal variation of the slope is primarily due to stellar or planetary effects, and to better understand what these effects may be.

Cepheid stars play a considerable role as extragalactic distances indicators, thanks to the simple empirical relation between their pulsation period and their luminosity. They overlap with that of secondary distance indicators, such as Type Ia supernovae, whose distance scale is tied to Cepheid luminosities. However, the Period-Luminosity (P-L) relation still lacks a calibration to better than 5%. Using an original combination of interferometric astrometry with optical and ultraviolet spectroscopy, we measured the geometrical distance d = 720.35+/-7.84 pc of the 3.33 d period Cepheid V1334 Cyg with an unprecedented accuracy of +/-1 %, providing the most accurate distance for a Cepheid. Placing this star in the P-L diagram provides an independent test of existing period-luminosity relations. We show that the secondary star has a significant impact on the integrated magnitude, particularly at visible wavelengths. Binarity in future high precision calibrations of the P-L relations is not negligible, at least in the short-period regime. Subtracting the companion flux leaves V1334 Cyg in marginal agreement with existing photometric-based P-L relations, indicating either an overall calibration bias or a significant intrinsic dispersion at a few percent level. Our work also enabled us to determine the dynamical masses of both components, M1 = 4.288 +/- 0.133 Msun (Cepheid) and M2 = 4.040 +/- 0.048 Msun (companion), providing the most accurate masses for a Galactic binary Cepheid system.

We present comprehensive models of Herbig Ae star, HD 142666, which aim to simultaneously explain its spectral energy distribution (SED) and near-infrared (NIR) interferometry. Our new sub-milliarcsecond resolution CHARA (CLASSIC and CLIMB) interferometric observations, supplemented with archival shorter baseline data from VLTI/PIONIER and the Keck Interferometer, are modeled using centro-symmetric geometric models and an axisymmetric radiative transfer code. CHARA's 330 m baselines enable us to place strong constraints on the viewing geometry, revealing a disk inclined at 58 degrees from face-on with a 160 degree major axis position angle. Disk models imposing vertical hydrostatic equilibrium provide poor fits to the SED. Models accounting for disk scale height inflation, possibly induced by turbulence associated with magneto-rotational instabilities, and invoking grain growth to >1 micron size in the disk rim are required to simultaneously reproduce the SED and measured visibility profile. However, visibility residuals for our best model fits to the SED indicate the presence of unexplained NIR emission, particularly along the apparent disk minor axis, while closure phase residuals indicate a more centro-symmetric emitting region. In addition, our inferred 58 degree disk inclination is inconsistent with a disk-based origin for the UX Ori-type variability exhibited by HD 142666. Additional complexity, unaccounted for in our models, is clearly present in the NIR-emitting region. We propose the disk is likely inclined toward a more edge-on orientation and/or an optically thick outflow component also contributes to the NIR circumstellar flux.

The Planet Formation Imager (PFI) is a near- and mid-infrared interferometer project with the driving science goal of imaging directly the key stages of planet formation, including the young proto-planets themselves. Here, we will present an update on the work of the Science Working Group (SWG), including new simulations of dust structures during the assembly phase of planet formation and quantitative detection efficiencies for accreting and non-accreting young exoplanets as a function of mass and age. We use these results to motivate two reference PFI designs consisting of a) twelve 3m telescopes with a maximum baseline of 1.2km focused on young exoplanet imaging and b) twelve 8m telescopes optimized for a wider range of young exoplanets and protoplanetary disk imaging out to the 150K H2O ice line. Armed with 4x8m telescopes, the ESO/VLTI can already detect young exoplanets in principle and projects such as MATISSE, Hi-5 and Heimdallr are important PFI pathfinders to make this possible. We also discuss the state of technology development needed to make PFI more affordable, including progress towards new designs for inexpensive, small field-of-view, large aperture telescopes and prospects for Cubesat-based space interferometry.

We present the design for MYSTIC, the Michigan Young STar Imager at CHARA. MYSTIC will be a K-band, cryogenic, 6-beam combiner for the Georgia State University CHARA telescope array. The design follows the image-plane combination scheme of the MIRC instrument where single-mode fibers bring starlight into a non-redundant fringe pattern to feed a spectrograph. Beams will be injected in polarization-maintaining fibers outside the cryogenic dewar and then be transported through a vacuum feedthrough into the ~220K cold volume where combination is achieved and the light is dispersed. We will use a C-RED One camera (First Light Imaging) based on the eAPD SAPHIRA detector to allow for near-photon-counting performance. We also intend to support a 4-telescope mode using a leftover integrated optics component designed for the VLTI-GRAVITY experiment, allowing better sensitivity for the faintest targets. Our primary science driver motivation is to image disks around young stars in order to better understand planet formation and how forming planets might influence disk structures.

The Planet Formation Imager (PFI, www.planetformationimager.org) is a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band combiner will be especially sensitive to thermal emission from young exoplanets (and their disks) with a high spectral resolution mode to probe the kinematics of CO and H2O gas. In this paper, we give an overview of the main science goals of PFI, define a baseline PFI architecture that can achieve those goals, point at remaining technical challenges, and suggest activities today that will help make the Planet Formation Imager facility a reality.

We characterize a near-infrared C-RED ONE camera from First Light Imaging (FLI). This camera uses a SAPHIRA electron avalanche photo-diode array (e-APD) from Leonardo (previously Selex). To do so, we developed a model of the signal distribution. This model allows a measurement of the gain and the Excess Noise Factor (ENF) independently of preexisting calibration such as the system gain. The results of this study show a gain which is 0.53 +/- 0.04 times the gain reported by the manufacturer. The measured ENF is 1.47 +/- 0.03 when we expected 1.25. For an avalanche gain of 60 and a frame rate larger than 100 Hz, the total noise can be lower than 1 e-/frame/pixel. The lowest dark current level is 90e-/s/pixel, in agreement with the expected H-band background passing through the camera window. These performance values provide a significant improvement compared to earlier-generation PICNIC camera and allowed us to improve the performance of the Michigan infrared combiner (MIRC) instrument at the Center for High Angular Resolution Astronomy (CHARA), as part of our MIRC-X instrumentation project.

Hi-5 is a high-contrast (or high dynamic range) infrared imager project for the VLTI. Its main goal is to characterize young extra-solar planetary systems and exozodiacal dust around southern main-sequence stars. In this paper, we present an update of the project and key technology pathways to improve the contrast achieved by the VLTI. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5~$\mu$m) and advanced fringe tracking strategies. We also address the strong exoplanet science case (young exoplanets, planet formation, and exozodiacal disks) offered by this wavelength regime as well as other possible science cases such as stellar physics (fundamental parameters and multiplicity) and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed.

In the era of high-angular resolution astronomical instrumentation, where long and very long baseline interferometers (constituted by many, $\sim$ 20 or more, telescopes) are expected to work not only in the millimeter and submillimeter domain, but also at near and mid infrared wavelengths (experiments such as the Planet Formation Imager, PFI, see Monnier et al. 2018 for an update on its design); any promising strategy to alleviate the costs of the individual telescopes involved needs to be explored. In a recent collaboration between engineers, experimental physicists and astronomers in Valparaiso, Chile, we are gaining expertise in the production of light carbon fiber polymer reinforced mirrors. The working principle consists in replicating a glass, or other substrate, mandrel surface with the mirrored adequate curvature, surface characteristics and general shape. Once the carbon fiber base has hardened, previous studies have shown that it can be coated (aluminum) using standard coating processes/techniques designed for glass-based mirrors. The resulting surface quality is highly dependent on the temperature and humidity control among other variables. Current efforts are focused on improving the smoothness of the resulting surfaces to meet near/mid infrared specifications, overcoming, among others, possible deteriorations derived from the replication process. In a second step, at the validation and quality control stage, the mirrors are characterized using simple/traditional tools like spherometers (down to micron precision), but also an optical bench with a Shack-Hartman wavefront sensor. This research line is developed in parallel with a more classical glass-based approach, and in both cases we are prototyping at the small scale of few tens of cms. We here present our progress on these two approaches.

We present the first results from MOPSS, The Michigan Optical Planetary Spectra Survey, aimed at creating a database of optical planetary transmission spectra all observed, reduced, and analyzed with a uniform method for the benefit of enabling comparative exoplanet studies. We discuss our methods, and present results for our first two targets observed with the Magellan Baade 6.5m telescope, one transit of the Hot Jupiter WASP-4b and two transits of the Hot Saturn WASP-52b. Both targets present flat, featureless spectra, corresponding to the presence of aerosols. We find that the cloud decks must begin no lower than 10$^{-4}$ bar for both planets. For WASP-52b, we also consider the effects of star spots on the transmission spectrum, including unocculted spots and spots on the stellar limb influencing the light curve limb darkening parameters. We discuss the usefulness of this program in the coming James Webb Space Telescope era.

Mass exchange and mass loss in close binaries can significantly affect their evolution, but a complete self-consistent theory of these processes is still to be developed. Processes such as radiative shielding due to a~hot-spot region, or a~hydrodynamical interaction of different parts of the gas stream have been studied previously. In order to test the respective predictions, it is necessary to carry out detailed observations of binaries undergoing the largescale mass exchange, especially for those that are in the rapid transfer phase. \bla is an archetype of such a system, having a long and rich observational history. Our goal for this first study is to quantitatively estimate the geometry and physical properties of the optically thick components, namely the Roche-lobe filling mass-losing star, and the accretion disk surrounding the mass-gaining star of \blae. A series of continuum visible and NIR spectro-interferometric observations by the NPOI, CHARA/MIRC and VEGA instruments covering the whole orbit of \bla acquired during a~two-week campaign in 2013 were complemented with \ubvr\ photometric observations acquired during a three-year monitoring of the system. We included NUV and FUV observations from OAO~A-2, IUE, and Voyager satellites.

Despite promising astrometric signals, to date there has been no success in direct imaging of a hypothesized third member of the Sirius system. Using the Clio instrument and MagAO adaptive optics system on the Magellan Clay 6.5 m telescope, we have obtained extensive imagery of Sirius through a vector apodizing phase plate (vAPP) coronagraph in a narrowband filter at 3.9 microns. The vAPP coronagraph and MagAO allow us to be sensitive to planets much less massive than the limits set by previous non-detections. However, analysis of these data presents challenges due to the target's brightness and unique characteristics of the instrument. We present a comparison of dimensionality reduction techniques to construct background illumination maps for the whole detector using the areas of the detector that are not dominated by starlight. Additionally, we describe a procedure for sub-pixel alignment of vAPP data using a physical-optics-based model of the coronagraphic PSF.

MIRC-X is an upgrade of the six-telescope infrared beam combiner at the CHARA telescope array, the world's largest baseline interferometer in the optical/infrared, located at the Mount Wilson Observatory in Los Angeles. The upgraded instrument features an ultra-low noise and fast frame rate infrared camera (SAPHIRA detector) based on e-APD technology. We report the MIRC-X sensitivity upgrade work and first light results in detail focusing on the detector characteristics and software architecture.

MIRC-X is a new beam combination instrument at the CHARA array that enables 6-telescope interferometric imaging on object classes that until now have been out of reach for milliarcsecond-resolution imaging. As part of an instrumentation effort lead by the University of Exeter and University of Michigan, we equipped the MIRC instrument with an ultra-low read-noise detector system and extended the wavelength range to the J and H-band. The first phase of the MIRC-X commissioning was successfully completed in June 2017. In 2018 we will commission polarisation control to improve the visibility calibration and implement a 'cross-talk resiliant' mode that will minimise visibility cross-talk and enable exoplanet searches using precision closure phases. Here we outline our key science drivers and give an overview about our commissioning timeline. We comment on operational aspects, such as remote observing, and the prospects of co-phased parallel operations with the upcoming MYSTIC combiner.

The development of high-contrast capabilities has long been recognized as one of the top priorities for the VLTI. As of today, the VLTI routinely achieves contrasts of a few 10$^{-3}$ in the near-infrared with PIONIER (H band) and GRAVITY (K band). Nulling interferometers in the northern hemisphere and non-redundant aperture masking experiments have, however, demonstrated that contrasts of at least a few 10$^{-4}$ are within reach using specific beam combination and data acquisition techniques. In this paper, we explore the possibility to reach similar or higher contrasts on the VLTI. After reviewing the state-of-the-art in high-contrast infrared interferometry, we discuss key features that made the success of other high-contrast interferometric instruments (e.g., integrated optics, nulling, closure phase, and statistical data reduction) and address possible avenues to improve the contrast of the VLTI by at least one order of magnitude. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5 $\mu$m), a sweet spot to image and characterize young extra-solar planetary systems. Finally, we address the science cases of a high-contrast VLTI imaging instrument and focus particularly on exoplanet science (young exoplanets, planet formation, and exozodiacal disks), stellar physics (fundamental parameters and multiplicity), and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed.

The next generation ground-based extremely large telescopes (ELTs) present incredible opportunities to discover and characterize diverse planetary systems, even potentially habitable worlds. Adaptive-optics assisted thermal-IR (3-14 micron) imaging is a powerful tool to study exoplanets with extant 6-12 meter telescopes. ELTs have the spatial resolution and sensitivity that offer an unparalleled expansion of the available discovery space. AO-assisted thermal-IR instruments on ELTs will be superior to JWST for high contrast imaging in the thermal-IR, and complementary to high contrast observations at shorter wavelengths, in space or with second-generation extreme AO instruments. With appropriate investments in instrumentation and pre-cursor observations, thermal-IR equipped ELTs could image the first terrestrial and super-earth planets around nearby stars, opening the door to characterization of potentially habitable planets from the ground and space.

High angular resolution observations of young stellar objects are required to study the inner astronomical units of protoplanetary disks in which the majority of planets form. As they evolve, gaps open up in the inner disk regions and the disks are fully dispersed within ~10 Myrs. MWC 614 is a pre-transitional object with a ~10au radius gap. We present a set of high angular resolution observations of this object including SPHERE/ZIMPOL polarimetric and coronagraphic images in the visible, KECK/NIRC2 near-infrared aperture masking observations and VLTI (AMBER, MIDI, and PIONIER) and CHARA (CLASSIC and CLIMB) long-baseline interferometry at infrared wavelengths. We find that all the observations are compatible with an inclined disk (i ~55deg at a position angle of ~20-30deg). The mid-infrared dataset confirms the disk inner rim to be at 12.3+/-0.4 au from the central star. We determined an upper mass limit of 0.34 Msun for a companion inside the cavity. Within the cavity, the near-infrared emission, usually associated with the dust sublimation region, is unusually extended (~10 au, 30 times larger than the theoretical sublimation radius) and indicates a high dust temperature (T~1800 K). As a possible result of companion-induced dust segregation, quantum heated dust grains could explain the extended near-infrared emission with this high temperature. Our observations confirm the peculiar state of this object where the inner disk has already been accreted onto the star exposing small particles inside the cavity to direct stellar radiation.

Combining visual and spectroscopic orbits of binary stars leads to a determination of the full 3D orbit, individual masses, and distance to the system. We present a full analysis of the evolved binary system $\delta$ Delphini using astrometric data from the MIRC and PAVO instruments on the CHARA long-baseline interferometer, 97 new spectra from the Fairborn Observatory, and 87 unpublished spectra from Lick Observatory. We determine the full set of orbital elements for $\delta$ Del, along with masses of $1.78 \pm 0.07$ $M_{\odot}$ and $1.62 \pm 0.07$ $M_{\odot}$ for each component, and a distance of $63.61 \pm 0.89$ pc. These results are important in two contexts: for testing stellar evolution models and defining the detection capabilities for future planet searches. We find that the evolutionary state of this system is puzzling, as our measured flux ratios, radii, and masses imply a $\sim$ 200 Myr age difference between the components using standard stellar evolution models. Possible explanations for this age discrepancy include mass transfer scenarios with a now ejected tertiary companion. For individual measurements taken over a span of 2 years we achieve $<10$ $\mu$-arcsecond precision on differential position with 10-minute observations. The high precision of our astrometric orbit suggests that exoplanet detection capabilities are within reach of MIRC at CHARA. We compute exoplanet detection limits around $\delta$ Del, and conclude that if this precision is extended to wider systems we should be able to detect most exoplanets $>2$ M$_{J}$ on orbits $>0.75$ AU around individual components of hot binary stars via differential astrometry.

Recent observational work has indicated that mechanisms for accretion and outflow in Herbig Ae/Be star-disk systems may differ from magnetospheric accretion (MA) as it is thought to occur in T Tauri star-disk systems. In this work, we assess the temporal evolution of spectral lines probing accretion and mass loss in Herbig Ae/Be systems and test for consistency with the MA paradigm. For two Herbig Ae/Be stars, HD 98922 (B9e) and V1295 Aql (A2e), we have gathered multi-epoch (~years) and high-cadence (~minutes) high-resolution optical spectra to probe a wide range of kinematic processes. Employing a line equivalent width evolution correlation metric introduced here, we identify species co-evolving (indicative of common line origin) via novel visualization. We interferometrically constrain often problematically degenerate parameters, inclination and inner disk radius, allowing us to focus on the structure of the wind, magnetosphere, and inner gaseous disk in radiative transfer models. Over all timescales sampled, the strongest variability occurs within the blueshifted absorption components of the Balmer series lines; the strength of variability increases with the cadence of the observations. Finally, high-resolution spectra allow us to probe substructure within the Balmer series' blueshifted absorption components: we observe static, low-velocity features and time-evolving features at higher velocities. Overall, we find the observed line morphologies and variability are inconsistent with a scaled-up T Tauri MA scenario. We suggest that as magnetic field structure and strength change dramatically with increasing stellar mass from T Tauri to Herbig Ae/Be stars, so too may accretion and outflow processes.

Herbig Ae/Be stars are young, pre-main-sequence stars that sample the transition in structure and evolution between low- and high-mass stars, providing a key test of accretion processes in higher-mass stars. Few Herbig Ae/Be stars have detected magnetic fields, calling into question whether the magnetospheric accretion paradigm developed for low-mass stars can be scaled to higher masses. We present He I 10830 \AA\ line profiles for 64 Herbig Ae/Be stars with a magnetic field measurement in order to test magnetospheric accretion in the physical regime where its efficacy remains uncertain. Of the 5 stars with a magnetic field detection, 1 shows redshifted absorption, indicative of infall, and 2 show blueshifted absorption, tracing mass outflow. The fraction of redshifted and blueshifted absorption profiles in the non-magnetic Herbig Ae/Be stars is remarkably similar, suggesting that the stellar magnetic field does not affect gas kinematics traced by He I 10830 \AA. Line profile morphology does not correlate with the luminosity, rotation rate, mass accretion rate, or disk inclination. Only the detection of a magnetic field and a nearly face-on disk inclination show a correlation (albeit for few sources). This provides further evidence for weaker dipoles and more complex field topologies as stars develop a radiative envelope. The small number of magnetic Herbig Ae/Be stars has already called into question whether magnetospheric accretion can be scaled to higher masses; accretion signatures are not substantially different in magnetic Herbig Ae/Be stars, casting further doubt that they accrete in the same manner as classical T Tauri stars.

The radial drift problem constitutes one of the most fundamental problems in planet formation theory, as it predicts particles to drift into the star before they are able to grow to planetesimal size. Dust-trapping vortices have been proposed as a possible solution to this problem, as they might be able to trap particles over millions of years, allowing them to grow beyond the radial drift barrier. Here, we present ALMA 0.04"-resolution imaging of the pre-transitional disk of V1247 Orionis that reveals an asymmetric ring as well as a sharply-confined crescent structure, resembling morphologies seen in theoretical models of vortex formation. The asymmetric ring (at 0.17"=54 au separation from the star) and the crescent (at 0.38"=120 au) seem smoothly connected through a one-armed spiral arm structure that has been found previously in scattered light. We propose a physical scenario with a planet orbiting at $\sim0.3$"$\approx$100 au, where the one-armed spiral arm detected in polarised light traces the accretion stream feeding the protoplanet. The dynamical influence of the planet clears the gap between the ring and the crescent and triggers two vortices that trap mm-sized particles, namely the crescent and the bright asymmetry seen in the ring. We conducted dedicated hydrodynamics simulations of a disk with an embedded planet, which results in similar spiral-arm morphologies as seen in our scattered light images. At the position of the spiral wake and the crescent we also observe $^{12}$CO (3-2) and H$^{12}$CO$^{+}$ (4-3) excess line emission, likely tracing the increased scale-height in these disk regions.

Nearby, active stars with relatively rapid rotation and large starspot structures offer the opportunity to compare interferometric, spectroscopic, and photometric imaging techniques. In this paper, we image a spotted star with three different methods for the first time. The giant primary star of the RS Canum Venaticorum binary $\sigma$ Geminorum ($\sigma$ Gem) was imaged for two epochs of interferometric, high-resolution spectroscopic, and photometric observations. The light curves from the reconstructions show good agreement with the observed light curves, supported by the longitudinally-consistent spot features on the different maps. However, there is strong disagreement in the spot latitudes across the methods.

Sunspots are cool areas caused by strong surface magnetic fields inhibiting convection. Moreover, strong magnetic fields can alter the average atmospheric structure, degrading our ability to measure stellar masses and ages. Stars more active than the Sun have more and stronger dark spots than in the solar case, including on the rotational pole itself. Doppler imaging, which has so far produced the most detailed images of surface structures on other stars than the Sun, cannot always distinguish the hemisphere in which the starspots are located, especially in the equatorial region and if the data quality is not optimal. This leads to problems in investigating the north-south distribution of starspot active latitudes (those latitudes with more spot activity), which are crucial constraints of dynamo theory. Polar spots, inferred only from Doppler tomography, could plausibly be observational artifacts, casting some doubt on their very existence. Here we report imaging of the old, magnetically-active star $\zeta$ Andromedae using long-baseline infrared interferometry. In our data, a dark polar spot is seen in each of two epochs, while lower-latitude spot structures in both hemispheres do not persist between observations revealing global starspot asymmetries. The north-south symmetry of active latitudes observed on the Sun is absent on $\zeta$ And, which hosts global spot patterns that cannot be produced by solar-type dynamos.

High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41.

Context. Red giant branch (RGB) stars are very bright objects in galaxies and are often used as standard candles. Interferometry is the ideal tool to characterize the dynamics and morphology of their atmospheres. Aims. We aim at precisely characterising the surface dynamics of a sample of RGB stars. Methods. We obtained interferometric observations for three RGB stars with the MIRC instrument mounted at the CHARA interfer- ometer. We looked for asymmetries on the stellar surfaces using limb-darkening models. Results. We measured the apparent diameters of HD197989 (Epsilon Cyg) = 4.61+-0.02 mas, HD189276 (HR7633) = 2.95+-0.01 mas, and HD161096 (Beta Oph) = 4.43+-0.01 mas. We detected departures from the centrosymmetric case for all three stars with the tendency of a greater effect for lower logg of the sample. We explored the causes of this signal and conclude that a possible explanation to the interferometric signal is the convection-related and/or the magnetic-related surface activity. However, it is necessary to monitor these stars with new observations, possibly coupled with spectroscopy, in order to firmly establish the cause.

In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AU above the midplane at a radial distance of 77 AU, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry.

High-mass multiples might form via fragmentation of self-gravitational disks or alternative scenarios such as disk-assisted capture. However, only few observational constraints exist on the architecture and disk structure of high-mass protobinaries and their accretion properties. Here we report the discovery of a close ($57.9\pm0.2$mas=170au) high-mass protobinary, IRAS17216-3801, where our VLTI/GRAVITY+AMBER near-infrared interferometry allows us to image the circumstellar disks around the individual components with 3 milliarcsecond resolution. We estimate the component masses to $\sim20$ and $\sim18 M_{\odot}$ and find that the radial intensity profiles can be reproduced with an irradiated disk model, where the inner regions are excavated of dust, likely tracing the dust sublimation region in these disks. The circumstellar disks are strongly misaligned with respect to the binary separation vector, which indicates that the tidal forces did not have time to realign the disks, pointing towards a young dynamical age of the system. We constrain the distribution of the Br$\gamma$ and CO-emitting gas using VLTI/GRAVITY spectro-interferometry and VLT/CRIRES spectro-astrometry and find that the secondary is accreting at a higher rate than the primary. VLT/NACO imaging shows $L'$-band emission on 3-4 times larger scales than the binary separation, matching the expected dynamical truncation radius for the circumbinary disk. The IRAS17216-3801 system is $\sim3\times$ more massive and $\sim5\times$ more compact than other high-mass multiplies imaged at infrared wavelengths and the first high-mass protobinary system where circumstellar and circumbinary dust disks could be spatially resolved. This opens exciting new opportunities for studying star-disk interactions and the role of multiplicity in high-mass star formation.

Magnetic fields are essential ingredients of many physical processes in the interiors and envelopes of cool stars. Yet their direct detection and characterisation is notoriously difficult, requiring high-quality observations and advanced analysis techniques. Significant progress has been recently achieved by several types of direct magnetic field studies on the surfaces of cool active stars. In particular, complementary techniques of the field topology mapping with polarisation data and total magnetic flux measurements from intensity spectra have been systematically applied to different classes of active stars leading to interesting and occasionally controversial results. In this paper we summarise the current status of direct magnetic field studies of cool stars, and investigations of surface inhomogeneities caused by the field, based on the material presented at the Cool Stars 19 splinter session.

Context. It is now generally accepted that the near-infrared excess of Herbig AeBe stars originates in the dust of a circumstellar disk. Aims. The aims of this article are to infer the radial and vertical structure of these disks at scales of order one au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6micron) PIONIER/VLTI interferometer. The largest baselines allowed us to resolve (at least partially) structures of a few tenths of an au at typical distances of a few hundred parsecs. Dedicated UBVRIJHK photometric measurements were also obtained. Spectral and 2D geometrical parameters are extracted via fits of a few simple models: ellipsoids and broadened rings with azimuthal modulation. Model bias is mitigated by parallel fits of physical disk models. Sample statistics were evaluated against similar statistics for the physical disk models to infer properties of the sample objects as a group. Results. We find that dust at the inner rim of the disk has a sublimation temperature Tsub~1800K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide delta_r>=0.5. A wide ring favors a rim that, on the star-facing side, looks more like a knife edge than a doughnut. The data are also compatible with a the combination of a narrow ring and an inner disk of unspecified nature inside the dust sublimation radius. The disk inner part has a thickness z/r~0.2, flaring to z/r~0.5 in the outer part. We confirm the known luminosity-radius relation; a simple physical model is consistent with both the mean luminosity-radius relation and the ring relative width; however, a significant spread around the mean relation is present. In some of the objects we find a halo component, fully resolved at the shortest interferometer spacing, that is related to the HAeBe class.

The Planet Formation Imager (PFI) is a project for a very large optical interferometer intended to obtain images of the planet formation process at scales as small as the Hill sphere of giant exoplanets. Its main science instruments will work in the thermal infrared but it will be cophased in the near infrared, where it requires also some capacity for scientific imaging. PFI imaging and resolution specifications imply an array of 12 to 20 apertures and baselines up to a few kilometers cophased at near infrared coherent magnitudes as large as 10. This paper discusses various cophasing architectures and the corresponding minimum diameter of individual apertures, which is the dominant element of PFI cost estimates. From a global analysis of the possible combinations of pairwise fringe sensors, we show that conventional approaches used in current interferometers imply the use of prohibitively large telescopes and we indicate the innovative strategies that would allow building PFI with affordable apertures smaller than 2 m in diameter. The approach with the best potential appears to be Hierarchical Fringe Tracking based on "two beams spatial filters" that cophase pairs of neighboring telescopes with all the efficiency of a two telescopes fringe tracker and transmit most of the flux as if it was produced by an unique single mode aperture to cophase pairs of pairs and then pairs of groups of apertures. We consider also the adaptation to PFI of more conventional approaches such as a combination of GRAVITY like fringe trackers or single or multiple chains of 2T fringe trackers.

We present interferometric observations of the sigma Orionis triple system using the CHARA Array, NPOI, and VLTI. Using these measurements, we spatially resolve the orbit of the close spectroscopic binary (Aa,Ab) for the first time and present a revised orbit for the wide pair (A,B). Combining the visual orbits with previously published radial velocity measurements and new radial velocities measured at CTIO, we derive dynamical masses for the three massive stars in the system of M_Aa = 16.99 +/- 0.20 Msun, M_Ab = 12.81 +/- 0.18 Msun, and M_B = 11.5 +/- 1.2 Msun. The inner and outer orbits in the triple are not coplanar, with a relative inclination of 120-127 deg. The orbital parallax provides a precise distance of 387.5 +/- 1.3 pc to the system. This is a significant improvement over previous estimates of the distance to the young sigma Orionis cluster.

Using the high-precision photometry from the Kepler space telescope, thousands of stars with stellar and planetary companions have been observed. The characterization of stars with companions is not always straightforward and can be contaminated by systematic and stellar influences on the light curves. Here, through a detailed analysis of starspots and eclipses, we identify KOI-1003 as a new, active RS CVn star---the first identified with data from Kepler. The Kepler light curve of this close binary system exhibits the system's primary transit, secondary eclipse, and starspot evolution of two persistent active longitudes. The near equality of the system's orbital and rotation periods indicates the orbit and primary star's rotation are nearly synchronized ($P_\mathrm{orb} = 8.360613\pm0.000003$ days; $P_\mathrm{rot} \sim 8.23$ days). By assuming the secondary star is on the main sequence, we suggest the system consists of a $1.45^{+0.11}_{-0.19} \ M_\odot$ subgiant primary and a $0.59^{+0.03}_{-0.04} \ M_\odot$ main-sequence companion. Our work gives a distance of $4400 \pm 600$ pc and an age of $t = 3.0^{-0.5}_{+2.0}$ Gyr, parameters which are discrepant with previous studies that included the star as a member of the open cluster NGC 6791.

The Planet Formation Imager (PFI) Project has formed a Technical Working Group (TWG) to explore possible facility architectures to meet the primary PFI science goal of imaging planet formation in situ in nearby star- forming regions. The goals of being sensitive to dust emission on solar system scales and resolving the Hill-sphere around forming giant planets can best be accomplished through sub-milliarcsecond imaging in the thermal infrared. Exploiting the 8-13 micron atmospheric window, a ground-based long-baseline interferometer with approximately 20 apertures including 10km baselines will have the necessary resolution to image structure down 0.1 milliarcseconds (0.014 AU) for T Tauri disks in Taurus. Even with large telescopes, this array will not have the sensitivity to directly track fringes in the mid-infrared for our prime targets and a fringe tracking system will be necessary in the near-infrared. While a heterodyne architecture using modern mid-IR laser comb technology remains a competitive option (especially for the intriguing 24 and 40{\mu}m atmospheric windows), the prioritization of 3-5{\mu}m observations of CO/H2O vibrotational levels by the PFI-Science Working Group (SWG) pushes the TWG to require vacuum pipe beam transport with potentially cooled optics. We present here a preliminary study of simulated L- and N-band PFI observations of a realistic 4-planet disk simulation, finding 21x2.5m PFI can easily detect the accreting protoplanets in both L and N-band but can see non-accreting planets only in L band. (abridged -- see PDF for full abstract)

The Planet Formation Imager (PFI) project aims to image the period of planet assembly directly, resolving structures as small as a giant planet's Hill sphere. These images will be required in order to determine the key mechanisms for planet formation at the time when processes of grain growth, protoplanet assembly, magnetic fields, disk/planet dynamical interactions and complex radiative transfer all interact - making some planetary systems habitable and others inhospitable. We will present the overall vision for the PFI concept, focusing on the key technologies and requirements that are needed to achieve the science goals. Based on these key requirements, we will define a cost envelope range for the design and highlight where the largest uncertainties lie at this conceptual stage.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

Transitional discs are a class of circumstellar discs around young stars with extensive clearing of dusty material within their inner regions on 10s of au scales. One of the primary candidates for this kind of clearing is the formation of planet(s) within the disc that then accrete or clear their immediate area as they migrate through the disc. Our sample included eight transitional discs. Using the Keck/NIRC2 instrument we utilised the Sparse Aperture Masking (SAM) interferometry technique to search for asymmetries indicative of ongoing planet formation. We searched for close-in companions using both model fitting and interferometric image reconstruction techniques. Using simulated data, we derived diagnostics that helped us to distinguish between point sources and extended asymmetric disc emission. In addition, we investigated the degeneracy between the contrast and separation that appear for marginally resolved companions. We found FP Tau to contain a previously unseen disc wall, and DM Tau, LkHa 330, and TW Hya to contain an asymmetric signal indicative of point source-like emission. We placed upper limits on the contrast of a companion in RXJ1842.9-3532 and V2246 Oph. We ruled the asymmetry signal in RXJ1615.3-3255 and V2062 Oph to be false positives. In the cases where our data indicated a potential companion we computed estimates for the value of $M_c \dot M_c$ and found values in the range of $10^{-5} - 10^{-3} M^2_J yr^{-1}$. We found significant asymmetries in four targets. Of these, three were consistent with companions. We resolved a previously unseen gap in the disc of FP Tau extending inwards from approximately 10 au.

During their formation phase stars gain most of their mass in violent episodic accretion events, such as observed in FU Orionis (FUor) and EXor stars. V346 Normae is a well-studied FUor that underwent a strong outburst beginning in $\sim1980$. Here, we report photometric and spectroscopic observations which show that the visual/near-infrared brightness has decreased dramatically between the 1990s and 2010 (${\Delta}R\approx10.9^{\rm m}$, ${\Delta}J\approx7.8^{\rm m}$, ${\Delta}K\approx5.8^{\rm m}$). The spectral properties of this fading event cannot be explained with variable extinction alone, but indicate a drop in accretion rate by 2-3 orders of magnitude, marking the first time that a member of the FUor class has been observed to switch to a very low accretion phase. Remarkably, in the last few years (2011-2015) V346 Nor has brightened again at all near-infrared wavelengths, indicating the onset of a new outburst event. The observed behaviour might be consistent with the clustered luminosity bursts that have been predicted by recent gravitational instability and fragmentation models for the early stages of protostellar evolution. Given V346 Nor's unique characteristics (concerning outburst duration, repetition frequency, and spectroscopic diagnostics), our results also highlight the need for revisiting the FUor/EXor classification scheme.

We report new CHARA/MIRC interferometric observations of the Cepheid archetype $\delta$ Cep, which aimed at detecting the newly discovered spectroscopic companion. We reached a maximum dynamic range $\Delta H $ = 6.4, 5.8, and 5.2 mag, respectively within the relative distance to the Cepheid $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. Our observations did not show strong evidence of a companion. We have a marginal detection at $3\sigma$ with a flux ratio of 0.21\%, but nothing convincing as we found other possible probable locations. We ruled out the presence of companion with a spectral type earlier than F0V, A1V and B9V, respectively for the previously cited ranges $r$. From our estimated sensitivity limits and the Cepheid light curve, we derived lower-limit magnitudes in the $H$ band for this possible companion to be $H_\mathrm{comp} > 9.15, 8.31$ and 7.77 mag, respectively for $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. We also found that to be consistent with the predicted orbital period, the companion has to be located at a projected separation $< 24$ mas with a spectral type later than a F0V star.

We present near and mid-infrared interferometric observations made with the Keck Interferometer Nuller and near-contemporaneous spectro-photometry from the IRTF of 11 well known young stellar objects, several observed for the first time in these spectral and spatial resolution regimes. With AU-level spatial resolution, we first establish characteristic sizes of the infrared emission using a simple geometrical model consisting of a hot inner rim and mid-infrared disk emission. We find a high degree of correlation between the stellar luminosity and the mid-infrared disk sizes after using near-infrared data to remove the contribution from the inner rim. We then use a semi-analytical physical model to also find that the very widely used "star + inner dust rim + flared disk" class of models strongly fails to reproduce the SED and spatially-resolved mid-infrared data simultaneously; specifically a more compact source of mid-infrared emission is required than results from the standard flared disk model. We explore the viability of a modification to the model whereby a second dust rim containing smaller dust grains is added, and find that the two-rim model leads to significantly improved fits in most cases. This complexity is largely missed when carrying out SED modelling alone, although detailed silicate feature fitting by McClure et al. 2013 recently came to a similar conclusion. As has been suggested recently by Menu et al. 2015, the difficulty in predicting mid-infrared sizes from the SED alone might hint at "transition disk"-like gaps in the inner AU; however, the relatively high correlation found in our mid-infrared disk size vs. stellar luminosity relation favors layered disk morphologies and points to missing disk model ingredients instead.

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 present new mid-infrared interferometric observations of the massive young stellar object IRAS 13481-6124, using VLTI/MIDI for spectrally-resolved, long-baseline measurements (projected baselines up to $\sim120$ m) and GSO/T-ReCS for aperture-masking interferometry in five narrow-band filters (projected baselines of $\sim1.8-6.4$ m) in the wavelength range of $7.5-13$ $\mu$m. We combine these measurements with previously-published interferometric observations in the $K$ and $N$ bands in order to assemble the largest collection of infrared interferometric observations for a massive YSO to date. Using a combination of geometric and radiative-transfer models, we confirm the detection at mid-infrared wavelengths of the disk previously inferred from near-infrared observations. We show that the outflow cavity is also detected at both near- and mid-infrared wavelengths, and in fact dominates the mid-infrared emission in terms of total flux. For the disk, we derive the inner radius ($\sim1.8$ mas or $\sim6.5$ AU at 3.6 kpc), temperature at the inner rim ($\sim1760$ K), inclination ($\sim48$ deg) and position angle ($\sim107$ deg). We determine that the mass of the disk cannot be constrained without high-resolution observations in the (sub-)millimeter regime or observations of the disk kinematics, and could be anywhere from $\sim10^{-3}$ to $20$ M$_\odot$. Finally, we discuss the prospects of interpreting the spectral energy distributions of deeply-embedded massive YSOs, and warn against attempting to infer disk properties from the SED.

At present, the principal limitation on the resolution and contrast of astronomical imaging instruments comes from aberrations in the optical path, which may be imposed by the Earth's turbulent atmosphere or by variations in the alignment and shape of the telescope optics. These errors can be corrected physically, with active and adaptive optics, and in post-processing of the resulting image. A recently-developed adaptive optics post-processing technique, called kernel phase interferometry, uses linear combinations of phases that are self-calibrating with respect to small errors, with the goal of constructing observables that are robust against the residual optical aberrations in otherwise well-corrected imaging systems. Here we present a direct comparison between kernel phase and the more established competing techniques, aperture masking interferometry, point spread function (PSF) fitting and bispectral analysis. We resolve the alpha Ophiuchi binary system near periastron, using the Palomar 200-Inch Telescope. This is the first case in which kernel phase has been used with a full aperture to resolve a system close to the diffraction limit with ground-based extreme adaptive optics observations. Excellent agreement in astrometric quantities is found between kernel phase and masking, and kernel phase significantly outperforms PSF fitting and bispectral analysis, demonstrating its viability as an alternative to conventional non-redundant masking under appropriate conditions.

High angular resolution images of IRC+10216 are presented in several near infrared wavelengths spanning more than 8 years. These maps have been reconstructed from interferometric observations obtained at both Keck and the VLT, and also from stellar occultations by the rings of Saturn observed with the Cassini spacecraft. The dynamic inner regions of the circumstellar environment are monitored over eight epochs ranging between January 2000 and July 2008. The system is shown to experience substantial evolution within this period including the fading of many previously reported persistent features, some of which had been identified as the stellar photosphere. These changes are discussed in context of existing models for the nature of the underlying star and the circumstellar environment. With access to these new images, we are able to report that none of the previously identified bright spots in fact contain the star, which is buried in its own dust and not directly visible in the near infrared.

Presented are the first interferometric images of cool starspots on the chromospherically active giant $\lambda$ Andromedae. These images represent the first model-independent images of cool starspots on a star other than the Sun to date. The interferometric observations, taken with the Michigan Infra-Red Combiner coupled to the Center for High Angular Resolution Astronomy Array, span 26 days from Aug 17$^{th}$, 2008 to Sep 24$^{th}$, 2011. The photometric time series acquired at Fairborn Observatory spanning Sep 20$^{th}$, 2008 to Jan 20$^{th}$, 2011 is also presented. The angular diameter and power law limb-darkening coefficient of this star are 2.759 $\pm$ 0.050 mas and 0.229 $\pm$ 0.111, respectively. Starspot properties are obtained from both modeled and SQUEEZE reconstructed images. The images from 2010 through 2011 show anywhere from one to four starspots. The measured properties of identical starspots identified in both the model and reconstructed images are within two $\sigma$ error bars in 51$\%$ of cases. The cadence in the data for the 2010 and 2011 data sets are sufficient to measure a stellar rotation period based on apparent starspot motion. This leads to estimates of the rotation period (P$_{2010}$ = 60 $\pm$ 13 days, P$_{2011}$ = 54.0 $\pm$ 7.6 days) that are consistent with the photometrically determined period of 54.8 days. In addition, the inclination and position angle of the rotation axis is computed for both the 2010 and 2011 data sets; values ($\bar{\Psi}$ = 21.5$\degree$, $\bar{\emph{i}}$ = 78.0$\degree$) for each are nearly identical between the two years. \end{abstract}

We report on a total of 106 nights of optical interferometric observations of the $\epsilon$ Aurigae system taken during the last 14 years by four beam combiners at three different interferometric facilities. This long sequence of data provides an ideal assessment of the system prior to, during, and after the recent 2009-2011 eclipse. We have reconstructed model-independent images from the 10 in-eclipse epochs which show that a disk-like object is indeed responsible for the eclipse. Using new 3D, time-dependent modeling software, we derive the properties of the F-star (diameter, limb darkening), determine previously unknown orbital elements ($\Omega$, $i$), and access the global structures of the optically thick portion of the eclipsing disk using both geometric models and approximations of astrophysically relevant density distributions. These models may be useful in future hydrodynamical modeling of the system. Lastly, we address several outstanding research questions including mid-eclipse brightening, possible shrinking of the F-type primary, and any warps or sub-features within the disk.

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.

Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3{\sigma} detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8{\sigma} using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.

A classical nova occurs when material accreting onto the surface of a white dwarf in a close binary system ignites in a thermonuclear runaway. Complex structures observed in the ejecta at late stages could result from interactions with the companion during the common envelope phase. Alternatively, the explosion could be intrinsically bipolar, resulting from a localized ignition on the surface of the white dwarf or as a consequence of rotational distortion. Studying the structure of novae during the earliest phases is challenging because of the high spatial resolution needed to measure their small sizes. Here we report near-infrared interferometric measurements of the angular size of Nova Delphini 2013, starting from one day after the explosion and continuing with extensive time coverage during the first 43 days. Changes in the apparent expansion rate can be explained by an explosion model consisting of an optically thick core surrounded by a diffuse envelope. The optical depth of the ejected material changes as it expands. We detect an ellipticity in the light distribution, suggesting a prolate or bipolar structure that develops as early as the second day. Combining the angular expansion rate with radial velocity measurements, we derive a geometric distance to the nova of 4.54 +/- 0.59 kpc from the Sun.

To measure the properties of both components of the RS CVn binary sigma Geminorum (sigma Gem), we directly detect the faint companion, measure the orbit, obtain model-independent masses and evolutionary histories, detect ellipsoidal variations of the primary caused by the gravity of the companion, and measure gravity darkening. We detect the companion with interferometric observations obtained with the Michigan InfraRed Combiner (MIRC) at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array with a primary-to-secondary H-band flux ratio of 270+/-70. A radial velocity curve of the companion was obtained with spectra from the Tillinghast Reflector Echelle Spectrograph (TRES) on the 1.5-m Tillinghast Reflector at Fred Lawrence Whipple Observatory (FLWO). We additionally use new observations from the Tennessee State University Automated Spectroscopic and Photometric Telescopes (AST and APT, respectively). From our orbit, we determine model-independent masses of the components (M_1 = 1.28 +/- 0.07 M_Sun, M_2 = 0.73 +/- 0.03 M_Sun), and estimate a system age of 5 -/+ 1 Gyr. An average of the 27-year APT light curve of sigma Gem folded over the orbital period (P = 19.6027 +/- 0.0005 days) reveals a quasi-sinusoidal signature, which has previously been attributed to active longitudes 180 deg apart on the surface of sigma Gem. With the component masses, diameters, and orbit, we find that the predicted light curve for ellipsoidal variations due to the primary star partially filling its Roche lobe potential matches well with the observed average light curve, offering a compelling alternative explanation to the active longitudes hypothesis. Measuring gravity darkening from the light curve gives beta < 0.1, a value slightly lower than that expected from recent theory.

The rapidly rotating Be star phi Persei was spun up by mass and angular momentum transfer from a now stripped-down, hot subdwarf companion. Here we present the first high angular resolution images of phi Persei made possible by new capabilities in longbaseline interferometry at near-IR and visible wavelengths. We observed phi Persei with the MIRC and VEGA instruments of the CHARA Array. Additional MIRC-only observations were performed to track the orbital motion of the companion, and these were fit together with new and existing radial velocity measurements of both stars to derive the complete orbital elements and distance. The hot subdwarf companion is clearly detected in the near-IR data at each epoch of observation with a flux contribution of 1.5% in the H band, and restricted fits indicate that its flux contribution rises to 3.3% in the visible. A new binary orbital solution is determined by combining the astrometric and radial velocity measurements. The derived stellar masses are 9.6+-0.3Msol and 1.2+-0.2Msol for the Be primary and subdwarf secondary, respectively. The inferred distance (186 +- 3 pc), kinematical properties, and evolutionary state are consistent with membership of phi Persei in the alpha Per cluster. From the cluster age we deduce significant constraints on the initial masses and evolutionary mass transfer processes that transformed the phi Persei binary system. The interferometric data place strong constraints on the Be disk elongation, orientation, and kinematics, and the disk angular momentum vector is coaligned with and has the same sense of rotation as the orbital angular momentum vector. The VEGA visible continuum data indicate an elongated shape for the Be star itself, due to the combined effects of rapid rotation, partial obscuration of the photosphere by the circumstellar disk, and flux from the bright inner disk.

Among the most fascinating and hotly-debated areas in contemporary astrophysics are the means by which planetary systems are assembled from the large rotating disks of gas and dust which attend a stellar birth. Although important work has already been, and is still being done both in theory and observation, a full understanding of the physics of planet formation can only be achieved by opening observational windows able to directly witness the process in action. The key requirement is then to probe planet-forming systems at the natural spatial scales over which material is being assembled. By definition, this is the so-called Hill Sphere which delineates the region of influence of a gravitating body within its surrounding environment. The Planet Formation Imager project (PFI) has crystallized around this challenging goal: to deliver resolved images of Hill-Sphere-sized structures within candidate planet-hosting disks in the nearest star-forming regions. In this contribution we outline the primary science case of PFI. For this purpose, we briefly review our knowledge about the planet-formation process and discuss recent observational results that have been obtained on the class of transition disks. Spectro-photometric and multi-wavelength interferometric studies of these systems revealed the presence of extended gaps and complex density inhomogeneities that might be triggered by orbiting planets. We present detailed 3-D radiation-hydrodynamic simulations of disks with single and multiple embedded planets, from which we compute synthetic images at near-infrared, mid-infrared, far-infrared, and sub-millimeter wavelengths, enabling a direct comparison of the signatures that are detectable with PFI and complementary facilities such as ALMA. From these simulations, we derive some preliminary specifications that will guide the array design and technology roadmap of the facility.

The Planet Formation Imager (PFI) is a future world facility that will image the process of planetary formation. It will have an angular resolution and sensitivity sufficient to resolve sub-Hill sphere structures around newly formed giant planets orbiting solar-type stars in nearby star formation regions. We present one concept for this design consisting of twenty-seven or more 4m telescopes with kilometric baselines feeding a mid-infrared spectrograph where starlight is mixed with a frequency-comb laser. Fringe tracking will be undertaken in H-band using a fiber-fed direct detection interferometer, meaning that all beam transport is done by communications band fibers. Although heterodyne interferometry typically has lower signal-to-noise than direct detection interferometry, it has an advantage for imaging fields of view with many resolution elements, because the signal in direct detection has to be split many ways while the signal in heterodyne interferometry can be amplified prior to combining every baseline pair. We compare the performance and cost envelope of this design to a comparable direct-detection design.

Complex non-linear and dynamic processes lie at the heart of the planet formation process. Through numerical simulation and basic observational constraints, the basics of planet formation are now coming into focus. High resolution imaging at a range of wavelengths will give us a glimpse into the past of our own solar system and enable a robust theoretical framework for predicting planetary system architectures around a range of stars surrounded by disks with a diversity of initial conditions. Only long-baseline interferometry can provide the needed angular resolution and wavelength coverage to reach these goals and from here we launch our planning efforts. The aim of the "Planet Formation Imager" (PFI) project is to develop the roadmap for the construction of a new near-/mid-infrared interferometric facility that will be optimized to unmask all the major stages of planet formation, from initial dust coagulation, gap formation, evolution of transition disks, mass accretion onto planetary embryos, and eventual disk dispersal. PFI will be able to detect the emission of the cooling, newly-formed planets themselves over the first 100 Myrs, opening up both spectral investigations and also providing a vibrant look into the early dynamical histories of planetary architectures. Here we introduce the Planet Formation Imager (PFI) Project (www.planetformationimager.org) and give initial thoughts on possible facility architectures and technical advances that will be needed to meet the challenging top-level science requirements.

One of the primary goals of exoplanet science is to find and characterize habitable planets, and direct imaging will play a key role in this effort. Though imaging a true Earth analog is likely out of reach from the ground, the coming generation of giant telescopes will find and characterize many planets in and near the habitable zones (HZs) of nearby stars. Radial velocity and transit searches indicate that such planets are common, but imaging them will require achieving extreme contrasts at very small angular separations, posing many challenges for adaptive optics (AO) system design. Giant planets in the HZ may even be within reach with the latest generation of high-contrast imagers for a handful of very nearby stars. Here we will review the definition of the HZ, and the characteristics of detectable planets there. We then review some of the ways that direct imaging in the HZ will be different from the typical exoplanet imaging survey today. Finally, we present preliminary results from our observations of the HZ of {\alpha} Centauri A with the Magellan AO system's VisAO and Clio2 cameras.

We aim at determining the masses of Cepheids in binary systems, as well as their geometric distances and the flux contribution of the companions. The combination of interferometry with spectroscopy will offer a unique and independent estimate of the Cepheid masses. Using long-baseline interferometry at visible and infrared wavelengths, it is possible to spatially resolve binary systems containing a Cepheid down to milliarcsecond separations. Based on the resulting visual orbit and radial velocities, we can then derive the fundamental parameters of these systems, particularly the masses of the components and the geometric distance. We therefore performed interferometric observations of the first-overtone mode Cepheid V1334 Cyg with the CHARA/MIRC combiner. We report the first detection of a Cepheid companion using long-baseline interferometry. We detect the signature of a companion orbiting V1334 Cyg at two epochs. We measure a flux ratio between the companion and the Cepheid f = 3.10+/-0.08%, giving an apparent magnitude mH = 8.47+/-0.15mag. The combination of interferometric and spectroscopic data have enabled the unique determination of the orbital elements: P = 1938.6+/-1.2 days, Tp = 2 443 616.1+/-7.3, a = 8.54+/-0.51mas, i = 124.7+/-1.8{\deg}, e = 0.190+/-0.013, {\omega} = 228.7+/-1.6{\deg}, and {\Omega} = 206.3+/-9.4{\deg}. We derive a minimal distance d ~ 691 pc, a minimum mass for both stars of 3.6 Msol, with a spectral type earlier than B5.5V for the companion star. Our measured flux ratio suggests that radial velocity detection of the companion using spectroscopy is within reach, and would provide an orbital parallax and model-free masses.

Optical interferometry is the only technique giving access to milli-arcsecond (mas) resolution at infrared wavelengths. For Cepheids, this is a powerful and unique tool to detect the orbiting companions and the circumstellar envelopes (CSE). CSEs are interesting because they might be used to trace the Cepheid evolution history, and more particularly they could impact the distance scale. Cepheids belonging to binary systems offer an unique opportunity to make progress in resolving the Cepheid mass discrepancy. The combination of spectroscopic and interferometric measurements will allow us to derive the orbital elements, distances, and dynamical masses. Here we focus on recent results using 2- to 6-telescopes beam combiners for the Cepheids X Sgr, T Mon and V1334 Cyg.

Classical Cepheid stars have been considered since more than a century as reliable tools to estimate distances in the universe thanks to their Period-Luminosity (P-L) relationship. Moreover, they are also powerful astrophysical laboratories, providing fundamental clues for studying the pulsation and evolution of intermediate-mass stars. When in binary systems, we can investigate the age and evolution of the Cepheid, estimate the mass and distance, and constrain theoretical models. However, most of the companions are located too close to the Cepheid (1-40 mas) to be spatially resolved with a 10-meter class telescope. The only way to spatially resolve such systems is to use long-baseline interferometry. Recently, we have started a unique and long-term interferometric program that aims at detecting and characterizing physical parameters of the Cepheid companions, with as main objectives the determination of accurate masses and geometric distances.

Two red supergiants of the Per OB1 association, RS Per and T Per, have been observed in H band using the MIRC instrument at the CHARA array. The data show clear evidence of departure from circular symmetry. We present here new techniques specially developed to analyze such cases, based on state-of-the-art statistical frameworks. The stellar surfaces are first modeled as limb-darkened discs based on SATLAS models that fit both MIRC interferometric data and publicly available spectrophotometric data. Bayesian model selection is then used to determine the most probable number of spots. The effective surface temperatures are also determined and give further support to the recently derived hotter temperature scales of red su- pergiants. The stellar surfaces are reconstructed by our model-independent imaging code SQUEEZE, making use of its novel regularizer based on Compressed Sensing theory. We find excellent agreement between the model-selection results and the reconstructions. Our results provide evidence for the presence of near-infrared spots representing about 3-5% of the stellar flux.

The Event Horizon Telescope (EHT) is a project to assemble a Very Long Baseline Interferometry (VLBI) network of mm wavelength dishes that can resolve strong field General Relativistic signatures near a supermassive black hole. As planned, the EHT will include enough dishes to enable imaging of the predicted black hole "shadow", a feature caused by severe light bending at the black hole boundary. The center of M87, a giant elliptical galaxy, presents one of the most interesting EHT targets as it exhibits a relativistic jet, offering the additional possibility of studying jet genesis on Schwarzschild radius scales. Fully relativistic models of the M87 jet that fit all existing observational constraints now allow horizon-scale images to be generated. We perform realistic VLBI simulations of M87 model images to examine detectability of the black shadow with the EHT, focusing on a sequence of model images with a changing jet mass load radius. When the jet is launched close to the black hole, the shadow is clearly visible both at 230 and 345 GHz. The EHT array with a resolution of 20-30$\mu$as resolution ($\sim$2-4 Schwarzschild radii) is able to image this feature independent of any theoretical models and we show that imaging methods used to process data from optical interferometers are applicable and effective for EHT data sets. We demonstrate that the EHT is also capable of tracing real-time structural changes on a few Schwarzschild radii scales, such as those implicated by VHE flaring activity of M87. While inclusion of ALMA in the EHT is critical for shadow imaging, generally the array is robust against loss of a station.

The emergence of optical interferometers with three and more telescopes allows image reconstruction of astronomical objects at the milliarcsecond scale. However, some objects contain components with very different spectral energy distributions (SED; i.e. different temperatures), which produces strong chromatic effects on the interferograms that have to be managed with care by image reconstruction algorithms. For example, the gray approximation for the image reconstruction process results in a degraded image if the total (u, v)-coverage given by the spectral supersynthesis is used. The relative flux contribution of the central object and an extended structure changes with wavelength for different temperatures. For young stellar objects, the known characteristics of the central object (i.e., stellar SED), or even the fit of the spectral index and the relative flux ratio, can be used to model the central star while reconstructing the image of the extended structure separately. Methods. We present a new method, called SPARCO (semi-parametric algorithm for the image reconstruction of chromatic objects), which describes the spectral characteristics of both the central object and the extended structure to consider them properly when reconstructing the image of the surrounding environment. We adapted two image-reconstruction codes (Macim, Squeeze, and MiRA) to implement this new prescription. SPARCO is applied using Macim, Squeeze and MiRA on a young stellar object model and also on literature data on HR5999 in the near-infrared with the VLTI. This method paves the way to improved aperture-synthesis imaging of several young stellar objects with existing datasets. More generally, the approach can be used on astrophysical sources with similar features such as active galactic nuclei, planetary nebulae, and asymptotic giant branch stars.

Imaging the direct light signal from a faint exoplanet against the overwhelming glare of its host star presents one of the fundamental challenges to modern astronomical instrumentation. Achieving sufficient signal-to-noise for detection by direct imaging is limited by three basic physical processes: aberration of the wavefronts (both instrumental and atmospheric), photon noise, and detector noise. In this paper, we advance a novel optical setup which synthesizes the advantages of two different techniques: nulling interferometry to mitigate photon noise, and closure phase to combat optical aberrations. Our design, which employs technology from integrated optics and photonics, is intended to combine the advantageous aspects of both a coronagraph and a non-redundant interferometer inside a single optical device. We show that such an instrument would have a dynamic range limited either by i) the readout noise (if perfect co-phasing), or ii) the photon noise due to stellar flux leakage (in the case of imperfect nulling). This concept is optimal when the readout noise is not the main limitation, ie, for space interferometry or for ground based observations of bright stellar hosts (apparent magnitude brighter than 10).

Cophasing six telescopes from the CHARA array, the CHARA-Michigan Phasetracker (CHAMP) and Michigan Infrared Combiner (MIRC) are pushing the frontiers of infrared long-baseline interferometric imaging in key scientific areas such as star- and planet-formation. Here we review our concepts and recent improvements on the CHAMP and MIRC control interfaces, which establish the communication to the real-time data recording & fringe tracking code, provide essential performance diagnostics, and assist the observer in the alignment and flux optimization procedure. For fringe detection and tracking with MIRC, we have developed a novel matrix approach, which provides predictions for the fringe positions based on cross-fringe information.

Aims: We aim at detecting and characterizing the main-sequence companion of the Cepheid AX Cir ($P_\mathrm{orb} \sim $ 18 yrs). The long-term objective is to estimate the mass of both components and the distance to the system. Methods: We used the PIONIER combiner at the VLT Interferometer to obtain the first interferometric measurements of the short-period Cepheid AX Cir and its orbiting component. Results: The companion is resolved by PIONIER at a projected separation $\rho = 29.2 \pm 0.2$ mas and projection angle $PA = 167.6 \pm 0.3^{\circ}$. We measured $H$-band flux ratios between the companion and the Cepheid of $0.90 \pm 0.10$ % and $0.75 \pm 0.17$ %, respectively at a pulsation phase for the Cepheid $\phi = 0.24$ and 0.48. The lower contrast at $\phi = 0.48$ is due to increased brightness of the Cepheid compared to the $\phi = 0.24$. This gives an average apparent magnitude $m\mathrm{_H (comp)} = 9.06 \pm 0.24$ mag. The limb-darkened angular diameter of the Cepheid at the two pulsation phases was measured to be $\theta_\mathrm{LD} = 0.839 \pm 0.023$ mas and $\theta_\mathrm{LD} = 0.742 \pm 0.020$ mas, respectively at $\phi = 0.24$ and 0.48. A lower limit on the total mass of the system was also derived based on our measured separation, we found $M_\mathrm{T} \geq 9.7 \pm 0.6 M_\odot$.

Binary post-AGB stars are interesting laboratories to study both the evolution of binaries as well as the structure of circumstellar disks. A multiwavelength high angular resolution study of the prototypical object 89 Herculis is performed with the aim of identifying and locating the different emission components seen in the SED. A large interferometric data set, collected over the past decade and covering optical and near-IR wavelengths, is analyzed with simple geometric models. Combining the interferometric constraints with the photometry and the optical spectra, we reassess the energy budget of the post-AGB star and its circumstellar environment. We report the first (direct) detection of a large (35-40%) optical circumstellar flux contribution and spatially resolve its emission region. Given this large amount of reprocessed and/or redistributed optical light, the fitted size of the emission region is rather compact and fits with(in) the inner rim of the circumbinary dust disk. This rim dominates our K band data through thermal emission and is rather compact, emitting significantly already at a radius of twice the orbital separation. We interpret the circumstellar optical flux as due to a scattering process, with the scatterers located in the extremely puffed-up inner rim of the disk and possibly also in a bipolar outflow seen pole-on. A non-LTE gaseous origin in an inner disk cannot be excluded but is considered highly unlikely. This direct detection of a significant amount of circumbinary light at optical wavelengths poses several significant questions regarding our understanding of both post-AGB binaries and the physics in their circumbinary disks. Although the identification of the source of emission/scattering remains inconclusive without further study on this and similar objects, the implications are manifold.

(Abridged) Dust is expected to be ubiquitous in extrasolar planetary systems owing to the dynamical activity of minor bodies. Inner dust populations are, however, still poorly known because of the high contrast and small angular separation with respect to their host star. We aim to determine the level of near-infrared exozodiacal dust emission around a sample of 42 nearby main sequence stars with spectral types ranging from A to K and to investigate its correlation with various stellar parameters and with the presence of cold dust belts. We use high-precision K-band visibilities obtained with the FLUOR interferometer on the shortest baseline of the CHARA array. The calibrated visibilities are compared with the expected visibility of the stellar photosphere to assess whether there is an additional, fully resolved circumstellar emission. Near-infrared circumstellar emission amounting to about 1% of the stellar flux is detected around 13 of our 42 target stars. Follow-up observations showed that one of them (eps Cep) is associated with a stellar companion, while another one was detected around what turned out to be a giant star (kap CrB). The remaining 11 excesses found around single main sequence stars are most probably associated with hot circumstellar dust, yielding an overall occurrence rate of 28+8-6% for our (biased) sample. We show that the occurrence rate of bright exozodiacal discs correlates with spectral type, K-band excesses being more frequent around A-type stars. It also correlates with the presence of detectable far-infrared excess emission in the case of solar-type stars. This study provides new insight into the phenomenon of bright exozodiacal discs, showing that hot dust populations are probably linked to outer dust reservoirs in the case of solar-type stars. For A-type stars, no clear conclusion can be made regarding the origin of the detected near-infrared excesses.

We present interferometric observations of 7 main-sequence and 3 giant stars with spectral types from B2 to F6 using the PAVO beam combiner at the CHARA array. We have directly determined the angular diameters for these objects with an average precision of 2.3%. We have also computed bolometric fluxes using available photometry in the visible and infrared wavelengths, as well as space-based ultraviolet spectroscopy. Combined with precise \textit{Hipparcos} parallaxes, we have derived a set of fundamental stellar properties including linear radius, luminosity and effective temperature. Fitting the latter to computed isochrone models, we have inferred masses and ages of the stars. The effective temperatures obtained are in good agreement (at a 3% level) with nearly-independent temperature estimations from spectroscopy. They validate recent sixth-order polynomial (B-V)-$T_\mathrm{eff}$ empirical relations \citep{Boyajian2012a}, but suggest that a more conservative third-order solution \citep{vanBelle2009} could adequately describe the (V-K)-$T_\mathrm{eff}$ relation for main-sequence stars of spectral type A0 and later. Finally, we have compared mass values obtained combining surface gravity with inferred stellar radius (\textit{gravity mass}) and as a result of the comparison of computed luminosity and temperature values with stellar evolutionary models (\textit{isochrone mass}). The strong discrepancy between isochrone and gravity mass obtained for one of the observed stars, $\gamma$\,Lyr, suggests that determination of the stellar atmosphere parameters should be revised.

We present the results of long-baseline optical interferometry observations using the Precision Astronomical Visual Observations (PAVO) beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array to measure the angular sizes of three bright Kepler stars: {\theta} Cygni, and both components of the binary system 16 Cygni. Supporting infrared observations were made with the Michigan Infrared Combiner (MIRC) and Classic beam combiner, also at the CHARA Array. We find limb-darkened angular diameters of 0.753+/-0.009 mas for {\theta} Cyg, 0.539+/-0.007 mas for 16 Cyg A and 0.490+/-0.006 mas for 16 Cyg B. The Kepler Mission has observed these stars with outstanding photometric precision, revealing the presence of solar-like oscillations. Due to the brightness of these stars the oscillations have exceptional signal-to-noise, allowing for detailed study through asteroseismology, and are well constrained by other observations. We have combined our interferometric diameters with Hipparcos parallaxes, spectrophotometric bolometric fluxes and the asteroseismic large frequency separation to measure linear radii ({\theta} Cyg: 1.48+/-0.02 Rsun, 16 Cyg A: 1.22+/-0.02 Rsun, 16 Cyg B: 1.12+/-0.02 Rsun), effective temperatures ({\theta} Cyg: 6749+/-44 K, 16 Cyg A: 5839+/-42 K, 16 Cyg B: 5809+/-39 K), and masses ({\theta} Cyg: 1.37+/-0.04 Msun, 16 Cyg A: 1.07+/-0.05 Msun, 16 Cyg B: 1.05+/-0.04 Msun) for each star with very little model dependence. The measurements presented here will provide strong constraints for future stellar modelling efforts.

Pre-transitional disks are protoplanetary disks with a gapped disk structure, potentially indicating the presence of young planets in these systems. In order to explore the structure of these objects and their gap-opening mechanism, we observed the pre-transitional disk V1247 Orionis using the Very Large Telescope Interferometer, the Keck Interferometer, Keck-II, Gemini South, and IRTF. This allows us spatially resolve the AU-scale disk structure from near- to mid-infrared wavelengths (1.5 to 13 {\mu}m), tracing material at different temperatures and over a wide range of stellocentric radii. Our observations reveal a narrow, optically-thick inner-disk component (located at 0.18 AU from the star) that is separated from the optically thick outer disk (radii >46 AU), providing unambiguous evidence for the existence of a gap in this pre-transitional disk. Surprisingly, we find that the gap region is filled with significant amounts of optically thin material with a carbon-dominated dust mineralogy. The presence of this optically thin gap material cannot be deduced solely from the spectral energy distribution, yet it is the dominant contributor at mid-infrared wavelengths. Furthermore, using Keck/NIRC2 aperture masking observations in the H, K', and L' band, we detect asymmetries in the brightness distribution on scales of about 15-40 AU, i.e. within the gap region. The detected asymmetries are highly significant, yet their amplitude and direction changes with wavelength, which is not consistent with a companion interpretation but indicates an inhomogeneous distribution of the gap material. We interpret this as strong evidence for the presence of complex density structures, possibly reflecting the dynamical interaction of the disk material with sub-stellar mass bodies that are responsible for the gap clearing.

We present the first high angular resolution observations in the nearinfrared H-band (1.6 microns) of the Luminous Blue Variable star P Cygni. We obtained six-telescope interferometric observations with the CHARA Array and the MIRC beam combiner. These show that the spatial flux distribution is larger than expected for the stellar photosphere. A two component model for the star (uniform disk) plus a halo (two-dimensional Gaussian) yields an excellent fit of the observations, and we suggest that the halo corresponds to flux emitted from the base of the stellar wind. This wind component contributes about 45% of the H-band flux and has an angular FWHM = 0.96 mas, compared to the predicted stellar diameter of 0.41 mas. We show several images reconstructed from the interferometric visibilities and closure phases, and they indicate a generally spherical geometry for the wind. We also obtained near-infrared spectrophotometry of P Cygni from which we derive the flux excess compared to a purely photospheric spectral energy distribution. The H-band flux excess matches that from the wind flux fraction derived from the two component fits to the interferometry. We find evidence of significant near-infrared flux variability over the period from 2006 to 2010 that appears similar to the variations in the H-alpha emission flux from the wind. Future interferometric observations may be capable of recording the spatial variations associated with temporal changes in the wind structure.

The Kepler target KIC 5110407, a K-type star, shows strong quasi-periodic light curve fluctuations likely arising from the formation and decay of spots on the stellar surface rotating with a period of 3.4693 days. Using an established light-curve inversion algorithm, we study the evolution of the surface features based on Kepler space telescope light curves over a period of two years (with a gap of .25 years). At virtually all epochs, we detect at least one large spot group on the surface causing a 1-10% flux modulation in the Kepler passband. By identifying and tracking spot groups over a range of inferred latitudes, we measured the surface differential rotation to be much smaller than that found for the Sun. We also searched for a correlation between the seventeen stellar flares that occurred during our observations and the orientation of the dominant surface spot at the time of each flare. No statistically-significant correlation was found except perhaps for the very brightest flares, suggesting most flares are associated with regions devoid of spots or spots too small to be clearly discerned using our reconstruction technique. While we may see hints of long-term changes in the spot characteristics and flare statistics within our current dataset, a longer baseline of observation will be needed to detect the existence of a magnetic cycle in KIC 5110407.

We report low resolution near infrared spectroscopic observations of the eruptive star FU Orionis using the Integral Field Spectrograph Project 1640 installed at the Palomar Hale telescope. This work focuses on elucidating the nature of the faint source, located 0.5" south of FU Ori, and identified in 2003 as FU Ori S. We first use our observations in conjunction with published data to demonstrate that the two stars are indeed physically associated and form a true binary pair. We then proceed to extract J and H band spectro-photometry using the damped LOCI algorithm, a reduction method tailored for high contrast science with IFS. This is the first communication reporting the high accuracy of this technique, pioneered by the Project 1640 team, on a faint astronomical source. We use our low resolution near infrared spectrum in conjunction with 10.2 micron interferometric data to constrain the infrared excess of FU Ori S. We then focus on estimating the bulk physical properties of FU Ori S. Our models lead to estimates of an object heavily reddened, A_V =8-12, with an effective temperature of ~ 4000-6500 K . Finally we put these results in the context of the FU Ori N-S system and argue that our analysis provides evidence that FU Ori S might be the more massive component of this binary system

Optical and infrared interferometers definitively established that the photometric standard Vega (alpha Lyrae) is a rapidly rotating star viewed nearly pole-on. Recent independent spectroscopic analyses could not reconcile the inferred inclination angle with the observed line profiles, preferring a larger inclination. In order to resolve this controversy, we observed Vega using the six-beam Michigan Infrared Combiner on the Center for High Angular Resolution Astronomy Array. With our greater angular resolution and dense (u,v)-coverage, we find Vega is rotating less rapidly and with a smaller gravity darkening coefficient than previous interferometric results. Our models are compatible with low photospheric macroturbulence and also consistent with the possible rotational period of ~0.71 days recently reported based on magnetic field observations. Our updated evolutionary analysis explicitly incorporates rapid rotation, finding Vega to have a mass of 2.15+0.10_-0.15 Msun and an age 700-75+150 Myrs, substantially older than previous estimates with errors dominated by lingering metallicity uncertainties (Z=0.006+0.003-0.002).

ASTRA (ASTrometric and phase-Referencing Astronomy) is an upgrade to the existing Keck Interferometer which aims at providing new self-phase referencing (high spectral resolution observation of YSOs), dual-field phase referencing (sensitive AGN observations), and astrometric (known exoplanetary systems characterization and galactic center general relativity in strong field regime) capabilities. With the first high spectral resolution mode now offered to the community, this contribution focuses on the progress of the dual field and astrometric modes.

While the importance of dusty asymptotic giant branch (AGB) stars to galactic chemical enrichment is widely recognised, a sophisticated understanding of the dust formation and wind-driving mechanisms has proven elusive due in part to the difficulty in spatially-resolving the dust formation regions themselves. We have observed twenty dust-enshrouded AGB stars as part of the Keck Aperture Masking Experiment, resolving all of them in multiple near-infrared bands between 1.5 microns and 3.1 microns. We find 45% of the targets to show measurable elongations that, when correcting for the greater distances of the targets, would correspond to significantly asymmetric dust shells on par with the well-known cases of IRC+10216 or CIT6. Using radiative transfer models, we find the sublimation temperature of 1130 +- 90 K and 1170 +- 60 K for silicates and amorphous carbon respectively, both somewhat lower than expected from laboratory measurements and vastly below temperatures inferred from the inner edge of YSO disks. The fact that O-rich and C-rich dust types showed the same sublimation temperature was surprising as well. For the most optically-thick shells (tau > 2 at 2.2 microns), the temperature profile of the inner dust shell is observed to change substantially, an effect we suggest could arise when individual dust clumps become optically-thick at the highest mass-loss rates.

We present the results of the fifth Interferometric Imaging Beauty Contest. The contest consists in blind imaging of test data sets derived from model sources and distributed in the OIFITS format. Two scenarios of imaging with CHARA/MIRC-6T were offered for reconstruction: imaging a T Tauri disc and imaging a spotted red supergiant. There were eight different teams competing this time: Monnier with the software package MACIM; Hofmann, Schertl and Weigelt with IRS; Thi\'ebaut and Soulez with MiRA ; Young with BSMEM; Mary and Vannier with MIROIRS; Millour and Vannier with independent BSMEM and MiRA entries; Rengaswamy with an original method; and Elias with the radio-astronomy package CASA. The contest model images, the data delivered to the contestants and the rules are described as well as the results of the image reconstruction obtained by each method. These results are discussed as well as the strengths and limitations of each algorithm.

Ground-based long baseline interferometers have long been limited in sensitivity by the short integration periods imposed by atmospheric turbulence. The first observation fainter than this limit was performed on January 22, 2011 when the Keck Interferometer observed a K=11.5 target, about one magnitude fainter than its K=10.3 limit. This observation was made possible by the Dual Field Phase Referencing instrument of the ASTRA project: simultaneously measuring the real-time effects of the atmosphere on a nearby bright guide star, and correcting for it on the faint target, integration time longer than the turbulence time scale are made possible. As a prelude to this demonstration, we first present the implementation of Dual Field Phase Referencing on the interferometer. We then detail its on-sky performance focusing on the accuracy of the turbulence correction, and on the resulting fringe contrast stability. We conclude with a presentation of early results obtained with Laser Guide Star AO and the interferometer.

We present simulations of Keck Interferometer ASTRA and VLTI GRAVITY observations of mock star fields in orbit within ~50 milliarcseconds of Sgr A*. Dual-field phase referencing techniques, as implemented on ASTRA and planned for GRAVITY, will provide the sensitivity to observe Sgr A* with infrared interferometers. Our results show an improvement in the confusion noise limit over current astrometric surveys, opening a window to study stellar sources in the region. Since the Keck Interferometer has only a single baseline, the improvement in the confusion limit depends on source position angles. The GRAVITY instrument will yield a more compact and symmetric PSF, providing an improvement in confusion noise which will not depend as strongly on position angle. Our Keck results show the ability to characterize the star field as containing zero, few, or many bright stellar sources. We are also able to detect and track a source down to mK~18 through the least confused regions of our field of view at a precision of ~200 microarcseconds along the baseline direction. This level of precision improves with source brightness. Our GRAVITY results show the potential to detect and track multiple sources in the field. GRAVITY will perform ~10 microarcsecond astrometry on a mK=16.3 source and ~200 microarcsecond astrometry on a mK=18.8 source in six hours of monitoring a crowded field. Monitoring the orbits of several stars will provide the ability to distinguish between multiple post-Newtonian orbital effects, including those due to an extended mass distribution around Sgr A* and to low-order General Relativistic effects. Early characterizations of the field by ASTRA including the possibility of a precise source detection, could provide valuable information for future GRAVITY implementation and observation.

Algol (Beta Per) is an extensively studied hierarchical triple system whose inner pair is a prototype semi-detached binary with mass transfer occurring from the sub-giant secondary to the main-sequence primary. We present here the results of our Algol observations made between 2006 and 2010 at the CHARA interferometer with the Michigan Infrared Combiner in the H band. The use of four telescopes with long baselines allows us to achieve better than 0.5 mas resolution and to unambiguously resolve the three stars. The inner and outer orbital elements, as well as the angular sizes and mass ratios for the three components are determined independently from previous studies. We report a significantly improved orbit for the inner stellar pair with the consequence of a 15% change in the primary mass compared to previous studies. We also determine the mutual inclination of the orbits to be much closer to perpendicularity than previously established. State-of-the-art image reconstruction algorithms are used to image the full triple system. In particular an image sequence of 55 distinct phases of the inner pair orbit is reconstructed, clearly showing the Roche-lobe-filling secondary revolving around the primary, with several epochs corresponding to the primary and secondary eclipses.

V921 Scorpii is a close binary system (separation 0.025") showing the B[e]-phenomenon. The system is surrounded by an enigmatic bipolar nebula, which might have been shaped by episodic mass-loss events, possibly triggered by dynamical interactions between the companion and the circumprimary disk (Kraus et al. 2012a). In this paper, we investigate the spatial structure and kinematics of the circumprimary disk, with the aim to obtain new insights into the still strongly debated evolutionary stage. For this purpose, we combine, for the first time, infrared spectro-interferometry (VLTI/AMBER, R=12,000) and spectro-astrometry (VLT/CRIRES, R=100,000), which allows us to study the AU-scale distribution of circumstellar gas and dust with an unprecedented velocity resolution of 3 km*s^-1. Using a model-independent photocenter analysis technique, we find that the Br-gamma-line emission rotates in the same plane as the dust disk. We can reproduce the wavelength-differential visibilities and phases and the double-peaked line profile using a Keplerian-rotating disk model. The derived mass of the central star is 5.4+/-0.4 M_sun*(d/1150 pc), which is considerably lower than expected from the spectral classification, suggesting that V921 Sco might be more distant (d approx 2kpc) than commonly assumed. Using the geometric information provided by our Br-gamma spectro-interferometric data and Paschen, Brackett, and Pfund line decrement measurements in 61 hydrogen recombination line transitions, we derive the density of the line-emitting gas (N_e=2...6*10^19 m^-3). Given that our measurements can be reproduced with a Keplerian velocity field without outflowing velocity component and the non-detection of age-indicating spectroscopic diagnostics, our study provides new evidence for the pre-main-sequence nature of V921 Sco.

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: 1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; 2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image reconstruction process was "blind" in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with 6 to 10 telescopes, should be well suited to provide images of complex sceneries.

We present spatially resolved observations of the canonical transition disk object TW Hya at 8.74 {\mu}m, 11.66 {\mu}m, and 18.30 {\mu}m, obtained with the T-ReCS instrument on the Gemini telescope. These observations are a result of a novel observing mode at Gemini that enables speckle imaging. Using this technique, we image our target with short enough exposure times to achieve diffraction limited images. We use Fourier techniques to reduce our data, which allows high-precision calibration of the instrumental point spread function. Our observations span two epochs and we present evidence for temporal variability at 11.66 {\mu}m in the disk of TW Hya. We show that previous models of TW Hya's disk from the literature are incompatible with our observations, and construct a model to explain the discrepancies. We detect marginal asymmetry in our data, most significantly at the shortest wavelengths. To explain our data, we require a model that includes an optically thin inner disk extending from 0.02 to 3.9 AU, an optically thick ring representing the outer disk wall at 3.9 AU and extending to 4.6 AU, and a hotter-than-disk-equilibrium source of emission located at ~3.5 AU.

We report the detection of thermal emission from the hot Jupiter WASP-3b in the KS band, using a newly developed guiding scheme for the WIRC instrument at the Palomar Hale 200in telescope. Our new guiding scheme has improved the telescope guiding precision by a factor of ~5-7, significantly reducing the correlated systematics in the measured light curves. This results in the detection of a secondary eclipse with depth of 0.181%\pm0.020% (9-{\sigma}) - a significant improvement in WIRC's photometric precision and a demonstration of the capability of Palomar/WIRC to produce high quality measurements of exoplanetary atmospheres. Our measured eclipse depth cannot be explained by model atmospheres with heat redistribution but favor a pure radiative equilibrium case with no redistribution across the surface of the planet. Our measurement also gives an eclipse phase center of 0.5045\pm0.0020, corresponding to an ecos{\omega} of 0.0070\pm0.0032. This result is consistent with a circular orbit, although it also suggests the planet's orbit might be slightly eccentric. The possible non-zero eccentricity provides insight into the tidal circularization process of the star-planet system, but also might have been caused by a second low-mass planet in the system, as suggested by a previous transit timing variation study. More secondary eclipse observations, especially at multiple wavelengths, are necessary to determine the temperature-pressure profile of the planetary atmosphere and shed light on its orbital eccentricity.

\gamma Cas is the prototypical classical Be star and is best known for its variable hard X-ray emission. To elucidate the reasons for this emission, we mounted a multiwavelength campaign in 2010 centered around 4 XMM observations. The observational techniques included long baseline optical interferometry (LBOI), monitoring by an Automated Photometric Telescope and Halpha observations. Because gamma Cas is also known to be in a binary, we measured Halpha radial velocities and redetermined its period as 203.55+/-0.2 days and an eccentricity near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric in 2010, has an inclination angle near 45^o, and a larger radius than previously reported. The Be star began an "outburst" at the beginning of our campaign, made visible by a disk brightening and reddening during our campaign. Our analyses of the new high resolution spectra disclosed many attributes found from spectra obtained in 2001 (Chandra) and 2004 (XMM). As well as a dominant hot 14 keV thermal component, these familiar ones included: (i) a fluorescent feature of Fe K stronger than observed at previous times, (ii) strong lines of N VII and Ne XI lines indicative of overabundances, and (iii) a subsolar Fe abundance from K-shell lines but a solar abundance from L-shell ions. We also found that 2 absorption columns are required to fit the continuum. While the first one maintained its historical average of 1X10^21 cm^-2, the second was very large and doubled to 7.4X10^23 cm^-2 during our X-ray observations. Although we found no clear relation between this column density and orbital phase, it correlates well with the disk brightening and reddening both in the 2010 and earlier observations. Thus, the inference from this study is that much (perhaps all?) of the X-ray emission from this source originates behind matter ejected by gamma Cas into our line of sight.

Astronomers usually need the highest angular resolution possible, but the blurring effect of diffraction imposes a fundamental limit on the image quality from any single telescope. Interferometry allows light collected at widely-separated telescopes to be combined in order to synthesize an aperture much larger than an individual telescope thereby improving angular resolution by orders of magnitude. Radio and millimeter wave astronomers depend on interferometry to achieve image quality on par with conventional visible and infrared telescopes. Interferometers at visible and infrared wavelengths extend angular resolution below the milli-arcsecond level to open up unique research areas in imaging stellar surfaces and circumstellar environments. In this chapter the basic principles of interferometry are reviewed with an emphasis on the common features for radio and optical observing. While many techniques are common to interferometers of all wavelengths, crucial differences are identified that will help new practitioners avoid unnecessary confusion and common pitfalls. Concepts essential for writing observing proposals and for planning observations are described, depending on the science wavelength, angular resolution, and field of view required. Atmospheric and ionospheric turbulence degrades the longest-baseline observations by significantly reducing the stability of interference fringes. Such instabilities represent a persistent challenge, and the basic techniques of phase-referencing and phase closure have been developed to deal with them. Synthesis imaging with large observing datasets has become a routine and straightforward process at radio observatories, but remains challenging for optical facilities. In this context the commonly-used image reconstruction algorithms CLEAN and MEM are presented. Lastly, a concise overview of current facilities is included as an appendix.

Belonging to the group of B[e] stars, V921 Scorpii is associated with a strong infrared excess and permitted and forbidden line emission, indicating the presence of low- and high-density circumstellar gas and dust. Many aspects of V921 Sco and other B[e] stars still remain mysterious, including their evolutionary state and the physical conditions resulting in the class-defining characteristics. In this paper, we employ VLTI/AMBER spectro-interferometry in order to reconstruct high-resolution (lambda/2B=0.0013") model-independent interferometric images for three wavelength bands around 1.65, 2.0, and 2.3 micrometer. In our images, we discover a close (25.0+/-0.8 milliarcsecond, corresponding to 29+/-0.9 AU at 1.15 kpc) companion around V921 Sco. Between two epochs in 2008 and 2009, we measure orbital motion of 7 degrees, implying an orbital period of about 35 years (for a circular orbit). Around the primary star, we detect a disk-like structure with indications for a radial temperature gradient. The polar axis of this AU-scale disk is aligned with the arcminute-scale bipolar nebula in which V921 Sco is embedded. Using Magellan/IMACS imaging, we detect multi-layered arc-shaped sub-structure in the nebula, suggesting episodic outflow activity from the system with a period of about 25 years, roughly matching the estimated orbital period of the companion. Our study supports the hypothesis that the B[e] phenomenon is related to dynamical interaction in a close binary system.

Wolf-Rayet stars represent one of the final stages of massive stellar evolution. Relatively little is known about this short-lived phase and we currently lack reliable mass, distance, and binarity determinations for a representative sample. Here we report the first visual orbit for WR 140(=HD193793), a WC7+O5 binary system known for its periodic dust production episodes triggered by intense colliding winds near periastron passage. The IOTA and CHARA interferometers resolved the pair of stars in each year from 2003--2009, covering most of the highly-eccentric, 7.9 year orbit. Combining our results with the recent improved double-line spectroscopic orbit of Fahed et al. (2011), we find the WR 140 system is located at a distance of 1.67 +/- 0.03 kpc, composed of a WR star with M_WR = 14.9 +/- 0.5 Msun and an O star with M_O = 35.9 +/- 1.3 Msun. Our precision orbit yields key parameters with uncertainties times 6 smaller than previous work and paves the way for detailed modeling of the system. Our newly measured flux ratios at the near-infrared H and Ks bands allow an SED decomposition and analysis of the component evolutionary states.

We report on the preliminary results of a three-year monitoring campaign of Galactic and Magellanic LBVs with both moderate and high resolution spectroscopy. We have collected more than 500 moderate-resolution spectra of 20 stars in the red portion of the optical spectrum, which includes the wind-sensitive transitions of H$\alpha$ and He I 5876 and 6678 Angstroms. We summarize our detailed study of 15 years of high resolution H$\alpha$ spectroscopy of the prototypical luminous blue variable, P Cygni. We report on the discovery of discrete absorption components in P Cygni's H$\alpha$ profile, and we discuss their relationship to structure in the wind. These results are compared to our recent high resolution interferometric imaging of the $H$-band emitting region surrounding the star. We discuss recent observations of $\eta$ Car, highlighting its unusual recovery from the previous periastron passage. Our results on HDE 326823 indicate that the star has a short period binary orbit, with Roche Lobe overflow onto an unseen massive companion.

We report a new detection of the H-band thermal emission of CoRoT-1b and two confirmation detections of the Ks-band thermal emission of WASP-12b at secondary eclipses. The H-band measurement of CoRoT-1b shows an eclipse depth of 0.145%\pm0.049% with a 3-{\sigma} percentile between 0.033% - 0.235%. This depth is consistent with the previous conclusions that the planet has an isother- mal region with inefficient heat transport from dayside to nightside, and has a dayside thermal inversion layer at high altitude. The two Ks band detections of WASP-12b show a joint eclipse depth of 0.299%\pm0.065%. This result agrees with the measurement of Croll & collaborators, providing independent confirmation of their measurement. The repeatability of the WASP-12b measurements also validates our data analysis method. Our measurements, in addition to a number of previous results made with other telescopes, demonstrate that ground-based observations are becoming widely available for characterization of atmospheres of hot Jupiters.

Direct detection of thermal emission from nearby hot Jupiters has greatly advanced our knowledge of extrasolar planets in recent years. Since hot Jupiter systems can be regarded as analogs of high contrast binaries, ground-based infrared long baseline interferometers have the potential to resolve them and detect their thermal emission with precision closure phase - a method that is immune to the systematic errors induced by the Earth's atmosphere. In this work, we present closure phase studies toward direct detection of nearby hot Jupiters using the CHARA interferometer array outfitted with the MIRC instrument. We carry out closure phase simulations and conduct a large number of observations for the best candidate {\upsion} And. Our experiments suggest the method is feasible with highly stable and precise closure phases. However, we also find much larger systematic errors than expected in the observations, most likely caused by dispersion across different wavelengths. We find that using higher spectral resolution modes (e.g., R=150) can significantly reduce the systematics. By combining all calibrators in an observing run together, we are able to roughly recalibrate the lower spectral resolution data, allowing us to obtain upper limits of the star-planet contrast ratios of {\upsion} And b across the H band. The data also allow us to get a refined stellar radius of 1.625\pm0.011 R\odot. Our best upper limit corresponds to a contrast ratio of 2.1\times10^3:1 with 90% confidence level at 1.52{\mu}m, suggesting that we are starting to have the capability of constraining atmospheric models of hot Jupiters with interferometry. With recent and upcoming improvements of CHARA/MIRC, the prospect of detecting emission from hot Jupiters with closure phases is promising.

Using CHARA and VLTI near-infrared spectro-interferometry with hectometric baseline lengths (up to 330m) and with high spectral resolution (up to 12000), we studied the gas distribution and kinematics around two classical Be stars. The combination of high spatial and spectral resolution achieved allows us to constrain the gas velocity field on scales of a few stellar radii and to obtain, for the first time in optical interferometry, a dynamical mass estimate using the position-velocity analysis technique known from radio astronomy. For our first target star, Beta Canis Minoris, we model the H+K-band continuum and Br-gamma line geometry with a near-critical rotating stellar photosphere and a geometrically thin equatorial disk. Testing different disk rotation laws, we find that the disk is in Keplerian rotation (v(r)=r^(-0.5+/-0.1)) and derive the disk position angle (140+/-1.7 deg), inclination (38.5+/-1 deg), and the mass of the central star (3.5+/-0.2 M_sun). As a second target star, we observed the prototypical Be star Zeta Tauri and spatially resolved the Br{\gamma} emission as well as nine transitions from the hydrogen Pfund series (Pf14-Pf22). Comparing the spatial origin of the different line transitions, we find that the Brackett (Br-gamma), Pfund (Pf14-17), and Balmer (H-alpha) lines originate from different stellocentric radii (R_cont < R_Pf < R_Br-gamma = R_H-alpha), which we can reproduce with an LTE line radiative transfer computation. Discussing different disk-formation scenarios, we conclude that our constraints are inconsistent with wind compression models predicting a strong outflowing velocity component, but support viscous decretion disk models, where the Keplerian-rotating disk is replenished with material from the near-critical rotating star.

Although debris discs have been detected around a significant number of main-sequence stars, only a few of them are known to harbour hot dust in their inner part where terrestrial planets may have formed. Thanks to infrared interferometric observations, it is possible to obtain a direct measurement of these regions, which are of prime importance for preparing future exo-Earth characterisation missions. In this context, we have resolved the exozodiacal dust disc around Vega with the help of infrared stellar interferometry and estimated the integrated H-band flux originating from the first few AUs of the debris disc. Using precise H-band interferometric measurements obtained with the 3-telescope IOTA/IONIC interferometer (Mount Hopkins, Arizona), thorough modelling of both interferometric data (squared visibility and closure phase) and spectral energy distribution was performed to constrain the nature of the near-infrared excess emission. The most straightforward scenario consists in a compact dust disc producing a thermal emission that is largely dominated by small grains located between 0.1 and 0.3 AU from Vega and accounting for 1.23 +/- 0.45% of the near-infrared stellar flux for our best-fit model. This flux ratio is shown to vary slightly with the geometry of the model used to fit our interferometric data (variations within +/-0.19%). Initially revealed by K-band CHARA/FLUOR observations, the presence of hot exozodiacal dust in the vicinity of Vega is confirmed by our H-band IOTA/IONIC measurements at the 3-sigma level. Whereas the origin of the dust is still uncertain, its presence and the possible connection with the outer disc suggest that the Vega system is currently undergoing major dynamical perturbations.

The CHARA Array is a six 1-m telescope optical and near infrared interferometer located at the Mount Wilson Observatory in southern California and operated by the Center for High Angular Resolution Astronomy of Georgia State University. The CHARA Array has been in regular scientific operation since 2005 and now has over 55 publications in the refereed literature, including two in Science and one in Nature. The Array now supports seven beam combiners ranging from 0.5 microns up to 2.3 microns and combing from 2 to 4 beams at a time. An upgrade to a full 6 beam combiner is now underway and fringes with all six telescopes were achieved soon after the meeting.We present some of the more recent results from the CHARA-Array.

Near-infrared interferometers have recently imaged a number of rapidly rotating A-type stars, finding levels of gravity darkening inconsistent with theoretical expectations. Here, we present new imaging of both a cooler star {\beta} Cas (F2IV) and a hotter one {\alpha} Leo (B7V) using the CHARA array and the MIRC instrument at the H band. Adopting a solid-body rotation model with a simple gravity darkening prescription, we modeled the stellar geometric properties and surface temperature distributions, confirming both stars are rapidly rotating and show gravity darkening anomalies. We estimate the masses and ages of these rapid rotators on L-Rpol and HR diagrams constructed for non-rotating stars by tracking their non-rotating equivalents. The unexpected fast rotation of the evolved subgiant {\beta} Cas offers a unique test of the stellar core-envelope coupling, revealing quite efficient coupling over the past ~ 0.5 Gyr. Lastly we summarize all our interferometric determinations of the gravity darkening coefficient for rapid rotators, finding none match the expectations from the widely used von Zeipel gravity darkening laws. Since the conditions of the von Zeipel law are known to be violated for rapidly rotating stars, we recommend using the empirically-derived {\beta} = 0.19 for such stars with radiation-dominated envelopes. Furthermore, we note that no paradigm exists for self-consistently modeling heavily gravity-darkened stars that show hot radiative poles with cool convective equators.

Young and close multiple systems are unique laboratories to probe the initial dynamical interactions between forming stellar systems and their dust and gas environment. Their study is a key building block to understanding the high frequency of main-sequence multiple systems. However, the number of detected spectroscopic young multiple systems that allow dynamical studies is limited. GW Orionis is one such system. It is one of the brightest young T Tauri stars and is surrounded by a massive disk. Our goal is to probe the GW Orionis multiplicity at angular scales at which we can spatially resolve the orbit. We used the IOTA/IONIC3 interferometer to probe the environment of GW Orionis with an astronomical unit resolution in 2003, 2004, and 2005. By measuring squared visibilities and closure phases with a good UV coverage we carry out the first image reconstruction of GW Ori from infrared long-baseline interferometry. We obtain the first infrared image of a T Tauri multiple system with astronomical unit resolution. We show that GW Orionis is a triple system, resolve for the first time the previously known inner pair (separation $\rho\sim$1.4 AU) and reveal a new more distant component (GW Ori C) with a projected separation of $\sim$8 AU with direct evidence of motion. Furthermore, the nearly equal (2:1) H-band flux ratio of the inner components suggests that either GW Ori B is undergoing a preferential accretion event that increases its disk luminosity or that the estimate of the masses has to be revisited in favour of a more equal mass-ratio system that is seen at lower inclination. Accretion disk models of GW Ori will need to be completely reconsidered because of this outer companion C and the unexpected brightness of companion B.

Integrated optic beam combiners offer many advantages over conventional bulk optic implementations for astronomical imaging. To date, integrated optic beam combiners have only been demonstrated at operating wavelengths below 4 microns. Operation in mid-infrared wavelength region, however, is highly desirable. In this paper, a theoretical design technique based on three coupled waveguides is developed to achieve fully achromatic, broadband, polarization-insensitive, lossless beam combining. This design may make it possible to achieve the very deep broadband nulls needed for exoplanet searching.

Despite a century of remarkable progress in understanding stellar interiors, we know surprisingly little about the inner workings of stars spinning near their critical limit. New interferometric imaging of these so-called ``rapid rotators'' combined with breakthroughs in asteroseismology promise to lift this veil and probe the strongly latitude-dependent photospheric characteristics and even reveal the internal angular momentum distribution of these luminous objects. Here, we report the first high precision photometry on the low-amplitude delta cuti variable star Rasalhague (alpha Oph, A5IV, 2.18 Msun, omega/omega_c~0.88) based on 30 continuous days of monitoring using the MOST satellite. We have identified 57+/-1 distinct pulsation modes above a stochastic granulation spectrum with a cutoff of ~26 cycles per day. Remarkably, we have also discovered that the fast rotation period of 14.5~hours modulates low-frequency modes (1-10 day periods) that we identify as a rich family of g-modes (|m| up to 7). The spacing of the g-modes is surprisingly linear considering Coriolis forces are expected to strongly distort the mode spectrum, suggesting we are seeing prograde ``equatorial Kelvin'' waves (modes l=m). We emphasize the unique aspects of Rasalhague motivating future detailed asteroseismic modeling -- a source with a precisely measured parallax distance, photospheric oblateness, latitude temperature structure, and whose low-mass companion provides an astrometric orbit for precise mass determinations.

The nearby star Alpha Oph (Ras Alhague) is a rapidly rotating A5IV star spinning at ~89% of its breakup velocity. This system has been imaged extensively by interferometric techniques, giving a precise geometric model of the star's oblateness and the resulting temperature variation on the stellar surface. Fortuitously, Alpha Oph has a previously known stellar companion, and characterization of the orbit provides an independent, dynamically-based check of both the host star and the companion mass. Such measurements are crucial to constrain models of such rapidly rotating stars. In this study, we combine eight years of Adaptive Optics imaging data from the Palomar, AEOS, and CFHT telescopes to derive an improved, astrometric characterization of the companion orbit. We also use photometry from these observations to derive a model-based estimate of the companion mass. A fit was performed on the photocenter motion of this system to extract a component mass ratio. We find masses of 2.40^{0.23}_{0.37} solar masses and 0.85^{0.06}_{0.04} solar masses for Alpha Oph A and Alpha Oph B, respectively. Previous orbital studies of this system found a mass too high for this system, inconsistent with stellar evolutionary calculations. Our measurements of the host star mass are more consistent with these evolutionary calculations, but with slightly higher uncertainties. In addition to the dynamically-derived masses, we use IJHK photometry to derive a model-based mass for Alpha Oph B, of 0.77 +/- 0.05 solar masses marginally consistent with the dynamical masses derived from our orbit. Our model fits predict a periastron passage on 2012 April 19, with the two components having a ~50 milliarcsec separation from March to May 2012. A modest amount of interferometric and radial velocity data during this period could provide a mass determination of this star at the few percent level.

We present interferometric observations of the Be star Zeta Tau obtained using the MIRC beam combiner at the CHARA Array. We resolved the disk during four epochs in 2007-2009. We fit the data with a geometric model to characterize the circumstellar disk as a skewed elliptical Gaussian and the central Be star as a uniform disk. The visibilities reveal a nearly edge-on disk with a FWHM major axis of ~ 1.8 mas in the H-band. The non-zero closure phases indicate an asymmetry within the disk. Interestingly, when combining our results with previously published interferometric observations of Zeta Tau, we find a correlation between the position angle of the disk and the spectroscopic V/R ratio, suggesting that the tilt of the disk is precessing. This work is part of a multi-year monitoring campaign to investigate the development and outward motion of asymmetric structures in the disks of Be stars.

We present observations and models for one of these MYSO candidates, NGC3603 IRS 9A. Our goal is to investigate with infrared interferometry the structure of IRS 9A on scales as small as 200AU, exploiting the fact that a cluster of O and B stars has blown away much of the obscuring foreground dust and gas. Observations in the N-band were carried out with the MIDI beam combiner attached to the VLTI. Additional interferometric observations which probe the structure of IRS 9A on larger scales were performed with an aperture mask installed in the T-ReCS instrument of Gemini South. The spectral energy distribution (SED) is constrained by the MIDI N-band spectrum and by data from the Spitzer Space Telescope. Our efforts to model the structure and SED of IRS 9A range from simple geometrical models of the brightness distribution to one- and two-dimensional radiative transfer computations. The target is resolved by T-ReCS, with an equivalent (elliptical) Gaussian width of 330mas by 280mas (2300 AU by 2000 AU). Despite this fact, a warm compact unresolved component was detected by MIDI which is possibly associated with the inner regions of a flattened dust distribution. Based on our interferometric data, no sign of multiplicity was found on scales between about 200AU and 700AU projected separation. A geometric model consisting of a warm (1000 K) ring (400 AU diameter) and a cool (140 K) large envelope provides a good fit to the data. No single model fitting all visibility and photometric data could be found, with disk models performing better than spherical models. While the data are clearly inconsistent with a spherical dust distribution they are insufficient to prove the existence of a disk but rather hint at a more complex dust distribution.

We present near-infrared H and K-band spectro-interferometric observations of the gaseous disk around the primary Be star in the delta Sco binary system, obtained in 2007 (between periastron passages in 2000 and 2011). Observations using the CHARA/MIRC instrument at H-band resolve an elongated disk with a Gaussian FWHM 1.18 x 0.91 mas. Using the Keck Interferometer, the source of the K-band continuum emission is only marginally spatially resolved, and consequently we estimate a relatively uncertain K-band continuum disk FWHM of 0.7 +/- 0.3 mas. Line emission on the other hand, He1 (2.0583 micron) and Br gamma (2.1657 micron), is clearly detected, with about 10% lower visibilities than those of the continuum. When taking into account the continuum/line flux ratio this translates into much larger sizes for the line emission regions: 2.2 +/- 0.4 mas and 1.9 +/- 0.3 mas for He1 and Br gamma respectively. Our KI data also reveal a relatively flat spectral differential phase response, ruling out significant off-center emission. We expect these new measurements will help constrain dynamical models being actively developed in order to explain the disk formation process in the delta Sco system and Be stars in general.

The ASTrometric and phase-Referenced Astronomy (ASTRA) project will provide phase referencing and astrometric observations at the Keck Interferometer, leading to enhanced sensitivity and the ability to monitor orbits at an accuracy level of 30-100 microarcseconds. Here we discuss recent scientific results from ASTRA, and describe new scientific programs that will begin in 2010-2011. We begin with results from the "self phase referencing" (SPR) mode of ASTRA, which uses continuum light to correct atmospheric phase variations and produce a phase-stabilized channel for spectroscopy. We have observed a number of protoplanetary disks using SPR and a grism providing a spectral dispersion of ~2000. In our data we spatially resolve emission from dust as well as gas. Hydrogen line emission is spectrally resolved, allowing differential phase measurements across the emission line that constrain the relative centroids of different velocity components at the 10 microarcsecond level. In the upcoming year, we will begin dual-field phase referencing (DFPR) measurements of the Galactic Center and a number of exoplanet systems. These observations will, in part, serve as precursors to astrometric monitoring of stellar orbits in the Galactic Center and stellar wobbles of exoplanet host stars. We describe the design of several scientific investigations capitalizing on the upcoming phase-referencing and astrometric capabilities of ASTRA.

We report on the successful science verification phase of a new observing mode at the Keck interferometer, which provides a line-spread function width and sampling of 150km/s at K'-band, at a current limiting magnitude of K'~7mag with spatial resolution of lam/2B ~2.7mas and a measured differential phase stability of unprecedented precision (3mrad at K=5mag, which represents 3uas on sky or a centroiding precision of 10^-3). The scientific potential of this mode is demonstrated by the presented observations of the circumstellar disk of the evolved Be-star 48Lib. In addition to indirect methods such as multi-wavelength spectroscopy and polaritmetry, the here described spectro-interferometric astrometry provides a new tool to directly constrain the radial density structure in the disk. We resolve for the first time several Pfund emission lines, in addition to BrGam, in a single interferometric spectrum, and with adequate spatial and spectral resolution and precision to analyze the radial disk structure in 48Lib. The data suggest that the continuum and Pf-emission originates in significantly more compact regions, inside of the BrGam emission zone. Thus, spectro-interferometric astrometry opens the opportunity to directly connect the different observed line profiles of BrGam and Pfund in the total and correlated flux to different disk radii. The gravitational potential of a rotationally flattened Be star is expected to induce a one-armed density perturbation in the circumstellar disk. Such a slowly rotating disk oscillation has been used to explain the well known periodic V/R spectral profile variability in these stars, as well as the observed V/R cycle phase shifts between different disk emission lines. The differential line properties and linear constraints set by our data lend support to the existence of a radius-dependent disk density perturbation.

We present here a new observational technique, Phase Closure Nulling (PCN), which has the potential to obtain very high contrast detection and spectroscopy of faint companions to bright stars. PCN consists in measuring closure phases of fully resolved objects with a baseline triplet where one of the baselines crosses a null of the object visibility function. For scenes dominated by the presence of a stellar disk, the correlated flux of the star around nulls is essentially canceled out, and in these regions the signature of fainter, unresolved, scene object(s) dominates the imaginary part of the visibility in particular the closure phase. We present here the basics of the PCN method, the initial proof-of-concept observation, the envisioned science cases and report about the first observing campaign made on VLTI/AMBER and CHARA/MIRC using this technique.

To understand how planetary systems form in the dusty disks around pre-main-sequence stars a detailed knowledge of the structure and evolution of these disks is required. While this is reasonably well understood for the regions of the disk beyond about 1 AU, the structure of these disks inward of 1 AU remains a puzzle. This is partly because it is very difficult to spatially resolve these regions with current telescopes. But it is also because the physics of this region, where the disk becomes so hot that the dust starts to evaporate, is poorly understood. With infrared interferometry it has become possible in recent years to directly spatially resolve the inner AU of protoplanetary disks, albeit in a somewhat limited way. These observations have partly confirmed current models of these regions, but also posed new questions and puzzles. Moreover, it has turned out that the numerical modeling of these regions is extremely challenging. In this review we give a rough overview of the history and recent developments in this exciting field of astrophysics.

We present near-infrared observations of T Tauri and Herbig Ae/Be stars with a spatial resolution of a few milli-arcseconds and a spectral resolution of ~2000. Our observations spatially resolve gas and dust in the inner regions of protoplanetary disks, and spectrally resolve broad-linewidth emission from the Brackett gamma transition of hydrogen gas. We use the technique of spectro-astrometry to determine centroids of different velocity components of this gaseous emission at a precision orders of magnitude better than the angular resolution. In all sources, we find the gaseous emission to be more compact than or distributed on similar spatial scales to the dust emission. We attempt to fit the data with models including both dust and Brackett gamma-emitting gas, and we consider both disk and infall/outflow morphologies for the gaseous matter. In most cases where we can distinguish between these two models, the data show a preference for infall/outflow models. In all cases, our data appear consistent with the presence of some gas at stellocentric radii of ~0.01 AU. Our findings support the hypothesis that Brackett gamma emission generally traces magnetospherically driven accretion and/or outflows in young star/disk systems.

Eclipses of the single-line spectroscopic binary star, epsilon Aurigae, provide an opportunity to study the poorly-defined companion. We used the MIRC beam combiner on the CHARA array to create interferometric images during eclipse ingress. Our results demonstrate that the eclipsing body is a dark disk that is opaque and tilted, and therefore exclude alternative models for the system. These data constrain the geometry and masses of the components, providing evidence that the F-star is not a massive supergiant star.

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.

This paper reports on H-band interferometric observations of Betelgeuse made at the three-telescope interferometer IOTA. We image Betelgeuse and its asymmetries to understand the spatial variation of the photosphere, including its diameter, limb darkening, effective temperature, surrounding brightness, and bright (or dark) star spots. We used different theoretical simulations of the photosphere and dusty environment to model the visibility data. We made images with parametric modeling and two image reconstruction algorithms: MIRA and WISARD. We measure an average limb-darkened diameter of 44.28 +/- 0.15 mas with linear and quadratic models and a Rosseland diameter of 45.03 +/- 0.12 mas with a MARCS model. These measurements lead us to derive an updated effective temperature of 3600 +/- 66 K. We detect a fully-resolved environment to which the silicate dust shell is likely to contribute. By using two imaging reconstruction algorithms, we unveiled two bright spots on the surface of Betelgeuse. One spot has a diameter of about 11 mas and accounts for about 8.5% of the total flux. The second one is unresolved (diameter < 9 mas) with 4.5% of the total flux. Resolved images of Betelgeuse in the H band are asymmetric at the level of a few percent. The MOLsphere is not detected in this wavelength range. The amount of measured limb-darkening is in good agreement with model predictions. The two spots imaged at the surface of the star are potential signatures of convective cells.

Integrated-optic, astronomical, two-beam and three-beam, interferometric combiners have been designed and fabricated for operation in the L band (3 - 4 microns) for the first time. The devices have been realized in titanium-indiffused, x-cut lithium niobate substrates, and on-chip electro-optic fringe scanning has been demonstrated. White light fringes were produced in the laboratory using the two-beam combiner integrated with an on-chip Y-splitter.

We present infrared interferometric imaging of the S-type Mira star Chi Cygni. The object was observed at four different epochs in 2005-2006 with the IOTA optical interferometer (H band). Images show up to 40% variation in the stellar diameter, as well as significant changes in the limb darkening and stellar inhomogeneities. Model fitting gave precise time-dependent values of the stellar diameter, and reveals presence and displacement of a warm molecular layer. The star radius, corrected for limb darkening, has a mean value of 12.1 mas and shows a 5.1mas amplitude pulsation. Minimum diameter was observed at phase 0.94+/-0.01. Maximum temperature was observed several days later at phase 1.02+/-0.02. We also show that combining the angular acceleration of the molecular layer with CO (Delta v = 3) radial velocity measurements yields a 5.9+/-1.5 mas parallax. The constant acceleration of the CO molecules -- during 80% of the pulsation cycle -- lead us to argument for a free-falling layer. The acceleration is compatible with a gravitational field produced by a 2.1(+1.5/-0.7) solar mass star. This last value is in agreement with fundamental mode pulsator models. We foresee increased development of techniques consisting in combining radial velocity with interferometric angular measurements, ultimately allowing total mapping of the speed, density, and position of the diverse species in pulsation driven atmospheres.

We present observations of the 51 Ophiuchi circumstellar disk made with the Keck interferometer operating in nulling mode at N-band. We model these data simultaneously with VLTI-MIDI visibility data and a Spitzer IRS spectrum using a variety of optically-thin dust cloud models and an edge-on optically-thick disk model. We find that single-component optically-thin disk models and optically-thick disk models are inadequate to reproduce the observations, but an optically-thin two-component disk model can reproduce all of the major spectral and interferometric features. Our preferred disk model consists of an inner disk of blackbody grains extending to ~4 AU and an outer disk of small silicate grains extending out to ~1200 AU. Our model is consistent with an inner "birth" disk of continually colliding parent bodies producing an extended envelope of ejected small grains. This picture resembles the disks around Vega, AU Microscopii, and Beta Pictoris, supporting the idea that 51 Ophiuchius may be a Beta Pictoris analog.

We present sub-milliarcseond resolution imaging and modeling of two nearby rapid rotators Alpha Cephei and Alpha Ophiuchi, obtained with the CHARA array - the largest optical/IR interferometer in the world. Incorporating a gravity darkening model, we are able to determine the inclination, the polar and equatorial radius and temperature, as well as the fractional rotation speed of the two stars with unprecedented precision. The polar and equatorial regions of the two stars have ~2000K temperature gradient, causing their apparent temperatures and luminosities to be dependent on their viewing angles. Our modeling allow us to determine the true effective temperatures and luminosities of Alpha Cep and Alpha Oph, permitting us to investigate their true locations on the H-R diagram. These properties in turn give us estimates of the masses and ages of the two stars within a few percent of error using stellar evolution models. Also, based on our gravity darkening modeling, we propose a new method to estimate the masses of single stars in a more direct way through Vsin(i) measurements and precise geometrical constraint. Lastly, we investigate the degeneracy between the inclination and the gravity darkening coefficient, which especially affects the modeling of Alpha Oph. Although incorporating Vsin(i) has lifted the degeneracy to some extent, higher resolution observations are still needed to further constrain the parameters independently.

The mid-infrared properties of pre-planetary disks are sensitive to the temperature and flaring profiles of disks for the regions where planet formation is expected to occur. In order to constrain theories of planet formation, we have carried out a mid-infrared (wavelength 10.7 microns) size survey of young stellar objects using the segmented Keck telescope in a novel configuration. We introduced a customized pattern of tilts to individual mirror segments to allow efficient sparse-aperture interferometry, allowing full aperture synthesis imaging with higher calibration precision than traditional imaging. In contrast to previous surveys on smaller telescopes and with poorer calibration precision, we find most objects in our sample are partially resolved. Here we present the main observational results of our survey of 5 embedded massive protostars, 25 Herbig Ae/Be stars, 3 T Tauri stars, 1 FU Ori system, and 5 emission-line objects of uncertain classification. The observed mid-infrared sizes do not obey the size-luminosity relation found at near-infrared wavelengths and a companion paper will provide further modelling analysis of this sample. In addition, we report imaging results for a few of the most resolved objects, including complex emission around embedded massive protostars, the photoevaporating circumbinary disk around MWC 361A, and the subarcsecond binaries T Tau, FU Ori and MWC 1080.

We present mid-IR observations from Gemini/TReCS that spatially resolve the dust emission around SR 21. The protoplanetary disk around SR 21 is believed to have a cleared gap extending from stellocentric radii of ~0.5 AU to ~20 AU, based on modeling of the observed spectral energy distribution. Our new observations resolve the dust emission, and our data are inconsistent with the previous model. We require the disk to be completely cleared within ~10 AU, without the hot inner disk spanning ~0.25-0.5 AU posited previously. To fit the SED and mid-IR imaging data together, we propose a disk model with a large inner hole, but with a warm companion--possibly surrounded by circumstellar material of its own--residing near the outer edge of the cleared region. We also discuss a model with a narrow ring included in a large cleared inner disk region, and argue that it is difficult to reconcile with the data.

We present the first high angular resolution observation of the B[e] star/X-ray transient object CI Cam, performed with the two-telescope Infrared Optical Telescope Array (IOTA), its upgraded three-telescope version (IOTA3T) and the Palomar Testbed Interferometer (PTI). Visibilities and closure phases were obtained using the IONIC-3 integrated optics beam combiner. CI Cam was observed in the near-infrared H and K spectral bands, wavelengths well suited to measure the size and study the geometry of the hot dust surrounding CI Cam. The analysis of the visibility data over an 8 year period from soon after the 1998 outburst to 2006 shows that the dust visibility has not changed over the years. The visibility data shows that CI Cam is elongated which confirms the disc-shape of the circumstellar environment and totally rules out the hypothesis of a spherical dust shell. Closure phase measurements show direct evidence of asymmetries in the circumstellar environment of CI Cam and we conclude that the dust surrounding CI Cam lies in an inhomogeneous disc seen at an angle. The near-infrared dust emission appears as an elliptical skewed Gaussian ring with a major axis a = 7.58 +/- 0.24 mas, an axis ratio r = 0.39 +/- 0.03 and a position angle theta = 35 +/- 2 deg.

We have detected asymmetry in the symbiotic star CH Cyg through the measurement of precision closure-phase with the IONIC beam combiner, at the IOTA interferometer. The position of the asymmetry changes with time and is correlated with the phase of the 2.1-yr period found in the radial velocity measurements for this star. We can model the time-dependent asymmetry either as the orbit of a low-mass companion around the M giant or as an asymmetric, 20% change in brightness across the M giant. We do not detect a change in the size of the star during a 3 year monitoring period neither with respect to time nor with respect to wavelength. We find a spherical dust-shell with an emission size of 2.2+/-0.1 D* FWHM around the M giant star. The star to dust flux ratio is estimated to be 11.63+/-0.3. While the most likely explanation for the 20% change in brightness is non-radial pulsation we argue that a low-mass companion in close orbit could be the physical cause of the pulsation. The combined effect of pulsation and low-mass companion could explain the behaviour revealed by the radial-velocity curves and the time-dependent asymmetry detected in the closure-phase data. If CH Cyg is a typical long secondary period variable then these variations could be explained by the effect of an orbiting low-mass companion on the primary star.

We have obtained H-band interferometric observations of three galactic red supergiant stars using the MIRC instrument on the CHARA array. The targets include AZ Cyg, a field RSG and two members of the Per OB1 association, RS Per and T Per. We find evidence of departures from circular symmetry in all cases, which can be modelled with the presence of hotspots. This is the first detection of these features in the $H$-band. The measured mean diameters and the spectral energy distributions were combined to determine effective temperatures. The results give further support to the recently derived hotter temperature scale of red supergiant stars by Levesque et al. (2005), which has been evoked to reconcile the empirically determined physical parameters and stellar evolutionary theories. We see a possible correlation between spottedness and mid-IR emission of the circumstellar dust, suggesting a connection between mass-loss and the mechanism that generates the spots.

The formation of planets is one of the major unsolved problems in modern astrophysics. Planets are believed to form out of the material in circumstellar disks known to exist around young stars, and which are a by-product of the star formation process. Therefore, the physical conditions in these disks - structure and composition as a function of stellocentric radius and vertical height, density and temperature profiles of each component - represent the initial conditions under which planets form. Clearly, a good understanding of disk structure and its time evolution are crucial to understanding planet formation, the evolution of young planetary systems (e.g. migration), and the recently discovered, and unanticipated, diversity of planetary architectures. However, the inner disk regions (interior to ~10 AU) most relevant in the context of planet formation are very poorly known, primarily because of observational challenges in spatially resolving this region. In this contribution we discuss opportunities for the next decade from spectrally and spatially resolved observations, and from direct imaging, using infrared long baseline interferometry.

The sensitivity and astrometry upgrade ASTRA of the Keck Interferometer is introduced. After a brief overview of the underlying interferometric principles, the technology and concepts of the upgrade are presented. The interferometric dual-field technology of ASTRA will provide the KI with the means to observe two objects simultaneously, and measure the distance between them with a precision eventually better than 100 uas. This astrometric functionality of ASTRA will add a unique observing tool to fields of astrophysical research as diverse as exo-planetary kinematics, binary astrometry, and the investigation of stars accelerated by the massive black hole in the center of the Milky Way as discussed in this contribution.

We present results from a spectro-interferometric study of the Miras o Cet, R Leo and W Hya obtained with the Keck Aperture Masking Experiment from 1998 Sep to 2002 Jul. The spectrally dispersed visibility data permit fitting with circularly symmetric brightness profiles such as a simple uniform disk. The stellar angular diameter obtained over up to ~ 450 spectral channels spaning the region 1.1-3.8 microns is presented. Use of a simple uniform disk brightness model facilitates comparison between epochs and with existing data and theoretical models. Strong size variations with wavelength were recorded for all stars, probing zones of H2O, CO, OH, and dust formation. Comparison with contemporaneous spectra extracted from our data show a strong anti-correlation between the observed angular diameter and flux. These variations consolidate the notion of a complex stellar atmosphere consisting of molecular shells with time-dependent densities and temperatures. Our findings are compared with existing data and pulsation models. The models were found to reproduce the functional form of the wavelength vs. angular diameter curve well, although some departures are noted in the 2.8-3.5 micron range.

We present infrared interferometric observations of the inner regions of two A-star debris disks, beta Leo and zeta Lep, using the FLUOR instrument at the CHARA interferometer on both short (30 m) and long (>200 m) baselines. For the target stars, the short baseline visibilities are lower than expected for the stellar photosphere alone, while those of a check star, delta Leo, are not. We interpret this visibility offset of a few percent as a near-infrared excess arising from dust grains which, due to the instrumental field of view, must be located within several AU of the central star. For beta Leo, the near-infrared excess producing grains are spatially distinct from the dust which produces the previously known mid-infrared excess. For zeta Lep, the near-infrared excess may be spatially associated with the mid-infrared excess producing material. We present simple geometric models which are consistent with the near and mid-infrared excess and show that for both objects, the near-infrared producing material is most consistent with a thin ring of dust near the sublimation radius with typical grain sizes smaller than the nominal radiation pressure blowout radius. Finally, we discuss possible origins of the near-infrared emitting dust in the context of debris disk evolution models.

We present comprehensive models for the Herbig Ae stars MWC275 and AB Aur that aim to explain their spectral energy distribution (from UV to millimeter) and long baseline interferometry (from near-infrared to millimeter) simultaneously. Data from the literature, combined with new mid-infrared (MIR) interferometry from the Keck Segment Tilting Experiment, are modeled using an axisymmetric Monte Carlo radiative transfer code. Models in which most of the near-infrared (NIR) emission arises from a dust rim fail to fit the NIR spectral energy distribution (SED) and sub-milli-arcsecond NIR CHARA interferometry. Following recent work, we include an additional gas emission component with similar size scale to the dust rim, inside the sublimation radius, to fit the NIR SED and long-baseline NIR interferometry on MWC275 and AB Aur. In the absence of shielding of star light by gas, we show that the gas-dust transition region in these YSOs will have to contain highly refractory dust, sublimating at ~1850K. Despite having nearly identical structure in the thermal NIR, the outer disks of MWC275 and AB Aur differ substantially. In contrast to the AB Aur disk, MWC275 lacks small grains in the disk atmosphere capable of producing significant 10-20micron emission beyond ~7AU, forcing the outer regions into the "shadow" of the inner disk

We present the first resolved images of the eclipsing binary Beta Lyrae, obtained with the CHARA Array interferometer and the MIRC combiner in the H band. The images clearly show the mass donor and the thick disk surrounding the mass gainer at all six epochs of observation. The donor is brighter and generally appears elongated in the images, the first direct detection of photospheric tidal distortion due to Roche-lobe filling. We also confirm expectations that the disk component is more elongated than the donor and is relatively fainter at this wavelength. Image analysis and model fitting for each epoch were used for calculating the first astrometric orbital solution for Beta Lyrae, yielding precise values for the orbital inclination and position angle. The derived semi-ma jor axis also allows us to estimate the distance of Beta Lyrae; however, systematic differences between the models and the images limit the accuracy of our distance estimate to about 15%. To address these issues, we will need a more physical, self-consistent model to account for all epochs as well as the multi-wavelength information from the eclipsing light curves.

We report the discovery of a circumbinary disk around the Herbig Ae/Be system v892 Tau. Our detailed mid-infrared images were made using segment-tilting interferometry on the Keck-1 Telescope and reveal an asymmetric disk inclined at ~60 degs with an inner hole diameter of 250 mas (35 AU), approximately 5X larger than the apparent separation of the binary components. In addition, we report a new measurement along the binary orbit using near-infrared Keck aperture masking, allowing a crude estimate of orbital parameters and the system mass for the first time. The size of the inner hole appears to be consistent with the minimum size prediction from tidal truncation theory, bearing a resemblance to the recently unmasked binary CoKu Tau/4. Our results have motivated a re-analysis of the system spectral energy distribution, concluding the luminosity of this system has been severely underestimated. With further study and monitoring, v892 Tau should prove a powerful testing ground for both predictions of dynamical models for disk-star interactions in young systems with gas-rich disks and for calibrations of pre-main-sequence tracks for intermediate-mass stars.

This paper is an H band interferometric examination of Arcturus, a star frequently used as a spatial and spectral calibrator. Using the IOTA 3 telescope interferometer, we performed spectro-interferometric observation (R~35) of Arcturus. Atmospheric models and prescriptions were fitted to the data to derive the brightness distribution of the photosphere. Image reconstruction was also obtained using two software algorithms: Wisard and Mira. An achromatic power law proved to be a good model of the brightness distribution, with a limb darkening compatible with the one derived from atmospheric model simulations using our Marcs model. A Rosseland diameter of 21.05+/-0.21 was derived, corresponding to an effective temperature of T_eff=4295+/-26 K. No companion was detected from the closure phases, with an upper limit on the brightness ratio of 8e-4 at 1AU. Dynamic range at such distance from the photosphere was established at 1.5e-4 (1sigma rms). An upper limit of 1.7e-3 was also derived for the level of brightness asymmetries present on the photosphere.

Using the longest optical-interferometeric baselines currently available, we have detected strong near-infrared (NIR) emission from inside the dust-destruction radius of Herbig Ae stars MWC275 and AB Aur. Our sub-milli-arcsecond resolution observations unambiguously place the emission between the dust-destruction radius and the magnetospheric co-rotation radius. We argue that this new component corresponds to hot gas inside the dust-sublimation radius, confirming recent claims based on spectrally-resolved interferometry and dust evaporation front modeling.

During the last few years, considerable effort has been directed towards large-scale (>> $1 Billion US) missions to detect and characterize earth-like planets around nearby stars, such as the Terrestrial Planet Finder Interferometer (TPF-I) and Darwin missions. However, technological and budgetary issues as well as shifting science priorities will likely prevent these missions from entering Phase A until the next decade. The secondary eclipse technique using the Spitzer Space Telescope has been used to directly measure the temperature and emission spectrum of extrasolar planets. However, only a small fraction of known extrasolar planets are in transiting orbits. Thus, a simplified nulling interferometer, which produces an artificial eclipse or occultation, and operates in the near- to mid-infrared (e.g. ~ 3 to 8 or 10 microns), can characterize the atmospheres of this much larger sample of the known but non-transiting exoplanets. Many other scientific problems can be addressed with a system like this, including imaging debris disks, active galactic nuclei, and low mass companions around nearby stars. We discuss the rationale for a probe-scale mission in the $600-800 Million range, which we name here as the Small Prototype Planet Finding Interferometer (SPPFI).

We report imaging observations of the symbotic long-period Mira variable R Aquarii (R Aqr) at near-infrared and radio wavelengths. The near-infrared observations were made with the IOTA imaging interferometer in three narrow-band filters centered at 1.51, 1.64, and 1.78 $\mu$m, which sample mainly water, continuum, and water features, respectively. Our near-infrared fringe visibility and closure phase data are analyzed using three models. (a) A uniform disk model with wavelength-dependent sizes fails to fit the visibility data, and is inconsistent with the closure phase data. (b) A three- component model, comprising a Mira star, water shell, and an off-axis point source, provide a good fit to all data. (c) A model generated by a constrained image reconstruction analysis provides more insight, suggesting that the water shell is highly non-uniform, i.e., clumpy. The VLBA observations of SiO masers in the outer molecular envelope show evidence of turbulence, with jet-like features containing velocity gradients.

Results from the most extensive study of the time-evolving dust structure around the prototype "Pinwheel" nebula WR 104 are presented. Encompassing 11 epochs in three near-infrared filter bandpasses, a homogeneous imaging data set spanning more than 6 years (or 10 orbits) is presented. Data were obtained from the highly successful Keck Aperture Masking Experiment, which can recover high fidelity images at extremely high angular resolutions, revealing the geometry of the plume with unprecedented precision. Inferred properties for the (unresolved) underlying binary and wind system are orbital period 241.5 +/- 0.5 days and angular outflow velocity of 0.28 +/- 0.02 mas/day. An optically thin cavity of angular size 13.3 +/- 1.4 mas was found to lie between the central binary and the onset of the spiral dust plume. Rotational motion of the wind system induced by the binary orbit is found to have important ramifications: entanglement of the winds results in strong shock activity far downstream from the nose of the bowshock. The far greater fraction of the winds participating in the collision may play a key role in gas compression and the nucleation of dust at large radii from the central binary and shock stagnation point. Investigation of the effects of radiative braking pointed towards significant modifications of the simple hydrostatic colliding wind geometry, extending the relevance of this phenomena to wider binary systems than previously considered. Limits placed on the maximum allowed orbital eccentricity of e < 0.06 argue strongly for a prehistory of tidal circularization in this system. Finally we discuss the implications of Earth's polar (i < 16 deg) vantage point onto a system likely to host supernova explosions at future epochs.

The angular diameters of six oxygen rich Mira-type long-period variables have been measured at various near-infrared (NIR) wavelengths using the aperture masking technique in an extensive observing program from 1997 Jan to 2004 Sep. These data sets span many pulsation cycles of the observed objects and represent the largest study of multi-wavelength, multi-epoch interferometric angular diameter measurements on Mira stars to date. The calibrated visibility data of o Cet, R Leo, R Cas, W Hya, chi Cyg and R Hya are fitted using a uniform disk brightness distribution model to facilitate comparison between epochs, wavelengths and with existing data and theoretical models. The variation of angular diameter as a function of wavelength and time are studied, and cyclic diameter variations are detected for all objects in our sample. These variations are believed to stem from time-dependent changes of density and temperature (and hence varying molecular opacities) in different layers of these stars. The similarities and differences in behaviour between these objects are analyzed and discussed in the context of existing theoretical models. Furthermore, we present time-dependent 3.08 micron angular diameter measurements, probing for the first time these zones of probable dust formation, which show unforeseen sizes and are consistently out of phase with other NIR layers shown in this study. The S-type Mira chi Cyg exhibits significantly different behaviour compared to the M-type Miras in this study.

We report long-baseline interferometric measurements of circumstellar dust around massive evolved stars with the MIDI instrument on the Very Large Telescope Interferometer and provide spectrally dispersed visibilities in the 8-13 micron wavelength band. We also present diffraction-limited observations at 10.7 micron on the Keck Telescope with baselines up to 8.7 m which explore larger scale structure. We have resolved the dust shells around the late type WC stars WR 106 and WR 95, and the enigmatic NaSt1 (formerly WR 122), suspected to have recently evolved from a Luminous Blue Variable (LBV) stage. For AG Car, the protoypical LBV in our sample, we marginally resolve structure close to the star, distinct from the well-studied detached nebula. The dust shells around the two WC stars show fairly constant size in the 8-13 micron MIDI band, with gaussian half-widths of ~ 25 to 40 mas. The compact dust we detect around NaSt1 and AG Car favors recent or ongoing dust formation. Using the measured visibilities, we build spherically symmetric radiative transfer models of the WC dust shells which enable detailed comparison with existing SED-based models. Our results indicate that the inner radii of the shells are within a few tens of AU from the stars. In addition, our models favor grain size distributions with large (~ 1 micron) dust grains. This proximity of the inner dust to the hot central star emphasizes the difficulty faced by current theories in forming dust in the hostile environment around WR stars. Although we detect no direct evidence for binarity for these objects, dust production in a colliding-wind interface in a binary system is a feasible mechanism in WR systems under these conditions.

Spatially resolving the surfaces of nearby stars promises to advance our knowledge of stellar physics. Using optical long-baseline interferometry, we present here a near-infrared image of the rapidly rotating hot star Altair with <1 milliarcsecond resolution. The image clearly reveals the strong effect of gravity darkening on the highly-distorted stellar photosphere. Standard models for a uniformly rotating star can not explain our results, requiring differential rotation, alternative gravity darkening laws, or both.

The Mira AB system is a nearby (~107 pc) example of a wind accreting binary star system. In this class of system, the wind from a mass-losing red giant star (Mira A) is accreted onto a companion (Mira B), as indicated by an accretion shock signature in spectra at ultraviolet and X-ray wavelengths. Using novel imaging techniques, we report the detection of emission at mid-infrared wavelengths between 9.7 and 18.3 $\mu$m from the vicinity of Mira B but with a peak at a radial position about 10 AU closer to the primary Mira A. We interpret the mid-infrared emission as the edge of an optically-thick accretion disk heated by Mira A. The discovery of this new class of accretion disk fed by M-giant mass loss implies a potential population of young planetary systems in white-dwarf binaries which has been little explored, despite being relatively common in the solar neighborhood.

Located in the Orion Trapezium cluster, Theta 1 Orionis C is one of the youngest and nearest high-mass stars (O5-O7) and also known to be a close binary system. Using new multi-epoch visual and near-infrared bispectrum speckle interferometric observations obtained at the BTA 6 m telescope, and IOTA near-infrared long-baseline interferometry, we trace the orbital motion of the Theta 1 Ori C components over the interval 1997.8 to 2005.9, covering a significant arc of the orbit. Besides fitting the relative position and the flux ratio, we apply aperture synthesis techniques to our IOTA data to reconstruct a model-independent image of the Theta 1 Ori C binary system. The orbital solutions suggest a high eccentricity (e approx. 0.91) and short-period (P approx. 10.9 yrs) orbit. As the current astrometric data only allows rather weak constraints on the total dynamical mass, we present the two best-fit orbits. From these orbital solutions one can be favoured, implying a system mass of 48 M_sun and a distance to the Trapezium cluster of 434 pc. When also taking the measured flux ratio and the derived location in the HR-diagram into account, we find good agreement for all observables, assuming a spectral type of O5.5 for Theta 1 Ori C1 (M=34.0 M_sun) and O9.5 for C2 (M=15.5 M_sun). We find indications that the companion C2 is massive itself, which makes it likely that its contribution to the intense UV radiation field of the Trapezium cluster is non-negligible. Furthermore, the high eccentricity of the preliminary orbit solution predicts a very small physical separation during periastron passage (approx. 1.5 AU, next passage around 2007.5), suggesting strong wind-wind interaction between the two O stars.

Dust-grain growth and settling are the first steps towards planet formation. An understanding of dust physics is therefore integral to a complete theory of the planet formation process. In this paper, we explore the possibility of using the dust evaporation front in YSO disks (`the inner rim') as a probe of the dust physics operating in circumstellar disks. The geometry of the rim depends sensitively on the composition and spatial distribution of dust. Using radiative transfer and hydrostatic equilibrium calculations we demonstrate that dust growth and settling can curve the evaporation front dramatically (from a cylindrical radius of about 0.5 AU in the disk mid-plane to 1.2 AU in the disk upper layers for an A0 star). We compute synthetic images and interferometric visibilities for our representative rim models and show that the current generation of near-IR long-baseline interferometers (VLTI, CHARA) can strongly constrain the dust properties of circumstellar disks, shedding light on the relatively poorly understood processes of grain growth, settling and turbulent mixing.

Lambda Virginis (LamVir) is a well-known double-lined spectroscopic Am binary with the interesting property that both stars are very similar in abundance but one is sharp-lined and the other is broad-lined. We present combined interferometric and spectroscopic studies of LamVir. The small scale of the LamVir orbit (~20 mas) is well resolved by the Infrared Optical Telescope Array (IOTA), allowing us to determine its elements as well as the physical properties of the components to high accuracy. The masses of the two stars are determined to be 1.897 Msun and 1.721 Msun, with 0.7% and 1.5% errors respectively, and the two stars are found to have the same temperature of 8280 +/- 200 K. The accurately determined properties of LamVir allow comparisons between observations and current stellar evolution models, and reasonable matches are found. The best-fit stellar model gives LamVir a subsolar metallicity of Z=0.0097, and an age of 935 Myr. The orbital and physical parameters of LamVir also allow us to study its tidal evolution time scales and status. Although currently atomic diffusion is considered to be the most plausible cause of the Am phenomenon, the issue is still being actively debated in the literature. With the present study of the properties and evolutionary status of LamVir, this system is an ideal candidate for further detailed abundance analyses that might shed more light on the source of the chemical anomalies in these A stars.

We report the results of a high angular resolution near-infrared survey of dusty Wolf-Rayet stars using the Keck-1 Telescope, including new multi-wavelength images of the pinwheel nebulae WR 98a, WR 104, and WR 112. Angular sizes were measured for an additional 8 dusty WR stars using aperture masking interferometry, allowing us to probe characteristics sizes down to ~20 milliarcseconds (~40 AU for typical sources). With angular sizes and specific fluxes, we can directly measure the wavelength-dependent surface brightness and size relations for our sample. We discovered tight correlations of these properties within our sample which could not be explained by simple spherically-symmetric dust shells or even the more realistic ``pinwheel nebula'' (3-D) radiative transfer model, when using optical constants of Zubko. While the tightly-correlated surface brightness relations we uncovered offer compelling indirect evidence of a shared and distinctive dust shell geometry amongst our sample, long-baseline interferometers should target the marginally-resolved objects in our sample in order to conclusively establish the presence or absence of the putative underyling colliding wind binaries thought to produce the dust shells around WC Wolf-Rayets.

The five enigmatic Cocoon stars after which the Quintuplet cluster was christened have puzzled astronomers since their discovery. Their extraordinary cool, featureless thermal spectra have been attributed to various stellar types from young to highly evolved, while their absolute luminosities places them among the supergiants. We present diffraction-limited images from the Keck 1 telescope which resolves this debate with the discovery of rotating spiral plumes characteristic of colliding-wind binary "pinwheel" nebulae. Such elegant spiral structures, found around high-luminosity Wolf-Rayet stars, have recently been implicated in the behavior of supernovae lightcurves in the radio and optical.

We have used the Palomar 200" Adaptive Optics (AO) system to directly detect the astrometric brown dwarf GJ 802B reported by Pravdo et al. 2005. This observation is achieved with a novel combination of aperture masking interferometry and AO. The dynamical masses are 0.175$\pm$0.021 M$_\odot$ and 0.064$\pm$0.032 M$_\odot$ for the primary and secondary respectively. The inferred absolute H band magnitude of GJ 802B is M$_H$=12.8 resulting in a model-dependent T$_\mathrm{eff}$ of 1850 $\pm$ 50K and mass range of 0.057--0.074 M$_\odot$.

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.

We have measured non-zero closure phases for about 29% of our sample of 56 nearby Asymptotic Giant Branch (AGB) stars, using the 3-telescope Infrared Optical Telescope Array (IOTA) interferometer at near-infrared wavelengths (H band) and with angular resolutions in the range 5-10 milliarcseconds. These nonzero closure phases can only be generated by asymmetric brightness distributions of the target stars or their surroundings. We discuss how these results were obtained, and how they might be interpreted in terms of structures on or near the target stars. We also report measured angular sizes and hypothesize that most Mira stars would show detectable asymmetry if observed with adequate angular resolution.

We report on the detection of localized off-center emission at 1-4 AU in the circumstellar environment of the young stellar object AB Aurigae. We used closure phase measurements in the near-infrared made at the long baseline interferometer IOTA, the first obtained on a young stellar object using this technique. When probing sub-AU scales, all closure phases are close to zero degrees, as expected given the previously-determined size of the AB Aurigae inner dust disk. However, a clear closure phase signal of -3.5 +/- 0.5 degrees is detected on one triangle containing relatively short baselines, requiring a high degree of non-point symmetry from emission at larger (AU-sized) scales in the disk. We have not identified any alternative explanation for these closure phase results and demonstrate that a ``disk hot spot'' model can fit our data. We speculate that such asymmetric near-infrared emission detected might arise as a result of localized viscous heating due to a gravitational instability in the AB Aurigae disk, or to the presence of a close stellar companion or accreting sub-stellar object.

Using the 3-telescope IOTA interferometer on Mt. Hopkins, we report results from the first near-infrared (lambda=1.65 mu) closure-phase survey of Young Stellar Objects (YSOs). These closure phases allow us to unambiguously detect departures from centrosymmetry (i.e., skew) in the emission pattern from YSO disks on the scale of ~4 milliarcseconds, expected from generic ``flared disk'' models. Six of fourteen targets showed small, yet statistically-significant, non-zero closure phases, with largest values from the young binary system MWC 361-A and the (pre-main sequence?) Be star HD 45677. Our observations are quite sensitive to the vertical structure of the inner disk and we confront the predictions of the ``puffed-up inner wall'' models of Dullemond, Dominik, and Natta (DDN). Our data support disks models with curved inner rims because the expected emission appear symmetrically-distributed around the star over a wide range of inclination angles. In contrast, our results are incompatible with the models possessing vertical inner walls because they predict extreme skewness (i.e., large closure phases) from the near-IR disk emission that is not seen in our data. In addition, we also present the discovery of mysterious H-band ``halos'' (~5-10% of light on scales 0.01-0.50 arcsec) around a few objects, a preliminary ``parametric imaging'' study for HD 45677, and the first astrometric orbit for the young binary MWC 361-A.

We review recent advances in our understanding of the innermost regions of the circumstellar environment around young stars, made possible by the technique of long baseline interferometry at infrared wavelengths. Near-infrared observations directly probe the location of the hottest dust. The characteristic sizes found are much larger than previously thought, and strongly correlate with the luminosity of the central young stars. This relation has motivated in part a new class of models of the inner disk structure. The first mid-infrared observations have probed disk emission over a larger range of scales, and spectrally resolved interferometry has for the first time revealed mineralogy gradients in the disk. These new measurements provide crucial information on the structure and physical properties of young circumstellar disks, as initial conditions for planet formation.

We present new K-band long baseline interferometer observations of three young stellar objects of the FU Orionis class, V1057 Cyg, V1515 Cyg and Z CMa-SE, obtained at the Keck Interferometer during its commissioning science period. The interferometer clearly resolves the source of near-infrared emission in all three objects. Using simple geometrical models we derive size scales (0.5-4.5 AU) for this emission. All three objects appear significantly more resolved than expected from simple models of accretion disks tuned to fit the broadband optical and infrared spectro-photometry. We explore variations in the key parameters that are able to lower the predicted visibility amplitudes to the measured levels, and conclude that accretion disks alone do not reproduce the spectral energy distributions and K-band visibilities simultaneously. We conclude that either disk models are inadequate to describe the near-infrared emission, or additional source components are needed. We hypothesize that large scale emission (10s of AU) in the interferometer field of view is responsible for the surprisingly low visibilities. This emission may arise in scattering by large envelopes believed to surround these objects.

We present observations of the T Tauri stars BP Tau, DG Tau, DI Tau, GM Aur, LkCa 15, RW Aur and V830 Tau, using long baseline infrared interferometry at K band (2.2 microns) from the Keck Interferometer. The target sources have a range of mass accretion rates and excess near-infrared emission. The interferometer is most sensitive to extended emission on characteristic size scales of 1 to 5 millarcseconds. All sources show evidence for resolved K band emission on these scales, although a few of the sources are marginally consistent with being unresolved. We calculate the infrared excess based on fitting stellar photosphere models to the optical photometry and estimate the physical size of the emission region using simple geometric models for the sources with a significant infrared excess. Assuming that the K band resolved emission traces the inner edge of the dust disk, we compare the measured characteristic sizes to predicted dust sublimation radii and find that the models require a range of dust sublimation temperatures and possibly optical depths within the inner rim to match the measured radii.

This paper describes the OI Exchange Format, a standard for exchanging calibrated data from optical (visible/infrared) stellar interferometers. The standard is based on the Flexible Image Transport System (FITS), and supports storage of the optical interferometric observables including squared visibility and closure phase -- data products not included in radio interferometry standards such as UV-FITS. The format has already gained the support of most currently-operating optical interferometer projects, including COAST, NPOI, IOTA, CHARA, VLTI, PTI, and the Keck Interferometer, and is endorsed by the IAU Working Group on Optical Interferometry. Software is available for reading, writing and merging OI Exchange Format files.

We present infrared aperture synthesis maps produced with the upgraded IOTA interferometer. Michelson interferograms on the close binary system Capella (Alpha Aur) were obtained in the H-band between 2002 November 12 and 16 using the IONIC3 beam combiner. With baselines of 15m < B < 38m, we were able to determine the relative position of the binary components with milliarcsecond (mas) precision and to track their movement along the approx. 14 degree arc covered by our observation run. We briefly describe the algorithms used for visibility and closure phase estimation. Three different Hybrid Mapping and Bispectrum Fitting techniques were implemented within one software framework and used to reconstruct the source brightness distribution. By dividing our data into subsets, the system could be mapped at three epochs, revealing the motion of the stars. The precise position of the binary components was also determined with model fits, which in addition revealed I_Aa/I_Ab=1.49 +/- 0.10 and apparent stellar uniform-disk (UD) diameters of Theta_Aa=8.9 +/- 0.6 mas and Theta_Ab=5.8 +/- 0.8 mas. To improve the u, v-plane coverage, we compensated this orbital motion by applying a rotation-compensating coordinate transformation. The resulting model-independent map with a beam size of 5.4 x 2.6 mas allows the resolution of the stellar surfaces of the Capella giants themselves.

We describe the fringe packet tracking system used to equalise the optical path lengths at the Infrared Optical Telescope Array (IOTA) interferometer. The measurement of closure phases requires obtaining fringes on three baselines simultaneously. This is accomplished using an algorithm based on double Fourier interferometry for obtaining the wavelength-dependent phase of the fringes and a group delay tracking algorithm for determining the position of the fringe packet. The comparison between data acquired with and without the fringe packet tracker shows about a factor 3 reduction of the error on the closure-phase measurement. The fringe packet tracker has been able so far to track fringes of signal-to-noise as low as 1.8 for stars as faint as mH=7.0.

High angular resolution images of IRC +10216 taken at various bandpasses within the near-infrared h, k & l bands are presented. The maps have the highest angular resolution yet recovered, and were reconstructed from interferometric measurements obtained at the Keck 1 telescope in 1997 December and 1998 April, forming a subset of a 7-epoch monitoring program presented earlier [paper 1]. Systematic changes with observing wavelength are found and discussed in context of present geometrical models for the circumstellar envelope. With these new high-resolution, multi-wavelength data and contemporaneous photometry, we also re-visit the hypothesis that the bright compact Core of the nebula (component `A') marks the location of the central carbon star. We find that directly measured properties of the Core (angular size, flux density, color temperature) are consistent with a reddened carbon star photosphere (line-of-sight tau=5.3).

We report the results of a sensitive K-band survey of Herbig Ae/Be disk sizes using the 85-m baseline Keck Interferometer. Targets were chosen to span the maximum range of stellar properties to probe the disk size dependence on luminosity and effective temperature. For most targets, the measured near-infrared sizes (ranging from 0.2 to 4 AU) support a simple disk model possessing a central optically-thin (dust-free) cavity, ringed by hot dust emitting at the expected sublimation temperatures (T_sub~1000-1500K). Furthermore, we find a tight correlation of disk size with source luminosity R propto L^(1/2) for Ae and late Be systems (valid over more than 2 decades in luminosity), confirming earlier suggestions based on lower-quality data. Interestingly, the inferred dust-free inner cavities of the highest luminosity sources (Herbig B0-B3 stars) are under-sized compared to predictions of the optically-thin cavity model, likely due to optically-thick gas within the inner AU.

We present the first spatially resolved observations of a sample of 23 Mira stars simultaneously measured in the near-infrared J, H and K' bands. The technique used was optical long baseline interferometry, and we present for each star visibility amplitude measurements as a function of wavelength. We also present characteristic sizes at each spectral band, obtained by fitting the measured visibilities to a simple uniform disk model. This approach reveals the general relation J diameter < H diameter < K' diameter.

We present Keck LWS images of the Orion BN/KL star forming region obtained in the first multi-wavelength study to have 0.3-0.5" resolution from 4.7 to 22 microns. The young stellar objects designated infrared source-n and radio source-I are believed to dominate the BN/KL region. We have detected extended emission from a probable accretion disk around source-n but infer a stellar luminosity on the order of only 2000 Lsun. Although source-I is believed to be more luminous, we do not detect an infrared counterpart even at the longest wavelengths. However, we resolve the closeby infrared source, IRc2, into an arc of knots ~1000 AU long at all wavelengths. Although the physical relation of source-I to IRc2 remains ambiguous, we suggest these sources mark a high density core (10^7-10^8 pc^-3 over 1000 AU) within the larger BN/KL star forming cluster. The high density may be a consequence of the core being young and heavily embedded. We suggest the energetics of the BN/KL region may be dominated by this cluster core rather than one or two individual sources.

We present radiative-transfer modelling of the dusty spiral Pinwheel Nebula observed around the Wolf-Rayet/OB-star binary WR104. The models are based on the three-dimensional radiative-transfer code TORUS, modified to include an adaptive mesh that allows us to adequately resolve both the inner spiral turns (sub-AU scales) and the outer regions of the nebula (distances of 10^4 AU from the central source). The spiral model provides a good fit to both the spectral energy distribution and Keck aperture masking interferometry, reproducing both the maximum entropy recovered images and the visibility curves. We deduce a dust creation rate of 8+-1 x 10^{-7} solar masses per year, corresponding to approximately 2% by mass of the carbon produced by the Wolf-Rayet star. Simultaneous modelling of the imaging and spectral data enables us to constrain both the opening-angle of the wind-wind collision interface and the dust grain size. We conclude that the dust grains in the inner part of the Pinwheel nebula are small (~100A), in agreement with theoretical predictions, although we cannot rule out the presence of larger grains (~1 micron) further from the central binary. The opening angle of the wind-wind collision interface appears to be about 40 degrees, in broad agreement with the wind parameters estimated for the central binary. We discuss the success and deficiencies of the model, and the likely benefits of applying similar techniques to the more the more complex nebulae observed around other WR/O star binaries.

We present first results of an experiment to combine data from Keck aperture masking and the Infrared-Optical Telescope Array (IOTA) to image the circumstellar environments of evolved stars with ~20 milliarcsecond resolution. The unique combination of excellent Fourier coverage at short baselines and high-quality long-baseline fringe data allows us to determine the location and clumpiness of the inner-most hot dust in the envelopes, and to measure the diameters of the underlying stars themselves. We find evidence for large-scale inhomogeneities in some dust shells and also significant deviations from uniform brightness for the photospheres of the most evolved M-stars. Deviations from spherically-symmetric mass loss in the red supergiant NML Cyg could be related to recent evidence for dynamically-important magnetic fields and/or stellar rotation. We point out that dust shell asymmetries, like those observed here, can qualitatively explain the difficulty recent workers have had in simultaneously fitting the broad-band spectral energy distributions and high-resolution spatial information, without invoking unusual dust properties or multiple distinct shells (from hypothetical ``superwinds''). This paper is the first to combine optical interferometry data from multiple facilities for imaging, and we discuss the challenges and potential for the future of this method, given current calibration and software limitations.

We report the first spatially-resolved observations of the spectroscopic binaries lambda Vir and WR 140, which includes the debut of aperture-synthesis imaging with the upgraded three-telescope IOTA interferometer. Using IONIC-3, a new integrated optics beam combiner capable of precise closure phase measurement, short observations were sufficient to extract the angular separation and orientation of each binary system and the component brightness ratio. Most notably, the underlying binary in the prototypical colliding-wind source WR 140 (WC7 + O4/5) was found to have a separation of ~13 milli-arcseconds with a position angle consistent with the images of the 2001 dust shell ejection only if the Wolf-Rayet star is fainter than the O star at 1.65 microns. We also highlight lambda Vir whose peculiar stellar properties of the Am star components will permit direct testing of current theories of tidal evolution when the full orbit is determined.

We report novel, high-angular resolution interferometric measurements that imply the near-infrared nuclear emission in NGC 4151 is unexpectedly compact. We have observed the nucleus of NGC 4151 at 2.2 microns using the two 10-meter Keck telescopes as an interferometer and find a marginally resolved source ~0.1 pc in diameter. Our measurements rule out models in which a majority of the K band nuclear emission is produced on scales larger than this size. The interpretation of our measurement most consistent with other observations is that the emission mainly originates directly in the central accretion disk. This implies that AGN unification models invoking hot, optically thick dust may not be applicable to NGC 4151.

Here I review the current state of the field of optical stellar interferometry, concentrating on ground-based work although a brief report of space interferometry missions is included. We pause both to reflect on decades of immense progress in the field as well as to prepare for a new generation of large interferometers just now being commissioned (most notably, the CHARA, Keck and VLT Interferometers). First, this review summarizes the basic principles behind stellar interferometry needed by the lay-physicist and general astronomer to understand the scientific potential as well as technical challenges of interferometry. Next, the basic design principles of practical interferometers are discussed, using the experience of past and existing facilities to illustrate important points. Here there is significant discussion of current trends in the field, including the new facilities under construction and advanced technologies being debuted. This decade has seen the influence of stellar interferometry extend beyond classical regimes of stellar diameters and binary orbits to new areas such as mapping the accretion disks around young stars, novel calibration of the Cepheid Period-Luminosity relation, and imaging of stellar surfaces. The third section is devoted to the major scientific results from interferometry, grouped into natural categories reflecting these current developments. Lastly, I consider the future of interferometry, highlighting the kinds of new science promised by the interferometers coming on-line in the next few years. I also discuss the longer-term future of optical interferometry, including the prospects for space interferometry and the possibilities of large-scale ground-based projects. Critical technological developments are still needed to make these projects attractive and affordable.

We present near-infrared high resolution observations of the young binary system Z Canis Majoris using the adaptive optics system at the Keck-II telescope. Both components are unresolved at 1.25 and 1.65 microns, although high dynamic range images reveal a previously unknown jet-like feature in the circumstellar environment. We argue that this feature probably arises from light scattered off the walls of a jet-blown cavity, and proper motion studies of this feature can probe the dynamics of the bipolar outflow. Potentially, the morphology of the dust-laden cavity walls offers a new probe of the momentum profile and collimation of bipolar winds from young stellar objects. We also derive high precision binary parameters, which when combined with historical data have allowed the first detection of orbital motion. Lastly, our observations confirm the high degree of flux variability in the system; the North-West binary component is dominant at H-band, in contrast to all previous observations.

Long-baseline optical interferometers can now detect and resolve hot dust emission thought to arise at the inner edge of circumstellar disks around young stellar objects (YSOs). We argue that the near-infrared sizes being measured are closely related to the radius at which dust is sublimated by the stellar radiation field. We consider how realistic dust optical properties and gas opacity dramatically affect the predicted location of this dust destruction radius, an exercise routinely done in other contexts but so far neglected in the analysis of near-infrared sizes of YSOs. We also present the accumulated literature of near-infrared YSO sizes in the form of a ``size-luminosity diagram'' and compare with theoretical expectations. We find evidence that large (>~ 1 micron) dust grains predominate in the inner disks of T Tauri and Herbig Ae/Be stars, under the assumption that the inner-most gaseous disks are optically-thin at visible wavelengths.

The Herbig Ae/Be star LkHalpha 101 has been imaged at high angular resolution at a number of wavelengths in the near-infrared (from 1 to 3 microns) using the Keck 1 Telescope, and also observed in the mid-infrared (11.15 microns) using the U.C. Berkeley Infrared Spatial Interferometer (ISI). The resolved circular disk with a central hole or cavity reported in Tuthill et al. (2001) is confirmed. This is consistent with an almost face-on view (inclination < 35 deg) onto a luminous pre- or early-main sequence object surrounded by a massive circumstellar disk. With a multiple-epoch study spanning almost four years, relative motion of the binary companion has been detected, together with evidence for changes in the brightness distribution of the central disk/star. The resolution of the LkHalpha 101 disk by ISI mid-infrared interferometry constitutes the first such measurement of a young stellar object in this wavelength region. The angular size was found to increase only slowly from 1.6 to 11.15 microns, inconsistent with standard power-law temperature profiles usually encountered in the literature, supporting instead models with a hot inner cavity and relatively rapid transition to a cool or tenuous outer disk. The radius of the dust-free inner cavity is consistent with a model of sublimation of dust in equilibrium with the stellar radiation field. Measurements from interferometry have been combined with published photometry enabling an investigation of the energetics and fundamental properties of this prototypical system.

The eccentric WR+O binary system WR 140 produces dust for a few months at intervals of 7.94 yrs coincident with periastron passage. We present the first resolved images of this dust shell, at binary phases ~0.039 and ~0.055, using aperture masking techniques on the Keck-I telescope to achieve diffraction-limited resolution. Proper motions of approximately 1.1 milliarcsecond per day were detected, implying a distance ~1.5 kpc from the known wind speed. The dust plume observed is not as simple as the ``pinwheel'' nebulae seen around other WR colliding wind binaries, indicating the orbital plane is highly inclined to our line-of-sight and/or the dust formation is very clumpy. Follow-up imaging in the mid-infrared and with adaptive optics is urgently required to track the dust motion further, necessary for unambiguously determining the orbital geometry which we only partially constrain here. With full knowledge of the orbital elements, these infrared images can be used to reconstruct the dust distribution along the colliding wind interface, providing a unique tool for probing the post-shock physical conditions of violent astrophysical flows.

A small number of dusty Wolf-Rayet stars have been resolved into pinwheel nebulae, defined by their ``rotating'' spiral dust shells observed in the infrared. This morphology is naturally explained by dust formation associated with colliding winds in a binary system. In order to confirm and further explore this hypothesis, we have observed the known pinwheel nebulae (WR 104 and WR 98a) as well as the suspected binary WR 112 at multiple radio wavelengths with the Very Large Array to search for non-thermal radio emission from colliding winds. The spectrum of each target is nearly flat between 5 and 22 GHz, consistent with the presence of non-thermal emission that is reduced at low frequencies by free-free absorption. This emission must lie outside the radio ``photosphere,'' leading us to estimate a lower limit to the physical size of the non-thermal emitting region that is larger than expected from current theory. Based on a radio and infrared comparison to WR 104 and 98a, we conclude that WR 112 is a likely candidate pinwheel nebula, but its temporal variability indicates an eccentric binary orbit or a pinwheel viewed nearly edge-on. A sensitive radio survey of IR-bright WRs would stringently test the hypothesis that colliding winds lie at the heart of_all_ dusty WR systems. We also discuss the effects of dust obscuration in the ultra-violet and how radio-determined mass-loss rates of pinwheel nebulae (and dusty WR stars in general) may be underestimated due to shadowing effects.

This paper discusses our ongoing efforts to characterize dust-enshrouded Wolf-Rayet (WR) stars in the radio and infrared. We have used the Very Large Array to measure the broadband radio spectrum of WR stars in suspected binary systems and discovered non-thermal emission, which is usually attributed to colliding winds. In addition, infrared imaging using aperture masking interferometry on the Keck-I telescope has resolved the dust shells around a number of WR stars with K-magnitudes brighter than ~6. Although this admittedly small study suffers from selection bias, we note that all the dust-enshrouded WR stars with radio detections show evidence for colliding winds, supporting the theory that wind compression in a binary system is necessary for efficient dust production. A consequence of this hypothesis is that virtually all WC8-10 stars must be in binaries, since most are dusty. Single-star and binary stellar evolution models will have to be modified to accommodate this observational result if confirmed.

Optical interferometers increasingly use single-mode fibers as spatial filters to convert varying wavefront distortion into intensity fluctuations which can be monitored for accurate calibration of fringe amplitudes. Here I propose using an asymmetric coupler to allow the photometric intensities of each telescope beam to be measured at the same time as the fringe visibility, but without the need for dedicated photometric outputs, which reduce the light throughput in the interferometric channels. In the read-noise limited case often encountered in the infrared, I show that a 53% improvement in signal-to-noise ratio for the visibility amplitude measurement is achievable, when compared to a balanced coupler setup with 50% photometric taps (e.g., the FLUOR experiment). In the Poisson-noise limit appropriate for visible light, the improvement is reduced to only ~8%. This scheme also reduces the cost and complexity of the beam combination since fewer components and detectors are required, and can be extended to more than two telescopes for "all-in-one" or pair-wise beam combination. Asymmetric beam combination can also be employed for monitoring scintillation and throughput variations in systems without spatial filtering.

A star forms when a cloud of dust and gas collapses. It is generally believed that this collapse first produces a flattened rotating disk, through which matter is fed onto the embryonic star at the center of the disk. When the temperature and density at the center of the star pass a critical threshold, thermonuclear fusion begins. The remaining disk, which can still contain up to 0.3 times the mass of the star, is then sculpted and eventually dissipated by the radiation and wind from the newborn star. Unfortunately this picture of the structure and evolution of the disk remains speculative because of the lack of morphological data of sufficient resolution and uncertainties regarding the underlying physical processes. Here we present resolved images of a young star, LkHa 101 in which the structure of the inner accretion disk is resolved. We find that the disk is almost face-on, with a central gap (or cavity) and a hot inner edge. The cavity is bigger than previous theoretical predictions, and we infer that the position of the inner edge is probably determined by sublimation of dust grains by direct stellar radiation, rather than by disk reprocessing or the viscous heating processes as usually assumed.

We present multi-epoch, diffraction-limited images of the nebula around the carbon star CIT 6 at 2.2 microns and 3.1 microns from aperture masking on the Keck-I telescope. The near-IR nebula is resolved into two main components, an elongated, bright feature showing time-variable asymmetry and a fainter component about 60 milliarcseconds away with a cooler color temperature. These images were precisely registered (~35 milliarcseconds) with respect to recent visible images from the Hubble Space Telescope (Trammell et al. 2000), which showed a bipolar structure in scattered light. The dominant near-IR feature is associated with the northern lobe of this scattering nebula, and the multi-wavelength dataset can be understood in terms of a bipolar dust shell around CIT 6. Variability of the near-IR morphology is qualitatively consistent with previously observed changes in red polarization, caused by varying illumination geometry due to non-uniform dust production. The blue emission morphology and polarization properties can not be explained by the above model alone, but require the presence of a wide binary companion in the vicinity of the southern polar lobe. The physical mechanisms responsible for the breaking of spherical symmetry around extreme carbon stars, such as CIT 6 and IRC+10216, remain uncertain.

Using the U.C. Berkeley Infrared Spatial Interferometer with an RF filterbank, the first interferometric observations of mid-infrared molecular absorption features of ammonia (NH_3) and silane (SiH_4) with very high spectral resolution (R ~ 100000) were made. Under the assumptions of spherical symmetry and uniform outflow, these new data permitted the molecular stratification around carbon star IRC+10216 and red supergiant VY CMa to be investigated. For IRC+10216, both ammonia and silane were found to form in the dusty outflow significantly beyond both the dust formation and gas acceleration zones. Specifically, ammonia was found to form before silane in a region of decaying gas turbulence (>~ 20 R_star), while the silane is produced in a region of relatively smooth gas flow much further from the star (>~ 80 R_star). The depletion of gas-phase SiS onto grains soon after dust formation may fuel silane-producing reactions on the grain surfaces. For VY CMa, a combination of interferometric and spectral observations suggest that NH_3 is forming near the termination of the gas acceleration phase in a region of high gas turbulence (~ 40 R_star).

The U. C. Berkeley Infrared Spatial Interferometer has measured the mid-infrared visibilities of the carbon star IRC+10216 and the red supergiant VY CMa. The dust shells around these sources have been previously shown to be time-variable, and these new data are used to probe the evolution of the dust shells on a decade time-scale, complementing contemporaneous studies at other wavelengths. Self-consistent, spherically-symmetric models at maximum and minimum light both show the inner radius of the IRC+10216 dust shell to be much larger (150 mas) than that expected from the dust condensation temperature, implying that dust production has slowed or stopped in recent years. Apparently, dust does not form every pulsational cycle (638 days), and these mid-infrared results are consistent with recent near-IR imaging which indicates little or no new dust production in the last three years (Tuthill et al 2000). Spherically symmetric models failed to fit recent VY CMa data, implying that emission from the inner dust shell is highly asymmetric and/or time-variable.

The U. C. Berkeley Infrared Spatial Interferometer has been outfitted with a filterbank system to allow interferometric observations of mid-infrared spectral lines with very high spectral resolution (R ~ 10^5). This paper describes the design, implementation, and performance of the matched 32-channel filterbank modules, and new spectral line observations of Mars and IRC+10216 are used to demonstrate their scientific capability. In addition, observing strategies are discussed for accurate calibration of fringe visibilities in spectral lines, despite strong atmospheric fluctuations encountered in the infrared.

The results of a 7-epoch interferometric imaging study, at wavelengths in the near-infrared K-band, of the carbon star IRC +10216 are presented. The use of non- and partially-redundant aperture masking techniques on the 10-m Keck-I telescope has allowed us to produce images of the innermost regions of the circumstellar dust envelope with unprecedented detail. With roughly twice the resolving power of previous work, the complex asymmetric structures reported within the central 0.5 arcsec (20 stellar radii) have been imaged at the size scale of the stellar disk itself (about 50 mas). A prominent dark lane at a position angle of approximately 120 deg is suggested to be an optically thick disk or torus of dust which could help to explain IRC +10216's well-known bipolarity at a position angle of 20 deg. Observations spanning more than a pulsational cycle (638 days) have revealed significant temporal evolution of the nebula, including the outward motion of bright knots and clumps. Registering these displacements against the compact bright core, which we tentatively identify as marking the location of the star, has allowed us to determine the apparent angular velocity at a number of points. The magnitudes of the proper motions were found to be in agreement with current estimates of the stellar distance and radial velocity. Higher outflow speeds were found for features with greater separation from the core. This is consistent with acceleration taking place over the region sampled by the measurements, however alternate interpretations are also presented. Although a number of changes of morphology were found, none were clearly interpreted as the condensation of new dust over the pulsation cycle.

We report the first use of Michelson interferometry on the Keck I telescope for diffraction-limited imaging in the near infrared JHK and L bands. By using an aperture mask located close to the f/25 secondary, the 10 m Keck primary mirror was transformed into a separate-element, multiple aperture interferometer. This has allowed diffraction-limited imaging of a large number of bright astrophysical targets, including the geometrically complex dust envelopes around a number of evolved stars. The successful restoration of these images, with dynamic ranges in excess of 200:1, highlights the significant capabilities of sparse aperture imaging as compared with more conventional filled-pupil speckle imaging for the class of bright targets considered here. In particular the enhancement of the signal-to-noise ratio of the Fourier data, precipitated by the reduction in atmospheric noise, allows high fidelity imaging of complex sources with small numbers of short-exposure images relative to speckle. Multi-epoch measurements confirm the reliability of this imaging technique and our whole dataset provides a powerful demonstration of the capabilities of aperture masking methods when utilized with the current generation of large-aperture telescopes. The relationship between these new results and recent advances in interferometry and adaptive optics is briefly discussed.

The results of a high-resolution interferometric campaign targeting the symbiotic long-period variable (LPV) R~Aqr are reported. With both near-infrared measurements on baselines out to 10m and mid-infrared data extending to 32m, we have been able to measure the characteristic sizes of regions from the photosphere of the LPV and its extended molecular atmosphere, out to the cooler circumstellar dust shell. The near-infrared data were taken using aperture masking interferometry on the Keck-I telescope and show R~Aqr to be partially resolved for wavelengths out to 2.2 microns but with a marked enlargement, possibly due to molecular opacity, at 3.1 microns. Mid-infrared interferometric measurements were obtained with the U.C. Berkeley Infrared Spatial Interferometer (ISI) operating at 11.15 microns from 1992 to 1999. Although this dataset is somewhat heterogeneous with incomplete coverage of the Fourier plane and sampling of the pulsation cycle, clear changes in the mid-infrared brightness distribution were observed, both as a function of position angle on the sky and as a function of pulsation phase. Spherically symmetric radiative transfer calculations of uniform-outflow dust shell models produce brightness distributions and spectra which partially explain the data, however limitations to this approximation are noted. Evidence for significant deviation from circular symmetry was found in the mid-infrared and more tentatively at 3.08 microns in the near-infrared, however no clear detection of binarity or of non-LPV elements in the symbiotic system is reported.

We present the first near-infrared images of dusty Wolf-Rayet star WR 98a. Aperture masking interferometry has been utilized to recover images at the diffraction-limit of the Keck-I telescope, ~<50 mas at 2.2 micron. Multi-epoch observations spanning about one year have resolved the dust shell into a ``pinwheel'' nebula, the second example of a new class of dust shell first discovered around WR 104 (Tuthill, Monnier, & Danchi 1999a). Interpreting the collimated dust outflow in terms of an interacting winds model, the binary orbital parameters and apparent wind speed are derived: a period of 565 +/- 50 days, a viewing angle of 35 +/- 6 degrees from the pole, and a wind speed of 99 +/- 23 mas/yr. This period is consistent with a possible ~588 day periodicity in the infrared light curve (Williams et al. 1995), linking the photometric variation to the binary orbit. Important implications for binary stellar evolution are discussed by identifying WR 104 and WR 98a as members of a class of massive, short-period binaries whose orbits were circularized during a previous red supergiant phase. The current component separation in each system is similar to the diameter of a red supergiant, indicating that the supergiant phase was likely terminated by Roche-lobe overflow, leading to the present Wolf-Rayet stage.

Wolf-Rayet (WR) stars are luminous massive blue stars thought to be immediate precursors to the supernova terminating their brief lives. The existence of dust shells around such stars has been enigmatic since their discovery some 30 years ago; the intense radiation field from the star should be inimical to dust survival. Although dust-creation models, including those involving interacting stellar winds from a companion star, have been put forward, high-resolution observations are required to understand this phenomena. Here we present resolved images of the dust outflow around Wolf-Rayet WR 104, obtained with novel imaging techniques, revealing detail on scales corresponding to about 40 AU at the star. Our maps show that the dust forms a spatially confined stream following precisely a linear (or Archimedian) spiral trajectory. Images taken at two separate epochs show a clear rotation with a period of 220 +/- 30 days. Taken together, these findings prove that a binary star is responsible for the creation of the circumstellar dust, while the spiral plume makes WR 104 the prototype of a new class of circumstellar nebulae unique to interacting wind systems.

Recent ground-based mid-infrared spectra of 29 late-type stars, most with substantial dust shells, are compared to ground-based spectra of these stars from the 1960s and 1970s and to IRAS-LRS spectra obtained in 1983. The spectra of about half the stars show no detectable changes, implying that their distributions of circumstellar material and associated dust grain properties have changed little over this time interval. However, many of the stars with strong silicate features showed marked changes. In nearly all cases the silicate peak has strengthened with respect to the underlying continuum, although there is one case (VY~CMa) in which the silicate feature has almost completely disappeared. This suggests that, in general, an oxygen-rich star experiences long periods of gradual silicate feature strengthening, punctuated by relatively rare periods when the feature weakens. We discuss various mechanisms for producing the changes, favoring the slow evolution of the intrinsic dust properties (i.e., the chemical composition or grain structure). Although most IRAS spectra agree well with ground-based spectra, there are a number of cases where they fall well outside the expected range of uncertainty. In almost all such cases the slopes of the red and blue LRS spectra do not match in their region of overlap.

This paper will review the technical progress of interferometric infrared observing techniques from the first 2-element interferometer 25 years ago to the 3+ element arrays now coming into service. To date, only the Infrared Spatial Interferometer (ISI) has published separate-element interferometric data on circumstellar dust shells in the infrared and many of these scientific results will be discussed. Speckle interferometry has also evolved significantly over the last few decades as slit-scanning techniques over single-pixel detectors have largely been replaced by fast-readout of large format detector arrays. Important near-infrared and mid-infrared results derived from speckle data will also be reviewed. Until recently, two-dimensional information about circumstellar dust distributions has been sorely lacking, hence detections of dust shell asymmetries have been difficult and uncertain. New speckle observations using modern, 10-m class telescopes have yielded surprising results, demonstrating the importance of accurate closure phase information in interpreting interferometric data. These discoveries hopefully precursor those to be made from closure-phase imaging with the new generation of separate-element, interferometric arrays.

We present new observations of the red supergiant VY CMa at 1.25 micron, 1.65 micron, 2.26 micron, 3.08 micron and 4.8 micron. Two complementary observational techniques were utilized: non-redundant aperture masking on the 10-m Keck-I telescope yielding images of the innermost regions at unprecedented resolution, and adaptive optics imaging on the ESO 3.6-m telescope at La Silla attaining extremely high (~10^5) peak-to-noise dynamic range over a wide field. For the first time the inner dust shell has been resolved in the near-infrared to reveal a one-sided extension of circumstellar emission within 0.1" (~15 R_star) of the star. The line-of-sight optical depths of the circumstellar dust shell at 1.65 micron, 2.26 micron, and 3.08 micron have been estimated to be 1.86 +/- 0.42, 0.85 +/- 0.20, and 0.44 +/- 0.11. These new results allow the bolometric luminosity of VY~CMa to be estimated independent of the dust shell geometry, yielding L_star ~ 2x10^5 L_sun. A variety of dust condensations, including a large scattering plume and a bow-shaped dust feature, were observed in the faint, extended nebula up to 4" from the central source. While the origin of the nebulous plume remains uncertain, a geometrical model is developed assuming the plume is produced by radially-driven dust grains forming at a rotating flow insertion point with a rotational period between 1200-4200 years, which is perhaps the stellar rotational period or the orbital period of an unseen companion.

New multi-epoch, mid-infrared (8-13 microns) spectrophotometric observations are presented for 30~late-type stars. The observations were collected over a four year period (1994-1997), permitting an investigation of the mid-infrared spectral shape as a function of the pulsation cycle (typically 1-2 years). The spectra of stars with little excess infrared emission and those with carbon-rich dust show the least spectral variability, while stars with evidence for dusty, oxygen-rich envelopes are most likely to show discernible variations in their spectral profile. Most significantly, a large fraction of variable stars with strong 9.7 micron emission features show clear spectral profile changes which repeat from one cycle to the next. The significant sharpening of the silicate feature near maximum light can not be fully explained by heating and cooling of the circumstellar dust shell during the pulsational cycle, suggesting that the dust optical properties themselves must also be varying. In addition, the appearance of a narrow emission feature near the silicate peak for a few stars may require the production of especially ``pure'' silicate dust near maximum light. The general narrowing of the silicate feature observed may reflect the evolution of the pre-existing dirty grains whose surface impurities have been evaporated off when the grain temperature rises preceding maximum light. An improved theory of dust formation which can explain the observed changes in the grain properties around a single, pulsating star may lead to a definitive explanation for the diversity of silicate emission profiles observed amongst oxygen-rich, late-type stars.

Measurements by the U.C. Berkeley Infrared Spatial Interferometer at 11.15 micron have yielded strong evidence for multiple dust shells and/or significant asymmetric dust emission around NML Cyg. New observations reported also include multiple 8-13 micron spectra taken from 1994-1995 and N band (10.2 micron) photometry from 1980-1992. These and past measurements are analyzed and fitted to a model of the dust distribution around NML Cyg. No spherically symmetric single dust shell model is found consistent with both near- and mid-infrared observations. However, a circularly symmetric maximum entropy reconstruction of the 11 micron brightness distribution suggests a double shell model for the dust distribution. Such a model, consisting of a geometrically thin shell of intermediate optical depth ($\tau_{11 micron} \sim 1.9$) plus an outer shell ($\tau_{11 micron} \sim 0.33$), is consistent not only with the 11 micron visibility data, but also with near-infrared speckle measurements, the broadband spectrum, and the 9.7 micron silicate feature. The outer shell, or large scale structure, is revealed only by long-baseline interferometry at 11 micron, being too cold ($\sim$ 400 K) to contribute in the near-infrared and having no unambiguous spectral signature in the mid-infrared. The optical constants of Ossenkopf, Henning, & Mathis (1992) proved superior to the Draine & Lee (1984) constants in fitting the detailed shape of the silicate feature and broadband spectrum for this object. Recent observations of H$_2$O maser emission around NML Cyg by Richards, Yates, & Cohen (1996) are consistent with the location of the two dust shells and provide further evidence for the two-shell model.