Free Access
Issue
A&A
Volume 625, May 2019
Article Number A9
Number of page(s) 7
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201834732
Published online 30 April 2019

© ESO 2019

1. Introduction

Understanding how the landscape of type-II-Plateau (II-P) supernova (SN) properties connects to the diversity of red-supergiant (RSG) star progenitors and their explosion is of great interest for astrophysics. Unfortunately discrepant estimates of the ejecta and progenitor masses are often obtained from radiation-hydrodynamics simulations that study the bolometric light-curve evolution1, from pre-explosion imaging, and from studies of host stellar populations.

While significant advances have been made in the modeling of the proto-neutron star phase leading to shock revival and explosion in massive star progenitors, there are still many unresolved issues. Although the ejecta properties in these studies are broadly consistent with inferred properties from type-II-P SN observations (Lentz et al. 2015; Müller et al. 2017; Glas et al. 2019; O’Connor & Couch 2018; Vartanyan et al. 2019), the simulations cannot predict which stars produce core-collapse SNe, or quantities such as the 56Ni mass and explosion energy. Complications also arise from uncertainties in the structure of the SN progenitor prior to core collapse. These theoretical studies would benefit from reliable observational inferences about SN and progenitor properties.

Using pre-explosion photometry, observations have allowed the detection of a progenitor RSG for a handful of type-II-P SNe (Van Dyk et al. 2003, 2012a,b; Smartt et al. 2004; Maund et al. 2005, 2013, 2014; Fraser et al. 2012, 2014; O’Neill et al. 2019). These works propose main sequence progenitor masses below about 17 M (and the lack of more massive RSG progenitors has been coined the “RSG” problem), although these inferences are often uncertain. The problem of inferring the progenitor luminosity is made difficult for example because of the frequent lack of multi-band photometry, the uncertain bolometric correction (Davies & Beasor 2018), or the uncertain reddening (Beasor & Davies 2016). Furthermore, these inferences are based on stellar evolution models, whose predictions depend on physical properties that are either not known a priori (e.g., initial rotation) or are hard to model (e.g., convection, overshooting; Arnett et al. 2015). These can impact the He-core mass and the progenitor luminosity at the time of explosion (Heger et al. 2000; Meynet & Maeder 2000; Hirschi et al. 2004).

Radiation hydrodynamics and radiative transfer can also be used to characterize the properties of type-II-P SNe. The most robust inference and the most easily done is for the 56Ni mass since full γ-ray trapping holds at the onset of the nebular phase. At that time, the bolometric luminosity equals the total power radiated from 56Ni and 56Co decay, and the mass of ejected 56Ni can be determined at an accuracy that depends only on the accuracy of the adopted distance, reddening, and amount of flux falling outside of the observed range (e.g., Woosley 1988; Hamuy 2003). Other inferences such as explosion energy, ejecta mass, and progenitor mass at collapse or on the main sequence require some modeling and are therefore more uncertain.

Radiation hydrodynamic simulations of type-II-P SNe have been performed for several decades. Early works identified the basic properties of the progenitors, in particular their large progenitor radii and massive H-rich envelopes compatible with RSG progenitors (Grassberg et al. 1971; Falk & Arnett 1977). Basic relations between ejecta and observed type-II-P SN properties were also drawn (Litvinova & Nadezhin 1985; Popov 1993). When applied to observed type-II-P SNe, these relations can however produce very perplexing numbers, such as ejecta masses of several 10 M or small supergiant radii atypical of RSG stars (see, e.g., Hamuy 2003). In many cases, these quantities conflict with predictions from stellar evolution and progenitor observations.

Studies dedicated to specific objects allow for a more refined analysis through the production of a tailored model (see, e.g., Turatto et al. 1998; Utrobin 2007; Chugai & Utrobin 2000; Bersten et al. 2011; Dessart et al. 2013; Lisakov et al. 2018; Morozova et al. 2018). But in this case, one notes a large disparity in inferred ejecta properties, in particular the ejecta mass, which then translates into an even larger disparity in the progenitor main sequence mass. These masses are often much larger than those inferred for the progenitor stars from pre-explosion imaging (Smartt 2009), or from the modeling of nebular-phase spectra (Maguire et al. 2012; Jerkstrand et al. 2012, 2014). Since mass is the most defining characteristic of a star, this discrepancy is a problem.

Type-II-Plateau SN light curves are however primarily sensitive to the H-rich envelope mass and not to the mass of the progenitor He-core (Dessart & Hillier 2011a; Dessart et al. 2013). Therefore, in many instances the inferred “ejecta” mass refers only to the H-rich envelope mass of the progenitor star at the time of explosion. This ambiguity arises because the words envelope mass and ejecta mass are used interchangeably in the literature, while the two words refer to two different entities. Stellar evolution indicates that their value may differ by a factor of ten.

To clarify this property, we use the public grid of massive star models from Woosley et al. (2002). Figure 1 shows a montage of density profiles at core collapse for massive stars evolved at solar metallicity with zero-age main sequence masses (MZAMS) of 11 up to 30 M. The left column shows the standard way of presenting such profiles, with the origin at the center of the star. As the initial mass MZAMS increases, the mass of the He core (MHe, c; shown in red, and corresponding to regions with a density greater than about 1 g cm−3) increases from about 3 M for the 11 M model up to 10 M for the 30 M model. It is in this hot and dense He core that nuclear fusion takes place. In the lower-density lower-temperature H-rich envelope, the composition deviates from its primordial value only through the effect of mixing (see, e.g., Davies & Dessart 2019). The mass of the H-rich envelope (MH, e) at core collapse is between 8 and 9 M for stars with MZAMS between 11 and 21 M, and decreases (not necessarily monotonically) as MZAMS increases further. This arises from an empirically based formulation of RSG mass-loss rates, which reflects the greater mass-loss rates inferred for more massive and more luminous RSGs (de Jager et al. 1988). This influence of mass loss on MH, e is what eventually turns more massive stars into H-deficient Wolf-Rayet stars (Maeder & Meynet 1987; Langer et al. 1994; Crowther 2007). While RSG star mass loss remains uncertain, this should only affect the value of MZAMS where MH, e starts decreasing with increasing MZAMS, and the actual value of MH, e for a given MZAMS2. The trend shown in Fig. 1 must hold.

thumbnail Fig. 1.

Density profiles vs. Lagrangian mass at the time of core collapse for a set of stellar models with an initial mass on the main sequence between 11 and 30 M (Woosley et al. 2002). The origin is at star center (left column) or at the He core edge (right column). The H-rich envelope is drawn with a thick black line and corresponds in all models to (extended) regions with a density below about 10−7 g cm−3. The He core is drawn with a thick red line and corresponds in all models to (compact) regions with a density above about 1 g cm−3. The vertical scale is squeezed but is suitable to reveal the evolution of the He core mass and the H-rich envelope mass for different ZAMS masses. Left panel: smooth progression in the He core mass with initial progenitor mass, and the varying progenitor mass at the onset of core collapse. Right panel: nonlinear evolution of the H-rich envelope mass with initial progenitor mass is emphasized. In these models, stars with ZAMS masses of 11–21 M contain similar H-rich envelope masses (8–9 M) at the onset of core collapse.

Because type-II-P SN light curves are primarily sensitive to the H-rich envelope mass, it is more instructive to show the density structure with respect to the edge of the He core (right column of Fig. 1). Ignoring the He core, the models with MZAMS below 21 M have similar density structures, while above the models have progressively lower values of MH, e. For a given ejecta kinetic energy Ekin imparted to the H-rich envelope, models with MZAMS below 21 M should have comparable light curves (modulo differences in R), while above 21 M, the light curve should be characterized by a faster post-breakout luminosity decline and a shorter photospheric phase (Bartunov & Blinnikov 1992). Furthermore, “ejecta” masses (i.e., MH, e values) inferred from type-II-P SN light curves should consequently be around 10 M since the RSG stars from 10 − 20 M stars on the main sequence are favored by the initial mass function. Masses of just a few M (which result from higher-mass progenitors) should be less frequently inferred (the reasoning here ignores binarity). This basic picture seems compatible with the observed diversity of type-II SN light curves (Anderson et al. 2014; Sanders et al. 2015). Furthermore, the large range of type-II SN ejecta masses (extending to large masses) published in the literature appears incongruous.

The goal of this paper is to demonstrate that the type-II-P SN ejecta masses inferred from light-curve modeling bear considerable uncertainty. One cannot infer the progenitor mass (either at core collapse or on the main sequence) using light-curve modeling. The problem is degeneracy since stars with widely different total masses may have the same MH, e at core collapse. To illustrate this property, we conduct a controlled experiment on a set of massive star models with different main sequence mass but the same H-rich envelope mass – our model set is very similar to that of Woosley et al. (2002) and thus representative of predictions for single star evolution. The next section presents the numerical approach. Section 3 then presents the SN radiative properties for each model and discusses the implications of the results.

2. Numerical approach

The progenitor models correspond to solar-metallicity nonrotating stars of 12, 15, 20, and 25 M on the main sequence. These were evolved with MESA version 7623 (Paxton et al. 2011, 2013, 2015) as part of earlier works (Lisakov et al. 2017, 2018). These stars are evolved as single stars, or in binary systems with a wide orbit: the possibility of mass exchange with a companion star is ignored. These simulations use standard MESA default parameters, except for a mixing length parameter of three. This choice is needed to produce more compact RSG stars at core collapse (see Dessart et al. 2013, and more recently Paxton et al. 2018). Because this produces RSG star models with a higher effective temperature, a slight enhancement on the default mass-loss rate parameter is needed in order to yield a mass-loss rate similar to that of RSGs with a lower effective temperature. The model properties are summarized in Table 1, while Fig. 2 shows the density structure at the onset of core collapse for each model. The models differ in He-core mass but have a similar H-rich envelope mass – their properties are similar to those of Woosley et al. (2002) for the corresponding ZAMS mass (Fig. 1). In radial space, the differences in He core properties are invisible since the He core has the size of the Sun while the surface radii cover approximately from 400 to 900 R. Because they make no distinction between the high-density He core and the low-density H-rich envelope, polytropic progenitor structures are unsuitable for light-curve modeling. To comply with the early-time observations of some type-II-P SNe (see, e.g., Yaron et al. 2017), about ∼0.2 M of atmospheric (static) material is added above R using a density scale height of about 0.2 R. This has however no relevance for the present discussion.

Table 1.

Properties of the model set.

thumbnail Fig. 2.

Density structure of the models m12, m15, m20, and m25 shown vs. Lagrangian mass (the origin is the outer edge of the He core). The differences in H-rich envelope density reflect the differences in surface radii (see Table 1).

The explosion is modeled with the 1D gray radiation-hydrodynamics code V1D in the usual manner (Livne 1993; Dessart et al. 2010). The explosion is triggered by means of a piston placed within the Si-rich shell, at a mass cut of 1.5–1.94 M. To compensate for the progenitor binding energy, the explosion energy is adjusted for each model to yield an asymptotic ejecta kinetic energy Ekin of ∼1.25 × 1051 erg. Although V1D treats explosive nucleosynthesis, the 56Ni mass is reset to be 0.04 M at 100 s after the piston trigger (the value in Table 1 deviates a little because of the remapping into CMFGEN). At 14 d after the piston trigger (when we remap into CMFGEN), the four models have essentially the same density structure in velocity space all the way down to the He-core material (whose outer location corresponds to the density jump; Fig. 3). The jump is at a velocity of 1000 − 2000 km s−1 greater in model m25 than in model m12, which reflects the differences in He-core mass between the models (see Dessart et al. 2010 for discussion). The temperature above that jump is lower in more compact progenitor models. The material that used to be in the He core has cooled considerably because it expanded by a factor of about 1000, while the H-rich envelope has expanded by a factor of about ten. This conspires to make the temperature very uniform throughout most of the ejecta. The relatively small volume occupied by the former He core stores a modest radiative energy. Instead, the large volume occupied by the shocked H-rich envelope stores most of the radiative energy that will be released during the high-brightness phase of a type-II-P SN. It also stores the bulk of the kinetic energy.

thumbnail Fig. 3.

Ejecta velocity, fractional kinetic energy Ek, frac, density, and temperature vs. mass at 14 d after the piston trigger. Ek, frac is defined as ∫V2dm/2Ekin, with the integration done inward. The x-axis origin is at the outermost ejecta layer in order to emphasize how similar the properties of the shocked H-rich envelope are in these different progenitor models.

At 14 d, the ejecta models computed by V1D are remapped into the nonlocal thermodynamic equilibrium (nonLTE) and time-dependent radiative-transfer code CMFGEN (Hillier & Dessart 2012). The evolution is then followed until 300 d using the standard procedure (see, e.g., Dessart et al. 2013).

There is at present no robust theory that predicts with certainty the ejecta kinetic energy and 56Ni mass that a given ZAMS mass should produce. However, the simulations of Sukhbold et al. (2016) suggest that stars in a broad mass range from about 12 to 25 M can yield an ejecta within 50% of 1.2 × 1051 erg and a 56Ni mass within a factor of two of 0.04 M (their is some scatter and some offset with MZAMS depending on the calibration used for their “thermal bomb”). Our choice is therefore realistic.

3. Results and discussion

Figure 4 summarizes the main results from the CMFGEN simulations for models m12, m15, m20, and m25. The bolometric light curves are similar. Models with larger progenitor radii are more luminous, as expected (see., e.g., Popov 1993; Kasen & Woosley 2009; Dessart et al. 2013). The V-band light curve and the V − I color evolution are also similar, with the bigger progenitors reddening later and being visually brighter (this is the same as for models s15 and s25 in Dessart & Hillier 2011a and for models m15mlt1 and m15mlt3 in Dessart et al. 2013). These color offsets are understood from the offset in photospheric temperature (also visible in the H ionization state), with the bigger progenitors remaining hotter at the photosphere for longer. Despite the large differences in progenitor and ejecta masses, the four models have a similar photospheric-phase duration because they have a similar H-rich envelope mass (Fig. 2). The models show the same qualitative and quantitative evolution in photospheric velocity, because the bulk of the kinetic energy is held up in what used to be the H-rich envelope. Throughout the photospheric phase, the H mass fraction at the photosphere is comparable in all models (i.e., the photosphere resides in the former H-rich envelope; Dessart & Hillier 2011a) which implies that the spectra should reveal lines from similar elements.

thumbnail Fig. 4.

Summary of results from the CMFGEN simulations of models m12, m15, m20, and m25, showing, from top to bottom, the bolometric light curve, the V-band light curve, the V − I color evolution, and the evolution of some photospheric properties (T3, ph: temperature in units of 1000 K; V3, ph: velocity in units of 1000 km s−1; XH, ph: hydrogen mass fraction; IH, ph: hydrogen ionization state).

While the light curves and colors show some differences, these are not sufficient to unambiguously determine the ejecta mass. The main driver behind the differences in photometric properties is the progenitor radius, not the ejecta mass. If we had used a different mixing length parameter for the four progenitor models (smaller in lower mass models), we could have produced progenitors with the same radius, and in that case the light curves would have been even more similar. Furthermore, because RSG observations probe the outer low-density fluff of the RSG atmosphere rather than the higher-density hydrostatic photosphere (as defined in stellar evolution models), surface radii are poorly constrained. Further, and contrary to a generally held belief, RSG radii are not predicted with any confidence by stellar-evolution models since R is controlled by convection and a prescribed mixing length (Maeder & Meynet 1987; Dessart et al. 2013).

Other factors also influence the light curves. For example, there is increasing evidence that the early-time light curves of many type-II-P SNe are influenced by circumstellar (CSM) interaction, which affects both the brightness and colors of the SN at early times (e.g., Yaron et al. 2017; Morozova et al. 2017, 2018; Dessart et al. 2017; Moriya et al. 2017; Förster et al. 2018). Clumping may also hasten the photosphere recession and boost the SN luminosity during the high-brightness phase (Dessart et al. 2018). Further, most simulations of core-collapse SNe lead to asymmetric explosions. This may affect both the light-curve evolution and colors. Observations of type-II-P SNe also show variable levels of polarization, and the degree of polarization is time-dependent (Leonard et al. 2006, 2012). The latter may indicate (although not necessarily) that the degree of asymmetry varies with time (Dessart & Hillier 2011b). For example, the core may be much more asymmetric than the massive H-rich envelope. The amount of 56Ni mixing is also known to affect the shape of the light curve (see, e.g., Kasen & Woosley 2009), and there are observational issues related to the distance of the SN, reddening and the inference of the bolometric light curve from a finite number of bands. Finally, we stress that even if one can infer the ejecta mass from the light curve, an extrapolation to the progenitor mass is fraught with biases because RSG mass-loss rates, and hence the progenitor mass-loss history, are poorly known (see, e.g., Fig. 1 in Meynet et al. 2015).

In nature, it is possible that the explosion of progenitors with very different ZAMS masses yields the same asymptotic ejecta kinetic energy (which is what we assumed here for our model set) and therefore naturally contributes to light-curve degeneracy. The power source for the explosion is gravitational contraction of the core and mass accretion onto it (see, e.g., Burrows & Goshy 1993). This process may release more energy in higher-mass progenitors, but this additional energy would be sapped by the greater binding energy of their He core. In the present set of models, to produce an ejecta kinetic energy of 1.25 × 1051 erg, the total energy released by the “thermal bomb” had to be larger in more massive progenitors. Specifically, it was 1.36 (m12), 1.39 (m15), 1.78 (m20), and 1.96 × 1051 erg (m25) – these values depend on the time left before the onset of collapse (e.g., the central density is not the same in all models). The trend nonetheless holds (see, e.g., Woosley et al. 2002).

Polytropes and nonevolutionary progenitor structures, designed with considerable freedom, conflict with the fundamental features of massive stars at death (primarily through the incorrect treatment of their core-halo structure). It is from these simulations that the largest discrepancies in progenitor masses arise (see also Sect. 4.1 in Morozova et al. 2018). Similarly, analytic scalings (e.g., Litvinova & Nadezhin 1985) are unable to yield a reliable ejecta mass since they are insensitive to the He core mass.

An extensive discussion about progenitor masses is provided by Dall’Ora et al. (2014). In order to examine the parameter space they considered both a semi-analytic approach and an approach based on direct hydrodynamical modeling. Their semi-analytic approach yields multiple minima in χ2 when comparing model predictions with observations. While not all solutions gave consistency with the inferred photospheric velocities, these studies do highlight degeneracies in inferring ejecta masses from light curves.

Morozova et al. (2018) have also performed a parameter study of type-II-P light curves. Unlike the study of Dall’Ora et al. (2014), their simulations are based on evolutionary models computed with KEPLER. With some exceptions, the 2D-χ2 plots of Morozova et al. (2018) generally show a broad range of ZAMS masses and radii that are consistent with their light curves. Given the degeneracies, assumptions in both the evolutionary models and in the light-curve modeling will significantly influence the results. Furthermore, Morozova et al. (2018) infer explosion energies and ejecta masses without using any spectral information, hence ignoring fundamental constraints on the expansion rate. For example, their ejecta kinetic energies are lower than standard, and with their adopted CSM mass, their ejecta models most likely fail to produce the broad Doppler-broadened lines seen in type-II SN spectra. Disregarding spectral constraints, light-curve modeling is even more subject to degeneracies.

While more difficult to model, spectra potentially provide another avenue to constrain ejecta and progenitor mass. Figure 5 shows the spectra for our model set at about 20, 62, and 155 d after explosion. At the first epoch, the difference is the greatest and results from the larger radii in more massive progenitors (which cause the brightness and color offset shown in Fig. 4). However, during the recombination phase (second epoch), the spectra are essentially identical. At the nebular epoch, the spectra are also very similar; they exhibit a strong Hα line, but with a stronger O I 6300 Å doublet line in more massive progenitors.

thumbnail Fig. 5.

Comparison of optical spectra for models m12, m15, m20, and m25 at about 20 d (top), 62 d (middle), and 155 d (bottom) after explosion. The color offset at early times results from the different progenitor radii. The greater O I 6300 Å line in more massive progenitors arises from their greater oxygen content.

A more reliable discriminant for progenitor mass may be sought from nebular-phase spectra. At such times, the stark contrast between a 12 and a 25 M progenitor can be revealed from the inspection of emission lines, and in particular O I 6300 Å. For the present models m12 to m25, the oxygen mass increases with the ZAMS mass (it is 0.16, 0.53, 1.33, and 2.33 M), and this conspires to produce a greater O I line strength in more massive progenitors, which is apparent here in the bottom panel of Fig. 5. Careful modeling is however needed to infer an accurate progenitor mass (Maguire et al. 2012; Jerkstrand et al. 2012, 2014).

Our results show that models from a 12, 15, 20, and 25 M star on the main sequence can yield a type-II-P SN with a similar plateau duration and spectral evolution if they have the same H-rich envelope mass and ejecta kinetic energy. Differences arise primarily from offsets in progenitor radius (affecting the early-time brightness and color) and nucleosynthetic yields (affecting the lines from intermediate-mass elements at nebular times). Because of the unconstrained mass contribution from the He core, light-curve modeling cannot yield a robust inference of the ejecta mass, nor of the progenitor mass at the time of explosion. Furthermore, stars with different main sequence mass can die with the same MH, e (which is function of the uncertain cumulative mass loss), meaning that the constraint of MH, e from light-curve modeling cannot be uniquely connected to a MZAMS. Overall, there may well be a “RSG problem” but light-curve modeling will not solve it.

One may argue that this position is pessimistic. For example, from their set of progenitor stars and explosion models, Sukhbold et al. (2016) obtain a trend in asymptotic ejecta Ekin and 56Ni mass, and we may eventually have a robust prediction of these quantities from ab-initio 3D simulations of the explosion mechanism. We may also eventually have a robust theory for convection, overshoot, and RSG mass loss, with an accurate way to constrain core rotation. Similarly, we may eventually be able to perform multi-dimensional radiative transfer in nonLTE, thereby accounting accurately for the treatment of ejecta asymmetry, clumping, chemical segregation, and mixing. When this time comes, we will be able to constrain ejecta and progenitor masses with great confidence. We will have great accuracy in our results and therefore the degeneracies discussed in this paper will be lifted. However, this is not the present situation and even the state of the art remains limited by our understanding of each of the above building blocks.


1

We refer loosely to this approach as “light-curve modeling”. In practice, it usually (but not always) includes the additional information from one spectral line to infer the evolution of the expansion rate at the photosphere.

2

These quantities should also vary with metallicity.

References

  1. Anderson, J. P., González-Gaitán, S., Hamuy, M., et al. 2014, ApJ, 786, 67 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arnett, W. D., Meakin, C., Viallet, M., et al. 2015, ApJ, 809, 30 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bartunov, O. S., & Blinnikov, S. I. 1992, Sov. Astron. Lett., 18, 43 [NASA ADS] [Google Scholar]
  4. Beasor, E. R., & Davies, B. 2016, MNRAS, 463, 1269 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bersten, M. C., Benvenuto, O., & Hamuy, M. 2011, ApJ, 729, 61 [NASA ADS] [CrossRef] [Google Scholar]
  6. Burrows, A., & Goshy, J. 1993, ApJ, 416, L75 [NASA ADS] [CrossRef] [Google Scholar]
  7. Chugai, N. N., & Utrobin, V. P. 2000, A&A, 354, 557 [NASA ADS] [Google Scholar]
  8. Crowther, P. A. 2007, ARA&A, 45, 177 [NASA ADS] [CrossRef] [Google Scholar]
  9. Dall’Ora, M., Botticella, M. T., Pumo, M. L., et al. 2014, ApJ, 787, 139 [NASA ADS] [CrossRef] [Google Scholar]
  10. Davies, B., & Beasor, E. R. 2018, MNRAS, 474, 2116 [Google Scholar]
  11. Davies, B., & Dessart, L. 2019, MNRAS, 483, 887 [NASA ADS] [CrossRef] [Google Scholar]
  12. de Jager, C., Nieuwenhuijzen, H., & van der Hucht, K. A. 1988, A&AS, 72, 259 [NASA ADS] [Google Scholar]
  13. Dessart, L., & Hillier, D. J. 2011a, MNRAS, 410, 1739 [NASA ADS] [Google Scholar]
  14. Dessart, L., & Hillier, D. J. 2011b, MNRAS, 415, 3497 [NASA ADS] [CrossRef] [Google Scholar]
  15. Dessart, L., Livne, E., & Waldman, R. 2010, MNRAS, 408, 827 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dessart, L., Hillier, D. J., Waldman, R., & Livne, E. 2013, MNRAS, 433, 1745 [NASA ADS] [CrossRef] [Google Scholar]
  17. Dessart, L., John Hillier, D., & Audit, E. 2017, A&A, 605, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Dessart, L., Hillier, D. J., & Wilk, K. D. 2018, A&A, 619, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Falk, S. W., & Arnett, W. D. 1977, ApJS, 33, 515 [NASA ADS] [CrossRef] [Google Scholar]
  20. Förster, F., Moriya, T. J., & Maureira, J. C. 2018, Nat. Astron., 2, 574 [NASA ADS] [CrossRef] [Google Scholar]
  21. Fraser, M., Maund, J. R., Smartt, S. J., et al. 2012, ApJ, 759, L13 [NASA ADS] [CrossRef] [Google Scholar]
  22. Fraser, M., Maund, J. R., Smartt, S. J., et al. 2014, MNRAS, 439, L56 [NASA ADS] [CrossRef] [Google Scholar]
  23. Glas, R., Just, O., Janka, H. T., & Obergaulinger, M. 2019, ApJ, 873, 45 [NASA ADS] [CrossRef] [Google Scholar]
  24. Grassberg, E. K., Imshennik, V. S., & Nadyozhin, D. K. 1971, Astrophys. Space Sci., 10, 28 [Google Scholar]
  25. Hamuy, M. 2003, ApJ, 582, 905 [NASA ADS] [CrossRef] [Google Scholar]
  26. Heger, A., Langer, N., & Woosley, S. E. 2000, ApJ, 528, 368 [NASA ADS] [CrossRef] [Google Scholar]
  27. Hillier, D. J., & Dessart, L. 2012, MNRAS, 424, 252 [NASA ADS] [CrossRef] [Google Scholar]
  28. Hirschi, R., Meynet, G., & Maeder, A. 2004, A&A, 425, 649 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Jerkstrand, A., Fransson, C., Maguire, K., et al. 2012, A&A, 546, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Jerkstrand, A., Smartt, S. J., Fraser, M., et al. 2014, MNRAS, 439, 3694 [NASA ADS] [CrossRef] [Google Scholar]
  31. Kasen, D., & Woosley, S. E. 2009, ApJ, 703, 2205 [NASA ADS] [CrossRef] [Google Scholar]
  32. Langer, N., Hamann, W. R., & Lennon, M. 1994, A&A, 290, 819 [NASA ADS] [Google Scholar]
  33. Lentz, E. J., Bruenn, S. W., Hix, W. R., et al. 2015, ApJ, 807, L31 [NASA ADS] [CrossRef] [Google Scholar]
  34. Leonard, D. C., Filippenko, A. V., Ganeshalingam, M., et al. 2006, Nature, 440, 505 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  35. Leonard, D. C., Pignata, G., & Dessart, L. 2012, ATel, 4033 [Google Scholar]
  36. Lisakov, S. M., Dessart, L., Hillier, D. J., Waldman, R., & Livne, E. 2017, MNRAS, 466, 34 [NASA ADS] [CrossRef] [Google Scholar]
  37. Lisakov, S. M., Dessart, L., Hillier, D. J., Waldman, R., & Livne, E. 2018, MNRAS, 473, 3863 [NASA ADS] [CrossRef] [Google Scholar]
  38. Litvinova, I. Y., & Nadezhin, D. K. 1985, Sov. Astron. Lett., 11, 145 [NASA ADS] [Google Scholar]
  39. Livne, E. 1993, ApJ, 412, 634 [NASA ADS] [CrossRef] [Google Scholar]
  40. Maeder, A., & Meynet, G. 1987, A&A, 182, 243 [NASA ADS] [Google Scholar]
  41. Maguire, K., Jerkstrand, A., Smartt, S. J., et al. 2012, MNRAS, 420, 3451 [NASA ADS] [CrossRef] [Google Scholar]
  42. Maund, J. R., Smartt, S. J., & Danziger, I. J. 2005, MNRAS, 364, L33 [NASA ADS] [CrossRef] [Google Scholar]
  43. Maund, J. R., Fraser, M., Smartt, S. J., et al. 2013, MNRAS, 431, L102 [NASA ADS] [CrossRef] [Google Scholar]
  44. Maund, J. R., Mattila, S., Ramirez-Ruiz, E., & Eldridge, J. J. 2014, MNRAS, 438, 1577 [NASA ADS] [CrossRef] [Google Scholar]
  45. Meynet, G., & Maeder, A. 2000, A&A, 361, 101 [NASA ADS] [Google Scholar]
  46. Meynet, G., Chomienne, V., Ekström, S., et al. 2015, A&A, 575, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Moriya, T. J., Yoon, S.-C., Gräfener, G., & Blinnikov, S. I. 2017, MNRAS, 469, L108 [NASA ADS] [CrossRef] [Google Scholar]
  48. Morozova, V., Piro, A. L., & Valenti, S. 2017, ApJ, 838, 28 [NASA ADS] [CrossRef] [Google Scholar]
  49. Morozova, V., Piro, A. L., & Valenti, S. 2018, ApJ, 858, 15 [NASA ADS] [CrossRef] [Google Scholar]
  50. Müller, B., Melson, T., Heger, A., & Janka, H.-T. 2017, MNRAS, 472, 491 [NASA ADS] [CrossRef] [Google Scholar]
  51. O’Connor, E. P., & Couch, S. M. 2018, ApJ, 865, 81 [NASA ADS] [CrossRef] [Google Scholar]
  52. O’Neill, D., Kotak, R., Fraser, M., et al. 2019, A&A, 622, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
  54. Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4 [NASA ADS] [CrossRef] [Google Scholar]
  55. Paxton, B., Marchant, P., Schwab, J., et al. 2015, ApJS, 220, 15 [Google Scholar]
  56. Paxton, B., Schwab, J., Bauer, E. B., et al. 2018, ApJS, 234, 34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Popov, D. V. 1993, ApJ, 414, 712 [Google Scholar]
  58. Sanders, N. E., Soderberg, A. M., Gezari, S., et al. 2015, ApJ, 799, 208 [NASA ADS] [CrossRef] [Google Scholar]
  59. Smartt, S. J. 2009, ARA&A, 47, 63 [NASA ADS] [CrossRef] [Google Scholar]
  60. Smartt, S. J., Maund, J. R., Hendry, M. A., et al. 2004, Science, 303, 499 [NASA ADS] [CrossRef] [Google Scholar]
  61. Sukhbold, T., Ertl, T., Woosley, S. E., Brown, J. M., & Janka, H.-T. 2016, ApJ, 821, 38 [Google Scholar]
  62. Turatto, M., Mazzali, P. A., Young, T. R., et al. 1998, ApJ, 498, L129 [NASA ADS] [CrossRef] [Google Scholar]
  63. Utrobin, V. P. 2007, A&A, 461, 233 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Van Dyk, S. D., Li, W., & Filippenko, A. V. 2003, PASP, 115, 1289 [NASA ADS] [CrossRef] [Google Scholar]
  65. Van Dyk, S. D., Davidge, T. J., Elias-Rosa, N., et al. 2012a, AJ, 143, 19 [NASA ADS] [CrossRef] [Google Scholar]
  66. Van Dyk, S. D., Cenko, S. B., Poznanski, D., et al. 2012b, ApJ, 756, 131 [NASA ADS] [CrossRef] [Google Scholar]
  67. Vartanyan, D., Burrows, A., Radice, D., Skinner, M. A., & Dolence, J. 2019, MNRAS, 482, 351 [NASA ADS] [CrossRef] [Google Scholar]
  68. Woosley, S. E. 1988, ApJ, 330, 218 [NASA ADS] [CrossRef] [Google Scholar]
  69. Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Rev. Mod. Phys., 74, 1015 [NASA ADS] [CrossRef] [Google Scholar]
  70. Yaron, O., Perley, D. A., Gal-Yam, A., et al. 2017, Nat. Phys., 13, 510 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Properties of the model set.

All Figures

thumbnail Fig. 1.

Density profiles vs. Lagrangian mass at the time of core collapse for a set of stellar models with an initial mass on the main sequence between 11 and 30 M (Woosley et al. 2002). The origin is at star center (left column) or at the He core edge (right column). The H-rich envelope is drawn with a thick black line and corresponds in all models to (extended) regions with a density below about 10−7 g cm−3. The He core is drawn with a thick red line and corresponds in all models to (compact) regions with a density above about 1 g cm−3. The vertical scale is squeezed but is suitable to reveal the evolution of the He core mass and the H-rich envelope mass for different ZAMS masses. Left panel: smooth progression in the He core mass with initial progenitor mass, and the varying progenitor mass at the onset of core collapse. Right panel: nonlinear evolution of the H-rich envelope mass with initial progenitor mass is emphasized. In these models, stars with ZAMS masses of 11–21 M contain similar H-rich envelope masses (8–9 M) at the onset of core collapse.

In the text
thumbnail Fig. 2.

Density structure of the models m12, m15, m20, and m25 shown vs. Lagrangian mass (the origin is the outer edge of the He core). The differences in H-rich envelope density reflect the differences in surface radii (see Table 1).

In the text
thumbnail Fig. 3.

Ejecta velocity, fractional kinetic energy Ek, frac, density, and temperature vs. mass at 14 d after the piston trigger. Ek, frac is defined as ∫V2dm/2Ekin, with the integration done inward. The x-axis origin is at the outermost ejecta layer in order to emphasize how similar the properties of the shocked H-rich envelope are in these different progenitor models.

In the text
thumbnail Fig. 4.

Summary of results from the CMFGEN simulations of models m12, m15, m20, and m25, showing, from top to bottom, the bolometric light curve, the V-band light curve, the V − I color evolution, and the evolution of some photospheric properties (T3, ph: temperature in units of 1000 K; V3, ph: velocity in units of 1000 km s−1; XH, ph: hydrogen mass fraction; IH, ph: hydrogen ionization state).

In the text
thumbnail Fig. 5.

Comparison of optical spectra for models m12, m15, m20, and m25 at about 20 d (top), 62 d (middle), and 155 d (bottom) after explosion. The color offset at early times results from the different progenitor radii. The greater O I 6300 Å line in more massive progenitors arises from their greater oxygen content.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.