We present a sample of 4,132 Mira-like variables (red variables with long periods and high amplitudes) in the Kilodegree Extremely Little Telescope (KELT) survey. Of these, 814 are new detections. We used 2MASS colors to identify candidate asymptotic giant branch (AGB) stars. We tested for photometric variability among the sample and used Lomb-Scargle to determine the periodicity of the variable sample. We selected variables with high amplitudes and strong periodic behavior using a Random Forest classifier. Of the sample of 4,132 Mira-like variables, we estimate that 70% are Miras, and 30% are semi-regular (SR) variables. We also adopt the method of using (W_{RP} - W_{K_s}) vs. (J - K_s) colors (Lebzelter et al. 2018) in distinguishing between O-rich and C-rich Miras and find it to be an improvement over 2MASS colors.

Time domain and multi-messenger astrophysics are growing and important modes of observational astronomy that will help define astrophysics in the 2020s. Significant effort is being put into developing the components of a follow-up system for dynamically turning survey alerts into data. This system consists of: 1) brokers that will aggregate, classify, and filter alerts; 2) Target Observation Managers (TOMs) for prioritizing targets and managing observations and data; and 3) observatory interfaces, schedulers, and facilities along with data reduction software and science archives. These efforts need continued community support and funding in order to complete and maintain them. Many of the efforts can be community open-source software projects but they will benefit from the leadership of professional software developers. The coordination should be done by institutions that are involved in the follow-up system such as the national observatories (e.g. LSST/Gemini/NOAO Mid-scale/Community Science and Data Center) or a new MMA institute. These tools will help the community to produce the most science from new facilities and will provide new capabilities for all users of the facilities that adopt them.

The field of asteroseismology has enjoyed a large swath of data coming from recent missions (e.g., CoRoT, Kepler, K2). This wealth of new data has allowed the field to expand beyond the previous limitation of a few extremely bright and evolved stars. Asteroseismology relies on accurate surface measurements for boundary conditions, but the predicted physical parameters in the Kepler Input Catalog (KIC) are unreliable for hot stars. We present stellar parameters of 25 candidate pulsating B star candidates in the Kepler field. We use blue optical spectra to measure the projected rotational velocity (V sini), effective temperature (Teff), and surface gravity (logg) using TLUSTY and Kurucz ATLAS9 model atmospheres. We find a large discrepancy between our spectroscopically derived parameters and those derived from photometry in the KIC and Gaia Data Release 2 (DR2). Using spectral energy distributions, we also measure the radii of these stars and later calculate the luminosities and masses. We find the extinctions (AV ) of these stars to be consistent with zero, which is expected for stars of high Galactic latitude.

The Carina Nebula is an active star forming region in the southern sky that is of particular interest due to the presence of a large number of massive stars in a wide array of evolutionary stages. Here we present the results of the spectroscopic analysis of 82 B-type stars and 33 O-type stars that were observed in 2013 and 2014. For 82 B-type stars without line blending, we fit model spectra from the Tlusty BSTAR2006 grid to the observed profiles of H{\gamma} and He {\lambda}{\lambda} 4026, 4388, and 4471 to measure the effective temperatures, surface gravities, and projected rotational velocities. We also measure the masses, ages, radii, bolometric luminosities, and distances of these stars. From the radial velocities measured in our sample, we find 31 single lined spectroscopic binary candidates. We find a high dispersion of radial velocities among our sample stars, and we argue that the Carina Nebula stellar population has not yet relaxed and become virialized.

In order to study the growth and evolution of circumstellar disks around classical Be stars, we analyze optical time-series photometry from the KELT survey with simultaneous infrared and visible spectroscopy from the APOGEE survey and BeSS database for a sample of 160 Galactic classical Be stars. The systems studied here show variability including transitions from a diskless to a disk-possessing state (and vice versa), and persistent disks that vary in strength, being replenished at either regularly or irregularly occurring intervals. We detect disk-building events (outbursts) in the light curves of 28\% of our sample. Outbursts are more commonly observed in early- (57\%), compared to mid- (27\%) and late-type (8\%) systems. A given system may show anywhere between 0 -- 40 individual outbursts in its light curve, with amplitudes ranging up to $\sim$0.5 mag and event durations between $\sim$2 -- 1000 days. We study how both the photometry and spectroscopy change together during active episodes of disk growth or dissipation, revealing details about the evolution of the circumstellar environment. We demonstrate that photometric activity is linked to changes in the inner disk, and show that, at least in some cases, the disk growth process is asymmetrical. Observational evidence of Be star disks both growing and clearing from the inside out is presented. The duration of disk buildup and dissipation phases are measured for 70 outbursts, and we find that the average outburst takes about twice as long to dissipate as it does to build up in optical photometry. Our analysis hints that dissipation of the inner disk occurs relatively slowly for late-type Be stars.

Be stars have generally been characterized by the emission lines in their spectra, and especially the time variability of those spectroscopic features. They are known to also exhibit photometric variability at multiple timescales, but have not been broadly compared and analyzed by that behavior. We have taken advantage of the advent of wide-field, long-baseline, and high-cadence photometric surveys that search for transiting exoplanets to perform a comprehensive analysis of brightness variations among a large number of known Be stars. The photometric data comes from the KELT transit survey, with a typical cadence of 30 minutes, baseline of up to ten years, photometric precision of about 1%, and coverage of about 60% of the sky. We analyze KELT light curves of 610 known Be stars in both the Northern and Southern hemispheres in an effort to study their variability. Consistent with other studies of Be star variability, we find most of the stars to be photometrically variable. We derive lower limits on the fraction of stars in our sample that exhibit features consistent with non-radial pulsations (25%), outbursts (36%), and long term trends in the circumstellar disk (37%), and show how these are correlated with spectral sub-type. Other types of variability, such as those owing to binarity, are also explored. Simultaneous spectroscopy for some of these systems from the Be Star Spectral Database (BeSS) allow us to better understand the physical causes for the observed variability, especially in cases of outbursts and changes in the disk.

The Carina star-forming region is one of the largest in the Galaxy, and its massive star population is still being unveiled. The large number of stars combined with high, and highly variable, interstellar extinction makes it inherently difficult to find OB stars in this type of young region. We present the results of a spectroscopic campaign to study the massive star population of the Carina Nebula, with the primary goal to confirm or reject previously identified Carina OB star candidates. A total of 141 known O- and B-type stars and 94 candidates were observed, of which 73 candidates had a high enough signal-to-noise ratio to classify. We find 23 new OB stars within the Carina Nebula, a 32% confirmation rate. One of the new OB stars has blended spectra and is suspected to be a double-lined spectroscopic binary (SB2). We also reclassify the spectral types of the known OB stars and discover nine new SB2s among this population. Finally, we discuss the spatial distribution of these new OB stars relative to known structures in the Carina NSebula.

We report on the results of the extensive multi-wavelength campaign from optical to GeV gamma-rays of the 2014 periastron passage of PSR B1259-63, which is a unique high-mass gamma-ray emitting binary system with a young pulsar companion. Observations demonstrate the stable nature of the post-periastron GeV flare and prove the coincidence of the flare with the start of rapid decay of the H$\alpha$ equivalent width, usually interpreted as a disruption of the Be stellar disk. Intensive X-ray observations reveal changes in the X-ray spectral behaviour happening at the moment of the GeV flare. We demonstrate that these changes can be naturally explained as a result of synchrotron cooling of monoenergetic relativistic electrons injected into the system during the GeV flare.

We report an updated analysis of the gamma-ray source AGL J2241+4454 that was detected as a brief two-day flare in 2010 by the AGILE satellite. The high-energy emission of AGL J2241+4454 has been attributed to the binary system HD 215227, which consists of a Be star being orbited by a black hole making it the first known Be-black hole binary system. We have analyzed the AGILE data and find a gamma-ray flux of $(1.8\pm0.7)\times10^{-6}$ ph cm$^{-2}$ s$^{-1}$, in agreement with the initial report. Additionally, we examined data from the Fermi LAT over several time intervals including the two day flare, the folded orbital phase, and the entire mission ($\sim$6-years). We do not detect AGL J2241+4454 over any of these time periods with Fermi and find upper limits of $1.1\times10^{-7}$ ph cm$^{-2}$ s$^{-1}$ and $5.2\times10^{-10}$ ph cm$^{-2}$ s$^{-1}$ for the flare and the full mission, respectively. We conclude that the HD 215227 Be-black hole binary is not a true gamma-ray binary as previous speculated. While analyzing the Fermi data of the AGL J2241+4454 region, we discovered a previously unknown gamma-ray source with average flux of $(13.56\pm0.02)\times10^{-8}$ ph cm$^{-2}$ s$^{-1}$ that is highly variable on monthly timescales. We associate this emission with the known quasar 87GB 215950.2+503417.

We report on broad multi-wavelength observations of the 2010-2011 periastron passage of the gamma-ray loud binary system PSR B1259-63. High resolution interferometric radio observations establish extended radio emission trailing the position of the pulsar. Observations with the Fermi Gamma-ray Space Telescope reveal GeV gamma-ray flaring activity of the system, reaching the spin-down luminosity of the pulsar, around 30 days after periastron. There are no clear signatures of variability at radio, X-ray and TeV energies at the time of the GeV flare. Variability around periastron in the H$\alpha$ emission line, can be interpreted as the gravitational interaction between the pulsar and the circumstellar disk. The equivalent width of the H$\alpha$ grows from a few days before periastron until a few days later, and decreases again between 18 and 46 days after periastron. In near infrared we observe the similar decrease of the equivalent width of Br$\gamma$ line between the 40th and 117th day after the periastron. For the idealized disk, the variability of the H$\alpha$ line represents the variability of the mass and size of the disk. We discuss possible physical relations between the state of the disk and GeV emission under assumption that GeV flare is directly related to the decrease of the disk size.

A rare group of high mass X-ray binaries (HMXBs) are known that also exhibit MeV, GeV, and/or TeV emission ("gamma-ray binaries"). Expanding the sample of gamma-ray binaries and identifying unknown Fermi sources are currently of great interest to the community. Based upon their positional coincidence with the unidentified Fermi sources 1FGL J1127.7-6244c and 1FGL J1808.5-1954c, the Be stars HD 99771 and HD 165783 have been proposed as gamma-ray binary candidates. During Fermi Cycle 4, we have performed multiwavelength observations of these sources using XMM-Newton and the CTIO 1.5m telescope. We do not confirm high energy emission from the Be stars. Here we examine other X-ray sources in the field of view that are potential counterparts to the Fermi sources.

We present a study of the B and Be star populations of the Double Cluster h and chi Persei. Blue optical spectroscopy is used to first measure projected rotational velocity, V sin i, effective surface temperature, T_eff, and surface gravity, log g, for B-type sample stars, while available Stromgren photometry is used to calculate T_eff and log g for the Be stars showing emission. In our sample of 104 objects for which we measured these stellar parameters, 28 are known or proposed Be stars. Of these Be stars, 22 show evidence of emission at the times of our observations, and furthermore, we find evidence in our data and the literature for at least 8 transient Be stars in the clusters. We find that the Be stars are not rotating near their critical velocity, contrary to the results of studies of similar open clusters. We compare the results of our analysis with other previous studies, and find that the cluster members are more evolved than found by Huang & Gies but still retain much of their initial rotational angular momentum.

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.

Gamma-ray binaries emit most of their radiated power beyond ~10 MeV. The non-thermal emission is thought to arise from the interaction of the relativistic wind of a rotation-powered pulsar with the stellar wind of its massive (O or Be) companion star. A powerful pulsar creates an extended cavity, filled with relativistic electrons, in the radiatively-driven wind of the massive star. As a result, the observed P Cyg profiles of UV resonant lines from the stellar wind should be different from those of single massive stars. We propose to use UV emission lines to detect and constrain the colliding wind region in gamma-ray binaries. We compute the expected orbital variability of P Cyg profiles depending upon the interaction geometry (set by the ratio of momentum fluxes from the winds) and the line-of-sight to the system. We predict little or no variability for the case of LS 5039 and PSR B1259-63, in agreement with currently available HST observations of LS 5039. However, variability between superior and inferior conjunction is expected in the case of LS I+61 303.

The gamma-ray binary LS I +61 303 is one of the brightest Fermi sources, with orbitally modulated emission across the electromagnetic spectrum. Here we present H-alpha spectra of LS I +61 303 that exhibit a dramatic emission burst shortly before apastron, observed as a redshifted shoulder in the line profile. A correlated burst in radio, X-ray, and GeV emission is observed at the same orbital phase. We interpret the source of the emission as a compact pulsar wind nebula that forms when a tidal mass stream from the Be circumstellar disk interacts with the relativistic pulsar wind. The H-alpha emission offers an important probe of the high energy emission morphology in this system.

The recently discovered gamma-ray binary 1FGL J1018.6-5856 has a proposed optical/near-infrared (OIR) counterpart 2MASS 10185560-5856459. We present Stromgren photometry of this star to investigate its photometric variability and measure the reddening and distance to the system. We find that the gamma-ray binary has E(B-V) = 1.34 +/- 0.04 and d = 5.4^+4.6_-2.1 kpc. While E(B-V) is consistent with X-ray observations of the neutral hydrogen column density, the distance is somewhat closer than some previous authors have suggested.

LS I +61 303 and LS 5039 are exceptionally rare examples of HMXBs with MeV-TeV emission, making them two of only five known or proposed "gamma-ray binaries". There has been disagreement within the literature over whether these systems are microquasars, with stellar winds accreting onto a compact object to produce high energy emission and relativistic jets, or whether their emission properties might be better explained by a relativistic pulsar wind colliding with the stellar wind. Here we present an attempt to detect radio pulsars in both systems with the Green Bank Telescope. The upper limits of flux density are between 4.1-14.5 uJy, and we discuss the null results of the search. Our spherically symmetric model of the wind of LS 5039 demonstrates that any pulsar emission will be strongly absorbed by the dense wind unless there is an evacuated region formed by a relativistic colliding wind shock. LS I +61 303 contains a rapidly rotating Be star whose wind is concentrated near the stellar equator. As long as the pulsar is not eclipsed by the circumstellar disk or viewed through the densest wind regions, detecting pulsed emission may be possible during part of the orbit.

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.

HD 259440 is a B0pe star that was proposed as the optical counterpart to the gamma-ray source HESS J0632+057. Here we present optical spectra of HD 259440 acquired to investigate the stellar parameters, the properties of the Be star disk, and evidence of binarity in this system. Emission from the H-alpha line shows evidence of a spiral density wave in the nearly edge-on disk. We find a best fit stellar effective temperature of 27500-30000 K and a log surface gravity of 3.75-4.0, although our fits are somewhat ambiguous due to scattered light from the circumstellar disk. We derive a mass of 13.2-19.0 M_sun and a radius of 6.0-9.6 R_sun. By fitting the spectral energy distribution, we find a distance between 1.1-1.7 kpc. We do not detect any significant radial velocity shifts in our data, ruling out orbital periods shorter than one month. If HD 259440 is a binary, it is likely a long period (> 100 d) system.

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

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.

We present blue optical spectra of 92 members of h and chi Per obtained with the WIYN telescope at Kitt Peak National Observatory. From these spectra, several stellar parameters were measured for the B-type stars, including V sin i, T_eff, log g_polar, M_star, and R_star. Stromgren photometry was used to measure T_eff and log g_polar for the Be stars. We also analyze photometric data of cluster members and discuss the near-to-mid IR excesses of Be stars.

Non-radial pulsations (NRPs) are a proposed mechanism for the formation of decretion disks around Be stars and are important tools to study the internal structure of stars. NGC 3766 has an unusually large fraction of transient Be stars, so it is an excellent location to study the formation mechanism of Be star disks. High resolution spectroscopy can reveal line profile variations from NRPs, allowing measurements of both the degree, l, and azimuthal order, m. However, spectroscopic studies require large amounts of time with large telescopes to achieve the necessary high S/N and time domain coverage. On the other hand, multi-color photometry can be performed more easily with small telescopes to measure l only. Here, we present representative light curves of Be stars and non-emitting B stars in NGC 3766 from the CTIO 0.9m telescope in an effort to study NRPs in this cluster.

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.

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.

The NSF's Astronomy and Astrophysics Postdoctoral Fellowship (AAPF) is exceptional among the available postdoctoral awards in Astronomy and Astrophysics. The fellowship is one of the few that allows postdoctoral researchers to pursue an original research program, of their own design, at the U.S. institution of their choice. However, what makes this fellowship truly unique is the ability of Fellows to lead an equally challenging, original educational program simultaneously. The legacy of this singular fellowship has been to encourage and advance leaders in the field who are equally as passionate about their own research as they are about sharing that research and their passion for astronomy with students and the public. In this positional paper we address the importance of fellowships like the AAPF to the astronomical profession by identifying the science and educational contributions that Fellows have made to the community. Further, we recommend that fellowships that encourage leading postdoctoral researchers to also become leaders in Astronomy education be continued and expanded.

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.

A growing number of observations indicate that magnetic fields are present among a small fraction of massive O- and B-type stars, yet the origin of these fields remains unclear. Here we present the results of a VLT/FORS1 spectropolarimetric survey of 15 B-type members of the open cluster NGC 3766. We have detected two magnetic B stars in the cluster, including one with a large field of nearly 2 kG, and we find marginal detections of two additional stars. There is no correlation between the observed longitudinal field strengths and the projected rotational velocity, suggesting that a dynamo origin for the fields is unlikely. We also use the Oblique Dipole Rotator model to simulate populations of magnetic stars with uniform or slightly varying magnetic flux on the ZAMS. None of the models successfully reproduces our observed range in B_l and the expected number of field detections, and we rule out a purely fossil origin for the observed fields.

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

As part of an ongoing study to determine the distribution of orbital parameters for massive binaries in the Cygnus OB2 Association, we present the orbital solutions for two new single-lined spectroscopic binaries, MT059 (O8V) & MT258 (O8V), and one double-lined eclipsing binary (Schulte 3). We also constrain the orbital elements of three additional double-lined systems (MT252, MT720, MT771). Periods for all systems range from 1.5-19 days and eccentricities range from 0-0.11. The six new OB binary systems bring the total number of multiple systems within the core region of Cyg OB2 to 11. The current sample does not show evidence for a ``twin-heavy'' binary distribution.

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.

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 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 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 report on high resolution Echelle spectroscopy of 20 giant stars in the Galactic old open clusters NGC 6791 obtained with Hydra at the WIYN telescope. High precision radial velocity allow us to isolate 15 {\it bona fide} cluster members. From 10 of them we derive a global [M/H]=+0.39$\pm$0.05. We therefore confirm that NGC 6791 is extremely metal rich, exhibits a few marginally sub-solar abundance ratios, and within the resolution of our spectra does not show evidences of spread in metal abundance. With these new data we re-derive the cluster fundamental parameters suggesting that it is about 8 Gyr old and 4.3 kpc far from the Sun. The combination of its chemical properties, age, position, and Galactic orbit hardly makes NGC 6791 a genuine Population I open cluster. We discuss possible interpretations of the cluster peculiarities suggesting that the cluster might be what remains of a much larger system, whose initial potential well could have been sufficient to produce high metallicity stars, and which has been depopulated by the tidal field of the Galaxy. Alternatively, its current properties may be explained by the perturbation of the Galactic bar on an object originated well inside the solar ring, where the metal enrichment had been very fast.

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

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

HD 16429 A is a triple star system consisting of a single-lined spectroscopic binary and a widely separated third component, previously identified via speckle interferometry. Here I present the first orbital elements for the unblended spectroscopic binary as well as estimates of the spectral types and relative flux contributions for each visible component based upon a Doppler tomographic reconstruction of their spectra. There are several stars around HD 16429 A, including the nearby Be X-ray binary and microquasar LS I +61 303, which all probably belong to a sub-cluster within the Cas OB6 association.

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.

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.

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.