Abstract We present initial results from a James Webb Space Telescope (JWST) survey of the youngest Galactic core-collapse supernova remnant, Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging ...mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI Medium Resolution Spectrograph integral field unit spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include (1) a weblike network of unshocked ejecta filaments resolved to ∼0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor’s oxygen layer with hot, high-entropy matter heated by neutrino interactions and radioactivity; (2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant’s interior pockmarked with small (∼1″) round holes formed by ≲0.″1 knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks; and (3) dozens of light echoes with angular sizes between ∼0.″1 and 1′ reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission (≲20 nJy at 3 μ m) from the neutron star in Cas A’s center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects.
Type IIb supernovae (SNe IIb) present a unique opportunity for investigating the evolutionary channels and mechanisms governing the evolution of stripped-envelope SN progenitors due to a variety of ...observational constraints. Comparison of these constraints with the full distribution of theoretical properties not only helps determine the prevalence of observed properties in nature, but can also reveal currently unobserved populations. In this follow-up paper, we use the large grid of models presented in Sravan et al. to derive distributions of single and binary SNe IIb progenitor properties and compare them to constraints from three independent observational probes: multiband SN light curves, direct progenitor detections, and X-ray/radio observations. Consistent with previous work, we find that while current observations exclude single stars as SN IIb progenitors, SN IIb progenitors in binaries can account for them. We also find that the distributions indicate the existence of an unobserved dominant population of binary SNe IIb at low metallicity that arise due to mass transfer initiated on the Hertzsprung Gap. In particular, our models indicate the existence of a group of highly stripped (envelope mass ∼0.1-0.2M☉) progenitors that are compact (<50R☉) and blue (Teff 105 K) with ∼104.5-105.5 L☉ and low-density circumstellar mediums. As discussed in Sravan et al., this group is necessary to account for SN IIb fractions and likely exist regardless of metallicity. The detection of the unobserved populations indicated by our models would support weak stellar winds and inefficient mass transfer in SN IIb progenitors.
We perform a Bayesian analysis of the mass distribution of stellar-mass black holes using the observed masses of 15 low-mass X-ray binary systems undergoing Roche lobe overflow and 5 high-mass, ...wind-fed X-ray binary systems. Using Markov Chain Monte Carlo calculations, we model the mass distribution both parametrically--as a power law, exponential, Gaussian, combination of two Gaussians, or log-normal distribution--and non-parametrically--as histograms with varying numbers of bins. We provide confidence bounds on the shape of the mass distribution in the context of each model and compare the models with each other by calculating their relative Bayesian evidence as supported by the measurements, taking into account the number of degrees of freedom of each model. The mass distribution of the low-mass systems is best fit by a power law, while the distribution of the combined sample is best fit by the exponential model. This difference indicates that the low-mass subsample is not consistent with being drawn from the distribution of the combined population. We examine the existence of a 'gap' between the most massive neutron stars and the least massive black holes by considering the value, M 1%, of the 1% quantile from each black hole mass distribution as the lower bound of black hole masses. Our analysis generates posterior distributions for M 1%; the best model (the power law) fitted to the low-mass systems has a distribution of lower bounds with M 1%>4.3 M with 90% confidence, while the best model (the exponential) fitted to all 20 systems has M 1%>4.5 M with 90% confidence. We conclude that our sample of black hole masses provides strong evidence of a gap between the maximum neutron star mass and the lower bound on black hole masses. Our results on the low-mass sample are in qualitative agreement with those of Ozel et al., although our broad model selection analysis more reliably reveals the best-fit quantitative description of the underlying mass distribution. The results on the combined sample of low- and high-mass systems are in qualitative agreement with Fryer & Kalogera, although the presence of a mass gap remains theoretically unexplained.
Abstract
This Letter presents the detection of a source at the position of the Type Ib/c supernova (SN) 2013ge more than four years after the radioactive component is expected to have faded. This ...source could mark the first post-SN direct detection of a surviving companion to a stripped-envelope Type Ib/c explosion. We test this hypothesis and find the shape of the source’s spectral energy distribution is most consistent with that of a B5 I supergiant. While binary models tend to predict OB-type stars for stripped-envelope companions, the location of the source on a color–magnitude diagram places it redward of its more likely position on the main sequence (MS). The source may be temporarily out of thermal equilibrium, or a cool and inflated non-MS companion, which is similar to the suggested companion of Type Ib SN 2019yvr that was constrained from pre-SN imaging. We also consider other possible physical scenarios for the source, including a fading SN, circumstellar shock interaction, line-of-sight coincidence, and an unresolved host star cluster, all of which will require future observations to more definitively rule out. Ultimately, the fraction of surviving companions (“binary fraction”) will provide necessary constraints on binary evolution models and the underlying physics.
Measurements of the growth rate of structures at
z
< 0.1 with peculiar velocity surveys have the potential of testing the validity of general relativity on cosmic scales. In this work, we present ...growth-rate measurements from realistic simulated sets of type-Ia supernovae (SNe Ia) from the Zwicky Transient Facility (ZTF). We describe our simulation methodology, the light-curve fitting, and peculiar velocity estimation. Using the maximum likelihood method, we derived constraints on
fσ
8
using only ZTF SN Ia peculiar velocities. We carefully tested the method and we quantified biases due to selection effects (photometric detection, spectroscopic follow-up for typing) on several independent realizations. We simulated the equivalent of 6 years of ZTF data, and considering an unbiased spectroscopically typed sample at
z
< 0.06, we obtained unbiased estimates of
fσ
8
with an average uncertainty of 19% precision. We also investigated the information gain in applying bias correction methods. Our results validate our framework, which can be used on real ZTF data.
Context: Type IIb supernovae (SNe) are important candidates to understand mechanisms that drive the stripping of stripped-envelope (SE) supernova (SN) progenitors. While binary interactions and their ...high incidence are generally cited to favor them as Type IIb SN progenitors, this idea has not been tested using models covering a broad parameter space. Aims: In this work we use single- and binary-star models at solar and low metallicities covering a broad parameter space to investigate the progenitors of and evolutionary pathways to Type IIb SNe. We also estimate theoretical Type IIb SN rates and make predictions for observable constraints. We also perform a case study on SN 2016gkg using Bayesian inference to derive the probability distributions of its progenitors using existing observational constraints. Methods: We use the largest database of self-consistently computed Type IIb progenitor models, statistical inference methods and multiple comparison methods to observations. As a result, our work provides the strongest constraints on Type IIb progenitors and progenitor channels to date. Results: We find that the parameter space for single-star SN IIb progenitors decreases and that for binary-star SN IIb progenitors increases with decreasing metallicity. We find that single and binary stars contribute roughly the same as Type IIb SNe at solar metallicity. Binary stars only dominate as progenitors at low metallicity. We also find that our models can account for less than half the observationally inferred rate for Type IIb SNe at high metallicity, making up < 4.5% of all core-collapse (CC) SNe. On the other hand, our models can account for the rates currently indicated by observations at low metallicity, making up 0.5–15% of all CC SNe. However, this requires low mass transfer efficiencies (~0.1) in the binaries. We find that potential binary star progenitors for SN 2016gkg have smaller pre-SN hydrogen-envelope and helium-core masses than potential single-star progenitors typically by 0.1 solar masses and 2 solar masses, respectively. We find that, a binary companion, if present, is a main-sequence or red-giant star. We demonstrate that the range of progenitor helium-core mass for SN 2016gkg inferred from observations could help improve constraints on the progenitor. We find that the probability that the progenitor of SN 2016gkg was a binary is 22% when we use constraints only on the progenitor luminosity and effective temperature. Imposing the range of pre-SN progenitor hydrogen-envelope mass and radius inferred from SN light-curves the probability the progenitor is a binary increases to 44%. However, there is no clear preference for a binary progenitor. Conclusions: We suggest that, at solar metallicity, the stellar wind mass-loss rates are lower than those currently used in most stellar evolution models. Lower mass-loss rates would widen the parameter space for binary Type IIb SNe at solar metallicity by allowing stars that initiate mass transfer earlier in their evolution to reach CC without getting stripped. Our analysis of SN 2016gkg demonstrates the importance of statistical inference methods to constrain progenitor channels. Our work indicates that to address the question of progenitors of SNe IIb we still need four pieces of information: (1) SN IIb rates as a function of metallicity, (2) better constraints on structural properties of SN IIb progenitors, (3) robust distributions for single- and binary- star properties, and (4) theoretical models for SN IIb progenitors at solar metallicity using the `correct' mass-loss prescription. Finally, our work highlights the importance of self-consistent broad parameter space modeling and statistical inference methods to constrain SN progenitor channels. Such methods will be especially important given the deluge of data expected with the imminent launch of Large Synoptic Survey Telescope (LSST).
We use cosmological simulations from the Feedback In Realistic Environments project, which implement a comprehensive set of stellar feedback processes, to study ultraviolet (UV) metal-line emission ...from the circum-galactic medium of high-redshift (z = 2–4) galaxies. Our simulations cover the halo mass range M
h ∼ 2 × 1011–8.5 × 1012 M⊙ at z = 2, representative of Lyman break galaxies. Of the transitions we analyse, the low-ionization C iii (977 Å) and Si iii (1207 Å) emission lines are the most luminous, with C iv (1548 Å) and Si iv (1394 Å) also showing interesting spatially extended structures. The more massive haloes are on average more UV-luminous. The UV metal-line emission from galactic haloes in our simulations arises primarily from collisionally ionized gas and is strongly time variable, with peak-to-trough variations of up to ∼2 dex. The peaks of UV metal-line luminosity correspond closely to massive and energetic mass outflow events, which follow bursts of star formation and inject sufficient energy into galactic haloes to power the metal-line emission. The strong time variability implies that even some relatively low-mass haloes may be detectable. Conversely, flux-limited samples will be biased towards haloes whose central galaxy has recently experienced a strong burst of star formation. Spatially extended UV metal-line emission around high-redshift galaxies should be detectable by current and upcoming integral field spectrographs such as the Multi Unit Spectroscopic Explorer on the Very Large Telescope and Keck Cosmic Web Imager.
We perform a Bayesian analysis of the mass distribution of stellar-mass black holes using the observed masses of 15 low-mass X-ray binary systems undergoing Roche lobe overflow and 5 high-mass, ...wind-fed X-ray binary systems. Using Markov Chain Monte Carlo calculations, we model the mass distribution both parametrically-as a power law, exponential, Gaussian, combination of two Gaussians, or log-normal distribution-and non-parametrically-as histograms with varying numbers of bins. We provide confidence bounds on the shape of the mass distribution in the context of each model and compare the models with each other by calculating their relative Bayesian evidence as supported by the measurements, taking into account the number of degrees of freedom of each model. The mass distribution of the low-mass systems is best fit by a power law, while the distribution of the combined sample is best fit by the exponential model. This difference indicates that the low-mass subsample is not consistent with being drawn from the distribution of the combined population. We examine the existence of a 'gap' between the most massive neutron stars and the least massive black holes by considering the value, M{sub 1%}, of the 1% quantile from each black hole mass distribution as the lower bound of black hole masses. Our analysis generates posterior distributions for M{sub 1%}; the best model (the power law) fitted to the low-mass systems has a distribution of lower bounds with M{sub 1%}>4.3 M{sub sun} with 90% confidence, while the best model (the exponential) fitted to all 20 systems has M{sub 1%}>4.5 M{sub sun} with 90% confidence. We conclude that our sample of black hole masses provides strong evidence of a gap between the maximum neutron star mass and the lower bound on black hole masses. Our results on the low-mass sample are in qualitative agreement with those of Ozel et al., although our broad model selection analysis more reliably reveals the best-fit quantitative description of the underlying mass distribution. The results on the combined sample of low- and high-mass systems are in qualitative agreement with Fryer and Kalogera, although the presence of a mass gap remains theoretically unexplained.
Astronomical surveys continue to provide unprecedented insights into the time-variable Universe and will remain the source of groundbreaking discoveries for years to come. However, their data ...throughput has overwhelmed the ability to manually synthesize alerts for devising and coordinating necessary follow-up with limited resources. The advent of Rubin Observatory, with alert volumes an order of magnitude higher at otherwise sparse cadence, presents an urgent need to overhaul existing human-centered protocols in favor of machine-directed infrastructure for conducting science inference and optimally planning expensive follow-up observations. We present the first implementation of autonomous real-time science-driven follow-up using value iteration to perform sequential experiment design. We demonstrate it for strategizing photometric augmentation of Zwicky Transient Facility Type Ia supernova light-curves given the goal of minimizing SALT2 parameter uncertainties. We find a median improvement of 2-6% for SALT2 parameters and 3-11% for photometric redshift with 2-7 additional data points in g, r and/or i compared to random augmentation. The augmentations are automatically strategized to complete gaps and for resolving phases with high constraining power (e.g. around peaks). We suggest that such a technique can deliver higher impact during the era of Rubin Observatory for precision cosmology at high redshift and can serve as the foundation for the development of general-purpose resource allocation systems.