Pulsational Pair-instability Supernovae Research Articles

Properties and Astrophysical Implications of the 150 M⊙ Binary Black Hole Merger GW190521

Abstract The gravitational-wave signal GW190521 is consistent with a binary black hole (BBH) merger source at redshift 0.8 with unusually high component masses, 85 − 14 + 21 M ⊙ and 66 − 18 + 17 M ⊙, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range 65–120 M ⊙. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger ( 142 − 16 + 28 M ⊙) classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular BBH coalescence, we detail the physical properties of GW190521’s source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of the coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be 0.13 − 0.11 + 0.30 Gpc − 3 yr − 1 . We discuss the astrophysical implications of GW190521 for stellar collapse and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescences, or via hierarchical mergers of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.

Predictions for the hydrogen-free ejecta of pulsational pair-instability supernovae

Present and upcoming time-domain astronomy efforts, in part driven by gravitational-wave follow-up campaigns, will unveil a variety of rare explosive transients in the sky. Here, we focus on pulsational pair-instability evolution, which can result in signatures that are observable with electromagnetic and gravitational waves. We simulated grids of bare helium stars to characterize the resulting black hole (BH) masses together with the ejecta composition, velocity, and thermal state. We find that the stars do not react “elastically” to the thermonuclear ignition in the core: there is not a one-to-one correspondence between pair-instability driven ignition and mass ejections, which causes ambiguity as to what is an observable pulse. In agreement with previous studies, we find that for initial helium core masses of 37.5 M⊙ ≲ MHe, init ≲ 41 M⊙, corresponding to carbon-oxygen core masses 27.5 M⊙ ≲ MCO ≲ 30.1 M⊙, the explosions are not strong enough to affect the surface. With increasing initial helium core mass, they become progressively stronger causing first large radial expansion (41 M⊙ ≲ MHe, init ≲ 42 M⊙, corresponding to 30.1 M⊙ ≲ MCO ≲ 30.8 M⊙) and, finally, also mass ejection episodes (for MHe, init ≳ 42 M⊙, or MCO ≳ 30.8 M⊙). The lowest mass helium core to be fully disrupted in a pair-instability supernova is MHe, init ≃ 80 M⊙, corresponding to MCO ≃ 55 M⊙. Models with MHe, init ≳ 200 M⊙ (MCO ≳ 114 M⊙) reach the photodisintegration regime, resulting in BHs with masses of MBH ≳ 125 M⊙. Although this is currently considered unlikely, if BHs from these models form via (weak) explosions, the previously-ejected material might be hit by the blast wave and convert kinetic energy into observable electromagnetic radiation. We characterize the hydrogen-free circumstellar material from the pulsational pair-instability of helium cores by simply assuming that the ejecta maintain a constant velocity after ejection. We find that our models produce helium-rich ejecta with mass of 10−3 M⊙ ≲ MCSM ≲ 40 M⊙, the larger values corresponding to the more massive progenitor stars. These ejecta are typically launched at a few thousand km s−1 and reach distances of ∼1012 − 1015 cm before the core-collapse of the star. The delays between mass ejection events and the final collapse span a wide and mass-dependent range (from subhour to 104 years), and the shells ejected can also collide with each other, powering supernova impostor events before the final core-collapse. The range of properties we find suggests a possible connection with (some) type Ibn supernovae.

An extremely energetic supernova from a very massive star in a dense medium

The interaction of a supernova with a circumstellar medium (CSM) can dramatically increase the emitted luminosity by converting kinetic energy to thermal energy. In 'superluminous' supernovae (SLSNe) of Type IIn -- named for narrow hydrogen lines in their spectra -- the integrated emission can reach $\sim 10^{51}$ erg, attainable by thermalising most of the kinetic energy of a conventional SN. A few transients in the centres of active galaxies have shown similar spectra and even larger energies, but are difficult to distinguish from accretion onto the supermassive black hole. Here we present a new event, SN2016aps, offset from the centre of a low-mass galaxy, that radiated $\gtrsim 5 \times 10^{51}$ erg, necessitating a hyper-energetic supernova explosion. We find a total (SN ejecta $+$ CSM) mass likely exceeding 50-100 M$_\odot$, with energy $\gtrsim 10^{52}$ erg, consistent with some models of pair-instability supernovae (PISNe) or pulsational PISNe -- theoretically-predicted thermonuclear explosions from helium cores $>50$ M$_\odot$. Independent of the explosion mechanism, this event demonstrates the existence of extremely energetic stellar explosions, detectable at very high redshifts, and provides insight into dense CSM formation in the most massive stars.

Pulsational Pair-instability Supernovae. II. Neutrino Signals from Pulsations and Their Detection by Terrestrial Neutrino Detectors

Abstract A Pulsational Pair-instability supernova (PPISN) evolves from a massive star with a mass ∼80–140 M ⊙ that develops electron–positron pair-instability after hydrostatic He-burning in the core has finished. In Leung et al. (Paper I), we examined the evolutionary tracks and the pulsational mass-loss history of this class of stars. In this paper, we analyze the thermodynamical history to explore the neutrino observables of PPISNe. We compute the neutrino light curves and spectra during pulsation. We then study the detailed neutrino emission profiles of these stars and estimate the expected neutrino detection count for different terrestrial neutrino detectors, including, e.g., KamLAND and Super-Kamiokande. Finally, we compare the neutrino pattern of PPISN with other types of supernovae based on a canonical 10 kt detector. The predicted neutrino signals can provide an early warning for telescopes to trace for the early time optical signals. The implications of neutrino physics on the expected detection are also discussed.

Pulsational Pair-instability Supernovae. I. Pre-collapse Evolution and Pulsational Mass Ejection

Abstract We calculate the evolution of massive stars, which undergo pulsational pair-instability (PPI) when the O-rich core is formed. The evolution from the main sequence through the onset of PPI is calculated for stars with initial masses of 80–140 M ⊙ and metallicities of Z = 10−3−1.0 Z ⊙. Because of mass loss, Z ≤ 0.5 Z ⊙ is necessary for stars to form He cores massive enough (i.e., mass >40 M ⊙) to undergo PPI. The hydrodynamical phase of evolution from PPI through the beginning of Fe-core collapse is calculated for He cores with masses of 40−62 M ⊙ and Z = 0. During PPI, electron–positron pair production causes a rapid contraction of the O-rich core, which triggers explosive O-burning and a pulsation of the core. We study the mass dependence of the pulsation dynamics, thermodynamics, and nucleosynthesis. The pulsations are stronger for more massive He cores and result in a large amount of mass ejection such as 3–13 M ⊙ for 40−62 M ⊙ He cores. These He cores eventually undergo Fe-core collapse. The 64 M ⊙ He core undergoes complete disruption and becomes a pair-instability supernova. The H-free circumstellar matter ejected around these He cores is massive enough to explain the observed light curve of Type I (H-free) superluminous supernovae with circumstellar interaction. We also note that the mass ejection sets the maximum mass of black holes (BHs) to be ∼50 M ⊙, which is consistent with the masses of BHs recently detected by VIRGO and aLIGO.

Predictions for the hydrogen-free ejecta of pulsational pair-instability supernovae

Predictions for the hydrogen-free ejecta of pulsational pair-instability supernovae

The Impact of Pair-instability Mass Loss on the Binary Black Hole Mass Distribution

A population of binary black hole mergers has now been observed in gravitational waves by Advanced LIGO and Virgo. The masses of these black holes appear to show evidence for a pile-up between $30$--$45$ $M_\odot$ and a cut-off above $\sim 45$ $M_\odot$. One possible explanation for such a pile-up and subsequent cut-off are pulsational pair-instability supernovae (PPISNe) and pair-instability supernovae (PISNe) in massive stars. We investigate the plausibility of this explanation in the context of isolated massive binaries. We study a population of massive binaries using the rapid population synthesis software COMPAS, incorporating models for PPISNe and PISNe. Our models predict a maximum black hole mass of $40$ $M_\odot$. We expect $\sim 10$\% of all binary black hole mergers at redshift z = 0 will include at least one component that went through a PPISN (with mass $30$--$40$ $M_\odot$), constituting $\sim 20$--$50$\% of binary black hole mergers observed during the first two observing runs of Advanced LIGO and Virgo. Empirical models based on fitting the gravitational-wave mass measurements to a combination of a power law and a Gaussian find a fraction too large to be associated with PPISNe in our models. The rates of PPISNe and PISNe track the low metallicity star formation rate, increasing out to redshift $z = 2$. These predictions may be tested both with future gravitational-wave observations and with observations of superluminous supernovae.

Pulsational Pair-instability Supernovae in Very Close Binaries

Abstract Pair-instability and pulsational pair-instability supernovae (PPISNe) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. If these BBHs are the product of binary evolution, then PPISNe could occur in very close binaries. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion and what impact these events have on the formation of merging BBHs through binary evolution. For this, we have computed a set of models of metal-poor (Z ⊙/10) single helium stars using the MESA software instrument. For PPISN progenitors with pre-pulse masses >50 M ⊙ we find that, after a pulse, heat deposited throughout the layers of the star that remain bound causes it to expand to more than 100 R ⊙ for periods of 102–104 yr depending on the mass of the progenitor. This results in long-lived phases of Roche lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated with PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final BH spin by ∼30%, induces eccentricities below the threshold of detectability of the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors.

SN 2016iet: The Pulsational or Pair Instability Explosion of a Low-metallicity Massive CO Core Embedded in a Dense Hydrogen-poor Circumstellar Medium

Abstract We present optical photometry and spectroscopy of SN 2016iet (=Gaia16bvd=PS17brq), an unprecedented Type I supernova (SN I) at z = 0.0676 with no obvious analog in the existing literature. SN 2016iet exhibits a peculiar light curve, with two roughly equal brightness peaks (≈ −19 mag) separated by about 100 days, and a subsequent slow decline by about 5 mag in 650 rest-frame days. The spectra are dominated by strong emission lines of calcium and oxygen, with a width of only 3400 km s−1, superposed on a strong blue continuum in the first year. There is no clear evidence for hydrogen or helium associated with the SN at any phase. The nebular spectra exhibit a ratio of , much larger than for core-collapse SNe and Type I superluminous SNe. We model the light curves with several potential energy sources: radioactive decay, a central engine, and ejecta–circumstellar medium (CSM) interaction. Regardless of the model, the inferred progenitor mass near the end of its life (i.e., the CO core mass) is ≳55 M ⊙ and potentially up to 120 M ⊙, clearly placing the event in the regime of pulsational pair instability supernovae (PPISNe) or pair instability supernovae (PISNe). The models of CSM interaction provide the most consistent explanation for the light curves and spectra, and require a CSM mass of ≈35 M ⊙ ejected in the final decade before explosion. We further find that SN 2016iet is located at an unusually large projected offset (16.5 kpc, 4.3 effective radii) from its low-metallicity dwarf host galaxy (Z ≈ 0.1 Z ⊙, L ≈ 0.02 L *, M ≈ 108.5 M ⊙), supporting the interpretation of a PPISN/PISN explosion. In our final spectrum at a phase of about 770 rest-frame days we detect weak and narrow Hα emission at the location of the SN, corresponding to a star formation rate of ≈3 × 10−4 M ⊙ yr−1, which is likely due to a dim underlying galaxy host or an H ii region. Despite the overall consistency of the SN and its unusual environment with PPISNe and PISNe, we find that the inferred properties of SN 2016iet challenge existing models of such events.

Fallback Accretion-powered Supernova Light Curves Based on a Neutrino-driven Explosion Simulation of a 40 M⊙ Star

Abstract We present synthetic light curves (LCs) of fallback-powered supernovae based on a neutrino-driven explosion of a 40 zero-metallicity star with significant fallback accretion onto a black hole that was previously simulated by Chan et al. until shock breakout. Here, we investigate the LC properties of the explosion after shock breakout for various fallback models. Without extra power from fallback accretion, the LC is that of an SN IIP with a plateau magnitude of around −14 mag and a plateau duration of 40 days. With extra power for the LC from fallback accretion, however, we find that the transient can be significantly more luminous. The LC shape can be SN 1987A-like or Type IIP-like, depending on the efficiency of the fallback engine. If the accretion disk forms soon after the collapse and more than 1% of the accretion energy is released as the central engine, fallback accretion-powered supernovae become as luminous as superluminous supernovae. We suggest that Type II superluminous supernovae with broad hydrogen features could be related to such hydrogen-rich supernovae powered by fallback accretion. In the future, such superluminous supernovae powered by fallback accretion might be found among the supernovae from the first stars in addition to pair-instability supernovae and pulsational pair-instability supernovae.

The Evolution of Massive Helium Stars, Including Mass Loss

Abstract The evolution of helium stars with initial masses in the range 1.6–120 is studied, including the effects of mass loss by winds. These stars are assumed to form in binary systems when their expanding hydrogenic envelopes are promptly lost just after helium ignition. Significant differences are found with single-star evolution, chiefly because the helium core loses mass during helium burning rather than gaining it from hydrogen shell burning. Consequently, presupernova stars for a given initial mass function have considerably smaller mass when they die and will be easier to explode. Even accounting for this difference, the helium stars with mass loss develop more centrally condensed cores that should explode more easily than their single-star counterparts. The production of low-mass black holes may be diminished. Helium stars with initial masses below 3.2 experience significant radius expansion after helium depletion, reaching blue supergiant proportions. This could trigger additional mass exchange or affect the light curve of the supernova. The most common black hole mass produced in binaries is estimated to be about 9 . A new maximum mass for black holes derived from pulsational pair-instability supernovae is derived, 46 , and a new potential gap at 10–12 is noted. Models pertinent to SN 2014ft are presented, and a library of presupernova models is generated.

Massive Stellar Mergers as Precursors of Hydrogen-rich Pulsational Pair Instability Supernovae

Abstract Interactions between massive stars in binaries are thought to be responsible for much of the observed diversity of supernovae. As surveys probe rarer populations of events, we should expect to see supernovae arising from increasingly uncommon progenitor channels. Here we examine a scenario in which massive stars merge after they have both formed a hydrogen-exhausted core. We suggest that this could produce stars that explode as pair-instability supernovae (PISNe) with significantly more hydrogen, at a given metallicity, than in single-star models with the same pre-explosion oxygen-rich core mass. We investigate the subset of those stellar mergers that later produce pulsational PISNe, and estimate that the rate of such post-merger, hydrogen-rich pulsational PISNe could approach a few in a thousand of all core-collapse supernovae. The nature and predicted rate of such hydrogen-rich pulsational PISNe are reminiscent of the very unusual supernova iPTF14hls. For plausible assumptions, PISNe from similar mergers might dominate the rate of PISNe in the local Universe.

Spin orientations of merging black holes formed from the evolution of stellar binaries

We study the expected spin misalignments of merging binary black holes (BHs) formed in isolation by combining state-of-the-art population-synthesis models with efficient post-Newtonian evolutions, thus tracking sources from stellar formation to gravitational-wave detection. We present extensive predictions of the properties of sources detectable by both current and future interferometers. We account for the fact that detectors are more sensitive to spinning BH binaries with suitable spin orientations and find that this significantly impacts the population of sources detectable by LIGO, while this is not the case for 3rd-generation detectors. We find that three formation pathways, differentiated by the order of core collapse and common-envelope phases, dominate the observed population, and that their relative importance critically depends on the recoils imparted to BHs at birth. Our models suggest that measurements of the "effective-spin" parameter $\chi_{\rm eff}$ will allow for powerful constraints. For instance, we find that the role of spin magnitudes and spin directions in $\chi_{\rm eff}$ can be largely disentangled, and that the symmetry of the effective-spin distribution is a robust indicator of the binary's formation history. Our predictions for individual spin directions and their precessional morphologies confirm and extend early toy models, while exploring substantially more realistic and broader sets of initial conditions. Our main conclusion is that specific subpopulations of BH binaries will exhibit distinctive precessional dynamics: these classes include (but are not limited to) sources where stellar tidal interactions act on sufficiently short timescales, and massive binaries produced in pulsational pair-instability supernovae. Measurements of BH spin orientations have enormous potential to constrain specific evolutionary processes in the lives of massive binary stars.

Models for the Unusual Supernova iPTF14hls

Abstract Supernova iPTF14hls maintained a bright, variable luminosity for more than 600 days, while lines of hydrogen and iron in its spectrum had different speeds but showed little evolution. Here, several varieties of models are explored for iPTF14hls-like events. They are based upon circumstellar medium (CSM) interaction in an ordinary supernova, pulsational pair-instability supernovae (PPISN), and magnetar formation. Each is able to explain the enduring emission and brightness of iPTF14hls but each has shortcomings when confronted with other observed characteristics. The PPISN model can, in some cases, produce a presupernova transient like the one observed at the site of iPTF14hls in 1954. It also offers a clear path to providing the necessary half solar mass of material at ∼5 × 1016 cm for CSM interaction to work and it can give an irregular light curve without invoking additional assumptions. It explains the 4000 km s−1 seen in the iron lines but without additional energy input it strains to explain the nearly constant 8000 km s−1 velocity seen in H α . Magnetar models can also explain most of the observed features but they give a smooth light curve and may be difficult to reconcile with the observation of slow-moving hydrogen at late times. The various models predict different spectral characteristics and a remnant that, today, could be a black hole, magnetar, or even a star. Further observations and calculations of radiation transport will narrow the range of possibilities.

Measuring the Binary Black Hole Mass Spectrum with an Astrophysically Motivated Parameterization

Abstract Gravitational-wave detections have revealed a previously unknown population of stellar mass black holes with masses above 20 M ⊙. These observations provide a new way to test models of stellar evolution for massive stars. By considering the astrophysical processes likely to determine the shape of the binary black hole mass spectrum, we construct a parameterized model to capture key spectral features that relate gravitational-wave data to theoretical stellar astrophysics. In particular, we model the signature of pulsational pair-instability supernovae, which are expected to cause all stars with initial mass 100 M ⊙ ≲ M ≲ 150 M ⊙ to form ∼40 M ⊙ black holes. This would cause a cutoff in the black hole mass spectrum along with an excess of black holes near 40 M ⊙. We carry out a simulated data study to illustrate some of the stellar physics that can be inferred using gravitational-wave measurements of binary black holes and demonstrate several such inferences that might be made in the near future. First, we measure the minimum and maximum stellar black hole mass. Second, we infer the presence of a peak due to pair-instability supernovae. Third, we measure the distribution of black hole mass ratios. Finally, we show how inadequate models of the black hole mass spectrum lead to biased estimates of the merger rate and the amplitude of the stochastic gravitational-wave background.

Modules for Experiments in Stellar Astrophysics (${\mathtt{M}}{\mathtt{E}}{\mathtt{S}}{\mathtt{A}}$): Convective Boundaries, Element Diffusion, and Massive Star Explosions

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact.

Very massive stars, pair-instability supernovae and intermediate-mass black holes with the sevn code

Understanding the link between massive ($\gtrsim 30$ M$_{\odot{}}$) stellar black holes (BHs) and their progenitor stars is a crucial step to interpret observations of gravitational-wave events. In this paper, we discuss the final fate of very massive stars (VMSs), with zero-age main sequence (ZAMS) mass $>150$ M$_{\odot}$, accounting for pulsational pair-instability supernovae (PPISNe) and for pair-instability supernovae (PISNe). We describe an updated version of our population synthesis code SEVN, in which we added stellar evolution tracks for VMSs with ZAMS mass up to $350 $M$_{\odot{}}$ and we included analytical prescriptions for PPISNe and PISNe. We use the new version of SEVN to study the BH mass spectrum at different metallicity $Z$, ranging from $Z=2.0\times 10^{-4}$ to $Z=2.0\times 10^{-2}$. The main effect of PPISNe and PISNe is to favour the formation of BHs in the mass range of the first gravitational-wave event (GW150914), while they prevent the formation of remnants with mass 60 - 120 M$_{\odot{}}$. In particular, we find that PPISNe significantly enhance mass loss of metal-poor ($Z\leq 2.0\times 10^{-3}$) stars with ZAMS mass $60\leq $M$_{\mathrm{ZAMS}}/$M$_{\odot{}}\leq 125$. In contrast, PISNe become effective only for moderately metal-poor ($Z<8.0\times 10^{-3}$) VMSs. VMSs with M$_{\rm ZAMS}\gtrsim{}220$ M$_\odot$ and $Z<10^{-3}$ do not undergo PISNe and form intermediate-mass BHs (IMBHs, with mass $\gtrsim 200 $M$_{\odot{}}$) via direct collapse.

Pulsational Pair-instability Supernovae

Abstract The final evolution of stars in the mass range 70–140 M ⊙ is explored. Depending upon their mass loss history and rotation rates, these stars will end their lives as pulsational pair-instability supernovae (PPISN) producing a great variety of observational transients with total durations ranging from weeks to millennia and luminosities from 1041 to over 1044 erg s−1. No nonrotating model radiates more than 5 × 10 50 erg of light or has a kinetic energy exceeding 5 × 10 51 erg, but greater energies are possible, in principle, in magnetar-powered explosions, which are explored. Many events resemble SNe Ibn, SNe Icn, and SNe IIn, and some potential observational counterparts are mentioned. Some PPISN can exist in a dormant state for extended periods, producing explosions millennia after their first violent pulse. These dormant supernovae contain bright Wolf–Rayet stars, possibly embedded in bright X-ray and radio sources. The relevance of PPISN to supernova impostors like Eta Carinae, to superluminous supernovae, and to sources of gravitational radiation is discussed. No black holes between 52 and 133 M ⊙ are expected from stellar evolution in close binaries.

Pulsational Pair-instability Model for Superluminous Supernova PTF12dam:Interaction and Radioactive Decay

Abstract Being a superluminous supernova, PTF12dam can be explained by a 56Ni-powered model, a magnetar-powered model, or an interaction model. We propose that PTF12dam is a pulsational pair-instability supernova, where the outer envelope of a progenitor is ejected during the pulsations. Thus, it is powered by a double energy source: radioactive decay of 56Ni and a radiative shock in a dense circumstellar medium. To describe multicolor light curves and spectra, we use radiation-hydrodynamics calculations of the STELLA code. We found that light curves are well described in the model with 40 M ⊙ ejecta and 20–40 M ⊙ circumstellar medium. The ejected 56Ni mass is about 6 M ⊙, which results from explosive nucleosynthesis with large explosion energy (2–3) × 1052 erg. In comparison with alternative scenarios of pair-instability supernova and magnetar-powered supernova, in the interaction model, all the observed main photometric characteristics are well reproduced: multicolor light curves, color temperatures, and photospheric velocities.

EXTREME SUPERNOVA MODELS FOR THE SUPER-LUMINOUS TRANSIENT ASASSN-15LH

ABSTRACT The recent discovery of the unprecedentedly super-luminous transient ASASSN-15lh (or SN 2015L) with its UV-bright secondary peak challenges all the power-input models that have been proposed for super-luminous supernovae. Here we examine some of the few viable interpretations of ASASSN-15lh in the context of a stellar explosion, involving combinations of one or more power inputs. We model the light curve of ASASSN-15lh with a hybrid model that includes contributions from magnetar spin-down energy and hydrogen-poor circumstellar interaction. We also investigate models of pure circumstellar interaction with a massive hydrogen-deficient shell and discuss the lack of interaction features in the observed spectra. We find that, as a supernova, ASASSN-15lh can be best modeled by the energetic core-collapse of an ∼40 M ⊙ star interacting with a hydrogen-poor shell of ∼20 M ⊙. The circumstellar shell and progenitor mass are consistent with a rapidly rotating pulsational pair-instability supernova progenitor as required for strong interaction following the final supernova explosion. Additional energy injection by a magnetar with an initial period of 1–2 ms and magnetic field of 0.1–1 × 1014 G may supply the excess luminosity required to overcome the deficit in single-component models, but this requires more fine-tuning and extreme parameters for the magnetar, as well as the assumption of efficient conversion of magnetar energy into radiation. We thus favor a single-input model where the reverse shock formed in a strong SN ejecta–circumstellar matter interaction following a very powerful core-collapse SN explosion can supply the luminosity needed to reproduce the late-time UV-bright plateau.