Core-collapse Supernovae Research Articles

Enhanced production of 60Fe in massive stars.

Massive stars are a major source of chemical elements in the cosmos, ejecting freshly produced nuclei through winds and core-collapse supernova explosions into the interstellar medium. Among the material ejected, long-lived radioisotopes, such as 60Fe (iron) and 26Al (aluminum), offer unique signs of active nucleosynthesis in our galaxy. There is a long-standing discrepancy between the observed 60Fe/26Al ratio by γ-ray telescopes and predictions from supernova models. This discrepancy has been attributed to uncertainties in the nuclear reaction networks producing 60Fe, and one reaction in particular, the neutron-capture on 59Fe. Here we present experimental results that provide a strong constraint on this reaction. We use these results to show that the production of 60Fe in massive stars is higher than previously thought, further increasing the discrepancy between observed and predicted 60Fe/26Al ratios. The persisting discrepancy can therefore not be attributed to nuclear uncertainties, and points to issues in massive-star models.

Preferential occurrence of fast radio bursts in massive star-forming galaxies.

Fast radio bursts (FRBs) are millisecond-duration events detected from beyond the Milky Way. FRB emission characteristics favour highly magnetized neutron stars, or magnetars, as the sources1, as evidenced by FRB-like bursts from a galactic magnetar2,3, and the star-forming nature of FRB host galaxies4,5. However, the processes that produce FRB sources remain unknown6. Although galactic magnetars are often linked to core-collapse supernovae (CCSNe)7, it is uncertain what determines which supernovae result in magnetars. The galactic environments of FRB sources can be used to investigate their progenitors. Here, we present the stellar population properties of 30 FRB host galaxies discovered by the Deep Synoptic Array (DSA-110). Our analysis shows a marked deficit of low-mass FRB hosts compared with the occurrence of star formation in the Universe, implying that FRBs are a biased tracer of star formation, preferentially selecting massive star-forming galaxies. This bias may be driven by galaxy metallicity, which is positively correlated with stellar mass8. Metal-rich environments may favour the formation of magnetar progenitors through stellar mergers9,10, as higher-metallicity stars are less compact and more likely to fill their Roche lobes, leading to unstable mass transfer. Although massive stars do not have convective interiors to generate strong magnetic fields by dynamo11, merger remnants are thought to have the requisite internal magnetic-field strengths to result in magnetars11,12. The preferential occurrence of FRBs in massive star-forming galaxies suggests that a core-collapse supernova of merger remnants preferentially forms magnetars.

Identifying Jittering Jet-shaped Ejecta in the Cygnus Loop Supernova Remnant

Analyzing images of the Cygnus Loop, a core-collapse supernova (CCSN) remnant, in different emission bands, we identify a point-symmetrical morphology composed of three symmetry axes that we attribute to shaping by three pairs of jets. The main jet axis has an elongated S shape, appearing as a faint narrow zone in visible and UV. We term it the S-shaped hose, and the structure of three symmetry lines the point-symmetric wind rose. The two other lines connect a protrusion (an ear or a bulge) with a hole on the opposite side of the center (a nozzle or a cavity), structures that we identify in the X-ray, UV, visible, IR, and/or radio images. There is a well-known blowout at the southern end of the S-shaped hose, and we identify a possible opposite blowout at the northern end of the S-shaped hose. The point-symmetrical morphology of the Cygnus Loop is according to the expectation of the jittering jets explosion mechanism (JJEM) of CCSNe, where several to few tens of pairs of jets with stochastically varying directions explode the star. The three pairs of jets that shaped the wind-rose structure of the Cygnus Loop are the last energetic pairs of this series of jets. Our study further supports the JJEM as the main explosion mechanism of CCSNe.

Generating tailored high frequency features in core collapse supernova gravitational wave signals applicable in LIGO interferometric studies

In this article, we introduce a methodology based on an analytical model of a damped harmonic oscillator subject to random forcing to generate transient gravitational wave signals. Such a model incorporates a simulated linear high-frequency component that mirrors the growing characteristic frequency over time observed in numerical simulations of core-collapse supernova gravitational wave signals. Unlike traditional numerical simulations, the method proposed in this study requires minimal computational resources, which makes it particularly advantageous for tasks such as data analysis, detection, and reconstruction of gravitational wave transients. To verify the physical accuracy of the generated signals, they are compared against the amplitude spectral of current LIGO interferometers and a 3D numerical simulation of a core-collapse supernova gravitational wave signal from the Andresen et al. 2017 model s15.nr. The results indicate that this approach is effective in generating scalable signals that align with LIGO interferometric data, offering potential utility in various gravitational wave transient investigations.

The next galactic supernova can uncover mass and couplings of particles decaying to neutrinos

Many particles predicted by extensions of the Standard Model feature interactions with neutrinos, e.g., Majoron-like bosons ϕ. If the mass of ϕ is larger than about 10 keV, they can be produced abundantly in the core of the next galactic core-collapse supernova through neutrino coalescence, and leave it with energies of around 100 MeV. Their subsequent decay to high-energy neutrinos and anti-neutrinos provides a distinctive signature at Earth. Ongoing and planned neutrino and dark matter experiments allow us to reconstruct the energy, flavor, and time of arrival of these high-energy neutrinos. For the first time, we show that these measurements can help pinpointing the mass of ϕ and its couplings to neutrinos of different flavor. Our results can be generalized in a straightforward manner to other hypothetical feebly interacting particles, like novel gauge bosons or heavy neutral leptons, that decay into neutrinos.

A census of the Sun’s ancestors and their contributions to the Solar System chemical composition

In this work, we compute the rates and numbers of different types of stars and phenomena (supernovae, novae, white dwarfs, merging neutron stars, black holes) that contributed to the chemical composition of the Solar System. During the Big Bang, only light elements formed, while all the heavy ones, from carbon to uranium and beyond, have since been created inside stars. Stars die and release the newly formed elements into the interstellar gas. This process is called ‘chemical evolution’. In particular, we analyse the death rates of stars of all masses, whether they die quiescently or explosively. These rates and total star numbers are computed in the context of a revised version of the two-infall model for the chemical evolution of the Milky Way, which reproduces the observed abundance patterns of several chemical species, the global solar metallicity, and the current gas, stellar, and total surface mass densities relatively well. We also compute the total number of stars ever born and still alive as well as the number of stars born up to the formation of the Solar System with mass and metallicity like those of the Sun. This latter number accounts for all the possible existing Solar systems that can host life in the solar vicinity. We conclude that, among all the stars (from 0.8 to 100 M⊙) that were born and died from the Big Bang up until the Solar System formation epoch and that contributed to its chemical composition, 93.00% were stars that died as single white dwarfs (without interacting significantly with a companion star) and originated in the mass range of 0.8–8 M⊙, while 5.24% were neutron stars and 0.73% were black holes, both originating from core-collapse supernovae (M > 8 M⊙); 0.64% were Type Ia supernovae and 0.40% were nova systems, both originating from the same mass range as the white dwarfs. The number of stars similar to the Sun born from the Big Bang up until the formation of the Solar System, with metallicity in the range 12+log(Fe/H)= 7.50 ± 0.04 dex, is ~31•107, and in particular our Sun is the ~2.61• 107-th star of this kind.

SDSS J0018-0939: A Clear Signature of Sub-Chandrasekhar Mass Type 1a Supernova

Very metal-poor (VMP) stars ([Fe/H] ≤ −2) that have subsolar values of [X/Fe] for α elements such as Mg, Si, and Ca, are referred to as α-poor VMP stars. They are quite rare among VMP stars and are thought to have formed from gas enriched predominantly by a single Type Ia supernova (SN 1a) in contrast to most VMP stars, which are α-enhanced and usually associated with core-collapse supernovae. The observed abundance pattern in such stars can provide a direct way to probe the nucleosynthesis in individual SN 1a. Although the abundance patterns in some α-poor VMP stars have been shown to be consistent with SN 1a ejecta, a clear nucleosynthetic signature for SN 1a resulting from the explosion of a near Chandrasekhar mass (near-M Ch) or a sub-Chandrasekhar mass (sub-M Ch) white dwarf, has not been unambiguously detected. We perform a detailed analysis of various formation channels of VMP stars and find that the α-poor VMP star SDSS J0018-0939, which was earlier reported as a star with potential pair-instability supernova origin, provides a near-smoking-gun signature of a sub-M Ch SN 1a resulting from He detonation. We find that compared to other α-poor VMP stars that were previously identified with SN 1a, SDSS J0018-0939 is the only star that has a clear and unambiguous signature of SN 1a. Interestingly, our results are consistent with constraints on SN1 a from recent galactic chemical evolution studies that indicate that sub-M Ch SN 1a accounts for ∼50%–75% of all SN 1a and is possibly the dominant channel in the early Galaxy.

Pulsar-wind nebulae meeting the circumstellar media of their progenitors

A significative fraction of high-mass stars sail away through the interstellar medium of the galaxies. Once they evolved and died via a core-collapse supernova, a magnetised, rotating neutron star (a pulsar) is usually left over. The immediate surroundings of the pulsar is the pulsar wind, which forms a nebula whose morphology is shaped by the supernova ejecta and channelled into the circumstellar medium of the progenitor star in the pre-supernova time. Irregular pulsar-wind nebulae display a large variety of radio appearances, screened by their interacting supernova blast wave, or harbour asymmetric up–down emission. Here, we present a series of 2.5-dimensional (2 dimensions for the scalar quantities plus a toroidal component for the vectors) black non-relativistic magneto-hydrodynamical simulations exploring the evolution of the pulsar-wind nebulae generated by a red supergiant and a Wolf-Rayet massive supernova progenitor, moving with Mach number $M=1$ and $M=2$ into the warm phase of the Galactic plane. black In such a simplified approach, the progenitor's direction of motion, the local ambient medium magnetic field, and the progenitor and pulsar axis of rotation, are all aligned; this restricted our study to peculiar pulsar-wind nebula of high-equatorial-energy flux. We find that the reverberation of the termination shock of the pulsar-wind nebulae, when sufficiently embedded into its dead stellar surroundings and interacting with the supernova ejecta, is asymmetric and differs greatly as a function of the past circumstellar evolution of its progenitor, which reflects into their projected radio synchrotron emission. This mechanism is particularly at work in the context of remnants involving slowly moving or very massive stars. We find that the mixing of material in plerionic core-collapse supernova remnants is strongly affected by the asymmetric reverberation in their pulsar-wind nebulae.

Lower-mass-gap Black Holes in Dense Star Clusters

The existence of compact stellar remnants in the mass range 2–5 M ⊙ has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5 M ⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte Carlo N-body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc−3 yr−1.

Convective vortices in collapsing stars

Abstract Recent studies show that non-radial structures arising from massive star shell convection play an important role in shaping core-collapse supernova explosions. During the collapse phase, convective vortices generate acoustic waves that interact with the supernova shock. This amplifies turbulence in the post-shock region, contributing to explosion. We study how various physical parameters influence the evolution of these convective vortices during stellar collapse using simplified simulations. We model the collapsing star with a transonic Bondi flow and represent convection as solenoidal velocity perturbations. Our results are consistent with previous studies, demonstrating that the peak perturbation amplitude scales linearly with the pre-collapse convective Mach number and inversely with the angular wavenumber of convection. While the shell radius and width primarily determine the timescale of accretion, they have little impact on the peak perturbation amplitudes. Finally, we show that when the convective Mach number is below approximately 0.2, the dynamics remain within the linear regime.

Mesiri: Mephisto Early Supernovae Ia Rapid Identifier

The early time observations of Type Ia supernovae (SNe Ia) play a crucial role in investigating and resolving longstanding questions about progenitor stars and the explosion mechanisms of these events. Colors of supernovae (SNe) in the initial days after the explosion can help differentiate between different types of SNe. However, the use of true color information to identify SNe Ia at the early-time explosion is still in its infancy. The Multi-channel Photometric Survey Telescope (Mephisto) is a photometric survey telescope equipped with three CCD cameras, capable of simultaneously imaging the same patch of sky in three bands (u, g, i or v, r, z), yielding real-time colors of astronomical objects. In this paper, we introduce a new time-series classification tool named Mephisto Early Supernovae Ia Rapid Identifier (Mesiri), which, for the first time, utilizes real-time color information to distinguish early-time SNe Ia from core-collapse supernovae. Mesiri is based on the deep learning approach and can achieve an accuracy of 96.75% ± 0.79%, and AUC of 98.87% ± 0.53% in case of single epoch random observation before the peak brightness. These values reach towards perfectness if additional data points on several night observations are considered. The classification with real-time color significantly outperforms that with pseudo-color, especially at the early time, i.e., with only a few points of observations. The BiLSTM architecture shows the best performance compared to others that have been tested in this work.

A Possible Formation Scenario of the Gaia BH1: Inner Binary Merger in Triple Systems

Based on astrometric measurements and spectral analysis from Gaia DR3, two quiescent black hole (BH) binaries, Gaia BH1 and BH2, have been identified. Their origins remain controversial, particularly for Gaia BH1. By considering a rapidly rotating (ω/ω crit = 0.8) and strongly magnetized (B 0 = 5000 G) merger product, we find that, at typical Galactic metallicity, the merger product can undergo efficient chemically homogeneous evolution. This results in the merger product having a significantly smaller radius during its evolution compared to that of a normally evolving massive star. Under the condition that the initial triple stability is satisfied, we use the Multiple Stellar Evolution code and the MESA code to identify an initial hierarchical triple that can evolve into Gaia BH1. It initially consists of three stars with masses of 9.03 M ⊙, 3.12 M ⊙, and 1 M ⊙, with inner and outer orbital periods of 2.21 days and 121.92 days, and inner and outer eccentricities of 0.41 and 0.45, respectively. This triple initially experiences triple evolution dynamics instability (TEDI) followed by Roche lobe overflow (RLOF). During RLOF, the inner orbit shrinks, and tidal effects gradually suppress the TEDI. Eventually, the inner binary undergoes a merger through contact (or collision). Finally, using models of rapidly rotating and strongly magnetic stars, along with standard core-collapse supernova (SN) or failed supernova (FSN) models, we find that a postmerger binary (PMB) consisting of an 12.11 M ⊙ merger product and a 1 M ⊙ companion star (originally an outer tertiary) can avoid RLOF. After an SN or FSN with a low ejected mass of ∼0.22 M ⊙ and a low kick velocity ( 46−33+25kms−1 or 9−8+16kms−1 ), the PMB can form Gaia BH1 in the Galactic disk.

Impacts of Black-hole-forming Supernova Explosions on the Diffuse Neutrino Background

The flux spectrum, event rate, and experimental sensitivity are investigated for the diffuse supernova (SN) neutrino background (DSNB), which originates from past stellar collapses and is also known as a supernova relic neutrino background. For this purpose, the contribution of collapses that lead to successful supernova explosion and black hole (BH) formation simultaneously, which are suggested to be a nonnegligible population from the perspective of Galactic chemical evolution, is taken into account. If the BH-forming SNe involve matter fallback onto the protoneutron star for the long term, their total emitted neutrino energy becomes much larger than that of ordinary SNe and failed SNe (BH formation without explosion). Then, in the case of the normal mass hierarchy in neutrino oscillations and with half of all core-collapse SNe being BH-forming SNe, the expected event rate according to the current DSNB model is enhanced by up to a factor of 2 due to the BH-forming SNe. While substantial uncertainties exist regarding the duration of the matter fallback, which determines the total amount of emitted neutrinos, and the fraction of BH-forming SNe, the operation time required to detect the DSNB at Hyper-Kamiokande would be reduced by such contribution in any case.

Can Supernovae from Runaway Stars Mimic the Signs of Absorbing “Supervirial” Gas?

The recent detection of large column density absorption lines from highly ionized gas in a few directions through the circumgalactic medium (CGM) of the Milky Way (MW) has been puzzling. The inferred temperature from these absorption lines far exceeds the virial temperature of the MW, and the column densities are also too large to be easily explained. In this paper, we propose a novel idea to explain these observations and claim that they may not have originated from the CGM, but from a totally different type of source, namely, stellar ejecta from supernovae (SNe) above the Galactic disk that happen to lie in the line of sight to the background quasars. About ∼20% of massive OB stars (progenitors of core-collapse supernovae) are known to be runaway stars that have high ejection velocities near the Galactic plane and can end up exploding as SNe above the Galactic disk. We show that the associated reverse shock in the supernova remnant in the early nonradiative phase can heat the ejecta to temperatures of ≳107 K and can naturally explain the observed high column density of ions in the observed “supervirial” phase along with α-enriched supersolar abundance that is typical of core-collapse supernovae. However, SNe from runaway stars has a covering fraction of ≲0.7% and thus can only explain the observations along limited sightlines.

Gravitational-wave and Gravitational-wave Memory Signatures of Core-collapse Supernovae

In this paper, we calculate the energy, signal-to-noise ratio (SNR), detection range, and angular anisotropy of the matter, matter memory, and neutrino memory gravitational-wave (GW) signatures of 21 three-dimensional initially nonrotating core-collapse supernova (CCSN) models carried to late times. We find that inferred energy, SNR, and detection range are angle-dependent quantities, and that the spread of possible energy, signal to noise, and detection ranges across all viewing angles generally increases with progenitor mass. When examining the low-frequency matter memory and neutrino memory components of the signal, we find that the neutrino memory is the most detectable component of a CCSN GW signal, and that DECIGO is best equipped to detect both matter memory and neutrino memory. Moreover, we find that the polarization angle between the h + and h × strains serves as a unique identifier of matter and neutrino memory. Finally, we develop a Galactic density- and stellar mass-weighted formalism to calculate the rate at which we can expect to detect CCSN GW signals with the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). When considering only the matter component of the signal, the aLIGO detection rate is around 65% of the total Galactic supernova rate, but increases to 90% when incorporating the neutrino memory component. We find that all future detectors (Einstein Telescope, Cosmic Explorer, DECIGO) will be able to detect CCSN GW signals from the entire Galaxy, and for the higher-mass progenitors even into the Local Group of galaxies.

Origin of LAMOST J1010+2358 Revisited

Signature from Pop III massive stars of 140– that end their lives as pair-instability supernovae (PISNe) are expected to be seen in very metal-poor (VMP) stars of . Although thousands of VMP stars have been discovered, the identification of a VMP star with a PISN signature has been elusive. Recently, the VMP star LAMOST J1010+2358 was claimed to be the first star with a clear PISN signature. A subsequent study showed that ejecta from low-mass core-collapse supernovae (CCSNe) can also fit the abundance pattern equally well and additional elements such as C and Al are required to differentiate the two sources. Follow-up observations of LAMOST J1010+2358 by two independent groups were able to detect both C and Al. Additionally, key odd elements such as Na and Sc were also detected whose abundances were found to be higher than the upper limits found in the original detection. We perform a detailed analysis of the newly observed abundance patterns by exploring various possible formation channels for VMP stars. We find that purely low-mass CCSN ejecta as well as the combination of CCSN and Type 1a SN ejecta can provide an excellent fit to the newly observed abundance pattern. Our results confirm earlier analysis that the newly observed abundance pattern is peculiar but has no signatures of PISN.

Measurement of the absolute efficiency of the X-ARAPUCA photon detector for the DUNE Far Detector 1

The DUNE far detector has been designed to detect photons and electrons generated by the charged products of the interaction of neutrinos with a massive liquid argon (LAr) target. The photon detection system (PDS) of the first DUNE far detector (FD1) is composed of 6000 photon detection units, named X-ARAPUCA. The detection of the prompt light pulse generated by the particle energy release in LAr will complement and boost the DUNE LAr Time Projection Chamber. It will improve the non-beam events tagging and enable at low energies the trigger and the calorimetry of the supernova neutrinos. The X-ARAPUCA is an assembly of several components. Its photon detection efficiency (PDE) depends on the design of the assembly, on the grade of the individual components and on their coupling. The X-ARAPUCA PDE is one of the leading parameters for the PDS sensitivity, that in turn determines the sensitivity of the DUNE for the detection of core-collapse supernova within the galaxy and for nucleon decay searches. In this work we present the final assessment of the absolute PDE of the FD1 X-ARAPUCA baseline design, measured in two laboratories with independent methods and setups. Preliminary results were reported in Palomares (JINST 18(02):C02064, https://doi.org/10.1088/1748-0221/18/02/C02064, 2023). One hundred sixty units of these X-ARAPUCA devices have been deployed in the NP04 facility at the CERN Neutrino Platform, the 1:20 scale FD1 prototype, and will be operated during the year 2024. The assessed value of the PDE is a key parameter both in the NP04 and in the DUNE analysis and reconstruction studies.

Late-time supernovae radio re-brightening in the VAST pilot survey

ABSTRACT We present our analysis of supernovae serendipitously found to be radio-bright several years after their optical discovery. We used recent observations from the Australian SKA Pathfinder (ASKAP) taken as part of the pilot Variables and Slow Transients and Rapid ASKAP Continuum Survey programmes. We identified 29 objects by cross-matching sources from these ASKAP observations with known core-collapse supernovae below a declination of $+40^{\circ }$ and with a redshift of $z\le 0.15$. Our results focus on eight cases that show potential late-time radio emission. These supernovae exhibit significantly greater amounts of radio emission than expected from the standard model of a single shockwave propagating through a spherical circumstellar medium, with a constant density structure produced by regular stellar mass-loss. We also discuss how we can learn from future ASKAP surveys about the circumstellar environments and emission mechanisms of supernovae that undergo late-time radio re-brightening. This pilot work tested and confirmed the potential of the Variables and Slow Transients survey to discover and study late-time supernova emission.

Quasithermal GeV neutrinos from neutron-loaded magnetized outflows in core-collapse supernovae: Spectra and light curves

Quasithermal GeV neutrinos from neutron-loaded magnetized outflows in core-collapse supernovae: Spectra and light curves

Insights into the Production of 44Ti and Nickel Isotopes in Core-collapse Supernovae

We report nucleosynthetic results for both 44Ti and nickel isotopes for 18 three-dimensional (3D) core-collapse supernova (CCSN) simulations extended to ∼20 s after bounce. We find that many of our long-term models are able to achieve 44Ti/56Ni ratios similar to that observed in Cassiopeia A, and modern supernova models can synthesize up to 2 × 10−4 M ⊙ of 44Ti. Neutrino-driven winds and the fact that there can be simultaneous accretion and explosion in 3D models of CCSNe play central roles in its production. We conclude that the 44Ti underproduction problem in previous CCSN models is no longer an issue. In addition, we discuss the production of both 57Ni and stable nickel/iron ratios and compare our results to observations of SN 1987A and the Crab.