Core-collapse Supernovae Research Articles

An Extensive Hubble Space Telescope Study of the Offset and Host Light Distributions of Type I Superluminous Supernovae

We present an extensive Hubble Space Telescope rest-frame UV imaging study of the locations of Type I superluminous supernovae (SLSNe) within their host galaxies. The sample includes 65 SLSNe with detected host galaxies in the redshift range z ≈ 0.05–2. Using precise astrometric matching with SN images, we determine the distributions of the physical and host-normalized offsets relative to the host centers, as well as the fractional flux distribution relative to the underlying UV light distributions. We find that the host-normalized offsets of SLSNe roughly track an exponential disk profile, but exhibit an overabundance of sources with large offsets of 1.5–4 times their hosts' half-light radii. The SLSNe normalized offsets are systematically larger than those of long gamma-ray bursts (LGRBs), and even Type Ib/c and Type II SNe. Furthermore, we find from a Monte Carlo procedure that about 37−8+6% of SLSNe occur in the dimmest regions of their host galaxies, with a median fractional flux value of 0.16, in stark contrast to LGRBs and Type Ib/c and Type II SNe. We do not detect any significant trends in the locations of SLSNe as a function of redshift, or as a function of explosion and magnetar engine parameters inferred from modeling of their optical light curves. The significant difference in SLSN locations compared to LGRBs (and normal core-collapse SNe) suggests that at least some of their progenitors follow a different evolutionary path. We speculate that SLSNe arise from massive runaway stars from disrupted binary systems, with velocities of ∼102 km s−1.

Low-energy Explosions in a Gravitational Field: Implications for Sub-energetic Supernovae and Fast X-Ray Transients

Observations and theory suggest that core-collapse supernovae can span a range of explosion energies, and when sub-energetic the shockwave initiating the explosion can decelerate to speeds comparable to the escape speed of the progenitor. In these cases, gravity will complicate the explosion hydrodynamics and conceivably cause the shock to stall at large radii within the progenitor star. To understand these unique properties of weak explosions, we develop a perturbative approach for modeling the propagation of an initially strong shock into a time-steady, infalling medium in the gravitational field of a compact object. This method writes the shock position and the post-shock velocity, density, and pressure as series solutions in the (time-dependent) ratio of the freefall speed to the shock speed, and predicts that the shock stalls within the progenitor if the explosion energy is below a critical value. We show that our model agrees very well with hydrodynamic simulations, and accurately predicts (for example) the time-dependent shock position and velocity and the radius at which the shock stalls. Our results have implications for black hole formation and the newly detected class of fast X-ray transients (FXTs). In particular, we propose that a “phantom shock breakout”—where the outer edge of the star falls through a stalled shock—can yield a burst of X-rays without a subsequent optical/UV signature, similar to FXTs. This model predicts the rise time of the X-ray burst, t d, and the mean photon energy, kT, are anticorrelated, approximately as T∝td−5/8 .

The origin of the X-ray emission from the non-starburst gas-rich luminous infrared galaxies Arp 302

ABSTRACT We present an analysis of the XMM–Newton observation of luminous infrared merging galaxies Arp 302 and a joint re-analysis of its Chandra observation. In particular, we focus on the more significant X-ray emitter of the pair, Arp 302N. Chandra detects significant soft X-ray emission from the hot gas in the star-forming region of Arp 302N spreading up to 12 kpc. We estimate the star formation rate of Arp 302N to be around 1–2 M⊙ yr−1 based on the X-ray luminosity of the star-forming region, similar to previous measurements at longer wavelengths. Chandra and XMM–Newton observations show evidence of a Si xiii emission line with 86 per cent confidence. Our best-fitting model infers a super-solar silicon abundance in the star-forming region, likely related to the past core-collapse supernovae in this galaxy. Similar silicon overabundance was reported in the circumstellar medium of core-collapse supernova remnants in our Galaxy. We also detect narrow Fe Kα and Fe Kβ (98.6 per cent confidence) emission lines as part of the active galactic nucleus (AGN) emission. Our best-fitting spectral model using mytorus indicates the evidence of a heavily obscured power-law emission with NH > 3 × 1024 cm−2 in addition to a weak, unobscured power-law emission. The scattering fraction of the unobscured power-law emission from Compton-thin materials is 0.7 per cent. All these spectral features suggest evidence of a Seyfert 2-like AGN in Arp 302N. The X-ray measurement of its AGN activity is consistent with the previous Spitzer measurement of the same object.

Constraining MeV-scale axionlike particles with Fermi-LAT observations of SN 2023ixf

The Fermi-LAT observations of SN 2023ixf, a type II supernova in the nearby Pinwheel Galaxy, Messier 101 (M101), presents us with an excellent opportunity to constrain MeV-scale axionlike particles (ALPs). By examining the photon decay signature from heavy ALPs that could be produced in the explosion, the existing constraints on the ALP-photon coupling can be improved, under optimistic assumptions, by up to a factor of ∼2 for masses ma≲3 MeV. Under very conservative assumptions, we find a bound that is slightly weaker than the existing ones for ma≲0.5 MeV. The exact reach of these searches depends mostly on properties of the SN progenitor. This study demonstrates the relevance of core-collapse supernovae, also beyond the Magellanic Clouds, as probes of fundamental physics. Published by the American Physical Society 2024

Collisional and fast neutrino flavor instabilities in two-dimensional core-collapse supernova simulation with Boltzmann neutrino transport

Collisional and fast neutrino flavor instabilities in two-dimensional core-collapse supernova simulation with Boltzmann neutrino transport

On the treatment of phenomenological turbulent effects in one-dimensional simulations of core-collapse supernovae

ABSTRACT We have developed a phenomenological turbulent model with one-dimensional (1D) simulation based on Reynolds decomposition. Using this method, we have systematically studied models with different effects of compression, mixing length parameters, and diffusion coefficient of internal energy, turbulence energy, and electron fraction. With employed turbulent effects, supernova explosion can be achieved in 1D geometry, which can mimic the evolution of shock in the 3D simulations. We found that enhancement of turbulent energy by compression affects the early shock evolution. The diffusion coefficients of internal energy and turbulent energy also affect the explodability. The smaller diffusion makes the shock revival faster. Our comparison between the two reveals that the diffusion coefficients of internal energy has a greater impact. These simulations would help understand the role of turbulence in core-collapse supernovae.

Detectability of neutrino-signal fluctuations induced by the hadron-quark phase transition in failing core-collapse supernovae

Detectability of neutrino-signal fluctuations induced by the hadron-quark phase transition in failing core-collapse supernovae

Real-time monitoring for the next core-collapse supernova in JUNO

The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neutrinos at the Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton liquid scintillator detector currently under construction in South China.The real-time monitoring system is designed to ensure both prompt alert speed and comprehensive coverage of progenitor stars. It incorporates prompt monitors on the electronic board as well as online monitors at the data acquisition stage.Assuming a false alert rate of 1 per year, this monitoring system exhibits sensitivity to pre-SN neutrinos up to a distance of approximately 1.6 (0.9) kiloparsecs and SN neutrinos up to about 370 (360) kiloparsecs for a progenitor mass of 30 solar masses, considering both normal and inverted mass ordering scenarios.The pointing ability of the CCSN is evaluated by analyzing the accumulated event anisotropy of inverse beta decay interactions from pre-SN or SN neutrinos. This, along with the early alert, can play a crucial role in facilitating follow-up multi-messenger observations of the next galactic or nearby extragalactic CCSN.

Keck Infrared Transient Survey. I. Survey Description and Data Release 1

We present the Keck Infrared Transient Survey, a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roman SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. Here, we describe the first release of data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients <17 mag in a red optical band (usually ZTF r or ATLAS o bands); a volume-limited sample including all transients within redshift z < 0.01 (D ≈ 50 Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete to z = 0.005. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events, luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction using PypeIt, which requires minimal human interaction to ensure reproducibility.

The Enge Split-Pole Spectrograph at the University of Notre Dame

Nuclear reactions play a crucial role in determining the nucleosynthesis that occurs in astrophysical events. The rates of many reactions that significantly impact certain nucleosynthesis processes can not be currently measured via direct means. These reactions must be constrained in another manner, such as determining the level energies and other structure properties of the compound nuclei. In order to measure level energies of nuclei relevant to nuclear astrophysics, the Enge split-pole spectrograph has been installed and commissioned at the University of Notre Dame’s Nuclear Science Laboratory. The first scientific measurement has also been performed. Structure properties of 58Cu were measured via the reaction 58Ni(3He,t)58Cu to provide the first experimental constraint of the 57Ni(p,γ)58Cu reaction rate, which impacts the production of of 44Ti, 57Fe, and 59Ni in core-collapse supernovae. Preliminary analysis of this measurement confirms the level energies of states in 58Cu that could lead to significant resonances in the 57Ni(p,γ)58Cu reaction rate, while suggesting the presence of additional states that have not been previously observed but could also lead to significant resonances.

Evolution and Final Fates of a Rotating 25M⊙ Pop III Star

In this manuscript, we present the one-dimensional stellar evolution of two rotating population III (Pop III) star models, each having a mass of 25M⊙ at the zero-age main-sequence (ZAMS). The slowly rotating model has an initial angular rotational velocity of 10% of the critical angular rotational velocity. In contrast, the rapidly rotating model has an initial angular rotational velocity of 70% of the critical angular rotational velocity. As an effect of rotationally enhanced mixing, we find that the rapidly rotating model suffers an enormous mass loss due to the deposition of a significant amount of CNO elements toward the surface after the main-sequence phase. We also display the simulated light curves as these models explode into core-collapse supernovae (CCSNe).

Search for Galactic Core-collapse Supernovae in a Decade of Data Taken with the IceCube Neutrino Observatory

The IceCube Neutrino Observatory has been continuously taking data to search for O(0.5–10) s long neutrino bursts since 2007. Even if a Galactic core-collapse supernova is optically obscured or collapses to a black hole instead of exploding, it will be detectable via the O(10) MeV neutrino burst emitted during the collapse. We discuss a search for such events covering the time between 2008 April 17 and 2019 December 31. Considering the average data taking and analysis uptime of 91.7% after all selection cuts, this is equivalent to 10.735 yr of continuous data taking. In order to test the most conservative neutrino production scenario, the selection cuts were optimized for a model based on an 8.8 solar mass progenitor collapsing to an O–Ne–Mg core. Conservative assumptions on the effects of neutrino oscillations in the exploding star were made. The final selection cut was set to ensure that the probability to detect such a supernova within the Milky Way exceeds 99%. No such neutrino burst was found in the data after performing a blind analysis. Hence, a 90% C.L. upper limit on the rate of core-collapse supernovae out to distances of ≈25 kpc was determined to be 0.23 yr−1. For the more distant Magellanic Clouds, only high neutrino luminosity supernovae will be detectable by IceCube, unless external information on the burst time is available. We determined a model-independent limit by parameterizing the dependence on the neutrino luminosity and the energy spectrum.

Spectropolarimetry of Type II supernovae

The explosion processes of supernovae (SNe) are imprinted in their explosion geometries. The recent discovery of several highly aspherical core-collapse SNe is significant, and studying these is regarded as being crucial in order to understand the underlying explosion mechanism. Here, we study the intrinsic polarization of 15 hydrogen-rich core-collapse SNe and explore the relation between polarization and the photometric and spectroscopic properties of these objects. Our sample shows diverse properties of the continuum polarization. Most SNe show a low degree of polarization at early phases but a sudden rise to ∼1% at certain points during the photospheric phase followed by a slow decline during the tail phase, with a constant polarization angle. The variation in the timing of peak polarization values implies diversity in the explosion geometry: some SNe have aspherical structures only in their helium cores, while in other SNe such structures reach out to a significant part of the outer hydrogen envelope with a common axis from the helium core to the hydrogen envelope. Other SNe show high polarization from early phases and a change in polarization angle around the middle of the photospheric phase. This implies that the ejecta are significantly aspherical out to the outermost layer and have multi-directional aspherical structures. Exceptionally, the Type IIL SN 2017ahn shows low polarization at both the photospheric and tail phases. Our results show that the timing of the polarization rise in Type IIP SNe is likely correlated with their brightness, velocity, and the amount of radioactive Ni produced: brighter SNe with faster ejecta velocity and a larger 56Ni mass have more extended aspherical explosion geometries. In particular, there is a clear correlation between the timing of the polarization rise and the explosion energy; that is, the explosion asphericity is proportional to the explosion energy. This implies that the development of a global aspherical structure, such as a jet, might be the key for the realisation of an energetic SN in the mechanism of SN explosions.

Supernova model discrimination with a kilotonne-scale Gd-H2O Cherenkov detector

Abstract The supernova model discrimination capabilities of the WATCHMAN detector concept are explored. This cylindrical kilotonne-scale water Cherenkov detector design has been developed to detect reactor antineutrinos through inverse β-decay for non-proliferation applications but also has the ability to observe antineutrino bursts of core-collapse supernovae within our galaxy. Detector configurations with sizes ranging from 16 m to 22 m tank diameter and 10% to 20% PMT coverage are used to compare the expected observable antineutrino spectra based on the Nakazato, Vartanyan and Warren supernova models. These spectra are then compared to each other with a fixed event count of 100 observed inverse β-decay events and a benchmark supernova at 10 kpc distance from Earth. By comparing the expected spectra, each detector configuration's ability to distinguish is evaluated. This analysis then demonstrates that the detector design is capable of meaningful event discrimination (90+% accuracy) with 100 observed supernova antineutrino events in most configurations. Furthermore, a larger tank configuration can maintain this performance at 10 kpc distance and above, indicating that overall target mass is the main factor for such a detector's discrimination capabilities. Finally, it is estimated that the detector design can provide early warning capability for supernova bursts for the entire Milky Way in all configurations.

Detectability of Late-time Supernova Neutrinos with Fallback Accretion onto Protoneutron Star

We investigate the late-time neutrino emission powered by fallback mass accretion onto a protoneutron star (PNS), using neutrino radiation-hydrodynamic simulations with full Boltzmann neutrino transport. We follow the time evolution of the accretion flow onto the PNS until the system reaches a quasi-steady state. A standing shock wave is commonly formed in the accretion flow, whereas the shock radius varies depending on the mass accretion rate and the PNS mass. A sharp increase in temperature emerges in the vicinity of the PNS (∼10 km), which characterizes neutrino emission. Both the neutrino luminosity and the average energy become higher with increasing mass accretion rate and PNS mass. The mean energy of the emitted neutrinos is in the range of 10 ≲ ϵ ≲ 20 MeV, which is higher than that estimated from PNS cooling models (≲10 MeV). Assuming a distance to core-collapse supernova of 10 kpc, we quantify neutrino event rates for Super-Kamiokande (Super-K) and Deep Underground Neutrino Experiment (DUNE). The estimated detection rates are well above the background, and their energy-dependent features are qualitatively different from those expected from PNS cooling models. Another notable feature is that the neutrino emission is strongly flavor dependent, exhibiting that the neutrino event rate hinges on the neutrino oscillation model. We estimate them in the case with the adiabatic Mikheev–Smirnov–Wolfenstein model, and show that the normal and inverted mass hierarchy offer a large number of neutrino detections in Super-K and DUNE, respectively. Hence the simultaneous observation with Super-K and DUNE of fallback neutrinos will provide a strong constraint on the neutrino mass hierarchy.

A New Insight into Electron Acceleration Properties from Theoretical Modeling of Double-peaked Radio Light Curves in Core-collapse Supernovae

It is recognized that some core-collapse supernovae (SNe) show a double-peaked radio light curve within a few years since the explosion. A shell of circumstellar medium (CSM) detached from the SN progenitor has been considered to play a viable role in characterizing such a rebrightening of radio emission. Here, we propose another mechanism that can give rise to the double-peaked radio light curve in core-collapse SNe. The key ingredient in the present work is to expand the model for the evolution of the synchrotron spectral energy distribution (SED) to a generic form, including fast and slow cooling regimes, as guided by the widely accepted modeling scheme of gamma-ray burst afterglows. We show that even without introducing an additional CSM shell, the radio light curve would show a double-peaked morphology when the system becomes optically thin to synchrotron self-absorption at the observational frequency during the fast cooling regime. We can observe this double-peaked feature if the transition from the fast cooling to slow cooling regime occurs during the typical observational timescale of SNe. This situation is realized when the minimum Lorentz factor of injected electrons is initially large enough for the nonthermal electrons’ SED to be discrete from the thermal distribution. We propose SN 2007bg as a special case of double-peaked radio SNe that can be possibly explained by the presented scenario. Our model can serve as a potential diagnostic for electron acceleration properties in SNe.

Search for Supernova Progenitor Stars with ZTF and LSST

The direct detection of core-collapse supernova (SN) progenitor stars is a powerful way of probing the last stages of stellar evolution. However, detections in archival Hubble Space Telescope images are limited to about one detection per year. Here, we explore whether we can increase the detection rate by using data from ground-based wide-field surveys. Due to crowding and atmospheric blurring, progenitor stars can typically not be identified in preexplosion images alone. Instead, we combine many pre-SN and late-time images to search for the disappearance of the progenitor star. As a proof of concept, we implement our search of ZTF data. For a few hundred images, we achieve limiting magnitudes of ∼23 mag in the g and r bands. However, no progenitor stars or long-lived outbursts are detected for 29 SNe within z ≤ 0.01, and the ZTF limits are typically several magnitudes less constraining than detected progenitors in the literature. Next, we estimate progenitor detection rates for the Legacy Survey of Space and Time (LSST) with the Vera C. Rubin telescope by simulating a population of nearby SNe. The background from bright host galaxies reduces the nominal LSST sensitivity by, on average, 0.4 mag. Over the 10 yr survey, we expect the detection of ∼50 red supergiant progenitors and several yellow and blue supergiants. The progenitors of Type Ib and Ic SNe will be detectable if they are brighter than −4.7 or −4.0 mag in the LSST i band, respectively. In addition, we expect the detection of hundreds of pre-SN outbursts depending on their brightness and duration.

Nonlinear Evolution of Fast Collective Neutrino Oscillations in Core-Collapse Supernovae

According to one of the most promising supernova theories, the neutrino-heating mechanism, neutrinos are responsible for transferring the energy released during the gravitational collapse of massive stars to their surroundings. If neutrino flavors are converted fast in the cores, the efficiency of neutrino heating is enhanced and can change the dynamics of the shock wave in supernovae. In this article, we investigate the dynamics of fast neutrino flavor conversions with collisions under energy-dependent treatment in detail. For the first time, we use a realistic initial condition, which is taken from the results of the self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry. We report that the neutrino flavor conversion will be significantly enhanced if the energy-dependent collision term is considered in the neutrino transport equation which has been ignored in previous studies. We present the preliminary results of our investigation here, which are consistent with other studies which were conducted under simple treatments. It is believed that such findings may have an impact on the explosion mechanism of core-collapse supernovae, nucleosynthesis and neutrino astronomy.

Visualizing convolutional neural network for classifying gravitational waves from core-collapse supernovae

Visualizing convolutional neural network for classifying gravitational waves from core-collapse supernovae

Universality of the neutrino collisional flavor instability in core-collapse supernovae

Universality of the neutrino collisional flavor instability in core-collapse supernovae