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.

The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (http://www.univie.ac.at/diana).

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.

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.

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.

Many nearby main-sequence stars have been searched for debris using the far-infrared Herschel satellite, within the DEBRIS, DUNES and Guaranteed-Time Key Projects. We discuss here 11 stars of spectral types A to M where the stellar inclination is known and can be compared to that of the spatially-resolved dust belts. The discs are found to be well aligned with the stellar equators, as in the case of the Sun's Kuiper belt, and unlike many close-in planets seen in transit surveys. The ensemble of stars here can be fitted with a star-disc tilt of ~<10 degrees. These results suggest that proposed mechanisms for tilting the star or disc in fact operate rarely. A few systems also host imaged planets, whose orbits at tens of AU are aligned with the debris discs, contrary to what might be expected in models where external perturbers induce tilts.

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

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

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.

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.

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.

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

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.

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.

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

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.

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.

We show star images obtained with a miniature ``densified pupil imaging interferometer'' also called a hyper-telescope. The formation of such images violates a ``golden rule of imaging interferometers'' which appeared to forbid the use of interferometric arrangements differing from a Fizeau interferometer. These produce useless images when the sub-apertures spacing is much wider than their size, owing to diffraction through the sub-apertures. The hyper-telescope arrangement solves these problems opening the way towards multi-kilometer imaging arrays in space. We experimentally obtain an intensity gain of 24 +- 3X when a densified-pupil interferometer is compared to an equivalent Fizeau-type interferometer and show images of the double star alpha Gem. The initial results presented confirm the possibility of directly obtaining high resolution and high dynamic range images in the recombined focal plane of a large interferometer if enough elements are used.

We have built an infrared beam combiner for the GI2T/REGAIN interferometer of the Observatoire de la Cote d'Azur. The beam combiner allows us to record spectrally dispersed Michelson interference fringes in the near-infrared J-, H- or K-bands. The beam combiner has the advantage that Michelson interferograms can simultaneously be recorded in about 128 different spectral channels. The tilt of the spectrally dispersed fringes is a measure of the instantaneous optical path difference. We present the optical design of the beam combiner and GI2T/REGAIN observations of the Mira star R Cas with this beam combiner in the spectral range of 2.00 micron - 2.18 micron (observations on 22 and 25 August 1999; variability phase 0.08; V-magnitude approx. 6; seven baselines between 12m and 24m; reference stars Vega and Beta Peg). The spectrograph of the beam combiner consists of an anamorphotic cylindrical lens system, an image plane slit, and a grism. A system of digital signal processors calculates the ensemble average power spectrum of the spectrally dispersed Michelson interferograms and the instantaneous optical path difference error in real time. From the observed R Cas visibilities at baselines 12.0m, 13.8m and 13.9m, a 2.1 micron uniform-disk diameter of 25.3mas +/-3.3mas was derived. The unusually high visibility values at baselines >16m show that the stellar surface of R Cas is more complex than previously assumed. The visibility values at baselines >16m can be explained by high-contrast surface structure on the stellar surface of R Cas or other types of unexpected center-to-limb variations. The R Cas observations were compared with theoretical Mira star models yielding a linear Rosseland radius of 276Rsun +/-66Rsun and an effective temperature of 2685K+/-238K for R Cas at phase 0.08.