arxiv: v1 [astro-ph.he] 20 Mar 2014

Transcription

1 Astronomy & Astrophysics manuscript no. pisnlc5 c ESO 216 November 16, 216 Observational properties of low redshift pair instability supernovae A. Kozyreva 1, S. Blinnikov 2,3, N. Langer 1, and S.-C. Yoon 4 1 Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, Bonn, Germany kozyreva@astro.uni-bonn.de 2 Institute for Theoretical and Experimental Physics, Bolshaya Cheremushkinskaya Str. 25, Moscow, Russia 3 Kavli Institute for the Physics and Mathematics of the Universe (WPI), The University of Tokyo, Kashiwanoha, Kashiwa, Chiba , Japan 4 Department of Physics & Astronomy, Seoul National University, Gwanak-ro 1, Gwanak-gu, , Seoul, South Korea arxiv: v1 [astro-ph.he] 2 Mar 214 Preprint online version: November 16, 216 ABSTRACT Context. So called superluminous supernovae have been recently discovered in the local Universe. It appears possible that some of them originate from stellar explosions induced by the pair instability mechanism. Recent stellar evolution models also predict pair instability supernovae from very massive stars at fairly high metallicities (i.e., Z.4). Aims. We provide supernova models and synthetic light curves for two progenitor models, a 15 M red-supergiant and a 25 M yellow-supergiant at a metallicity of Z =.1, for which the evolution from the main sequence to collapse, and the initiation of the pair instability supernova (PISN) itself, has been previously computed in a realistic and self-consistent way. Methods. We are using the radiation hydrodynamics code STELLA to describe the supernova evolution of both models over a time frame of about 5 days. Results. We describe the shock-breakout phases of both supernovae which are characterized by a higher luminosity, a longer duration and a lower effective temperature than those of ordinary Type IIP supernovae. We derive the bolometric as well as the U, B, V, R and I light curves of our pair instability supernova models, which show a long-lasting plateau phase with maxima at M bol 19.3 mag and 21.3 mag for our lower and higher mass model, respectively. While we do not produce synthetic spectra, we also describe the photospheric composition and velocity as function of time. Conclusions. We conclude that the light curve of the explosion of our initially 15 M star resembles those of relatively bright type IIP supernovae, whereas its photospheric velocity at early times is somewhat smaller. Its 56 Ni mass of.4 M also falls well into the range found in ordinary core collapse supernovae. The light curve and photospheric velocity of our 25 M models has a striking resemblance with that of the superluminous SN 27bi, strengthening its interpretation as pair instability supernova. We conclude that pair instability supernovae may occur more frequently in the local universe than previously assumed. Key words. stars: massive stars: evolution stars: supernovae: super-luminous supernovae supernovae: pair instability supernovae supernovae: general 1. Introduction The final fate of very massive stars with initial masses between approximately 14 M and 26 M has been studied in many papers (Fowler & Hoyle 1964; Bisnovatyi-Kogan & Kazhdan 1967; Rakavy & Shaviv 1967; Barkat et al. 1967; Fraley 1968; El Eid & Langer 1986; ger et al. 23). Such massive stars undergo the dynamical instability caused by the creation of electron-positron pairs in oxygen cores, if they can retain their oxygen core masses high enough ( 6 M ) until carbon exhaustion at the center (ger et al. 23). This leads to explosive oxygen burning that eventually causes complete disruption of the star without leaving a compact remnant behind. Theoretical models predict that these pair instability supernovae (PISNe) can be much more energetic and luminous than ordinary SNe. Explosion energies of up to 1 53 erg and masses of radioactive nickel up to 4 M (e.g., ger & Woosley 22) are found to be achieved depending on the progenitor mass. The corresponding light curves are characterized by a longduration of several hundreds of days and luminosities of up to erg s 1 (Scannapieco et al. 25; Kasen et al. 211; Dessart et al. 213; Whalen et al. 213a,c). It is believed that PISNe are particularly relevant to the first generations of stars in the early Universe. Theoretical studies indicate that a significant fraction of the first stars would be massive enough to be potential progenitors of PISNe, mainly because of the lack of efficient coolants in the starforming regions in the early Universe (e.g., Bromm et al. 1999; Nakamura & Umemura 21; Abel et al. 22; Omukai & Palla 23; O Shea & Norman 26; Ohkubo et al. 29). Such very massive stars in the early Universe would not lose much mass during the pre-supernova evolutionary stages in favor of PISN production, because metal-free massive stars are not supposed to have strong line-driven winds (Krtička & Kubát 26) and because they are expected to be stable against pulsations (Baraffe et al. 21). Several numerical studies have been therefore presented to discuss the nature and detectability of PISNe in the early Universe (Scannapieco et al. 25; Kasen et al. 211; Dessart et al. 213; Whalen et al. 213a). These studies considered a variety of PISN progenitors including red supergiants (RSG), blue supergiants (BSG) and pure helium stars, but the considered metallicities of these progenitor models were limited to zero or very small values (Z 1 4 ). However, several recent studies indicate that PISNe are likely to occur not only in the early Universe, but also in the local 1

2 Table 1: Characteristics of the PISN models: name of the model, initial metallicity Z, initial/final mass (M ini /M f ) in solar masses, radius of the star at the onset of the radiative calculation R in solar radii, explosion (E expl ) and kinetic (E kin ) energy in foe (= 1 51 erg), specific energy (E/M) in units [1 5 erg/m ], the velocity at the outer edge of the 56 Ni-rich layer in km s 1 (v ni,max ), bulk yields of the isotopes in the ejecta in solar masses (hydrogen 1 H, helium 4, carbon 12 C, oxygen 16 O, silicon 28 Si, nickel 56 Ni). The names of the models starting with R indicate red supergiants. The names starting with B are blue supergiants. Models labeled.d (Dessart et al. 213) and labeled.k (Kasen et al. 211) are given for comparison. Name Z M in /M f R E expl E kin E/M v ni,max 1 H 4 12 C 16 O 28 Si 56 Ni (M ) (R ) (foe) (foe) (km/s) (M ) (M ) (M ) (M ) (M ) (M ) / R15.K / no data R19.D / R25.K / B19.D / B21.D / / B25.K 25/ Universe where metallicity is systematically higher than in the environments at high redshift. Crowther et al. (21) found several very massive stars (VMS) with initial masses of 15 M M 32 M in the Large Magellanic Clouds (LMC), which are potential PISN progenitors. The final fate of such VMS stars is critically determined by mass loss (e.g., Vink et al. 211). Given the strong metallicity dependence of the stellar winds mass loss rate (e.g., Mokiem et al. 27), it is generally believed that VMSs can not retain enough mass to produce PISNe for high metallicity. However, Langer et al. (27) point out that the production of PISNe does not necessarily require extremely metal poor environments, although low metallicity is still preferred. Using stellar evolution models, they argue that the metallicity threshold below which PISNe may occur (Z PISN ) can be as high as Z /3 (Langer 29). More recently Yusof et al. (213) drew a similar conclusion. The recent discovery of several superluminous SNe (SLSNe) that cannot be easily explained by usual core-collapse and/or interaction supernovae (see Gal-Yam 212a, for a recent review) also provides evidence for PISN in the local Universe. In particular, the observed properties of SN 27bi and SN seem to imply a large amount of radioactive nickel in these supernovae (i.e, more than 3 M of 56 Ni), for which PISN explosion gives one of the best explanations (Gal-Yam et al. 29; Young et al. 21; Cooke et al. 212). Given that some other possibilities like very energetic core-collapse explosion or magnetar-driven SN have also been suggested to explain SLSNe (Moriya et al. 21; Woosley 21; Kasen & Bildsten 21) and that future observational SN surveys will discover more diverse SLSN events in nearby galaxies, detailed studies on the observable properties of PISNe in the local Universe are still needed to have a solid conclusion on the association of SLSNe and PISNe (e.g., Dessart et al. 212b). Langer et al. (27) presented PISN progenitor models with initial masses of 15 M and 25 M, which roughly represent the low- and high-mass ends of PISN progenitors, at a metallicity of Z=.1 to discuss the possible event rate of PISNe in the local Universe. Kozyreva et al. (214) investigated the explosions and consequent nucleosynthesis of these models to discuss the implications of PISNe for the chemical evolution. In this paper, we explore these models further to investigate their observable properties using the radiation hydrodynamics code STELLA (Blinnikov et al. 26). These models together with those from rzig et al. (199) and Whalen et al. (213b) are among the highest metallicity PISN light curve models available in the literature, and therefore useful to identify PISN events at low redshift. We explain the method of calculations in Section 2 and discuss the produced light curves in Section 3. In Section 4, we compare our results with other synthetic light curves for PISNe and with observational light curves of usual and unusual core collapse SNe (CCSNe). We conclude our work in Section Evolutionary models and light curves modeling 2.1. Description of the evolutionary models For simulating the light curves of PISNe we use the results of evolutionary calculations produced with the Binary Evolution Code (BEC, Langer et al. 27; Kozyreva et al. 214). These are initially 15 M and 25 M models at a metallicity of Z = 1 3, with an initial rotation velocity of 1 km s 1. The evolution of the stars was calculated all the way from the zero age main sequence to the thermonuclear explosion caused by the pair instability. To follow the nucleosynthesis, Kozyreva et al. (214) used the Torch nuclear network developed by Timmes (1998, 1999) using 2 isotopes. Our progenitor models are red and yellow supergiants for the 15 M and 25 M models, respectively (Langer et al. 27; Kozyreva et al. 214). The main characteristics of the models are summarized in Table 1. Our stellar models lose a large fraction of their initial mass during core hydrogen and core helium burning due to line-driven winds (Langer et al. see 27, for details and Kudritzki et al. 1989; Vink et al. 21; Vink & de Koter 25; Yoon et al. 21). The average mass-loss rate is about M yr 1. By the onset of the pair instability our 15 M star has lost 56 M and our 25 M star 81 M. We show the pre-supernova evolution of stellar mass in Figure 3. The stellar evolutionary models were calculated with the assumption of semi-convection using a large semi-convective mixing parameter (α SEM = 1), and without convective core overshooting. The mixing length parameter was chosen to be 1.5 (Yoon et al. 26). Note that we neglected any convective mixing during the pair instability phase (Kozyreva et al. 214). The input models for the STELLA calculations already contain the shock wave generated at the boundary between the oxygen and helium shells emerging from the interaction between the rapidly expanding inner core and the nearly static outer 2

3 6 26 log ρ, g cm Stellar mass,solar masses Hydrogen burning lium burning Mass, solar masses Fig. 1: The density structure of the 15 M and 25 M PISN progenitor models at the onset of the STELLA calculations Time before the explosion Fig. 3: The pre-supernova evolution of the stellar mass of our models due to stellar wind mass loss. 1 3 helium-hydrogen envelope (Kozyreva et al. 214). The physics of the pair instability explosion mechanism is well understood and smoothly reproduced by evolutionary simulations without the need of invoking artificial assumptions. This is contrasted to other SN types that often need free parameters for supernova modeling, such as the critical density for the transition from deflagration to detonation for type Ia supernovae and the mass cut, degree of chemical mixing, and explosion energy for corecollapse supernovae. In Figures 1 and 2, we show the density and the chemical structure of our progenitor models at the onset of the STELLA calculations. At this moment the matter below the helium shell expands nearly homologously (v r). In Figure 2, we truncated the outer hydrogen-helium envelope which moves at very low velocity compared to the fast moving inner ejecta. We also do not show in this figure the bottom of the hydrogen-helium envelope where the shock is located and causes a strong discontinuities. The ejecta in the 15 M model contain.4 M of 56 Ni in the innermost region, moving at a velocity of less than 5 km s 1. The 25 M model contains 19 M of nickel over an extended region with km s 1 v 5 km s 1, which corresponds to M M r 5 M. We do not apply any mixing in our models. Later we discuss the impact of this chemical structure on light curves. There is a small amount of helium in the innermost region of the 25 M star as seen in Figure 2. This helium was produced by photo-disintegration of heavy elements during the explosive burning. We summarize the main properties of the PISNe in Table 1 along with other PISN models that we use for comparison Simulation of theoretical light curves and SEDs 1 For simulating the hydrodynamic evolution during the pair instability explosion we use the one-dimensional (1D) multigroup radiation Lagrangian implicit hydrodynamics code STELLA (Blinnikov et al. 26 and references therein). The code solves the hydrodynamic equations coupled with the radiative transfer equations without assuming radiative equilibrium. The STELLA code uses one temperature for the matter and has no specific 1 SED spectral energy distribution. temperature for radiation. The non-steady radiative transfer is solved for each of the one hundred frequency groups in all radial zones. The energy groups are uniformly distributed in logarithmic scale (i.e. in geometric progression) between the maximum wavelength of 5, Å and the minimum wavelength of 1 Å. The light curves are computed by integration of fluxes calculated with the STELLA code with UBVRI Bessel filter function using logarithmic interpolation. The opacity is calculated for each frequency group taking about 16 spectral lines into account according to Kurucz data (Kurucz 1991). The opacity also includes photoionization, free-free absorption and electron scattering assuming local thermodynamical equilibrium. Because of the large velocity gradient the opacity is calculated accounting for the effect of an expanding medium following Friend & Castor (1983) and Eastman & Pinto (1993). Local thermodynamic equilibrium is assumed in the plasma allowing to use the Boltzmann-Saha s distribution for ionization and level populations. This is needed for determining absorption, scattering and emission coefficients. Gamma-ray transfer is calculated using the one-group approximation for the non-local deposition of energy from the radioactive decay (Ambwani & Sutherland 1988). The code treats strong discontinuities (shock propagation) with an artificial viscosity term. The STELLA code does not follow any nuclear reactions. To map the BEC models into the STELLA code we remeshed the models. The original number of zones was reduced from 1931 BEC zones to 242 STELLA zones for the 15 M model and from 222 BEC zones to 276 STELLA zones for the 25 M model. We reduce the number of isotopes from 2 (from hydrogen to germanium) to 16 (H,, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Fe, Co, Ni). The envelopes of the progenitor models are optically thick and photons are not able to leave them on a short time scale. Therefore, STELLA radiation transport calculations are computationally very expensive for an input model in which a shock is located far from the photosphere. On the other hand, using a more evolved input model in which the shock front is closer to the photosphere causes more numerical instabilities just behind the shock because the BEC code does not include artificial vis- 3

4 Isotope Mass fraction Si S ArCa Ni O Mg Ne Velocity, 1 km/s C H Isotope Mass fraction Ni ArCa Si S O Velocity, 1 km/s Mg Ne C H Fig. 2: The chemical structure of the exploding 15 M and 25 M stars at metallicity Z = 1 3 at the onset of the STELLA calculations. Absolute magnitude Mbol R 2 U+2 B+1 V Mbol I Absolute magnitude Mbol R 2 I 1 B+1 U V Mbol Fig. 4: Bolometric and multiband (U, B, V, I, R) light curves for 15 M and 25 M PISNe at metallicity Z= 1 3. U, B, V, I, R magnitudes are plotted with a shift of+2,+1,, -1, -2 magnitudes, respectively. cosity. We therefore compared two different cases to test how the choice of the starting model affects the final outcome: one in which the shock front is located at the bottom of the helium shell and the other in which the shock is propagating in the middle of the helium envelope. We find that the final solutions differ only by less than 1% in the terms of bolometric luminosity. In the following we use the model in which the shock reached the bottom of hydrogen-helium envelope. This defines the point t= for all the figures. Note that the zero point in time (t = ) is not the same as the time of the explosion. Nevertheless, the time between the onset of the pair instability explosion (namely the highest central density) and the beginning of the light curve calculations is 1 hour for 15 M model and 44 seconds for 25 M model. These are the time intervals that the shock takes to propagate through the shallow pure helium layer to reach the bottom of hydrogenhelium envelope. These intervals are relatively small compared to the light curve evolution (months), therefore, one can consider the time in all figures as approximate time since the onset of the pair instability explosion. 3. Results In Figure 4, we plot the bolometric and U, B, V, R, I band light curves for the 15 M and 25 M models. reafter, in figures and tables we designate the models as Model and Model respectively. The U, B, V, R, I magnitudes are plotted with a shift of+2,+1,, -1, -2 magnitudes, respectively. The bolometric curves follow closely the V-band light curves, except for early time where it needs ultraviolet bolometric corrections, and for late time where it needs infrared bolometric corrections (Bersten & Hamuy 29). We summarize the details about shock breakout events (duration, maximum luminosity, effective and color temperatures, spectral wavelength peak, spectral 4

5 O O Si log Mass fraction at photosphere 1 2 H C Ne Mg Si S log Mass fraction at photosphere 1 2 H C Mg Ne Ar Ca S Fe Co Ni Fig. 5: Evolution of the composition at the receding electron-scattering photosphere for the 15 M and 25 M PISN models. Bolometric magnitude Time, hours Fig. 6: Shock breakout events for 15 M and 25 M PISN explosions. energy peak) and light curves in comparison to other radiative simulations of PISN explosion (Kasen et al. 211; Dessart et al. 213) in Table 2. For the discussion on the light curves below, we also show the evolution of the chemical composition at the photosphere which recedes along the decreasing Lagrangian mass coordinate, in Figure 5. The photosphere is defined as the mass zone where the electron-scattering optical depth in B-band turns above 2/3. This consideration is valid while matter and radiation temperatures are comparable at an accuracy better than 1%. The estimated period during which such conditions are present is approximately the first 35 days for Model and 25 days for Model. At later times, the opacity is dominated by line opacity and non-local thermodynamical equilibrium consideration is required. Soon after the explosion, the shock wave forms when the fast expanding inner part of the star encounters the slowly moving helium layer and the hydrogen-helium envelope. Because the whole envelope is optically thick it takes a long time for photons generated during the explosive burning to reach the surface (optically transparent medium). Therefore, firstly, during about 6 days for Model and 12 hours for Model the shock passes through the entire envelope (see e.g. Shigeyama et al. 1987; Blinnikov et al. 2). The shock wave reaches the surface, and the envelope matter behind the shock is heated and ionized. When the shock emerges at the surface the energetic shock breakout flash of ultraviolet and X-ray emission appears with a duration of a few hours. The shock breakout is shown in the Figure 6. We draw some important conclusions about the shock breakout in Section The 15 M model Soon after the shock breakout, the temperature drops rapidly due to adiabatic cooling. Later the recombination losses become comparable to adiabatic cooling (Grasberg et al. 1971). The recession of the photosphere along the Lagrangian mass coordinate is compensated by an overall expansion of the envelope. The combination of expansion and cooling provides the condition for only slight variations of the luminosity for some time (Grasberg & Nadezhin 1976; Imshennik & Nadezhin 1989). This emerges as a plateau phase in the light curve. The plateau phase lasts for about 1 days for Model. The light curve shape of Model looks similar to that of ordinary type IIP supernovae. This fact is also noticed by Scannapieco et al. (25) and Kasen et al. (211). The absolute V-band magnitude becomes about M V = 19 mag at the visual maximum, which is about 1 times brighter than an average SN IIP. However, the peak luminosities of type IIP SNe vary by a factor of 1, and that of our model is still contained in this range. Its high luminosity results from the relatively high supernova energy (about 1 foe) and the large radius (3394 R ) of the progenitor. Once the photosphere recedes to the region below the hydrogen-rich envelope, the luminosity decreases rapidly and then becomes governed by the radioactive decay of 56 Co. It is interesting to consider the effective temperature evolution which we present in Figure 7. Note that the recombination effectively slows the adiabatic cooling during the expansion. The recombination starts playing a role when recombination 5

6 Temperature, K Tcolor Teff Tcolor Teff Fig. 7: Evolution of the effective and color temperature of 15 M and 25 M PISN models. The solid and dashed curves present effective and color temperatures correspondingly. The thin line is located at 5 K and approximately corresponds to hydrogen recombination temperature. radiation becomes comparable to adiabatic cooling. This happens when the effective temperature is approximately 1 K (Imshenik & Nadezhin 1965), i.e. well before the establishment of a recombination cooling wave. The effective temperature of an ordinary SN II remains nearly constant ( 55 K) after the establishment of a recombination cooling wave. This phase begins about 2 days after the explosion for an ordinary SN II (Bersten & Hamuy 29; Dessart & Hillier 211). This delay of 2 days depends mostly on the progenitor radius (Grasberg et al. 1971; Imshennik & Utrobin 1977). Due to the very large progenitor radius of 15 M PISN (34 R ) the onset of the recombination cooling wave is delayed up to day 1 in this case. For about 2 days thereafter, when the photosphere is still inside the hydrogen-helium layer (Figure 5), the recombination cooling wave is established for a while and the effective temperature is kept at the hydrogen recombination level. Although the mass fraction ratio between hydrogen and helium is 4:1, the number of hydrogen atoms is about equal to the number of helium atoms and the relative contribution to the electron density from hydrogen stays high. The plateau phase of this model mostly corresponds to the phase of rapidly evolving photospheric temperature. This is different from the case of ordinary SN IIP (see e.g. Bersten & Hamuy 29). In Section 4.3 we provide the comparison of the color temperature with those of typical and bright SNe IIP, and with SN 29kf which was bright in the near-ultraviolet (NUV) range at early time The 25 M model The light curve looks very different for our Model (Figure 4). After the shock breakout, the bolometric magnitude drops to about M 18 mag during the first 5 days. The luminosity decreases due to the adiabatic expansion more rapidly than in Model, because the progenitor of Model is more compact (i.e., R = 745 R, compared to 3394 R in Model ). A precursor-like event happens at around day 2 in U, B, V-bands at magnitude M 18.5 mag. Such precursor could be perceived as a separate hydrogen-rich SN being discovered long before reaching the maximum luminosity. A rebrightening occurs thereafter, as the energy from the radioactive decay of nickel and cobalt in the ejecta diffuses out (Dessart et al. 21). The photosphere recedes below the hydrogen-rich envelope starting at about day 175 (Figure 5), but unlike in the case of Model, the luminosity keeps increasing until about 22 days, reaching M V = 21 mag at the visual maximum. At later times, the light curve is largely governed by the radioactive decay of 56 Co. A hump-like feature on top of the main maximum phase is shown around day 22. At this time the photosphere leaves the oxygen-rich layer and moves down to the silicon-rich layers. The receding front gradually encounters a bubble of diffusing photons generated by the radioactive decay. It is interesting to note that the photosphere is located in the inner layers where there is no hydrogen and helium well before the luminosity reaches its peak, as mentioned above (Figures 4 and 5). Based on our results we can not predict whether hydrogen and helium lines are expected in the spectrum at the time of maximum luminosity due to excitation by radiation above the photosphere. However, Dessart & Hillier (211) showed that the Balmer lines disappear after the photosphere leaves the hydrogen-rich shell of ejecta in the case of type II plateau SN and there are no Balmer-continuum photons at later time. The spectral models for PISNe by Dessart et al. (213) also indicate that hydrogen lines do not appear once the photosphere moves down to the hydrogen-free core. Non-thermal excitation of hydrogen can happen if there is some degree of mixing of cobalt and nickel into layers containing hydrogen (Li et al. 212). However it is expected that PISNe do not experience such mixing during the explosion according to the recent numerical simulations of Joggerst & Whalen (211). Therefore, we expect that PISNe like Model will not have any hydrogen lines once the photosphere recedes to the hydrogen-free core. For the same reason, we do not expect helium lines in the optical bands either, because strong mixing of helium and radioactive nickel is needed to excite helium lines as shown in Dessart et al. (212a). This means that the Model would appear as a type II SN initially but look like a SN Ic from about day 175, well before the luminosity reaches its peak value. 4. Discussion 4.1. Comparison with other theoretical PISN light curves In this section we compare our results with other PISN models from red or blue supergiant progenitors. These include PISN models carried out with the Monte Carlo radiation transport code SEDONA (labeled.k in figures, tables; Kasen et al. 26, 211) and with the non-lte radiative-transfer code CMFGEN (labeled.d ; Dessart et al. 213). Light curves from the SEDONA simulations are integrated from spectra provided by D. Kasen and artificially smoothed. The noisy appearance is a consequence of the statistical approach of Monte Carlo transport calculations. The advantage of the STELLA calculations over the SEDONA and the non-lte CMFGEN codes is the capability to simulate the radiation-hydrodynamical evolution from the onset of the explosion. nce we reproduce the emergence of the shock breakout event and the hydrodynamical evolution during the pair instability explosion. Both, SEDONA and CMFGEN, carry out the radiative transfer calculations once the homologous expansion is reached. The codes have their advantages 6

7 45 46 log Bolometric luminosity R19.D R25.K R15.K log Bolometric luminosity B21.D B25.K B19.D Time, Days Fig. 8: Bolometric luminosity for the red supergiant models (blue curve), R15.K (green curve), R19.D (red curve), and R25.K (cyan curve) Time, Days Fig. 1: Bolometric luminosity for the Model (blue curve), Model B19.D (green curve), Model B21.D (red curve), and Model B25.K (cyan curve). Photospheric Velocity, 1 km/s, R15.K R25.K R19.D Time, Days Fig. 9: The photospheric velocities for the same models as shown in Figure 8. Labels and colours have identical meanings. for spectral computations well described in a number of papers (Kasen et al. 26; Woosley et al. 27; Dessart & Hillier 21, 211) Red supergiants Red supergiant (RSG) progenitors are common for producing SNe IIP (Grasberg et al. 1971; Smartt 29). Typically the light curve from the explosion of a red supergiant is characterized by a pronounced plateau phase lasting for about 1 days, followed by a radioactive tail (Barbon et al. 1979). In this section we compare our lower mass PISN Model with other PISNe produced by RSG progenitors. In Figure 8, we plot the theoretical light curves for different RSG PISN models: Model, Model R15.K, Model R15.D, and Model R25.K (see Table 1 for details). As already discussed by Kasen et al. (211) and Dessart et al. (213), more massive progenitors result in higher luminosities and broader light curves. In particular, heating due to radioactive decay of nickel and cobalt leads to a rebrightening in Model R19.D and Model R25.K which produced large amounts of 56 Ni (2.63 M and 37.9 M, respectively), while it is not seen in Model R15.K and Model which have only.7 M and.4 M of 56 Ni, respectively. The tail of each light curve is powered by cobalt decay and its luminosity is directly proportional to the amount of radioactive nickel generated during the explosion. This dependence is clearly visible in Figure 8. Our Model is directly comparable to Model R15.K, because it has the same initial progenitor mass, a similar explosion energy and a similar nickel mass (Table 1). Compared to Model, Model R15.K has a 5 times lower metallicity, which results in several important differences in the progenitor properties: 1. The ejecta mass is much lower for Model (94 M ) than that of Model R15.K (143 M ). 2. The hydrogen-helium envelope masses are 29 M in Model and 71 M in Model R15.K, respectively. 3. Because of stronger mass loss, the hydrogen mass fraction in the envelope of Model is much lower (X H.2) than that of Model R15.K (X H.7). 4. The radius at the pre-supernova stage is significantly larger for Model (R = 3394 R ) than for Model R15.K (R=2314 R ). As a consequence the overall plateau duration is shorter in Model than in Model R15.K, because of the smaller hydrogen envelope mass and lower mass fraction of hydrogen in the envelope. Moreover, nickel was additionally mixed into outer layers in Model R15.K which in turn caused extra nickel heating during plateau phase. Figure 9 plots the photospheric velocities of the RSG models. The Model and Model R15.K have a phase where the photospheric velocity remains nearly constant, which results from the interplay between the recession of the photospheric front and the expansion of the ejecta. For the models with higher E/M ratio, i.e. Model R19.D and Model R25.K, the photo- 7

8 Photospheric Velocity, 1 km/s, B19.D B21.D B25.K Time, Days Fig. 11: The photospheric velocities for the same models as shown in Figure 1. Labels and colours have identical meanings. spheric velocities stay at a high level for a longer time. This is the result of a stronger explosion and hence a more powerful shock. More energetic photons keep the medium ionised for longer and the photosphere remains in hydrogen-rich regions (larger Lagrangian mass coordinate) for a longer time (up to around day 25 and day 3 for Model R19.D and Model R25.K, respectively). The photospheric velocity in Model is very low at the moment of shock breakout (4 km/s) because the shock spends a large fraction of its energy on ionizing the medium while traveling along the extremely extended envelope. As showed in Grasberg et al. (1971), Falk & Arnett (1977), Blinnikov & Tolstov (211) and Tolstov et al. (213), extended progenitors result in lower photospheric velocities than more compact ones Yellow and blue supergiants In Figure 1, we plot light curves for PISNe from more compact progenitors than RSGs: Model, Model B19.D 2, Model B21.D, and Model B25.K. These PISNe produce light curves reminiscent to that of SN 1987A: the luminosity decreases rapidly in the beginning, but the supernova rebrightenes as the thermalized photons from the radioactive decay of nickel and cobalt diffuse out. The bolometric luminosity of our yellow supergiant Model is higher than that of the other blue supergiant models during the initial phase (t 5 days), because of the larger radius of the progenitor (e.g. Popov 1993). In Figure 11, we show the evolution of the photospheric velocity for the compact progenitor models. In all of these models, the photosphere recedes very rapidly through the outer hydrogen-helium layer after shock breakout. In Model, Model B21.D and Model B25.K, where the amounts of nickel exceed 19 M, the reversion of the photospheric velocity occurs when the recombination and cooling wave encounters the expanding bubble of diffusing photons generated by nickel 2 Blue supergiant B19.D is the same evolutionary model as red supergiant R19.D but with truncated hydrogen atmosphere (about 3 M of outer hydrogen layer is cut; Dessart et al. 213). radioactive decay. These photons ionize the just recombined medium and push the photospheric front to outer shells. A larger envelope mass leads to a broader phase of this reverse photosphere motion. The sudden drop in the velocity of Model occurs around day 175. At this time the photosphere leaves the hydrogen-helium envelope (Figure 5) and moves rapidly through hotter layers of oxygen, neon, carbon and magnesium heated by diffusing nickel photons. There is a lack of such a sharp drop in other calculations because their input chemical structure was smeared (Kasen et al. 211; Dessart et al. 213) to mimic hydrodynamical mixing happened during the explosion. However, it was shown by a number of studies that the degree of mixing for the inner regions containing radioactive nickel is not so prominent (Joggerst & Whalen 211; Chen et al. 212; Chatzopoulos et al. 213; Chen et al. 214). Mixing is more efficient in the oxygen layer where the shock emerges and above it due to the propagation of the reverse shock. Red supergiants exhibit a higher degree of mixing compared to more compact blue supergiants, similar to red and blue supergiant core collapse SNe (Joggerst et al. 29) The chemical structure during the coasting phase Figure 12 shows the ejecta structure of Model and Model at coasting phase. In fact, the coasting phase begins around day 1 for both our models. In Figure 12, we plot the chemical composition of the ejecta at day 95 and day 12, respectively. The plots demonstrate what would be the degree of Doppler broadening for the spectral lines of the given elements 1 days after the explosion and later. The oxygen shell in Model, which contains small amounts of carbon, neon, magnesium and silicon, expands at a velocity of 1 22 km s 1. The inner silicon-sulfur shell moves at velocities below 2 km s 1. Model shows systematically higher velocities. The oxygen-rich layers travel at 4 5 km s 1, the silicon-sulfur shell at about 4 km s 1, and the nickel-rich layer moves at 39 km s 1. In the Model the structure has some peculiarity at the sulfur-silicon layer moving at velocity around 4 km s 1. The reason is the hydrodynamical effect of nickel heating (Arnett 1988; Woosley 1988; Arnett et al. 1989; Pinto & Eastman 2; Woosley et al. 27). The additional radiation pressure produced by nickel heating provides an acceleration of the innermost 3 M of ejecta. The effect increases the velocity field in the inner part of the ejecta by a few hundred km s 1 and decreases the density. The overall changes in density and velocity do not exceed 1% on a relative scale. This effect has no influence on the light curve shape around the luminosity peak (Woosley et al. 27) because the photosphere retains far from the innermost region at this time Comparison with observed SNe Relatively low-mass PISNe and type IIP SNe Our Model is particularly relevant for the identification of PISNe in the local Universe. This is because SNe from the low-mass end of the PISN regime like Model are expected to be more abundant than those from the high-mass end like our Model. In Figure 13, therefore, we show the first 17 days of light curves for our synthetic PISN models and compare it to several observed plateau SNe. The data for these particular SNe are taken from the Sternberg Astronomical 8

9 Table 2: Shock breakout and plateau-maximum phase characteristics. Shock breakout duration (defined as a full width at halfmaximum), bolometric peak luminosity, effective and color temperature, spectral wavelength peak in Å, spectral energy peak in kev. The photospheric phase begins after relaxation from the shock breakout and is limited by the transition to the radioactive decay tail. Model and Model are those of 15 M and 25 M PISNe simulated in the frame of present study. Labeled R19.D, B19.D and B21.D are models from (Dessart et al. 213) and R15.K and R25.K are models simulated by (Kasen et al. 211). shock breakout photospheric phase duration, L, T eff /T color, λ max, E p, duration, L, (hours) (erg/s) (1 3 K) (Å) (kev) (days) (erg/s) / R15.K / R19.D R25.K / B19.D B21.D / B25.K / Mass fraction Mass fraction Si S O Ar Ca Fe Co Ni Ne C Mg Velocity, 1 km/s Fe Co Ni Ar Ca Si Velocity, 1 km/s O S Ne Mg C Fig. 12: The chemical structure of ejecta of 15 M and 25 M PISNe at day 95 and day 12 correspondingly. H H Absolute V magnitude H 29kf 1999em 1992am Fig. 13: The absolute V-band light curve for 15 M PISN (blue solid line) is shown together with the absolute V-band magnitudes for typical plateau SN 1999em (green circles) and three bright plateau SNe 1992H (magenta circles), 1992am (cyan squares), 29kf (red squares). Institute Supernova Light Curve Catalogue (Tsvetkov et al. 21; Tsvetkov & Pavlyuk 213) and are compiled from the original data (Tsvetkov 1994; Clocchiatti et al. 1996; Hamuy 21; Botticella et al. 21). Model has a plateau phase during the first 115 days which is not unusual for SNe IIP. Although this model has a higher envelope mass (29 M ) than those of typical SNe IIP progenitors (< 1 M ), the helium mass fraction is very high (about 8%). At the same time it is expected that the envelopes of typical SN IIP progenitors have helium mass fractions of 35% 5% depending on their initial masses (Woosley et al. 22; Langer 212). lium recombines at higher temperatures to reduce the electron scattering opacity and the hydrogen recombination front recedes more rapidly with a higher fraction of helium in the envelope (Kasen & Woosley 29). This makes our Model to have a fairly short duration of the plateau phase compared 9

10 12 Color Temperature, 1 4 K em 29kf Fig. 14: Color temperature evolution of 15 M PISN is shown together with those of typical SN IIP SN 1999em and NUV-bright SN 29kf. SN 1999em and SN 29kf data are taken from Bersten & Hamuy (29); Botticella et al. (21). SN 29kf data are shifted in time by 15 days. Photospheric Velocity, 1 km/s, SN IIP: 1999br 1999em 1999gi 23bn 23hn Fig. 15: The photospheric velocity for our 15 M PISN model (blue solid line) is shown together with observational data for several SNe IIP (Jones 28). The zero point for the observed data is the first spectroscopic observation. to high redshift RSG PISN models (cf. Figure 8). The plateau duration of the light curve for our Model is comparable to those of typical SNe IIP, despite the relatively high envelope mass. During the plateau phase, the V-band magnitude varies by 1-2 magnitudes (Barbon et al. 1979; Hamuy 23). Compared to the typical plateau supernova SN 1999em, the plateau luminosity of the Model is higher by 2-3 magnitudes, but comparable to those of the three bright SNe IIP (SN 1992H, 1992am and 29kf, cf. Figure 13). The estimated nickel masses for these SNe IIP are comparable or higher than that in Model : M Ni =.58 M,.75M,.36M for 1999em, 1992H, 1992am (Clocchiatti et al. 1996; Hamuy 23; Nadyozhin 23; Bersten et al. 211), respectively, and limited by.4 M for 29kf (Botticella et al. 21). Our Model has a very large initial radius of 3394 R (see Table 1), while those of ordinary SNe IIP progenitors have radii of less than 1 R. Due to this difference our PISN explosion and SN IIP explosions have a different appearance of the shock breakout event. The shock breakout duration is longer for a larger progenitor (Tominaga et al. 211). The color temperature is higher for a SN IIP shock breakout than that for our PISN (T col R 1/2 ) while the peak luminosity is similar depending mostly on the explosion energy (Tolstov et al. 213). Since our PISN shock breakout is redder the spectral maximum occurs at a longer wavelength (see Table 2). Unfortunately the shock breakout of local SNe II is very difficult to detect because it appears as a ultraviolet/x-ray burst lasting only minutes to hours (Calzavara & Matzner 24). The large progenitor radius of our 15 M PISN model has consequences for the photospheric temperature evolution. In Figure 14, we compare the color temperature of our 15 M PISN with those of the typical SN IIP SN 1999em and the near UV-bright SN 29kf. The color temperature evolution of our model is very different to that of SN IIP 1999em, but similar to that of SN IIP SN 29kf, which had a high NUV excess at early time. Nevertheless, the NUV light curve of our model does not show such a high luminosity as SN 29kf. The maximum NUV luminosity reaches 2.5 mag and 22 mag for our 15 M PISN model and for SN 29kf, respectively (Botticella et al. 21). The high NUV luminosity of 29kf is explained by the interaction of the SN shock with a dense stellar wind preceding the SN explosion (Moriya et al. 211). The higher temperature is the direct consequence of the ultraviolet excess. However, in case of our low mass PISN model the shock breakout event occurs and the high temperature is related to the relaxation of the highly excited medium. Any interaction of the SN ejecta with the progenitor wind is neglected in our model. Another possible way to distinguish a SN IIP from a PISN explosion is to check the photospheric velocities. In Figure 15 we show the photospheric velocity for our 15 M PISN model along with those of a few ordinary plateau SNe taken from Jones (28). The estimate for the photospheric velocity is based on the H β absorption line (Jones et al. 29). The distinct property of our model is the low photospheric velocity at earlier time compared to the maximum photospheric velocities of SNe IIP (cf., Grasberg et al. 1971; Young 24). The photospheric velocities at later time are similar for both SNe IIP and our low mass PISN model SLSNe linked to type Ic SNe As discussed in Section 3 the photosphere of our 25 M PISN model during its peak luminosity phase recedes to the bottom of the hydrogen-rich envelope (see also Kasen et al. 211; Dessart et al. 213). Therefore, one could classify it as a SN Ic, like SN 27bi, if it was discovered during the maximum or post-maximum phase. In Figure 16, therefore, we compare the synthetic bolometric light curve of our 25 M model with the quasi-bolometric light curves of some superluminous type Ic SNe: SN 21gx, PS1-1awh, PS1-1ky PS1-1bzj, SN 27bi, and PTF1nmn. From Figure 16 (left panel) it is clear that Model can not explain several of the unusually luminous SNe Ic, because of the very broad light curve of our model. Unambiguous evidence for a pair instability explosion of an initially high mass 1

11 log Bolometric Luminosity SN21gx slsn Ic PS1 1awh slsn Ic PS1 1ky slsn Ic PS1 1bzj slsn Ic SN27bi slsn Ic Absolute magnitude Mbol PTF1nmn R band R band SN27bi R band Fig. 16: Left panel: The absolute bolometric light curve for our 25 M PISN model (blue solid curve) is shown together with quasi-bolometric and bolometric light curves of superluminous type Ic SN 21gx (red, Inserra et al. 213), PS1-1awh and PS1-1ky (green and black, Chomiuk et al. 211) PS1-1bzj (magenta, Lunnan et al. 213), and PISN candidate SN 27bi (yellow, Young et al. 21). Right panel: The absolute R-band and bolometric light curve for our 25 M PISN model (green and blue) with superimposed R-band data for SLSN 27bi (cyan) and PTF-1nmn (red). Data are taken from Gal-Yam (212b). For both plots the observed curves are shifted in time by 2 days to coincide with the maximum phase of the theoretical curve. progenitor would be the observation of an extremely long rise to the maximum phase (Benetti et al. 213). Therefore, some attempts are made to complete the light curves with retrospective detection of data points before the maximum phase(nicholl et al. 213). The synthetic light curve of our 25 M PISN model shows 2 days of rise while many SNe Ibc and SLSNe demonstrate significantly shorter rise time less than 4 days (see e.g. Modjaz et al. 29; Drout et al. 211; Chomiuk et al. 211). A good example of a SLSN with a reliable long-lasting rise is PTF1nmn. The data for this particular SLSN and for SN 27bi are shown in the right panel of Figure 16. Our Model agrees well with the broad light curves of PTF1nmn and SN 27bi. The photospheric velocity of our high mass PISN model around the luminosity peak is smaller than typical velocities of luminous SNe Ibc (Figure 17). The reason for the low velocity is the high ejecta mass of our model. This may be another criterion for distinguishing luminous PISNe powered by large amounts of nickel from superluminous CCSNe. However, the photospheric velocities measured for SN 27bi demonstrate that the ejecta of this particular SLSN moves at a lower velocity than those of other SLSNe, which is in good agreement with our model. This renders precise spectroscopic observations important to shed light on this question. Being discovered around or after its maximum, SN 27bi resembles some other SLSNe which show a short rise to their peak luminosity (Nicholl et al. 213). Particularly, this rules out the pair instability origin of these SNe. Nevertheless, we conclude that our Model agrees with observed properties (light curve, photospheric velocity) of SLSN 27bi well. Therefore, SN 27bi might emerge from pair instability explosion of very massive star with initial mass above 2 M. 5. Conclusions We carried out simulations of shock breakouts and light curves of pair instability supernovae using two evolutionary models of 15 M and 25 M at metallicity Z = 1 3 (Langer et al. 27; Kozyreva et al. 214). We used the radiation hydrodynamics code STELLA for this purpose (Blinnikov et al. 26). The considered metallicity (Z= 1 3 ) is among the highest of PISN models that have been so far presented in the literature (rzig et al. 199; Whalen et al. 213b). Therefore, our models may serve as useful references for future studies on PISNe observed in the local Universe, as well as in the early Universe. From our qualitative comparison to ordinary core collapse SNe we conclude that it is difficult to distinguish low mass pair instability explosions from hydrogen-rich core collapse explosions. The photometric and spectroscopic observations, including X-ray and ultraviolet (for detection of shock breakout events), should be very detailed from the earliest epoch to help shedding light on this. The increasing number of SN surveys allows to increase the number of discovered SNe and detailed data from the very early epoch after explosion, especially those missions which have the short cadences (e.g. PTF). Given the low-mass preference of the stellar initial mass function, a large fraction of PISNe that will be observed in the local Universe could resemble our 15 M model, which represents PISNe from the low-mass end of the PISN regime. These PISNe are predicted to have the following characteristics: 1. The progenitors are likely to be red-supergiants having very extended envelopes (R 3 R ), if they can retain some fraction of the hydrogen envelopes by the time of explosion. Our 15 M model has the final mass of 94 M and the envelope mass is 29 M, which is significantly smaller than in the corresponding case of zero or extremely low metallicity ( 7 M ). The hydrogen mass fraction in the envelope is only about