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.

We present optical spectrophotometry of the red supergiant Betelgeuse from 2020 February 15, during its recent unprecedented dimming episode. By comparing this spectrum to stellar atmosphere models for cool supergiants, as well as spectrophotometry of other Milky Way red supergiants, we conclude that Betelgeuse has a current effective temperature of 3600 +/- 25 K. While this is slightly cooler than previous measurements taken prior to Betelgeuse's recent lightcurve evolution, this drop in effective temperature is insufficient to explain Betelgeuse's recent optical dimming. We propose that episodic mass loss and an increase in the amount of large-grain circumstellar dust along our sightline to Betelgeuse is the most likely explanation for its recent photometric evolution.

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 present the detection of the first candidate colliding-wind binary (CWB) in M33, located in the giant H II region NGC 604. The source was first identified in archival {\it Chandra} imaging as a relatively soft X-ray point source, with the likely primary star determined from precise astrometric alignment between archival {\it Hubble Space Telescope} and {\it Chandra} imaging. The candidate primary star in the CWB is classified for the first time in this work as a carbon-rich Wolf-Rayet star with a likely O star companion based on spectroscopy obtained from Gemini-North. We model the X-ray spectrum using {\it Chandra} and {\it XMM-Newton} observations, and find the CWB is well-fit as a $\sim$ 1 keV thermal plasma with a median unabsorbed luminosity in the 0.5--2.0 keV band of $L_{\rm X}$ $\sim$ 3 $\times$ 10$^{35}$ erg s$^{-1}$, making this source among the brightest of CWBs observed to date. We present a long term light curve for the candidate CWB from archival {\it Chandra} and {\it XMM-Newton} observations, and discuss the constraints placed on the binary by this light curve, as well as the X-ray luminosity at maximum. Finally, we compare this candidate CWB in M33 to other well-studied, bright CWBs in the Galaxy and Magellanic Clouds, such as $\eta$ Car.

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 investigate the red supergiant (RSG) population of M31, obtaining radial velocities of 255 stars. These data substantiate membership of our photometrically-selected sample, demonstrating that Galactic foreground stars and extragalactic RSGs can be distinguished on the basis of B-V, V-R two-color diagrams. In addition, we use these spectra to measure effective temperatures and assign spectral types, deriving physical properties for 192 RSGs. Comparison with the solar-metallicity Geneva evolutionary tracks indicates astonishingly good agreement. The most luminous RSGs in M31 are likely evolved from 25-30 Mo stars, while the vast majority evolved from stars with initial masses of 20 Mo or less. There is an interesting bifurcation in the distribution of RSGs with effective temperatures that increases with higher luminosities, with one sequence consisting of early K-type supergiants, and with the other consisting of M-type supergiants that become later (cooler) with increasing luminosities. This separation is only partially reflected in the evolutionary tracks, although that might be due to the mis-match in metallicities between the solar Geneva models and the higher-than-solar metallicity of M31. As the luminosities increase the median spectral type also increases; i.e., the higher mass RSGs spend more time at cooler temperatures than do those of lower luminosities, a result which is new to this study. Finally we discuss what would be needed observationally to successfully build a luminosity function that could be used to constrain the mass-loss rates of RSGs as our Geneva colleagues have suggested.

We investigate Wolf-Rayet (WR) stars as a source of feedback contributing to the removal of natal material in the early evolution of massive star clusters. Despite previous work suggesting that massive star clusters clear out their natal material before the massive stars evolve into the WR phase, WR stars have been detected in several emerging massive star clusters. These detections suggest that the timescale for clusters to emerge can be at least as long as the time required to produce WR stars (a few million years), and could also indicate that WR stars may be providing the tipping point in the combined feedback processes that drive a massive star cluster to emerge. We explore the potential overlap between the emerging phase and the WR phase with an observational survey to search for WR stars in emerging massive star clusters hosting WR stars. We select candidate emerging massive star clusters from known radio continuum sources with thermal emission and obtain optical spectra with the 4m Mayall Telescope at Kitt Peak National Observatory and the 6.5m MMT. We identify 21 sources with significantly detected WR signatures, which we term "emerging WR clusters." WR features are detected in $\sim$50% of the radio-selected sample, and thus we find that WR stars are commonly present in massive star clusters currently emerging. The observed extinctions and ages suggest that clusters without WR detections remain embedded for longer periods of time, and may indicate that WR stars can aid, and therefore accelerate, the emergence process.

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.

A significant percentage of OB stars are runaways, so we should expect a similar percentage of their evolved descendants to also be runaways. However, recognizing such stars presents its own set of challenges, as these older, more evolved stars will have drifted further from their birthplace, and thus their velocities might not be obviously peculiar. Several Galactic red supergiants (RSGs) have been described as likely runaways, based upon the existence of bow shocks, including Betelgeuse. Here we announce the discovery of a runaway RSG in M31, based upon a 300 km s$^{-1}$ discrepancy with M31's kinematics. The star is found about 21' (4.6 kpc) from the plane of the disk, but this separation is consistent with its velocity and likely age ($\sim$10 Myr). The star, J004330.06+405258.4, is an M2 I, with $M_V=-5.7$, $\log L/L_\odot$=4.76, an effective temperature of 3700 K, and an inferred mass of 12-15$M_\odot$. The star may be a high-mass analog of the hypervelocity stars, given that its peculiar space velocity is probably 400-450 km s$^{-1}$, comparable to the escape speed from M31's disk.

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.

The feedback from massive stars is important to super star cluster (SSC) evolution and the timescales on which it occurs. SSCs form embedded in thick material, and eventually, the cluster is cleared out and revealed at optical wavelengths -- however, this transition is not well understood. We are investigating this critical SSC evolutionary transition with a multi-wavelength observational campaign. Although previously thought to appear after the cluster has fully removed embedding natal material, we have found that SSCs may host large populations of Wolf-Rayet stars. These evolved stars provide ionization and mechanical feedback that we hypothesize is the tipping point in the combined feedback processes that drive a SSC to emerge. Utilizing optical spectra obtained with the 4m Mayall Telescope at Kitt Peak National Observatory and the 6.5m MMT, we have compiled a sample of embedded SSCs that are likely undergoing this short-lived evolutionary phase and in which we confirm the presence of Wolf-Rayet stars. Early results suggest that WRs may accelerate the cluster emergence.

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.

The post main-sequence evolution of massive stars is very sensitive to many parameters of the stellar models. Key parameters are the mixing processes, the metallicity, the mass-loss rate and the effect of a close companion. We study how the red supergiant lifetimes, the tracks in the Hertzsprung-Russel diagram (HRD), the positions in this diagram of the pre-supernova progenitor as well as the structure of the stars at that time change for various mass-loss rates during the red supergiant phase (RSG), and for two different initial rotation velocities. The surface abundances of RSGs are much more sensitive to rotation than to the mass-loss rates during that phase. A change of the RSG mass-loss rate has a strong impact on the RSG lifetimes and therefore on the luminosity function of RSGs. At solar metallicity, the enhanced mass-loss rate models do produce significant changes on the populations of blue, yellow and red supergiants. When extended blue loops or blue ward excursions are produced by enhanced mass-loss, the models predict that a majority of blue (yellow) supergiants are post RSG objects. These post RSG stars are predicted to show much smaller surface rotational velocities than similar blue supergiants on their first crossing of the HR gap. The position in the HRD of the end point of the evolution depends on the mass of the hydrogen envelope. More precisely, whenever, at the pre-supernova stage, the H-rich envelope contains more than about 5\% of the initial mass, the star is a red supergiant, and whenever the H-rich envelope contains less than 1\% of the total mass the star is a blue supergiant. For intermediate situations, intermediate colors/effective temperatures are obtained. Yellow progenitors for core collapse supernovae can be explained by the enhanced mass-loss rate models, while the red progenitors are better fitted by the standard mass-loss rate models.

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.

Thorne-Zytkow objects (TZOs) are a theoretical class of star in which a compact neutron star is surrounded by a large, diffuse envelope. Supergiant TZOs are predicted to be almost identical in appearance to red supergiants (RSGs). The best features that can be used at present to distinguish TZOs from the general RSG population are the unusually strong heavy-element and Li lines present in their spectra, products of the star's fully convective envelope linking the photosphere with the extraordinarily hot burning region in the vicinity of the neutron star core. Here we present our discovery of a TZO candidate in the Small Magellanic Cloud. It is the first star to display the distinctive chemical profile of anomalous element enhancements thought to be unique to TZOs. The positive detection of a TZO will provide the first direct evidence for a completely new model of stellar interiors, a theoretically predicted fate for massive binary systems, and never-before-seen nucleosynthesis processes that would offer a new channel for Li and heavy-element production in our universe.

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.

We present results of a high angular resolution survey of massive OB stars in the Cygnus OB2 association that we conducted with the Fine Guidance Sensor 1R (FGS1r) on the Hubble Space Telescope. FGS1r is able to resolve binary systems with a magnitude difference delta-V < 4 down to separations as small as 0.01 arcsec. The sample includes 58 of the brighter members of Cyg OB2, one of the closest examples of an environment containing a large number of very young and massive stars. We resolved binary companions for 12 targets and confirmed the triple nature of one other target, and we offer evidence of marginally resolved companions for two additional stars. We confirm the binary nature of 11 of these systems from complementary adaptive optics imaging observations. The overall binary frequency in our study is 22% to 26% corresponding to orbital periods ranging from 20 - 20,000 years. When combined with the known short-period spectroscopic binaries, the results supports the hypothesis that the binary fraction among massive stars is > 60%. One of the new discoveries is a companion to the hypergiant star MT 304 = Cyg OB2-12, and future measurements of orbital motion should provide mass estimates for this very luminous star.

The star-forming galaxies of the Local Group act as our laboratories for testing massive star evolutionary models. In this review, I briefly summarize what we believe we know about massive star evolution, and the connection between OB stars, Luminous Blue Variables, yellow supergiants, red supergiants, and Wolf-Rayet stars. The difficulties and recent successes in identifying these various types of massive stars in the neighboring galaxies of the Local Group will be discussed.

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.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

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.

We report the detection of oxygen-rich circumstellar envelopes in stars of the nearby (700 kpc) starburst galaxy IC 10. The star formation history and the chemical environment of this galaxy makes it an ideal target to observe dust production by high-mass stars in a low-metallicity environment. The goal of this study is to identify oxygen-rich stars in IC 10 and to constrain their nature between asymptotic giant branch stars (AGBs), red supergiants (RSGs), and other infrared bright sources. We examine the mass-loss rate of the stars and compare to results obtained for the Magellanic Clouds. Our objectives are to (1) assess whether RSGs can be significant dust producers in IC 10, and (2), solve the discrepancy between the star formation history of IC 10 and the relatively low number of RSGs detected in the optical. We search for silicate dust in emission by using the spectral map observed with the Infrared Spectrograph on board the Spitzer Space Telescope. The optical (UBVRI) and infrared (JHK, Spitzer/IRAC and Spitzer/MIPS) photometry is used to assert the membership of the stars to IC 10 and disentangle between AGBs and RSGs. Radiative models are used to infer mass-loss rates and stellar luminosities. The luminosity and colors of at least 9 silicate emission sources are consistent with stars within IC 10. Furthermore, the photometry of 2 of these sources is consistent with RSGs. We derive dust mass-loss rates similar to the values found in the Magellanic Clouds. Accounting for the sample completeness, RSGs are not important contributors to the dust mass budget in IC 10.

This is the third paper in a series aiming at the analysis of nitrogen abundances in O-type stars, to enable further constraints on the early evolution of massive stars. We provide first theoretical predictions for the NIV4058/NIII4640 emission line ratio in dependence of various parameters, and confront them with results from the analysis of a sample of early LMC/SMC O-stars. Stellar and wind parameters are determined by line profile fitting of H/He/N lines, exploiting the helium and nitrogen ionization balance. Corresponding synthetic spectra are calculated using the NLTE atmospheric code FASTWIND. Though there is a monotonic relationship between the emission line ratio and Teff, all other parameters being equal, theoretical predictions indicate additional dependencies, most notably, on the nitrogen abundance. These basic predictions are confirmed by results from atmospheric code CMFGEN. The effective temperatures for the earliest O-stars, inferred from the nitrogen ionization balance, are partly considerably hotter than indicated by previous studies. Consistent with earlier results, effective temperatures increase from supergiants to dwarfs for all spectral types in the LMC. The relation between observed NIV4058/NIII4640 emission line ratio and Teff, for a given luminosity class, turned out to be quite monotonic for our sample stars, and fairly consistent with our model predictions. The scatter within a spectral sub-type is mainly produced by abundance effects. Our findings suggest that the Walborn et al. (2002) classification scheme is able to provide a meaningful relation between spectral type and Teff, provided that it is possible to discriminate for the luminosity class. This might be difficult to achieve in low-Z environments such as the SMC, owing to rather low wind-strengths. According to our predictions, the major bias of the classification scheme is due to nitrogen content.

We used a recent grid of stellar models computed with and without rotation to make predictions concerning the WR populations and the frequency of different types of core-collapse SNe. Current rotating models were checked to provide good fits to the following features: solar luminosity and radius at the solar age, main-sequence width, red-giant and red-supergiant (RSG) positions in the HRD, surface abundances, and rotational velocities. Rotating stellar models predict that about half of the observed WR stars and at least half of the type Ibc SNe may be produced through the single-star evolution channel. Rotation increases the duration of the WNL and WNC phases, while reducing those of the WNE and WC phases, as was already shown in previous works. Rotation increases the frequency of type Ic SNe. The upper mass limit for type II-P SNe is \sim 19.0 MSun for the non rotating models and \sim 16.8 MSun for the rotating ones. Both values agree with observations. Moreover, present rotating models provide a very good fit to the progenitor of SN 2008ax. We discuss future directions of research for further improving the agreement between the models and the observations. We conclude that the mass-loss rates in the WNL and RSG phases are probably underestimated at present. We show that up to an initial mass of 40 M\odot, a surface magnetic field inferior to about 200 G may be sufficient to produce some braking. Much lower values are needed at the red supergiant stage. We suggest that the presence/absence of any magnetic braking effect may play a key role in questions regarding rotation rates of young pulsars and the evolution leading to LGRBs.

We present moderate-resolution spectroscopic observations of red supergiants (RSGs) in the low-metallicity Local Group galaxies NGC 6822 (Z = 0.4Zsun) and Wolf-Lundmark-Melotte (WLM; Z = 0.1Zsun). By combining these observations with reduction techniques for multislit data reduction and flux calibration, we are able to analyze spectroscopic data of 16 RSGs in NGC 6822 and spectrophotometric data of 11 RSGs in WLM. Using these observations we determine spectral types for these massive stars, comparing them to Milky Way and Magellanic Clouds RSGs and thus extending observational evidence of the abundance-dependent shift of RSG spectral types to lower metallicities. In addition, we have uncovered two RSGs with unusually late spectral types (J000158.14-152332.2 in WLM, with a spectral type of M3 I, and J194453.46-144552.6 in NGC 6822, with a spectral type of M4.5 I) and a third RSG (J194449.96-144333.5 in NGC 6822) whose spectral type has varied from a M2.5 in 1997 to a K5 in 2008. All three of these stars could potentially be members of a recently-discovered class of extreme RSG variables.

Yellow and red supergiants are evolved massive stars whose numbers and locations on the HR diagram can provide a stringent test for models of massive star evolution. Previous studies have found large discrepancies between the relative number of yellow supergiants observed as a function of mass and those predicted by evolutionary models, while a disagreement between the predicted and observed locations of red supergiants on the HR diagram was only recently resolved. Here we extend these studies by examining the yellow and red supergiant populations of M33. Unfortunately, identifying these stars is difficult as this portion of the color-magnitude diagram is heavily contaminated by foreground dwarfs. We identify the red supergiants through a combination of radial velocities and a two-color surface gravity discriminant and, after re-characterizing the rotation curve of M33 with our newly selected red supergiants, we identify the yellow supergiants through a combination of radial velocities and the strength of the OI $\lambda$7774 triplet. We examine ~1300 spectra in total and identify 121 yellow supergiants (a sample which is unbiased in luminosity above log(L/L\odot) ~ 4.8) and 189 red supergiants. After placing these objects on the HR diagram, we find that the latest generation of Geneva evolutionary tracks show excellent agreement with the observed locations of our red and yellow supergiants, the observed relative number of yellow supergiants with mass and the observed red supergiant upper mass limit. These models therefore represent a drastic improvement over previous generations.

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.

To explore the properties of extreme outer stellar disks, we obtained ultra-deep V and GALEX UV images of 4 dwarf irregular galaxies and one Blue Compact Dwarf galaxy and ultra-deep B images of 3 of these. Our V-band surface photometry extends to 29.5 magnitudes arcsec^-2. We convert the FUV and V-band photometry, along with Halpha photometry, into radial star formation rate profiles that are sensitive to timescales from 10 Myrs to the lifetime of the galaxy. We also compare the stellar distributions, surface brightness profiles, and star formation rate profiles to HI-line emission maps, gas surface density profiles, and gas kinematics. Our data lead us to two general observations: First, the exponential disks in these irregular galaxies are extraordinarily regular. The stellar disks continue to decline exponentially as far as our measurements extend. In spite of lumpiness in the distribution of young stars and HI distributions and kinematics that have significant unordered motions, sporadic processes that have built the disks-star formation, radial movement of stars, and perhaps even perturbations from the outside-have, nevertheless, conspired to produce standard disk profiles. Second, there is a remarkable continuity of star formation throughout these disks over time. In four out of five of our galaxies the star formation rate in the outer disk measured from the FUV tracks that determined from the V-band, to within factors of 5, requiring star formation at a fairly steady rate over the galaxy's lifetime. Yet, the HI surface density profiles generally decline with radius more shallowly than the stellar light, and the gas is marginally gravitationally stable against collapse into clouds. Outer stellar disks are challenging our concepts of star formation and disk growth and provide a critical environment in which to understand processes that mold galaxy disks.

Spectroscopy is one of the most important tools that an astronomer has for studying the universe. This chapter begins by discussing the basics, including the different types of optical spectrographs, with extension to the ultraviolet and the near-infrared. Emphasis is given to the fundamentals of how spectrographs are used, and the trade-offs involved in designing an observational experiment. It then covers observing and reduction techniques, noting that some of the standard practices of flat-fielding often actually degrade the quality of the data rather than improve it. Although the focus is on point sources, spatially resolved spectroscopy of extended sources is also briefly discussed. Discussion of differential extinction, the impact of crowding, multi-object techniques, optimal extractions, flat-fielding considerations, and determining radial velocities and velocity dispersions provide the spectroscopist with the fundamentals needed to obtain the best data. Finally the chapter combines the previous material by providing some examples of real-life observing experiences with several typical instruments.

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.

We discuss, in the context of the single star scenario, the nature of the progenitors of Red Supergiants (RSG), of Luminous Blue Variables (LBV) and of Wolf-Rayet (WR) stars. These three different populations correspond to evolved phases of Main Sequence (MS) OB stars. Axial rotation and mass loss have a great influence on massive star evolution in general and more specifically on the durations of these different phases. Moderate rotation and mass loss, during the MS phase, favor the evolution towards the RSG stage. Fast rotation and strong mass loss during the MS phase, in contrast, prevent the star from becoming a RSG and allow the star to pass directly from the OB star phase into the WR phase. Mass loss during the RSG stage may make the star evolve back in the blue part of the HR diagram. We argue that such an evolution may be more common than presently accounted for in stellar models. This might be the reason for the lack of type IIP SNe with RSG progenitors having initial masses between 18 and 30 M$_\odot$. The LBVs do appear as a possible transition phase between O and WR stars or between WNL and WNE stars. Fast rotation and/or strong mass loss during the Main-Sequence phase prevent the formation of LBV stars. The mechanisms driving the very strong ejections shown by LBV stars are still unknown. (abridged)

We have carried out a narrow-band survey of the Local Group galaxy, M33, in the HeII4686 emission line, to identify HeII nebulae in this galaxy. With spectroscopic follow-up observations, we confirm three of seven candidate objects, including identification of two new HeII nebulae, BCLMP651, HBW673. We also obtain spectra of associated ionizing stars for all the HII regions, identifying two new WN stars. We demonstrate that the ionizing source for the known HeII nebula, MA 1, is consistent with being the early-type WN star MC8 (M33-WR14), by carrying out a combined stellar and nebular analysis of MC8 and MA1. We were unable to identify the helium ionizing sources for HBW 673 and BCLMP 651, which do not appear to be Wolf-Rayet stars. According to the [OIII]5007/Hbeta vs [NII]6584/Halpha diagnostic diagram, excitation mechanisms apart from hot stellar continuum are needed to account for the nebular emission in HBW 673, which appears to have no stellar source at all.

We have obtained UBVRI images with the Kitt Peak and Cerro Tololo 4-m telescopes and Mosaic cameras of seven dwarfs in (or near) the Local Group, all of which have known evidence of recent star formation: IC10, NGC 6822, WLM, Sextans B, Sextans A, Pegasus,and Phoenix. We construct color-magnitude diagrams (CMDs) of these systems, as well as neighboring regions that can be used to evaluate the degree of foreground contamination by stars in the Milky Way. Inter-comparison of these CMDs with those of M31, M33, the LMC, and the SMC permits us to determine improved reddening values for a typical OB star found within these galaxies. All of the CMDs reveal a strong or modest number of blue supergiants. All but Pegasus and Phoenix also show the clear presence of red supergiants in the CMD, although IC10 appears to be deficient in these objects given its large WR population. The bright stars of intermediate color in the CMD are badly contaminated by foreground stars (30-100%), and considerable spectroscopy is needed before statistics on the yellow supergiants in these systems will be known. This study is intended to serve both as the impetus and "finding charts" for further space-based imaging, and for many spectroscopic programs at large aperture.

There are various ways of measuring the slope of the upper end of the IMF. Arguably the most direct of these is to place stars on the H-R diagram and compare their positions with stellar evolutionary models. Even so, the masses one infers from this depend upon the exact methodology used. I briefly discusssome of the caveats and go through a brief error analysis. I conclude that the current data suggest that the IMF slopes are the same to within the errors. Similarly the determination of the upper mass "limit" is dependent upon how well one can determine the masses of the most massive stars within a cluster. The recent finding by Crowther et al (2010) invalidates the claim that there is a 150Mo upper limit to the IMF, but this is really not surprising given the weakness of the previous evidence.

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.

When we look at a nearby galaxy, we see a mixture of foreground stars and bona fide extragalactic stars. I will describe what we need to do to get meaningful statistics on the massive star populations across the H-R diagram. Such a census provides the means of a very powerful test of massive star evolutionary theory.

The yellow supergiant content of nearby galaxies can provide a critical test of stellar evolution theory, bridging the gap between the hot, massive stars and the cool red supergiants. But, this region of the color-magnitude diagram is dominated by foreground contamination, requiring membership to somehow be determined. Fortunately, the large negative systemic velocity of M31, coupled to its high rotation rate, provides the means for separating the contaminating foreground dwarfs from the bona fide yellow supergiants within M31. Using the MMT, we obtained spectra of about 2900 stars, selected using the color and magnitude range to be yellow supergiants. Comparing the velocities to that of M31's rotation curve, we identified 54 certain, and 66 probable yellow supergiants from among the sea of foreground dwarfs. We find excellent agreement between the location of yellow supergiants in the H-R diagram and that predicted by the latest Geneva evolutionary tracks which include rotation. However, the relative number of yellow supergiants seen as a function of mass varies from that predicted by the models by a factor of more than 10, in the sense that more high mass yellow supergiants are predicted than are actually observed. Comparing the total number of yellow supergiants with masses above 20Mo with the estimated number of unevolved O stars indicates that the duration of the yellow supergiant phase is about 3000 years. This is consistent with what the 12Mo and 15Mo evolutionary tracks predict, but disagrees with the 20,000-80,000 year time scales predicted by the models for higher masses.

Red supergiants are a short-lived stage in the evolution of moderately massive stars (10-25Mo), and as such their location in the H-R diagram provides an exacting test of stellar evolutionary models. Since massive star evolution is strongly affected by the amount of mass-loss a star suffers, and since the mass-loss rates depend upon metallicity, it is highly desirable to study the physical properties of these stars in galaxies of various metallicities. Here we identify a sample of red supergiants in M31 (the most metal-rich of the Local Group galaxies) and derive their physical properties by fitting MARCS atmosphere models to moderate resolution optical spectroscopy, and from V-K photometry.

WOH G64 is an unusual red supergiant (RSG) in the Large Magellanic Cloud (LMC), with a number of properties that set it apart from the rest of the LMC RSG population, including a thick circumstellar dust torus, an unusually late spectral type, maser activity, and nebular emission lines. Its reported physical properties are also extreme, including the largest radius for any star known and an effective temperature that is much cooler than other RSGs in the LMC, both of which are at variance with stellar evolutionary theory. We fit moderate-resolution optical spectrophotometry of WOH G64 with the MARCS stellar atmosphere models, determining an effective temperature of 3400 +/- 25 K. We obtain a similar result from the star's broadband V - K colors. With this effective temperature, and taking into account the flux contribution from the aysmmetric circumstellar dust envelope, we calculate log(L/L_sun) = 5.45 +/- 0.05 for WOH G64, quite similar to the luminosity reported by Ohnaka and collaborators based on their radiative transfer modeling of the star's dust torus. We determine a radius of R/R_sun = 1540, bringing the size of WOH G64 and its position on the H-R diagram into agreement with the largest known Galactic RSGs, although it is still extreme for the LMC. In addition, we use the Ca II triplet absorption feature to determine a radial velocity of 294 +/- 2 km/s for the star; this is the same radial velocity as the rotating gas in the LMC's disk, which confirms its membership in the LMC and precludes it from being an unusual Galactic halo giant. Finally, we describe the star's unusual nebula emission spectrum; the gas is nitrogen-rich and shock-heated, and displays a radial velocity that is significantly more positive than the star itself by 50 km/s.

Most stars are born in rich young stellar clusters (YSCs) embedded in giant molecular clouds. The most massive stars live out their short lives there, profoundly influencing their natal environments by ionizing HII regions, inflating wind-blown bubbles, and soon exploding as supernovae. Thousands of lower-mass pre-main sequence stars accompany the massive stars, and the expanding HII regions paradoxically trigger new star formation as they destroy their natal clouds. While this schematic picture is established, our understanding of the complex astrophysical processes involved in clustered star formation have only just begun to be elucidated. The technologies are challenging, requiring both high spatial resolution and wide fields at wavelengths that penetrate obscuring molecular material and remove contaminating Galactic field stars. We outline several important projects for the coming decade: the IMFs and structures of YSCs; triggered star formation around YSC; the fate of OB winds; the stellar populations of Infrared Dark Clouds; the most massive star clusters in the Galaxy; tracing star formation throughout the Galactic Disk; the Galactic Center region and YSCs in the Magellanic Clouds. Programmatic recommendations include: developing a 30m-class adaptive optics infrared telescope; support for high-resolution and wide field X-ray telescopes; large-aperture sub-millimeter and far-infrared telescopes; multi-object infrared spectrographs; and both numerical and analytical theory.

In order to better determine the physical properties of hot, massive stars as a function of metallicity, we obtained very high SNR optical spectra of 26 O and early B stars in the Magellanic Clouds. These allow accurate modeling even in cases where the He I 4471 line has an equivalent width of only a few tens of mA. The spectra were modeled with FASTWIND, with good fits obtained for 18 stars; the remainder show signatures of being binaries. We include stars in common to recent studies to investigate possible systematic differences. The "automatic" FASTWIND modeling method of Mokiem and collaborators produced temperatures 1100 K hotter on the average, presumably due to the different emphasis given to various temperature-sensitive lines. More significant, however, is that the automatic method always produced some "best" answer, even for stars we identify as composite (binaries). The temperatures found by the TLUSTY/CMFGEN modeling of Bouret, Heap, and collaborators yielded temperatures 1000 K cooler than ours, on average. Significant outliers were due either to real differences in the data (some of the Bouret/Heap data were contaminated by moonlight continua) or the fact we could detect the HeI line needed to better constrain the temperature. Our new data agrees well with the effective temperature scale we presented previously. We confirm that the "Of" emission-lines do not track luminosity classes in the exact same manner as in Milky Way stars. We revisit the the issue of the "mass discrepancy", finding that some of the stars in our sample do have spectroscopic masses that are significantly smaller than those derived from stellar evolutionary models. We do not find that the size of the mass discrepancy is simply related to either effective temperature or surface gravity.

We investigate the massive star content of NGC 3603, the closest known giant H II region. We have obtained spectra of 26 stars in the central cluster using the Baade 6.5-m telescope (Magellan I). Of these 26 stars, 16 had no previous spectroscopy. We also obtained photometry of all of the stars with previous or new spectroscopy, primarily using archival HST ACS/HRC images. We use these data to derive an improved distance to the cluster, and to construct an H-R diagram for discussing the masses and ages of the massive star content of this cluster.

Red supergiants (RSGs) are an evolved stage in the life of intermediate massive stars (than than 25 solar masses). For many years, their location in the H-R diagram was at variance with the evolutionary models. Using the MARCS stellar atmospheres, we have determined new effective temperatures and bolometric luminosities for RSGs in the Milky Way, LMC, and SMC, and our work has resulted in much better agreement with the evolutionary models. We have also found evidence of significant visual extinction due to circumstellar dust. Although in the Milky Way the RSGs contribute only a small fraction (than than 1 percent) of the dust to the interstellar medium (ISM), in starburst galaxies or galaxies at large look-back times, we expect that RSGs may be the main dust source. We are in the process of extending this work now to RSGs of higher and lower metallicities using the galaxies M31 and WLM.

We present the results of optical wavelength observations of the unusual SMC eclipsing binary system HD 5980 obtained in 1999 and 2004--2005. Radial velocity curves for the erupting LBV/WR object (star A) and its close WR-like companion (star B) are obtained by deblending the variable emission-line profiles of N IV and N V lines under the simplistic assumption that these lines originate primarily in the winds of star A and star B. The derived masses M_A=58--79 Mo and M_B=51--67 Mo, are more consistent with the stars' location near the top of the HRD than previous estimates. The presence of a wind-wind interaction region is inferred from the orbital phase-dependent behavior of He I P Cygni absorption components. The emission-line intensities continued with the declining trend previously seen in UV spectra. The behavior of the photospheric absorption lines is consistent with the results of Schweickhardt (2002) who concludes that the third object in the combined spectrum, star C, is also a binary system with P(starC)~96.5 days, e=0.83. The data used in this paper will be made publicly available for further analysis.

We describe a search for H-alpha emission-lined stars in M31, M33, and seven dwarfs in or near the Local Group (IC 10, NGC 6822, WLM, Sextans B, Sextans A, Pegasus and the Phoenix dwarf) using interference filter imaging with the KPNO and CTIO 4-m telescope and Mosaic cameras. The survey is aimed primarily at identifying new Luminous Blue Variables (LBVs) from their spectroscopic similarity to known LBVs, avoiding the bias towards photometric variability, which may require centuries to manifest itself if LBVs go through long quiescent periods. Followup spectroscopy with WIYN confirms that our survey detected a wealth of stars whose spectra are similar to the known LBVs. We "classify" the spectra of known LBVs, and compare these to the spectra of the new LBV candidates. We demonstrate spectacular spectral variability for several of the new LBV candidates, such as AM2, previously classified as a Wolf-Rayet star, which now shows FeI, FeII and Balmer emission lines but neither the NIII 4634,42 nor HeII 4686 emission that it did in 1982. Profound spectral changes are also noted for other suspected and known LBVs. Several of the LBV candidates also show >0.5 mag changes in V over the past 10-20 years. The number of known or suspected LBVs is now 24 in M31, 37 in M33, 1 in NGC 6822, and 3 in IC 10. We estimate that the total number of LBVs in M31 and M33 may be several hundred, in contrast to the 8 known historically through large-scale photometric variability. This has significant implications for the time scale of the LBV phase. We also identify a few new WRs and peculiar emission-lined objects.

Red supergiants (RSGs) are a He-burning phase in the evolution of moderately high mass stars (10-25 solar masses). The evolution of these stars, particularly at low metallicities, is still poorly understood. The latest-type RSGs in the Magellanic Clouds are cooler than the current evolutionary tracks allow, occupying the region to the right of the Hayashi limit where stars are no longer in hydrodynamic equilibrium. We have discovered four Cloud RSGs in this region that display remarkably similar unusual behavior. All of them show considerable variations in their V magnitudes and effective temperatures (and spectral types). Two of these stars, HV 11423 and [M2002] SMC 055188, have been observed in an M4.5 I state, considerably later and cooler than any other supergiant in the SMC. These stars suffer dramatic physical changes on timescales of months - when they are at their warmest, they are also brighter, more luminous, and show an increased amount of extinction. This variable extinction is characteristic of the effects of circumstellar dust, and can be connected with sporadic dust production from these stars in their cooler states. We suggest that these unusual properties are indicative of an unstable (and short-lived) evolutionary phase not previously associated with RSGs, and consider the implications such behavior could have for our understanding of the latest stages of massive star evolution in low-metallicity environments.

Red supergiants (RSGs) are an evolved stage in the life of intermediate massive stars (<25Mo). For many years their location in the H-R diagram was at variance with the evolutionary models. Using the MARCS stellar atmosphere models, we have determined new effective temperatures and bolometric luminosities for RSGs in the Milky Way, LMC, and SMC, and our work has resulted in much better agreement with the evolutionary models. We have also found evidence of significant visual extinction due to circumstellar dust. Although in the Milky Way the RSGs contribute only a small fraction (<1%) of the dust to the interstellar medium (ISM), in starburst galaxies or galaxies at large look-back times, we expect that RSGs may be the main dust source. We are in the process of extending this work now to RSGs of higher and lower metallicities using the galaxies M31 and WLM.

We have identified seven red supergiants (RSGs) in the Large Magellanic Cloud (LMC) and four RSGs in the Small Magellanic Cloud (SMC), all of which have spectral types that are considerably later than the average type observed in their parent galaxy. Using moderate-resolution optical spectrophotometry and the MARCS stellar atmosphere models, we determine their physical properties and place them on the H-R diagram for comparison with the predictions of current stellar evolutionary tracks. The radial velocities of these stars suggest that they are likely all members of the Clouds rather than foreground dwarfs or halo giants. Their locations in the H-R diagram also show us that those stars are cooler than the current evolutionary tracks allow, appearing to the right of the Hayashi limit, a region in which stars are no longer in hydrodynamic equilibrium. These stars exhibit considerable variability in their V magnitudes, and three of these stars also show changes in their effective temperatures (and spectral types) on the time-scales of months. One of these stars, [M2002] SMC 055188, was caught in an M4.5 I state, as late as that seen in HV 11423 at its recent extreme: considerable later, and cooler, than any other supergiant in the SMC. In addition, we find evidence of variable extinction due to circumstellar dust and changes in the stars' luminosities, also consistent with our recent findings for HV 11423 - when these stars are hotter they are also dustier and more luminous. We suggest that these stars have unusual properties because they are in an unstable (and short-lived) evolutionary phase.

We call attention to the fact that one of the brightest red supergiants in the SMC has recently changed its spectral type from K0-1 I (December 2004) to M4 I (December 2005) and back to K0-1 I (September 2006). An archival spectrum from the Very Large Telescope reveals that the star was even cooler (M4.5-M5 I) in December 2001. By contrast, the star was observed to be an M0 I in both October 1978 and October 1979. The M4-5 I spectral types is by far the latest type seen for an SMC supergiant, and its temperature in that state places it well beyond the Hayashi limit into a region of the H-R diagram where the star should not be in hydrostatic equilibrium. The star is variable by nearly 2 mag in V, but essentially constant in K. Our modeling of its spectral energy distribution shows that the visual extinction has varied during this time, but that the star has remained essentially constant in bolometric luminosity. We suggest that the star is currently undergoing a period of intense instability, with its effective temperature changing from 4300 K to 3300 K on the time-scale of months. It has one of the highest 12-micron fluxes of any RSG in the SMC, and we suggest that the variability at V is due primarily to changes in effective temperature, and secondly, due to changes in the local extinction due to creation and dissipation of circumstellar dust. We speculate that the star may be nearing the end of its life.

We present a combined analysis of the kinematics of the Large Magellanic Cloud through its HI gas, carbon stars, and red supergiant stars. After correcting the line-of-sight velocities for the recent accurate measurement of the LMC's space motion, we find that each kinematic tracer clearly defines a flat rotation curve with similar shape but different amplitude for each tracer: 61 km/s for the carbon stars, 80 km/s for the HI, and 107 km/s for the red supergiants. Previously identified tidal HI features are seen to harbor numerous carbon stars, with the tidally disturbed stars comprising 7-15% of the total sample. This discovery implies that we cannot depend on the carbon star sample alone to construct a reliable model of the LMC's gravitational potential. We also find red supergiants with peculiar kinematics, but their association with tidal features is unclear, and may instead be interacting with supergiant HI shell SGS4. In addition, although the local velocity dispersion of the red supergiants is small, ~8 km/s, their velocity dispersion about the carbon star rotation solution is 17 km/s, equal to the velocity dispersion of the carbon stars themselves. We thus appear to be witnessing the tidal heating of the LMC's stellar disk.

The star VY CMa is a late-type M supergiant with many peculiarities, mostly related to the intense circumstellar environment due to the star's high mass-loss rate. Claims have been made that would suggest this star is considerably more luminous (L = 5 x 10^5 Lo) and larger (R=2800 Ro) than other Galactic red supergiants (RSGs). Indeed, such a location in the H-R diagram would be well in the "Hayashi forbidden zone" where stars cannot be in hydrostatic equilibrium. These extraordinary properties, however, rest upon an assumed effective temperature of 2800-3000 K, far cooler than recent work have shown RSGs to be. To obtain a better estimate, we fit newly obtained spectrophotometry in the optical and NIR with the same MARCS models used for our recent determination of the physical properties of other RSGs; we also use $V-K$ and $V-J$ from the literature to derive an effective temperatures. We find that the star likely has a temperature of 3650 K, a luminosity L = 6 x 10^4 Lo, and a radius of 600Ro. These values are consistent with VY CMa being an ordinary evolved 15 Mo RSG, and agree well with the Geneva evolutionary tracks. We find that the circumstellar dust region has a temperature of 760 K, and an effective radius of approximately 130 AU, if spherical geometry is assumed for the latter. What causes this star to have such a high mass-loss, and large variations in brightness (but with little change in color), remains a mystery at present, although we speculate that perhaps this star (and NML Cyg) are simply normal RSGs caught during an unusually unstable time.

We discuss H and K observations of three fields in the bulge and disk of M31 obtained with the Altair adaptive optics system and NIRI instrument on Gemini North. These are the highest resolution and deepest near-infrared observations obtained to date of the inner regions of M31, and demonstrate the promise of ground-based adaptive optics for studying the crowded regions of nearby galaxies. We have combined our observations with previously published HST/NICMOS observations of nine M31 fields and have derived the coarse star formation histories of M31's bulge and inner disk. From fits to the M_K luminosity functions, we find the stellar population mix to be dominated by old, nearly solar-metallicity stars. The old populations, which we define as having age >~6 Gyr, indeed dominate the star formation histories at all radii independent of the relative contributions of bulge and disk stars. Although all of our fields contain some bulge contribution, our results suggest that there is no age difference between the bulge and disk to the limit of our precision.

We use the MARCS stellar atmosphere to derive the physical properties of 36 red supergiants (RSGs) in the LMC, and 39 RSGs in the SMC using moderate-resolution optical spectrophotometry (4000-9000A) and broad-band colors (V-R, V-K). The results from the dereddened V-R colors are in good agreement with those derived from the spectrophotometry, but the dereddened V-K colors give temperatures that are 3-4% warmer for the SMC data, with the LMC and Milky Way showing a smaller but similar effect. We conclude that this discrepancy is due to the limitations of 1D models. Our newly derived effective temperatures and bolometric luminosities bring the Magellanic Cloud RSGs into good agreement with stellar evolutionary models that include the effects of rotation. A typical M2~I in the SMC is about 150 K cooler than its Galactic counterpart; one in the LMC is about 50 K cooler. This is in the sense expected due to the lower chemical abundances in the SMC and LMC, although it is not sufficient to explain the shift in average RSG spectral type seen between the SMC, LMC, and Milky Way. Instead, that is due primarily to the change in Hayashi limit with metallicity, as first proposed by Elias et al. (1985). Finally, our study confirms that many RSGs in the Magellanic Clouds are significantly more reddened than OB stars, consistent with our recent findings for Galactic stars that circumstellar dust may contribute several magnitudes of extra visual extinction.

We present UBVRI photometry obtained from Mosaic images of M31 and M33 using the KPNO 4-m telescope. The survey covers 2.2 sq degrees of M31, and 0.8 sq degrees of M33, chosen so as to include all of the regions currently active in forming massive stars. The catalog contains 371,781 and 146,622 stars in M31 and M33, respectively, where every star has a counterpart (at least) in B, V, and R. We compare our photometry to previous studies. We provide cross references to the stars confirmed as members by spectroscopy, and compare the location of these to the complete set in color-magnitude diagrams. While follow-up spectroscopy is needed for many projects, we demonstrate the success of our photometry in being able to distinguish M31/M33 members from foreground Galactic stars. We also present the results of newly obtained spectroscopy, which identifies 34 newly confirmed members, including B-A supergiants, the earliest O star known in M31, and two new Luminous Blue Variable candidates whose spectra are similar to that of P Cygni.

We present spectroscopy and discuss the photometric history of a previously obscure star in M31. The spectrum of the star is an extremely close match to that of P Cygni, one of the archetypes of Luminous Blue Variables (LBVs). The star has not shown much variability over the past 40 years ($<0.2$ mag), although small-scale (0.05 mag) variations over a year appear to be real. Nevertheless, the presence of a sub-arcsecond extension around the star is indicative of a past outburst, and from the nebula's size (0.5 pc diameter) we estimate the outburst took place roughly 2000 yrs ago. P Cygni itself exhibits a similar photometric behavior, and has a similar nebula (0.2 pc diameter). We argue that this may be more typical behavior for LBVs than commonly assumed. The star's location in the HR diagram offers substantial support for stellar evolutionary models that include the effects of rotation, as the star is just at a juncture in the evolutionary track of a 85Mo star. The star is likely in a transition from an O star to a late-type WN Wolf-Rayet.

We present new spectroscopic observations of 13 H II regions in the Local Group spiral galaxy M33. The regions observed range from 1 to 7 kpc in distance from the nucleus. Of the 13 H II regions observed, the [O III] 4363 Angstrom line was detected in six regions. Electron temperatures were thus able to be determined directly from the spectra using the [O III] 4959,5007 A/4363 A line ratio. Based on these temperature measurements, oxygen and neon abundances and their radial gradients were calculated. For neon, a gradient of -0.016 +/- 0.017 dex/kpc was computed, which agrees with the Ne/H gradient derived previously from ISO spectra. A gradient of -0.012 +/- 0.011 dex/kpc was computed for O/H, much shallower than was derived in previous studies. The newly calculated O/H and Ne/H gradients are in much better agreement with each other, as expected from predictions of stellar nucleosynthesis. We examine the correlation between the WC/WN ratio and metallicity, and find that the new M33 abundances do not impact the observed correlation significantly. We also identify two new He II-emitting H II regions in M33, the first to be discovered in a spiral galaxy other than the Milky Way. In both cases the nebular He II emission is not associated with Wolf-Rayet stars. Therefore, caution is warranted in interpreting the relationship between nebular He II emission and Wolf-Rayet stars when both are observed in the integrated spectrum of an H II region.

Deriving the physical properties of red supergiants (RSGs) depends upon accurate corrections for reddening by dust. We use our recent modeling of the optical spectra of RSGs to address this topic. We find: (1) Previous broad-band studies have underestimated the correction for extinction in the visible, and hence the luminousities, if derived from V. (2) A significant fraction of RSGs in Galactic OB associations and clusters show up to several magnitudes of excess visual extinction compared to OB stars in the same regions; we argue that this is likely due to circumstellar dust around the RSGs. (3) RSGs contribute dust grains at the rate of $3 \times 10^{-8} M_\odot$ yr$^{-1}$ kpc$^{-2}$ in the solar neighborhood, comparable to what we estimate for late-type WC Wolf-Rayet stars, $1 \times 10^{-7} M_\odot$ yr$^{-1}$ kpc$^{-2}$. In the solar neighborhood this represents only a few percent of the dust production (which is dominated by low-mass AGBs), but we note that in low-metallicity starbursts, dust production by RSGs would likely dominate over other sources.

We analyze the optical and UV spectra of an additional sample of 20 Magellanic Cloud O stars, and draw conclusions from the complete sample of 40 stars. We find (1) The SMC O3-7 dwarfs are about 4000 K cooler than their Galactic counterparts; this is in the sense expected from the decreased signficiance of line-blanketing and wind-blanketing at lower metallicities. The difference decreases with later types, becoming negligible by B0 V. A similar difference is found for the supergiants. (2) The wind momentum of these stars scales with luminosity and metallicty in the ways predicted by radiatively-driven wind theory. (3) A comparison of the masses derived from spectroscopy with those derived from stellar evolutionary theory shows a significant discrepancy for stars hotter than 45000, although good agreement is found for cooler stars. (4) For the hottest O stars (O2-3.5) neither the NIII/NIV ratio, nor even the HeI/HeII ratio, does a good job of predicting the effective temperature by itself. Instead, a full analysis is needed to derive physical parameters. Thus there are O3.5V stars which are as hot or hotter than stars classified as O2V. (5) The two stars with the most discordant radial velocities in our sample happen to be O3 "field stars". This provides the first compelling observational evidence that the "field" O stars in the Magellanic Clouds may be runaway OB stars, ejected from their birth place.

We use moderate-resolution optical spectrophotometry and the new MARCS stellar atmosphere models to determine the effective temperatures of 74 Galactic red supergiants. From these we find a new effective temperature scale that is significantly warmer than those in the literature. We show that this temperature scale, along with the newly derived bolometric corrections, gives much better agreement between our red supergiants and stellar evolutionary tracks. This agreement provides an independent verification of our new temperature scale. The combination of effective temperature and bolometric luminosities allows us to calculate stellar radii; the coolest and most luminous stars have radii of roughly 1500 solar radii (7 AU), in excellent accordance with the largest stellar radii predicted from current evolutionary theory. We find that similar results are obtained for the effective temperatures and bolometric luminosities using only the de-reddened V-K colors, providing a powerful demonstration of the self-consistency of the MARCS models.

HIZSS3 is an H I source in the Zone of Avoidance. Its radio characteristics are consistent with it being nearby (~ 1.8 Mpc) low-mass dwarf irregular (dIm) galaxy. It contains a modest H II region, but its stellar population has not been resolved by optical observations. MMT spectra of the H II region are used to derive the line-of-sight extinction (E(B-V) = 1.41 +/- 0.04) and gas metallicity (log O/H + 12 ~ 7.8). ESO VLT near-IR images clearly reveal a resolved stellar population. Narrow-band Pbeta images of the H II region are used in combination with previously published Halpha data to obtain an independent line-of-sight extinction estimate: E(B-V) = 1.32 +/- 0.04. The adopted foreground extinction is E(B-V) = 1.36 +/- 0.06. Based on the K-band luminosity function and K,J-K color-magnitude diagram, the apparent magnitude and color of the tip of the red-giant branch (TRGB) are derived. In turn, these parameters are combined with the adopted foreground extinction to estimate the distance (1.69 +/- 0.07 Mpc) and mean red giant branch metallicity ([Fe/H] = -0.5 +/- 0.1). As an ensemble, these new observations significantly strengthen the conclusion that HIZSS3 is a newly discovered low-mass dIm lurking behind the Milky Way in the outskirts of the Local Group.

We report on very soft X-ray sources (VSSs) in M31. In a survey which was most sensitive to soft sources in four 8'x8' regions covered by Chandra's ACIS-S S3 CCD, we find 33 VSSs that appear to belong to M31. Fifteen VSSs have spectral characteristics mirroring the supersoft X-ray sources studied in the Magellanic Cloud and Milky Way ($k T_{eff} \leq 100$ eV); we therefore call these ``classical'' supersoft sources, or simply supersoft sources (SSSs). Eighteen VSSs may either have small (< 10%) hard components, or slightly higher effective temperatures (but still < 350 eV). We refer to these VSSs as quasisoft sources (QSSs). While hot white dwarf models may apply to SSSs, the effective temperatures of QSSs are too high, unless, e.g., the radiation emanates from only a small portion of surface. Two of the SSSs were first detected and identified as such through ROSAT observations. One SSS and one QSS may be identified with symbiotics, and 2 SSSs with supernova remnants. Both SSSs and QSSs in the disk are found near star-forming regions, possibly indicating that they are young. VSSs in the outer disk and halo are likely to be old systems; in these regions, there are more QSSs than SSSs, which is opposite to what is found in fields closer to the galaxy center. The largest density of bright VSSs is in the bulge; some of the bulge sources are close enough to the nucleus to be remnants of the tidal disruption of a giant by the massive central black hole. By using Chandra data in combination with ROSAT and XMM observations, we find most VSSs to be highly variable, fading from or brightening toward detectability on time scales of months. There is evidence for VSSs with low luminosities ($\sim 10^{36}$ erg s$^{-1}$).

From a survey of the 3400 A region in the earliest O-type spectra, we have found that two of the four O2 giants observed in the LMC have the OIV lines there stronger than the NIV, while the other two have the opposite. A SMC counterpart also has the NIV stronger than OIV. Inspection of the blue spectra of these stars shows that the former pair have weaker N lines in all ionization states (III, IV, V) present as well as lines of CIV lambda 4658, while the latter three have stronger N lines and greater He/H. Space ultraviolet observations of two of the N-strong stars show NV wind profiles substantially stronger than those of CIV, while in the N-weak stars the CIV features are equal to or stronger than the NV. The N-strong stars are now reclassified as ON2 III(f*), newly defining that category. These characteristics strongly suggest a larger fraction of processed material in the atmospheres of the ON2 stars, which we confirm by modeling the optical spectra. In the context of current models, it is in turn implied that the ON2 stars are in a more advanced evolutionary state than the others, and/or that they had higher initial rotational velocities. The recent formulation of the effects of rotation on massive stellar evolution introduces an additional fundamental parameter, which the CNO abundances are in principle able to constrain. We present some illustrative comparisons with current Geneva evolutionary models for rotating massive stars. It is possible that these very hot, nitrogen-rich objects are products of homogeneous evolution. Our results will provide motivation for further physical modeling of the atmospheres and evolutionary histories of the most massive hot stars.

We have obtained HST and ground-based observations of a sample of 20 O-type stars in the LMC and SMC, including six of the hottest massive stars known (subtypes O2-3) in the R136 cluster. In general, these data include (a) the HST UV spectra in order to measure the terminal velocities of the stellar winds, (b) high signal-to-noise, blue-optical data where the primary temperature- and gravity-sensitive photospheric lines are found, and (c) nebular-free H-alpha profiles, which provide the mass-loss rates. The line-blanketed non-LTE atmosphere code FASTWIND was then used to determine the physical parameters of this sample of stars. We find good agreement between the synthetic line profiles for the hydrogen, He I, and He II lines in the majority of the stars we analyzed; the three exceptions show evidence of being incipiently resolved spectroscopic binaries or otherwise spectral composites. One such system is apparently an O3 V+O3 V eclipsing binary, and a follow-up radial velocity study is planned to obtain Keplerian masses. Although we did not use them to constrain the fits, good agreement is also found for the He I $\lambda 3187$ and He II $\lambda 3203$ lines in the near-UV, which we plan to exploit in future studies. Our effective temperatures are compared to those recently obtained by Repolust, Puls & Herrero for a sample of Galactic stars using the same techniques. We find that the Magellanic Cloud sample is 3,000-4,000$^\circ$K hotter than their Galactic counterparts for the early through mid-O's. These higher temperatures are the consequence of a decreased importance of wind emission, wind blanketing, and metal-line blanketing at lower metallicities.

We present HST/STIS Ultraviolet spectra of the six known Ofpe/WN9 stars (``slash stars'') in M33. These stars were selected for showing the characteristics of the Ofpe/WN9 class from previous optical ground-based spectroscopy. The UV spectra are rich in wind lines, whose strength and terminal velocity vary greatly among our target sample. We analyse the STIS spectra with non-LTE, line blanketed, spherical models with hydrodynamics, computed with the WM-basic code. We find C to be underabundant and N overabundant, respect to the solar values, with a ratio (by mass) of C/N between 0.02 to 0.9 across the sample. Some stars show very conspicuous wind lines (P Cygni profiles), while two stars have extremely weak winds. The mass-loss rates thus vary greatly across the sample. The mass-loss rates of the hottest stars are lower than typical values of WNL stars, but higher than expected for normal population I massive stars. There is indication that the mass-loss rates may be variable in time. The C/N ratio, and the other physical parameters derived by the spectral modeling (Teff, Lbol, mass), are consistent with evolutionary calculations for objects with moderately high initial masses (approx 30-50 solar masses), evolving towards the WNL stage through an enhanced mass-loss phase.

We investigate the red supergiant (RSG) content of the SMC and LMC using multi-object spectroscopy on a sample of red stars previously identified by {\it BVR} CCD photometry. We obtained high accuracy ($<1$ km s$^{-1}$) radial velocities for 118 red stars seen towards the SMC and 167 red stars seen towards the LMC, confirming most of these (89% and 95%, respectively) as red supergiants (RSGs). Spectral types were also determined for most of these RSGs. We find that the distribution of spectral types is skewed towards earlier type at lower metallicities: the average (median) spectral type is K5-7 I in the SMC, M1 I in the LMC, and M2 I in the Milky Way. We argue that RSGs in the Magellanic Clouds are 100deg (LMC) and 300deg (SMC) cooler than Galactic RSGs of the same spectral type. We compare the distribution of RSGs in the H-R diagram to that of various stellar evolutionary models; we find that none of the models produce RSGs as cool and luminous as what is actually observed. In all of our H-R diagrams, however, there is an elegant sequence of decreasing effective temperatures with increasing luminosities; explaining this will be an important test of future stellar evolutionary models.

We report the discovery of a relatively faint (V=15.5) early-type WN star in the SMC. The line strength and width of He II lambda 4686 emission is similar to that of the other SMC WNs, and the presense of N V lambda 4603,19 emission (coupled with the lack of N III) suggests this star is of spectral type WN3-4.5, and thus is similar in type to the other SMC WRs. Also like the other SMC WN stars, an early-type absorption spectrum is weakly present. The absolute magnitude is comparable to that of other (single) Galactic early-type WNs. The star is located in the Hodge 53 OB association, which is also the home of two other SMC WNs. This star, which we designate SMC-WR12, was actually detected at a high significance level in an earlier interference-filter survey, but the wrong star was observed as part of a spectroscopic followup, and this case of mistaken identity resulted in its Wolf-Rayet nature not being recognized until now.

We have obtained VLA and optical follow-up observations of the low-velocity H I source HIZSS 3 discovered by Henning et al. (2000) and Rivers (2000) in a survey for nearby galaxies hidden by the disk of the Milky Way. Its radio characteristics are consistent with this being a nearby (~1.8 Mpc) low-mass dwarf irregular galaxy (dIm). Our optical imaging failed to reveal a resolved stellar population, but did detect an extended H-alpha emission region. The location of the H-alpha source is coincident with a partially-resolved H I cloud in the 21-cm map. Spectroscopy confirms that the H$\alpha$ source has a similar radial velocity to that of the H I emission at this location, and thus we have identified an optical counterpart. The H$\alpha$ emission (100 pc in diameter and with a luminosity of $1.4\times10^{38}$ ergs s$^{-1}$) is characteristic of a single H II region containing a modest population of OB stars. The galaxy's radial velocity and distance from the solar apex suggests that it is not a Local Group member, although a more accurate distance is needed to be certain. The properties of HIZSS 3 are comparable to those of GR 8, a nearby dIm with a modest amount of current star formation. Further observations are needed to characterize its stellar population, determine the chemical abundances, and obtain a more reliable distance estimate.

The galaxies of the Local Group that are currently forming stars can serve as our laboratories for understanding star formation and the evolution of massive stars. In this talk I will summarize what I think we've learned about these topics over the past few decades of research, and briefly mention what I think needs to happen next.

IC10 is the nearest starburst galaxy, as revealed both by its Halpha surface brightness and the large number of Wolf-Rayet stars (WRs) per unit area. The relative number of known WC- to WN-type WRs has been thought to be unusually high (~2), unexpected for IC10's metallicity. In this Letter we report the first results of a new and deeper survey for WRs in IC10. We sucessfully detected all of the spectroscopically known WRs, and based upon comparisons with a neighboring control field, estimate that the total number of WRs in IC10 is about 100. We present spectroscopic confirmation of two of our WR candidates, both of which are of WN type. Our photometric survey predicts that the actual WC/WN ratio is ~0.3. This makes the WC/WN ratio of IC 10 consistent with that expected for its metallicity, but greatly increases the already unusually high number of WRs, resulting in a surface density that is about 20 times higher than in the LMC. If the majority of these candidates are spectroscopically confirmed, IC10 must have an exceptional population of high mass stars.

The h and Chi Per "double cluster" is examined using wide-field (0.98 deg x 0.98 deg) CCD UBV imaging supplemented by optical spectra of several hundred of the brightest stars. Restricting our analysis to near the cluster nuclei, we find identical reddenings (E(B-V)=0.56+/-0.01), distance moduli (11.85+/-0.05), and ages (12.8+/-1.0 Myr) for the two clusters. In addition, we find an IMF slope for each of the cluster nuclei that is quite normal for high-mass stars, Gamma=-1.3+/-0.2, indistinguishable from a Salpeter value. We derive masses of 3700 M_Sun (h) and 2800 M_Sun (Chi) integrating the PDMF from 1 to 120 M_Sun. There is evidence of mild mass segregation within the cluster cores. Our data are consistent with the stars having formed at a single epoch; claims to the contrary are very likely due to the inclusion of the substantial population of early-type stars located at similar distances in the Perseus spiral arm, in addition to contamination by G and K giants at various distances. We discuss the uniqueness of the double cluster, citing other examples of such structures in the literature, but concluding that the nearly identical nature of the two cluster cores is unusual. We fail to settle the long-standing controversy regarding whether or not the double cluster is the core of the Per OB1 association, and argue that this may be unanswerable with current techniques. We also emphasize the need for further work on the pre-main sequence population of this nearby and highly interesting region.

We present a study of optical spectra of the Wolf--Rayet star AzV 336a (= SMC WR7) in the Small Magellanic Cloud. Our study is based on data obtained at several Observatories between 1988 and 2001. We find SMC WR7 to be a double lined WN+O6 spectroscopic binary with an orbital period of 19.56 days. The radial velocities of the He absorption lines of the O6 component and the strong He{\sc ii} emission at $\lambda$4686\AA of the WN component describe antiphased orbital motions. However, they show a small phase shift of $\sim$ 1 day. We discuss possible explanations for this phase shift. The amplitude of the radial velocity variations of He {\sc ii} emission is twice that of the absorption lines. The binary components have fairly high minimum masses, $\sim$ 18 $\modot$ and 34 $\modot$ for the WN and O6 components, respectively.

We present our catalog of U, B, V, and I stellar photometry of the central 18 sq. degree area of the Small Magellanic Cloud. We combine our data with the 2MASS and DENIS catalogs to provide, when available, U through K_S data for stars. Internal and external astrometric and photometric tests using existing optical photometry (U, B, and V from Massey's bright star catalog; B, V, and I from the microlensing database of OGLE and I from the near-infrared sky survey DENIS) are used to determine the observational uncertainties and identify systematic errors. We fit stellar atmosphere models to the optical data to check the consistency of the photometry for individual stars across the passbands and to estimate the line-of-sight extinction. Finally, we use the estimated line-of-sight extinctions to produce an extinction map across the Small Magellanic Cloud and investigate the nature of extinction as a function of stellar population.

We present photometry and a preliminary interpretation of a UBVR survey of the Large and Small Magellanic Clouds. We determine improved values for the relative number of blue and red supergiants. We also compare the relative number of Red Supergiants (RSGs) and Wolf-Rayet stars, demonstrating a strong, tight correlation with metallicity, and reinvestigate the initial mass function slope of massive stars found in the field.

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.

We report on the detection of optically thick free-free radio sources in the galaxies M33, NGC 253, and NGC 6946 using data in the literature. We interpret these sources as being young, embedded star birth regions, which are likely to be clusters of ultracompact HII regions. All 35 of the sources presented in this article have positive radio spectral indices alpha>0 suggesting an optically thick thermal bremsstrahlung emission arising in the HII region surrounding hot stars. Energy requirements indicate a range of a several to >500 O7V star equivalents powering each HII region. Assuming a Salpeter IMF, this corresponds to integrated stellar masses of 0.1--60,000 Msun. For roughly half of the sources in our sample, there is no obvious optical counterpart, giving further support for their deeply embedded nature. Their luminosities and radio spectral energy distributions are consistent with HII regions having electron densities from 1500 cm^-3 to 15000 cm^-3 and radii of 1 - 7 pc. We suggest that the less luminous of these sources are extragalactic ultracompact HII region complexes, those of intermediate luminosity are similar to W49 in the Galaxy, while the brightest will be counterparts to 30 Doradus. These objects constitute the lower mass range of extragalactic ``ultradense HII regions'' which we argue are the youngest stages of massive star cluster formation yet observed. This sample is beginning to fill in the continuum of objects between small associations of ultracompact HII regions and the massive extragalactic clusters that may evolve into globular clusters.

The current understanding of mechanical feedback is reviewed by evaluating the standard, adiabatic model for shell formation and evolution. This model is relevant to phenomena ranging from individual stellar-wind bubbles to galactic superwinds, forming the basis for our understanding of the multiphase ISM, IGM, and galactic evolutionary processes. Although significant discrepancies between the model and observation have been identified, to date there are none that require a fundamental revision. A variety of evidence, ranging over three orders of magnitude in spatial scale, is broadly consistent with the standard model. This includes kinematics of individual objects, observations of hot gas, the size distribution of HI shells, and outflow rates from starburst galaxies. However, some of the most pressing issues relating to shell evolution are still outstanding and obstruct efforts to resolve key questions like the fate of the hot gas.

Using photometry from the Ultraviolet Imaging Telescope (UIT) and photometry and spectroscopy from three ground-based optical datasets we have analyzed the stellar content of OB associations and field areas in and around the regions N 79, N 81, N 83, and N 94 in the LMC. We compare data for the OB association Lucke-Hodge 2 (LH 2) to determine how strongly the initial mass function (IMF) may depend on different photometric reductions and calibrations. We also correct for the background contribution of field stars, showing the importance of correcting for field star contamination in determinations of the IMF of star formation regions. It is possible that even in the case of an universal IMF, the variability of the density of background stars could be the dominant factor creating the differences between calculated IMFs for OB associations. We have also combined the UIT data with the Magellanic Cloud Photometric Survey to study the distribution of the candidate O-type stars in the field. We find a significant fraction, roughly half, of the candidate O-type stars are found in field regions, far from any obvious OB associations. These stars are greater than 2 arcmin (30 pc) from the boundaries of existing OB associations in the region, which is a distance greater than most O-type stars with typical dispersion velocities will travel in their lifetimes. The origin of these massive field stars (either as runaways, members of low-density star-forming regions, or examples of isolated massive star formation) will have to be determined by further observations and analysis.

We conducted an extensive search for Wolf-Rayet stars (W-Rs) in the SMC, using the same interference filter imaging techniques that have proved successful in finding W-Rs in more distant members of the Local Group. Photometry of some 1.6 million stellar images resulted in some 20 good candidates, which we then examined spectroscopically. Two of these indeed proved to be newly found W-Rs, bringing the total known in the SMC from 9 to 11. Other finds included previously unknown Of-type stars (one as early as O5f?p)),the recovery of the Luminous Blue Variable S18, and the discovery of a previously unknown SMC symbiotic star. More important, however, is the fact that there does not exist a significant number of W-Rs waiting to be discovered in the SMC. The number of W-Rs in the SMC is a factor of 3 lower than in the LMC (per unit luminosity), and we argue this is the result of the SMC's low metallicity on the evolution of the most massive stars.

Here we examine 12 Galactic clusters and OB associations containing Wolf-Rayet stars (W-Rs) and/or Luminous Blue Variables (LBVs), in order to assess the progentor masses of these evolved massive stars. We find that in the Milky Way, early-type WN stars evolve from stars with masses as low as 20Mo. Some early-type WN stars are found in clusters with very high turn-off masses, supporting our earlier results that the "WNE" class of W-Rs is a common stage in the evolution of stars with a wide range of mases. In contrast, the late-type WNs (particularly those of WN7 class) and the LBVs are only found in clusters with the very highest turn-off masses, suggesting that only the most massive stars evolve to these objects. An examination of the environment of the archetype LBV Eta Car reveals that this star is coeval with the rest of the Tr14/16 association, suggesting that its LBV properties is a normal consequence of evolution, and has little to do with the possibility that it is a binary.

We combine new CCD UBV photometry and spectroscopy with that from the literature to investigate 19 Magellanic Cloud OB associations that contain Wolf-Rayet (WR) and other types of evolved massive stars. Our spectroscopy reveals a wealth of newly identified interesting objects, including early O-type supergiants, a high mass double-lined binary in the SMC, and, in the LMC, a newly confirmed LBV (R 85), a newly discovered WR star (Sk-69 194), and a newly found luminous B[e] star (LH85-10). We use the data to evaluate the extent that each association is coeval, and for the ones with little age spread, determine the turn-off masses of the unevolved stars, as a means of finding the progenitor masses of LBVs and WRs. Our data is also used to find the (minimum) bolometric corrections that apply to these evolved objects. We find that WRs in the SMC come from only the highest mass stars (>70Mo), while in the LMC the WNEs come from a wide range of masses---possibly all stars with masses >30Mo pass though a WNE stage at LMC metallicities. In both the SMC and LMC the WCs appear to come from the same mass range as do the WNs, a finding that has important implications for interpreting the WC/WN ratio among Local Group galaxies. LBVs appear to come from the very highest mass stars. The Ofpe/WN9 stars, sometimes implicated in the "LBV phenonomon" come from very modest mass stars, 25-35Mo, suggesting a different origin for their variability. The BCs of WRs are considerably more negative than those of O-type stars, with BC(WNE)=-6.0 mag and BC(WC4)=-5.5mag; more modest values are found for the Ofpe/WN9 stars (-2 to -4 mag). These findings are consistent with recent analysis done with Hillier's "standard model".

We use the number of papers published in 1998 and 1999 to test the hypothesis that the queue observing mode at WIYN leads to a significantly higher scientific throughput than classical mode observing. We use the papers published from the 4-m, and papers published from the non-queue WIYN time as controls, requiring only that the data be obtained after 1996 August 1, at which time the WIYN queue was in its third full semester of operation, and the WIYN instruments functional and stable. The number of papers published from the queue data is actually 1.5 times smaller (on a per night basis) than from the 4-m, and roughly comparable to (but lower than) the number published from non-queue WIYN time. Thus neither comparison offers any support for the hypothesis that queue leads to a higher scientific throughput. The number of papers is relatively small, but the statistics are sufficiently robust to reject the possibility that queue observing at WIYN leads to a factor of 1.5 enhancement in publication rate with a 99.3% confidence in comparison to the 4-m, and with an 89.9% confidence in comparison with non-queue WIYN time. We consider several explanations, and urge other observatories planning to employ the queue mode include some controls to provide an objective evaluation of its success.