The stellar mass-luminosity relation is poorly constrained by observations for high mass stars. We describe our program to find eclipsing massive binaries in the Magellanic Clouds using photometry of regions rich in massive stars, and our spectroscopic follow-up to obtain radial velocities and orbits. Our photometric campaign identified 48 early-type periodic variables, of which only 15 (31%) were found as part of the microlensing surveys. Spectroscopy is now complete for 17 of these systems, and in this paper we present analysis of the first two, LMC 172231 and ST2-28, simple detached systems of late-type O dwarfs of relatively modest masses. Our orbit analysis yields very precise masses (2%) and we use tomography to separate the components and determine effective temperatures by model fitting, necessary for determining accurate (0.05-0.07 dex) bolometric luminosities in combination with the light-curve analysis. Our approach allows more precise comparisons with evolutionary theory than previously possible. To our considerable surprise, we find a small, but significant, systematic discrepancy: all of the stars are slightly under-massive, by typically 11% (or over-luminous by 0.2 dex) compared to that predicted by the evolutionary models. We examine our approach for systematic problems, but find no satisfactory explanation. The discrepancy is in the same sense as the long-discussed and elusive discrepancy between the masses measured from stellar atmosphere analysis with the stellar evolutionary models, and might suggest that either increased rotation or convective overshooting is needed in the models. Additional systems will be discussed in future papers of this series, and will hopefully confirm or refute this trend.

We present projected rotational velocity values for 97 Galactic, 55 SMC, and 106 LMC O-B type stars from archival FUSE observations. The evolved and unevolved samples from each environment are compared through the Kolmogorov-Smirnov test to determine if the distribution of equatorial rotational velocities is metallicity dependent for these massive objects. Stellar interior models predict that massive stars with SMC metallicity will have significantly reduced angular momentum loss on the main sequence compared to their Galactic counterparts. Our results find some support for this prediction but also show that even at Galactic metallicity, evolved and unevolved massive stars have fairly similar fractions of stars with large V sin i values. Macroturbulent broadening that is present in the spectral features of Galactic evolved massive stars is lower in the LMC and SMC samples. This suggests the processes that lead to macroturbulence are dependent upon metallicity.

We present the results of a Doppler tomographic reconstruction of the UV spectra of the double-lined, O binary HD 100213 based on observations made with the International Ultraviolet Explorer (IUE). We used cross-correlation methods to obtain radial velocities, confirm the orbital elements, estimate the UV flux ratio, and determine projected rotational velocities. The individual component spectra are classified as O7 V + O8 V using UV criteria defined by Penny, Gies, & Bagnuolo. We present a model fit of the eclipsing light curve from observations from the Hipparcos satellite and published observations of Andersen & Gronbech. We derive an orbital inclination, i=77.7deg+/-1.0deg. This analysis indicates that both stars are currently experiencing Roche lobe overflow (RLOF), which confirms earlier results that this is one of only a few massive contact binaries. Our derived masses, Mp/Msolar=16.7+/-0.4 and Ms/Msolar=10.4+/-0.4, are significantly lower than those computed from the Doppler shifts of lines in the optical spectrum. We suggest that the difference occurs because mutual irradiation decreases the upper atmospheric temperature gradient in the inward-facing hemispheres of both stars, which makes lower excitation lines appear weaker there and shifts their center of light away from the center of mass. We compare the current state of HD 100213 with predicted outcomes of massive close binary evolutionary models, and suggest that the system is currently in a very slow case AA mass transfer stage.

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.

Recent theoretical predictions for low metallicity massive stars predict that these stars should have drastically reduced equatorial winds (mass loss) while on the main sequence, and as such should retain most of their angular momentum. Observations of both the Be/(B+Be) ratio and the blue-to-red supergiant ratio appear to have a metallicity dependence that may be caused by high rotational velocities. We have analyzed 39 archival Hubble Space Telescope Imaging Spectrograph (STIS), high resolution, ultraviolet spectra of O-type stars in the Magellanic Clouds to determine their projected rotational velocities V sin i. Our methodology is based on a previous study of the projected rotational velocities of Galactic O-type stars using International Ultraviolet Explorer (IUE) Short Wavelength Prime (SWP) Camera high dispersion spectra, which resulted in a catalog of V sin i values for 177 O stars. Here we present complementary V sin i values for 21 Large Magellanic Cloud and 22 Small Magellanic Cloud O-type stars based on STIS and IUE UV spectroscopy. The distribution of V sin i values for O type stars in the Magellanic Clouds is compared to that of Galactic O type stars. Despite the theoretical predictions and indirect observational evidence for high rotation, the O type stars in the Magellanic Clouds do not appear to rotate faster than their Galactic counterparts.

We present the first orbital elements for the massive close binary, HD 101131, one of the brightest objects in the young, open cluster, IC 2944. This system is a double-lined spectroscopic binary in an elliptical orbit with a period of 9.64659 +/- 0.00012 days. It is a young system of unevolved stars (approximately 2 million years old) that are well within their critical Roche surfaces. We use a Doppler tomography algorithm to reconstruct the individual component optical spectra, and we apply well known criteria to arrive at classifications of O6.5 V((f)) and O8.5 V for the primary and secondary, respectively. We compare the reconstructed spectra of the components to single star, spectrum standards to determine a flux ratio of f_2/f_1 = 0.55 +/- 0.08 in the V-band. Both components are rotating faster than synchronously. We estimate the temperatures and luminosities of the components from the observed spectral classifications, composite V magnitude, and cluster distance modulus. The lower limits on the masses derived from the orbital elements and the lack of eclipses are 25 solar masses and 14 solar masses for the primary and secondary, respectively. These limits are consistent with the larger masses estimated from the positions of the stars in the Hertzsprung-Russell diagram and evolutionary tracks for single stars.

We present the first orbital elements for the massive close binary, HD 115071, a double-lined spectroscopic binary in a circular orbit with a period of 2.73135 +/- 0.00003 days. The orbital semiamplitudes indicate a mass ratio of M_2/M_1 = 0.58 +/- 0.02 and yet the stars have similar luminosities. We used a Doppler tomography algorithm to reconstruct the individual component optical spectra, and we applied well known criteria to arrive at classifications of O9.5 V and B0.2 III for the primary and secondary, respectively. We present models of the Hipparcos light curve of the ellipsoidal variations caused by the tidal distortion of the secondary, and the best fit model for a Roche-filling secondary occurs for an inclination of i = 48.7 +/- 2.1 degrees. The resulting masses are 11.6 +/- 1.1 and 6.7 +/- 0.7 solar masses for the primary and secondary, respectively, so that both stars are very overluminous for their mass. The system is one of only a few known semi-detached, Algol-type binaries that contain O-stars. We suggest that the binary has recently emerged from extensive mass transfer (possibly through a delayed contact and common envelope process).

Double-lined spectroscopic orbital elements have recently been found for the central binary in the massive triple, delta Orionis A based on radial velocities from cross-correlation techniques applied to IUE high dispersion spectra and He I 6678 spectra obtained at Kitt Peak. The primary and secondary velocity amplitudes were found to be 94.9 +/- 0.6 km/s and 186 +/- 9 km/s respectively. Tomographic reconstructions of the primary and secondary stars' spectra confirm the O9.5 II classification of the primary and indicate a B0.5 III type for the secondary. The widths of the UV cross-correlation functions are used to estimate the projected rotational velocities, Vsin i = 157 +/- 6 km/s and 138 +/- 16 km/s for the primary and secondary, respectively implying that both stars rotate faster than their orbital motion. We used the spectroscopic results to make a constrained fit of the Hipparcos light curve of this eclipsing binary, and the model fits limit the inclination to the range between 67 and 77 degrees. The i = 67 degrees solution, which corresponds to a near Roche-filling configuration, results in a primary mass of 11.2 solar masses and a secondary mass of 5.6 solar masses, both of which are substantially below the expected masses for stars of their luminosity. This binary may have experienced a mass ratio reversal caused by Case A Roche lobe overflow, or the system may have suffered extensive mass loss through a binary interaction, perhaps during a common envelope phase, in which most of the primary's mass was lost from the system rather than transferred to the secondary.

We present radial velocity and photometry for four early-type, massive double-lined spectroscopic binaries in the R136 cluster. Three of these systems are eclipsing, allowing orbital inclinations to be determined. One of these systems, R136-38 (O3 V + O6 V), has one of the highest masses ever measured, 57 Mo, for the primary. Comparison of our masses with those derived from standard evolutionary tracks shows excellent agreement. We also identify five other light variables in the R136 cluster which are worthy of follow-up study.