In recent years, the viability of the pair-instability supernova (PISN) scenario for explaining superluminous supernovae has all but disappeared except for a few slowly-evolving examples. However, PISN are not predicted to be superluminous throughout the bulk of their mass range. In fact, it is more likely that the first PISN we see (if we have not seen one already) will not be superluminous. Here, we present hydrodynamic simulations of PISNe for four stellar models with unique envelope properties spanning the PISN mass range. In addition, we compute synthetic light curves for comparison with current and future observations. We also investigate, in the context of our most massive model, the prospect of mixing in the supernova ejecta alleviating discrepancies between current PISN models and the remaining superluminous candidate events. To this end, we present the first published 3D hydrodynamic simulations of PISNe. After achieving convergence between 1D, 2D, and 3D simulations we examine mixing in the supernova ejecta and its affect on the bolometric light curve. We observe slight deviations from spherical symmetry which increase with the number of dimensions. We find no significant effects on the bolometric light curve, however we conclude that mixing between the silicon and oxygen rich layers caused by the Rayleigh-Taylor instability may affect spectra.

The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500 M$_{\odot}$ Z $\sim$ 0.1 - 0.4 Z$_{\odot}$ stars. These stars lose up to 80% of their mass to strong line-driven winds and explode as bare He cores. We find that their light curves and spectra are quite different from those of Population III pair-instability explosions, which therefore cannot be used as templates for low-redshift events. Although non-zero metallicity pair-instability supernovae are generally dimmer than their Population III counterparts, in some cases they will be bright enough to be detected at the earliest epochs at which they can occur, the formation of the first galaxies at $z \sim$ 10 - 15. Others can masquerade as dim, short duration supernovae that are only visible in the local universe and that under the right conditions could be hidden in a wide variety of supernova classes. We also report for the first time that some pair-instability explosions can create black holes with masses of $\sim$ 100 M$_{\odot}$.

(shortened) We provide a grid of single star models covering a mass range from 0.8 to 120 Msun with an initial metallicity Z = 0.002 with and without rotation. We discuss the impact of a change in the metallicity by comparing the current tracks with models computed with exactly the same physical ingredients but with a metallicity Z = 0.014 (solar). We show that the width of the main-sequence (MS) band in the upper part of the Hertzsprung-Russell diagram (HRD), for luminosity above log(L/Lsun) > 5.5, is very sensitive to rotational mixing. Strong mixing significantly reduces the MS width. We confirm, but here for the first time on the whole mass range, that surface enrichments are stronger at low metallicity provided that comparisons are made for equivalent initial mass, rotation and evolutionary stage. We show that the enhancement factor due to a lowering of the metallicity (all other factors kept constant) increases when the initial mass decreases. Present models predict an upper luminosity for the red supergiants (RSG) of log (L/Lsun) around 5.5 at Z = 0.002 in agreement with the observed upper limit of RSG in the Small Magellanic Cloud. We show that models using shear diffusion coefficient calibrated to reproduce the surface enrichments observed for MS B-type stars at Z = 0.014 can also reproduce the stronger enrichments observed at low metallicity. In the framework of the present models, we discuss the factors governing the timescale of the first crossing of the Hertzsprung gap after the MS phase. We show that any process favouring a deep localisation of the H-burning shell (steep gradient at the border of the H-burning convective core, low CNO content) and/or the low opacity of the H-rich envelope favour a blue position in the HRD for the whole or at least a significant fraction of the core He-burning phase.

There is observational evidence that supports the existence of Very Massive Stars in the local universe. First, very massive stars (Mini<=320 M) have been observed in the Large Magellanic Cloud . Second, there are observed SNe that bear the characteristics of Pair Creation Supernovae which have very massive stars as progenitors. The most promising candidate to date is SN2007bi. In order to investigate the evolution and fate of nearby very massive stars, we calculated a new grid of models for such objects, for solar, LMC and SMC metallicities, which covers the initial mass range from 120 to 500M. Both rotating and non-rotating models were calculated using the Geneva stellar evolution code and evolved until at least the end of helium burning and for most models until oxygen burning. Since very massive stars have very large convective cores during the Main-Sequence phase, their evolution is not so much affected by rotational mixing, but more by mass loss through stellar winds. Their evolution is never far from a homogeneous evolution even without rotational mixing. All the VMS, at all the metallicities studied here, end their life as WC(WO) type stars. At solar metallicity, none of our models is expected to explode as a PCSN. At the metallicity of the LMC, only stars more massive than 300 M are expected to explode as PCSNe. At the SMC metallicity, the mass range for the PCSN progenitors is much larger and comprises stars with initial masses between about 100 and 290 M . All VMS stars in the metallicity range studied here produce either a type Ib or a type Ic SN. We estimate that the progenitor of SN2007bi, assuming a SMC metallicity, had an initial mass between 160 and 175 M . None of models presented in this grid produce GRBs or magnetars. They lose too much angular momentum by mass loss or avoid the formation of a BH by producing a completely disruptive PCSN.

We use contemporary evolutionary models for Very Massive Stars (VMS) to assess whether the Eddington limit constrains the upper stellar mass limit. We also consider the interplay between mass and age for the wind properties and spectral morphology of VMS, with reference to the recently modified classification scheme for O2-3.5If*/WN stars. Finally, the death of VMS in the local universe is considered in the context of pair instability supernovae.

We recently determined the mass of the most massive star known to the date, R136a1 with a mass at birth 320 times the mass of our sun, as well as the mass of several other stars that are more massive than 150 M. Such massive stars (~150-300 M) may end their life as pair-instability supernovae (PISN) if they retain enough mass until they die. We have calculated a grid of stellar evolution models in order to investigate the impact of mass loss and rotation on the evolution and fate of these very massive stars. As mass loss is very strong at solar metallicity, our models predict that most of the very massive stars will die as type Ic SNe. Only slowly and non-rotating stars at metallicities below that of the LMC might retain enough mass to produce a PISN. This would mean that the first stellar generations might have produced PISN although their chemical signature is not observed in extremely metal poor stars in the halo of our galaxy.

Spectroscopic analyses of H-rich WN5-6 stars within the young star clusters NGC 3603 and R136 are presented, using archival HST & VLT spectroscopy, & high spatial resolution near-IR photometry. We derive high T* for the WN stars in NGC 3603 (T*~42+/-2 kK) & R136 (T*~53+/-3 kK) plus clumping-corrected dM/dt ~ 2-5x10^-5 Msun/yr which closely agree with theoretical predictions. These stars make a disproportionate contribution to the global budget of their host clusters. R136a1 alone supplies ~7% of N(LyC) of the entire 30 Dor region. Comparisons with stellar models calculated for the main-sequence evolution of 85-500 Msun suggest ages of ~1.5 Myr & M_init in the range 105 - 170 Msun for 3 systems in NGC 3603, plus 165-320 Msun for 4 stars in R136. Our high stellar masses are supported by dynamical mass determinations for the components of NGC 3603 A1. We consider the predicted L_X of the R136 stars if they were close, colliding wind binaries. R136c is consistent with a colliding wind binary system. However, short period, colliding wind systems are excluded for R136a WN stars if mass ratios are of order unity. Widely separated systems would have been expected to harden owing to early dynamical encounters with other massive stars in such a dense environment. From simulated star clusters, whose constituents are randomly sampled from the Kroupa IMF, both clusters are consistent with a tentative upper mass limit of ~300 Msun. The Arches cluster is either too old, exhibits a deficiency of very massive stars, or more likely stellar masses have been underestimated - M_init for the most luminous stars in the Arches cluster approach 200 Msun according to contemporary stellar & photometric results. The potential for stars greatly exceeding 150 Msun within metal-poor galaxies suggests that such pair-instability SNe could occur within the local universe, as has been claimed for SN 2007bi (abridged).