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

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

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

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.

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

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.

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.

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

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.

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.

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

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

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

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 present velocity-resolved reverberation results for five active galactic nuclei. We recovered velocity-delay maps using the maximum-entropy method for four objects: Mrk 335, Mrk 1501, 3C120, and PG2130+099. For the fifth, Mrk 6, we were only able to measure mean time delays in different velocity bins of the Hbeta emission line. The four velocity-delay maps show unique dynamical signatures for each object. For 3C120, the Balmer lines show kinematic signatures consistent with both an inclined disk and infalling gas, but the HeII 4686 emission line is suggestive only of inflow. The Balmer lines in Mrk 335, Mrk 1501, and PG 2130+099 show signs of infalling gas, but the HeII emission in Mrk 335 is consistent with an inclined disk. We also see tentative evidence of combined virial motion and infalling gas from the velocity-binned analysis of Mrk 6. The maps for 3C120 and Mrk 335 are two of the most clearly defined velocity-delay maps to date. These maps constitute a large increase in the number of objects for which we have resolved velocity-delay maps and provide evidence supporting the reliability of reverberation-based black hole mass measurements.

The last part of SpS5 dealt with the circumstellar environment. Structures are indeed found around several types of massive stars, such as blue and red supergiants, as well as WRs and LBVs. As shown in the last years, the potential of IR for their study is twofold: first, IR can help discover many previously unknown nebulae, leading to the identification of new massive stars as their progenitors; second, IR can help characterize the nebular features. Current and new IR facilities thus pave the way to a better understanding of the feedback from massive stars.

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 the results from a detailed analysis of photometric and spectrophotometric data on five Seyfert 1 galaxies observed as a part of a recent reverberation mapping program. The data were collected at several observatories over a 140-day span beginning in 2010 August and ending in 2011 January. We obtained high sampling-rate light curves for Mrk 335, Mrk 1501, 3C120, Mrk 6, and PG2130+099, from which we have measured the time lag between variations in the 5100 Angstrom continuum and the H-beta broad emission line. We then used these measurements to calculate the mass of the supermassive black hole at the center of each of these galaxies. Our new measurements substantially improve previous measurements of MBH and the size of the broad line-emitting region for four sources and add a measurement for one new object. Our new measurements are consistent with photoionization physics regulating the location of the broad line region in active galactic nuclei.

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.

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

We present the first results from a detailed analysis of photometric and spectrophotometric data on the narrow-line Seyfert 1 galaxy Mrk 335, collected over a 120-day span in the fall of 2010. From these data we measure the lag in the He II 4686 broad emission line relative to the optical continuum to be 2.7 \pm 0.6 days and the lag in the H\beta 4861 broad emission line to be 13.9 \pm 0.9 days. Combined with the line width, the He II lag yields a black hole mass, MBH = (2.6 \pm 0.8)\times 10^7 Msun. This measurement is consistent with measurements made using the H\beta 4861 line, suggesting that the He II emission originates in the same region as H\beta, but at a much smaller radius. This constitutes the first robust lag measurement for a high-ionization line in a narrow-line Seyfert 1 galaxy.

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

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

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

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

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.

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.

A technique for calculation of variability coherence at different timescales performed directly in the time domain is introduced. Simulations are made to compare the coherence spectrum derived by the time domain technique with the coherence function by Fourier analysis. The results indicate that in comparison with the Fourier analysis the time domain technique is more powerful for revealing signal coherence in noisy data. We apply the time domain technique to the real data of the black hole binaries Cygnus X-1 and GX 339-4 and compare the results with their Fourier coherence spectra.