Supernova Rate Research Articles

BPASS stellar evolution models incorporating α-enhanced composition – I. Single star models from 0.1 to 316 M⊙

Abstract Stellar evolution modelling is fundamental to many areas of astrophysics including stellar populations in both nearby and distant galaxies. It is heavily influenced by chemical composition. Observations of distant galaxies and nucleosynthesis calculations show that α-process elements are enriched faster than iron group elements. We present a dense grid of single-star models calculated using the bpass stellar evolution code and covering masses (0.1 ≤ M/M⊙≤316), metallicity mass fractions (10−5 ≤ Z ≤ 0.04) and α-to-iron abundance ratios (−0.2 ≤$[\alpha /\rm {Fe}]$≤+0.6). By comparing Solar-scaled models to ones enriched in α-process elements, we find that stellar radii, surface temperatures, Main Sequence lifetimes, supernova progenitor properties and supernova rates are all sensitive to changes in $[\alpha /\rm {Fe}]$. Lifetimes of low-mass stars differ by up to 0.4 dex, while surface temperatures of massive stars at the end of the Main Sequence also differ by around 0.4 dex. Inferred supernova rates when $[\rm {Fe}/\rm {H}]$ is unknown can be highly uncertain. Models with different $[\alpha /\rm {Fe}]$ but comparable iron abundances show smaller variations, indicating that while iron primarily defines the course of evolution; α-enhancement nonetheless has an impact of up to 0.1 dex on stellar properties. Such changes are small for individual stars, but have a large cumulative effect when considering an entire stellar population as demonstrated by isochrone fitting to nearby clusters. Changes in radii and lifetimes have further consequences for a stellar population including binary stars, as they influence the timing, nature and occurrence rate of mass transfer events.

Nebular and Nonthermal Radio Emissions for Young Stellar Populations with PARSEC v1.2s

Abstract In this paper, we compute, by means of the recently and thoroughly updated PARSE v1.2 s database of stellar nonrotating evolutionary tracks, the integrated stellar spectra, the ionizing photon budget, and the supernovae rates of young simple stellar populations (SSPs), for five metallicities between 0.0001 and 0.02 and four choices of stellar initial mass function (IMF) upper mass limits between 40 M ⊙ and 350 M ⊙. Using the photo-ionization code CLOUDY, we compute, at this same range of metallicities and limits, the intensities of some selected recombination and collisionally excited lines as a function of the age of the SSP. We account for the electron temperature dependence on IMF upper mass limit and metallicity while computing the thermal radio emission component, and also accounted for recent advances in core-collapse supernova explosion models while computing the nonthermal radio emission component. We self-consistently add the emission lines, nebular continuum, and nonthermal radio emission to the original SSP integrated photospheric spectra. Finally, from the resulting new suite of SSPs, we provide a consistent set of analytical relations between star formation rate (SFR) and ultraviolet, optical, and thermal radio luminosities that can be used to convert attenuation-corrected and dust-unaffected luminosities to SFR estimates. In a forthcoming paper, we will use our new SSP libraries as input to the state-of-the-art radiative transfer model GRAphites and SILicates to test the overall performance of these SSPs in reproducing the observed spectral energy distribution of young star-forming galaxies.

The Impact of Host-galaxy Properties on Supernova Classification with Hierarchical Labels

With the advent of the Vera C. Rubin Observatory, the discovery rate of supernovae (SNe) will surpass the rate of SNe with real time spectroscopic follow-up by 3 orders of magnitude. Accurate photometric classifiers are essential to both select interesting events for follow-up in real time and for archival population-level studies. In this work, we investigate the impact of observable host-galaxy information on the classification of SNe, both with and without additional light-curve and redshift information. We find that host-galaxy information alone can successfully isolate relatively pure (>90%) samples of Type Ia SNe with or without redshift information. With redshift information, we can additionally produce somewhat pure (>70%) samples of Type II SNe and superluminous SNe. Additionally with redshift information, host-galaxy properties do not significantly improve the accuracy of SN classification when paired with complete light curves. In the absence of redshift information, however, galaxy properties significantly increase the accuracy of photometric classification. As a part of this analysis, we present the first formal application of a new objective function, the weighted hierarchical cross entropy, to the problem of SN classification. This objective function more naturally accounts for the hierarchical nature of SN classes and, more broadly, transients. Finally, we present a new set of SN classifications for the Pan-STARRS Medium Deep Survey of SNe that lack spectroscopic redshift, increasing the full photometric sample to >4400 events.

SN 1885A and Supernova Remnants in the Centre of M31 with LOFAR

We present the first LOFAR image of the center of M31 at a frequency of 150 MHz. We clearly detect three supernova remnants, which, along with archival VLA data at 3 GHz and other published radio and X-ray data, allows us to characterize them in detail. Our observations also allow us to obtain upper limits of the historical SN 1885A, which is undetected even at a low frequency of 150 MHz. From analytical modeling, we find that SN 1885A will stay in its free-expansion phase for at least another couple of centuries. We find an upper limit of n H ≲ 0.04 cm−3 for the interstellar medium of SN 1885A, and that the SN ejecta density is not shallower than ∝r −9 (on average). From the 2.6σ tentative detection in X-ray, our analysis shows that nonthermal emission is expected to dominate the SN 1885A emission. Comparing our results with those on G1.9+0.3, we find that it is likely that the asymmetries in G1.9+0.3 make it a more efficient radio and X-ray emitter than SN 1885A. For Braun 80, 95, and 101, the other remnants in this region, we estimate ages of 5200, 8100, and 13,100 yr, and shock speeds of 1150, 880, and 660 km s−1, respectively. Based on this, the supernova rate in the central 0.5 kpc × 0.6 kpc of M31 is at least one per ∼3000 yr. We estimate radio spectral indices of −0.66 ± 0.05, −0.37 ± 0.03, and −0.50 ± 0.03 for the remnants, respectively, which match fairly well with previous studies.

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.

Using dust to constrain dark matter models

Abstract In this paper, we use hydrodynamic zoom-in simulations of Milky Way-type haloes to explore using dust as an observational tracer to discriminate between cold and warm dark matter universes. Comparing a cold and 3.5keV warm dark matter particle model, we tune the efficiency of galaxy formation in our simulations using a variable supernova rate to create Milky Way systems with similar satellite galaxy populations while keeping all other simulation parameters the same. Cold dark matter, having more substructure, requires a higher supernova efficiency than warm dark matter to achieve the same satellite galaxy number. These different supernova efficiencies create different dust distributions around their host galaxies, which we generate by post-processing the simulation output with the powderday codebase. Analysing the resulting dust in each simulation, we find ∼4.5 times more dust in our cold dark matter Milky Way halos compared with warm dark matter. The distribution of dust out to R200c is then explored, revealing that the warm dark matter simulations are noticeably less concentrated than their cold dark matter counterparts, although differences in substructure complicate the comparison. Our results indicate that dust is a possible unique probe to test theories of dark matter.

Constraining the Initial Mass Function via Stellar Transients

The stellar initial mass function (IMF) represents a fundamental quantity in astrophysics and cosmology describing the mass distribution of stars from low mass all the way up to massive and very massive stars. It is intimately linked to a wide variety of topics, including stellar and binary evolution, galaxy evolution, chemical enrichment, and cosmological reionization. Nonetheless, the IMF still remains highly uncertain. In this work, we aim to determine the IMF with a novel approach based on the observed rates of transients of stellar origin. We parametrize the IMF with a simple but flexible Larson shape, and insert it into a parametric model for the cosmic UV luminosity density, local stellar mass density, type Ia supernova (SN Ia), core-collapse supernova (CCSN), and long gamma-ray burst (LGRB) rates as a function of redshift. We constrain our free parameters by matching the model predictions to a set of empirical determinations for the corresponding quantities via a Bayesian Markov Chain Monte Carlo method. Remarkably, we are able to provide an independent IMF determination with a characteristic mass mc=0.10−0.08+0.24M⊙ and high-mass slope ξ=−2.53−0.27+0.24 that are in accordance with the widely used IMF parameterizations (e.g., Salpeter, Kroupa, Chabrier). Moreover, the adoption of an up-to-date recipe for the cosmic metallicity evolution allows us to constrain the maximum metallicity of LGRB progenitors to Zmax=0.12−0.05+0.29Z⊙. We also find which progenitor fraction actually leads to SN Ia or LGRB emission (e.g., due to binary interaction or jet-launching conditions), put constraints on the CCSN and LGRB progenitor mass ranges, and test the IMF universality. These results show the potential of this kind of approach for studying the IMF, its putative evolution with the galactic environment and cosmic history, and the properties of SN Ia, CCSN, and LGRB progenitors, especially considering the wealth of data incoming in the future.

The Scale of Stellar Yields: Implications of the Measured Mean Iron Yield of Core Collapse Supernovae

The scale of α-element yields is difficult to predict from theory because of uncertainties in massive star evolution, supernova physics, and black hole formation, and it is difficult to constrain empirically because the impact of higher yields can be compensated by greater metal loss in galactic winds. We use a recent measurement of the mean iron yield of core collapse supernovae (CCSN) by Rodriguez et al., y¯Fecc=0.058±0.007M⊙ , to infer the scale of α-element yields by assuming that the plateau of [α/Fe] abundance ratios observed in low-metallicity stars represents the yield ratio of CCSN. For a plateau at [α/Fe]cc = 0.45, we find that the population-averaged yields of O and Mg are about equal to the solar abundance of these elements, logyOcc/ZO,⊙=logyMgcc/ZMg,⊙=−0.01±0.1 , where yXcc is the mass of element X produced by massive stars per unit mass of star formation. The inferred O and Fe yields agree with predictions of the Sukhbold et al. CCSN models assuming their Z9.6+N20 neutrino-driven engine, a scenario in which many progenitors with M < 40M ⊙ implode to black holes rather than exploding. The yields are lower than assumed in many models of the galaxy mass–metallicity relation, reducing the level of outflows needed to match observed abundances. Our one-zone chemical evolution models with η=Ṁout/Ṁ*≈0.6 evolve to solar metallicity at late times. By further requiring that models reach [α/Fe] ≈ 0 at late times, we infer a Hubble-time integrated Type Ia supernova rate of 1.1×10−3M⊙−1 , compatible with estimates from supernova surveys.

Long-term evolutionary links between the isolated neutron star populations

ABSTRACT We have investigated the evolutionary connections of the isolated neutron star (NS) populations including radio pulsars (RPs), anomalous X-ray pulsars (AXPs), soft gamma repeaters (SGRs), dim isolated NSs (XDINs), ‘high-magnetic field’ RPs (‘HBRPs’), central compact objects (CCOs), rotating radio transients (RRATs), and long-period pulsars (LPPs) in the fallback disc model. The model can reproduce these NS families as a natural outcome of different initial conditions (initial period, disc mass, and dipole moment, μ) with a continuous μ distribution in the $\sim 10^{27} - 5 \times 10^{30}$ G cm$^3$ range. Results of our simulations can be summarized as follows: (1) A fraction of ‘HBRPs’ with relatively high μ evolve into the persistent AXP/SGR properties, and subsequently become LPPs. (2) Persistent AXP/SGRs do not have evolutionary links with CCOs, XDINs, and RRATs. (3) For a wide range of μ, most RRATs evolve passing through RP or ‘HBRP’ properties during their early evolutionary phases. (4) A fraction of RRATs which have the highest estimated birth rate seem to be the progenitors of XDINs. (5) LPPs, whose existence was predicted by the fallback disc model, are the sources evolving in the late stage of evolution before the discs become inactive. These results provide concrete support to the ideas proposing evolutionary connections between the NS families to account for the ‘birth rate problem’, the discrepancy between the cumulative birth rate estimated for these systems and the core-collapse supernova rate.

A targeted search for strongly lensed supernovae with the Las Cumbres Observatory

ABSTRACT Gravitationally lensed supernovae (glSNe) are of interest for time delay cosmology and SN physics. However, glSNe detections are rare, owing to the intrinsic rarity of SN explosions, the necessity of alignment with a foreground lens, and the relatively short window of detectability. We present the Las Cumbres Observatory Lensed Supernova Search, LCOLSS, a targeted survey designed for detecting glSNe in known strong lensing systems. Using cadenced $r^\prime$-band imaging, LCOLSS targeted 114 galaxy–galaxy lensing systems with high expected SN rates, based on estimated star formation rates. No plausible glSN was detected by LCOLSS during our two year observing program. We carry out an analysis here to measure a detection efficiency for these observations. We then perform Monte Carlo simulations using the predicted supernova rates to determine the expected number of glSN detections. The results of the simulation suggest an expected number of detections and 68 per cent Poisson confidence intervals, $N_{\mathrm{ SN}} = 0.20, [0,2.1]$, $N_{\mathrm{ Ia}} = 0.08, [0,2.0]$, $N_{\mathrm{ CC}} = 0.12, [0,2.0]$, for all SNe, Type Ia SNe, and core-collapse (CC) SNe, respectively. These results are broadly consistent with the absence of a detection in our survey. Analysis of the survey strategy can provide insights for future efforts to develop targeted glSN discovery programs, especially considering the large anticipated yields of upcoming surveys.

The cosmic rate of pair-instability supernovae

Abstract Pair-instability supernovae (PISNe) have crucial implications for many astrophysical topics, including the search for very massive stars, the black hole mass spectrum, and galaxy chemical enrichment. To this end, we need to understand where PISNe are across cosmic time, and what are their favourable galactic environments. We present a new determination of the PISN rate as a function of redshift, obtained by combining up-to-date stellar evolution tracks from the PARSEC and FRANEC codes, with an up-to-date semi-empirical determination of the star formation rate and metallicity evolution of star-forming galaxies throughout cosmic history. We find the PISN rate to exhibit a huge dependence on the model assumptions, including the criterion to identify stars unstable to pair production, and the upper limit of the stellar initial mass function. Remarkably, the interplay between the maximum metallicity at which stars explode as PISNe, and the dispersion of the galaxy metallicity distribution, dominates the uncertainties, causing a ∼ seven-orders-of-magnitude PISN rate range. Furthermore, we show a comparison with the core-collapse supernova rate, and study the properties of the favourable PISN host galaxies. According to our results, the main contribution to the PISN rate comes from metallicities between ∼10−3 and 10−2, against the common assumption that views very-low-metallicity, Population III stars as exclusive or dominant PISN progenitors. The strong dependencies we find offer the opportunity to constrain stellar and galaxy evolution models based on possible future (or the lack of) PISN observations.

Linking Transients to their Host Galaxies: II. A Comparison of Host Galaxy Properties and Rate Dependencies across Supernova Types

Abstract We use the latest dataset of supernova (SN) host galaxies to investigate how the host properties – stellar mass, star formation rate, metallicity, absolute magnitude, and colour – differ across SN types, with redshift-driven selection effects controlled. SN Ib and Ic host galaxies, on average, are more massive, metal-rich, and redder than SN II hosts. For subtypes, SN Ibn and Ic-BL have bluer hosts than their normal SN Ib and Ic siblings; SN IIb has consistent host properties with SN Ib, while hosts of SN IIn are more metal-rich than those of SN II. Hydrogen-deficient superluminous supernovae feature bluer and lower luminosity hosts than most subtypes of core-collapse supernova (CC SN). Assuming simple proportionality of CC SN rates and host star formation rates (SFRs) does not recover the observed mean host properties; either a population of long-lived progenitors or a metallicity-dependent SN production efficiency better reproduces the observed host properties. Assuming the latter case, the rates of SN II are insensitive to host metallicity, but the rates of SN Ib and Ic are substantially enhanced in metal-rich hosts by a factor of ∼10 per dex increase in metallicity. Hosts of SN Ia are diverse in their observed properties; subtypes including SN Ia-91T, Ia-02cx, and Ia-CSM prefer star-forming hosts, while subtypes like SN Ia-91bg and Ca-rich prefer quiescent hosts. The rates of SN Ia-91T, Ia-02cx, and Ia-CSM are closely dependent on, or even proportional to, their host SFRs, indicating relatively short-lived progenitors. Conversely, the rates of SN Ia-91bg and Ca-rich transients are proportional to the total stellar mass, favoring long-lived progenitors.

Forecast of strongly lensed supernovae rates in the China Space Station Telescope surveys

Strong gravitationally lensed supernovae (SNe) are a powerful probe for cosmology and stellar physics. The relative time delays between lensed SN images provide an independent way of measuring a fundamental cosmological parameter --- the Hubble constant ---, the value of which is currently under debate. The time delays also serve as a ``time machine'', offering a unique opportunity to capture the extremely early phase of the SN explosion, which can be used to constrain the SN progenitor and explosion mechanism. Although there are only a handful of strongly lensed SN discoveries so far, which greatly hinders scientific applications, the sample size is expected to grow substantially with next-generation surveys. In this work, we investigate the capability of detecting strongly lensed SNe with the China Space Station Telescope (CSST), a two-meter space telescope to be launched around 2026. Through Monte Carlo simulations, we predict that CSST can detect 1008.53 and 51.78 strongly lensed SNe from its Wide Field Survey (WFS, covering 17,500 deg$^2$) and Deep Field Survey (DFS, covering 400 deg$^2$) over the course of ten years. In both surveys, about 35&lt;!PCT!&gt; of the events involve Type Ia SNe as the background sources. Our results suggest that the WFS and DFS of CSST, although not designed or optimized for discovering transients, can still make a great contribution to the strongly lensed SNe studies.

Cosmic Type Ia supernova rate and constraints on supernova Ia progenitors

Type Ia supernovae play a key role in the evolution of galaxies by polluting the interstellar medium with a fraction of iron peak elements larger than that released in the core-collapse supernova events. Their light curve, moreover, is widely used in cosmological studies as it constitutes a reliable distance indicator on extragalactic scales. Among the mechanisms proposed to explain the Type Ia supernovae (SNe), the single- and double-degenerate channels are thought to be the dominant ones, which implies a different distribution of time delays between the progenitor formation and the explosion. In this paper, we aim to determine the dominant mechanism by comparing a compilation of Type Ia SN rates with those computed from various cosmic star-formation histories coupled with different delay-time distribution functions. We also evaluate the relative contributions of both channels. By using a least-squares fitting procedure, we modeled the observations of Type Ia SN rates assuming different combinations of three recent cosmic star-formation rates and seven delay-time distributions. The goodness of these fits are statistically quantified by the $ test. For two of the three cosmic star-formation rates, the single degenerate scenario provides the most accurate explanation for the observations, while a combination of 34&lt;!PCT!&gt; single-degenerate- and 66&lt;!PCT!&gt; double-degenerate delay-time distributions is more plausible for the remaining tested cosmic star-formation rates. Though dependent on the assumed cosmic star-formation rate, we find arguments in favor of the single-degenerate model. From the theoretic point of view, at least sim 34&lt;!PCT!&gt; of the Type Ia SN must have been produced through the single-degenerate channel to account for the observations. The wide, double-degenerate mechanism slightly under-predicts the observations at redshift $z 1$, unless the cosmic SFR flattens in that regime. On the contrary, although the purely close double-degenerate scenario can be ruled out, we cannot rule out a mixed scenario with single- and double-degenerate progenitors.

Fate of supernova progenitors in massive binary systems

ABSTRACT How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate systematically how binary interactions influence core–collapse progenitors and their fates. Binary evolution simulations are performed to survey the parameter space of supernova progenitors in solar metallicity binary systems and to delineate major evolutionary paths. We first study fixed binary mass ratios ($q=M_2/M_1$ = 0.5, 0.7, and 0.9) to elucidate the impacts of initial mass and initial separation on the outcomes, treating separately Type Ibc supernova, Type II supernova, accretion-induced collapse (AIC), rapidly rotating supernova (Ibc-R), black hole formation, and long gamma ray burst (long GRB). We then conduct 12 binary population synthesis model calculations, varying the initial condition distributions and binary evolution parameters, to estimate various supernova fractions. We obtain a Milky Way supernova rate $R_{\rm SN} = (1.78$–$2.47) \times 10^{-2} \, {\rm yr}^{-1}$ which is consistent with observations. We find the rates of AIC, Ibc-R, and long GRB to be $\sim 1/100$ the rate of regular supernovae. Our estimated long GRB rates are higher than the observed long GRB rate and close to the low luminosity GRB rate, although care must be taken considering our models are computed with solar metallicity. Furthering binary modelling and improving the inputs one by one will enable more detailed studies of these and other transients associated with massive stars.

Luminosity Functions of the Host Galaxies of Supernova

We present the luminosity functions and stellar mass functions of supernova (SN) host galaxies and test if they differ from the functions of normal field galaxies. We utilize homogeneous samples consisting of 273 SNe Ia (z ≤ 0.3) and 44 core-collapse (CC) SNe (z ≤ 0.1) from the Sloan Digital Sky Survey II Supernova Survey and the high-signal-to-noise-ratio photometry of galaxies from the Hyper Suprime-Cam Subaru Strategic Program. SN hosts are classified into star-forming and passive galaxy groups based on the spectral energy distribution fitting. We find that the SN host luminosity functions and stellar mass functions deviate from those of normal field galaxies. Star-forming galaxies dominate the low-mass end of the SN Ia host mass function, while passive galaxies dominate the high-mass end. CC SNe are predominantly hosted by star-forming galaxies. In addition, intermediate-mass hosts produce CC SNe with the highest efficiency, while the efficiency of producing SNe Ia monotonically increases as the hosts become more massive. Furthermore, we derive the pseudo mass normalized SN rates (pSNuM) based on the mass functions. We find that the star-forming component of pSNuM Ia is less sensitive to the changes in stellar mass, in comparison with the total rate. The behavior of pSNuM CC suggests that the CC rate is proportional to the star-forming rate.

Modelling the galaxy radio continuum from star formation and active galactic nuclei in the Shark semi-analytic model

ABSTRACT We present a model of radio continuum emission associated with star formation (SF) and active galactic nuclei (AGNs) implemented in the Shark semi-analytic model of galaxy formation. SF emission includes free-free and synchrotron emission, which depend on the free-electron density and the rate of core-collapse supernovae with a minor contribution from supernova remnants, respectively. AGN emission is modelled based on the jet production rate, which depends on the black hole mass, accretion rate, and spin, and includes synchrotron self-absorption. Shark reproduces radio luminosity functions (RLFs) at $1.4\, \rm GHz$ and $150\, \rm MHz$ for 0 ≤ z ≤ 4, and scaling relations between radio luminosity, star formation rate, and infrared luminosity of galaxies in the local and distant universe in good agreement with observations. The model also reproduces observed number counts of radio sources from 150 MHz to 8.4 GHz to within a factor of 2 on average, though larger discrepancies are seen at the very bright fluxes at higher frequencies. We use this model to understand how the radio continuum emission from radio-quiet AGNs can affect the measured RLFs of galaxies. We find current methods to exclude AGNs from observational samples result in large fractions of radio-quiet AGNs contaminating the ‘star-forming galaxies’ selection and a brighter end to the resulting RLFs.We investigate how this affects the infrared-radio correlation (IRRC) and show that AGN contamination can lead to evolution of the IRRC with redshift. Without this contamination, our model predicts a redshift- and stellar mass-independent IRRC, except at the dwarf-galaxy regime.

Cosmic-Ray Acceleration of Galactic Outflows in Multiphase Gas

We investigate the dynamical interaction between cosmic rays (CRs) and the multiphase interstellar medium (ISM) using numerical magnetohydrodynamic (MHD) simulations with a two-moment CR solver and TIGRESS simulations of star-forming galactic disks. We previously studied the transport of CRs within TIGRESS outputs using a “postprocessing” approach, and we now assess the effects of the MHD backreaction to CR pressure. We confirm our previous conclusion that there are three quite different regimes of CR transport in multiphase ISM gas, while also finding that simulations with “live MHD” predict a smoother CR pressure distribution. The CR pressure near the midplane is comparable to other pressure components in the gas, but the scale height of CRs is far larger. Next, with a goal of understanding the role of CRs in driving galactic outflows, we conduct a set of controlled simulations of the extraplanar region above z = 500 pc, with imposed boundary conditions flowing from the midplane into this region. We explore a range of thermal and kinematic properties for the injected thermal gas, encompassing both hot, fast-moving outflows, and cooler, slower-moving outflows. The boundary conditions for CR energy density and flux are scaled from the supernova rate in the underlying TIGRESS model. Our simulations reveal that CRs efficiently accelerate extraplanar material if the latter is mostly warm/warm-hot gas, in which CRs stream at the Alfvén speed, and the effective sound speed increases as density decreases. In contrast, CRs have very little effect on fast, hot outflows where the Alfvén speed is small, even when the injected CR momentum flux exceeds the injected MHD momentum flux.

Supernova rates and luminosity functions from ASAS-SN I: 2014–2017 Type Ia SNe and their subtypes

ABSTRACT We present the volumetric rates and luminosity functions (LFs) of Type Ia supernovae (SNe Ia) from the V-band All-Sky Automated Survey for Supernovae (ASAS-SN) catalogues spanning discovery dates from UTC 2014 January 26 to UTC 2017 December 29. Our standard sample consists of 404 SNe Ia with $m_{\mathrm{{\it V},peak}} \lt 17\, \mathrm{mag}$ and Galactic latitude |b| &amp;gt; 15°. Our results are both statistically more precise and systematically more robust than previous studies due to the large sample size and high spectroscopic completeness. We make completeness corrections based on both the apparent and absolute magnitudes by simulating the detection of SNe Ia in ASAS-SN light curves. We find a total volumetric rate for all subtypes of $R_{\mathrm{tot}} = 2.28^{+0.20}_{-0.20} \times 10^{4}\, \mathrm{yr}^{-1}\, \mathrm{Gpc}^{-3}\, h^{3}_{70}$ for $M_{\mathrm{{\it V},peak}} \lt -16.5\, \mathrm{mag}$ ($R_{\mathrm{tot}} = 1.91^{+0.12}_{-0.12} \times 10^{4}\, \mathrm{yr}^{-1}\, \mathrm{Gpc}^{-3}\, h^{3}_{70}$ for $M_{\mathrm{{\it V},peak}} \lt -17.5\, \mathrm{mag}$) at the median redshift of our sample, zmed = 0.024. This is in agreement (1σ) with the local volumetric rates found by previous studies. We also compile LFs for the entire sample as well as for subtypes of SNe Ia for the first time. The major subtypes with more than one SN include Ia-91bg, Ia-91T, Ia-CSM, and Ia-03fg with total rates of $R_{\mathrm{Ia-91bg}} = 1.4^{+0.5}_{-0.5} \times 10^{3}\, \mathrm{yr}^{-1}\, \mathrm{Gpc}^{-3}\, h^{3}_{70}$, $R_{\mathrm{Ia-91T}} = 8.5^{+1.6}_{-1.7} \times 10^{2}\, \mathrm{yr}^{-1}\, \mathrm{Gpc}^{-3}\, h^{3}_{70}$, $R_{\mathrm{Ia-CSM}} = 10^{+7}_{-7}\, \mathrm{yr}^{-1}\, \mathrm{Gpc}^{-3}\, h^{3}_{70}$, and $R_{\mathrm{Ia-03fg}} = 30^{+20}_{-20}\, \mathrm{yr}^{-1}\, \mathrm{Gpc}^{-3}\, h^{3}_{70}$, respectively. We estimate a mean host extinction of $E(V-r) \approx 0.2\, \mathrm{mag}$ based on the shift between our V band and the Zwicky Transient Facility r-band LFs.

The JWST Discovery of the Triply Imaged Type Ia “Supernova H0pe” and Observations of the Galaxy Cluster PLCK G165.7+67.0

A Type Ia supernova (SN) at z = 1.78 was discovered in James Webb Space Telescope Near Infrared Camera imaging of the galaxy cluster PLCK G165.7+67.0 (G165; z = 0.35). The SN is situated 1.5–2 kpc from the host-galaxy nucleus and appears in three different locations as a result of gravitational lensing by G165. These data can yield a value for Hubble’s constant using time delays from this multiply imaged SN Ia that we call “SN H0pe.” Over the cluster, we identified 21 image multiplicities, confirmed five of them using the Near-Infrared Spectrograph, and constructed a new lens model that gives a total mass within 600 kpc of (2.6 ± 0.3) × 1014 M ⊙. The photometry uncovered a galaxy overdensity coincident with the SN host galaxy. NIRSpec confirmed six member galaxies, four of which surround the SN host galaxy with relative velocity ≲900 km s−1 and projected physical extent ≲33 kpc. This compact galaxy group is dominated by the SN host galaxy, which has a stellar mass of (5.0 ± 0.1) × 1011 M ⊙. The group members have specific star formation rates of 2–260 Gyr−1 derived from the Hα-line fluxes corrected for stellar absorption, dust extinction, and slit losses. Another group centered on a strongly lensed dusty star-forming galaxy is at z = 2.24. The total (unobscured and obscured) SFR of this second galaxy group is estimated to be (≳ 100 M ⊙ yr−1), which translates to a supernova rate of ∼1 SNe yr−1, suggesting that regular monitoring of this cluster may yield additional SNe.