Search

The cores of stars are the cosmic furnaces where light elements are fused into heavier nuclei. The fusion of hydrogen to helium initially powers all stars. The ashes of the fusion reactions are then predicted to serve as fuel in a series of stages, eventually transforming massive stars into a structure of concentric shells. These are composed of natal hydrogen on the outside, and consecutively heavier compositions inside, predicted to be dominated by helium, carbon/oxygen, oxygen/neon/magnesium, and oxygen/silicon/sulphur. Silicon and sulphur are fused into inert iron, leading to the collapse of the core and either a supernova explosion or the direct formation of a black hole. Stripped stars, where the outer hydrogen layer has been removed and the internal He-rich layer (in Wolf-Rayet WN stars) or even the C/O layer below it (in Wolf-Rayet WC/WO stars) are exposed, provide evidence for this shell structure, and the cosmic element production mechanism it reflects. The types of supernova explosions that arise from stripped stars embedded in shells of circumstellar material (most notably Type Ibn supernovae from stars with outer He layers, and Type Icn supernovae from stars with outer C/O layers) confirm this scenario. However, direct evidence for the most interior shells, which are responsible for the production of elements heavier than oxygen, is lacking. Here, we report the discovery of the first-of-its-kind supernova arising from a star peculiarly stripped all the way to the silicon and sulphur-rich internal layer. Whereas the concentric shell structure of massive stars is not under debate, it is the first time that such a thick, massive silicon and sulphur-rich shell, expelled by the progenitor shortly before the SN explosion, has been directly revealed., Comment: 48 pages, 12 figures and 10 tables. Submitted to a high-impact journal. The reduced spectra and photometry will be made available via the journal webpage and the WISeREP archive after the acceptance of the paper

Supernova (SN) 1987A offers a unique opportunity to study how a spatially resolved SN evolves into a young supernova remnant (SNR). We present and analyze Hubble Space Telescope (HST) imaging observations of SN 1987A obtained in 2022 and compare them with HST observations from 2009 to 2021. These observations allow us to follow the evolution of the equatorial ring (ER), the rapidly expanding ejecta, and emission from the center over a wide range in wavelength from 2000 to 11 000 AA. The ER has continued to fade since it reached its maximum ~8200 days after the explosion. In contrast, the ejecta brightened until day ~11000 before their emission levelled off; the west side brightened more than the east side, which we attribute to the stronger X-ray emission by the ER on that side. The asymmetric ejecta expand homologously in all filters, which are dominated by various emission lines from hydrogen, calcium, and iron. From this overall similarity, we infer the ejecta are chemically well-mixed on large scales. The exception is the diffuse morphology observed in the UV filters dominated by emission from the Mg II resonance lines that get scattered before escaping. The 2022 observations do not show any sign of the compact object that was inferred from highly-ionized emission near the remnant's center observed with JWST. We determine an upper limit on the flux from a compact central source in the [O III] HST image. The non-detection of this line indicates that the S and Ar lines observed with JWST originate from the O free inner Si - S - Ar rich zone and/or that the observed [O III] flux is strongly affected by dust scattering.

We present initial results from a JWST survey of the youngest Galactic core-collapse supernova remnant Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI/MRS IFU spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include: 1) a web-like network of unshocked ejecta filaments resolved to 0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor's oxygen layer with hot, high-entropy matter heated by neutrino interactions and radioactivity, 2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant's interior pockmarked with small (approximately one arcsecond) round holes formed by knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks, 3) dozens of light echoes with angular sizes between 0.1 arcsecond to 1 arcminute reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission from the neutron star in Cas A's center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects., Comment: 27 pages, 10 figures, now published in ApJL

JWST NIRCam observations at 1.5-4.5 $\mu$m have provided broad and narrow band imaging of the evolving remnant of SN 1987A with unparalleled sensitivity and spatial resolution. Comparing with previous marginally spatially resolved Spitzer IRAC observations from 2004-2019 confirms that the emission arises from the circumstellar equatorial ring (ER), and the current brightness at 3.6 and 4.5 $\mu$m was accurately predicted by extrapolation of the declining brightness tracked by IRAC. Despite the regular light curve, the NIRCam observations clearly reveal that much of this emission is from a newly developing outer portion of the ER. Spots in the outer ER tend to lie at position angles in between the well-known ER hotspots. We show that the bulk of the emission in the field can be represented by 5 standard spectral energy distributions (SEDs), each with a distinct origin and spatial distribution. This spectral decomposition provides a powerful technique for distinguishing overlapping emission from the circumstellar medium (CSM) and the supernova (SN) ejecta, excited by the forward and reverse shocks respectively., Comment: Accepted for publication in ApJ. 16 pages, 12 figures. 2 animations not included here

We use the light curve and spectral synthesis code JEKYLL to calculate a set of macroscopically mixed Type IIb supernova (SN) models, which are compared to both previously published and new late-phase observations of SN 2020acat. The models differ in the initial mass, the radial mixing and expansion of the radioactive material, and the properties of the hydrogen envelope. The best match to the photospheric and nebular spectra and lightcurves of SN 2020acat is found for a model with an initial mass of 17 solar masses, strong radial mixing and expansion of the radioactive material, and a 0.1 solar mass hydrogen envelope with a low hydrogen mass-fraction of 0.27. The most interesting result is that strong expansion of the clumps containing radioactive material seems to be required to fit the observations of SN 2020acat both in the diffusion phase and the nebular phase. These "Ni bubbles" are expected to expand due to heating from radioactive decays, but the degree of expansion is poorly constrained. Without strong expansion there is a tension between the diffusion phase and the subsequent evolution, and models that fit the nebular phase produce a diffusion peak that is too broad. The diffusion phase lightcurve is sensitive to the expansion of the "Ni bubbles", as the resulting Swiss-cheese-like geometry decreases the effective opacity and therefore the diffusion time. This effect has not been taken into account in previous lightcurve modelling of stripped-envelope SNe, which may lead to a systematic underestimate of their ejecta masses. It should be emphasized, though, that JEKYLL is limited to a geometry that is spherically symmetric on average, and large-scale asymmetries may also play a role. The relatively high initial mass found for the progenitor of SN 2020acat places it at the upper end of the mass distribution of Type IIb SN progenitors, and a single star origin can not be excluded., Comment: Accepted for publication by Astronomy and Astrophysics

Abridged - Stars with ZAMS masses between 140 and $260 M_\odot$ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN2018ibb is a H-poor SLSN at $z=0.166$ that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the NIR with 2-10m class telescopes. SN2018ibb radiated $>3\times10^{51} \rm erg$ during its evolution, and its bolometric light curve reached $>2\times10^{44} \rm erg\,s^{-1}$ at peak. The long-lasting rise of $>93$ rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source ($^{56}$Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions for their photometric and spectroscopic properties. SN2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25-44 $M_\odot$ of freshly nucleosynthesised $^{56}$Ni, pointing to the explosion of a metal-poor star with a He-core mass of 120-130 $M_\odot$ at the time of death. This interpretation is also supported by the tentative detection of [Co II]$\lambda$1.025$\mu$m, which has never been observed in any other PISN candidate or SLSN before. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN2018ibb by far the best candidate for being a PISN, to date., Comment: Accepted in A&A, the revised version includes a PISN rate estimate and an additional test with PISN models. 47 pages, main text 41 pages, 38 figures, 16 Tables

A growing number of supernovae (SNe) are now known to exhibit evidence for significant interaction with a dense, pre-existing, circumstellar medium (CSM). SNe Ibn comprise one such class that can be characterised by both rapidly evolving light curves and persistent narrow He I lines. The origin of such a dense CSM in these systems remains a pressing question, specifically concerning the progenitor system and mass-loss mechanism. In this paper, we present multi-wavelength data of the Type Ibn SN 2020nxt, including $HST$/STIS ultraviolet spectra. We fit the data with recently updated CMFGEN models designed to handle configurations for SNe Ibn. The UV coverage yields strong constraints on the energetics and, when combined with the CMFGEN models, offer new insight on potential progenitor systems. We find the most successful model is a $\lesssim4 {\rm M}_\odot$ helium star that lost its $\sim 1\,{\rm M}_\odot$ He-rich envelope in the years preceding core collapse. We also consider viable alternatives, such as a He white dwarf merger. Ultimately, we conclude at least some SNe Ibn do not arise from single, massive ($>30 {\rm M}_\odot$) Wolf-Rayet-like stars., Comment: 17 pages, 13 figures, 1 table, submitted to MNRAS

We present a sample of 14 hydrogen-rich superluminous supernovae (SLSNe II) from the Zwicky Transient Facility (ZTF) between 2018 and 2020. We include all classified SLSNe with peaks $M_{g}<-20$ mag and with observed \emph{broad} but not narrow Balmer emission, corresponding to roughly 20 per cent of all hydrogen-rich SLSNe in ZTF phase I. We examine the light curves and spectra of SLSNe II and attempt to constrain their power source using light-curve models. The brightest events are photometrically and spectroscopically similar to the prototypical SN 2008es, while others are found spectroscopically more reminiscent of non-superluminous SNe II, especially SNe II-L. $^{56}$Ni decay as the primary power source is ruled out. Light-curve models generally cannot distinguish between circumstellar interaction (CSI) and a magnetar central engine, but an excess of ultraviolet (UV) emission signifying CSI is seen in most of the SNe with UV data, at a wide range of photometric properties. Simultaneously, the broad H$\alpha$ profiles of the brightest SLSNe II can be explained through electron scattering in a symmetric circumstellar medium (CSM). In other SLSNe II without narrow lines, the CSM may be confined and wholly overrun by the ejecta. CSI, possibly involving mass lost in recent eruptions, is implied to be the dominant power source in most SLSNe II, and the diversity in properties is likely the result of different mass loss histories. Based on their radiated energy, an additional power source may be required for the brightest SLSNe II, however -- possibly a central engine combined with CSI., Comment: 27 pages including appendices; 14 figures. Revised version, resubmitted to MNRAS after the referee's comments

We use the new NLTE lightcurve and spectral synthesis code JEKYLL to evolve a macroscopically mixed ejecta model of a type IIb Supernova (SN) originating from a star with an initial mass of 12 solar masses through the photospheric and nebular phase. We compare to SN 2011dh, and find that both the spectra and the lightcurves are well reproduced. Our work further strengthens the evidence that this SN originated from a star with an initial mass of about 12 solar masses that had lost all but tiny (<0.1 solar masses) fraction of its hydrogen envelope. We also investigate the effects of the macroscopic mixing by comparing macroscopically and microscopically mixed models, and by varying the clumping geometry. In the photospheric phase, we find strong effects on the effective opacity in the macroscopically mixed regions, which affects the model lightcurves. The diffusion peak is considerably narrower in the macroscopically mixed case, and differs strongly if the radioactive material in the helium envelope is allowed to expand more than in our standard model. The effect is mainly geometrical, and is driven by the expansion of the clumps containing radioactive material. These findings has implications for lightcurve modelling of stripped-envelope SNe in general, and the effect would increase the estimated ejecta masses. In the nebular phase, we find strong effects on the collisional cooling rates in the macroscopically mixed regions, which affects lines driven by collisional cooling, in particular the [Ca II] 7291, 7323 and [O I] 6300, 6364 lines. As these lines are often used for mass determinations, it highlights the importance of how the calcium- and oxygen-rich material is mixed. As shown in this and earlier work, both NLTE and macroscopic mixing are essential ingredients to accurately model the lightcurves and spectra of Type IIb SNe throughout their evolution., Comment: Accepted for publication by Astronomy and Astrophysics

We present photometric and spectroscopic observations of the Type IIn supernova SN 2019zrk (also known as ZTF20aacbyec). The SN shows a $\gtrsim$ 100 day precursor, with a slow rise, followed by a rapid rise to M $\sim -19.2$ in the $r$ and $g$ bands. The post-peak light-curve decline is well fit with an exponential decay with a timescale of $\sim 39$ days, but it shows prominent undulations, with an amplitude of $\sim 1$ mag. Both the light curve and spectra are dominated by an interaction with a dense circumstellar medium (CSM), probably from previous mass ejections. The spectra evolve from a scattering-dominated Type IIn spectrum to a spectrum with strong P-Cygni absorptions. The expansion velocity is high, $\sim 16,000$ km s$^{-1}$, even in the last spectra. The last spectrum $\sim 110$ days after the main eruption reveals no evidence for advanced nucleosynthesis. From analysis of the spectra and light curves, we estimate the mass-loss rate to be $\sim 4 \times 10^{-2}$ M$_\odot$ yr$^{-1}$ for a CSM velocity of 100 km s$^{-1}$, and a CSM mass of $\gtrsim 1$ M$_\odot$. We find strong similarities for both the precursor, general light curve, and spectral evolution with SN 2009ip and similar SNe, although SN 2019zrk displays a brighter peak magnitude. Different scenarios for the nature of the 09ip-class of SNe, based on pulsational pair instability eruptions, wave heating, and mergers, are discussed. }, Comment: 18 pages, 12 figures. Astronomy & Astrophysics, in press

We present spectroscopy of the ejecta of SN 1987A in 2017 and 2018 from the Hubble Space Telescope and the Very Large Telescope, covering the wavelength range between $1150$ and $10000$ {\AA}. At 31 years, this is the first epoch with coverage over the ultraviolet-to-near-infrared range since 1995. We create velocity maps of the ejecta in the H$\alpha$, Mg II $\lambda\lambda2796,2804$ and [O I] $\lambda\lambda6302,6366$ (vacuum) emission lines and study their morphology. All three lines have a similar morphology, but Mg II is blueshifted by $\sim$1000 km s$^{-1}$ relative to the others and stronger in the northwest. We also study the evolution of the line fluxes, finding a brightening by a factor of $\sim$9 since 1999 in Mg II, while the other line fluxes are similar in 1999 and 2018. We discuss implications for the power sources of emission lines at late times: thermal excitation due to heating by the X-rays from the ejecta-ring interaction is found to dominate the ultraviolet Mg II lines, while the infrared Mg II doublet is powered mainly by Ly$\alpha$ fluorescence. The X-ray deposition is calculated based on merger models of SN 1987A. Far-ultraviolet emission lines of H$_2$ are not detected. Finally, we examine the combined spectrum of recently-discovered hotspots outside the equatorial ring. Their unresolved Balmer emission lines close to zero velocity are consistent with the interaction of fast ejecta and a clumpy, slowly moving outflow. A clump of emission in this spectrum, south of the equatorial ring at $\sim$1500 km s$^{-1}$, is likely associated with the reverse shock., Comment: 19 pages, 17 figures. Accepted for publication in MNRAS

Supernova 1987A offers a unique opportunity to study an evolving supernova in unprecedented detail over several decades. The X-ray emission is dominated by interactions between the ejecta and the circumstellar medium, primarily the equatorial ring (ER). We analyze 3.3 Ms of NuSTAR data obtained between 2012 and 2020, and two decades of XMM-Newton data. Since ${\sim}$2013, the flux below 2 keV has declined, the 3-8 keV flux has increased, but has started to flatten, and the emission above 10 keV has remained nearly constant. The spectra are well described by a model with three thermal shock components. Two components at 0.3 and 0.9 keV are associated with dense clumps in the ER, and a 4 keV component may be a combination of emission from diffuse gas in the ER and the surrounding low-density H II region. We disfavor models that involve non-thermal X-ray emission and place constraints on non-thermal components, but cannot firmly exclude an underlying power law. Radioactive lines show a $^{44}$Ti redshift of $670^{+520}_{-380}$ km s$^{-1}$, $^{44}$Ti mass of $1.73_{-0.29}^{+0.27}\times{}10^{-4}$ solar masses, and $^{55}$Fe mass of $<4.2\times{}10^{-4}$ solar masses. The 35-65 keV luminosity limit on the compact object is $2\times{}10^{34}$ erg s$^{-1}$, and $<15$% of the 10-20 keV flux is pulsed. Considering previous limits, we conclude that there are currently no indications of a compact object, aside from a possible hint of dust heated by a neutron star in recent ALMA images., Comment: 39 pages, 15 figures, 12 tables, no pulsar. Accepted for publication in ApJ. Updated to match the accepted version. We have made a number of clarifications and significantly expanded the discussion about instrumental uncertainties in Appendices C and D. The scientific conclusions remain unchanged

We carry out a comprehensive analysis of supernova \snem\ covering a period of 19 years since discovery. SN 2001em is the oldest supernova known to have undergone a metamorphosis from a stripped envelope to an interacting supernova. An early spectrum indicates it exploded as a Type Ib supernova. Later, the ejecta caught up with a dense circumstellar H-shell, ejected a few thousand years before the explosion, triggering interaction between the supernova ejecta and the dense shell, producing radio, X-ray and H$\alpha$ emission. We use archival data with the Very Large Array in radio bands and with Chandra, XMM-NEWTON and SWIFT-XRT in X-ray bands, along with published H$\alpha$ measurements. We combine these data with our low radio frequency observations with the Giant Metrewave Radio Telescope at two epochs covering three frequencies. While the observations missed the phase when the shock entered the dense shell, the X-rays indicate that the shock came out of the dense shell at around 1750 days. The data suggest a forward shock origin of the X-ray emission. Radio data show a spectral inversion at late epochs ($> 5000$\,d) at around 3 GHz, which mimics the properties of the central absorbed component seen in SN 1986J. A possible explanation for this component is that the progenitor of \snem\ was a massive binary system which had undergone a period of common envelope evolution. The hydrogen envelope from the \snem\ progenitor may have been lost as a result of binary interaction., Comment: Accepted for publication in ApJ

The Type IIn supernova (SN) 2005ip is one of the most well-studied and long-lasting examples of a SN interacting with its circumstellar environment. The optical light curve plateaued at a nearly constant level for more than five years, suggesting ongoing shock interaction with an extended and clumpy circumstellar medium (CSM). Here we present continued observations of the SN from $\sim 1000-5000$ days post-explosion at all wavelengths, including X-ray, ultraviolet, near-infrared, and mid-infrared. The UV spectra probe the pre-explosion mass loss and show evidence for CNO processing. From the bolometric light curve, we find that the total radiated energy is in excess of $10^{50}$ erg, the progenitor star's pre-explosion mass-loss rate was $\gtrsim 1 \times 10^{-2}\,{\rm M_{\odot}~ yr}^{-1}$, and the total mass lost shortly before explosion was $\gtrsim 1\,{\rm M_\odot}$, though the mass lost could have been considerably larger depending on the efficiency for the conversion of kinetic energy to radiation. The ultraviolet through near-infrared spectrum is characterised by two high density components, one with narrow high-ionisation lines, and one with broader low-ionisation H I, He I, [O I], Mg II, and Fe II lines. The rich Fe II spectrum is strongly affected by Ly$\alpha$ fluorescence, consistent with spectral modeling. Both the Balmer and He I lines indicate a decreasing CSM density during the late interaction period. We find similarities to SN 1988Z, which shows a comparable change in spectrum at around the same time during its very slow decline. These results suggest that, at long last, the shock interaction in SN 2005ip may finally be on the decline., Comment: 16 pages, 9 figures, 9 tables. Accepted to MNRAS

We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$\alpha$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit., Comment: 32 pages, containing 19 figures and three appendices

During the first few hundred days after the explosion, core-collapse supernovae (SNe) emit down-scattered X-rays and gamma-rays originating from radioactive line emissions, primarily from the $^{56}$Ni $\rightarrow$ $^{56}$Co $\rightarrow$ $^{56}$Fe chain. We use SN models based on three-dimensional neutrino-driven explosion simulations of single stars and mergers to compute this emission and compare the predictions with observations of SN 1987A. A number of models are clearly excluded, showing that high-energy emission is a powerful way of discriminating between models. The best models are almost consistent with the observations, but differences that cannot be matched by a suitable choice of viewing angle are evident. Therefore, our self-consistent models suggest that neutrino-driven explosions are able to produce, in principle, sufficient mixing, although remaining discrepancies may require small changes to the progenitor structures. The soft X-ray cutoff is primarily determined by the metallicity of the progenitor envelope. The main effect of asymmetries is to vary the flux level by a factor of ${\sim}$3. For the more asymmetric models, the shapes of the light curves also change. In addition to the models of SN 1987A, we investigate two models of Type II-P SNe and one model of a stripped-envelope Type IIb SN. The Type II-P models have similar observables as the models of SN 1987A, but the stripped-envelope SN model is significantly more luminous and evolves faster. Finally, we make simple predictions for future observations of nearby SNe., Comment: 27 pages, 12 figures. Accepted for publication in ApJ. Updated to match the accepted version; added two paragraphs in Section 2 on differences between the models/the level of mixing and rephrased the NuSTAR predictions slightly

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising ($1.4\pm0.1$ mag/hr) and luminous ($M_{g,\mathrm{peak}}=-20$ mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The high peak luminosity ($L_{\mathrm{bol}} \gtrsim 3 \times 10^{44}$ erg $\mathrm{sec}^{-1}$), the short rise time ($t_{\mathrm{rise}}= 3$ days in $g$-band), and the blue colors at peak ($g-r\sim-0.4$) all resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature ($T_{\mathrm{eff}}\gtrsim40,000$ K) spectra of a stripped-envelope SN. A retrospective search revealed luminous ($M_g \sim M_r \approx -14\,$mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release $E_\mathrm{\gamma,iso}<4.9 \times 10^{48}$ erg, a limit on X-ray emission $L_{\mathrm{X}} < 10^{40}\,$erg sec$^{-1}$, and a limit on radio emission $\nu L_\nu \lesssim 10^{37}\,$erg sec$^{-1}$. Taken together, we find that the early ($<10\,$days) data are best explained by shock breakout in a massive shell of dense circumstellar material ($0.02\,M_\odot$) at large radii ($3 \times 10^{14}\,$cm) that was ejected in eruptive pre-explosion mass-loss episodes. The late-time ($>10$ days) light curve requires an additional energy source, which could be the radioactive decay of Ni-56., Comment: Accepted to ApJ on 1 Oct 2019. In this version, we made minor corrections and removed extraneous references

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions., Comment: 14 pages including 3 figures

The material expelled by core-collapse supernova (SN) explosions absorbs X-rays from the central regions. We use SN models based on three-dimensional neutrino-driven explosions to estimate optical depths to the center of the explosion, compare different progenitor models, and investigate the effects of explosion asymmetries. The optical depths below 2 keV for progenitors with a remaining hydrogen envelope are expected to be high during the first century after the explosion due to photoabsorption. A typical optical depth is $100 t_4^{-2} E^{-2}$, where $t_4$ is the time since the explosion in units of 10 000 days (${\sim}$27 years) and $E$ the energy in units of keV. Compton scattering dominates above 50 keV, but the scattering depth is lower and reaches unity already at ${\sim}$1000 days at 1 MeV. The optical depths are approximately an order of magnitude lower for hydrogen-stripped progenitors. The metallicity of the SN ejecta is much higher than in the interstellar medium, which enhances photoabsorption and makes absorption edges stronger. These results are applicable to young SN remnants in general, but we explore the effects on observations of SN 1987A and the compact object in Cas A in detail. For SN 1987A, the absorption is high and the X-ray upper limits of ${\sim}$100 Lsun on a compact object are approximately an order of magnitude less constraining than previous estimates using other absorption models. The details are presented in an accompanying paper. For the central compact object in Cas A, we find no significant effects of our more detailed absorption model on the inferred surface temperature., Comment: 21 pages, 7 figures. Accepted for publication in ApJ. Updated to match accepted version; added Section 2.5 on asymmetries and discussion on homologous expansion in preamble of Section 4

Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun ($6\times 10^{-29}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in optical if our line-of-sight is free of ejecta dust, and $10^{36}$ erg s$^{-1}$ ($2\times 10^{-30}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in 2-10 keV X-rays. Our X-ray limits are an order of magnitude less constraining than previous limits because we use a more realistic ejecta absorption model based on three-dimensional neutrino-driven SN explosion models (presented in an accompanying article). The allowed bolometric luminosity of the compact object is 22 Lsun if our line-of-sight is free of ejecta dust, or 138 Lsun if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star to <4-8 MK and do not exclude models, which typically are in the range 3-4 MK. For the simplest accretion model, the accretion rate for an efficiency $\eta$ is limited to $< 10^{-11} \eta^{-1}$ Msun yr$^{-1}$, which excludes most predictions. For pulsar activity modeled by a rotating magnetic dipole in vacuum, the limit on the magnetic field strength ($B$) for a given spin period ($P$) is $B < 10^{14} P^2$ G s$^{-2}$. By combining information about radiation reprocessing and geometry, it is likely that the compact object is a dust-obscured thermally-emitting neutron star, which may appear as a region of higher-temperature ejecta dust emission., Comment: 35 pages, 10 figures. Accepted for publication in ApJ. Updated to match accepted version; slightly shortened and clarified the discussion

SN 2017dio shows both spectral characteristics of a type-Ic supernova (SN) and signs of a hydrogen-rich circumstellar medium (CSM). Prominent, narrow emission lines of H and He are superposed on the continuum. Subsequent evolution revealed that the SN ejecta are interacting with the CSM. The initial SN Ic identification was confirmed by removing the CSM interaction component from the spectrum and comparing with known SNe Ic, and reversely, adding a CSM interaction component to the spectra of known SNe Ic and comparing them to SN 2017dio. Excellent agreement was obtained with both procedures, reinforcing the SN Ic classification. The light curve constrains the pre-interaction SN Ic peak absolute magnitude to be around $M_g = -17.6$ mag. No evidence of significant extinction is found, ruling out a brighter luminosity required by a SN Ia classification. These pieces of evidence support the view that SN 2017dio is a SN Ic, and therefore the first firm case of a SN Ic with signatures of hydrogen-rich CSM in the early spectrum. The CSM is unlikely to have been shaped by steady-state stellar winds. The mass loss of the progenitor star must have been intense, $\dot{M} \sim 0.02$ $(\epsilon_{H\alpha}/0.01)^{-1}$ $(v_\textrm{wind}/500$ km s$^{-1}$) $(v_\textrm{shock}/10 000$ km s$^{-1})^{-3}$ $M_\odot$~yr$^{-1}$, peaking at a few decades before the SN. Such a high mass loss rate might have been experienced by the progenitor through eruptions or binary stripping., Comment: ApJL, 11 pages, 5 figures (accepted 22 Jan 2018)

It has become clear during the last decades that the interaction between the supernova ejecta and the circumstellar medium is playing a major role both for the observational properties of the supernova and for understanding the evolution of the progenitor star leading up to the explosion. In addition, it provides an opportunity to understand the shock physics connected to both thermal and non-thermal processes, including relativistic particle acceleration, radiation processes and the hydrodynamics of shock waves. This chapter has an emphasis on the information we can get from radio and X-ray observations, but also their connection to observations in the optical and ultraviolet. We first review the different physical processes involved in circumstellar interaction, including hydrodynamics, thermal X-ray emission, acceleration of relativistic particles and non-emission processes in the radio and X-ray ranges. Finally, we discuss applications of these to different types of supernovae., Comment: 63 pages, 18 figures. Submitted as Chapter for the 'Supernova Handbook' to be published by Springer (eds. A. Alsabti and P. Murdin). Major update of the review in 'Supernovae and Gamma-Ray Bursters', Ed. K. Weiler., Springer Lecture Notes in Physics, vol. 598, p.171-194 Springer (arXiv:astro-ph/0110060). Comments welcome

Both CO and SiO have been observed at early and late phases in SN 1987A. H_2 was predicted to form at roughly the same time as these molecules, but was not detected at early epochs. Here we report the detection of NIR lines from H_2 at 2.12 mu and 2.40 mu in VLT/SINFONI spectra obtained between days 6489 and 10,120. The emission is concentrated to the core of the supernova in contrast to H-alpha and approximately coincides with the [Si I]/[Fe II] emission detected previously in the ejecta. Different excitation mechanisms and power sources of the emission are discussed. From the nearly constant H_2 luminosities we favour excitation resulting from the 44Ti decay., Comment: 5 pages 3 figures, accepted for ApJ Letters

The observational effects of the 'Infrared Catastrophe' are discussed in view of the very late observations of the Type Ia SN 2011fe. Our model spectra at 1000d take non-local radiative transfer into account, and find that this has a crucial impact on the spectral formation. Although rapid cooling of the ejecta to a few 100 K occurs also in these models, the late-time optical/NIR flux is brighter by 1-2 magnitudes due to redistribution of UV emissivity, resulting from non-thermal excitation and ionization. This effect brings models into better agreement with late-time observations of SN 2011fe and other Type Ia supernovae, and offers a solution to the long standing discrepancy between models and observations. The models show that spectral formation shifts from Fe II and Fe III at 300d to Fe I at 1000d, which explains the apparent wavelength shifts seen in SN2011fe. We discuss effects of time dependence and energy input from 57Co, finding both to be important at 1000d., Comment: 5 pages, 5 figures. Accepted for publication in ApJ Letters

Searches for circumstellar material around Type Ia supernovae (SNe Ia) are one of the most powerful tests of the nature of SN Ia progenitors, and radio observations provide a particularly sensitive probe of this material. Here we report radio observations for SNe Ia and their lower-luminosity thermonuclear cousins. We present the largest, most sensitive, and spectroscopically diverse study of prompt (delta t <~ 1 yr) radio observations of 85 thermonuclear SNe, including 25 obtained by our team with the unprecedented depth of the Karl G. Jansky Very Large Array. With these observations, SN 2012cg joins SN 2011fe and SN 2014J as a SN Ia with remarkably deep radio limits and excellent temporal coverage (six epochs, spanning 5--216 days after explosion, yielding Mdot/v_w <~ 5 x 10^-9 M_sun/yr / (100 km/s), assuming epsilon_B = 0.1 and epsilon_e = 0.1). All observations yield non-detections, placing strong constraints on the presence of circumstellar material. We present analytical models for the temporal and spectral evolution of prompt radio emission from thermonuclear SNe as expected from interaction with either wind-stratified or uniform density media. These models allow us to constrain the progenitor mass loss rates, with limits ranging from Mdot <~ 10^-9--10^-4 M_sun/yr, assuming a wind velocity v_w=100 km/s. We compare our radio constraints with measurements of Galactic symbiotic binaries to conclude that <~10% of thermonuclear SNe have red giant companions., Comment: Submitted to ApJ

We present all X-ray and radio observations of the Type IIn supernova SN 2010jl. The X-ray observations cover a period up to day 1500 with Chandra, XMM-Newton, NuSTAR and Swift-XRT. The Chandra observations after 2012 June, the XMM-Newton observation in 2013 November, and most of the Swift-XRT observations until 2014 December are presented for the first time. All the spectra can be fitted by an absorbed hot thermal model except for \chandra spectra on 2011 October and 2012 June when an additional component is needed. Although the origin of this component is uncertain, it is spatially coincident with the supernova and occurs when there are changes to the supernova spectrum in the energy range close to that of the extra component, indicating that the emission is related to the supernova. The X-ray light curve shows an initial plateau followed by a steep drop starting at day $\sim 300$. We attribute the drop to a decrease in the circumstellar density. The column density to the X-ray emission drops rapidly with time, showing that the absorption is in the vicinity of the supernova. We also present Very Large Array radio observations of SN 2010jl. Radio emission was detected from \sn from day 570 onwards. The radio light curves and spectra suggest that the radio luminosity was close to its maximum at the first detection. The velocity of the shocked ejecta derived assuming synchrotron self absorption is much less than that estimated from the optical and X-ray observations, suggesting that free-free absorption dominates., Comment: Accepted for publication in The Astrophysical Journal. 21 pages (in ApJ format), 6 tables, 11 figures, Typos corrected

We present imaging and spectroscopic observations with HST and VLT of the ring of SN 1987A from 1994 to 2014. After an almost exponential increase of the shocked emission from the hotspots up to day ~8,000 (~2009), both this and the unshocked emission are now fading. From the radial positions of the hotspots we see an acceleration of these up to 500-1000 km/s, consistent with the highest spectroscopic shock velocities from the radiative shocks. In the most recent observations (2013 and 2014), we find several new hotspots outside the inner ring, excited by either X-rays from the shocks or by direct shock interaction. All of these observations indicate that the interaction with the supernova ejecta is now gradually dissolving the hotspots. We predict, based on the observed decay, that the inner ring will be destroyed by ~2025., Comment: 5 pages, 4 figures, accepted for ApJ Letters

We present radio and X-ray observations of the nearby Type IIb Supernova 2013df in NGC4414 from 10 to 250 days after the explosion. The radio emission showed a peculiar soft-to-hard spectral evolution. We present a model in which inverse Compton cooling of synchrotron emitting electrons can account for the observed spectral and light curve evolution. A significant mass loss rate, $\dot{M} \approx 8 \times 10^{-5}\,\rm M_{\odot}/yr$ for a wind velocity of 10 km/s, is estimated from the detailed modeling of radio and X-ray emission, which are primarily due to synchrotron and bremsstrahlung, respectively. We show that SN 2013df is similar to SN 1993J in various ways. The shock wave speed of SN 2013df was found to be average among the radio supernovae; $v_{sh}/c \sim 0.07$. We did not find any significant deviation from smooth decline in the light curve of SN 2013df. One of the main results of our self-consistent multiband modeling is the significant deviation from energy equipartition between magnetic fields and relativistic electrons behind the shock. We estimate $\epsilon_{e} = 200 \epsilon_{B}$. In general for Type IIb SNe, we find that the presence of bright optical cooling envelope emission is linked with free-free radio absorption and bright thermal X-ray emission. This finding suggests that more extended progenitors, similar to that of SN 2013df, suffer from substantial mass loss in the years before the supernova., Comment: 15 pages, 7 figures, 1 table; Submitted to The Astrophysical Journal

We present new {\it Hubble Space Telescope} images of high-velocity H-$\alpha$ and Lyman-$\alpha$ emission in the outer debris of SN~1987A. The H-$\alpha$ images are dominated by emission from hydrogen atoms crossing the reverse shock. For the first time we observe emission from the reverse shock surface well above and below the equatorial ring, suggesting a bipolar or conical structure perpendicular to the ring plane. Using the H$\alpha$ imaging, we measure the mass flux of hydrogen atoms crossing the reverse shock front, in the velocity intervals ($-$7,500~$<$~$V_{obs}$~$<$~$-$2,800 km s$^{-1}$) and (1,000~$<$~$V_{obs}$~$<$~7,500 km s$^{-1}$), $\dot{M_{H}}$ = 1.2~$\times$~10$^{-3}$ M$_{\odot}$ yr$^{-1}$. We also present the first Lyman-$\alpha$ imaging of the whole remnant and new $Chandra$ X-ray observations. Comparing the spatial distribution of the Lyman-$\alpha$ and X-ray emission, we observe that the majority of the high-velocity Lyman-$\alpha$ emission originates interior to the equatorial ring. The observed Lyman-$\alpha$/H-$\alpha$ photon ratio, $\langle$$R(L\alpha / H\alpha)$$\rangle$ $\approx$~17, is significantly higher than the theoretically predicted ratio of $\approx$ 5 for neutral atoms crossing the reverse shock front. We attribute this excess to Lyman-$\alpha$ emission produced by X-ray heating of the outer debris. The spatial orientation of the Lyman-$\alpha$ and X-ray emission suggests that X-ray heating of the outer debris is the dominant Lyman-$\alpha$ production mechanism in SN 1987A at this phase in its evolution., Comment: 6 pages, 5 figures. ApJL - accepted

SN 2012ec is a Type IIP supernova (SN) with a progenitor detection and comprehensive photospheric-phase observational coverage. Here, we present Very Large Telescope and PESSTO observations of this SN in the nebular phase. We model the nebular [O I] 6300, 6364 lines and find their strength to suggest a progenitor main-sequence mass of 13-15 Msun. SN 2012ec is unique among hydrogen-rich SNe in showing a distinct and unblended line of stable nickel [Ni II] 7378. This line is produced by 58Ni, a nuclear burning ash whose abundance is a sensitive tracer of explosive burning conditions. Using spectral synthesis modelling, we use the relative strengths of [Ni II] 7378 and [Fe II] 7155 (the progenitor of which is 56Ni) to derive a Ni/Fe production ratio of 0.20pm0.07 (by mass), which is a factor 3.4pm1.2 times the solar value. High production of stable nickel is confirmed by a strong [Ni II] 1.939 micron line. This is the third reported case of a core-collapse supernova producing a Ni/Fe ratio far above the solar value, which has implications for core-collapse explosion theory and galactic chemical evolution models., Comment: Published version