We present a narrow-band imaging survey of the Large Magellanic Cloud, designed to isolate the C II $\lambda\lambda$7231, 7236 emission lines in objects as faint as $m_{\lambda7400}\sim18$. The work is motivated by the recent serendipitous discovery in the LMC of the first confirmed extragalactic [WC11] star, whose spectrum is dominated by C II emission, and the realization that the number of such objects is currently largely unconstrained. The survey, which imaged $\sim$50$~$deg$^2$ using on-band and off-band filters, will significantly increase the total census of these rare stars. In addition, each new LMC [WC] star has a known luminosity, a quantity quite uncertain in the Galactic sample. Multiple known C II emitters were easily recovered, validating the survey design. We find 38 new C II emission candidates; spectroscopy of the complete sample will be needed to ascertain their nature. In a preliminary spectroscopic reconnaissance, we observed three candidates, finding C II emission in each. One is a new [WC11]. Another shows both the narrow C II emission lines characteristic of a [WC11], but also broad emission of C IV, O V, and He II characteristic of a much hotter [WC4] star; we speculate that this is a binary [WC]. The third object shows weak C II emission, but the spectrum is dominated by a dense thicket of strong absorption lines, including numerous O II transitions. We conclude it is likely an unusual hot, hydrogen-poor post-AGB star, possibly in transition from [WC] to white dwarf. Even lacking a complete spectroscopic program, we can infer that late [WC] stars do not dominate the central stars of LMC planetary nebulae, and that the detected C II emitters are largely of an old population.

The mass-loss rates of red supergiant stars (RSGs) are poorly constrained by direct measurements, and yet the subsequent evolution of these stars depends critically on how much mass is lost during the RSG phase. In 2012 the Geneva evolutionary group updated their mass-loss prescription for RSGs with the result that a 20 solar mass star now loses 10x more mass during the RSG phase than in the older models. Thus, higher mass RSGs evolve back through a second yellow supergiant phase rather than exploding as Type II-P supernovae, in accord with recent observations (the so-called "RSG Problem"). Still, even much larger mass-loss rates during the RSG phase cannot be ruled out by direct measurements of their current dust-production rates, as such mass-loss is episodic. Here we test the models by deriving a luminosity function for RSGs in the nearby spiral galaxy M31 which is sensitive to the total mass loss during the RSG phase. We carefully separate RSGs from asymptotic giant branch stars in the color-magnitude diagram following the recent method exploited by Yang and collaborators in their Small Magellanic Cloud studies. Comparing our resulting luminosity function to that predicted by the evolutionary models shows that the new prescription for RSG mass-loss does an excellent job of matching the observations, and we can readily rule out significantly larger values.

Wolf-Rayet stars (WRs) represent the end of a massive star's life as it is about to turn into a supernova. Obtaining complete samples of such stars across a large range of metallicities poses observational challenges, but presents us with an exacting way to test current stellar evolutionary theories. A technique we have developed and refined involves interference filter imaging combined with image subtraction and crowded-field photometry. This helps us address one of the most controversial topics in current massive star research: the relative importance of binarity in the evolution of massive stars and formation of WRs. Here we discuss the current state of the field, including how the observed WR populations match with the predictions of both single and binary star evolutionary models. We end with what we believe are the most important next steps in WR research.

We report the discovery of 11 newly found quasars behind the stellar disks of the spiral galaxies M31 and M33 in the fields covered by the Local Group Galaxy Survey. Their redshifts range from 0.37 to 2.15. Most are X-ray, UV, and IR sources. We also report the discovery of 5 normal background galaxies. Most of these objects were observed owing to their anomalous colors, as part of a program (reported elsewhere) to confirm spectroscopically candidate red supergiant plus B star binaries; others were discovered as part of our identification of early-type massive stars based upon their optical colors. There are 15 previously known quasars in the same fields, for a grand total of 26, 15 behind M31 and 11 behind M33. Of these, only eight were discovered as part of surveys for quasars; the rest were found accidentally. The quasars are well distributed in the M31 and M33 fields, except for the inner regions, and have the potential for being good probes of the interstellar medium in these stellar disks, as well as serving as zero-point calibrators for Gaia parallaxes.

The percentage of massive main sequence OB stars in binary systems is thought to be as high as 100%. However, very few Galactic binary red supergiants (RSGs) have been identified, despite the fact that these stars are the evolved descendants of OB stars. As shown in our recent paper, binary RSGs will likely have B-type companions, as dictated by stellar evolution considerations. Such a system will have a very unique photometric signature due to the shape of the spectral energy distribution. Using photometric cutoffs it should therefore be possible to detect candidate RSG+B star binary systems. Here we present our spectoscopic follow-up observations of such candidates. Out of our initial list of 280 candidates in M31 and M33, we observed 149 and confirmed 63 as newly discovered RSG+B star binary systems. Additional spectra of four candidate systems in the Small Magellanic Cloud confirmed all of them as new RSG+B star binaries including the first known RSG+Be star system. By fitting BSTAR06 and MARCS model atmospheres to the newly-obtained spectra we place estimates on the temperatures and subtypes of both the B stars and RSGs. Overall, we have found 87 new RSG+B star binary systems in M31, M33 and the Small and Large Magellanic Clouds. Our future studies are aimed at determining the binary fraction of RSGs.

In the standard view of massive star evolution, luminous blue variables (LBVs) are transitional objects between the most massive O-type stars and Wolf-Rayet (WR) stars. With short lifetimes, these stars should all be found near one another. A recent study of LBVs in the Large Magellanic Cloud (LMC) found instead that LBVs are considerably more isolated than either O-type stars or WRs, with a distribution intermediate between that of the WRs and red supergiants (RSGs). A similar study, using a more restricted sample of LBVs, reached the opposite conclusion. Both studies relied upon the distance to the nearest spectroscopically identified O-type star to define the degree of isolation. However, our knowledge of the spectroscopic content of the LMC is quite spotty. Here we re-examine the issue using carefully defined photometric criteria to select the highest mass unevolved stars ("bright blue stars," or BBSs), using spatially complete photometric catalogs of the LMC, M31, and M33. Our study finds that the LBVs are no more isolated than BBSs or WRs. This result holds no matter which sample of LBVs we employ. A statistical test shows that we can rule out the LBVs having the same distribution as the RSGs, which are about 2x more isolated. We demonstrate the robustness of our results using the second-closest neighbor. Furthermore, the majority of LBVs in the LMC are found in or near OB associations as are the BBS and WRs; the RSGs are not. We conclude that the spatial distribution of LBVs therefore is consistent with the standard picture of massive star evolution.

New evidence provided by the Gaia satellite places the location of the runaway star J01020100-7122208 in the halo of the Milky Way (MW) rather than in the Small Magellanic Cloud as previously thought. We conduct a reanalysis of the star's physical and kinematic properties, which indicates that the star may be an even more extraordinary find than previously reported. The star is a 180 Myr old 3-4 Mo G5-8 bright giant, with an effective temperature of 4800+/-100 K, a metallicity of {Fe/H]=-0.5, and a luminosity of log L/Lo=2.70+/-0.20 dex. A comparison with evolutionary tracks identifies the star as being in a giant or early asymptotic giant branch stage. The proper motion, combined with the previously known radial velocity, yields a total Galactocentric space velocity of 296 km/s. The star is currently located 6.4 kpc below the plane of the Milky Way, but our analysis of its orbit shows it passed through the disk ~25 Myr ago. The star's metallicity and age argue against it being native to the halo, and we suggest that the star was likely ejected from the disk. We discuss several ejection mechanisms, and conclude that the most likely scenario is ejection by the Milky Way's central black hole based upon our analysis of the star's orbit. The identification of the large radial velocity of J01020100-7122208 came about as a happenstance of it being seen in projection with the SMC, and we suggest that many similar objects may be revealed in Gaia data.

With the exception of a few well-known and studied systems, the binary population of red supergiants (RSGs) remains relatively uncharacterized. Famous systems such as VV Cep, 31 Cyg and zeta Aur contain RSG + B star binaries and here we explore whether B stars are the main type of companion we expect from an evolutionary point of view. Using the Geneva evolutionary models we find that this is indeed the case. However, few such systems are known, and we use model spectra to determine how easy such binaries would be to detect observationally. We find that it should be quite difficult to hide a B-type companion given a reasonable signal-to-noise in the optical / blue portion of the spectrum. We next examine spectra of Magellanic Cloud RSGs and newly acquired spectra of Galactic RSGs looking for new systems and refining our conclusions about what types of stars could be hidden in the spectra. Finally, we develop a set of photometric criteria that can help select likely binaries in the future without the overhead of large periodic or spectroscopic surveys.

We summarize the results of our four year survey searching for Wolf-Rayet (WR) stars in the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). Over the course of this survey we've discovered 15 new WRs and 12 Of-type stars. In this last year we discovered two rare Of-type stars: an O6.5f?p and an O6nfp in addition to the two new Of?p stars discovered in our first year and the three Onfp stars discovered in our second and third years. However, even more exciting was our discovery of a new type of WR, ones we are calling WN3/O3s due to their spectroscopic signatures. We describe the completeness limits of our survey and demonstrate that we are sensitive to weak-lined WRs several magnitudes fainter than any we have discovered, arguing that there is not a population of fainter WRs waiting to be discovered. We discuss the nature of the WN3/O3s, summarizing the results of our extensive spectroscopy and modeling. We also examine the important claim made by others that the WN3/O3s are isolated compared to other massive stars. We find that if we use a more complete sample of reference massive stars, the WN3/O3s show the same spatial distribution as other early WNs, consistent with a common origin. Finally, we use this opportunity to present the "Fifth Catalog of LMC Wolf-Rayet Stars," which includes revised coordinates and updated spectral types for all 154 known LMC WRs.

We recently discovered a yellow supergiant (YSG) in the Small Magellanic Cloud (SMC) with a heliocentric radial velocity of ~300 km/s which is much larger than expected for a star in its location in the SMC. This is the first runaway YSG ever discovered and only the second evolved runaway star discovered in a different galaxy than the Milky Way. We classify the star as G5-8I, and use de-reddened broad-band colors with model atmospheres to determine an effective temperature of 4700+/-250K, consistent with what is expected from its spectral type. The star's luminosity is then L/Lo ~ 4.2+/-0.1, consistent with it being a ~30Myr 9Mo star according to the Geneva evolution models. The star is currently located in the outer portion of the SMC's body, but if the star's transverse peculiar velocity is similar to its peculiar radial velocity, in 10Myr the star would have moved 1.6 degrees across the disk of the SMC, and could easily have been born in one of the SMC's star-forming regions. Based on its large radial velocity, we suggest it originated in a binary system where the primary exploded as a supernovae thus flinging the runaway star out into space. Such stars may provide an important mechanism for the dispersal of heavier elements in galaxies given the large percentage of massive stars that are runaways. In the future we hope to look into additional evolved runaway stars that were discovered as part of our other past surveys.

The nearby galaxies of the Local Group can act as our laboratories in helping to bridge the gap between theory and observations. In this review we will describe the complications of identifying samples of OB stars, yellow and red supergiants, and Wolf-Rayet stars, and what we have so far learned from these studies.

As part of a search for Wolf-Rayet (WR) stars in the Magellanic Clouds, we have discovered a new type of WR star in the Large Magellanic Cloud (LMC). These stars have both strong emission lines, as well as He ii and Balmer absorption lines and spectroscopically resemble a WN3 and O3V binary pair. However, they are visually too faint to be WN3+O3V binary systems. We have found nine of these WN3/O3s, making up ~6% of the population of LMC WRs. Using CMFGEN, we have successfully modeled their spectra as single stars and have compared the physical parameters with those of more typical LMC WNs. Their temperatures are around 100,000 K, a bit hotter than the majority of WN stars (by around 10,000 K) although a few hotter WNs are known. The abundances are what you would expect for CNO equilibrium. However, most anomalous are their mass-loss rates which are more like that of an O-type star than a WN star. While their evolutionary status is uncertain, their low mass-loss rates and wind velocities suggest that they are not products of homogeneous evolution. It is possible instead that these stars represent an intermediate stage between O stars and WNs. Since WN3/O3 stars are unknown in the Milky Way, we suspect that their formation depends upon metallicity, and we are investigating this further by a deep survey in M33, which posses a metallicity gradient.

We summarize here recent work in identifying and characterizing red supergiants (RSGs) in the galaxies of the Local Group.

As part of a multi-year survey for Wolf-Rayet stars in the Magellanic Clouds, we have discovered a new type of Wolf-Rayet star with both strong emission and absorption. While one might initially classify these stars as WN3+O3V binaries based on their spectra, such a pairing is unlikely given their faint visual magnitudes. Spectral modeling suggests effective temperatures and bolometric luminosities similar to those of other early-type LMC WNs but with mass-loss rates that are three to five times lower than expected. They additionally retain a significant amount of hydrogen, with nitrogen at its CNO-equilibrium value (10x enhanced). Their evolutionary status remains an open question. Here we discuss why these stars did not evolve through quasi-homogeneous evolution. Instead we suggest that based on a link with long-duration gamma ray bursts, they may form in lower metallicity environments. A new survey in M33, which has a large metallicity gradient, is underway.

For the past three years we have been conducting a survey for WR stars in the Large and Small Magellanic Clouds (LMC, SMC). Our previous work has resulted in the discovery of a new type of WR star in the LMC, which we are calling WN3/O3. These stars have the emission-line properties of a WN3 star (strong N V but no N IV), plus the absorption-line properties of an O3 star (Balmer hydrogen plus Pickering He II but no He I). Yet these stars are 15x fainter than an O3 V star would be by itself, ruling out these being WN3+O3 binaries. Here we report the discovery of two more members of this class, bringing the total number of these objects to 10, 6.5% of the LMC's total WR population. The optical spectra of nine of these WN3/O3s are virtually indistinguishable from each other, but one of the newly found stars is significantly different, showing a lower excitation emission and absorption spectrum (WN4/O4-ish). In addition, we have newly classified three unusual Of-type stars, including one with a strong C III 4650 line, and two rapidly rotating "Oef" stars. We also "rediscovered" a low mass x-ray binary, RX J0513.9-6951, and demonstrate its spectral variability. Finally, we discuss the spectra of ten low priority WR candidates that turned out not to have He II emission. These include both a Be star and a B[e] star.

We describe our spectroscopic follow-up to the Local Group Galaxy Survey (LGGS) photometry of M31 and M33. We have obtained new spectroscopy of 1895 stars, allowing us to classify 1496 of them for the first time. Our study has identified many foreground stars, and established membership for hundreds of early- and mid-type supergiants. We have also found 9 new candidate Luminous Blue Variables and a previously unrecognized Wolf-Rayet star. We republish the LGGS M31 and M33 catalogs with improved coordinates and including spectroscopy from the literature and our new results. The spectroscopy in this paper is responsible for the vast majority of the stellar classifications in these two nearby spiral neighbors. The most luminous (and hence massive) of the stars in our sample are early-type B supergiants, as expected; the more massive O stars will be fainter visually, and thus mostly remain unobserved so far. The majority of the unevolved stars in our sample are in the 20-40Mo range.

As part of our ongoing Wolf-Rayet (WR) Magellanic Cloud survey, we have discovered 13 new WRs. However, the most exciting outcome of our survey is not the number of new WRs, but their unique characteristics. Eight of our discoveries appear to belong to an entirely new class of WRs. While one might naively classify these stars as WN3+O3V binaries, such a pairing is unlikely. Preliminary CMFGEN modeling suggests physical parameters similar to early-type WNs in the Large Magellanic Cloud except with mass-loss rates three to five times lower and slightly higher temperatures. The evolution status of these stars remains an open question.

We summarize past and current surveys for Wolf-Rayet stars among the Local Group galaxies, emphasizing both the how and the why. Such studies are invaluable for helping us learn about massive star evolution, and for providing sensitive tests of the stellar evolution models. But for such surveys to be useful, the completeness limits must be well understood. We illustrate that point in this review by following the "evolution" of the observed WC/WN ratio in nearby galaxies. We end by examining our new survey for WR stars in the Magellanic Clouds, which has revealed a new type of WN star, never before seen.

The numbers and types of evolved massive stars found in nearby galaxies provide an exacting test of stellar evolution models. Because of their proximity and rich massive star populations, the Magellanic Clouds have long served as the linchpins for such studies. Yet the continued accidental discoveries of Wolf-Rayet (WR) stars in these systems demonstrate that our knowledge is not as complete as usually assumed. Therefore, we undertook a multi-year survey for WRs in the Magellanic Clouds. Our results from our first year (reported previously) confirmed nine new LMC WRs. Of these, six were of a type never before recognized, with WN3-type emission combined with O3-type absorption features. Yet these stars are 2-3 magnitudes too faint to be WN3+O3 V binaries. Here we report on the second year of our survey, including the discovery of four more WRs, two of which are also WN3/O3s, plus two "slash" WRs. This brings the total of LMC WRs known to 152, 13 (8.2%) of which were found by our survey, which is now 60% complete. We find that the spatial distribution of the WN3/O3s are similar to that of other WRs in the LMC, suggesting that they are descended from the same progenitors. We call attention to the fact that five of the 12 known SMC WRs may in fact be similar WN3/O3s rather than the binaries they have often assumed to be. We also discuss our other discoveries: a newly found Onfp-type star, and a peculiar emission-line object. Finally, we consider the completeness limits of our survey.

Our new survey for Wolf-Rayet stars in the Magellanic Clouds is only 15% complete but has already found 9 new WRs in the LMC. This suggests that the total WR population in the LMC may be underestimated by 10-40%. Eight of the nine are WNs, demonstrating that the "observed" WC to WN ratio is too large, and is biased towards WC stars. The ninth is another rare WO star, the second we have found in the LMC in the past two years. Five (and possibly six) of the 8 WNs are of a new class of WRs, which pose a significant challenge to our understanding. Naively we would classify these stars as "WN3+O3V," but there are several reasons why such a pairing is unlikely, not the least of which is that the absolute visual magnitudes of these stars are faint, with M_V~ -2.3 to -3.1. We have performed a preliminary analysis with CMFGEN, and we find that (despite the faint visual magnitudes) the bolometric luminosities of these stars are normal for early-type WNs. Our fitting suggests that these stars are evolved, with significantly enriched N and He. Their effective temperatures are also normal for early-type WNs. What is unusual about these stars is that they have a surprisingly small mass-loss rate compared to other early-type WNs. How these stars got to be the way they are (single star evolution? binary evolution?) remains an open question. For now, we are designating this class as WN3/O3, in analogy to the late-type WN "slash" stars.

Massive star evolutionary models generally predict the correct ratio of WC-type and WN-type Wolf-Rayet stars at low metallicities, but underestimate the ratio at higher (solar and above) metallicities. One possible explanation for this failure is perhaps single-star models are not sufficient and Roche-lobe overflow in close binaries is necessary to produce the "extra" WC stars at higher metallicities. However, this would require the frequency of close massive binaries to be metallicity dependent. Here we test this hypothesis by searching for close Wolf-Rayet binaries in the high metallicity environments of M31 and the center of M33 as well as in the lower metallicity environments of the middle and outer regions of M33. After identifying ~100 Wolf-Rayet binaries based on radial velocity variations, we conclude that the close binary frequency of Wolf-Rayets is not metallicity dependent and thus other factors must be responsible for the overabundance of WC stars at high metallicities. However, our initial identifications and observations of these close binaries have already been put to good use as we are currently observing additional epochs for eventual orbit and mass determinations.

Over the years, directed surveys and incidental spectroscopy have identified 12 Wolf-Rayet (WR) stars in the SMC and 139 in the LMC, numbers which are often described as "essentially} complete." Yet, new WRs are discovered in the LMC almost yearly. We have therefore initiated a new survey of both Magellanic Clouds using the same interference-filter imaging technique previously applied to M31 and M33. We report on our first observing season, in which we have successfully surveyed ~15$% of our intended area of the SMC and LMC. Spectroscopy has confirmed 9 newly found WRs in the LMC (a 6% increase), including one of WO-type, only the third known in that galaxy and the second to be discovered recently. The other eight are WN3 stars that include an absorption component. In two, the absorption is likely from an O-type companion, but the other six are quite unusual. Five would be classified naively as "WN3+O3 V," but such a pairing is unlikely given the rarity of O3 stars, the short duration of this phase (which is incommensurate with the evolution of a companion to a WN star), and because these stars are considerably fainter than O3 V stars. The sixth star may also fall into this category. CMFGEN modeling suggests these stars are hot, bolometrically luminous, and N-rich like other WN3 stars, but lack the strong winds that characterize WNs. Finally, we discuss two rare Of?p stars and four Of supergiants we found, and propose that the B[e] star HD 38489 may have a WN companion.

The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified.

The model atmosphere programs FASTWIND and CMFGEN are both elegantly designed to perform non-LTE analyses of the spectra of hot massive stars, and include sphericity and mass-loss. The two codes differ primarily in their approach towards line blanketing, with CMFGEN treating all of the lines in the co-moving frame and FASTWIND taking an approximate approach which speeds up execution times considerably. Although both have been extensively used to model the spectra of O-type stars, no studies have used the codes to independently model the same spectra of the same stars and compare the derived physical properties. We perform this task on ten O-type stars in the Magellanic Clouds. For the late-type O supergiants, both CMFGEN and FASTWIND have trouble fitting some of the He I lines, and we discuss causes and cures. We find that there is no difference in the average effective temperatures found by the two codes for the stars in our sample, although the dispersion is large, due primarily to the various difficulties each code has with He I. The surface gravities determined using FASTWIND are systematically lower by 0.12 dex compared to CMFGEN, a result we attribute to the better treatment of electron scattering by CMFGEN. This has implications for the interpretation of the origin of the so-called mass discrepancy, as the masses derived by FASTWIND are on average lower than inferred from stellar evolutionary models, while those found by CMFGEN are in better agreement.

While observing OB stars within the most crowded regions of the Large Magellanic Cloud, we happened upon a new Wolf-Rayet star in Lucke-Hodge 41, the rich OB association that contains S Doradus and numerous other massive stars. At first glance the spectrum resembled that of a WC4 star, but closer examination showed strong OVI 3811, 34 lines, leading us to classify it as a WO4. This is only the second known WO in the LMC, and the first known WO4 (the other being a WO3). This rarity is to be expected due to these stars' short lifespans as they represent the most advanced evolutionary stage in a massive star's lifetime before exploding as SNe. This discovery shows that while the majority of WRs within the LMC have been discovered, there may be a few WRs left to be found.

Wolf-Rayet stars are evolved massive stars, and the relative number of WC-type and WN-type WRs should vary with the metallicity of the host galaxy, providing a sensitive test of stellar evolutionary theory. However, past studies of the WR content of M31 have been biased towards detecting WC stars, as their emission line signatures are much stronger than those of WNs. Here we present the results of a survey covering all of M31's optical disk (2.2 deg^2), with sufficient sensitivity to detect the weaker-lined WN-types. We identify 107 newly found WR stars, mostly of WN-type. This brings the total number of spectroscopically confirmed WRs in M31 to 154, a number we argue is complete to ~95%, except in regions of unusually high reddening. This number is consistent with what we expect from the integrated Halpha luminosity compared to that of M33. The majority of these WRs formed in OB associations around the Population I ring, although 5% are truly isolated. Both the relative number of WC to WN-type stars as well as the WC subtype distribution suggest that most WRs exist in environments with higher-than-solar metallicities, which is consistent with studies of M31's metallicity. Although the WC-to-WN ratio we find for M31 is much lower than that found by previous studies, it is still higher than what the Geneva evolutionary models predict. This may suggest that Roche-lobe overflow produces the excess of WC stars observed at high metallicity, or that the assumed rotational velocities in the models are too high.

Due to their transitionary nature, yellow supergiants provide a critical challenge for evolutionary modeling. Previous studies within M31 and the SMC show that the Geneva evolutionary models do a poor job at predicting the lifetimes of these short-lived stars. Here we extend this study to the LMC while also investigating the galaxy's red supergiant content. This task is complicated by contamination by Galactic foreground stars that color and magnitude criteria alone cannot weed out. Therefore, we use proper motions and the LMC's large systemic radial velocity (\sim278 km/s) to separate out these foreground dwarfs. After observing nearly 2,000 stars, we identified 317 probable yellow supergiants, 6 possible yellow supergiants and 505 probable red supergiants. Foreground contamination of our yellow supergiant sample was \sim80%, while that of the the red supergiant sample was only 3%. By placing the yellow supergiants on the H-R diagram and comparing them against the evolutionary tracks, we find that new Geneva evolutionary models do an exemplary job at predicting both the locations and the lifetimes of these transitory objects.

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.

Wolf-Rayet stars (WRs) are evolved massive stars, and the relative number of WC-type and WN-type WRs should vary with metallicity, providing a sensitive test of stellar evolutionary theory. The observed WC/WN ratio is much higher than that predicted by theory in some galaxies but this could be due to observational incompleteness for WN-types, which have weaker lines. Previous studies of M33's WR content show a galactocentric gradient in the relative numbers of WCs and WNs, but only small regions have been surveyed with sufficient sensitivity to detect all of the WNs. Here we present a sensitive survey for WRs covering all of M33, finding 55 new WRs, mostly of WN type. Our spectroscopy also improves the spectral types of many previously known WRs, establishing in one case that the star is actually a background quasar. The total number of spectroscopically confirmed WRs in M33 is 206, a number we argue is complete to approximately 5%, with most WRs residing in OB associations, although approximately 2% are truly isolated. The WC/WN ratio in the central regions (<2 kpc) of M33 is much higher than that predicted by the current Geneva evolutionary models, while the WC/WN ratios in the outer regions are in good accord, as are the values in the SMC and LMC. The WC/WN ratio and the WC subtype distribution both argue that the oxygen abundance gradient in M33 is significantly larger than found by some recent studies, but are consistent with the two-component model proposed by Magrini et al.

New absolute spectrophotometry of the Kitt Peak night sky has been obtained in 2009/10, which we compare to previously published data obtained in 1988 and 1999, allowing us to look for changes over the past two decades. A comparison of the data between 1988, 1999 and 2009/10 reveals that the sky brightness of Kitt Peak has stayed remarkably constant over the past 20 years. Compared to 1988, the 2009/10 data show no change in the sky brightness at Zenith though, as expected, the sky glow has increased most dramatically in the direction of Tucson. Comparisons between the 1999 and 2009/10 data suggest that the sky has actually decreased in brightness compared to 10 years ago. However, the older data were both taken during times of increased solar activity. When we correct the measurements for the solar irradiance fluctuations, we find that compared to 20 years ago, the sky is ~0.1 magnitude brighter at Zenith and ~0.3 magnitudes brighter towards Tucson. But even after these corrections, we still find that the sky over Kitt Peak is comparable to what it was 10 years ago at Zenith and ~0.1 magnitude darker towards Tucson. This suggests that the strengthened lighting ordinances Tucson and Pima County established in the early 2000s have been quite effective. With some care, the Kitt Peak night sky will remain this dark many years into the future.

The yellow supergiant content of nearby galaxies provides a critical test of massive star evolutionary theory. While these stars are the brightest in a galaxy, they are difficult to identify because a large number of foreground Milky Way stars have similar colors and magnitudes. We previously conducted a census of yellow supergiants within M31 and found that the evolutionary tracks predict a yellow supergiant duration an order of magnitude longer than we observed. Here we turn our attention to the SMC, where the metallicity is 10x lower than that of M31, which is important as metallicity strongly affects massive star evolution. The SMC's large radial velocity (~160 km/s) allows us to separate members from foreground stars. Observations of ~500 candidates yielded 176 near-certain SMC supergiants, 16 possible SMC supergiants, along with 306 foreground stars and provide good relative numbers of yellow supergiants down to 12Mo. Of the 176 near-certain SMC supergiants, the kinematics predicted by the Besancon model of the Milky Way suggest a foreground contamination of >4%. After placing the SMC supergiants on the H-R diagram and comparing our results to the Geneva evolutionary tracks, we find results similar to those of the M31 study: while the locations of the stars on the H-R diagram match the locations of evolutionary tracks well, the models over-predict the yellow supergiant lifetime by a factor of ten. Uncertainties about the mass-loss rates on the main-sequence thus cannot be the primary problem with the models.