Neutron star (NS) mergers are known to produce heavy elements through rapid neutron capture (r-process) nucleosynthesis. Actinides are expected to be created solely by the r-process in the most neutron-rich environments. Confirming if NS mergers provide the requisite conditions for actinide creation is therefore central to determining their origin in the Universe. Actinide signatures in kilonova (KN) spectra may yield an answer, provided adequate models are available in order to interpret observational data. In this study, we investigate actinide signatures in neutron-rich merger ejecta. We use three ejecta models with different compositions and radioactive power, generated by nucleosynthesis calculations using the same initial electron fraction (Ye=0.15) but with different nuclear physics inputs and thermodynamic expansion history. These are evolved from 10 to 100 d after merger using the sumo non-local thermodynamic equilibrium (NLTE) radiative transfer code. We highlight how uncertainties in nuclear properties, as well as choices in thermodynamic trajectory, may yield entirely different outputs for equal values of Ye. We consider an actinide-free model and two actinide-rich models, and find that the emergent spectra and light-curve evolution are significantly different depending on the amount of actinides present, and the overall decay properties of the models. We also present potential key actinide spectral signatures, of which doubly ionized 89Ac and 90Th may be particularly interesting as spectral indicators of actinide presence in KN ejecta.

Stars are initially powered by the fusion of hydrogen to helium. These ashes serve as fuel in a series of stages1, 2–3, transforming massive stars into a structure of shells. These are composed of natal hydrogen on the outside and consecutively heavier compositions inside, predicted to be dominated by He, C/O, O/Ne/Mg and O/Si/S (refs. 4,5). Silicon and sulfur are fused into iron, leading to the collapse of the core and either a supernova explosion or the formation of a black hole6, 7, 8–9. Stripped stars, in which the outer hydrogen layer has been removed and the internal He-rich or even the C/O layer below it is exposed10, provide evidence for this shell structure and the cosmic element production mechanism it reflects. The supernova types that arise from stripped stars embedded in shells of circumstellar material (CSM) confirm this scenario11, 12, 13, 14–15. However, direct evidence for the most interior shells, which are responsible for producing elements heavier than oxygen, is lacking. Here we report the discovery of the supernova (SN) 2021yfj resulting from a star stripped to its O/Si/S-rich layer. We directly observe a thick, massive Si/S-rich shell, expelled by the progenitor shortly before the supernova explosion. Exposing such an inner stellar layer is theoretically challenging and probably requires a rarely observed mass-loss mechanism. This rare supernova event reveals advanced stages of stellar evolution, forming heavier elements, including silicon, sulfur and argon, than those detected on the surface of any known class of massive stars.

e present supernova (SN) 2023ufx, a unique Type IIP SN with the shortest known plateau duration (t_PT ∼ 47 days), a luminous V-band peak (M_V = −​​​​​​18.42 ± 0.08 mag), and a rapid early decline rate (s1 = 3.47 ± 0.09 mag (50 days)⁻¹). By comparing observed photometry to a hydrodynamic MESA+STELLA model grid, we constrain the progenitor to be a massive red supergiant with M_ZAMS ∼ 19–25 M_⊙. Independent comparisons with nebular spectral models also suggest an initial He-core mass of ∼6 M_⊙, and thus a massive progenitor. For a Type IIP, SN 2023ufx produced an unusually high amount of nickel (⁵⁶Ni) ∼0.14 ± 0.02 M_⊙, during the explosion. We find that the short plateau duration in SN 2023ufx can be explained with the presence of a small hydrogen envelope (M_Henv ∼ 1.2 M_⊙), suggesting partial stripping of the progenitor. About ∼0.09 M_⊙ of circumstellar material through mass loss from late-time stellar evolution of the progenitor is needed to fit the early time (≲10 days) pseudo-bolometric light curve. Nebular line diagnostics of broad and multipeak components of [O i] λλ6300, 6364, Hα, and [Ca ii] λλ7291, 7323 suggest that the explosion of SN 2023ufx could be inherently asymmetric, preferentially ejecting material along our line of sight.

To investigate spectra of kilonovae in the nonlocal thermal equilibrium phase (t ≳ 1 week), we perform atomic calculations for dielectronic recombination (DR) rates for the light r-process elements Se (Z = 34), Rb (Z = 37), Sr (Z = 38), Y (Z = 39), and Zr (Z = 40) using the HULLAC code. For the different elements, our results for the DR rate coefficients for recombining from the ionization states of II to I, III to II, and IV to III vary between 2 × 10⁻¹²–5 × 10⁻¹¹ cm³ s⁻¹, 10⁻¹³–5 × 10⁻¹¹ cm³ s⁻¹, and 2 × 10⁻¹⁵–10⁻¹¹ cm³ s⁻¹, respectively, at a temperature of T = 10,000 K. Using this new atomic data (DR), we study the impact on kilonova model spectra at phases of t = 10 days and t = 25 days after the merger using the spectral synthesis code SUMO. Compared to models using the previous treatment of recombination as a constant rate, the new models show significant changes in ionization and temperature, and, correspondingly, in emergent spectra. With the new rates, we find that Zr (Z = 40) plays a yet more dominant role in kilonova spectra for light r-process compositions. Furthermore, we show that previously predicted mid-infrared (e.g., [Se III] 4.55 μm) and optical (e.g., Rb I 7802, 7949 Å) lines weaken in the new model. Instead, [Se I] 5.03 μm emerges as a signature. These results demonstrate the importance of considering the detailed microphysics for modelling and interpreting the late-time kilonova spectra.

Aims. Stripped envelope (SE) supernovae are explosions of stars that have somehow lost most of their outer envelopes. We present the discovery and analyse the observations of the Type Ib supernova 2019odp (a.k.a. ZTF19abqwtfu) covering epochs within days of the explosion to late nebular phases at 360 d post-explosion.Methods. Our observations include an extensive set of photometric observations and low- to medium-resolution spectroscopic observations, both covering the complete observable time range. We analysed the data using analytic models for the recombination cooling emission of the early excess emission and the diffusion of the peak light curve. We expanded on existing methods to derive oxygen mass estimates from nebular phase spectroscopy, and briefly discuss progenitor models based on this analysis.Results. Our spectroscopic observations confirm the presence of He in the supernova ejecta and we thus (re)classify SN 2019odp as a Type Ib supernova. From the pseudo-bolometric light curve, we estimate a high ejecta mass of Mej ∼ 4 − 7 M⊙. The high ejecta mass, large nebular [O I]/[Ca II] line flux ratio (1.2 − 1.9), and an oxygen mass above ⪆0.5 M⊙ point towards a progenitor with a pre-explosion mass higher than 18 M⊙. Whereas a majority of analysed SE supernovae in the literature seem to have low ejecta masses, indicating stripping in a binary star system, SN 2019odp instead has parameters that are consistent with an origin in a single massive star. The compact nature of the progenitor (≲10 R⊙) suggests that a Wolf-Rayet star is the progenitor.

Stars with zero-age main sequence masses between 140 and 260 M_⊙ 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. SN 2018ibb is a hydrogen-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 near-infrared (NIR) with 2–10 m class telescopes. SN 2018ibb radiated > 3 × 10⁵¹ erg during its evolution, and its bolometric light curve reached > 2 × 10⁴⁴ erg s⁻¹ at its 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 (⁵⁶Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions as to their photometric and spectroscopic properties. SN 2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, and potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25–44 M_⊙ of freshly nucleosynthesised ⁵⁶Ni, pointing to the explosion of a metal-poor star with a helium core mass of 120–130 M_⊙ at the time of death. This interpretation is also supported by the tentative detection of [Co II] λ 1.025 μm, which has never been observed in any other PISN candidate or SLSN before. We observe a significant excess in the blue part of the optical spectrum during the nebular phase, which is in tension with predictions of existing PISN models. However, we have compelling observational evidence for an eruptive mass-loss episode of the progenitor of SN 2018ibb shortly before the explosion, and our dataset reveals that the interaction of the SN ejecta with this oxygen-rich circumstellar material contributed to the observed emission. That may explain this specific discrepancy with PISN models. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN 2018ibb by far the best candidate for being a PISN, to date.

SN 2019neq was a very fast evolving superluminous supernova. At a redshift z = 0.1059, its peak absolute magnitude was −21.5 ± 0.2 mag in g band. In this work, we present data and analysis from an extensive spectrophotometric follow-up campaign using multiple observational facilities. Thanks to a nebular spectrum of SN 2019neq, we investigated some of the properties of the host galaxy at the location of SN 2019neq and found that its metallicity and specific star formation rate are in a good agreement with those usually measured for SLSNe-I hosts. We then discuss the plausibility of the magnetar and the circumstellar interaction scenarios to explain the observed light curves, and interpret a nebular spectrum of SN 2019neq using published SUMO radiative-transfer models. The results of our analysis suggest that the spin-down radiation of a millisecond magnetar with a magnetic field B ≃ 6×10¹⁴ G could boost the luminosity of SN 2019neq.

Understanding the explosion mechanism and hydrodynamic evolution of core-collapse supernovae (SNe) is a long-standing quest in astronomy. The asymmetries caused by the explosion are encoded into the line profiles which appear in the nebular phase of the SN evolution – with particularly clean imprints in He star explosions. Here, we carry out nine different supernova simulations of He-core progenitors, exploding them in 3D with parametrically varied neutrino luminosities using the PROMETHEUS-HOTB code, hydrodynamically evolving the models to the homologous phase. We then compute nebular phase spectra with the 3D Non-Local Thermodynamic Equilibrium spectral synthesis code EXTRASS (EXplosive TRAnsient Spectral Simulator). We study how line widths and shifts depend on progenitor mass, explosion energy, and viewing angle. We compare the predicted line profile properties against a large set of Type Ib observations, and discuss the degree to which current neutrino-driven explosions can match observationally inferred asymmetries. With self-consistent 3D modelling – circumventing the difficulties of representing 56Ni mixing and clumping accurately in 1D models – we find that neither low-mass He cores exploding with high energies nor high-mass cores exploding with low energies contribute to the Type Ib SN population. Models which have line profile widths in agreement with this population give sufficiently large centroid shifts for calcium emission lines. Calcium is more strongly affected by explosion asymmetries connected to the neutron star kicks than oxygen and magnesium. Lastly, we turn to the near-infrared spectra from our models to investigate the potential of using this regime to look for the presence of He in the nebular phase.

We present JWST observations of the Crab Nebula, the iconic remnant of the historical SN 1054. The observations include NIRCam and MIRI imaging mosaics plus MIRI/MRS spectra that probe two select locations within the ejecta filaments. We derive a high-resolution map of dust emission and show that the grains are concentrated in the innermost, high-density filaments. These dense filaments coincide with multiple synchrotron bays around the periphery of the Crab's pulsar wind nebula (PWN). We measure synchrotron spectral index changes in small-scale features within the PWN’s torus region, including the well-known knot and wisp structures. The index variations are consistent with Doppler boosting of emission from particles with a broken power-law distribution, providing the first direct evidence that the curvature in the particle injection spectrum is tied to the acceleration mechanism at the termination shock. We detect multiple nickel and iron lines in the ejecta filaments and use photoionization models to derive nickel-to-iron abundance ratios that are a factor of 3-8 higher than the solar ratio. We also find that the previously reported order-of-magnitude higher Ni/Fe values from optical data are consistent with the lower values from JWST when we reanalyze the optical emission using updated atomic data and account for local extinction from dust. We discuss the implications of our results for understanding the nature of the explosion that produced the Crab Nebula and conclude that the observational properties are most consistent with a low-mass Fe core-collapse supernova, even though an electron-capture explosion cannot be ruled out.

Multipeaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility. Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (−17.4 mag) and long (∼100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew's spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak's photospheric phase, before Hα resurfaces again during the nebular phase. The nebular lines ([O i], [Ca ii], Mg i], Hα) exhibit a double-peaked structure that hints toward a clumpy or nonspherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than that of normal SESNe as well as requiring a very high ⁵⁶Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak and also powers the second peak of the light curve through interaction of the SN with the circumstellar medium.