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.

The rapidly rotating Be star phi Persei was spun up by mass and angular momentum transfer from a now stripped-down, hot subdwarf companion. Here we present the first high angular resolution images of phi Persei made possible by new capabilities in longbaseline interferometry at near-IR and visible wavelengths. We observed phi Persei with the MIRC and VEGA instruments of the CHARA Array. Additional MIRC-only observations were performed to track the orbital motion of the companion, and these were fit together with new and existing radial velocity measurements of both stars to derive the complete orbital elements and distance. The hot subdwarf companion is clearly detected in the near-IR data at each epoch of observation with a flux contribution of 1.5% in the H band, and restricted fits indicate that its flux contribution rises to 3.3% in the visible. A new binary orbital solution is determined by combining the astrometric and radial velocity measurements. The derived stellar masses are 9.6+-0.3Msol and 1.2+-0.2Msol for the Be primary and subdwarf secondary, respectively. The inferred distance (186 +- 3 pc), kinematical properties, and evolutionary state are consistent with membership of phi Persei in the alpha Per cluster. From the cluster age we deduce significant constraints on the initial masses and evolutionary mass transfer processes that transformed the phi Persei binary system. The interferometric data place strong constraints on the Be disk elongation, orientation, and kinematics, and the disk angular momentum vector is coaligned with and has the same sense of rotation as the orbital angular momentum vector. The VEGA visible continuum data indicate an elongated shape for the Be star itself, due to the combined effects of rapid rotation, partial obscuration of the photosphere by the circumstellar disk, and flux from the bright inner disk.

We report 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 the discovery of \geq 100 MeV {\gamma} rays from the binary system PSR B1259-63/LS 2883 using the Large Area Telescope (LAT) on board Fermi. The system comprises a radio pulsar in orbit around a Be star. We report on LAT observations from near apastron to ~ 60 days after the time of periastron, tp, on 2010 December 15. No {\gamma}-ray emission was detected from this source when it was far from periastron. Faint {\gamma}-ray emission appeared as the pulsar approached periastron. At ~ tp + 30d, the \geq 100 MeV {\gamma}-ray flux increased over a period of a few days to a peak flux 20-30 times that seen during the pre-periastron period, but with a softer spectrum. For the following month, it was seen to be variable on daily time scales, but remained at ~ 1 - 4 \times 10^-6 cm^-2 s^-1 before starting to fade at ~ tp + 57d. The total {\gamma}-ray luminosity observed during this period is comparable to the spin-down power of the pulsar. Simultaneous radio and X-ray observations of the source showed no corresponding dramatic changes in radio and X-ray flux between the pre-periastron and post-periastron flares. We discuss possible explanations for the observed {\gamma}-ray-only flaring of the source.

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.

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.

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.

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.

LS I +61 303 and LS 5039 are two of only a handful of known high mass X-ray binaries (HMXBs) that exhibit very high energy emission in the MeV-TeV range, and these "gamma-ray binaries" are of renewed interest due to the recent launch of the Fermi Gamma-ray Space Telescope. Here we present new radial velocities of both systems based on recent red and blue optical spectra. Both systems have somewhat discrepant orbital solutions available in the literature, and our new measurements result in improved orbital elements and resolve the disagreements. The improved geometry of each orbit will aid in studies of the high energy emission region near each source.

Multi-wavelength studies at radio, infrared, optical, X-ray, and TeV wavelengths have discovered probable counterparts to many Galactic sources of GeV emission detected by EGRET. These include pulsar wind nebulae, high mass X-ray binaries, and mixed morphology supernova remnants. Here we provide an overview of the observational properties of Galactic sources which emit across 19 orders of magnitude in energy. We also present new observations of several sources.

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

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.

The production of runaway massive binaries offers key insights into the evolution of close binary stars and open clusters. The stars HD 14633 and HD 15137 are rare examples of such runaway systems, and in this work we investigate the mechanism by which they were ejected from their parent open cluster, NGC 654. We discuss observational characteristics that can be used to distinguish supernova ejected systems from those ejected by dynamical interactions, and we present the results of a new radio pulsar search of these systems as well as estimates of their predicted X-ray flux assuming that each binary contains a compact object. Since neither pulsars nor X-ray emission are observed in these systems, we cannot conclude that these binaries contain compact companions. We also consider whether they may have been ejected by dynamical interactions in the dense environment where they formed, and our simulations of four-body interactions suggest that a dynamical origin is possible but unlikely. We recommend further X-ray observations that will conclusively identify whether HD 14633 or HD 15137 contain neutron stars.

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