Radial velocities of massive stars in Cyg OB2

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, then 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.7AU).