The star system HR 6819 was recently proposed by Rivinius et al. as the site of the nearest example of a stellar mass black hole. Their spectra show evidence of two components: a B3 III star in a 40 d orbit and a stationary B-emission line star. Based upon the orbital mass function and the lack of evidence of a spectral component with reflex orbital motion, Rivinius et al. suggested that HR 6819 is a triple system consisting of an inner B3 III star plus black hole binary orbited by a distant Be star. Here we present an alternative model based upon an examination of the H-alpha emission line in their spectra. We show that the emission line displays the small reflex motion expected for the companion in the inner orbit and shows variations in profile shape related to orbital phase. This indicates that HR 6819 is a binary system consisting of a massive Be star and a low mass companion that is the stripped down remnant of a former mass donor star in a mass transfer binary.

We present the visual orbits of two long period spectroscopic binary stars, HD 8374 and HD 24546, using interferometric observations acquired with the CHARA Array and the Palomar Testbed Interferometer. We also obtained new radial velocities from echelle spectra using the APO 3.5 m and Fairborn 2.0 m telescopes. By combining the visual and spectroscopic observations, we solve for the full, three-dimensional orbits and determine the stellar masses and distances to within 3% uncertainty. We then estimate the effective temperature and radius of each component star through Doppler tomography and spectral energy distribution analyses, in order to compare the observed stellar parameters to the predictions of stellar evolution models. For HD 8374, we find masses of M1 = 1.636 +/- 0.050 Msun and M2 = 1.587 +/- 0.049 Msun, radii of R1 = 1.84 +/- 0.05 Rsun and R2 = 1.66 +/- 0.12 Rsun, temperatures of Teff1 = 7280 +/- 110 K and Teff2 = 7280 +/- 120 K, and an estimated age of 1.0 Gyr. For HD 24546, we find masses of M1 = 1.434 +/- 0.014 Msun and M2 = 1.409 +/- 0.014 Msun, radii of R1 = 1.67 +/- 0.06 Rsun and R2 = 1.60 +/- 0.10 Rsun, temperatures of Teff1 = 6790 +/- 120 K and Teff2 = 6770 +/- 90 K, and an estimated age of 1.4 Gyr. HD 24546 is therefore too old to be a member of the Hyades cluster, despite its physical proximity to the group.

Rapid rotation is a fundamental characteristic of classical Be stars and a crucial property allowing for the formation of their circumstellar disks. Past evolution in a mass and angular momentum transferring binary system offers a plausible solution to how Be stars attained their fast rotation. Although the subdwarf remnants of mass donors in such systems should exist in abundance, only a few have been confirmed due to tight observational constraints. An indirect method of detecting otherwise hidden companions is offered by their effect on the outer parts of Be star disks, which are expected to be disrupted or truncated. In the context of the IR and radio continuum excess radiation originating in the disk, the disk truncation can be revealed by a turndown in the spectral energy distribution due to reduced radio flux levels. In this work we search for signs of spectral turndown in a sample of 57 classical Be stars with radio data, which include new data for 23 stars and the longest wavelength detections so far ($\lambda \approx$ 10\,cm) for 2 stars. We confidently detect the turndown for all 26 stars with sufficient data coverage (20 of which are not known to have close binary companions). For the remaining 31 stars, data are inconclusive as to whether the turndown is present or not. The analysis suggests that many if not all Be stars have close companions influencing their outer disks. If confirmed to be subdwarf companions, the mass transfer spin-up scenario might explain the existence of the vast majority of classical Be stars.

Rapid rotation in massive stars imposes a latitudinal variation in the mass loss from radiatively driven winds that can lead to enhanced mass loss at the poles (with little angular momentum loss) and/or equator (with maximal angular momentum loss). Here we present an examination of the stellar wind lines of the two O-type stars with the fastest known equatorial velocities, VFTS 102 ($V\sin i = 610 \pm 30$ km/s; O9:Vnnne+) and VFTS 285 ($V\sin i = 609 \pm 29$ km/s; O7.5 Vnnn) in the Large Magellanic Cloud. Ultraviolet spectra of both stars were obtained with the Hubble Space Telescope Cosmic Origins Spectrograph. The spectrum of VFTS 285 displays a fast outflow in N V and a much slower wind in Si IV, and we argue that there is a two-wind regime in which mass loss is strong at the poles (fast and tenuous wind) but dominant at the equator (slow and dense winds). These ions and wind lines are not present in the spectrum of the cooler star VFTS 102, but the double-peaked H$\alpha$ emission in its spectrum implies equatorial mass loss into a circumstellar disk. The results suggest that in the fastest rotating O-stars, most mass is lost as an equatorial outflow, promoting angular momentum loss that contributes to a spin down over time.

We present the visual orbit of the double-lined eclipsing binary, HD 185912, from long baseline interferometry with the CHARA Array. We also obtain echelle spectra from the Apache Point observatory to update the spectroscopic orbital solution and analyze new photometry from Burggraaff et al. to model the eclipses. By combining the spectroscopic and visual orbital solutions, we find component masses of M1 = 1.361 +/- 0.004 Msun and M2 = 1.331 +/- 0.004 Msun, and a distance of d = 40.75 +/- 0.30 pc from orbital parallax. From the light curve solution, we find component radii of R1 = 1.348 +/- 0.016 Rsun and R2 = 1.322 +/- 0.016 Rsun. By comparing these observed parameters to stellar evolution models, we find that HD 185912 is a young system near the zero age main sequence with an estimated age of 500 Myr.

Super-high spatial resolution observations in the infrared are now enabling major advances in our understanding of supermassive black hole systems at the centers of galaxies. Infrared interferometry, reaching resolutions of milliarcseconds to sub-milliarcseconds, is drastically changing our view of the central structure from a static to a very dynamic one by spatially resolving to the pc-scale. We are also starting to measure the dynamical structure of fast moving gas clouds around active supermassive black holes at a scale of less than a light year. With further improvements on resolution and sensitivity, we will be able to directly image the exact site of the black hole's feedback to its host galaxy, and quantify the effect of such interaction processes. Near-future high angular resolution studies will definitely advance our mass determinations for these black holes, and we might even witness the existence of binary black hole systems at the center of galaxies.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

Direct determination of fundamental stellar parameters has many profound and wide-ranging impacts throughout astrophysics. These determinations are rooted in high angular resolution observations. In particular, as long-baseline optical interferometry has matured over the past decade, increasingly large survey samples are serving to empirically ground the basic parameters of these building blocks of the universe. True imaging and improved parametric fitting are becoming routinely available, an essential component of fully characterizing stars, stellar environments, and planets these stars may host.

Evolved stars dominate galactic spectra, enrich the galactic medium, expand to change their planetary systems, eject winds of a complex nature, produce spectacular nebulae and illuminate them, and transfer material between binary companions. While doing this, they fill the HR diagram with diagnostic loops that write the story of late stellar evolution. Evolved stars sometimes release unfathomable amounts of energy in neutrinos, light, kinetic flow, and gravitational waves. During these late-life times, stars evolve complexly, with expansion, convection, mixing, pulsation, mass loss. Some processes have virtually no spatial symmetries, and are poorly addressed with low-resolution measurements and analysis. Even a "simple" question as how to model mass loss resists solution. However, new methods offer increasingly diagnostic tools. Astrometry reveals populations and groupings. Pulsations/oscillations support study of stellar interiors. Optical/radio interferometry enable 2-3d imagery of atmospheres and shells. Bright stars with rich molecular spectra and velocity fields are a ripe opportunity for imaging with high spatial and spectral resolution, giving insight into the physics and modeling of later stellar evolution.

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.

We present the visual orbit of the double-lined spectroscopic binary HD 224355 from interferometric observations with the CHARA Array, as well as an updated spectroscopic analysis using echelle spectra from the Apache Point Observatory 3.5m telescope. By combining the visual and spectroscopic orbital solutions, we find the binary components to have masses of M1 = 1.626 +/- 0.005 Msun and M2 = 1.608 +/- 0.005 Msun, and a distance of d = 63.98 +/- 0.26 pc. Using the distance and the component angular diameters found by fitting spectrophotometry from the literature to spectral energy distribution models, we estimate the stellar radii to be R1 = 2.65 +/- 0.21 Rsun and R2 = 2.47 +/- 0.23 Rsun. We then compare these observed fundamental parameters to the predictions of stellar evolution models, finding that both components are evolved towards the end of the main sequence with an estimated age of 1.9 Gyr.

We present interferometric observations of six O-type stars that were made with the Precision Astronomical Visible Observations (PAVO) beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array. The observations include multiple brackets for three targets, $\lambda$~Ori~A, $\zeta$~Oph, and 10~Lac, but there are only preliminary, single observations of the other three stars, $\xi$~Per, $\alpha$~Cam, and $\zeta$~Ori~A. The stellar angular diameters range from 0.55 milliarcsec for $\zeta$~Ori~A down to 0.11 mas for 10~Lac, the smallest star yet resolved with the CHARA Array. The rotational oblateness of the rapidly rotating star $\zeta$ Oph is directly measured for the first time. We assembled ultraviolet to infrared flux measurements for these stars, and then derived angular diameters and reddening estimates using model atmospheres and an effective temperature set by published results from analysis of the line spectrum. The model-based angular diameters are in good agreement with observed angular diameters. We also present estimates for the effective temperatures of these stars derived by setting the interferometric angular size and fitting the spectrophotometry.

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.

We present the discovery of a gap near $M_{G}\approx$10 in the main sequence on the Hertzsprung-Russell Diagram (HRD) based on measurements presented in Gaia Data Release 2 (DR2). Using an observational form of the HRD with $M_{G}$ representing luminosity and $G_{BP}-G_{RP}$ representing temperature, the gap presents a diagonal feature that dips toward lower luminosities at redder colors. The gap is seen in samples extracted from DR2 with various distances, and is not unique to the {\it Gaia} photometry --- it also appears when using near-IR photometry ($J-K_{s}$ vs $M_{K_{s}}$). The gap is very narrow ($\sim$0.05 mag) and is near the luminosity-temperature regime where M dwarf stars transition from partially to fully convective, i.e., near spectral type M3.0V. This gap provides a new feature in the H-R Diagram that hints at an underlying astrophysical cause and we propose that it is linked to the onset of full convection in M dwarfs.

Prompted by peculiar spectroscopic variability observed in SDSS/APOGEE $H$-band spectra, we monitored the Be star HD 55606 using optical spectroscopy and found that it is an exotic double-lined spectroscopic binary (SB2) consisting of a Be star and a hot, compact companion that is probably an OB subdwarf (sdOB) star. Motion of the sdOB star is traced by its impact on the strong He~I lines, observed as radial velocity ($V_{r}$) variable, double-peaked emission profiles with narrow central absorption cores. Weak He II 4686 {\AA} absorption associated with the companion star is detected in most spectra. Use of the emission peaks of low-ionization emission lines to trace the Be star $V_{r}$ and the He I lines to trace the companion star $V_{r}$ yields a circular orbital solution with a 93.8-day period and masses of $M_{\rm Be}=6.2$ $M_{\rm \odot}$ and $M_{\rm sdOB}=0.9$ $M_{\rm \odot}$ in the case of $i=80^{\circ}$. HD 55606 exhibits a variety of phase-locked variability, including the development of shell lines twice per orbit. The shell phases coincide with variation in the double emission peak separations, and both forms of variability are likely caused by a two-armed spiral density perturbation in the Be disk. The intensity ratios of the double emission peaks are also phase-locked, possibly indicating heating by the sdOB star of the side of the Be disk facing it. HD 55606 is a new member of the growing sample of Be+sdOB binaries, in which the Be star's rapid rotation and ability to form a disk can be attributed to past mass transfer.

Eclipsing binaries are important tools for studying stellar evolution and stellar interiors. Their accurate fundamental parameters are used to test evolutionary models, and systems showing apsidal motion can also be used to test the model's internal structure predictions. For this purpose, we present a photometric and spectroscopic analysis of the eclipsing binary BW Aquarii, an evolved F-type binary with slow apsidal motion. We model the K2 C3 light curve using the Eclipsing Light Curve code to determine several orbital and stellar parameters, as well as measure the eclipse times to determine updated apsidal motion parameters for the system. Furthermore, we obtain high-resolution spectra of BW Aqr using the CHIRON echelle spectrograph on the CTIO 1.5m for radial velocity analysis. We then reconstruct the spectra of each component using Doppler tomography in order to determine the atmospheric parameters. We find that both components of BW Aqr are late F-type stars with M1 = 1.365 +/- 0.008 Msun, M2 = 1.483 +/- 0.009 Msun, and R1 = 1.782 +/- 0.021 Rsun, R2 = 2.053 +/- 0.020 Rsun. We then compare these results to the predictions of several stellar evolution models, finding that the models cannot reproduce the observed properties of both components at the same age.

Stripped-envelope stars (SESs) form in binary systems after losing mass through Roche-lobe overflow. They bear astrophysical significance as sources of UV and ionizing radiation in older stellar populations and, if sufficiently massive, as stripped supernova progenitors. Binary evolutionary models predict them to be common, but only a handful of subdwarfs (i.e., SESs) with B-type companions are known. This could be the result of observational biases hindering detection, or an incorrect understanding of binary evolution. We reanalyze the well-studied post-interaction binary phi Persei. Recently, new data improved the orbital solution of the system, which contains a ~1.2 Msun SES and a rapidly rotating ~9.6 Msun Be star. We compare with an extensive grid of evolutionary models using a Bayesian approach and find initial masses of the progenitor of 7.2+/-0.4 Msun for the SES and 3.8+/-0.4 Msun for the Be star. The system must have evolved through near-conservative mass transfer. These findings are consistent with earlier studies. The age we obtain, 57+/-9 Myr, is in excellent agreement with the age of the alpha Persei cluster. We note that neither star was initially massive enough to produce a core-collapse supernova, but mass exchange pushed the Be star above the mass threshold. We find that the subdwarf is overluminous for its mass by almost an order of magnitude, compared to the expectations for a helium core burning star. We can only reconcile this if the subdwarf is in a late phase of helium shell burning, which lasts only 2-3% of the total lifetime as a subdwarf. This could imply that up to ~50 less evolved, dimmer subdwarfs exist for each system similar to phi Persei. Our findings can be interpreted as a strong indication that a substantial population of SESs indeed exists, but has so far evaded detection because of observational biases and lack of large-scale systematic searches.

There is growing evidence that some Be stars were spun up through mass transfer in a close binary system, leaving the former mass donor star as a hot, stripped-down object. There are five known cases of Be stars with hot subdwarf (sdO) companions that were discovered through International Ultraviolet Explorer (IUE) spectroscopy. Here we expand the search for Be+sdO candidates using archival FUV spectra from IUE. We collected IUE spectra for 264 stars and formed cross-correlation functions (CCFs) with a model spectrum for a hot subdwarf. Twelve new candidate Be+sdO systems were found, and eight of these display radial velocity variations associated with orbital motion. The new plus known Be+sdO systems have Be stars with spectral subtypes of B0 to B3, and the lack of later-type systems is surprising given the large number of cooler B-stars in our sample. We discuss explanations for the observed number and spectral type distribution of the Be+sdO systems, and we argue that there are probably many Be systems with stripped companions that are too faint for detection through our analysis.

We report the discovery of a post-mass transfer Gamma Doradus/Delta Scuti hybrid pulsator in the eclipsing binary KIC~9592855. This binary has a circular orbit, an orbital period of 1.2 days, and contains two stars of almost identical masses ($M_1=1.72M_{\odot}, M_2=1.71M_{\odot}$). However, the cooler secondary star is more evolved ($R_2=1.96R_{\odot}$) while the hotter primary is still on the zero-age-main-sequence ($R_1=1.53R_{\odot}$). Coeval models from single star evolution cannot explain the observed masses and radii, and binary evolution with mass-transfer needs to be invoked. After subtracting the binary light curve, the Fourier spectrum shows low-order pressure-mode pulsations, and more dominantly, a cluster of low-frequency gravity modes at about $2$ day$^{-1}$. These g-modes are nearly equally-spaced in period, and the period spacing pattern has a negative slope. We identify these g-modes as prograde dipole modes and find that they stem from the secondary star. The frequency range of unstable p-modes also agrees with that of the secondary. We derive the internal rotation rate of the convective core and the asymptotic period spacing from the observed g-modes. The resulting values suggest that the core and envelope rotate nearly uniformly, i.e., their rotation rates are both similar to the orbital frequency of this synchronized binary.

Eclipsing binaries are vital for directly determining stellar parameters without reliance on models or scaling relations. Spectroscopically derived parameters of detached and semi-detached binaries allow us to determine component masses that can inform theories of stellar and binary evolution. Here we present moderate resolution ground-based spectra of stars in close binary systems with and without (detected) tertiary companions observed by NASA's Kepler mission and analyzed for eclipse timing variations. We obtain radial velocities and spectroscopic orbits for five single-lined and 35 double-lined systems, and confirm one false positive eclipsing binary. For the double-lined spectroscopic binaries we also determine individual component masses and examine the mass ratio M_2/M_1 distribution, which is dominated by binaries with like-mass pairs and semi-detached classical Algol systems that have undergone mass transfer. Finally, we constrain the mass of the tertiary component for five double-lined binaries with previously detected companions.

The most powerful tests of stellar models come from the brightest stars in the sky, for which complementary techniques, such as astrometry, asteroseismology, spectroscopy, and interferometry can be combined. The K2 Mission is providing a unique opportunity to obtain high-precision photometric time series for bright stars along the ecliptic. However, bright targets require a large number of pixels to capture the entirety of the stellar flux, and bandwidth restrictions limit the number and brightness of stars that can be observed. To overcome this, we have developed a new photometric technique, that we call halo photometry, to observe very bright stars using a limited number of pixels. Halo photometry is simple, fast and does not require extensive pixel allocation, and will allow us to use K2 and other photometric missions, such as TESS, to observe very bright stars for asteroseismology and to search for transiting exoplanets. We apply this method to the seven brightest stars in the Pleiades open cluster. Each star exhibits variability; six of the stars show what are most-likely slowly pulsating B-star (SPB) pulsations, with amplitudes ranging from 20 to 2000 ppm. For the star Maia, we demonstrate the utility of combining K2 photometry with spectroscopy and interferometry to show that it is not a 'Maia variable', and to establish that its variability is caused by rotational modulation of a large chemical spot on a 10 d time scale.

We used archival International Ultraviolet Explorer (IUE) high-dispersion, short wavelength spectra data to search for evidence of the spectra of hot subdwarf companions of six rapidly rotating Be stars in binary systems. We searched for the signature of a hot companion through an analysis of the cross-correlation functions of observed and model spectra that were separated into primary and secondary components using a Doppler tomography algorithm and adopted spectroscopic orbital solutions. A positive detection of the flux from a hot companion was made for the reconstructed secondary cross-correlation function of just one target, 60 Cygni (B1 Ve). We estimate that the companion of the Be star in 60 Cygni has $T_{\rm eff} = 42 \pm 4$ kK, mass ratio $M_2/M_1 = 0.15 \pm 0.02$, and monochromatic flux ratio ${f_2}/{f_1} = 0.034 \pm 0.002$ in the spectral region near 1525 \AA. If the companions of the other target Be stars are also hot, then they must be faint and contribute less than $\approx 1\%$ of the UV flux ($<0.6 \%$ in the case of $\gamma$ Cas). We also discuss in an appendix a shell episode of Pleione (28 Tau) recorded in the IUE spectra.

KIC~8262223 is an eclipsing binary with a short orbital period ($P=1.61$ d). The {\it Kepler} light curves are of Algol-type and display deep and partial eclipses, ellipsoidal variations, and pulsations of Delta Scuti type. We analyzed the {\it Kepler} photometric data, complemented by phase-resolved spectra from the R-C Spectrograph on the 4-meter Mayall telescope at Kitt Peak National Observatory and determined the fundamental parameters of this system. The low mass and oversized secondary ($M_2=0.20M_{\odot}$, $R_2=1.31R_{\odot}$) is the remnant of the donor star that transferred most of its mass to the gainer, and now the primary star. The current primary star is thus not a normal $\delta$ Scuti star but the result of mass accretion from a lower mass progenitor. We discuss the possible evolutionary history and demonstrate with the MESA evolution code that the system can be understood as the result of non-conservative binary evolution similar to that for the formation of EL CVn type binaries. The pulsations of the primary star can be explained as radial and non-radial pressure modes. The equilibrium models from single star evolutionary tracks can match the observed mass and radius ($M_1=1.94M_{\odot}$, $R_1=1.67R_{\odot}$) but the predicted unstable modes associated with these models differ somewhat from those observed. This work presents a preliminary asteroseismic analysis of the abnormal Delta Scuti pulsators, and we discuss the need for better theoretical understanding of such pulsating mass gaining stars.

KIC 3230227 is a short period ($P\approx 7.0$ days) eclipsing binary with a very eccentric orbit ($e=0.6$). From combined analysis of radial velocities and {\it Kepler} light curves, this system is found to be composed of two A-type stars, with masses of $M_1=1.84\pm 0.18M_{\odot}$, $M_2=1.73\pm 0.17M_{\odot}$ and radii of $R_1=2.01\pm 0.09R_{\odot}$, $R_2=1.68\pm 0.08 R_{\odot}$ for the primary and secondary, respectively. In addition to an eclipse, the binary light curve shows a brightening and dimming near periastron, making this a somewhat rare eclipsing heartbeat star system. After removing the binary light curve model, more than ten pulsational frequencies are present in the Fourier spectrum of the residuals, and most of them are integer multiples of the orbital frequency. These pulsations are tidally driven, and both the amplitudes and phases are in agreement with predictions from linear tidal theory for $l=2, m=-2$ prograde modes.

We use 612 single stars with previously published trigonometric parallaxes placing them within 25 pc to evaluate parallaxes released in Gaia's first data release (DR1). We find that the Gaia parallaxes are, on average, $0.24 \pm 0.02$ mas smaller than the weighted mean trigonometric parallax values for these stars in the solar neighborhood. We also find that the offset changes with distance out to 100 pc, in the sense that the closer the star, the larger the offset. We find no systematic trends in the parallax offsets with stellar $V$ magnitude, $V-K$ color, or proper motion. We do find that the offset is roughly twice as large for stars south of the ecliptic compared to those that are north.

We present a photometric and light curve analysis of an eccentric eclipsing binary in the K2 Campaign 0 field that resides in Sh 2-252E, a young star cluster embedded in an H II region. We describe a spectroscopic investigation of the three brightest stars in the crowded aperture to identify which is the binary system. We find that none of these stars are components of the eclipsing binary system, which must be one of the fainter nearby stars. These bright cluster members all have remarkable spectra: Sh 2-252a (EPIC 202062176) is a B0.5 V star with razor sharp absorption lines, Sh 2-252b is a Herbig A0 star with disk-like emission lines, and Sh 2-252c is a pre-main sequence star with very red color.

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 a spectroscopic investigation of the Be+sdO binary system HR 2142 that is based upon large sets of ultraviolet observations from the International Ultraviolet Explorer and ground-based H-alpha observations. We measured radial velocities for the Be star component from these spectra, and computed a revised orbit. In order to search for the spectral signature of the hot subdwarf, we cross-correlated the short wavelength end of each IUE spectrum with a model hot star spectrum, and then we used the predicted Doppler shifts of the subdwarf to shift-and-add all the cross-correlation functions to the frame of the subdwarf. This merged function shows the weak signal from the spectral lines of the hot star, and a best fit is obtained with a mass ratio M_2/M_1 = 0.07 +/- 0.02, companion temperature T_{eff} > 43 +/- 5 kK, projected rotational velocity V sin i < 30 km/s, and a monochromatic flux ratio near 1170 Angstroms of f_2/f_1 > 0.009 +/- 0.001. This hot subdwarf creates a one-armed spiral, tidal wake in the disk of the Be star, and we present a circumbinary disk model that can explain the occurrence of shell absorption lines by gas enhancements that occur where gas crossing the gap created by the subdwarf strikes the disk boundaries. The faint companion of HR 2142 may be representative of a significant fraction of Be stars with undetected former mass donor companion stars.

We report on interferometric observations with the CHARA Array of two classical Wolf-Rayet stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H-band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems ($f_{\rm WR, 137} = 0.59\pm0.04$; $f_{\rm WR, 138} = 0.67\pm0.01$), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modeling of each system. Our spectral modeling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edge-on. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane.

KIC 9851944 is a short period ($P=2.16$ days) eclipsing binary in the {\it Kepler} field of view. By combining the analysis of {\it Kepler} photometry and phase resolved spectra from Kitt Peak National Observatory and Lowell Observatory, we determine the atmospheric and physical parameters of both stars. The two components have very different radii ($2.27R_{\odot}$, $3.19R_{\odot}$) but close masses ($1.76 M_{\odot}$, $1.79M_{\odot}$) and effective temperatures ($7026$K, $6902$K), indicating different evolutionary stages. The hotter primary is still on the main sequence (MS), while the cooler and larger secondary star has evolved to post-MS, burning hydrogen in a shell. A comparison with coeval evolutionary models shows that it requires solar metallicity and a higher mass ratio to fit the radii and temperatures of both stars simultaneously. Both components show $\delta$ Scuti type pulsations which we interpret as p-modes and p and g mixed modes. After a close examination of the evolution of $\delta$ Scuti pulsational frequencies, we make a comparison of the observed frequencies with those calculated from MESA/GYRE.

Eclipsing binaries serve as a valuable source of stellar masses and radii that inform stellar evolutionary models and provide insight into additional astrophysical processes. The exquisite light curves generated by space-based missions such as Kepler offer the most stringent tests to date. We use the Kepler light curve of the 4.8-day eclipsing binary KIC 5739896 with ground based optical spectra to derive fundamental parameters for the system. We reconstruct the component spectra to determine the individual atmospheric parameters, and model the Kepler photometry with the binary synthesis code ELC to obtain accurate masses and radii. The two components of KIC 5738698 are F-type stars with M1 = 1.39+/-0.04M, M2 = 1.34+/-0.06M, and R1 = 1.84+/-0.03R, R2 = 1.72+/-0.03R. We also report a small eccentricity (e < 0.0017) and unusual albedo values that are required to match the detailed shape of the Kepler light curve. Comparisons with evolutionary models indicate an approximate age of 2.3 Gyr for the system.

Many short-period binary stars have distant orbiting companions that have played a role in driving the binary components into close separation. Indirect detection of a tertiary star is possible by measuring apparent changes in eclipse times of eclipsing binaries as the binary orbits the common center of mass. Here we present an analysis of the eclipse timings of 41 eclipsing binaries observed throughout the NASA Kepler mission of long duration and precise photometry. This subset of binaries is characterized by relatively deep and frequent eclipses of both stellar components. We present preliminary orbital elements for seven probable triple stars among this sample, and we discuss apparent period changes in seven additional eclipsing binaries that may be related to motion about a tertiary in a long period orbit. The results will be used in ongoing investigations of the spectra and light curves of these binaries for further evidence of the presence of third stars.

MWC 314 is a bright candidate luminous blue variable that resides in a fairly close binary system, with an orbital period of 60.753$\pm$0.003 d. We observed MWC 314 with a combination of optical spectroscopy, broad-band ground- and space-based photometry, as well as with long baseline, near-infrared interferometry. We have revised the single-lined spectroscopic orbit and explored the photometric variability. The orbital light curve displays two minima each orbit that can be partially explained in terms of the tidal distortion of the primary that occurs around the time of periastron. The emission lines in the system are often double-peaked and stationary in their kinematics, indicative of a circumbinary disc. We find that the stellar wind or circumbinary disc is partially resolved in the K\prime-band with the longest baselines of the CHARA Array. From this analysis, we provide a simple, qualitative model in an attempt to explain the observations. From the assumption of Roche Lobe overflow and tidal synchronisation at periastron, we estimate the component masses to be M1 $\approx 5$ M$_\odot$ and M2$\approx 15$ M$_\odot$, which indicates a mass of the LBV that is extremely low. In addition to the orbital modulation, we discovered two pulsational modes with the MOST satellite. These modes are easily supported by a low-mass hydrogen-poor star, but cannot be easily supported by a star with the parameters of an LBV. The combination of these results provides evidence that the primary star was likely never a normal LBV, but rather is the product of binary interactions. As such, this system presents opportunities for studying mass-transfer and binary evolution with many observational techniques.

We report on high-resolution spectroscopy of the 2009.0 spectroscopic event of $\eta$ Carinae collected via SMARTS observations using the CTIO 1.5 m telescope and echelle spectrograph. Our observations were made almost every night over a two-month interval around the photometric minimum of $\eta$ Car associated with the periastron passage of a hot companion. The photoionizing flux of the companion and heating related to colliding winds causes large changes in the wind properties of the massive primary star. Here we present an analysis of temporal variations in a sample of spectral lines that are clearly formed in the wind of the primary star. These lines are affected by a changing illumination of the flux of the secondary star during the periastron passage. We document the sudden onset of blue-shifted absorption that occurred in most of the lines near or slightly after periastron, and we argue that these absorption components are seen when we view the relatively undisturbed wind of the foreground primary star. We present time series measurements of the net equivalent width of the wind lines and of the radial velocities of the absorption trough minima and the emission peak midpoints. Most lines decrease in emission strength around periastron, and those high excitation lines formed close to the primary exhibit a red-ward velocity excursion. We show how these trends can be explained using an illuminated hemisphere model that is based on the idea that the emission originates primarily from the side of the primary facing the hot companion.

The multiple star system of delta Orionis is one of the closest examples of a system containing a luminous O-type, bright giant star (component Aa1). It is often used as a spectral-type standard and has the highest observed X-ray flux of any hot-star binary. The main component Aa1 is orbited by two lower mass stars, faint Aa2 in a 5.7 day eclipsing binary, and Ab, an astrometric companion with an estimated period of 346 years. Generally the flux from all three stars is recorded in ground-based spectroscopy, and the spectral decomposition of the components has proved difficult. Here we present HST/STIS ultraviolet spectroscopy of delta Ori A that provides us with spatially separated spectra of Aa and Ab for the first time. We measured radial velocities for Aa1 and Ab in two observations made near the velocity extrema of Aa1. We show tentative evidence for the detection of the Aa2 component in cross-correlation functions of the observed and model spectra. We discuss the appearance of the UV spectra of Aa1 and Ab with reference to model spectra. Both stars have similar effective temperatures, but Ab is fainter and is a rapid rotator. The results will help in the interpretation of ground-based spectroscopy and in understanding the physical and evolutionary parameters of these massive stars.

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.

KOI-81 is a totally eclipsing binary discovered by the Kepler mission that consists of a rapidly rotating B-type star and a small, hot companion. The system was forged through large scale mass transfer that stripped the mass donor of its envelope and spun up the mass gainer star. We present an analysis of UV spectra of KOI-81 that were obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope that reveal for the first time the spectral features of the faint, hot companion. We present a double-lined spectroscopic orbit for the system that yields mass estimates of 2.92 M_sun and 0.19 M_sun for the B-star and hot subdwarf, respectively. We used a Doppler tomography algorithm to reconstruct the UV spectra of the components, and a comparison of the reconstructed and model spectra yields effective temperatures of 12 kK and 19 - 27 kK for the B-star and hot companion, respectively. The B-star is pulsating, and we identified a number of peaks in the Fourier transform of the light curve, including one that may indicate an equatorial rotation period of 11.5~hours. The B-star has an equatorial velocity that is 74% of the critical velocity where centrifugal and gravitational accelerations balance at the equator, and we fit the transit light curve by calculating a rotationally distorted model for the photosphere of the B-star.

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 present the results of an all-sky survey made with the Fine Guidance Sensor on Hubble Space Telescope to search for angularly resolved binary systems among the massive stars. The sample of 224 stars is comprised mainly of Galactic O- and B-type stars and Luminous Blue Variables, plus a few luminous stars in the Large Magellanic Cloud. The FGS TRANS mode observations are sensitive to detection of companions with an angular separation between 0."01 and 1."0 and brighter than $\triangle m = 5$. The FGS observations resolved 52 binary and 6 triple star systems and detected partially resolved binaries in 7 additional targets (43 of these are new detections). These numbers yield a companion detection frequency of 29% for the FGS survey. We also gathered literature results on the numbers of close spectroscopic binaries and wider astrometric binaries among the sample, and we present estimates of the frequency of multiple systems and the companion frequency for subsets of stars residing in clusters and associations, field stars, and runaway stars. These results confirm the high multiplicity fraction, especially among massive stars in clusters and associations. We show that the period distribution is approximately flat in increments of log P. We identify a number of systems of potential interest for long term orbital determinations, and we note the importance of some of these companions for the interpretation of the radial velocities and light curves of close binaries that have third companions.

We present basic observational data and association membership analysis for 45 young and active low-mass stellar systems from the ongoing RECONS photometry and astrometry program at the Cerro Tololo Inter-American Observatory. Most of these systems have saturated X-ray emission (log(Lx/Lbol) > -3.5) based on X-ray fluxes from the ROSAT All-Sky Survey, and many are significantly more luminous than main-sequence stars of comparable color. We present parallaxes and proper motions, Johnson-Kron-Cousins VRI photometry, and multiplicity observations from the CTIOPI program on the CTIO 0.9m telescope. To this we add low-resolution optical spectroscopy and line measurements from the CTIO 1.5m telescope, and interferometric binary measurements from the Hubble Space Telescope Fine Guidance Sensors. We also incorporate data from published sources: JHKs photometry from the 2MASS point source catalog; X-ray data from the ROSAT All-Sky Survey; and radial velocities from literature sources. Within the sample of 45 systems, we identify 21 candidate low-mass pre-main-sequence members of nearby associations, including members of beta Pictoris, TW Hydrae, Argus, AB Doradus, two ambiguous 30 Myr old systems, and one object that may be a member of the Ursa Major moving group. Of the 21 candidate young systems, 14 are newly identified as a result of this work, and six of those are within 25 parsecs of the Sun.

We use near-infrared interferometric data coupled with trigonometric parallax values and spectral energy distribution fitting to directly determine stellar radii, effective temperatures, and luminosities for the exoplanet host stars 61 Vir, $\rho$ CrB, GJ 176, GJ 614, GJ 649, GJ 876, HD 1461, HD 7924, HD 33564, HD 107383, and HD 210702. Three of these targets are M dwarfs. Statistical uncertainties in the stellar radii and effective temperatures range from 0.5% -- 5% and from 0.2% -- 2%, respectively. For eight of these targets, this work presents the first directly determined values of radius and temperature; for the other three, we provide updates to their properties. The stellar fundamental parameters are used to estimate stellar mass and calculate the location and extent of each system's circumstellar habitable zone. Two of these systems have planets that spend at least parts of their respective orbits in the system habitable zone: two of GJ 876's four planets and the planet that orbits HD 33564. We find that our value for GJ 876's stellar radius is more than 20% larger than previous estimates and frequently used values in the astronomical literature.

We present results of a high angular resolution survey of massive OB stars in the Cygnus OB2 association that we conducted with the Fine Guidance Sensor 1R (FGS1r) on the Hubble Space Telescope. FGS1r is able to resolve binary systems with a magnitude difference delta-V < 4 down to separations as small as 0.01 arcsec. The sample includes 58 of the brighter members of Cyg OB2, one of the closest examples of an environment containing a large number of very young and massive stars. We resolved binary companions for 12 targets and confirmed the triple nature of one other target, and we offer evidence of marginally resolved companions for two additional stars. We confirm the binary nature of 11 of these systems from complementary adaptive optics imaging observations. The overall binary frequency in our study is 22% to 26% corresponding to orbital periods ranging from 20 - 20,000 years. When combined with the known short-period spectroscopic binaries, the results supports the hypothesis that the binary fraction among massive stars is > 60%. One of the new discoveries is a companion to the hypergiant star MT 304 = Cyg OB2-12, and future measurements of orbital motion should provide mass estimates for this very luminous star.

We report the discovery of a nearby M5.0V dwarf at 24.6 pc, SCR 1848-6855, that is orbited by an unusual companion causing an astrometric perturbation of more than 200 milli-arcseconds. This is by far the largest perturbation found to date among more than 700 targets observed during our long-term astrometry/photometry program at the CTIO 0.9-m telescope. We present here a suite of astrometric, photometric, and spectroscopic observations of this high proper motion (~1.3"/yr) system in an effort to reveal the nature of this unusual binary. The measured near-UV and optical U band fluxes exceed those expected for comparable M5.0V stars, and excess flux is also detected in the spectral range 4000A--7000A. The elusive companion has been detected in HST-STIS+MAMA images at 1820A and 2700A, and our analysis shows that it is probably a rare, cool, white dwarf with T = 4600-5500K. Given the long-term astrometric coverage, the prospects for an accurate mass determination are excellent, although as yet we can only provide limits on the unusual companion's mass.

KIC 9406652 is a remarkable variable star in the Kepler field of view that shows both very rapid oscillations and long term outbursts in its light curve. We present an analysis of the light curve over quarters 1 to 15 and new spectroscopy that indicates that the object is a cataclysmic variable with an orbital period of 6.108 hours. However, an even stronger signal appears in the light curve periodogram for a shorter period of 5.753 hours, and we argue that this corresponds to the modulation of flux from the hot spot region in a tilted, precessing disk surrounding the white dwarf star. We present a preliminary orbital solution from radial velocity measurements of features from the accretion disk and the photosphere of the companion. We use a Doppler tomography algorithm to reconstruct the disk and companion spectra, and we also consider how these components contribute to the object's spectral energy distribution from ultraviolet to infrared wavelengths. This target offers us a remarkable opportunity to investigate disk processes during the high mass transfer stage of evolution in cataclysmic variables.

Based on CHARA Array measurements, we present the angular diameters of 23 nearby, main- sequence stars, ranging from spectral type A7 to K0, five of which are exoplanet host stars. We derive linear radii, effective temperatures, and absolute luminosities of the stars using HIPPARCOS parallaxes and measured bolometric fluxes. The new data are combined with previously published values to create an Angular Diameter Anthology of measured angular diameters to main-sequence stars (luminosity class V and IV). This compilation consists of 125 stars with diameter uncertainties of less than 5%, ranging in spectral types from A to M. The large quantity of empirical data are used to derive color-temperature relations to an assortment of color indices in the Johnson (BVRIJHK), Cousins (RI), Kron (RI), Sloan (griz), and WISE (W3W4) photometric systems. These relations have an average standard deviation of ~3% and are valid for stars with spectral types A0 to M4. To derive even more accurate relations for Sun-like stars, we also determined these temperature relations omitting early-type stars (Teff > 6750 K) that may have biased luminosity estimates because of rapid rotation; for this subset the dispersion is only ~2.5%. We find effective temperatures in agreement within a couple percent for the interferometrically characterized sample of main sequence stars compared to those derived via the infrared-flux method and spectroscopic analysis.

R144 is a WN6h star in the 30 Doradus region. It is suspected to be a binary because of its high luminosity and its strong X-ray flux, but no periodicity could be established so far. Here, we present new Xshooter multi-epoch spectroscopy of R144 obtained at the ESO Very Large Telescope (VLT). We detect variability in position and/or shape of all the spectral lines. We measure radial velocity variations with an amplitude larger than 250 km/s in NIV and NV lines. Furthermore, the NIII and NV line Doppler shifts are anti-correlated and the NIV lines show a double-peaked profile on six of our seven epochs. We thus conclude that R144 is a double-lined spectroscopic binary. Possible orbital periods range from 2 to 6 months, although a period up to one year is allowed if the orbit is highly eccentric. We estimate the spectral types of the components to be WN5-6h and WN6-7h, respectively. The high luminosity of the system (log Lbol/Lsun ~ 6.8) suggests a present-day total mass content in the range of about 200 to 300 Msun, depending on the evolutionary stage of the components. This makes R144 the most massive binary identified so far, with a total mass content at birth possibly as large as 400 Msun. We briefly discuss the presence of such a massive object 60 pc away from the R136 cluster core in the context of star formation and stellar dynamics.

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 report on new radial velocity measurements of massive stars that are either suspected binaries or lacking prior observations. This is part of a survey to identify and characterize spectroscopic binaries among O-type stars with the goal of comparing the binary fraction of field and runaway stars with those in clusters and associations. We present orbits for HDE 308813, HD 152147, HD 164536, BD-16 4826 and HDE 229232, Galactic O-type stars exhibiting single-lined spectroscopic variation. By fitting model spectra to our observed spectra we obtain estimates for effective temperature, surface gravity, and rotational velocity. We compute orbital periods and velocity semiamplitudes for each system and note the lack of photometric variation for any system. These binaries probably appear single-lined because the companions are faint and because their orbital Doppler shifts are small compared to the width of the rotationally broadened lines of the primary.

We present interferometric diameter measurements of 21 K- and M- dwarfs made with the CHARA Array. This sample is enhanced by literature radii measurements to form a data set of 33 K-M dwarfs with diameters measured to better than 5%. For all 33 stars, we compute absolute luminosities, linear radii, and effective temperatures (Teff). We develop empirical relations for \simK0 to M4 main- sequence stars between the stellar Teff, radius, and luminosity to broad-band color indices and metallicity. These relations are valid for metallicities between [Fe/H] = -0.5 and +0.1 dex, and are accurate to ~2%, ~5%, and ~4% for Teff, radius, and luminosity, respectively. Our results show that it is necessary to use metallicity dependent transformations to convert colors into stellar Teffs, radii, and luminosities. We find no sensitivity to metallicity on relations between global stellar properties, e.g., Teff-radius and Teff-luminosity. Robust examinations of single star Teffs and radii compared to evolutionary model predictions on the luminosity-Teff and luminosity-radius planes reveals that models overestimate the Teffs of stars with Teff < 5000 K by ~3%, and underestimate the radii of stars with radii < 0.7 R\odot by ~5%. These conclusions additionally suggest that the models overestimate the effects that the stellar metallicity may have on the astrophysical properties of an object. By comparing the interferometrically measured radii for single stars to those of eclipsing binaries, we find that single and binary star radii are consistent. However, the literature Teffs for binary stars are systematically lower compared to Teffs of single stars by ~ 200 to 300 K. Lastly, we present a empirically determined HR diagram for a total of 74 nearby, main-sequence, A- to M-type stars, and define regions of habitability for the potential existence of sub-stellar mass companions in each system. [abridged]

We present the first results of a Kepler survey of 41 eclipsing binaries that we undertook to search for third star companions. Such tertiaries will periodically alter the eclipse timings through light travel time and dynamical effects. We discuss the prevalence of starspots and pulsation among these binaries and how these phenomena influence the eclipse times. There is no evidence of short period companions (P < 700 d) among this sample, but we do find evidence for long term timing variations in 14 targets (34%). We argue that this finding is consistent with the presence of tertiary companions among a significant fraction of the targets, especially if many have orbits measured in decades. This result supports the idea that the formation of close binaries involves the deposition of angular momentum into the orbital motion of a third star.

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

The stellar mass-luminosity relation is poorly constrained by observations for high mass stars. We describe our program to find eclipsing massive binaries in the Magellanic Clouds using photometry of regions rich in massive stars, and our spectroscopic follow-up to obtain radial velocities and orbits. Our photometric campaign identified 48 early-type periodic variables, of which only 15 (31%) were found as part of the microlensing surveys. Spectroscopy is now complete for 17 of these systems, and in this paper we present analysis of the first two, LMC 172231 and ST2-28, simple detached systems of late-type O dwarfs of relatively modest masses. Our orbit analysis yields very precise masses (2%) and we use tomography to separate the components and determine effective temperatures by model fitting, necessary for determining accurate (0.05-0.07 dex) bolometric luminosities in combination with the light-curve analysis. Our approach allows more precise comparisons with evolutionary theory than previously possible. To our considerable surprise, we find a small, but significant, systematic discrepancy: all of the stars are slightly under-massive, by typically 11% (or over-luminous by 0.2 dex) compared to that predicted by the evolutionary models. We examine our approach for systematic problems, but find no satisfactory explanation. The discrepancy is in the same sense as the long-discussed and elusive discrepancy between the masses measured from stellar atmosphere analysis with the stellar evolutionary models, and might suggest that either increased rotation or convective overshooting is needed in the models. Additional systems will be discussed in future papers of this series, and will hopefully confirm or refute this trend.

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.

The hot star HDE 326823 is a candidate transition-phase object that is evolving into a nitrogen-enriched Wolf-Rayet star. It is also a known low-amplitude, photometric variable with a 6.123 d period. We present new, high and moderate resolution spectroscopy of HDE 326823, and we show that the absorption lines show coherent Doppler shifts with this period while the emission lines display little or no velocity variation. We interpret the absorption line shifts as the orbital motion of the apparently brighter star in a close, interacting binary. We argue that this star is losing mass to a mass gainer star hidden in a thick accretion torus and to a circumbinary disk that is the source of the emission lines. HDE 326823 probably belongs to a class of objects that produce short-period WR+O binaries.

We report on the preliminary results of a three-year monitoring campaign of Galactic and Magellanic LBVs with both moderate and high resolution spectroscopy. We have collected more than 500 moderate-resolution spectra of 20 stars in the red portion of the optical spectrum, which includes the wind-sensitive transitions of H$\alpha$ and He I 5876 and 6678 Angstroms. We summarize our detailed study of 15 years of high resolution H$\alpha$ spectroscopy of the prototypical luminous blue variable, P Cygni. We report on the discovery of discrete absorption components in P Cygni's H$\alpha$ profile, and we discuss their relationship to structure in the wind. These results are compared to our recent high resolution interferometric imaging of the $H$-band emitting region surrounding the star. We discuss recent observations of $\eta$ Car, highlighting its unusual recovery from the previous periastron passage. Our results on HDE 326823 indicate that the star has a short period binary orbit, with Roche Lobe overflow onto an unseen massive companion.

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.

The star HD 193322 is a remarkable multiple system of massive stars that lies at the heart of the cluster Collinder 419. Here we report on new spectroscopic observations and radial velocities of the narrow-lined component Ab1 that we use to determine its orbital motion around a close companion Ab2 ($P = 312$ d) and around a distant third star Aa ($P = 35$ y).We have also obtained long baseline interferometry of the target in the $K^\prime$-band with the CHARA Array that we use in two ways. First, we combine published speckle interferometric measurements with CHARA separated fringe packet measurements to improve the visual orbit for the wide Aa,Ab binary. Second, we use measurements of the fringe packet from Aa to calibrate the visibility of the fringes of the Ab1,Ab2 binary, and we analyze these fringe visibilities to determine the visual orbit of the close system. The two most massive stars, Aa and Ab1, have masses of approximately 21 and $23 M_\odot$, respectively, and their spectral line broadening indicates that they represent extremes of fast and slow projected rotational velocity, respectively.

We present a catalog of 56 candidate intermediate mass eclipsing binary systems extracted from the 3rd data release of the All Sky Automated Survey. We gather pertinent observational data and derive orbital properties, including ephemerides, for these systems as a prelude to anticipated spectroscopic observations. We find that 37 of the 56, or ~66% of the systems are not identified in the Simbad Astronomical Database as known binaries. As a specific example, we show spectroscopic data obtained for the system HI Mon (B0 V + B0.5 V) observed at key orbital phases based on the computed ephemeris, and we present a combined spectroscopic and photometric solution for the system and give stellar parameters for each component.

P Cygni is a prototype of the Luminous Blue Variables (or S Doradus variables), and the star displays photometric and emission line variability on a timescale of years (known as the "short S Doradus phase" variations). Here we present new high resolution H-alpha spectroscopy of P Cyg that we combine with earlier spectra and concurrent V-band photometry to document the emission and continuum flux variations over a 24 y time span. We show that the emission and continuum fluxes vary in concert on timescales of 1.6 y and longer, but differ on shorter timescales. The H-alpha profile shape also varies on the photometric timescales, and we describe the observed co-variations of the emission peak and absorption trough properties. We argue that the episodes of photometric and emission brightening are caused by increases in the size of the emission region that are related to variations in wind mass loss rate and outflow speed. We find evidence of blueward accelerating, Discrete Absorption Components (DACs) in the absorption trough of the H-alpha profile, and these features have slower accelerations and longer durations than those observed in other lines. The DAC strengths also appear to vary on the photometric timescales, and we suggest that the propagation of the DAC-related wind structures is closely related to changes in the overall wind mass loss rate and velocity.

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.

LS I +61 303 is an exceptionally rare example of a high mass X-ray binary (HMXB) that also exhibits MeV-TeV emission, making it one of only a handful of "gamma-ray binaries". Here we present H-alpha spectra that show strong variability during the 26.5 day orbital period and over decadal time scales. We detect evidence of a spiral density wave in the Be circumstellar disk over part of the orbit. The H-alpha line profile also exhibits a dramatic emission burst shortly before apastron, observed as a redshifted shoulder in the line profile, as the compact source moves almost directly away from the observer. We investigate several possible origins for this red shoulder, including an accretion disk, mass transfer stream, and a compact pulsar wind nebula that forms via a shock between the Be star's wind and the relativistic pulsar wind.

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.

Two recent observing campaigns provide us with moderate dispersion spectra of more than 230 cluster and 370 field B stars. Combining them and the spectra of the B stars from our previous investigations ($\sim$430 cluster and $\sim$100 field B stars) yields a large, homogeneous sample for studying the rotational properties of B stars. We derive the projected rotational velocity $V\sin i$, effective temperature, gravity, mass, and critical rotation speed $V_{\rm crit}$ for each star. We find that the average $V\sin i$ is significantly lower among field stars because they are systematically more evolved and spun down than their cluster counterparts. The rotational distribution functions of $V_{\rm eq}/V_{\rm crit}$ for the least evolved B stars show that lower mass B stars are born with a larger proportion of rapid rotators than higher mass B stars. However, the upper limit of $V_{\rm eq}/V_{\rm crit}$ that may separate normal B stars from emission line Be stars (where rotation promotes mass loss into a circumstellar disk) is smaller among the higher mass B stars. We compare the evolutionary trends of rotation (measured according to the polar gravity of the star) with recent models that treat internal mixing. The spin-down rates observed in the high mass subset ($\sim 9 M_\odot$) agree with predictions, but the rates are larger for the low mass group ($\sim 3 M_\odot$). The faster spin down in the low mass B stars matches well with the predictions based on conservation of angular momentum in individual spherical shells. Our results suggest the fastest rotators (that probably correspond to the emission line Be stars) are probably formed by evolutionary spin up (for the more massive stars) and by mass transfer in binaries (for the full range of B star masses).

We present the results of a comprehensive assessment of companions to solar-type stars. A sample of 454 stars, including the Sun, was selected from the Hipparcos catalog with {\pi} > 40 mas, {\sigma}_{\pi}/{\pi} < 0.05, 0.5 < B - V < 1.0 (~ F6-K3), and constrained by absolute magnitude and color to exclude evolved stars. New observational aspects of this work include surveys for (1) very close companions with long-baseline interferometry at the CHARA Array, (2) close companions with speckle interferometry, and (3) wide proper motion companions identified by blinking multi-epoch archival images. In addition, we include the results from extensive radial-velocity monitoring programs and evaluate companion information from various catalogs. The overall observed fractions of single, double, triple, and higher order systems are 56% \pm 2%, 33% \pm 2%, 8% \pm 1%, and 3% \pm 1%, respectively, counting all confirmed stellar and brown dwarf companions. Our completeness analysis indicates that only a few undiscovered companions remain in this well-studied sample, implying that the majority (54% \pm 2%) of solar-type stars are single, in contrast to the results of prior multiplicity studies. The orbital-period distribution of companions is unimodal and roughly log-normal with a peak of about 300 years. The period-eccentricity relation shows a roughly flat distribution beyond the expected circularization for periods below 12 days. The mass-ratio distribution shows a preference for like-mass pairs, which occur more frequently in relatively close pairs. The fraction of planet hosts among single, binary, and multiple systems are statistically indistinguishable, suggesting that planets are as likely to form around single stars as they are around components of binary or multiple systems with sufficiently wide separations.

The young open cluster Collinder 419 surrounds the massive O star, HD 193322, that is itself a remarkable multiple star system containing at least four components. Here we present a discussion of the cluster distance based upon new spectral classifications of the brighter members, UBV photometry, and an analysis of astrometric and photometric data from the UCAC3 and 2MASS catalogs. We determine an average cluster reddening of E(B-V)=0.37 +- 0.05 mag and a cluster distance of 741 +- 36 pc. The cluster probably contains some very young stars that may include a reddened M3 III star, IRAS~20161+4035.

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.

We report on H-alpha spectroscopy of the 2009.0 spectroscopic event of eta Carinae collected via SMARTS observations using the CTIO 1.5 m telescope and echelle spectrograph. Our observations were made almost every night over a two month interval around the predicted minimum of eta Car. We observed a significant fading of the line emission that reached a minimum seven days after the X-ray minimum. About 17 d prior to the H-alpha flux minimum, the H-alpha profile exhibited the emergence of a broad, P Cygni type, absorption component (near a Doppler shift of -500 km/s) and a narrow absorption component (near -144 km/s and probably associated with intervening gas from the Little Homunculus Nebula). All these features were observed during the last event in 2003.5 and are probably related to the close periastron passage of the companion. We argue that these variations are consistent with qualitative expectations about changes in the primary star's stellar wind that result from the wind-wind collision with a massive binary companion and from atmospheric eclipses of the companion.

HD 15137 is an intriguing runaway O-type binary system that offers a rare opportunity to explore the mechanism by which it was ejected from the open cluster of its birth. Here we present recent blue optical spectra of HD 15137 and derive a new orbital solution for the spectroscopic binary and physical parameters of the O star primary. We also present the first XMM-Newton observations of the system. Fits of the EPIC spectra indicate soft, thermal X-ray emission consistent with an isolated O star. Upper limits on the undetected hard X-ray emission place limits on the emission from a proposed compact companion in the system, and we rule out a quiescent neutron star in the propellor regime or a weakly accreting neutron star. An unevolved secondary companion is also not detected in our optical spectra of the binary, and it is difficult to conclude that a gravitational interaction could have ejected this runaway binary with a low mass optical star. HD 15137 may contain an elusive neutron star in the ejector regime or a quiescent black hole with conditions unfavorable for accretion at the time of our observations.

We report on the discovery of a high velocity B star, HD 69686. We estimate its space velocity, distance, surface temperature, gravity, and age. With these data, we are able to reconstruct the trajectory of the star and to trace it back to its birthplace. We use evolutionary tracks for single stars to estimate that HD 69686 was born 73 Myr ago in the outer part of our Galaxy ($r \sim 12$ kpc) at a position well below the Galactic plane ($z \sim -1.8$ kpc), a very unusual birthplace for a B star. Along the star's projected path in the sky, we also find about 12 other stars having similar proper motions, and their photometry data suggest that they are located at the same distance as HD 69686 and probably have the same age. We speculate on the origin of this group by star formation in a high velocity cloud or as a Galactic merger fragment.

We present an examination of high resolution, ultraviolet spectroscopy from Hubble Space Telescope of the photospheric spectrum of the O-supergiant in the massive X-ray binary HD 226868 = Cyg X-1. We analyzed this and ground-based optical spectra to determine the effective temperature and gravity of the O9.7 Iab supergiant. Using non-local thermodynamic equilibrium (non-LTE), line blanketed, plane parallel models from the TLUSTY grid, we obtain T_eff = 28.0 +/- 2.5kK and log g > 3.00 +/- 0.25, both lower than in previous studies. The optical spectrum is best fit with models that have enriched He and N abundances. We fit the model spectral energy distribution for this temperature and gravity to the UV, optical, and IR fluxes to determine the angular size of and extinction towards the binary. The angular size then yields relations for the stellar radius and luminosity as a function of distance. By assuming that the supergiant rotates synchronously with the orbit, we can use the radius - distance relation to find mass estimates for both the supergiant and black hole as a function of the distance and the ratio of stellar to Roche radius. Fits of the orbital light curve yield an additional constraint that limits the solutions in the mass plane. Our results indicate masses of 23^{+8}_{-6} M_sun for the supergiant and 11^{+5}_{-3} M_sun for the black hole.

We recently discovered a large number of highly active Be stars in the open cluster NGC 3766, making it an excellent location to study the formation mechanism of Be star disks. To explore whether similar disk appearances and/or disappearances are common among the Be stars in other open clusters, we present here multiple epochs of H-alpha spectroscopy for 296 stars in eight open clusters. We identify 12 new transient Be stars and confirm 17 additional Be stars with relatively stable disks. By comparing the H-alpha equivalent widths to the photometric y - H-alpha colors, we present a method to estimate the strength of the H-alpha emission when spectroscopy is not available. For a subset of 128 stars in four open clusters, we also use blue optical spectroscopy and available Stromgren photometry to measure their projected rotational velocities, effective temperatures, and polar surface gravities. We combine our Be star detections from these four clusters to investigate physical differences between the transient Be stars, stable Be stars, and normal B-type stars with no line emission. Both types of Be stars are faster rotating populations than normal B-type stars, and we find no significant physical differences between the transient and stable Be stars in our sample.

We present projected rotational velocity values for 97 Galactic, 55 SMC, and 106 LMC O-B type stars from archival FUSE observations. The evolved and unevolved samples from each environment are compared through the Kolmogorov-Smirnov test to determine if the distribution of equatorial rotational velocities is metallicity dependent for these massive objects. Stellar interior models predict that massive stars with SMC metallicity will have significantly reduced angular momentum loss on the main sequence compared to their Galactic counterparts. Our results find some support for this prediction but also show that even at Galactic metallicity, evolved and unevolved massive stars have fairly similar fractions of stars with large V sin i values. Macroturbulent broadening that is present in the spectral features of Galactic evolved massive stars is lower in the LMC and SMC samples. This suggests the processes that lead to macroturbulence are dependent upon metallicity.

We present the results of a Doppler tomographic reconstruction of the UV spectra of the double-lined, O binary HD 100213 based on observations made with the International Ultraviolet Explorer (IUE). We used cross-correlation methods to obtain radial velocities, confirm the orbital elements, estimate the UV flux ratio, and determine projected rotational velocities. The individual component spectra are classified as O7 V + O8 V using UV criteria defined by Penny, Gies, & Bagnuolo. We present a model fit of the eclipsing light curve from observations from the Hipparcos satellite and published observations of Andersen & Gronbech. We derive an orbital inclination, i=77.7deg+/-1.0deg. This analysis indicates that both stars are currently experiencing Roche lobe overflow (RLOF), which confirms earlier results that this is one of only a few massive contact binaries. Our derived masses, Mp/Msolar=16.7+/-0.4 and Ms/Msolar=10.4+/-0.4, are significantly lower than those computed from the Doppler shifts of lines in the optical spectrum. We suggest that the difference occurs because mutual irradiation decreases the upper atmospheric temperature gradient in the inward-facing hemispheres of both stars, which makes lower excitation lines appear weaker there and shifts their center of light away from the center of mass. We compare the current state of HD 100213 with predicted outcomes of massive close binary evolutionary models, and suggest that the system is currently in a very slow case AA mass transfer stage.

We present a spectroscopic and photometric analysis of the remarkable massive binary system, BD +36 4063. We argue that the visible ON star is undergoing a rapid mass transfer episode that results in a thick and opaque disk that surrounds and renders invisible its massive companion. A comparison of the projected rotational velocity and the orbital semiamplitude of the visible star indicates a mass ratio near unity. Models for conservative mass transfer show that the equal mass state occurs at the point of minimum separation, and we argue that BD +36 4063 may represent the first system identified at this rapid and rare stage of massive binary evolution.

We present the results of a speckle interferometric survey of Galactic massive stars that complements and expands upon a similar survey made over a decade ago. The speckle observations were made with the KPNO and CTIO 4 m telescopes and USNO speckle camera, and they are sensitive to the detection of binaries in the angular separation regime between 0.03" and 5" with relatively bright companions (Delta V < 3). We report on the discovery of companions to 14 OB stars. In total we resolved companions of 41 of 385 O-stars (11%), 4 of 37 Wolf-Rayet stars (11%), and 89 of 139 B-stars (64%; an enriched visual binary sample that we selected for future orbital determinations). We made a statistical analysis of the binary frequency among the subsample that are listed in the Galactic O Star Catalog by compiling published data on other visual companions detected through adaptive optics studies and/or noted in the Washington Double Star Catalog and by collecting published information on radial velocities and spectroscopic binaries. We find that the binary frequency is much higher among O-stars in clusters and associations compared to the numbers for field and runaway O-stars, consistent with predictions for the ejection processes for runaway stars. We present a first orbit for the O-star Delta Orionis, a linear solution of the close, apparently optical, companion of the O-star Iota Orionis, and an improved orbit of the Be star Delta Scorpii. Finally, we list astrometric data for another 249 resolved and 221 unresolved targets that are lower mass stars that we observed for various other science programs.

We present angular diameters of the Hyades giants, gamma, delta^1, epsilon, and theta^1 Tau from interferometric measurements with the CHARA Array. Our errors in the limb-darkened angular diameters for these stars are all less than 2%, and in combination with additional observable quantities, we determine the effective temperatures, linear radii and absolute luminosities for each of these stars. Additionally, stellar masses are inferred from model isochrones to determine the surface gravities. These data show that a new calibration of effective temperatures with errors well under 100K is now possible from interferometric angular diameters of stars.

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.

We present a radial velocity study of the rapidly rotating B-star Regulus that indicates the star is a single-lined spectroscopic binary. The orbital period (40.11 d) and probable semimajor axis (0.35 AU) are large enough that the system is not interacting at present. However, the mass function suggests that the secondary has a low mass (M_2 > 0.30 M_sun), and we argue that the companion may be a white dwarf. Such a star would be the remnant of a former mass donor that was the source of the large spin angular momentum of Regulus itself.

The rapid rotation of Be stars may be caused in some cases by past mass and angular momentum accretion in an interacting binary in which the mass donor is currently viewed as a small, hot subdwarf stripped of its outer envelope. Here we report on the spectroscopic detection of such a subdwarf in the Be binary system FY Canis Majoris from the analysis of data acquired by the IUE spacecraft and KPNO Coude Feed Telescope over the course of 16 and 21 years, respectively. We present a double-lined spectroscopic orbit for the binary based upon radial velocities from the IUE spectra and use the orbital solutions with a Doppler tomography algorithm to reconstruct the components' UV spectra. The subdwarf is hot (T_eff = 45+/-5 kK) and has a mass of about 1.3 M_sun and a radius of about 0.6 R_sun. It contributes about 4% as much flux as the Be star does in the FUV. We also present observations of the H-alpha and He I 6678 emission features that are formed in the circumstellar disk of the Be star. Orbital flux and velocity variations in the He I 6678 profile indicate that much of the emission forms along the disk rim facing the hot subdwarf where the disk is probably heated by the incident radiation from the subdwarf. A study of the FUV infall shell lines discovered in the 1980s confirms their episodic presence but reveals that they tend to be found around both quadrature phases, unlike the pattern in Algol binaries. Phase-dependent variations in the UV N V doublet suggest the presence of a N-enhanced wind from the subdwarf and a possible shock-interaction region between the stars where the subdwarf's wind collides with the disk of the Be star.

We present the results of an adaptive optics survey for faint companions among Galactic O-type star systems (V < 8) using the Advanced Electro-Optical System (AEOS) 3.6-meter telescope on Haleakala. We surveyed these O star systems in I-band, typically being able to detect a companion with a magnitude difference of about 6 in the projected separation range of 0.5 to 1.0 arcseconds, and of about 9.5 in the range of 1.0 to 5.0 arcseconds. In the course of the survey, we discovered 40 new companions among 31 of the 116 objects examined and made astrometric and differential magnitude measurements of 24 additional known pairs, several of them being confirmation detections. We present new astrometric orbits for two binaries, BU 1032AB (WDS 05387-0236; sig Ori AB) and SEE 322 (WDS 17158-3344; HD 155889AB). We lack magnitude differences for other filter bands, so it is difficult to determine physical from line-of-sight companions, but we present empirical arguments for the limiting magnitude difference where field contamination is significant.

We present the results of a spectroscopic investigation of 108 nearby field B-stars. We derive their key stellar parameters, $V \sin i$, $T_{\rm eff}$, $\log g$, and $\log g_{\rm polar}$, using the same methods that we used in our previous cluster B-star survey. By comparing the results of the field and the cluster samples, we find that the main reason for the overall slower rotation of the field sample is that it contains a larger fraction of older stars than found in the (mainly young) cluster sample.

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 results from an optical spectroscopic investigation of the massive binary system [L72] LH~54-425 in the LH 54 OB association in the Large Magellanic Cloud. We revise the ephemeris of [L72] LH 54-425 and find an orbital period of 2.247409 +/- 0.000010 days. We find spectral types of O3 V for the primary and O5 V for the secondary. We made a combined solution of the radial velocities and previously published V-band photometry to determine the inclination for two system configurations, i = 52 degrees for the configuration of the secondary star being more tidally distorted and i = 55 degrees for the primary as the more tidally distorted star. We argue that the latter case is more probable, and this solution yields masses and radii of M_1 = 47 +/- 2 M_Sun and R_1 = 11.4 +/- 0.1 R_Sun for the primary, and M_2 = 28 +/- 1 M_Sun and R_2 = 8.1 +/- 0.1 R_Sun for the secondary. Our analysis places LH 54-425 amongst the most massive stars known. Based on the position of the two stars plotted on a theoretical HR diagram, we find the age of the system to be about 1.5 Myr.

We present results from Hubble Space Telescope UV spectroscopy of the massive X-ray binary system, HD226868 = Cyg X-1. The spectra were obtained at both orbital conjunction phases in two separate runs in 2002 and 2003 when the system was in the X-ray high/soft state. The stellar wind lines suffer large reductions in strength when the black hole is in the foreground due to the X-ray ionization of the wind ions. A comparison of HST and archival IUE spectra shows that similar photoionization effects occur in both the X-ray states. We constructed model UV wind line profiles assuming that X-ray ionization occurs everywhere in the wind except the zone where the supergiant blocks the X-ray flux. The good match between the observed and model profiles indicates that the wind ionization extends to near to the hemisphere of the supergiant facing the X-ray source. The H-alpha emission strength is generally lower in the high/soft state compared to the low/hard state, but the He II 4686 emission is relatively constant between states. The results suggest that mass transfer in Cyg X-1 is dominated by a focused wind flow that peaks along the axis joining the stars and that the stellar wind contribution is shut down by X-ray photoionization effects. The strong stellar wind from the shadowed side of the supergiant will stall when Coriolis deflection brings the gas into the region of X-ray illumination. This stalled gas component may be overtaken by the orbital motion of the black hole and act to inhibit accretion from the focused wind. The variations in the strength of the shadow wind component may then lead to accretion rate changes that ultimately determine the X-ray state.

The microquasar SS 433 is an interacting massive binary consisting of an evolved mass donor and a compact companion that ejects relativistic jets. The mass donor was previously identified through spectroscopic observations of absorption lines in the blue part of the spectrum that showed Doppler shifts associated with orbital motion and strength variations related to the orbital modulation of the star-to-disk flux ratio and to disk obscuration. However, subsequent observations revealed other absorption features that lacked these properties and that were probably formed in the disk gas outflow. We present here follow-up observations of SS 433 at orbital and precession phases identical to those from several previous studies with the goals of confirming the detection of the mass donor spectrum and providing more reliable masses for the two system components. We show that the absorption features present as well as those previously observed almost certainly belong to the mass donor star, and we find revised masses of 12.3 +/- 3.3 and 4.3 +/- 0.8 M(sun) for the mass donor and compact object, respectively.

(Shortened) Ultraviolet observations of the black hole X-ray binary Cygnus X-1 were obtained using the STIS on HSTubble. We detect P Cygni line features show strong, broad absorption components when the X-ray source is behind the companion star and noticeably weaker absorption when the X-ray source is between us and the companion star. We fit the P Cygni profiles using the SEI method applied to a spherically symmetric stellar wind subject to X-ray photoionization from the black hole. The Si IV doublet provides the most reliable estimates of the parameters of the wind and X-ray illumination. The velocity $v$ increases with radius $r$ according to $v=v_\infty(1-r_\star/r)^\beta$, with$\beta\approx0.75$ and $v_\infty\approx1420$ km s$^{-1}$.The microturbulent velocity was $\approx160$ km s$^{-1}$. Our fit implies a ratio of X-ray luminosity to wind mass-loss rate of L$_{X,38}/\dot M_{-6} \approx 0.33$, measured at $\dot M_{-6}$ = 4.8. Our models determine parameters that may be used to estimate the accretion rate onto the black hole and independently predict the X-ray luminosity. Our predicted L$_x$ matches that determined by contemporaneous RXTE ASM remarkably well, but is a factor of 3 lower than the rate according to Bondi-Hoyle-Littleton spherical wind accretion. We suggest that some of the energy of accretion may go into powering a jet.

We present multiple epochs of H-alpha spectroscopy for 47 members of the open cluster NGC 3766 to investigate the long term variability of its Be stars. Sixteen of the stars in this sample are Be stars, including one new discovery. Of these, we observe an unprecedented 11 Be stars that undergo disk appearances and/or near disappearances in our H-alpha spectra, making this the most variable population of Be stars known to date. NGC 3766 is therefore an excellent location to study the formation mechanism of Be star disks. From blue optical spectra of 38 cluster members and existing Stromgren photometry of the cluster, we also measure rotational velocities, effective temperatures, and polar surface gravities to investigate the physical and evolutionary factors that may contribute to the Be phenomenon. Our analysis also provides improvements to the reddening and distance of NGC 3766, and we find E(B-V) = 0.22 +/- 0.03 and (V-M_V)_0 = 11.6 +/- 0.2, respectively. The Be stars are not associated with a particular stage of main-sequence evolution, but they are a population of rapidly rotating stars with a velocity distribution generally consistent with rotation at 70-80% of the critical velocity, although systematic effects probably underestimate the true rotational velocities so that the rotation is much closer to critical. Our measurements of the changing disk sizes are consistent with the idea that transitory, nonradial pulsations contribute to the formation of these highly variable disks.

We present extensive optical spectroscopy of the early-type magnetic star HD 191612 (O6.5f?pe-O8fp). The Balmer and HeI lines show strongly variable emission which is highly reproducible on a well-determined 538-d period. Metal lines and HeII absorptions (including many selective emission lines but excluding He II 4686A emission) are essentially constant in line strength, but are variable in velocity, establishing a double-lined binary orbit with P(orb) = 1542d, e=0.45. We conduct a model-atmosphere analysis of the primary, and find that the system is consistent with a O8: giant with a B1: main-sequence secondary. Since the periodic 538-d changes are unrelated to orbital motion, rotational modulation of a magnetically constrained plasma is strongly favoured as the most likely underlying `clock'. An upper limit on the equatorial rotation is consistent with this hypothesis, but is too weak to provide a strong constraint.

The massive interacting binary RY Scuti is an important representative of an active mass-transferring system that is changing before our eyes and which may be an example of the formation of a Wolf-Rayet star through tidal stripping. Utilizing new and previously published spectra, we present examples of how a number of illustrative absorption and emission features vary during the binary orbit. We identify spectral features associated with each component, calculate a new, double-lined spectroscopic binary orbit, and find masses of 7.1 +/- 1.2 M_sun for the bright supergiant and 30.0 +/- 2.1 M_sun for the hidden massive companion. Through tomographic reconstruction of the component spectra from the composite spectra, we confirm the O9.7 Ibpe spectral class of the bright supergiant and discover a B0.5 I spectrum associated with the hidden massive companion; however, we suggest that the latter is actually the spectrum of the photosphere of the accretion torus immediately surrounding the massive companion. We describe the complex nature of the mass loss flows from the system in the context of recent hydrodynamical models for beta Lyr, leading us to conclude RY Scuti has matter leaving the system in two ways: 1) a bipolar outflow from winds generated by the hidden massive companion, and 2) mass from the bright O9.7 Ibpe supergiant flowing from the region near the L2 point to fill out a large, dense circumbinary disk. This circumbinary disk (radius ~ 1 AU) may feed the surrounding double-toroidal nebula (radius ~ 2000 AU).

We present new results from a radial velocity study of six bright OB stars with little or no prior measurements. One of these, HD 45314, may be a long-period binary, but the velocity variations of this Be star may be related to changes in its circumstellar disk. Significant velocity variations were also found for HD 60848 (possibly related to nonradial pulsations) and HD 61827 (related to wind variations). The other three targets, HD 46150, HD 54879, and HD 206183, are constant velocity objects, but we note that HD 54879 has H$\alpha$ emission that may originate from a binary companion. We illustrate the average red spectrum of each target.

We present the results from an optical spectroscopic analysis of the massive stars HD 37366 and HD 54662. We find that HD 37366 is a double-lined spectroscopic binary with a period of 31.8187 +/- 0.0004 days, and HD 54662 is also a double lined binary with a much longer period of 557.8 +/- 0.3 days. The primary of HD 37366 is classified as O9.5 V, and it contributes approximately two-thirds of the optical flux. The less luminous secondary is a broad-lined, early B-type main-sequence star. Tomographic reconstruction of the individual spectra of HD 37366 reveals absorption lines present in each component, enabling us to constrain the nature of the secondary and physical characteristics of both stars. Tomographic reconstruction was not possible for HD 54662; however, we do present mean spectra from our observations that show that the secondary component is approximately half as bright as the primary. The observed spectral energy distributions (SEDs) were fit with model SEDs and galactic reddening curves to determine the angular sizes of the stars. By assuming radii appropriate for their classifications, we determine distance ranges of 1.4 - 1.9 and 1.2 - 1.5 kpc for HD 37366 and HD 54662, respectively.

We present results from an H-alpha monitoring campaign of the Be X-ray binary systems HDE 245770 = A 0535+26 and X Per. We use the H-alpha equivalent widths together with adopted values of the Be star effective temperature, disk inclination, and disk outer boundary to determine the half-maximum emission radius of the disk as a function of time. The observations of HDE 245770 document the rapid spectral variability that apparently accompanied the regeneration of a new circumstellar disk. This disk grew rapidly during the years 1998 - 2000, but then slowed in growth in subsequent years. The outer disk radius is probably truncated by resonances between the disk gas and neutron star orbital periods. Two recent X-ray outbursts appear to coincide with the largest disk half-maximum emission radius attained over the last decade. Our observations of X Per indicate that its circumstellar disk has recently grown to near record proportions, and concurrently the system has dramatically increased in X-ray flux, presumably the result of enhanced mass accretion from the disk. We find that the H-alpha half-maximum emission radius of the disk surrounding X Per reached a size about six times larger than the stellar radius, a value, however, that is well below the minimum separation between the Be star and neutron star. We suggest that spiral arms excited by tidal interaction at periastron may help lift disk gas out to radii where accretion by the neutron star companion becomes more effective.

We report on a three year spectroscopic monitoring program of the H-alpha emission in the massive X-ray binary LS I +65 010 = 2S 0114+650, which consists of a B-supergiant and a slowly rotating X-ray pulsar. We present revised orbital elements that yield a period of P=11.5983 +/- 0.0006 d and confirm that the orbit has a non-zero eccentricity e=0.18 +/- 0.05. The H-alpha emission profile is formed in the base of the wind of the B-supergiant primary, and we show how this spectral feature varies on timescales that are probably related to the rotational period of the B-supergiant. We also examine the X-ray fluxes from the Rossi X-ray Timing Explorer All-Sky Monitor instrument, and we show that the X-ray orbital light curve has a maximum at periastron and a minimum at the inferior conjunction of the B-supergiant. We also show that the wind emission strength and the high energy X-ray flux appear to vary in tandem on timescales of approximately a year.

We present the results of an H-alpha monitoring campaign on the BeXRB and microquasar system LS I +61 303. We use radial velocity measurements of HeI lines in our spectra to re-evaluate the orbital elements and to better establish the time of periastron. We list equivalent widths and other parameters for the H-alpha emission line and discuss the orbital phase related variations observed. We call attention to a dramatic episode of emission weakening that occurred in less than a day that probably resulted from exposure to a transient source of ionizing radiation. We argue that the increase in H-alpha and X-ray emission following periastron probably results from the creation of an extended density wave in the disk created by tidal forces. We also discuss estimates of the size of the disk from the H-alpha equivalent width measurements, and we suggest that the disk radius from the average equivalent width corresponds to a resonant truncation radius of the disk while the maximum equivalent width corresponds to a radius limited by the separation of the stars at periastron. We note that a nearby faint companion is probably an unrelated foreground object.

Identifying binaries among runaway O- and B-type stars offers valuable insight into the evolution of open clusters and close binary stars. Here we present a spectroscopic investigation of 12 known or suspected binaries among field and runaway OB stars. We find new orbital solutions for five single-lined spectroscopic binaries (HD 1976, HD 14633, HD 15137, HD 37737, and HD 52533), and we classify two stars thought to be binaries (HD 30614 and HD 188001) as single stars. In addition, we reinvestigate their runaway status using our new radial velocity data with the UCAC2 proper motion catalogs. Seven stars in our study appear to have been ejected from their birthplaces, and at least three of these runaways are spectroscopic binaries and are of great interest for future study.

We present numerical models of the circumstellar disks of Be stars, and we describe the resulting synthetic H-alpha emission lines and maps of the wavelength-integrated emission flux projected onto the sky. We demonstrate that there are monotonic relationships between the emission line equivalent width and the ratio of the angular half-width at half maximum of the projected disk major axis to the radius of the star. These relationships depend mainly upon the temperatures of the disk and star, the inclination of the disk normal to the line of sight, and the adopted outer boundary for the disk radius. We show that the predicted H-alpha disk radii are consistent with those observed directly through long baseline interferometry of nearby Be stars (especially once allowance is made for disk truncation in binaries and for dilution of the observed H-alpha equivalent width by continuum disk flux in the V-band).

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.

We derive the effective temperatures and gravities of 461 OB stars in 19 young clusters by fitting the H-gamma profile in their spectra. We use synthetic model profiles for rotating stars to develop a method to estimate the polar gravity for these stars, which we argue is a useful indicator of their evolutionary status. We combine these results with projected rotational velocity measurements obtained in a previous paper on these same open clusters. We find that the more massive B-stars experience a spin down as predicted by the theories for the evolution of rotating stars. Furthermore, we find that the members of binary stars also experience a marked spin down with advanced evolutionary state due to tidal interactions. We also derive non-LTE-corrected helium abundances for most of the sample by fitting the He I 4026, 4387, 4471 lines. A large number of helium peculiar stars are found among cooler stars with Teff < 23000 K. The analysis of the high mass stars (8.5 solar masses < M < 16 solar masses) shows that the helium enrichment process progresses through the main sequence (MS) phase and is greater among the faster rotators. This discovery supports the theoretical claim that rotationally induced internal mixing is the main cause of surface chemical anomalies that appear during the MS phase. The lower mass stars appear to have slower rotation rates among the low gravity objects, and they have a large proportion of helium peculiar stars. We suggest that both properties are due to their youth. The low gravity stars are probably pre-main sequence objects that will spin up as they contract. These young objects very likely host a remnant magnetic field from their natal cloud, and these strong fields sculpt out surface regions with unusual chemical abundances.

Open clusters offer us the means to study stellar properties in samples with well-defined ages and initial chemical composition. Here we present a survey of projected rotational velocities for a large sample of mainly B-type stars in young clusters to study the time evolution of the rotational properties of massive stars. The survey is based upon moderate resolution spectra made with the WIYN 3.5 m and CTIO 4 m telescopes and Hydra multi-object spectrographs, and the target stars are members of 19 young open clusters with an age range of approximately 6 to 73 Myr. We made fits of the observed lines He I 4026, 4387, 4471 and Mg II 4481 using model theoretical profiles to find projected rotational velocities for a total of 496 OB stars. We find that there are fewer slow rotators among the cluster B-type stars relative to nearby B stars in the field. We present evidence consistent with the idea that the more massive B stars (M > 9 solar masses) spin down during their main sequence phase. However, we also find that the rotational velocity distribution appears to show an increase in the numbers of rapid rotators among clusters with ages of 10 Myr and higher. These rapid rotators appear to be distributed between the zero age and terminal age main sequence locations in the Hertzsprung-Russell diagram, and thus only a minority of them can be explained as the result of a spin up at the terminal age main sequence due to core contraction. We suggest instead that some of these rapid rotators may have been spun up through mass transfer in close binary systems.

We present new results from a study of high quality, red spectra of the massive binary star system HD 1383 (B0.5 Ib + B0.5 Ib). We determined radial velocities and revised orbital elements (P = 20.28184 +/- 0.0002 d) and made Doppler tomographic reconstructions of the component spectra. A comparison of these with model spectra from non-LTE, line blanketed atmospheres indicates that both stars have almost identical masses (M_2/M_1 = 1.020 +/- 0.014), temperatures (T_eff = 28000 +/- 1000 K), gravities (log g = 3.25 +/- 0.25), and projected rotational velocities (V sin i < 30 km/s). We investigate a number of constraints on the radii and masses of the stars based upon the absence of eclipses, surface gravity, stellar wind terminal velocity, and probable location in the Perseus spiral arm of the Galaxy, and these indicate a range in probable radius and mass of R/R_sun = 14 - 20 and M/M_sun = 16 - 35, respectively. These values are consistent with model evolutionary masses for single stars of this temperature and gravity. Both stars are much smaller than their respective Roche radii, so the system is probably in a pre-contact stage of evolution. A fit of the system's spectral energy distribution yields a reddening of E(B-V)=0.55 +/- 0.05 and a ratio of total-to-selective extinction of R=2.97 +/- 0.15. We find no evidence of H-alpha emission from colliding stellar winds, which is probably the consequence of the low gas densities in the colliding winds zone.

We have measured the angular diameters of six M dwarfs with the CHARA Array, a long-baseline optical interferometer located at Mount Wilson Observatory. Spectral types range from M1.0 V to M3.0 V and linear radii from 0.38 to 0.69 Rsun. These results are consistent with the seven other M-dwarf radii measurements from optical interferometry and with those for sixteen stars in eclipsing binary systems. We compare all directly measured M dwarf radii to model predictions and find that current models underestimate the true stellar radii by up to 15-20%. The differences are small among the metal-poor stars but become significantly larger with increasing metallicity. This suggests that theoretical models for low mass stars may be missing some opacity source that alters the computed stellar radii.

We present the results of time-resolved spectroscopy of 13 O-type stars in the Cas OB6 stellar association. We conducted a survey for radial velocity variability in search of binary systems, which are expected to be plentiful in young OB associations. Here we report the discovery of two new single-lined binaries, and we present new orbital elements for three double-lined binaries (including one in the multiple star system HD 17505). One of the double-lined systems is the eclipsing binary system DN Cas, and we present a preliminary light curve analysis that yields the system inclination, masses, and radii. We compare the spectra of the single stars and the individual components of the binary stars with model synthetic spectra to estimate the stellar effective temperatures, gravities, and projected rotational velocities. We also make fits of the spectral energy distributions to derive E(B-V), R=A_V/E(B-V), and angular diameter. A distance of 1.9 kpc yields radii that are consistent with evolutionary models. We find that 7 of 14 systems with spectroscopic data are probable binaries, consistent with the high binary frequency found for other massive stars in clusters and associations.

Be stars are a class of rapidly rotating B stars with circumstellar disks that cause Balmer and other line emission. There are three possible reasons for the rapid rotation of Be stars: they may have been born as rapid rotators, spun up by binary mass transfer, or spun up during the main-sequence (MS) evolution of B stars. To test the various formation scenarios, we have conducted a photometric survey of 55 open clusters in the southern sky. Of these, five clusters are probably not physically associated groups and our results for two other clusters are not reliable, but we identify 52 definite Be stars and an additional 129 Be candidates in the remaining clusters. We use our results to examine the age and evolutionary dependence of the Be phenomenon. We find an overall increase in the fraction of Be stars with age until 100 Myr, and Be stars are most common among the brightest, most massive B-type stars above the zero-age MS (ZAMS). We show that a spin-up phase at the terminal-age MS (TAMS) cannot produce the observed distribution of Be stars, but up to 73% of the Be stars detected may have been spun-up by binary mass transfer. Most of the remaining Be stars were likely rapid rotators at birth. Previous studies have suggested that low metallicity and high cluster density may also favor Be star formation. Our results indicate a possible increase in the fraction of Be stars with increasing cluster distance from the Galactic center (in environments of decreasing metallicity). However, the trend is not significant and could be ruled out due to the intrinsic scatter in our data. We also find no relationship between the fraction of Be stars and cluster density.

We present a calculation of the Sun's motion through the Milky Way Galaxy over the last 500 million years. The integration is based upon estimates of the Sun's current position and speed from measurements with Hipparcos and upon a realistic model for the Galactic gravitational potential. We estimate the times of the Sun's past spiral arm crossings for a range in assumed values of the spiral pattern angular speed. We find that for a difference between the mean solar and pattern speed of Omega_Sun - Omega_p = 11.9 +/- 0.7 km/s/kpc the Sun has traversed four spiral arms at times that appear to correspond well with long duration cold periods on Earth. This supports the idea that extended exposure to the higher cosmic ray flux associated with spiral arms can lead to increased cloud cover and long ice age epochs on Earth.

We are conducting a high angular resolution imaging survey of Galactic Wolf-Rayet stars using the Wide Field and Planetary Camera 2 aboard the Hubble Space Telescope. We have found a small stellar cluster associated with the faint, close pair WR 38 and WR 38a. We present astrometric measurements and photometry in the wide-band F336W (U), F439W (B), and F555W (V) filter system for these cluster and nearby stars. We compare their colors and magnitudes with calibrated model results for reddened stars to identify seven probable main sequence members of the cluster. A least-squares fit of the colors and magnitudes of this set yields a cluster reddening of E(B-V) = 1.45 +/- 0.14 mag and a distance of 3.7^{+3.8}_{-1.2} kpc. We discuss the relationship of this cluster to other objects along the line of sight, and we argue that the distance probably lies in the range 5 - 8 kpc (but is not as great as 14.5 kpc distance recently advocated by Shorlin, Turner, & Pedreros). At a distance of 8 kpc, the cluster would reside in a dense region of the Carina spiral arm, close to a giant molecular cloud and the starburst cluster NGC 3603.

We report on K-band interferometric observations of the bright, rapidly rotating star Regulus (type B7 V) made with the CHARA Array on Mount Wilson, California. Through a combination of interferometric and spectroscopic measurements, we have determined for Regulus the equatorial and polar diameters and temperatures, the rotational velocity and period, the inclination and position angle of the spin axis, and the gravity darkening coefficient. These first results from the CHARA Array provide the first interferometric measurement of gravity darkening in a rapidly rotating star and represent the first detection of gravity darkening in a star that is not a member of an eclipsing binary system.

We describe a technique to identify Be stars in open clusters using Stromgren b, y, and narrow-band Halpha photometry. We first identify the B-type stars of the cluster using a theoretical isochrone fit to the (b-y, y) color-magnitude diagram. The strongest Be stars are easily identified in a (b-y, y-Halpha) color-color diagram, but those with weaker Halpha emission (classified as possible Be star detections) may be confused with evolved or foreground stars. Here we present such photometry plus Halpha spectroscopy of members of the cluster NGC 3766 to demonstrate the accuracy of our technique. Statistical results on the relative numbers of Be and B-type stars in additional clusters will be presented in a future paper.

We present results from a radial velocity study of two runaway O-type stars, HD 14633 (ON8.5 V) and HD 15137 (O9.5 III(n)). We find that HD 14633 is a single-lined spectroscopic binary with an orbital period of 15.4083 days. The second target HD 15137 is a radial velocity variable and a possible single-lined spectroscopic binary with a period close to 1 month. Both binaries have large eccentricity, small semiamplitude, and a small mass function. We show the trajectories of the stars in the sky based upon an integration of motion in the Galactic potential, and we suggest that both stars were ejected from the vicinity of the open cluster NGC 654 in the Perseus spiral arm. The binary orbital parameters and runaway velocities are consistent with the idea that both these stars were ejected by supernova explosions in binaries and that they host neutron star companions. We find that the time-of-flight since ejection is longer than the predicted evolutionary timescales for the stars, which may indicate that they have a lower mass than normally associated with their spectral classifications and/or that their lives have been extended through rapid rotation.

Recent theoretical predictions for low metallicity massive stars predict that these stars should have drastically reduced equatorial winds (mass loss) while on the main sequence, and as such should retain most of their angular momentum. Observations of both the Be/(B+Be) ratio and the blue-to-red supergiant ratio appear to have a metallicity dependence that may be caused by high rotational velocities. We have analyzed 39 archival Hubble Space Telescope Imaging Spectrograph (STIS), high resolution, ultraviolet spectra of O-type stars in the Magellanic Clouds to determine their projected rotational velocities V sin i. Our methodology is based on a previous study of the projected rotational velocities of Galactic O-type stars using International Ultraviolet Explorer (IUE) Short Wavelength Prime (SWP) Camera high dispersion spectra, which resulted in a catalog of V sin i values for 177 O stars. Here we present complementary V sin i values for 21 Large Magellanic Cloud and 22 Small Magellanic Cloud O-type stars based on STIS and IUE UV spectroscopy. The distribution of V sin i values for O type stars in the Magellanic Clouds is compared to that of Galactic O type stars. Despite the theoretical predictions and indirect observational evidence for high rotation, the O type stars in the Magellanic Clouds do not appear to rotate faster than their Galactic counterparts.

We present an investigation of H-alpha emission line variations observed in the massive Algol binary, RY Per. We give new radial velocity data for the secondary based upon our optical spectra and for the primary based upon high dispersion UV spectra. We present revised orbital elements and an estimate of the primary's projected rotational velocity (which indicates that the primary is rotating 7 times faster than synchronous). We use a Doppler tomography algorithm to reconstruct the individual primary and secondary spectra in the region of H-alpha, and we subtract the latter from each of our observations to obtain profiles of the primary and its disk alone. Our H-alpha observations of RY Per show that the mass gaining primary is surrounded by a persistent but time variable accretion disk. The profile that is observed outside-of-eclipse has weak, double-peaked emission flanking a deep central absorption, and we find that these properties can be reproduced by a disk model that includes the absorption of photospheric light by the band of the disk seen in projection against the face of the star. We developed a new method to reconstruct the disk surface density distribution from the ensemble of H-alpha profiles observed around the orbit, and this method accounts for the effects of disk occultation by the stellar components, the obscuration of the primary by the disk, and flux contributions from optically thick disk elements. The resulting surface density distribution is elongated along the axis joining the stars, in the same way as seen in hydrodynamical simulations of gas flows that strike the mass gainer near trailing edge of the star. This type of gas stream configuration is optimal for the transfer of angular momentum, and we show that rapid rotation is found in other Algols that have passed through a similar stage.

We present spectroscopy of the microquasar SS 433 obtained near primary eclipse and disk precessional phase Psi = 0.0, when the accretion disk is expected to be most ``face-on''. The likelihood of observing the spectrum of the mass donor is maximized at this combination of orbital and precessional phases since the donor is in the foreground and above the extended disk believed to be present in the system. The spectra were obtained over four different runs centered on these special phases. The blue spectra show clear evidence of absorption features consistent with a classification of A3-7 I. The behavior of the observed lines indicates an origin in the mass donor. The observed radial velocity variations are in anti-phase to the disk, the absorption lines strengthen at mid-eclipse when the donor star is expected to contribute its maximum percentage of the total flux, and the line widths are consistent with lines created in an A supergiant photosphere. We discuss and cast doubt on the possibility that these lines represent a shell spectrum rather than the mass donor itself. We re-evaluate the mass ratio of the system and derive masses of 10.9 +/- 3.1 Msun and 2.9 +/- 0.7 Msun for the mass donor and compact object plus disk, respectively. We suggest that the compact object is a low mass black hole. In addition, we review the behavior of the observed emission lines from both the disk/wind and high velocity jets.

We present an investigation of new optical and ultraviolet spectra of the mass donor star in the massive X-ray binary LS 5039. The optical band spectral line strengths indicate that the atmosphere is N-rich and C-poor, and we classify the stellar spectrum as type ON6.5 V((f)). The N-strong and C-weak pattern is also found in the stellar wind P Cygni lines of N V 1240 and C IV 1550. We suggest that the N-enrichment may result from internal mixing if the O-star was born as a rapid rotator, or the O-star may have accreted N-rich gas prior to a common-envelope interaction with the progenitor of the supernova. We re-evaluated the orbital elements to find an orbital period of P=4.4267 +/- 0.0010 d. We compared the spectral line profiles with new non-LTE, line-blanketed model spectra, from which we derive an effective temperature T_eff = 37.5 +/- 1.7 kK, gravity log g = 4.0 +/- 0.1, and projected rotational velocity V sin i = 140 +/- 8 km/s. We fit the UV, optical, and IR flux distribution using a model spectrum and extinction law with parameters E(B-V)= 1.28 +/- 0.02 and R= 3.18 +/- 0.07. We confirm the co-variability of the observed X-ray flux and stellar wind mass loss rate derived from the H-alpha profile, which supports the wind accretion scenario for the X-ray production in LS 5039. Wind accretion models indicate that the compact companion has a mass M_X/M_sun = 1.4 +/- 0.4, consistent with its identification as a neutron star. The observed eccentricity and runaway velocity of the binary can only be reconciled if the neutron star received a modest kick velocity due to a slight asymmetry in the supernova explosion (during which >5 solar masses was ejected).

We present the results of a spectroscopic monitoring program (from 1998 to 2002) of the H-alpha emission strength in HDE 226868, the optical counterpart of the black hole binary, Cyg X-1. The H-alpha profiles consist of (1) a P Cygni component associated with the wind of the supergiant, (2) emission components that attain high velocity at the conjunctions and that probably form in enhanced outflows both towards and away from the black hole, and (3) an emission component that moves in anti-phase with the supergiant's motion. We argue that the third component forms in accreted gas near the black hole, and the radial velocity curve of the emission is consistent with a mass ratio of M_X / M_opt = 0.36 +/- 0.05. We find that there is a general anti-correlation between the H-alpha emission strength and X-ray flux in the sense that when the H-alpha emission is strong (W_\lambda < -0.5 Angstroms) the X-ray flux is weaker and the spectrum harder. On the other hand, there is no correlation between H-alpha emission strength and X-ray flux when H-alpha is weak. During the low/hard X-ray state, the strong wind is fast and the accretion rate is relatively low, while in the high/soft state the weaker, highly ionized wind attains only a moderate velocity and the accretion rate increases. We argue that the X-ray transitions from the normal low/hard to the rare high/soft state are triggered by episodes of decreased mass loss rate in the supergiant donor star.

We present results from a short series of blue, moderate resolution spectra of the microquasar binary, SS 433. The observations were made at a time optimized to find the spectrum of the donor star, i.e., when the donor was in the foreground and well above the plane of the obscuring disk. In addition to the well-known stationary and jet emission lines, we find evidence of a weak absorption spectrum that resembles that of an A-type evolved star. These lines display radial velocity shifts opposite to those associated with the disk surrounding the compact star, and they appear strongest when the disk is maximally eclipsed. All these properties suggest that these absorption lines form in the atmosphere of the hitherto unseen mass donor star in SS 433. The radial velocity shifts observed are consistent with a mass ratio M_X / M_O = 0.57 +/- 0.11 and masses of M_O = (19 +/- 7) M_sun and M_X = (11 +/- 5) M_sun. These results indicate that the system consists of an evolved, massive donor and a black hole mass gainer.

We present the first orbital elements for the massive close binary, HD 101131, one of the brightest objects in the young, open cluster, IC 2944. This system is a double-lined spectroscopic binary in an elliptical orbit with a period of 9.64659 +/- 0.00012 days. It is a young system of unevolved stars (approximately 2 million years old) that are well within their critical Roche surfaces. We use a Doppler tomography algorithm to reconstruct the individual component optical spectra, and we apply well known criteria to arrive at classifications of O6.5 V((f)) and O8.5 V for the primary and secondary, respectively. We compare the reconstructed spectra of the components to single star, spectrum standards to determine a flux ratio of f_2/f_1 = 0.55 +/- 0.08 in the V-band. Both components are rotating faster than synchronously. We estimate the temperatures and luminosities of the components from the observed spectral classifications, composite V magnitude, and cluster distance modulus. The lower limits on the masses derived from the orbital elements and the lack of eclipses are 25 solar masses and 14 solar masses for the primary and secondary, respectively. These limits are consistent with the larger masses estimated from the positions of the stars in the Hertzsprung-Russell diagram and evolutionary tracks for single stars.

We present the first orbital elements for the massive close binary, HD 115071, a double-lined spectroscopic binary in a circular orbit with a period of 2.73135 +/- 0.00003 days. The orbital semiamplitudes indicate a mass ratio of M_2/M_1 = 0.58 +/- 0.02 and yet the stars have similar luminosities. We used a Doppler tomography algorithm to reconstruct the individual component optical spectra, and we applied well known criteria to arrive at classifications of O9.5 V and B0.2 III for the primary and secondary, respectively. We present models of the Hipparcos light curve of the ellipsoidal variations caused by the tidal distortion of the secondary, and the best fit model for a Roche-filling secondary occurs for an inclination of i = 48.7 +/- 2.1 degrees. The resulting masses are 11.6 +/- 1.1 and 6.7 +/- 0.7 solar masses for the primary and secondary, respectively, so that both stars are very overluminous for their mass. The system is one of only a few known semi-detached, Algol-type binaries that contain O-stars. We suggest that the binary has recently emerged from extensive mass transfer (possibly through a delayed contact and common envelope process).

There is much evidence to suggest that stellar wind capture, rather than Roche lobe overflow, serves as the accretion mechanism onto the compact secondary object in the massive X-ray binary LS 5039. The lack of significant emission combined with only a modest X-ray flux provide observational evidence that no large-scale mass transfer is occurring (consistent with our estimate of the radius of the O6.5 V((f)) optical star that is smaller than its critical Roche radius). Here we determine the mass loss rate of the optical star from the broad, residual emission in the H-alpha profile. Using a stellar wind accretion model for a range in assumed primary mass, we compute the predicted X-ray luminosity for the system. We compare our results to the observed X-ray luminosity to determine the mass of the compact object for each case. The companion appears to be a neutron star with a mass between 1 and 3 solar masses. With our new constraints on the masses of both components, we discuss their implications on the evolution of the system before and after the supernova event that created the compact companion. The binary experienced significant mass loss during the supernova, and we find that the predictions for the resulting runaway velocity agree well with the observed peculiar space velocity. LS 5039 may be the fastest runaway object among known massive X-ray binaries.

Double-lined spectroscopic orbital elements have recently been found for the central binary in the massive triple, delta Orionis A based on radial velocities from cross-correlation techniques applied to IUE high dispersion spectra and He I 6678 spectra obtained at Kitt Peak. The primary and secondary velocity amplitudes were found to be 94.9 +/- 0.6 km/s and 186 +/- 9 km/s respectively. Tomographic reconstructions of the primary and secondary stars' spectra confirm the O9.5 II classification of the primary and indicate a B0.5 III type for the secondary. The widths of the UV cross-correlation functions are used to estimate the projected rotational velocities, Vsin i = 157 +/- 6 km/s and 138 +/- 16 km/s for the primary and secondary, respectively implying that both stars rotate faster than their orbital motion. We used the spectroscopic results to make a constrained fit of the Hipparcos light curve of this eclipsing binary, and the model fits limit the inclination to the range between 67 and 77 degrees. The i = 67 degrees solution, which corresponds to a near Roche-filling configuration, results in a primary mass of 11.2 solar masses and a secondary mass of 5.6 solar masses, both of which are substantially below the expected masses for stars of their luminosity. This binary may have experienced a mass ratio reversal caused by Case A Roche lobe overflow, or the system may have suffered extensive mass loss through a binary interaction, perhaps during a common envelope phase, in which most of the primary's mass was lost from the system rather than transferred to the secondary.

We present results from new optical and UV spectroscopy of the unusual binary system SS 433, and we discuss the relationship of the particular spectral components we observe to the properties of the binary. (1) The continuum spectrum which we associate with flux from the super-Eddington accretion disk and the dense part of its wind. (2) H-alpha moving components which are formed far from the binary orbital plane in the relativistic jets. (3) H-alpha and He I "stationary" emission lines which we suggest are formed in the disk wind in a volume larger than the dimensions of the binary. (4) A weak "stationary" emission feature we identify as a C II 7231,7236 blend that attains maximum radial velocity at the orbital quadrature of disk recession. (5) Absorption and emission features from outflowing clumps in the disk wind (seen most clearly in an episode of blue-shifted Na I emission). (6) We found no clear evidence of the absorption line spectrum of the optical star, although we point out the presence of He I absorption features (blended with the stationary emission) with the expected radial velocity trend at the orbital and precessional phases when the star might best be seen. (7) A rich interstellar absorption spectrum of diffuse interstellar bands. The results suggest that the binary is embedded in an expanding thick disk (detected in recent radio observations) which is fed by the wind from the super-Eddington accretion disk.