Search

We report constraints on a variety of non-standard cosmological models using the full 5-yr photometrically classified type Ia supernova sample from the Dark Energy Survey (DES-SN5YR). Both Akaike Information Criterion (AIC) and Suspiciousness calculations find no strong evidence for or against any of the non-standard models we explore. When combined with external probes, the AIC and Suspiciousness agree that 11 of the 15 models are moderately preferred over Flat- |$\Lambda$| CDM suggesting additional flexibility in our cosmological models may be required beyond the cosmological constant. We also provide a detailed discussion of all cosmological assumptions that appear in the DES supernova cosmology analyses, evaluate their impact, and provide guidance on using the DES Hubble diagram to test non-standard models. An approximate cosmological model, used to perform bias corrections to the data holds the biggest potential for harbouring cosmological assumptions. We show that even if the approximate cosmological model is constructed with a matter density shifted by |$\Delta \Omega _{\rm m}\sim 0.2$| from the true matter density of a simulated data set the bias that arises is subdominant to statistical uncertainties. Nevertheless, we present and validate a methodology to reduce this bias. [ABSTRACT FROM AUTHOR]

We use the latest data set of supernova (SN) host galaxies to investigate how the host properties – stellar mass, star formation rate (SFR), 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 SNe feature bluer and lower luminosity hosts than most subtypes of core-collapse supernova (CC SN). Assuming simple proportionality of CC SN rates and host 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 |$\sim 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, favouring long-lived progenitors. [ABSTRACT FROM AUTHOR]

The future of time-domain optical astronomy relies on the development of techniques and software capable of handling a rising amount of data and gradually complementing, or replacing if necessary, real observations. Next-generation surveys, like the Large Synoptic Survey Telescope, will open the door to the new era of optical astrophysics, creating, at the same time, a deficiency in spectroscopic data necessary to confirm the nature of each event and to fully recover the parametric space. In this framework, we developed Core collApse Supernovae parameTers estimatOR (CASTOR), a novel software for data analysis. CASTOR combines Gaussian process and other machine-learning techniques to build time-series templates of synthetic spectra and to estimate parameters of core-collapse supernovae (CCSNe) for which only multiband photometry is available. Techniques to build templates are fully data driven and non-parametric through empirical and robust models, and rely on the direct comparison with a training set of 111 CCSNe from the literature. Furthermore, CASTOR employees the real photometric data and the reconstructed synthetic spectra of an event to estimate parameters that belong to the supernova ejecta, to the stellar progenitor and to the event itself, in a rapid and user-friendly framework. In this work, we provide a demonstration of how CASTOR works, studying available data from SN 2015ap and comparing our results with those available in literature. [ABSTRACT FROM AUTHOR]

Time-domain astrophysics continues to grow rapidly, with the inception of new surveys drastically increasing data volumes. Democratized, distributed approaches to training sets for machine learning classifiers are crucial to make the most of this torrent of discovery – with citizen science approaches proving effective at meeting these requirements. In this paper, we describe the creation of and the initial results from the Kilonova Seekers citizen science project, built to find transient phenomena from the GOTO telescopes in near real-time. Kilonova Seekers launched in 2023 July and received over 600 000 classifications from approximately 2000 volunteers over the course of the LIGO-Virgo-KAGRA O4a observing run. During this time, the project has yielded 20 discoveries, generated a 'gold-standard' training set of 17 682 detections for augmenting deep-learned classifiers, and measured the performance and biases of Zooniverse volunteers on real-bogus classification. This project will continue throughout the lifetime of GOTO, pushing candidates at ever-greater cadence, and directly facilitate the next-generation classification algorithms currently in development. [ABSTRACT FROM AUTHOR]

Nitrogen is produced by CNO-cycling in massive stars, and can be ejected in significant amounts in supernova explosions. While in H-rich SNe, its [N  ii ] 6548, 6583 emission becomes obscured by strong H |$\alpha$|⁠ , in explosions of He stars, this nitrogen emission becomes more visible. We here explore the formation of this line, using the sumo code to compute spectra for a grid of 1D models with parametrized mixing informed from new 2D simulations. Because the mass fraction of nitrogen in the ejecta decreases with larger He-core masses, as more of the He/N zone gets processed by shell helium burning and is lost to winds, the [N  ii ] luminosity relative to the overall optical flux probes the He-core mass. By comparing to large samples of data, we find that low-mass He cores (⁠|$M_{\rm preSN}\lesssim \ 3\ \mathrm{ M}_\odot$|⁠) are exclusively associated with Type IIb SNe, with the exception of Type Ib SN 2007Y. Seeing no strong nitrogen emission in other Type Ib SNe, the implication is either an origin from low-mass stars with the He/N layer (but not the He/C) layer peeled away, or from higher mass He cores. We also see no clear nitrogen emission in Type Ic SNe. We discuss the diagnostic potential of this new line metric, and also dependencies on mass-loss rate and metallicity. [ABSTRACT FROM AUTHOR]

We investigate impacts of stellar rotation and magnetic fields on black hole (BH) formation and its subsequent explosive activities, by conducting axisymmetric radiation-magnetohydrodynamics simulations of gravitational collapse of a 70  |$\mathrm{M}_\odot$| star with two-moment multi energy neutrino transport in full general relativity for the first time. Due to its dense stellar structure, all models cannot avoid the eventual BH formation even though a strongly magnetized model experiences the so-called magnetorotational explosion prior to the BH formation. One intriguing phenomenon observed in the strongly magnetized model is the formation of a relativistic jet in the post-BH formation. The relativistic jet is the outcome of a combination of strong magnetic fields and low-density materials above the BH. The jet further enhances the explosion energy beyond |$\sim 10^{52}$|  erg, which is well exceeding the gravitational overburden ahead of the shock. Our self-consistent supernova models demonstrate that rotating magnetized massive stars at the high-mass end of supernova progenitors could be a potential candidate of hypernova and long gamma-ray burst progenitors. [ABSTRACT FROM AUTHOR]

We report on Chandra X-ray observations of SN 2016jae and SN 2018cqj, both low-luminosity Type Ia supernova (SN) that showed the presence of an H line in their early optical spectrum. No X-ray emission is detected at the location of either SN. Upper limits to the luminosity of up to 2 |$\times 10^{40}$|  erg s |$^{-1}$| are calculated for each SN, depending on the assumed spectral model, temperature, and column density. This luminosity is comparable to that of another low-luminosity Type Ia SN, SN 2018fhw, that was observed with Chandra. It is generally lower than upper limits calculated for Type Ia-CSM SNe observed in X-rays, and also below that of SN 2012ca, the only Type Ia-CSM SN to have been detected in X-rays. Comparisons are made to other Type Ia SN with an H line observed in X-rays. The observations suggest that while the density into which the SN is expanding may have been high at the time the H |$\alpha$| line was detected, it had decreased considerably by the time of X-ray observations. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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

Neutron stars (NSs) are observed with high space velocities and elliptical orbits in binaries. The magnitude of these effects points to natal kicks that originate from asymmetries during the supernova (SN) explosions. Using a growing set of long-time 3D SN simulations with the Prometheus-Vertex code, we explore the interplay of NS kicks that are induced by asymmetric neutrino emission and by asymmetric mass ejection. Anisotropic neutrino emission can arise from a large-amplitude dipolar convection asymmetry inside the proto-NS (PNS) termed LESA (Lepton-number Emission Self-sustained Asymmetry) and from aspherical accretion downflows around the PNS, which can lead to anisotropic neutrino emission (absorption/scattering) with a neutrino-induced NS kick roughly opposite to (aligned with) the kick by asymmetric mass ejection. In massive progenitors, hydrodynamic kicks can reach up to more than 1300 km s−1, whereas our calculated neutrino kicks reach (55–140) km s−1 (estimated upper bounds of (170–265) km s−1) and only ∼(10–50) km s−1, if LESA is the main cause of asymmetric neutrino emission. Therefore, hydrodynamic NS kicks dominate in explosions of high-mass progenitors, whereas LESA-induced neutrino kicks dominate for NSs born in low-energy SNe of the lowest-mass progenitors, when these explode nearly spherically. Our models suggest that the Crab pulsar with its velocity of ∼160 km s−1, if born in the low-energy explosion of a low-mass, single-star progenitor, should have received a hydrodynamic kick in a considerably asymmetric explosion. Black holes, if formed by the collapse of short-lived PNSs and solely kicked by anisotropic neutrino emission, obtain velocities of only some km s−1. [ABSTRACT FROM AUTHOR]

Millisecond pulsars (MSPs) are a kind of radio pulsars with short spin periods, playing a key role in many aspects of stellar astrophysics. In recent years, some more MSPs with wide orbits (⁠|$\gt 30\, \rm d$|⁠) have been discovered, but their origin is still highly unclear. In this work, according to an adiabatic power-law assumption for the mass-transfer process, we carried out a large number of complete binary evolution computations for the formation of MSPs with wide orbits through the iron core-collapse supernova (CCSN) channel, in which a neutron star originating from a CCSN accretes matter from a red-giant star and spun up to millisecond periods. We found that this channel can form the observed MSPs with wide orbits in the range of |$30\!-\!1200\, {\rm d}$|⁠ , in which the white dwarf (WD) companions have masses in the range of |$0.28\!-\!0.55\, \rm M_{\odot }$|⁠. We also found that almost all the observed MSPs can be reproduced by this channel in the WD companion mass versus orbital period diagram. We estimate that the Galactic numbers of the resulting MSPs from the CCSN channel are in the range of |$\sim 4.8\!-\!8.5\times 10^{5}$|⁠. Compared with the accretion-induced collapse channel, the CCSN channel provides a main way to produce MSPs with wide orbits. [ABSTRACT FROM AUTHOR]

We present the first discoveries of the double-lined double white dwarf (DBL) survey that targets overluminous sources with respect to the canonical white dwarf cooling sequence according to a set of well-defined criteria. The primary goal of the DBL survey is to identify compact double white dwarf binary star systems from a unique spectral detection of both stars, which then enables a precise quantification of the atmospheric parameters and radial velocity variability of a system. Our search of 117 candidates that were randomly selected from a magnitude-limited sample of 399 yielded a 29 per cent detection efficiency with 34 systems exhibiting a double-lined signature. A further 38 systems show strong evidence of being single-lined or potentially DBL binaries and seven single-lined sources from the full observed sample are radial velocity variable. The 45 remaining candidates appear as a single WD with no companion or a non-DA white dwarf, bringing the efficiency of detecting binaries to 62 per cent. Atmospheric fitting of all double-lined systems reveals a large fraction that have two similar mass components that combine to a total mass of 1.0–1.3  |$\mathrm{M}_\odot$| – a class of double white dwarf binaries that may undergo a sub-Chandrasekhar mass type Ia detonation or merge to form a massive O/Ne WD, although orbital periods are required to infer on which time-scales. One double-lined system located 49 pc away, WDJ181058.67+311940.94, is super-Chandrasekhar mass, making it the second such double white dwarf binary to be discovered. [ABSTRACT FROM AUTHOR]

The role of type Ia supernovae (SN Ia), mainly the delay time distributions (DTDs) determined by the binary systems, and the yields of elements created by different explosion mechanisms, are studied by using the MulChem chemical evolution model applied to our Galaxy. We explored 15 DTDs and 12 tables of elemental yields produced by different SN Ia explosion mechanisms, doing a total of 180 models. Chemical abundances for |$\alpha$| -elements (O, Mg, Si, and Ca) and Fe derived from these models are compared with recent solar region observational data of |$\alpha$| -elements over Fe relative abundances, [X/Fe], as a function of [Fe/H] and age. A multidimensional maximum-likelihood analysis shows that 52 models are able to fit all these data sets simultaneously, considering the 1 |$\sigma$| confidence level. The combination of STROLG1 DTD from Strolger et al. (2020) and LN20181 SN Ia yields from Leung & Nomoto (2018) provides the best fit. The exponential model with very prompt events is a possible DTD, but a combination of several channels is more probable. The SN Ia yields that include MCh or Near MCh correspond to 39 (75 per cent) of the 52 best models. Regarding the DTD, 31 (60 per cent) of the 52 most probable models correspond to the SD scenario, while the remaining 21 (40 per cent) are based on the DD scenario. Our results also show that the relatively large dispersion of the observational data may be explained by the stellar migration from other radial regions, and/or perhaps a combination of DTDs and explosion channels. [ABSTRACT FROM AUTHOR]

We present the results of a series of 3D special relativistic hydrodynamic simulations of a gamma-ray burst (GRB) jet in a massive circumstellar medium (CSM) surrounding the progenitor star. Our simulations reproduce the jet morphology transitioning from a well-collimated state to a thermal pressure-driven state for a range of CSM masses and outer radii. The jet–CSM interaction redistributes the jet energy to materials expanding into a wide solid angle and results in a quasi-spherical ejecta with four-velocities from |$\Gamma \beta \simeq 0.1$| to 10. The mass and kinetic energy of the ejecta with velocities faster than |$0.1c$| are typically of the order of |$0.1\, M_{\odot }$| and |$10^{51}\:\mbox{erg}$| with only a weak dependence on the CSM mass and radius for the explored CSM parameter ranges. We find that the numerically obtained density structure of the mildly relativistic ejecta is remarkably universal. The radial density profile is well approximated as a power-law function of the radial velocity with an index of |$-5$|⁠ , |$\rho \propto v^{-5}$|⁠ , in agreement with our previous simulations and other studies, as well as those suggested from recent studies on early-phase spectra of supernovae associated with GRBs. Such fast ejecta rapidly becomes transparent following its expansion. Gradually releasing the trapped thermal photons, the ejecta gives rise to bright UV–optical emission within |${\sim} 1\:$| d. We discuss the potential link of the relativistic ejecta resulting from jet–CSM interaction to GRB-associated supernovae as well as fast and blue optical transients. [ABSTRACT FROM AUTHOR]

We report isotope data for C, N, Al, Si, and S of 33 presolar SiC and Si3N4 grains (0.3–1.6  |$\mu$| m) of Type X, C, D, and N from the Murchison CM2 meteorite of likely core-collapse supernova (CCSN) origin which we discuss together with data of six SiC X grains from an earlier study. The isotope data are discussed in the context of hydrogen ingestion supernova (SN) models. We have modified previously used ad-hoc mixing schemes in that we considered (i) heterogeneous H ingestion into the He shell of the pre-SN star, (ii) a variable C-N fractionation for the condensation of SiC grains in the SN ejecta, and (iii) smaller mass units for better fine-tuning. With our modified ad-hoc mixing approach over small scales (0.2–0.4 M⊙), with major contributions from the O-rich O/nova zone, we find remarkably good fits (within a few per cent) for 12C/13C, 26Al/27Al, and 29Si/28Si ratios. The 14N/15N ratio of SiC grains can be well matched if variable C-N fractionation is considered. However, the Si3N4 isotope data point to overproduction of 15N in hydrogen ingestion CCSN models and lower C-N fractionation during SiC condensation than applied here. Our ad-hoc mixing approach based on current CCSN models suggests that the O-rich O/nova zone, which uniquely combines explosive H- and He-burning signatures, is favourable for SiC and Si3N4 formation. The effective range of C/O abundance variations in the He shell triggered by H ingestion events in the massive star progenitor is currently not well constrained and needs further investigation. [ABSTRACT FROM AUTHOR]

Context. Supernova (SN) 2023ixf was discovered in the galaxy M 101 in May 2023. Its proximity provided the scientific community an extremely valuable opportunity to study the characteristics of the SN and its progenitor. A point source detected on archival images and hydrodynamical modeling of the bolometric light curve have been used to constrain the former star's properties. There is a significant variation in the published results regarding the initial mass of the progenitor. Nebular spectroscopy can be used to enhance our understanding of the SN and its progenitor. Aims. We determined the SN progenitor mass by studying the first published nebular spectrum, taken 259 days after the explosion. Methods. We analyzed the nebular spectrum taken with GMOS at the Gemini North Telescope. We identified typical emission lines, such as [O I], Hα, and [Ca II], among others. Some species' line profiles show broad and narrow components, indicating two ejecta velocities and an asymmetric ejecta. We inferred the progenitor mass of SN 2023ixf by comparing its spectra with synthetic spectra and by measuring the forbidden oxygen doublet flux. Results. Based on the flux ratio and the direct comparison with spectra models, the progenitor star of SN 2023ixf had a MZAMS between 12 and 15 M⊙. We find that using the [O I] doublet flux provides a less tight constraint on the progenitor mass. Our results agree with those from hydrodynamical modeling of the early light curve and pre-explosion image estimates that point to a relatively low-mass progenitor. [ABSTRACT FROM AUTHOR]

Aims. Recent observations of core-collapse supernovae (SNe) suggest aspherical explosions. Globally, aspherical structures in SN explosions are thought to encode information regarding the underlying explosion mechanism. However, the exact explosion geometries from the inner cores to the outer envelopes are poorly understood. Methods. Here, we present photometric, spectroscopic, and polarimetric observations of the Type IIP SN 2021yja and discuss its explosion geometry in comparison to those of other Type IIP SNe that show large-scale aspherical structures in their hydrogen envelopes (SNe 2012aw, 2013ej and 2017gmr). Results. During the plateau phase, SNe 2012aw and 2021yja exhibit high continuum polarization characterized by two components with perpendicular polarization angles. This behavior can be interpreted as being due to a bipolar explosion, where the SN ejecta is composed of a polar (energetic) component and an equatorial (bulk) component. In such a bipolar explosion, an aspherical axis created by the polar ejecta would dominate at early phases, while the perpendicular axis along the equatorial ejecta would emerge at late phases after the photosphere in the polar ejecta has receded. Our interpretation of the explosions in SNe 2012aw and 2021yja as bipolar is also supported by other observational properties, including the time evolution of the line velocities and the line shapes in the nebular spectra. The polarization of other Type IIP SNe that show large-scale aspherical structures in the hydrogen envelope (SNe 2013ej and 2017gmr) is also consistent with the bipolar-explosion scenario, although this is not conclusive. [ABSTRACT FROM AUTHOR]

Aims. In this work we extended an energy-integrated neutrino transport method to facilitate efficient, yet precise, modeling of compact astrophysical objects. We particularly focus on core-collapse supernovae. Methods. We implemented a gray neutrino-transport framework from the literature into FLASH and performed a detailed evaluation of its accuracy in core-collapse supernova simulations. Based on comparisons with results from simulations using energy-dependent neutrino transport, we incorporated several improvements to the original scheme. Results. Our analysis shows that our gray neutrino transport method successfully reproduces key aspects from more complex energy-dependent transport across a variety of progenitors and equations of state. We find both qualitative and reasonable quantitative agreement with multi-group M1 transport simulations. However, the gray scheme tends to slightly favor shock revival. In terms of gravitational wave and neutrino signals, there is a good alignment with the energy-dependent transport, although we find 15–30% discrepancies in the average energy and luminosity of heavy-lepton neutrinos. Simulations using the gray transport are around four times faster than those using energy-dependent transport. [ABSTRACT FROM AUTHOR]

Narrow absorption features in nearby supernova (SN) spectra are a powerful diagnostic of the slow-moving material in the line of sight: they are extensively used to infer dust extinction from the host galaxies, and they can also serve in the detection of circumstellar material originating from the SN progenitor and present in the vicinity of the explosion. Despite their wide use, very few studies have examined the biases of the methods to characterize narrow lines, and not many statistical analyses exist. This is the first paper of a series in which we present a statistical analysis of narrow lines of SN spectra of various resolutions. We developed a robust automated methodology to measure the equivalent width (EW) and velocity of narrow absorption lines from intervening material in the line of sight of SNe, including Na I D, Ca II H&K, K I, and diffuse interstellar bands. We carefully studied systematic biases in heterogeneous spectra from the literature by simulating different signal-to-noise, spectral resolution, size and orientation of the slit, and we present the real capabilities and limitations of using low- and mid-resolution spectra to study these lines. In particular, we find that the measurement of the EW of the narrow lines in low-resolution spectra is highly affected by the evolving broad P-Cygni profiles of the SN ejecta, both for core-collapse and type Ia SNe, inducing a conspicuous apparent evolution. Such pervading non-physical evolution of narrow lines might lead to wrong conclusions on the line-of-sight material, for example concerning circumstellar material ejected from the SN progenitors. We thus present an easy way to detect and exclude those cases to obtain more robust and reliable measurements. Finally, after considering all possible effects, we analysed the temporal evolution of the narrow features in a large sample of nearby SNe to detect any possible variation in their EWs over time. We find no time evolution of the narrow line features in our large sample for all SN types. [ABSTRACT FROM AUTHOR]

Context. The analysis of non-variable stars is generally neglected in the literature. However, such objects are needed for many calibration processes and for testing pulsational models. The photometric time series of the Kepler satellite mission still stand as the most accurate data available today and are excellently suited to the search for non-variable stars. Aims. We analysed all long-cadence light curves for stars not reported as a variable so far from the Kepler satellite mission. Using the known characteristics and flaws of these data sets, we defined three different frequency ranges where we searched for non-variability. Methods. We used the Lomb–Scargle periodogram and the false-alarm probability (FAP) to analyse the cleaned data sets of 138 451 light curves. We then used log FAP ≥ −2 to define a star as 'non-variable' in the ranges below 0.1 c/d, 0.1 to 2.0 c/d, and 2.0 to 25.0 c/d, respectively. Furthermore, we also calculated the standard deviation of the mean light curve to obtain another parameter. Results. In total, we found 14 154 stars that fulfil the set criteria. These objects are mostly cooler than the 7000 K populating the whole main sequence (MS) to the red giant branch (RGB). [ABSTRACT FROM AUTHOR]

We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M⊙ star with spectral neutrino transport. Our focus is to study the gravitational-wave (GW) signatures from the magnetorotationally (MR)-driven models. By parametrically changing the initial angular velocity and the strength of the magnetic fields in the core, we compute four models. Among our models, only those with cores having an initial magnetic field strength of 1012 G and rotation rates of 1 or 2 rad s−1 produce MHD jets. Seen from the direction perpendicular to the rotational axis, a characteristic waveform is obtained exhibiting a monotonic time increase in the wave amplitude. As previously identified, this stems from the propagating MHD outflows along the axis. We show that the GW amplitude from anisotropic neutrino emission becomes more than one order-of-magnitude bigger than that from the matter contribution, whereas seen from the rotational axis, both of the two components are in the same order-of-magnitudes. Due to the memory effect, the frequency of the neutrino GW from our full-fledged 3D-MHD models is in the range less than ∼10 Hz. Toward the future GW detection for a Galactic core-collapse supernova, if driven by the MR mechanism, the planned next-generation detector as DECIGO is urgently needed to catch the low-frequency signals. [ABSTRACT FROM AUTHOR]

Multiple explosion mechanisms have been proposed to explain type Ia supernovae (SNe Ia). Empirical modelling tools have also been developed that allow for fast, customized modelling of individual SNe and direct comparisons between observations and explosion model predictions. Such tools have provided useful insights, but the subjective nature with which empirical modelling is performed makes it difficult to obtain robust constraints on the explosion physics or expand studies to large populations of objects. Machine learning accelerated tools have therefore begun to gain traction. In this paper, we present riddler , a framework for automated fitting of SNe Ia spectral sequences up to shortly after maximum light. We train a series of neural networks on realistic ejecta profiles predicted by the W7 and N100 explosion models to emulate full radiative transfer simulations and apply nested sampling to determine the best-fitting model parameters for multiple spectra of a given SN simultaneously. We show that riddler is able to accurately recover the parameters of input spectra and use it to fit observations of two well-studied SNe Ia. We also investigate the impact of different weighting schemes when performing quantitative spectral fitting and show that best-fitting models and parameters are highly dependent on the assumed weighting schemes and priors. As spectroscopic samples of SNe Ia continue to grow, automated spectral fitting tools such as riddler will become increasingly important to maximise the physical constraints that can be gained in a quantitative and consistent manner. [ABSTRACT FROM AUTHOR]

Nuclear reactions are key to our understanding of stellar evolution, particularly the |$^{12}{\rm C}(\alpha ,\gamma)^{16}{\rm O}\,$| rate, which is known to significantly influence the lower and upper ends of the black hole (BH) mass distribution due to pair-instability supernovae (PISNe). However, these reaction rates have not been sufficiently determined. We use the mesa stellar evolution code to explore the impact of uncertainty in the |$^{12}{\rm C}(\alpha ,\gamma)^{16}{\rm O}\,$| rate on PISN explosions, focusing on nucleosynthesis and explosion energy by considering the high resolution of the initial mass. Our findings show that the mass of synthesized radioactive nickel (56Ni) and the explosion energy increase with |$^{12}{\rm C}(\alpha ,\gamma)^{16}{\rm O}\,$| rate for the same initial mass, except in the high-mass edge region. With a high (about twice the starlib standard value) rate, the maximum amount of nickel produced falls below 70 M⊙, while with a low rate (about half of the standard value) it increases up to 83.9 M⊙. These results highlight that carbon 'preheating' plays a crucial role in PISNe by determining core concentration when a star initiates expansion. Our results also suggest that the onset of the expansion, which means the end of compression, competes with collapse caused by helium photodisintegration, and the maximum mass that can lead to an explosion depends on the |$^{12}{\rm C}(\alpha ,\gamma)^{16}{\rm O}\,$| reaction rate. [ABSTRACT FROM AUTHOR]

We apply the hierarchical probabilistic spectral energy distribution (SED) model bayesn to analyse a sample of 475 type Ia supernovae (0.015 < z < 0.4) from Foundation, DES3YR and PS1MD to investigate the properties of dust in their host galaxies. We jointly infer the dust law RV population distributions at the SED level in high- and low-mass galaxies simultaneously with dust-independent, intrinsic differences. We find an intrinsic mass step of −0.049 ± 0.016 mag, at a significance of 3.1σ, when allowing for a constant intrinsic, achromatic magnitude offset. We additionally apply a model allowing for time- and wavelength-dependent intrinsic differences between SNe Ia in different mass bins, finding ∼2σ differences in magnitude and colour around peak and 4.5σ differences at later times. These intrinsic differences are inferred simultaneously with a difference in population mean RV of ∼2σ significance, demonstrating that both intrinsic and extrinsic differences may play a role in causing the host galaxy mass step. We also consider a model which allows the mean of the RV distribution to linearly evolve with redshift but find no evidence for any evolution – we infer the gradient of this relation η R  = −0.38 ± 0.70. In addition, we discuss in brief a new, GPU-accelerated python implementation of bayesn suitable for application to large surveys which is publicly available and can be used for future cosmological analyses; this code can be found here: https://github.com/bayesn/bayesn. [ABSTRACT FROM AUTHOR]

Transient surveys routinely detect supernovae (SNe) without obvious host galaxies. To understand the demographics of these 'hostless' SNe and to constrain the possible host properties, we identify 161 SNe reported to the Transient Name Server since 2016 that do not have hosts catalogued from pre-explosion wide-field galaxy surveys. Using forced aperture photometry, we detect excess flux around only 56 of these SNe. Both thermonuclear and core-collapse (CC) SNe are present in our sample. Compared to flux-limited SNe samples with known hosts, superluminous supernovae (SLSNe), particularly hydrogen-deficient SLSNe, are over-represented here relative to all other SNe types; among CC SNe, there is also a higher fraction of interacting SNe than non-interacting. On the low-luminosity side, seven SNe have host absolute magnitude upper limits fainter than Mg  = −12, about 1 per cent of the Small Magellanic Cloud's luminosity; the faintest limits are close to the luminosity of globular clusters or ultra-faint dwarf galaxies (Mg ≃ −8). Fitting multiband forced photometry, 11 SNe have host stellar masses <106 M⊙ assuming quiescent hosts, and 13 SNe have host stellar masses <105 M⊙ assuming star-forming hosts. The spatial distribution of hostless SNe indicates that the majority are not associated with known galaxy groups and clusters, ruling out intracluster stellar light as the primary contributor of such SNe. Hostless Type Ia SNe tend to be more luminous and slow-fading than SNe Ia with known host galaxies, implying a hidden population of low-mass and star-forming hosts. We conclude that any undetected host galaxies are likely star-forming dwarfs in the field. [ABSTRACT FROM AUTHOR]

It is widely accepted that the width–luminosity relation used to standardize normal Type Ia supernovae (SNe Ia) breaks down in underluminous, 1991bg-like SNe Ia. This breakdown may be due to the choice of parameter used as a stand-in for the width of the SN Ia light curve. Using the colour stretch parameter s BV instead of older parameters resolves this issue. Here, I assemble a sample of 14 nearby 1991bg-like SNe Ia from the literature, all of which have independent host-galaxy distance moduli and little to no reddening. I use Gaussian process regression to fit the light curves of these SNe in U / u, B, V, g, R / r, I / i , and H , and measure their peak absolute magnitudes. I find statistically significant (>5σ confidence level in the optical and >4σ in the near-infrared) correlations between the peak absolute magnitudes of the 1991bg-like SNe Ia and their s BV values in the range 0.2 < s BV < 0.6. These correlations are broadly consistent with fits to s BV < 0.7 SNe Ia with preliminary B - and V -band peak absolute magnitudes from the Carnegie Supernova Project and significantly inconsistent with similar fits to normal and transitional SNe Ia (with 0.7 < s BV < 1.1). The underluminous width–luminosity relation shown here needs to be properly calibrated with a homogeneous sample of 1991bg-like SNe Ia, after which it could be used as a rung on a new cosmological distance ladder. With surface-brightness fluctuations (or another non-Cepheid method) used to calibrate distances to nearby 1991bg-like SNe Ia, such a ladder could produce an independent measurement of the Hubble–Lemaître constant, H 0. [ABSTRACT FROM AUTHOR]

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 | > 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, z med = 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. [ABSTRACT FROM AUTHOR]

Electron-capture-supernovae (EC-SNe) provide an alternative channel for producing neutron stars (NSs). They play an important role in the formation of double NS (DNS) systems and the chemical evolution of galaxies, and contribute to the NS mass distribution in observations. It is generally believed that EC-SNe originate from e -captures on |$\rm ^{24}Mg$| and |$\rm ^{20}Ne$| in the massive degenerate oxygen–neon (ONe) cores with masses close to the Chandrasekhar limit (M Ch). However, the origin of EC-SNe is still uncertain. In this paper, we systematically studied the EC-SNe in NS + He star systems by considering the explosive oxygen burning that may occur in the near- M Ch ONe core. We provided the initial parameter spaces for producing EC-SNe in the initial orbital period − initial He star mass (log |$P_{\rm orb}^{\rm i}-M_{\rm He}^{\rm i}$|⁠) diagram, and found that both |$M_{\rm He}^{\rm i}$| and minimum |$P_{\rm orb}^{\rm i}$| for EC-SNe increase with metallicity. Then, by considering NS kicks added to the newborn NS, we investigated the properties of the formed DNS systems after the He star companions collapse into NSs, such as the orbital periods, eccentricities, and spin periods of recycle pulsars (P spin), etc. The results show that most of the observed DNS systems can be produced by NS kicks of |$\lesssim$| 50 km s−1. In addition, we found that NSs could accrete more material if the residual H envelope on the He star companions is considered, which can form the mildly recycled pulsars (P spin ∼ 20 ms) in DNS systems. [ABSTRACT FROM AUTHOR]

We use rest-frame optical and near-infrared (NIR) observations of 42 Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project at low- z and 37 from the RAISIN (SNIA in the IR) Survey at high- z to investigate correlations between SN Ia host galaxy dust, host mass, and redshift. This is the first time the SN Ia host galaxy dust extinction law at high- z has been estimated using combined optical and rest-frame NIR data (YJ band). We use the BayeSN hierarchical model to leverage the data's wide rest-frame wavelength range (extending to ∼1.0–1.2 μ m for the RAISIN sample at 0.2 ≲ z ≲ 0.6). By contrasting the RAISIN and Carnegie Supernova Project (CSP) data, we constrain the population distributions of the host dust R V parameter for both redshift ranges. We place a limit on the difference in population mean R V between RAISIN and CSP of −1.16 < Δμ(R V) < 1.38 with 95 per cent posterior probability. For RAISIN we estimate μ(R V) = 2.58 ± 0.57, and constrain the population standard deviation to σ(R V) < 0.90 [2.42] at the 68 [95] per cent level. Given that we are only able to constrain the size of the low- to high- z shift in μ(R V) to ≲1.4 – which could still propagate to a substantial bias in the equation-of-state parameter w – these and other recent results motivate continued effort to obtain rest-frame NIR data at low- and high-redshifts (e.g. using the Roman Space Telescope). [ABSTRACT FROM AUTHOR]

We present gausSN , a Bayesian semiparametric Gaussian Process (GP) model for time-delay estimation with resolved systems of gravitationally lensed supernovae (glSNe). gausSN models the underlying light curve non-parametrically using a GP. Without assuming a template light curve for each SN type, gausSN fits for the time delays of all images using data in any number of wavelength filters simultaneously. We also introduce a novel time-varying magnification model to capture the effects of microlensing alongside time-delay estimation. In this analysis, we model the time-varying relative magnification as a sigmoid function, as well as a constant for comparison to existing time-delay estimation approaches. We demonstrate that gausSN provides robust time-delay estimates for simulations of glSNe from the Nancy Grace Roman Space Telescope and the Vera C. Rubin Observatory's Legacy Survey of Space and Time (Rubin-LSST). We find that up to 43.6 per cent of time-delay estimates from Roman and 52.9 per cent from Rubin-LSST have fractional errors of less than 5  per cent. We then apply gausSN to SN Refsdal and find the time delay for the fifth image is consistent with the original analysis, regardless of microlensing treatment. Therefore, gausSN maintains the level of precision and accuracy achieved by existing time-delay extraction methods with fewer assumptions about the underlying shape of the light curve than template-based approaches, while incorporating microlensing into the statistical error budget. gausSN is scalable for time-delay cosmography analyses given current projections of glSNe discovery rates from Rubin-LSST and Roman. [ABSTRACT FROM AUTHOR]

Type Iax supernovae (SNe Iax) are proposed to arise from deflagrations of Chandrasekhar mass white dwarfs (WDs). Previous deflagration simulations have achieved good agreement with the light curves and spectra of intermediate-luminosity and bright SNe Iax. However, the model light curves decline too quickly after peak, particularly in red optical and near-infrared (NIR) bands. Deflagration models with a variety of ignition configurations do not fully unbind the WD, leaving a remnant polluted with 56Ni. Emission from such a remnant may contribute to the luminosity of SNe Iax. Here we investigate the impact of adding a central energy source, assuming instantaneous powering by 56Ni decay in the remnant, in radiative transfer calculations of deflagration models. Including the remnant contribution improves agreement with the light curves of SNe Iax, particularly due to the slower post-maximum decline of the models. Spectroscopic agreement is also improved, with intermediate-luminosity and faint models showing greatest improvement. We adopt the full remnant 56Ni mass predicted for bright models, but good agreement with intermediate-luminosity and faint SNe Iax is only possible for remnant 56Ni masses significantly lower than those predicted. This may indicate that some of the 56Ni decay energy in the remnant does not contribute to the radiative luminosity but instead drives mass ejection, or that escape of energy from the remnant is significantly delayed. Future work should investigate the structure of remnants predicted by deflagration models and the potential roles of winds and delayed energy escape, as well as extend radiative transfer simulations to late times. [ABSTRACT FROM AUTHOR]

Recent studies of the neutrino-driven wind from proto-neutron stars have indicated that the wind is likely proton rich for much of its lifetime, and the high flux of neutrinos can induce νp-process nucleosynthesis allowing for the formation of heavy elements. It has also been shown that gravito-acoustic waves, generated by convection within the proto-neutron star, can significantly alter the dynamics and nucleosynthesis in the wind. Therefore, we present a study of the effects of convection-driven waves on the nucleosynthesis in proton-rich neutrino-driven winds, focusing on the νp-process. We find that wave effects can strongly impact νp-process nucleosynthesis even at wave luminosities a factor of 10−5 smaller than the total neutrino luminosity. The momentum flux of the waves accelerates the wind, reducing the net neutrino heating and the persistent neutron abundance created by p (⁠|$\bar{\nu }_\mathrm{ e},e^+$|⁠), which impedes νp-process nucleosynthesis. However, this effect is generally counteracted by the effects of waves on seed nucleus formation, as the acceleration of the wind and the heating that occurs as these waves shock both favour a more α-rich environment with very little heavy seed nucleus formation. Overall, higher wave luminosities correlate (albeit non-monotonically) with heavier element νp-processing, up to A ≈ 200 in some cases. At very high wave luminosities (≳10−3 L ν), early shock heating by the waves disrupts α recombination, and drives a suppressed, fast-outflow r-process proceeding up to A ≈ 200. This occurs despite an assumed neutrino spectrum that predicts a proton-rich wind with equilibrium Y e = 0.6. [ABSTRACT FROM AUTHOR]

Over the last two decades there have been considerable advances in modelling the spectra of massive stars and supernovae (SNe). Despite this progress, there are still numerous uncertainties that affect the accuracy of models. For massive stars, convection, instabilities, clumping, and our inability to model stellar winds self-consistently likely introduce systematic errors into our analyses. For SNe, and particularly for core-collapse SNe, departures from spherical symmetry strongly affect observed spectra and need to be taken into account. There are also issues with clumping, and mixing processes (both in the progenitor and the SN explosion) that need to be resolved. For both massive stars and SNe, the accuracy and availability of atomic data continues to be an ongoing issue influencing analyses. [ABSTRACT FROM AUTHOR]

Long gamma-ray bursts (LGRBs) and superluminous supernovae (SLSNe) are expected to result from massive star deaths. However to date, there has been no direct observational measurement of their cloud collapse timescales nor progenitor lifetimes to help constrain their mass. Our analyses of z 2 LGRB afterglow spectra and Hubble Space Telescope images find a higher fraction of host galaxies that are interacting, have a close companion, and/or may have experienced a recent galaxy 'fly by' as compared to the general z 2 galaxy population. A smaller set of z 2 SLSNe suggests a similar result. Under the hypothesis that galaxy interactions induce cloud collapse and star formation near their closest approach, we explore measurements of the host and companion galaxy velocities and separations at the time of the LRGB/SLSN event as a direct physical means to measure the timescale of cloud collapse plus progenitor star lifetime. [ABSTRACT FROM AUTHOR]

Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries. [ABSTRACT FROM AUTHOR]

Local HII environment metallicities of 65 supernovae (SNe), obtained with INT/IDS, have been determined using the N2 and O3N2 strong emission line methods. Resulting cumulative distribution functions reveal a narrower distribution for Ib SNe (standard deviation σ ∼ 0.06 dex) compared to Ic and IIP distributions (σ ∼ 0.15 dex). This narrow distribution of Ib SNe is confirmed with an extended dataset using data from Galbany et al. (2018). Statistical tests confirm a statistically significant difference between the Ib and II-P metallicity distributions with < 5% probability that they result from the same progenitors. This narrow distribution suggests a lack of Type Ib SNe in low metallicity environments and points towards single star progenitors for these Type Ib SNe, rather than binaries. It also suggests that single massive stars at low metallicity are not commonly able to produce helium-rich Type Ib supernovae. [ABSTRACT FROM AUTHOR]

Supernova properties in radio strongly depend on their circumstellar environment and they are an important probe to investigate the mass loss of supernova progenitors. Recently, core-collapse supernova observations in radio have been assembled and the rise time and peak luminosity distribution of core-collapse supernovae in radio has been obtained. In this talk, we will discuss the constraints on the mass-loss prescriptions of red supergiants obtained from the assembled radio properties of Type II supernovae. We take a couple of mass-loss prescriptions for red supergiants, calculate the rise time and peak luminosity distribution based on them, and compare the results with the observed distribution. We found that the widely spread radio rise time and peak luminosity distribution of Type II supernovae can only be explained by mass-loss prescriptions having strong dependence on the luminosity. Red supergiant mass-loss prescriptions should have steep luminosity dependence in the supernova progenitor range. [ABSTRACT FROM AUTHOR]

Despite the important role mass loss in the red supergiant phase plays in controlling stellar evolution and massive stars' final supernova fates, a theoretical explanation of the mechanism driving this mass loss has been elusive. In this contribution we present a recent breakthrough [Kee et∼al., 2021] showing that turbulent pressure alone is sufficient to markedly extend the atmospheres of red supergiants and allow a wind to be launched. The resulting theory provides a fully analytic prescription for red supergiant mass-loss rates. Moreover, the theoretical mass-loss rates computed from observationally inferred turbulent velocities are in overall good agreement with observationally inferred red supergiant mass loss. A particularly interesting aspect of this theory is that it is not sensitive to metallicity, providing important implications for stellar evolution and the so-called "red-supergiant problem" for supernova progenitors in various environments. [ABSTRACT FROM AUTHOR]

Recent observations of Type II supernovae have revealed that their red-supergiant progenitors lose a significant amount of mass during the last years of their evolution. However, because it is difficult to discover supernovae within days of explosion, the diversity of mass loss in red supergiants has not yet been fully mapped. This talk presented the case of SN 2021yja, which was serendipitously imaged within hours of explosion and observed with a sub-day cadence during its rise to peak. From the exceptionally long plateau period and the high nickel mass, we infer a relatively massive red-supergiant progenitor star. However, archival imaging from the Hubble Space Telescope places a stringent upper limit of on its progenitor mass. We discuss these conflicting constraints in the context of the larger sample of exploding red supergiants. Our analysis helps illuminate the poorly understood mechanism(s) behind red-supergiant mass loss. [ABSTRACT FROM AUTHOR]

Over the last decade, evidence has accumulated that massive stars do not typically evolve in isolation but instead follow a tumultuous journey with a companion star on their way to core collapse. While Roche-lobe overflow appears instrumental for the production of a large fraction of Type Ib and Ic supernovae (SNe), variations in the initial orbital period, Pinit, of massive interacting binaries may also produce a wide diversity of case B, BC, or C systems, with pre-SN stars endowed from minute to massive H-rich envelopes. Focusing here on the explosion of the primary donor star, originally 12.6 M⊙, we used radiation hydrodynamics and nonlocal thermodynamic equilibrium time-dependent radiative transfer to document the gas and radiation properties of such SNe, covering Types Ib, IIb, II-L, and II-P. Variations in Pinit are the root cause of the wide diversity of our SN light curves, which present single-peak, double-peak, fast-declining, or plateau-like morphologies in the V band. The different ejecta structures, expansion rates, and relative abundances (e.g., H, He, and 56Ni) can lead to a great deal of diversity in terms of spectral line shapes (absorption versus emission strength and width) and evolution. We emphasize that Hα is a key tracer of these modulations, and that He I 7065 Å is an enduring optical diagnostic for the presence of He. Our grid of simulations fares well against representative Type Ib, IIb, and II-P SNe, but interaction with circumstellar material, which is ignored in this work, is likely at the origin of the tension between our Type II-L SN models and observations (e.g., of SN 2006Y). Remaining discrepancies in the rise time to bolometric maximum of our models call for a proper account of both small-scale and large-scale structures in core-collapse SN ejecta. Discrepant Type II-P SN models, with a high plateau brightness but small spectral line widths, can be fixed by adopting more compact red-supergiant star progenitors. [ABSTRACT FROM AUTHOR]

We present an optical photometric and spectroscopic analysis of the fast-declining hydrogen-rich Type II supernova (SN) 2019nyk. The light curve properties of SN 2019nyk align well with those of other fast-declining Type II SNe, such as SNe 2013by and 2014G. SN 2019nyk exhibits a peak absolute magnitude of −18.09 ± 0.17 mag in the V band, followed by a rapid decline at 2.84  ±  0.03 mag (100 d)−1 during the recombination phase. The early spectra of SN 2019nyk exhibit high-ionisation emission features as well as narrow H Balmer lines, persisting until 4.1 d since explosion, indicating the presence of circumstellar material (CSM) in close proximity. A comparison of these features with other Type II SNe displaying an early interaction reveals similarities between these features and those observed in SNe 2014G and 2023ixf. We also compared the early spectra to literature models, estimating a mass-loss rate of the order of 10−3M⊙ yr−1. Radiation hydrodynamical modelling of the light curve also suggests the mass loss from the progenitor within a short period prior to explosion, totalling 0.16 M⊙ of material within 2900 R⊙ of the progenitor. Furthermore, light curve modelling infers a zero-age main sequence mass of 15 M⊙ for the progenitor, a progenitor radius of 1031 R⊙, and an explosion energy of 1.1 × 1051 erg. [ABSTRACT FROM AUTHOR]

SN 2020zbf is a hydrogen-poor superluminous supernova (SLSN) at z = 0.1947 that shows conspicuous C II features at early times, in contrast to the majority of H-poor SLSNe. Its peak magnitude is Mg = −21.2 mag and its rise time (≲26.4 days from first light) places SN 2020zbf among the fastest rising type I SLSNe. We used spectra taken from ultraviolet (UV) to near-infrared wavelengths to identify spectral features. We paid particular attention to the C II lines as they present distinctive characteristics when compared to other events. We also analyzed UV and optical photometric data and modeled the light curves considering three different powering mechanisms: radioactive decay of 56Ni, magnetar spin-down, and circumstellar medium (CSM) interaction. The spectra of SN 2020zbf match the model spectra of a C-rich low-mass magnetar-powered supernova model well. This is consistent with our light curve modeling, which supports a magnetar-powered event with an ejecta mass Mej = 1.5 M⊙. However, we cannot discard the CSM-interaction model as it may also reproduce the observed features. The interaction with H-poor, carbon-oxygen CSM near peak light could explain the presence of C II emission lines. A short plateau in the light curve around 35–45 days after peak, in combination with the presence of an emission line at 6580 Å, can also be interpreted as being due to a late interaction with an extended H-rich CSM. Both the magnetar and CSM-interaction models of SN 2020zbf indicate that the progenitor mass at the time of explosion is between 2 and 5 M⊙. Modeling the spectral energy distribution of the host galaxy reveals a host mass of 108.7M⊙, a star formation rate of 0.24−0.12+0.41M⊙ yr−1, and a metallicity of ∼0.4 Z⊙. [ABSTRACT FROM AUTHOR]

Context. The features in the light curves and spectra of many Type I and Type II supernovae (SNe) can be understood by assuming an interaction of the SN ejecta with circumstellar matter (CSM) surrounding the progenitor star. This suggests that many massive stars may undergo various degrees of envelope stripping shortly before exploding, and may therefore produce a considerable diversity in their pre-explosion CSM properties. Aims. We explore a generic set of about 100 detailed massive binary evolution models in order to characterize the amount of envelope stripping and the expected CSM configurations. Methods. Our binary models were computed with the MESA stellar evolution code, considering an initial primary star mass of 12.6 M⊙ and secondaries with initial masses of between ∼12 M⊙ and ∼1.3 M⊙, and focus on initial orbital periods above ∼500 d. We compute these models up to the time of iron core collapse in the primary. Results. Our models exhibit varying degrees of stripping due to mass transfer, resulting in SN progenitor models ranging from fully stripped helium stars to stars that have not been stripped at all. We find that Roche lobe overflow often leads to incomplete stripping of the mass donor, resulting in a large variety of pre-SN envelope masses. In many of our models, the red supergiant (RSG) donor stars undergo core collapse during Roche lobe overflow, with mass transfer and therefore system mass-loss rates of up to 0.01 M⊙ yr−1 at that time. The corresponding CSM densities are similar to those inferred for Type IIn SNe, such as SN 1998S. In other cases, the mass transfer becomes unstable, leading to a common-envelope phase at such late time that the mass donor explodes before the common envelope is fully ejected or the system has merged. We argue that this may cause significant pre-SN variability, as witnessed for example in SN 2020tlf. Other models suggest a common-envelope ejection just centuries before core collapse, which may lead to the strongest interactions, as observed in superluminous Type IIn SNe, such as SN 1994W and SN 2006gy. Conclusions. Wide massive binaries exhibit properties that may not only explain the diverse envelope stripping inferred in Type Ib, IIb, IIL, and IIP SNe, but also offer a natural framework to understand a broad range of hydrogen-rich interacting SNe. On the other hand, the flash features observed in many Type IIP SNe, such as SN 2013fs, may indicate that RSG atmospheres are more extended than currently assumed; this could enhance the parameter space for wide binary interaction. [ABSTRACT FROM AUTHOR]

We present photometry and spectroscopy of the slowly evolving superluminous Type IIn supernova (SN) 2015da. SN 2015da is extraordinary for its very high peak luminosity, and also for sustaining a high luminosity for several years. Even at 8 yr after explosion, SN 2015da remains as luminous as the peak of a normal SN II-P. The total radiated energy integrated over this time period (with no bolometric correction) is at least |$1.6 \times 10^{51}$| erg (or 1.6 FOE). Including a mild bolometric correction, adding kinetic energy of the expanding cold dense shell of swept-up circumstellar material (CSM), and accounting for asymmetry, the total explosion kinetic energy was likely 5–10 FOE. Powering the light curve with CSM interaction requires an energetic explosion and 20 M |$_{\odot }$| of H-rich CSM, which in turn implies a massive progenitor system |$\gt $| 30 M |$_{\odot }$|⁠. Narrow P Cyg features show steady CSM expansion at 90 km s |$^{-1}$|⁠ , requiring a high average mass-loss rate of |$\sim$| 0.1 M |$_{\odot }$|  yr |$^{-1}$| sustained for two centuries before explosion (although ramping up toward explosion time). No current theoretical model for single-star pre-SN mass-loss can account for this. The slow CSM, combined with broad wings of H  |$\alpha$| indicating H-rich material in the unshocked ejecta, disfavours a pulsational pair instability model for the pre-SN mass-loss. Instead, violent pre-SN binary interaction is a likely culprit. Finally, SN 2015da exhibits the characteristic asymmetric blueshift in its emission lines from shortly after peak until the present epoch, adding another well-studied superluminous SNe IIn with unambiguous evidence of post-shock dust formation. [ABSTRACT FROM AUTHOR]

We present an alternative model of unusual type-IIP SN 2018gj. Despite the short plateau and early gamma-ray escape seeming to favor low-mass ejecta, our hydrodynamic model requires a large ejected mass (≈23 M ⊙ ). The high ejecta velocity, we find from hydrogen lines in early spectra, is among the crucial constraints on the hydrodynamic model. We recover the wind density that rules out a notable contribution of the circumstellar interaction to the bolometric luminosity. The early radioactive gamma-ray escape is found to be due to the high velocity of 56Ni, whereas the asymmetry of the H α emission is attributed to the asymmetry of the 56Ni ejecta. The available sample of type-IIP supernovae studied hydrodynamically in a uniform way indicates that the asymmetry of the 56Ni ejecta is probably their intrinsic property. Hydrogen lines in the early spectra of SN 2018gi and SN 2020jfo are found to imply a clumpy structure of the outer ejecta. With two already known similar cases of SN 2008in and SN 2012A we speculate that the clumpiness of the outer ejecta is inherent to type-IIP supernovae related to the red supergiant explosion. [ABSTRACT FROM AUTHOR]

The environmental dependence of Type Ia supernova (SN Ia) luminosities is well established, and efforts are being made to find its origin. Previous studies typically use the currently observed status of the host galaxy. However, given the delay time between the birth of the progenitor star and the SN Ia explosion, the currently observed status may differ from the birth environment of the SN Ia progenitor star. In this paper, employing the chemical evolution and accurately determined stellar population properties of 44 early-type host galaxies, we, for the first time, estimate the SN Ia progenitor star birth environment, specifically [Fe/H]Birth and [α/Fe]Birth. We show that [α/Fe]Birth has a |$30.4^{\text{+10.6}}_{-10.1}{{\ \rm per\ cent}}$| wider range than the currently observed [α/Fe]Current, while the range of [Fe/H]Birth is not statistically different (⁠|$17.9^{\text{+26.0}}_{-27.1}{{\ \rm per\ cent}}$|⁠) to that of [Fe/H]Current. The birth and current environments of [Fe/H] and [α/Fe] are sampled from different populations (p -values of the Kolmogorov–Smirnov test <0.01). We find that light-curve fit parameters are insensitive to [Fe/H]Birth (<0.9σ for the non-zero slope), while a linear trend is observed with Hubble residuals (HRs) at the 2.4σ significance level. With [α/Fe]Birth, no linear trends (<1.1σ) are observed. Interestingly, we find that [α/Fe]Birth clearly splits the SN Ia sample into two groups: SN Ia exploded in [α/Fe]Birth-rich or [α/Fe]Birth-poor environments. SNe Ia exploded in different [α/Fe]Birth groups have different weighted-means of light-curve shape parameters: 0.81 ± 0.33 (2.5σ). They are thought to be drawn from different populations (p -value = 0.01). Regarding SN Ia colour and HRs, there is no difference (<1.0σ) in the weighted-means and distribution (p -value > 0.27) of each [α/Fe]Birth group. [ABSTRACT FROM AUTHOR]

Observations of Type Ia supernovae (SNe Ia) reveal diversity, even within assumed subcategories. Here, the composition of the peculiar iPTF16abc (SN 2016bln) is derived by modelling a time series of optical spectra. iPTF16abc's early spectra combine traits of SNe 1999aa and 1991T known for weak Si  ii λ 6355 and prominent Fe  iii features. However, it differs with weak early Fe  iii lines, and persistent C  ii lines post-peak. It also exhibits a weak Ca  ii H&K feature aligning it with SN 1991T, an observation supported by their bolometric light curves. The early attenuation of Fe  iii results from abundance effect. The weakening of the Si  ii λ 6355 line, stems from silicon depletion in the outer shells, a characteristic shared by both SNe 1999aa and 1991T, indicating a common explosion mechanism that terminates nuclear burning at around 12 000 km s−1 unseen in normal events. Beneath a thin layer of intermediate mass elements (IMEs) with a total mass of 0.18 M⊙, extends a 56Ni rich shell totaling 0.76 M⊙ and generating a bolometric luminosity as high as L peak = 1.60 ± 0.1 × 1043 ergs s−1. Inner layers, typical of SNe Ia, hold neutron-rich elements, (54Fe and 58Ni), totaling 0.20 M⊙. Stable iron, exceeding solar abundance, and carbon, coexist in the outermost layers, challenging existing explosion models. The presence of carbon down to v ≈ 9000 km s−1, totalling ∼0.01 M⊙ unprecedented in this class, links iPTF16abc to SN 2003fg-like events. The retention of 91T-like traits in iPTF16abc underscores its importance in understanding the diversity of SNe Ia. [ABSTRACT FROM AUTHOR]

The production of heavy elements is one of the main by-products of the explosive end of massive stars. A long sought goal is finding differentiated patterns in the nucleosynthesis yields, which could permit identifying a number of properties of the explosive core. Among them, the traces of the magnetic field topology are particularly important for extreme supernova (SN) explosions, most likely hosted by magnetorotational effects. We investigate the nucleosynthesis of five state-of-the-art magnetohydrodynamic models with fast rotation that have been previously calculated in full 3D and that involve an accurate neutrino transport (M1). One of the models does not contain any magnetic field and synthesizes elements around the iron group, in agreement with other CC-SNe models in literature. All other models host a strong magnetic field of the same intensity, but with different topology. For the first time, we investigate the nucleosynthesis of MR-SNe models with a quadrupolar magnetic field and a 90° tilted dipole. We obtain a large variety of ejecta compositions reaching from iron nuclei to nuclei up to the third r -process peak. We assess the robustness of our results by considering the impact of different nuclear physics uncertainties such as different nuclear masses, β−-decays and β−-delayed neutron emission probabilities, neutrino reactions, fission, and a feedback of nuclear energy on the temperature. We find that the qualitative results do not change with different nuclear physics input. The properties of the explosion dynamics and the magnetic field configuration are the dominant factors determining the ejecta composition. [ABSTRACT FROM AUTHOR]

Dust absorption is invoked in a number of contexts for hiding a star that has survived some sort of transient event from view. Dust formed in a transient is expanding away from the star and, in spherical models, the mass and energy budgets implied by a high optical depth at late times make such models untenable. Concentrating the dust in a disc or torus can in principle hide a source from an equatorial observer using less mass and so delay this problem. However, using axisymmetric dust radiation transfer models with a range of equatorial dust concentrations, we find that this is quite difficult to achieve in practice. The polar optical depth must be either low or high to avoid scattering optical photons to equatorial observers. Most of the emission remains at wavelengths easily observed by JWST. The equatorial brightness can be significantly suppressed for very discy configurations with little polar optical depth – but only by a factor of ∼2 for polar optical depths of τp = 1 and ∼5 for τp = 0.1 even for a very high optical depth disc (τe = 1000) viewed edge-on. It is particularly difficult to hide a source with silicate dusts because the absorption feature near 10 µm frequently leads to the emission being concentrated just bluewards of the feature, near 8 µm. [ABSTRACT FROM AUTHOR]

We present multi-epoch spectropolarimetry and spectra for a sample of 14 Type IIn supernovae (SNe IIn). We find that after correcting for likely interstellar polarization, SNe IIn commonly show intrinsic continuum polarization of 1–3 per cent at the time of peak optical luminosity, although a few show weaker or negligible polarization. While some SNe IIn have even stronger polarization at early times, their polarization tends to drop smoothly over several hundred days after peak. We find a tendency for the intrinsic polarization to be stronger at bluer wavelengths, especially at early times. While polarization from an electron scattering region is expected to be grey, scattering of SN light by dusty circumstellar material (CSM) may induce such a wavelength-dependent polarization. For most SNe IIn, changes in polarization degree and wavelength dependence are not accompanied by changes in the position angle, requiring that asymmetric pre-SN mass loss had a persistent geometry. While 2–3 per cent polarization is typical, about 30 per cent of SNe IIn have very low or undetected polarization. Under the simplifying assumption that all SN IIn progenitors have axisymmetric CSM (i.e. disc/torus/bipolar), then the distribution of polarization values we observe is consistent with similarly asymmetric CSM seen from a distribution of random viewing angles. This asymmetry has very important implications for understanding the origin of pre-SN mass loss in SNe IIn, suggesting that it was shaped by binary interaction. [ABSTRACT FROM AUTHOR]