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.

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.

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

We present the first statistical analysis of exoplanet direct imaging surveys combining adaptive optics imaging at small separations with deep seeing-limited observations at large separations allowing us to study the entire orbital separation domain from 5 to 5000~au simultaneously. Our sample of 344 stars includes only confirmed members of nearby young associations and is based on all AO direct-imaging detection limits readily available online, with addition of our own previous seeing limited surveys. Assuming that the companion distribution in mass and semi-major axis follows a power law distribution and adding a dependence on the mass of the host star, such as $d^2n\propto fM^{\alpha}a^{\beta} (M_\star/M_{\odot})^{\gamma}$d$ M $d$a$, we constrain the parameters to obtained $\alpha=-0.18^{+0.77}_{-0.65}$, $\beta=-1.43^{+0.23}_{-0.24}$, and $\gamma=0.62^{+0.56}_{-0.50}$,at a 68\% confidence level, and we obtain $f=0.11^{+0.11}_{-0.05}$, for the overall planet occurrence rate for companions with masses between 1 to 20~\mj\ in the range 5--5000~au. Thus, we find that occurrence of companions is negatively correlated with semi-major axis and companion mass (marginally) but is positively correlated with the stellar host mass. Our inferred mass distribution is in good agreement with other distributions found previously from direct imaging surveys for planets and brown dwarfs, but is shallower as a function of mass than the distributions inferred by radial velocity surveys of gas giants in the 1--3\,au range. This may suggest that planets at these wide and very-wide separations represent the low-mass tail of the brown dwarfs and stellar companion distribution rather than an extension of the distribution of the inner planets.

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.

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.

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

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u-v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction.

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

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, 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 report results from the Wide-orbit Exoplanet search with InfraRed Direct imaging (WEIRD), a survey designed to search for Jupiter-like companions on very wide orbits (1000 to 5000 AU) around young stars ($<$120 Myr) that are known members of moving groups in the solar neighborhood ($<$70 pc). Sharing the same age, distance, and metallicity as their host while being on large enough orbits to be studied as "isolated" objects make such companions prime targets for spectroscopic observations and valuable benchmark objects for exoplanet atmosphere models. The search strategy is based on deep imaging in multiple bands across the near-infrared domain. For all 177 objects of our sample, $z_{ab}^\prime$, $J$, [3.6] and [4.5] images were obtained with CFHT/MegaCam, GEMINI/GMOS, CFHT/WIRCam, GEMINI/Flamingos-2, and $Spitzer$/IRAC. Using this set of 4 images per target, we searched for sources with red $z_{ab}^\prime$ and $[3.6]-[4.5]$ colors, typically reaching good completeness down to 2Mjup companions, while going down to 1Mjup for some targets, at separations of $1000-5000$ AU. The search yielded 4 candidate companions with the expected colors, but they were all rejected through follow-up proper motion observations. Our results constrain the occurrence of 1-13 Mjup planetary-mass companions on orbits with a semi-major axis between 1000 and 5000 AU at less than 0.03, with a 95\% confidence level.

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.

The ability to observe extrasolar planets transiting their stars has profoundly changed our understanding of these planetary systems. However, these measurements depend on how well we understand the properties of the host star, such as radius, luminosity and limb darkening. Traditionally, limb darkening is treated as a parameterization in the analysis, but these simple parameterizations are not accurate representations of actual center-to-limb intensity variations (CLIV) to the precision needed for interpreting these transit observations. This effect leads to systematic errors for the measured planetary radii and corresponding measured spectral features. We compute synthetic planetary transits using model stellar atmosphere CLIV and corresponding best-fit limb-darkening laws for a grid spherically symmetric model stellar atmospheres. From these light curves we measure the differences in flux as a function of the star's effective temperature, gravity, mass, and the inclination of the planet's orbit.

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.

Context. Red supergiant stars are one of the latest stages in the evolution of massive stars. Their photospheric convection may play an important role in the launching mechanism of their mass loss; however, its characteristics and dynamics are still poorly constrained. Aims. By observing red supergiant stars with near infrared interferometry at different epochs, we expect to reveal the evolution of bright convective features on their stellar surface. Methods. We observed the M2Iab-Ib red supergiant star CE Tau with the VLTI/PIONIER instrument in the H band at two different epochs separated by one month. Results. We derive the angular diameter of the star and basic stellar parameters, and reconstruct two reliable images of its H-band photosphere. The contrast of the convective pattern of the reconstructed images is $5 \pm 1 \%$ and $6 \pm 1 \%$ for our two epochs of observation. Conclusions. The stellar photosphere shows few changes between the two epochs. The contrast of the convective pattern is below the average contrast variations obtained on 30 randomly chosen snapshots of the best matching 3D radiative hydrodynamics simulation: $23 \pm 1~\%$ for the original simulation images and $16 \pm 1~\%$ for the maps degraded to the reconstruction resolution. We offer two hypotheses to explain this observation. CE Tau may be experiencing a quiet convective activity episode or it could be a consequence of its warmer effective temperature (hence its smaller radius) compared to the simulation.

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.

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.

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.

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.

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.

We present near-infrared and optical emission-line and stellar kinematics of the Seyfert 2 galaxy Mrk 573 using the Near-Infrared Field Spectrograph (NIFS) at Gemini North and Dual Imaging Spectrograph (DIS) at Apache Point Observatory, respectively. By obtaining full kinematic maps of the infrared ionized and molecular gas and stellar kinematics in a 700 x 2100 pc^2 circumnuclear region of Mrk 573, we find that kinematics within the Narrow-Line Region (NLR) are largely due to a combination of both rotation and in situ acceleration of material originating in the host disk. Combining these observations with large-scale, optical long-slit spectroscopy that traces ionized gas emission out to several kpcs, we find that rotation kinematics dominate the majority of the gas. We find that outflowing gas extends to distances less than 1 kpc, suggesting that outflows in Seyfert galaxies may not be powerful enough to evacuate their entire bulges.

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)

One of the great challenges in understanding stars is measuring their masses. The best methods for measuring stellar masses include binary interaction, asteroseismology and stellar evolution models, but these methods are not ideal for red giant and supergiant stars. In this work, we propose a novel method for inferring stellar masses of evolved red giant and supergiant stars using interferometric and spectrophotometric observations combined with spherical model stellar atmospheres to measure what we call the stellar mass index, defined as the ratio between the stellar radius and mass. The method is based on the correlation between different measurements of angular diameter, used as a proxy for atmospheric extension, and fundamental stellar parameters. For a given star, spectrophotometry measures the Rosseland angular diameter while interferometric observations generally probe a larger limb-darkened angular diameter. The ratio of these two angular diameters is proportional to the relative extension of the stellar atmosphere, which is strongly correlated to the star's effective temperature, radius and mass. We show that these correlations are strong and can lead to precise measurements of stellar masses.

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.

We report the discovery of 14 low-mass binary systems containing mid-M to mid-L dwarf companions with separations larger than 250 AU. We also report the independent discovery of 9 other systems with similar characteristics that were recently discovered in other studies. We have identified these systems by searching for common proper motion sources in the vicinity of known high proper motion stars, based on a cross-correlation of wide area near-infrared surveys (2MASS, SDSS, and SIMP). An astrometric follow-up, for common proper motion confirmation, was made with SIMON and/or CPAPIR at the OMM 1.6 m and CTIO 1.5 m telescopes for all the candidates identified. A spectroscopic follow-up was also made with GMOS or GNIRS at Gemini to determine the spectral types of 11 of our newly identified companions and 10 of our primaries. Statistical arguments are provided to show that all of the systems we report here are very likely to be physical binaries. One of the new systems reported features a brown dwarf companion: LSPM J1259+1001 (M5) has an L4.5 (2M1259+1001) companion at about 340 AU. This brown dwarf was previously unknown. Seven other systems have a companion of spectral type L0-L1 at a separation in the 250-7500 AU range. Our sample includes 14 systems with a mass ratio below 0.3.

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.

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.

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.

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.

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.

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

We present a new method based on a Bayesian analysis to identify new members of nearby young kinematic groups. The analysis minimally takes into account the position, proper motion, magnitude and color of a star, but other observables can be readily added (e.g. radial velocity, distance). We use this method to find new young low-mass stars in the \beta Pictoris (\beta PMG) and AB Doradus (ABDMG) moving groups and in the TW Hydrae (TWA), Tucana-Horologium (THA), Columba, Carina and Argus associations. Starting from a sample of 758 mid-KM (K5V-M5V) stars showing youth indicators such as H\alpha\ and X-ray emission, our analysis yields 215 new highly probable low-mass members of the kinematic groups analyzed. One is in TWA, 37 in \beta PMG, 17 in THA, 20 in Columba, 6 in Carina, 50 in Argus, 33 in ABDMG, and the remaining 51 candidates are likely young but have an ambiguous membership to more than one association. The false alarm rate for new candidates is estimated to be 5% for \beta PMG and TWA, 10% for THA, Columba, Carina and Argus, and 14% for ABDMG. Our analysis confirms the membership of 58 stars proposed in the literature. Firm membership confirmation of our new candidates will require measurement of their radial velocity (predicted by our analysis), parallax and lithium 6708 {\AA} equivalent width. We have initiated these follow-up observations for a number of candidates and we have identified two stars (2MASSJ0111+1526, 2MASSJ0524-1601) as very strong candidate members of the \beta PMG and one strong candidate member (2MASSJ0533-5117) of the THA; these three stars have radial velocity measurements confirming their membership and lithium detections consistent with young age. Finally, we proposed that six stars should be considered as new bona fide members of \beta PMG and ABDMG, one of which being first identified in this work, the others being known candidates from the literature.

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.

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.

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.

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.

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.

We present the results of the fourth Optical/IR Interferometry Imaging Beauty Contest. The contest consists of blind imaging of test data sets derived from model sources and distributed in the OI-FITS format. The test data consists of spectral data sets on an object "observed" in the infrared with spectral resolution. There were 4 different algorithms competing this time: BSMEM the Bispectrum Maximum Entropy Method by Young, Baron & Buscher; RPR the Recursive Phase Reconstruction by Rengaswamy; SQUEEZE a Markov Chain Monte Carlo algorithm by Baron, Monnier & Kloppenborg; and, WISARD the Weak-phase Interferometric Sample Alternating Reconstruction Device by Vannier & Mugnier. The contest model image, 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.

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.

Context. The star HD 87643, exhibiting the "B[e] phenomenon", has one of the most extreme infrared excesses for this object class. It harbours a large amount of both hot and cold dust, and is surrounded by an extended reflection nebula. Aims. One of our major goals was to investigate the presence of a companion in HD87643. In addition, the presence of close dusty material was tested through a combination of multi-wavelength high spatial 5Aresolution observations. Methods. We observed HD 87643 with high spatial resolution techniques, using the near-IR AMBER/VLTI interferometer with baselines ranging from 60 m to 130 m and the mid-IR MIDI/VLTI interferometer with baselines ranging from 25 m to 65 m. These observations are complemented by NACO/VLT adaptive-optics-corrected images in the K and L-bands, ESO-2.2m optical Wide-Field Imager large-scale images in the B, V and R-bands, Results. We report the direct detection of a companion to HD 87643 by means of image synthesis using the AMBER/VLTI instrument. The presence of the companion is confirmed by the MIDI and NACO data, although with a lower confidence. The companion is separated by ~ 34 mas with a roughly north-south orientation. The period must be large (several tens of years) and hence the orbital parameters are not determined yet. Binarity with high eccentricity might be the key to interpreting the extreme characteristics of this system, namely a dusty circumstellar envelope around the primary, a compact dust nebulosity around the binary system and a complex extended nebula witnessing past violent ejections.

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.

Fringe tracking in interferometers is typically analyzed with the implicit assumption that there is a single phase associated with each telescope in the array. If the telescopes have apertures significantly larger than r0 and only partial adaptive optics correction, then the phase measured by a fringe sensor may differ significantly from the "piston" component of the aperture phase. In some cases, speckle noise will cause "branch points" in the measured phase as a function of time, causing large and sudden jumps in the phase. We present simulations showing these effects in order to understand their implications for the design of fringe tracking algorithms.

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.