The `Red Supergiant Problem' describes the claim that the brightest Red Supergiant (RSG) progenitors to type II-P supernovae are significantly fainter than RSGs in the field. This mismatch has been interpreted by several authors as being a manifestation of the mass threshold for the production of black holes (BHs), such that stars with initial masses above a cutoff of $M_{\rm hi}=17$M$_\odot$ and below 25$M_\odot$ will die as RSGs, but with no visible SN explosion as the BH is formed. However, we have previously cautioned that this cutoff is more likely to be higher and has large uncertainties ($M_{\rm hi}=19^{+4}_{-2}M_{\odot}$), meaning that the statistical significance of the RSG Problem is less than $2\sigma$. Recently, Kochanek (2020) has claimed that our work is statistically flawed, and with their analysis has argued that the upper mass cutoff is as low as $M_{\rm hi} = 15.7 \pm 0.8M_\odot$, giving the RSG Problem a significance of $>10\sigma$. In this letter, we show that Kochanek's low cutoff is caused by a statistical misinterpretation, and the associated fit to the progenitor mass spectrum can be ruled out at the 99.6\% confidence level. Once this problem is remedied, Kochanek's best fit becomes $M_{\rm hi} =19^{+4}_{-2}M_{\odot}$, in excellent agreement with our work. Finally, we argue that, in the search for a RSG `vanishing' as it collapses directly to a BH, any such survey would have to operate for decades before the absence of any such detection became statistically significant.

We use the Binary Population and Spectral Synthesis (BPASS) models to test the recent suggestion that red supergiants can provide an accurate age estimate of a co-eval stellar population that is unaffected by interacting binary stars. Ages are estimated by using both the minimum luminosity red supergiant and the mean luminosity of red supergiants in a cluster. We test these methods on a number of observed star clusters and find our results in agreement with previous estimates. Importantly we find the difference between the ages derived from stellar population models with and without a realistic population of interacting binary stars is only a few 100,000 years at most. We find that the mean luminosity of red supergiants in a cluster is the best method to determine the age of a cluster because it is based o the entire red supergiant population rather than using only the least luminous red supergiant.

Evolutionary models have shown the substantial effect that strong mass-loss rates ($\dot{M}$) can have on the fate of massive stars. Red supergiant (RSG) mass-loss is poorly understood theoretically, and so stellar models rely on purely empirical \mdot-luminosity relations to calculate evolution. Empirical prescriptions usually scale with luminosity and effective temperature, but $\dot{M}$ should also depend on the current mass and hence the surface gravity of the star, yielding more than one possible $\dot{M}$ for the same position on the Hertzsprung-Russell diagram. One can solve this degeneracy by measuring $\dot{M}$ for RSGs that reside in clusters, where age and initial mass ($M_{\rm init}$) are known. In this paper we derive $\dot{M}$ values and luminosities for RSGs in two clusters, NGC 2004 and RSGC1. Using newly derived $M_{\rm init}$ measurements, we combine the results with those of clusters with a range of ages and derive an $M_{\rm init}$-dependent $\dot{M}$-prescription. When comparing this new prescription to the treatment of mass-loss currently implemented in evolutionary models, we find models drastically over-predict the total mass-loss, by up to a factor of 20. Importantly, the most massive RSGs experience the largest downward revision in their mass-loss rates, drastically changing the impact of wind mass-loss on their evolution. Our results suggest that for most initial masses of RSG progenitors, quiescent mass-loss during the RSG phase is not effective at removing a significant fraction of the H-envelope prior to core-collapse, and we discuss the implications of this for stellar evolution and observations of SNe and SN progenitors.

By comparing the properties of Red Supergiant (RSG) supernova progenitors to those of field RSGs, it has been claimed that there is an absence of progenitors with luminosities $L$ above $\log(L/L_\odot) > 5.2$. This is in tension with the empirical upper luminosity limit of RSGs at $\log(L/L_\odot) = 5.5$, a result known as the `Red Supergiant Problem'. This has been interpreted as evidence for an upper mass threshold for the formation of black-holes. In this paper, we compare the observed luminosities of RSG SN progenitors with the observed RSG $L$-distribution in the Magellanic Clouds. Our results indicate that the absence of bright SN II-P/L progenitors in the current sample can be explained at least in part by the steepness of the $L$-distribution and a small sample size, and that the statistical significance of the Red Supergiant Problem is between 1-2$\sigma$ . Secondly, we model the luminosity distribution of II-P/L progenitors as a simple power-law with an upper and lower cutoff, and find an upper luminosity limit of $\log(L_{\rm hi}/L_\odot) = 5.20^{+0.17}_{-0.11}$ (68\% confidence), though this increases to $\sim$5.3 if one fixes the power-law slope to be that expected from theoretical arguments. Again, the results point to the significance of the RSG Problem being within $\sim 2 \sigma$. Under the assumption that all progenitors are the result of single-star evolution, this corresponds to an upper mass limit for the parent distribution of $M_{\rm hi} = 19.2{\rm M_\odot}$, $\pm1.3 {\rm M_\odot (systematic)}$, $^{+4.5}_{-2.3} {\rm M_\odot}$ (random) (68\% confidence limits).

Galactic, young massive star clusters are approximately coeval aggregates of stars, close enough to resolve the individual stars, massive enough to have produced large numbers of massive stars, and young enough for these stars to be in a pre-supernova state. As such these objects represent powerful natural laboratories in which to study the evolution of massive stars. To be used in this way, it is crucial that accurate and precise distances are known, since this affects both the inferred luminosities of the cluster members and the age estimate for the cluster itself. Here we present distance estimates for three star clusters rich in Red Supergiants ($\chi$ Per, NGC 7419 and Westerlund 1) based on their average astrometric parallaxes $\bar{\pi}$ in Gaia Data Release 2, where the measurement of $\bar{\pi}$ is obtained from a proper-motion screened sample of spectroscopically-confirmed cluster members. We determine distances of $d=2.25^{+0.16}_{-0.14}$kpc, $d=3.00^{+0.35}_{-0.29}$kpc, and $d=3.87^{+0.95}_{-0.64}$kpc for the three clusters respectively. We find that the dominant source of error is that in Gaia's zero-point parallax offset $\pi_{\rm ZP}$, and we argue that more precise distances cannot be determined without an improved characterization of this quantity.

There is growing evidence that star clusters can no longer be considered simple stellar populations (SSPs). Intermediate and old age clusters are often found to have extended main sequence turn-offs (eMSTOs) which are difficult to explain with single age isochrones, an effect attributed to rotation. In this paper, we provide the first characterisation of this effect in young (<20Myr) clusters. We determine ages for 4 young massive clusters (2 LMC, 2 Galactic) by three different methods: using the brightest single turn-off (TO) star; using the luminosity function (LF) of the TO; and by using the lowest $L_{\rm bol}$ red supergiant (RSG). The age found using the cluster TO is consistently younger than the age found using the lowest RSG $L_{\rm bol}$. Under the assumption that the lowest luminosity RSG age is the `true' age, we argue that the eMSTOs of these clusters cannot be explained solely by rotation or unresolved binaries. We speculate that the most luminous stars above the TO are massive blue straggler stars formed via binary interaction, either as mass gainers or merger products. Therefore, using the cluster TO method to infer ages and initial masses of post-main sequence stars such as Wolf-Rayet stars, luminous blue variables and RSGs, will result in ages inferred being too young and masses too high.

It has been argued in the literature that the star HV~2112 in the Small Magellanic Cloud is the first known example of a T\.ZO, a Red Supergiant with a degenerate neutron core. This claim is based on the star having a high luminosity ($\log (L/L_\odot)$~> 5), an extremely cool effective temperature, and a surface enriched in in lithium, calcium and various $irp$-process elements. In this paper we re-examine this evidence, and present new measurements of the stellar properties. By compiling archival photometry from blue to mid-IR for HV~2112 and integrating under it's spectral energy distribution we find a bolometric luminosity in the range of $\log (L/L_\odot)$=4.70-4.91, lower than that found in previous work and comparable to bright asymptotic giant branch (AGB) stars. We compare a VLT+XSHOOTER spectrum of HV~2112 to other late type, luminous SMC stars, finding no evidence for enhancements in Rb, Ca or K, though there does seem to be an enrichment in Li. We therefore conclude that a much more likely explanation for HV~2112 is that it is an intermediate mass($\sim$5M$_\odot$) AGB star. However, from our sample of comparison stars we identify a new T\.ZO candidate, HV~11417, which seems to be enriched in Rb but for which we cannot determine a Li abundance.

The empirical upper luminosity boundary $L_{\rm max}$ of cool supergiants, often referred to as the Humphreys-Davidson limit, is thought to encode information on the general mass-loss behaviour of massive stars. Further, it delineates the boundary at which single stars will end their lives stripped of their hydrogen-rich envelope, which in turn is a key factor in the relative rates of Type-II to Type-Ibc supernovae from single star channels. In this paper we have revisited the issue of $L_{\rm max}$ by studying the luminosity distributions of cool supergiants (SGs) in the Large and Small Magellanic Clouds (LMC/SMC). We assemble samples of cool SGs in each galaxy which are highly-complete above $\log L/L_{\odot}$=5.0, and determine their spectral energy distributions from the optical to the mid-infrared using modern multi-wavelength survey data. We show that in both cases $L_{\rm max}$ appears to be lower than previously quoted, and is in the region of $\log L/L_{\odot}$=5.5. There is no evidence for $L_{\rm max}$ being higher in the SMC than in the LMC, as would be expected if metallicity-dependent winds were the dominant factor in the stripping of stellar envelopes. We also show that $L_{\rm max}$ aligns with the lowest luminosity of single nitrogen-rich Wolf-Rayet stars, indicating of a change in evolutionary sequence for stars above a critical mass. From population synthesis analysis we show that the Geneva evolutionary models greatly over-predict the numbers of cool SGs in the SMC. We also argue that the trend of earlier average spectral types of cool SGs in lower metallicity environments represents a genuine shift to hotter temperatures. Finally, we use our new bolometric luminosity measurements to provide updated bolometric corrections for cool supergiants.

The fate of massive stars with initial masses >8M$_\odot$ depends largely on the mass-loss rate (\mdot ) in the end stages of their lives. Red supergiants (RSGs) are the direct progenitors to Type II-P core collapse supernovae (SN), but there is uncertainty regarding the scale and impact of any mass-loss during this phase. Here we used near and mid-IR photometry and the radiative transfer code DUSTY to determine luminosity and \mdot\ values for the RSGs in two Galactic clusters (NGC 7419 and $\chi$ Per) where the RSGs are all of similar initial mass ($M_{\rm initial}$$\sim$16M$_\odot$), allowing us to study how \mdot\ changes with time along an evolutionary sequence. We find a clear, tight correlation between luminosity and \mdot\ suggesting the scatter seen in studies of field stars is caused by stars of similar luminosity being of different initial masses. From our results we estimate how much mass a 16M$_\odot$ star would lose during the RSG phase, finding a star of this mass would lose a total of 0.61$^{+0.92}_{-0.31}$M$_\odot$. This is much less than expected for \mdot\ prescriptions currently used in evolutionary models.

There are a growing number of nearby SNe for which the progenitor star is detected in archival pre-explosion imaging. From these images it is possible to measure the progenitor's brightness a few years before explosion, and ultimately estimate its initial mass. Previous work has shown that II-P and II-L supernovae (SNe) have Red Supergiant (RSG) progenitors, and that the range of initial masses for these progenitors seems to be limited to $<$17M$_\odot$. This is in contrast with the cutoff of 25-30M$_\odot$ predicted by evolutionary models, a result which is termed the 'Red Supergiant Problem'. Here we investigate one particular source of systematic error present in converting pre-explosion photometry into an initial mass, that of the bolometric correction (BC) used to convert a single-band flux into a bolometric luminosity. We show, using star clusters, that RSGs evolve to later spectral types as they approach SN, which in turn causes the BC to become larger. Failure to account for this results in a systematic underestimate of a star's luminosity, and hence its initial mass. Using our empirically motivated BCs we reappraise the II-P and II-L SNe that have their progenitors detected in pre-explosion imaging. Fitting an initial mass function to these updated masses results in an increased upper mass cutoff of $M_{\rm hi}$=$19.0^{+2.5}_{-1.3}$M$_\odot$, with a 95% upper confidence limit of $<$27M$_\odot$. Accounting for finite sample size effects and systematic uncertainties in the mass-luminosity relationship raises the cutoff to $M_{\rm hi}$=25M$_\odot$ ($<$33M$_\odot$, 95% confidence). We therefore conclude that there is currently no strong evidence for `missing' high mass progenitors to core-collapse SNe.

We present an abundance analysis of seven super-star clusters in the disk of M83. The near-infrared spectra of these clusters are dominated by Red Supergiants, and the spectral similarity in the J-band of such stars at uniform metallicity means that the integrated light from the clusters may be analysed using the same tools as those applied to single stars. Using data from VLT/KMOS we estimate metallicities for each cluster in the sample. We find that the abundance gradient in the inner regions of M83 is flat, with a central metallicity of [Z] = 0.21$\pm$0.11 relative to a Solar value of $Z_\odot$=0.014, which is in excellent agreement with the results from an analysis of luminous hot stars in the same regions. Compiling this latest study with our other recent work, we construct a mass-metallicity relation for nearby galaxies based entirely on the analysis of RSGs. We find excellent agreement with the other stellar-based technique, that of blue supergiants, as well as with temperature-sensitive (`auroral' or `direct') \hii-region studies. Of all the HII-region strong-line calibrations, those which are empirically calibrated to direct-method studies (N2 and O3N2) provide the most consistent results.

With red supergiants (RSGs) predicted to end their lives as Type IIP core collapse supernova (CCSN), their behaviour before explosion needs to be fully understood. Mass loss rates govern RSG evolution towards SN and have strong implications on the appearance of the resulting explosion. To study how the mass-loss rates change with the evolution of the star, we have measured the amount of circumstellar material around 19 RSGs in a coeval cluster. Our study has shown that mass loss rates ramp up throughout the lifetime of an RSG, with more evolved stars having mass loss rates a factor of 40 higher than early stage RSGs. Interestingly, we have also found evidence for an increase in circumstellar extinction throughout the RSG lifetime, meaning the most evolved stars are most severely affected. We find that, were the most evolved RSGs in NGC2100 to go SN, this extra extinction would cause the progenitor's initial mass to be underestimated by up to 9M$_\odot$.

The mass loss rates of red supergiants (RSGs) govern their evolution towards supernova and dictate the appearance of the resulting explosion. To study how mass-loss rates change with evolution we measure the mass-loss rates (\mdot) and extinctions of 19 red supergiants in the young massive cluster NGC2100 in the Large Magellanic Cloud. By targeting stars in a coeval cluster we can study the mass-loss rate evolution whilst keeping the variables of mass and metallicity fixed. Mass-loss rates were determined by fitting DUSTY models to mid-IR photometry from WISE and Spitzer/IRAC. We find that the \mdot\ in red supergiants increases as the star evolves, and is well described by \mdot\ prescription of de Jager, used widely in stellar evolution calculations. We find the extinction caused by the warm dust is negligible, meaning the warm circumstellar material of the inner wind cannot explain the higher levels of extinction found in the RSGs compared to other cluster stars. We discuss the implications of this work in terms of supernova progenitors and stellar evolution theory. We argue there is little justification for substantially increasing the \mdot\ during the RSG phase, as has been suggested recently in order to explain the absence of high mass Type IIP supernova progenitors. We also argue that an increase in reddening towards the end of the RSG phase, as observed for the two most evolved cluster stars, may provide a solution to the red supergiant problem.