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

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.

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.

We report on a high angular resolution survey of circumstellar disks around 24 northern sky Be stars. The K-band continuum survey was made using the CHARA Array long baseline interferometer (baselines of 30 to 331 m). The interferometric visibilities were corrected for the flux contribution of stellar companions in those cases where the Be star is a member of a known binary or multiple system. For those targets with good uv coverage, we used a four parameter Gaussian elliptical disk model to fit the visibilities and to determine the axial ratio, position angle, K-band photospheric flux contribution, and angular diameter of the disk major axis. For the other targets with relatively limited uv coverage, we constrained the axial ratio, inclination angle, and or disk position angle where necessary in order to resolve the degeneracy between possible model solutions. We also made fits of the ultraviolet and infrared spectral energy distributions to estimate the stellar angular diameter and infrared flux excess of each target. The mean ratio of the disk diameter (measured in K-band emission) to stellar diameter (from SED modeling) is 4.4 among the 14 cases where we reliably resolved the disk emission, a value which is generally lower than the disk size ratio measured in the higher opacity Halpha emission line. We estimated the equatorial rotational velocity from the projected rotational velocity and disk inclination for 12 stars, and most of these stars rotate close to or at the critical rotational velocity.

We report on the detection of a hot subdwarf component in the Be binary system, 59 Cygni. The spectral signature is found in cross-correlation functions of photospheric model spectra with far-ultraviolet spectra obtained by the International Ultraviolet Explorer Satellite, and we used radial velocities from the cross-correlation functions to determine a double-lined spectroscopic orbit. The individual spectra of the binary components were extracted using a Doppler tomography algorithm. The flux of the system is dominated by the rapidly rotating Be star. However, the subdwarf contributes approximately 4% of the UV flux, and its spectrum bears a strong resemblance to that of the hot sdO star BD +75 325. Based upon the appearance of the UV spectrum and the orbital elements, we present estimates for the stellar masses, radii, and temperatures. The presence of the hot companion causes excess emission from the outer part of the Be disk facing the companion. We present a set of red spectra that show the orbital phase variations of the He I 6678 emission formed in the heated region of the disk, which probably occurs near the disk outer boundary. 59 Cygni, FY Canis Majoris, and phi Persei comprise the known set of Be binaries with detected hot evolved companions, which are the stripped down remains of mass transfer. Their properties demonstrate that some fraction of Be stars were spun up through angular momentum transfer by Roche lobe overflow.

We have executed a survey of nearby, main sequence A, F, and G-type stars with the CHARA Array, successfully measuring the angular diameters of fortyfour stars with an average precision of ~ 1.5%. We present new measures of the bolometric flux, which in turn leads to an empirical determination of the effective temperature for the stars observed. In addition, these CHARA-determined temperatures, radii, and luminosities are fit to Yonsei-Yale model isochrones to constrain the masses and ages of the stars. These results are compared to indirect estimates of these quantities obtained by collecting photometry of the stars and applying them to model atmospheres and evolutionary isochrones. We find that for most cases, the models overestimate the effective temperature by ~ 1.5-4%, when compared to our directly measured values. The overestimated temperatures and underestimated radii in these works appear to cause an additional offset in the star's surface gravity measurements, which consequently yield higher masses and younger ages, in particular for stars with masses greater than ~1.3 Msun. Additionally, we compare our measurements to a large sample of eclipsing binary stars, and excellent agreement is seen within both data sets. Finally, we present temperature relations with respect to (B-V) and (V-K) color as well as spectral type showing that calibration of effective temperatures with errors ~ 1% is now possible from interferometric angular diameters of stars.

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 an analysis of the near-infrared continuum emission from the circumstellar gas disks of Be stars using a radiative transfer code for a parametrized version of the viscous decretion disk model. This isothermal gas model creates predicted images that we use to estimate the HWHM emission radius along the major axis of the projected disk and the spatially integrated flux excess at wavelengths of 1.7, 2.1, 4.8, 9, and 18 ?m. We discuss in detail the effect of the disk base density, inclination angle, stellar effective temperature, and other physical parameters on the derived disk sizes and color excesses. We calculate color excess estimates relative to the stellar V -band flux for a sample of 130 Be stars using photometry from 2MASS and the AKARI infrared camera all-sky survey. The color excess relations from our models make a good match of the observed color excesses of Be stars. We also present our results on the projected size of the disk as a function of wavelength for the classical Be star ? Tauri, and we show that the model predictions are consistent with interferometric observations in the H, K', and 12 \mu m bands.

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.

The emission-line Be star HD 215227 lies within the positional error circle of the newly identified gamma-ray source AGL J2241+4454. We present new blue spectra of the star, and we point out the morphological and variability similarities to other Be binaries. An analysis of the available optical photometry indicates a variation with a period of 60.37 +/- 0.04 d, which may correspond to an orbital modulation of the flux from the disk surrounding the Be star. The distance to the star of 2.6 kpc and its relatively large Galactic latitude suggest that the binary was ejected from the plane by a supernova explosion that created the neutron star or black hole companion. The binary and runaway properties of HD 215227 make it an attractive candidate as the optical counterpart of AGL J2241+4454 and as a new member of the small class of gamma-ray emitting binaries.

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.

We present spectrophotometric data from 0.4 to 4.2 microns for bright, northern sky, Be stars and several other types of massive stars. Our goal is to use these data with ongoing, high angular resolution, interferometric observations to model the density structure and sky orientation of the gas surrounding these stars. We also present a montage of the H-alpha and near-infrared emission lines that form in Be star disks. We find that a simplified measurement of the IR excess flux appears to be correlated with the strength of emission lines from high level transitions of hydrogen. This suggests that the near-IR continuum and upper level line fluxes both form in the inner part of the disk, close to the star.

Using the longest baselines of the CHARA Array, we have measured the angular diameter of the G5 V subdwarf $\mu$ Cas A, the first such determination for a halo population star. We compare this result to new diameters for the higher metallicity K0 V stars, $\sigma$ Dra and HR 511, and find that the metal-poor star, $\mu$ Cas A, has an effective temperature ($T_{\rm eff}=5297\pm32$ K), radius ($R=0.791\pm0.008 R_{\rm \odot}$), and absolute luminosity ($L=0.442\pm0.014 L_{\rm \odot}$) comparable to the other two stars with later spectral types. We show that stellar models show a discrepancy in the predicted temperature and radius for $\mu$ Cas A, and we discuss these results and how they provide a key to understanding the fundamental relationships for stars with low metallicity.

We present the first K'-band, long-baseline interferometric observations of the northern Be stars gamma Cas, phi Per, zeta Tau, and kappa Dra. The measurements were made with multiple telescope pairs of the CHARA Array interferometer, and in every case the observations indicate that the circumstellar disks of the targets are resolved. We fit the interferometric visibilities with predictions from a simple disk model that assumes an isothermal gas in Keplerian rotation. We derive fits of the four model parameters (disk base density, radial density exponent, disk normal inclination, and position angle) for each of the targets. The resulting densities are in broad agreement with prior studies of the IR excess flux and the resulting orientations generally agree with those from interferometric H-alpha and continuum polarimetric observations. We find that the angular size of the K' disk emission is smaller than that determined for the H-alpha emission, and we argue that the difference is the result of a larger H-alpha opacity and the relatively larger neutral hydrogen fraction with increasing disk radius. All the targets are known binaries with faint companions, and we find that companions appear to influence the interferometric visibilities in the cases of phi Per and kappa Dra. We also present contemporaneous observations of the H-alpha, H-gamma, and Br-gamma emission lines. Synthetic model profiles of these lines that are based on the same disk inclination and radial density exponent as derived from the CHARA Array observations match the observed emission line strength if the disk base density is reduced by approximately 1.7 dex.