We investigate the effects of mass loss during the main-sequence (MS) and post-MS phases on the final black hole (BH) masses of massive stars. We compute solar metallicity Geneva stellar evolution models of an 85 $M_\odot$ star with mass-loss rate ($M_\odot$) prescriptions for MS and post-MS phases. Such models could lead to massive BHs such as the recently detected 70 $M_\odot$ BH in the LB-1 system. Based on the observational constraints for $M_\odot$ of luminous stars, we discuss two possible scenarios that could produce such a massive BH at high metallicity. First, if the progenitor of LB-1 evolved from the observed population of WNh stars, we show that its average mass loss rate during its post-MS evolution was less than $1\,\times10^{-5}\,M_\odot/$yr. However, this is lower than the typical observed mass-loss rates of LBVs. A second possibility is that the progenitor evolved from a yet undetected population of 80-85~\msun\ stars with strong surface magnetic fields, which could quench mass loss during the MS evolution. In this case, the average mass-loss rate during the post-MS luminous blue variable (LBV) phase has to be less than $5\,\times10^{-5}\,M_\odot/$yr. This value is still low, considering that LBVs such as AG Carinae typically have average mass-loss rates from quiescent stellar winds close or above this level. LBVs directly collapsing to massive BHs are apparently at odds with the evidence that LBVs are the direct progenitors of some supernovae (SNe). To reconcile this, we suggest that LBVs from single stars (or mergers that fully mixed) have large cores and form BHs, while binary LBVs may have smaller cores and larger envelopes (either mergers or mass gainers that do not fully mix) and could produce a SN.

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.