The difficulty of inferring progenitor masses from Type II-Plateau supernova light curves
Luc Dessart, D. John Hillier

TL;DR
This study demonstrates that Type II-P supernova light curves are degenerate with respect to progenitor mass, making it difficult to infer initial stellar masses solely from light curve data due to similar observable properties across different progenitor masses.
Contribution
The paper provides a controlled experiment showing that supernova light curves and spectra are degenerate across different progenitor masses with similar envelope masses, challenging the reliability of mass inference from light curves.
Findings
Light curves are similar for different progenitor masses with comparable envelope masses.
Spectral features are degenerate and influenced by progenitor radius and oxygen yield.
Light curve modeling cannot reliably determine progenitor mass due to degeneracies.
Abstract
Much controversy surrounds the inferred progenitor masses of Type II-Plateau (II-P) supernovae (SNe). The debate is nourished by the discrepant results from radiation-hydrodynamics simulations, from pre-explosion imaging, and from studies of host stellar populations. Here, we present a controlled experiment using four solar metallicity models with zero-age main-sequence masses of 12, 15, 20, and 25Msun. Because of the effects of core burning and surface mass loss, these models reach core collapse as red-supergiant (RSG) stars with a similar H-rich envelope mass of 8 to 9Msun but with final masses in the range 11 to 16Msun. We explode the progenitors using a thermal bomb, adjusting the energy deposition to yield an asymptotic ejecta kinetic energy of 1.25 x 10^51 erg and an initial 56Ni mass of 0.04Msun. The resulting SNe produce similar photometric and spectroscopic properties from 10…
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8| (i) | ||||||||
|---|---|---|---|---|---|---|---|---|
| [M⊙] | [M⊙] | [M⊙] | [R⊙] | [erg] | [M⊙] | [M⊙] | [M⊙] | |
| m12 | 11.21 | 3.26 | 7.95 | 406 | 1.27(51) | 9.79 | 0.16 | 0.046 |
| m15 | 13.12 | 4.56 | 8.56 | 589 | 1.28(51) | 11.57 | 0.53 | 0.041 |
| m20 | 14.86 | 6.81 | 8.05 | 843 | 1.25(51) | 13.18 | 1.33 | 0.041 |
| m25 | 15.96 | 8.59 | 7.37 | 872 | 1.21(51) | 13.93 | 2.33 | 0.042 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
11institutetext: Unidad Mixta Internacional Franco-Chilena de Astronomía (CNRS, UMI 3386), Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile 22institutetext: Department of Physics and Astronomy & Pittsburgh Particle Physics, Astrophysics, and Cosmology Center (PITT PACC), University of Pittsburgh, 3941 O’Hara Street, Pittsburgh, PA 15260, USA.
The difficulty of inferring progenitor masses from Type II-Plateau supernova light curves
Luc Dessart The difficulty of inferring progenitor masses from Type II-Plateau supernova light curvesThe difficulty of inferring progenitor masses from Type II-Plateau supernova light curves
D. John Hillier The difficulty of inferring progenitor masses from Type II-Plateau supernova light curvesThe difficulty of inferring progenitor masses from Type II-Plateau supernova light curves
(Received; accepted)
Much controversy surrounds the inferred progenitor masses of Type II-Plateau (II-P) supernovae (SNe). The debate is nourished by the discrepant results from radiation-hydrodynamics simulations, from pre-explosion imaging, and from studies of host stellar populations. Here, we present a controlled experiment using four solar metallicity models with zero-age main-sequence masses of 12, 15, 20, and 25 M*⊙. Because of the effects of core burning and surface mass loss, these models reach core collapse as red-supergiant (RSG) stars with a similar H-rich envelope mass of 8 to 9 M⊙* but with final masses in the range to M*⊙. We explode the progenitors using a thermal bomb, adjusting the energy deposition to yield an asymptotic ejecta kinetic energy of erg and an initial i mass of 0.04 M⊙*. The resulting SNe produce similar photometric and spectroscopic properties from 10 to 200 d. The spectral characteristics are degenerate. The scatter in early-time color results from the range in progenitor radii, while the differences in late-time spectra reflect the larger oxygen yields in more massive progenitors. Because the progenitors have a comparable H-rich envelope mass, the photospheric phase duration is comparable for all models; the difference in He-core mass is invisible. As different main-sequence masses can produce progenitors with a similar H-rich envelope mass, light curve modeling cannot provide a robust and unique solution for the ejecta mass of Type II-P SNe. The numerous uncertainties in massive star evolution and wind mass loss also prevent a robust association with a main-sequence star mass. Light curve modeling can at best propose compatibility.
Key Words.:
** radiative transfer – radiation: dynamics – supernovae: general **
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 evolution111We will 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., 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 SN II-P observations (Lentz et al., 2015; Müller et al., 2017; Glas et al., 2018; O’Connor & Couch, 2018; Vartanyan et al., 2019), the simulations cannot predict which stars produce core collapse SNe, and quantities such as the i 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 having 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 SNe II-P (Van Dyk et al., 2003, 2012b, 2012a; 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 either not known a priori (e.g., initial rotation) or 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 SNe II-P. The most robust inference and the most easily done is for the i 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 i and o decay, and the mass of ejected i 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, 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 done for several decades. Early works identified the basic properties of the progenitors, in particular their large progenitor radii and massive H-rich envelope compatible with RSG progenitors (Grassberg et al., 1971; Falk & Arnett, 1977). Basic relations between ejecta and observed SN II-P properties were also drawn (Litvinova & Nadezhin, 1985; Popov, 1993). When applied to observed SNe II-P, 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.
SN II-P 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). Hence, 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 inter-changeably 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 utilize the public grid of massive star models from Woosley et al. (2002). Fig. 1 shows a montage of density profiles at core collapse for massive stars evolved at solar metallicity with zero-age main sequence masses () 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 increases, the mass of the He core (; 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 2018). The mass of the H-rich envelope () at core collapse is between 8 and 9 M*⊙* for stars with between 11 and 21 M*⊙*, and decreases (not necessarily monotonically) as 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 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 where starts going down with increasing , and the actual value of for a given .222These quantities should also vary with metallicity. The trend shown in Fig. 1 must hold.
Because SN II-P 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 below 21 M*⊙* have similar density structures, while above the models have progressively lower values of . For a given ejecta kinetic energy imparted to the H-rich envelope, models with below 21 M*⊙* should have comparable light curves (modulo differences in ), 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., values) inferred from SN II-P light curves should consequently be around 10 M⊙* since the RSG stars from 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). And 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 SN II-P 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 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 non-rotating 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 unsuited for light curve modeling. To comply with the early-time observations of some SNe II-P (see, e.g., Yaron et al. 2017), about 0.2 M⊙* of atmospheric (static) material is added above using a density scale height of about 0.2 . This has however no relevance for the present discussion.
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 to 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 of erg. Although v1d treats explosive nucleosynthesis, the i 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 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 10. 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 SN II-P. It also stores the bulk of the kinetic energy.
At 14 d, the ejecta models computed by v1d are remapped into the non local thermodynamic equilibrium 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 i 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 erg and a i mass within a factor of two of 0.04 M*⊙* (their is some scatter and some offset with depending on the calibration used for their “thermal bomb”). Our choice is thus not unrealistic.
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 -band light curve and the 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.
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 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 i 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 yield the same asymptotic ejecta kinetic energy (which is what we assumed here for our model set) and therefore naturally contribute 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 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 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 non-evolutionary 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 Section 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 when comparing model predictions with observations. While not all gave consistency with the inferred photospheric velocities, it does 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- 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.
A more reliable discriminant for progenitor mass may thus 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 SN II-P with a similar plateau duration and spectral evolution if they have the same H-rich envelope mass and ejecta kinetic energy. Differences arise then 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. And stars with different main-sequence mass can die with the same (which is function of the uncertain cumulative mass loss), so that the constraint of from light curve modeling cannot be uniquely connected to a . 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 and i mass, and we may eventually have a robust prediction of these quantities from ab-initio three-dimensional 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 non-LTE, 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 so 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.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Anderson et al. (2014) Anderson, J. P., González-Gaitán, S., Hamuy, M., et al. 2014, Ap J, 786, 67
- 2Arnett et al. (2015) Arnett, W. D., Meakin, C., Viallet, M., et al. 2015, Ap J, 809, 30
- 3Bartunov & Blinnikov (1992) Bartunov, O. S. & Blinnikov, S. I. 1992, Soviet Astronomy Letters, 18, 43
- 4Beasor & Davies (2016) Beasor, E. R. & Davies, B. 2016, MNRAS, 463, 1269
- 5Bersten et al. (2011) Bersten, M. C., Benvenuto, O., & Hamuy, M. 2011, Ap J, 729, 61
- 6Burrows & Goshy (1993) Burrows, A. & Goshy, J. 1993, Ap JL, 416, L 75
- 7Chugai & Utrobin (2000) Chugai, N. N. & Utrobin, V. P. 2000, A&A, 354, 557
- 8Crowther (2007) Crowther, P. A. 2007, ARA&A, 45, 177
