We investigate a suspected very massive star in one of the most metal-poor dwarf galaxies, PHL~293B. Excitingly, we find the sudden disappearance of the stellar signatures from our 2019 spectra, in particular the broad H lines with P~Cygni profiles that have been associated with a massive luminous blue variable (LBV) star. Such features are absent from our spectra obtained in 2019 with the ESPRESSO and X-shooter instruments of the ESO's VLT. We compute radiative transfer models using CMFGEN that fit the observed spectrum of the LBV and are consistent with ground-based and archival Hubble Space Telescope photometry. Our models show that during 2001--2011 the LBV had a luminosity $L_* = 2.5-3.5 \times 10^6 ~L_{\odot}$, a mass-loss rate $\dot{M} = 0.005-0.020 ~M_{\odot}$~yr$^{-1}$, a wind velocity of 1000~km~s$^{-1}$, and effective and stellar temperatures of $T_\mathrm{eff} = 6000-6800$~K and $T_\mathrm{*}=9500-15000$~K. These stellar properties indicate an eruptive state. We consider two main hypotheses for the absence of the broad emission components from the spectra obtained since 2011. One possibility is that we are seeing the end of an LBV eruption of a surviving star, with a mild drop in luminosity, a shift to hotter effective temperatures, and some dust obscuration. Alternatively, the LBV could have collapsed to a massive black hole without the production of a bright supernova.

We compute an extensive set of early-time spectra of supernovae interacting with circumstellar material using the radiative transfer code CMFGEN. Our models are applicable to events observed from 1 to a few days after explosion. Using these models, we constrain the progenitor and explosion properties of a sample of 17 observed interacting supernovae at early-times. Because massive stars have strong mass loss, these spectra provide valuable information about supernova progenitors, such as mass-loss rates, wind velocities, and surface abundances. We show that these events span a wide range of explosion and progenitor properties, exhibiting supernova luminosities in the 1e8 to 1e12 Lsun range, temperatures from 10000 to 60000 K, progenitor mass-loss rates from a few 1e-4 up to 1 Msun/yr, wind velocities from 100 to 800 km/s, and surface abundances from solar-like to H-depleted. Our results suggest that many progenitors of supernovae interacting with circumstellar material have significantly increased mass-loss before explosion compared to what massive stars show during the rest of their lifetimes. We also infer a lack of correlation between surface abundances and mass-loss rates. This may point to the pre-explosion mass-loss mechanism being independent of stellar mass. We find that the majority of these events have CNO-processed surface abundances. In the single star scenario this points to a preference towards high-mass RSGs as progenitors of interacting SNe, while binary evolution could impact this conclusion. Our models are publicly available and readily applicable to analyze results from ongoing and future large scale surveys such as the Zwicky Transient Factory.

The effects of rotation on stellar evolution are particularly important at low metallicity, when mass loss by stellar winds diminishes and the surface enrichment due to rotational mixing becomes relatively more pronounced than at high metallicities. Here we investigate the impact of rotation and metallicity on stellar evolution. Using a similar physics as in our previous large grids of models at Z=0.002 and Z=0.014, we compute stellar evolution models with the Geneva code for rotating and nonrotating stars with initial masses (Mini) between 1.7 and 120 Msun and Z=0.0004 (1/35 solar). This is comparable to the metallicities of the most metal poor galaxies observed so far, such as I Zw 18. Concerning massive stars, both rotating and nonrotating models spend most of their core-helium burning phase with an effective temperature higher than 8000 K. Stars become red supergiants only at the end of their lifetimes, and few RSGs are expected. Our models predict very few to no classical Wolf-Rayet stars as a results of weak stellar winds at low metallicity. The most massive stars end their lifetimes as luminous blue supergiants or luminous blue variables, a feature that is not predicted by models with higher metallicities. Interestingly, due to the behavior of the intermediate convective zone, the mass domain of stars producing pair-instability supernovae is smaller at Z=0.0004 than at Z=0.002. We find that during the main sequence phase, the ratio between nitrogen and carbon abundances (N/C) remains unchanged for nonrotating models. However, N/C increases by factors of 10-20 in rotating models at the end of the MS. Cepheids coming from stars with Mini > 4-6 Msun are beyond the core helium burning phase and spend little time in the instability strip. Since they would evolve towards cooler effective temperatures, these Cepheids should show an increase of the pulsation period as a function of age.

Observational surveys are now able to detect an increasing number of transients, such as core-collapse supernovae (SN) and powerful non-terminal outbursts (SN impostors). Dedicated spectroscopic facilities can follow up these events shortly after detection. Here we investigate the properties of these explosions at early times. We use the radiative transfer code CMFGEN to build an extensive library of spectra simulating the interaction of supernovae and their progenitor's winds/circumstellar medium (CSM). We consider a range of progenitor mass-loss rates (Mdot = 5e-4 to 1e-2 Msun/yr), abundances (solar, CNO-processed, and He-rich), and SN luminosities (L = 1.9e8 to 2.5e10 Lsun). The models simulate events ~1 day after explosion, and we assume a fixed location of the shock front as Rin=8.6e13 cm. We show that the large range of massive star properties at the pre-SN stage causes a diversity of early-time interacting SN and impostors. We identify three main classes of early-time spectra consisting of relatively high-ionisation (e.g. Ovi), medium-ionisation (e.g. Ciii), and low-ionisation lines (e.g. Feii/iii). They are regulated by L and the CSM density. Given a progenitor wind velocity Vinf, our models also place a lower limit of Mdot > 5e-4 (Vinf/150 km/s) Msun/yr for detection of CSM interaction signatures in observed spectra. Early-time SN spectra should provide clear constraints on progenitors by measuring H, He, and CNO abundances if the progenitors come from single stars. The connections are less clear considering the effects of binary evolution. Yet, our models provide a clear path for linking the final stages of massive stars to their post-explosion spectra at early times, and guiding future observational follow-up of transients with facilities such as the Zwicky Transient Facility.

In this paper we analyse the pre-explosion spectrum of SN2015bh by performing radiative transfer simulations using the CMFGEN code. This object has attracted significant attention due to its remarkable similarity to SN2009ip in both its pre- and post-explosion behaviour. They seem to belong to a class of events for which the fate as a genuine core-collapse supernova or a non-terminal explosion is still under debate. Our CMFGEN models suggest that the progenitor of SN2015bh had an effective temperature between 8700 and 10000 K, luminosity in the range ~ 1.8-4.74e6 Lsun, contained at least 25% H in mass at the surface, and half-solar Fe abundances. The results also show that the progenitor of SN 2015bh generated an extended wind with a mass-loss rate of ~ 6e-4 to 1.5e-3 Msun/yr and a velocity of 1000 km/s. We determined that the wind extended to at least 2.57e14 cm and lasted for at least 30 days prior to the observations, releasing 5e-5 Msun into the circumstellar medium. In analogy to 2009ip, we propose that this is the material that the explosive ejecta could interact at late epochs, perhaps producing observable signatures that can be probed with future observations. We conclude that the progenitor of SN 2015bh was most likely a warm luminous blue variable of at least 35 Msun before the explosion. Considering the high wind velocity, we cannot exclude the possibility that the progenitor was a Wolf-Rayet star that inflated just before the 2013 eruption, similar to HD5980 during its 1994 episode. If the star survived, late-time spectroscopy may reveal either a similar LBV or a Wolf-Rayet star, depending on the mass of the H envelope before the explosion. If the star exploded as a genuine SN, 2015bh would be a remarkable case of a successful explosion after black-hole formation in a star with a possible minimum mass 35 Msun at the pre-SN stage.