Both rotation and interactions with binary companions can significantly affect massive star evolution, altering interior and surface abundances, mass loss rates and mechanisms, observed temperatures and luminosities, and their ultimate core-collapse fates. The Geneva and BPASS stellar evolution codes include detailed treatments of rotation and binary evolutionary effects, respectively, and can illustrate the impact of these phenomena on massive stars and stellar populations. However, a direct comparison of these two widely-used codes is vital if we hope to use their predictions for interpreting observations. In particular, rotating and binary models will predict different young stellar populations, impacting the outputs of stellar population synthesis (SPS) and the resulting interpretation of large massive star samples based on commonly-used tools such as star count ratios. Here we compare the Geneva and BPASS evolutionary models, using an interpolated SPS scheme introduced in our previous work and a novel Bayesian framework to present the first in-depth direct comparison of massive stellar populations produced from single, rotating, and binary non-rotating evolution models. We calculate both models' predicted values of star count ratios and compare the results to observations of massive stars in Westerlund 1, $h + \chi$ Persei, and both Magellanic Clouds. We also consider the limitations of both the observations and the models, and how to quantitatively include observational completeness limits in SPS models. We demonstrate that the methods presented here, when combined with robust stellar evolutionary models, offer a potential means of estimating the physical properties of massive stars in large stellar populations.

We present the first results from a study of TESS Sector 1 and 2 light curves for eight evolved massive stars in the LMC: six yellow supergiants (YSGs) and two luminous blue variables (LBVs), including S Doradus. We use an iterative prewhitening procedure to characterize the short-timescale variability in all eight stars. The periodogram of one of the YSGs, HD 269953, displays multiple strong peaks at higher frequencies than its fellows. While the field surrounding HD 269953 is quite crowded, it is the brightest star in the region, and has infrared colors indicating it is dusty. We suggest HD 269953 may be in a post-red supergiant evolutionary phase. We find a signal with a period of $\sim5$ days for the LBV HD 269582. The periodogram of S Doradus shows a complicated structure, with peaks below frequencies of 1.5 cycles per day. We fit the shape of the background noise of all eight light curves, and find a red noise component in all of them. However, the power law slope of the red noise and the timescale over which coherent structures arise changes from star to star. Our results highlight the potential for studying evolved massive stars with TESS.

Populations of massive stars are directly reflective of the physics of stellar evolution. Counting subtypes of massive stars and ratios of massive stars in different evolutionary states have been used ubiquitously as diagnostics of age and metallicity effects. While the binary fraction of massive stars is significant, inferences are often based upon models incorporating only single-star evolution. In this work, we utilize custom synthetic stellar populations from the Binary Population and Stellar Synthesis (BPASS) code to determine the effect of stellar binaries on number count ratios of different evolutionary stages in both young massive clusters and galaxies with massive stellar populations. We find that many ratios are degenerate in metallicity, age, and/or binary fraction. We develop diagnostic plots using these stellar count ratios to help break this degeneracy, and use these plots to compare our predictions to observed data in the Milky Way and the Local Group. These data suggest a possible correlation between the massive star binary fraction and metallicity. We also examine the robustness of our predictions in samples with varying levels of completeness. We find including binaries and imposing a completeness limit can both introduce $\gtrsim0.1$ dex changes in inferred ages. Our results highlight the impact that binary evolution channels can have on the massive star population.

We have obtained near-simultaneous Swift/XRT imaging and Gemini GMOS spectroscopy for the ultraluminous X-ray source (ULX) NGC~300 ULX-1 (formerly designated SN~2010da). The observed X-ray emission is consistent with an inhomogeneous wind that partially obscures a central, bright inner accretion disk. We simultaneously fit eleven 0.3-10 keV spectra obtained over a $\sim$1 year time period (2016 April to 2017 July) using the same partial covering model, and find that although the covering fraction varies significantly (from 78% to consistent with 0%), the unabsorbed luminosity remains essentially constant across all observations ($2-6\times10^{39}$ erg s$^{-1}$). A relatively high 0.3-10 keV fractional variability amplitude ($F_{\rm var}$) of $\sim$30% is observed in all eleven observations. Optical spectra from Gemini exhibit numerous emission lines (e.g., H$\alpha$, H$\beta$, He II $\lambda$4686) which suggest that the neutron star primary is photoionizing material in the immediate vicinity of the binary. We compare the He II $\lambda$4686 line luminosity ($\sim7-9\times10^{35}$ erg s$^{-1}$) to the contemporaneous soft X-ray emission and find the X-ray emission is broadly consistent with the observed He II line luminosity. The combination of our X-ray observations and optical spectroscopy suggest that geometric beaming effects in the ULX-1 system are minimal, making ULX-1 one of only a few bona fide ULXs to be powered by accretion onto a neutron star.

We identify an object previously thought to be a star in the disk of M31, J0045+41, as a background $z\approx0.215$ AGN seen through a low-absorption region of M31. We present moderate resolution spectroscopy of J0045+41 obtained using GMOS at Gemini-North. The spectrum contains features attributable to the host galaxy. We model the spectrum to estimate the AGN contribution, from which we estimate the luminosity and virial mass of the central engine. Residuals to our fit reveal a blue-shifted component to the broad H$\alpha$ and H$\beta$ at a relative velocity of $\sim4800$ km s$^{-1}$. We also detect \ion{Na}{1} absorption in the Milky Way restframe. We search for evidence of periodicity using $g$-band photometry from the Palomar Transient Factory and find evidence for multiple periodicities ranging from $\sim80-350$ days. Two of the detected periods are in a 1:4 ratio, which is identical to the predictions of hydrodynamical simulations of binary supermassive black hole systems. If these signals arise due to such a system, J0045+41 is well within the gravitational wave regime. We calculate the time until inspiral due to gravitational radiation, assuming reasonable values of the mass ratio of the two black holes. We discuss the implications of our findings and forthcoming work to identify other such interlopers in the light of upcoming photometric surveys such as the Zwicky Transient Facility (ZTF) or the Large Synoptic Survey Telescope (LSST) projects.

Thanks to incredible advances in instrumentation, surveys like the Sloan Digital Sky Survey have been able to find and catalog billions of objects, ranging from local M dwarfs to distant quasars. Machine learning algorithms have greatly aided in the effort to classify these objects; however, there are regimes where these algorithms fail, where interesting oddities may be found. We present here an X-ray bright quasar misidentified as a red supergiant/X-ray binary, and a subsequent search of the SDSS quasar catalog for X-ray bright stars misidentified as quasars.