We present the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), a spatially complete survey of uniformly selected field OB stars that covers the entire star-forming body of the SMC. Using the IMACS multislit spectrograph and MIKE echelle spectrograph on the Magellan telescopes, we obtained spectra of 374 early-type field stars that are at least 28 pc from any other OB candidates. We also obtained spectra of an additional 23 field stars in the SMC bar identified from slightly different photometric criteria. Here, we present the observational catalog of stars in the RIOTS4 survey, including spectral classifications and radial velocities. For three multi-slit fields covering 8% of our sample, we carried out monitoring observations over 9-16 epochs to study binarity, finding a spectroscopic, massive binary frequency of at least $\sim$60% in this subsample. Classical Oe/Be stars represent a large fraction of RIOTS4 (42%), occurring at much higher frequency than in the Galaxy, consistent with expectation at low metallicity. RIOTS4 confirmed a steep upper IMF in the field, apparently caused by the inability of the most massive stars to form in the smallest clusters. Our survey also yields evidence for in-situ field OB star formation, and properties of field emission-line star populations, including sgB[e] stars and classical Oe/Be stars. We also discuss the radial velocity distribution and its relation to SMC kinematics and runaway stars. RIOTS4 presents a first quantitative characterization of field OB stars in an external galaxy, including the contributions of sparse, but normal, star formation; runaway stars; and candidate isolated star formation.

The Cygnus OB2 Association is one of the nearest and largest collections of massive stars in the Galaxy. Situated at the heart of the "Cygnus X" complex of star-forming regions and molecular clouds, its distance has proven elusive owing to the ambiguous nature of kinematic distances along this $\ell\simeq80$ degree sightline and the heavy, patchy extinction. In an effort to refine the three-dimensional geometry of key Cygnus~X constituents, we have measured distances to four eclipsing double-lined OB-type spectroscopic binaries that are probable members of Cyg~OB2. We find distances of $1.33\pm0.17$, $1.32\pm0.07$, $1.44\pm0.18$, and $1.32\pm0.13$ kpc toward MT91~372, MT91~696, CPR2002~A36, and Schulte~3 respectively. We adopt a weighted average distance of 1.33$\pm$0.06~kpc. This agrees well with spectrophotometric estimates for the Association as a whole and with parallax measurements of protostellar masers in the surrounding interstellar clouds, thereby linking the ongoing star formation in these clouds with Cyg~OB2. We also identify Schulte 3C (O9.5V), a 4" visual companion to the 4.75 day binary Schulte~3(A+B), as a previously unrecognized Association member.

We analyze orbital solutions for 48 massive multiple-star systems in the Cygnus OB2 Association, 23 of which are newly presented here, to find that the observed distribution of orbital periods is approximately uniform in log P for P<45 d, but it is not scale-free. Inflections in the cumulative distribution near 6 d, 14, d, and 45 d, suggest key physical scales of ~0.2, ~0.4, and ~1 A.U. where yet-to-be-identified phenomena create distinct features. No single power law provides a statistically compelling prescription, but if features are ignored, a power law with exponent beta = -0.22 provides a crude approximation over P=1.4 -- 2000 d, as does a piece-wise linear function with a break near 45 d. The cumulative period distribution flattens at P > 45 d, even after correction for completeness, indicating either a lower binary fraction or a shift toward low-mass companions. A high degree of similarity (91% likelihood) between the Cyg OB2 period distribution and that of other surveys suggests that the binary properties at P<25 d are determined by local physics of disk/clump fragmentation and are relatively insensitive to environmental and evolutionary factors. Fully 30% of the unbiased parent sample is a binary with period P < 45 d. Completeness corrections imply a binary fraction near 55% for P < 5000 d. The observed distribution of mass ratios 0.2<q<1 is consistent with uniform, while the observed distribution of eccentricities 0.1<e<0.6 is consistent with uniform plus an excess of e ~ 0 systems. We identify six stars, all supergiants, that exhibit aperiodic velocity variations of ~30 km/s attributed to atmospheric fluctuations.

Massive binary stars may constitute a substantial fraction of progenitors to supernovae and gamma-ray bursts, and the distribution of their orbital characteristics holds clues to the formation process of massive stars. As a contribution to securing statistics on OB-type binaries, we report the discovery and orbital parameters for five new systems as part of the Cygnus OB2 Radial Velocity Survey. Four of the new systems (MT070, MT174, MT267, and MT734 (a.k.a. VI Cygni #11) are single-lined spectroscopic binaries while one (MT103) is a double-lined system (B1V+B2V). MT070 is noteworthy as the longest period system yet measured in Cyg OB2, with P=6.2 yr. The other four systems have periods ranging between 4 and 73 days. MT174 is noteworthy for having a probable mass ratio q<0.1, making it a candidate progenitor to a low-mass X-ray binary. These measurements bring the total number of massive binaries in Cyg OB2 to 25, the most currently known in any single cluster or association.

This work provides a statistical analysis of the massive star binary characteristics in the Cygnus OB2 Association using radial velocity information of 114 B3-O3 primary stars and orbital properties for the 24 known binaries. We compare these data to a series of Monte Carlo simulations to infer the intrinsic binary fraction and distributions of mass ratios, periods, and eccentricities. We model the distribution of mass ratio, log-period, and eccentricity as power-laws and find best fitting indices of alpha=0.1+/-0.5, beta=0.2+/-0.4, and gamma=-0.6+/-0.3 respectively. These distributions indicate a preference for massive companions, short periods, and low eccentricities. Our analysis indicates that the binary fraction of the cluster is 44+/-8% if all binary systems are (artificially) assumed to have P<1000 days; if the power-law period distribution is extrapolated to 10^4 years, a plausible upper limit for bound systems, the binary fraction is ~90+/-10%. Of these binary (or higher order) systems, ~45% will have companions close enough to interact during pre- or post-main-sequence evolution (semi-major axis ~/<4.7 AU). The period distribution for P<27 days is not well reproduced by any single power-law owing to an excess of systems with periods around 3-5 days (0.08-0.31 AU) and a relative shortage of systems with periods around 7-14 days (0.14-0.62 AU). We explore the idea that these longer-period systems evolved to produce the observed excess of short-period systems. The best fitting binary parameters imply that secondaries generate, on average, ~16% of the V-band light in young massive populations. This means that photometrically based distance measurements for young massive clusters & associations will be systematically low by ~8% (0.16 mag in the distance modulus) if the luminous contributions of unresolved secondaries are not taken into account.

We report the discovery and orbital solutions for two new OB binaries in the Cygnus OB2 Association, MT311 (B2V+B3V) and MT605 (B0.5V+B2.5:V). We also identify the system MT429 as a probable triple system consisting of a tight eclipsing 2.97 day B3V+B6V pair and a B0V at a projected separation of 138 AU. We further provide the first spectroscopic orbital solutions to the eclipsing, double-lined, O-star binary MT696 (O9.5V+B1:V), the double-lined, early B binary MT720 (B0-1V+B1-2V), and the double-lined, O-star binary MT771 (O7V+O9V). These systems exhibit orbital periods between 1.5 days and 12.3 days, with the majority having P<6 days. The two new binary discoveries and six spectroscopic solutions bring the total number of known massive binaries in the central region of the Cygnus OB2 Association to 20, with all but two having full orbital solutions.

After the first detection of its binary nature, the spectroscopic monitoring of the non-thermal radio emitter Cyg OB2 #9 (P=2.4yrs) has continued, doubling the number of available spectra of the star. Since the discovery paper of 2008, a second periastron passage has occurred in February 2009. Using a variety of techniques, the radial velocities could be estimated and a first, preliminary orbital solution was derived from the HeI5876 line. The mass ratio appears close to unity and the eccentricity is large, 0.7--0.75. X-ray data from 2004 and 2007 are also analyzed in quest of peculiarities linked to binarity. The observations reveal no large overluminosity nor strong hardness, but it must be noted that the high-energy data were taken after the periastron passage, at a time where colliding wind emission may be low. Some unusual X-ray variability is however detected, with a 10% flux decrease between 2004 and 2007. To clarify their origin and find a more obvious signature of the wind-wind collision, additional data, taken at periastron and close to it, are needed.

We use mid-IR images from the Spitzer Cygnus~X Legacy Survey to search for stellar bowshocks, a signature of early type "runaway" stars with high space velocities. We identify ten arc-shaped nebulae containing centrally located stars as candidate bowshocks. New spectroscopic observations of five stars show that all are late O to early B dwarfs. Our morphologically selected sample of bowshock candidates encompasses diverse physical phenomena. Three of the stars appear to be pre-main-sequence objects on the basis of rising SEDs in the mid-IR, and their nebulae may be photon-dominated regions (PDRs). Four objects have ambiguous classification. These may be partial dust shells or bubbles. We conclude that three of the objects are probable bowshocks, based on their morphological similarity to analytic prescriptions. Their nebular morphologies reveal no systematic pattern of orientations that might indicate either a population of stars ejected from or large-scale hydrodynamic outflows from Cyg OB2. The fraction of runaways among OB stars near Cyg OB2 identified either by radial velocity or bowshock techniques is ~0.5%, much smaller than the 8% estimated among field OB stars. We also obtained a heliocentric radial velocity for the previously known bowshock star, BD+43\degr3654, of -66.2+/-9.4 km/s, solidifying its runaway status and implying a space velocity of 77+/-10 km/s. We use the principles of momentum-driven bowshocks to arrive at a novel method for estimating stellar mass loss rates. Derived mass loss rates range between 10^-7 and few x10^-6 solar masses/yr for the three O5V -- ~B2V stars identified as generating bowshocks. These values are at the upper range of, but broadly consistent with, estimates from other methods. (Abridged)

We present the orbital solutions for four OB spectroscopic binaries, MT145, GSC 03161-00815, 2MASS J20294666+4105083, and Schulte 73, and the partial orbital solution to the B spectroscopic binary, MT372, as part of an ongoing study to determine the distribution of orbital parameters for massive binaries in the Cygnus OB2 Association. MT145 is a new, single-lined, moderately eccentric (e=0.291+/-0.009) spectroscopic binary with period of 25.140+/-0.008 days. GSC 03161-00815 is a slightly eccentric (e=0.10+/-0.01), eclipsing, interacting and double-lined spectroscopic binary with a period of 4.674+/-0.004 days. 2MASS J20294666+4105083 is a moderately eccentric (e=0.273+/-0.002) double-lined spectroscopic binary with a period of 2.884+/-0.001 days. Schulte 73 is a slightly eccentric (e=0.169+/-0.009), double-lined spectroscopic binary with a period of 17.28+/-0.03 days and the first "twin" in our survey with a mass ratio of q=0.99+/-0.02. MT372 is a single-lined, eclipsing system with a period of 2.228 days and low eccentricity (e~0). Of the now 18 known OB binaries in Cyg OB2, 14 have periods and mass ratios. Emerging evidence also shows that the distribution of log(P) is flat and consistent with Oepik's Law.

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 conducted a radial velocity survey of the Cygnus OB2 Association over a 6 year (1999 - 2005) time interval to search for massive close binaries. During this time we obtained 1139 spectra on 146 OB stars to measure mean systemic radial velocities and radial velocity variations. We spectroscopically identify 73 new OB stars for the first time, the majority of which are likely to be Association members. Spectroscopic evidence is also presented for a B3Iae classification and temperature class variation (B3 - B8) on the order of 1 year for Cygnus OB2 No. 12. Calculations of the intial mass function with the current spectroscopic sample yield Gamma = -2.2 +/- 0.1. Of the 120 stars with the most reliable data, 36 are probable and 9 are possible single-lined spectroscopic binaries. We also identify 3 new and 8 candidate double-lined spectroscopic binaries. These data imply a lower limit on the massive binary fraction of 30% - 42%. The calculated velocity dispersion for Cygnus OB2 is 2.44 +/- km/s, which is typical of open clusters. No runaway OB stars were found.