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 n
ot 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 s
tatistics 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.
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 bowsh
ocks. 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+43degr3654, 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)
