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.

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

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.

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.

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.

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

Infrared interferometry of supergiant and Mira stars has recently been reinterpreted as revealing the presence of deep molecular layers. Empirical models for a photosphere surrounded by a simple molecular layer or envelope have led to a consistent interpretation of previously inconsistent data. The stellar photospheres are found to be smaller than previously understood, and the molecular layer is much higher and denser than predicted by hydrostatic equilibrium. However, the analysis was based on spatial observations with medium-band optical filters, which mixed the visibilities of different spatial structures. This paper reports spatial interferometry with narrow spectral bands, isolating near-continuum and strong molecular features, obtained for the supergiant mu Cep. The measurements confirm strong variation of apparent diameter across the K-band. A layer model shows that a stellar photosphere of angular diameter 14.11+/-0.60 mas is surrounded by a molecular layer of diameter 18.56+/-0.26 mas, with an optical thickness varying from nearly zero at 2.15 microns to >1 at 2.39 microns. Although mu Cep and alpha Ori have a similar spectral type, interferometry shows that they differ in their radiative properties. Comparison with previous broad-band measurements shows the importance of narrow spectral bands. The molecular layer or envelope appears to be a common feature of cool supergiants.

We present near-IR interferometric measurements of the K1.5 giant Arcturus (alpha Bootis), obtained at the IOTA interferometer with the FLUOR instrument, in four narrow filters with central wavelengths ranging from 2.03 to 2.39 micron. These observations were expected to allow us to quantify the wavelength dependence of the diameter of a typical K giant. They are compared to predictions from both plane-parallel and spherical model atmospheres. Unexpectedly, neither can explain the observed visibilities. We show how these data suggest the presence of a companion, in accordance with the Hipparcos data on this star, and discuss this solution with respect to Arcturus' single star status.

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.

The tip-tilt correction system at the Infrared Optical Telescope Array (IOTA) has been upgraded with a new star tracker camera. The camera features a backside-illuminated CCD chip offering doubled overall quantum efficiency and a four times higher system gain compared to the previous system. Tests carried out to characterize the new system showed a higher system gain with a lower read-out noise electron level. Shorter read-out cycle times now allow to compensate tip-tilt fluctuations so that their error imposed on visibility measurements becomes comparable to, and even smaller than, that of higher-order aberrations.

We report the observations in the K band of the red supergiant star alpha Orionis and of the bright giant star alpha Herculis with the FLUOR beamcombiner at the IOTA interferometer. The high quality of the data allows us to estimate limb-darkening and derive precise diameters in the K band which combined with bolometric fluxes yield effective temperatures. In the case of Betelgeuse, data collected at high spatial frequency although sparse are compatible with circular symmetry and there is no clear evidence for departure from circular symmetry. We have combined the K band data with interferometric measurements in the L band and at 11.15 micron. The full set of data can be explained if a 2055 K layer with optical depths $\tau_{K}=0.060\pm0.003$, $\tau_{L}=0.026\pm0.002$ and $\tau_{11.15\mu m}=2.33\pm0.23$ is added 0.33 $R_{\star}$ above the photosphere providing a first consistent view of the star in this range of wavelengths. This layer provides a consistent explanation for at least three otherwise puzzling observations: the wavelength variation of apparent diameter, the dramatic difference in limb darkening between the two supergiant stars, and the previously noted reduced effective temperature of supergiants with respect to giants of the same spectral type. Each of these may be simply understood as an artifact due to not accounting for the presence of the upper layer in the data analysis. This consistent picture can be considered strong support for the presence of a sphere of warm water vapor, proposed by Tsuji (2000) when interpreting the spectra of strong molecular lines.

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 present K'-band observations of five Mira stars with the IOTA interferometer. The interferograms were obtained with the FLUOR fiber optics beam combiner, which provides high-accuracy visibility measurements in spite of time-variable atmospheric conditions. For the M-type Miras X Oph, R Aql, RU Her, R Ser, and the C-type Mira V CrB we derived the uniform-disk diameters 11.7mas, 10.9mas, 8.4mas, 8.1mas, and 7.9mas (+/- 0.3mas), respectively. Simultaneous photometric observations yielded the bolometric fluxes. The derived angular Rosseland radii and the bolometric fluxes allowed the determination of effective temperatures. For instance, the effective temperature of R Aql was determined to be 2970 +/- 110 K. A linear Rosseland radius for R Aql of (250 +100/-60) Rsun was derived from the angular Rosseland radius of 5.5mas +/- 0.2mas and the HIPPARCOS parallax of 4.73mas +/- 1.19mas. The observations were compared with theoretical Mira star models of Bessel et al. (1996) and Hofmann et al. (1998). The effective temperatures of the M-type Miras and the linear radius of R Aql indicate fundamental mode pulsation.

Cepheids are the primary distance indicators for extragalactic astronomy and therefore are of very high astrophysical interest. Unfortunately, they are rare stars, situated very far from Earth.Though they are supergiants, their typical angular diameter is only a few milliarcseconds, making them very challenging targets even for long-baseline interferometers. We report observations that were obtained in the K prime band (2-2.3 microns), on the Cepheid Zeta Geminorum with the FLUOR beam combiner, installed at the IOTA interferometer. The mean uniform disk angular diameter was measured to be 1.64 +0.14 -0.16 mas. Pulsational variations are not detected at a significant statistical level, but future observations with longer baselines should allow a much better estimation of their amplitude. The distance to Zeta Gem is evaluated using Baade-Wesselink diameter determinations, giving a distance of 502 +/- 88 pc.

We present K-band observations of five Mira stars with the IOTA interferometer. The interferograms were obtained with the FLUOR fiber optics beam combiner which provides high-accuracy visibility measurements in spite of time-variable atmospheric conditions. For the Mira stars X Oph, R Aql, RU Her, R Ser, and V CrB we derived the uniform-disk diameters 11.7 mas, 10.9 mas, 8.4 mas, 8.1 mas, and 7.9 mas (+/-0.3 mas), respectively. Simultaneous photometric observations yielded the bolometric fluxes. The derived angular Rosseland radii and the bolometric fluxes allowed the determination of effective temperatures. For instance, the effective temperature of R Aql was determined to be 3072 K +/- 161 K. A Rosseland radius for R Aql of 250 Rsun +/- 63 Rsun was derived from the angular Rosseland radius of 5.5 mas +/- 0.2 mas and the HIPPARCOS parallax of 4.73 mas +/- 1.19 mas. The observations were compared with theoretical Mira star models of Bessel, Scholz and Wood (1996) and Hofmann, Scholz and Wood (1998).