Search

We report on radio follow-up observations of the nearby Type II supernova, SN 2023ixf, spanning from 1.7 to 269.9 days after the explosion, conducted using three very long baseline interferometers (VLBIs), which are the Japanese VLBI Network (JVN), the VLBI Exploration of Radio Astrometry (VERA), and the Korean VLBI Network (KVN). In three observation epochs (152.3, 206.1, and 269.9 days), we detected emission at the 6.9 and 8.4 GHz bands, with a flux density of $\sim 5$ mJy. The flux density reached a peak at around 206.1 days, which is longer than the timescale to reach the peak observed in typical Type II supernovae. Based on the analytical model of radio emission, our late-time detections were inferred to be due to the decreasing optical depth. In this case, the mass-loss rate of the progenitor is estimated to have increased from $\sim 10^{-6} - 10^{-5}\, M_{\odot}\,{\rm yr^{-1}}$ to $\sim 10^{-4}\, M_{\odot}\,{\rm yr^{-1}}$ between 28 and 6 years before the explosion. Our radio constraints are also consistent with the mass-loss rate to produce a confined circumstellar medium proposed by previous studies, which suggest that the mass-loss rate increased from $\sim 10^{-4}\, M_{\odot}\,{\rm yr^{-1}}$ to $\gtrsim 10^{-2}\, M_{\odot}\,{\rm yr^{-1}}$ in the last few years before the explosion., Comment: 12 pages, 3 figures, 3 tables. Accepted for publication in ApJ

Type Ia supernovae arise from the thermonuclear explosions of white dwarfs in multiple star systems. A rare sub-class of SNe Ia exhibit signatures of interaction with circumstellar material (CSM), allowing for direct constraints on companion material. While most known events show evidence for dense nearby CSM identified via peak-light spectroscopy (as SNe Ia-CSM), targeted late-time searches have revealed a handful of cases exhibiting delayed CSM interaction with detached shells. Here, we present the first all-sky search for late CSM interaction in SNe Ia using a new image-subtraction pipeline for mid-infrared data from the NEOWISE space telescope. Analyzing a sample of $\approx 8500$ SNe Ia, we report evidence for late-time mid-infrared brightening in six previously overlooked events spanning sub-types SNe Iax, normal SNe Ia, SNe Ia-91T and super-Chandra SNe Ia. Our systematic search doubles the known sample, and suggests that $\gtrsim 0.1$% of SNe Ia exhibit mid-infrared signatures of delayed CSM interaction. The mid-infrared light curves ubiquitously indicate the presence of multiple (or extended) detached CSM shells located at $\gtrsim 10^{16}-10^{17}$ cm, containing $10^{-4}-10^{-2}$ $M_\odot$ of dust, with some sources showing evidence for new dust formation, likely within the cold, dense shell of the ejecta. We do not detect interaction signatures in spectroscopic and radio follow-up; however, the limits are largely consistent with previously confirmed events given the sensitivity and observation phase. Our results highlight that CSM interaction is more prevalent than previously estimated from optical and ultraviolet searches, and that mid-infrared synoptic surveys provide a unique window into this phenomenon., Comment: 16 pages, 5 figures, 1 table. Comments welcome!

Recent high-cadence transient surveys have discovered rapid transients whose light curve timescales are shorter than those of typical supernovae. In this paper, we present a systematic search for rapid transients at medium-high redshifts among 3381 supernova candidates obtained from the Subaru HSC-SSP transient survey. We developed a machine learning classifier to classify the supernova candidates into four types (Type Ia, Ibc, II supernovae, and rapid transients) based on the features derived from the light curves. By applying this classifier to the 3381 supernova candidates and by further applying the quality cut, we selected 14 rapid transient samples. They are located at a wide range of redshifts ($0.34 \leq z \leq 1.85$) and show a wide range of the peak absolute magnitude ($-17 \geq M \geq -22$). The event rate of the rapid transients is estimated to be $\sim 6\times10^3~\rm{events~yr^{-1}~Gpc^{-3}}$ at $z \sim 0.74$, which corresponds to about $2$ $\%$ of the event rate of normal core-collapse supernovae at the similar redshift. Based on the luminosity and color evolution, we selected two candidates of Type Ibn supernovae at $z\sim0.75$. The event rate of Type Ibn SN candidates is more than 1 $\%$ of Type Ib SN rate at the same redshift, suggesting that this fraction of massive stars at this redshift range eruptively ejects their He-rich envelope just before the explosions. Also, two objects at $z=1.37$ and 1.85 show high luminosities comparable to superluminous supernovae. Their event rate is about 10-25 $\%$ of superluminous supernovae at $z\sim 2$., Comment: 17 pages, 13 figures, accepted for publication in ApJ

Interacting binary star systems play a critical role in many areas of astrophysics. One interesting example of a binary merger product are Thorne-\.Zytkow Objects (T\.ZOs), stars that look like red supergiants but contain neutron stars at their cores. T\.ZOs were theorized nearly five decades ago, and significant work has gone into understanding the physics of their formation, evolution, and stability. Several searches for T\.ZO candidates have also been carried out. Whether or not T\.ZOs could even exist or if they would be stable after formation has also been investigated. Understanding the existence and possible prevalence of T\.ZOs would have important effects on our understanding of binary evolution, stellar mergers, and inform binary population synthesis models. In this chapter, we review the formation channels, evolution and structure, final fates, and observable signatures of T\.ZOs, as well as candidates in the literature, from the inception of T\.ZO theory to recent progress in the field., Comment: This is a pre-print of a chapter for the Encyclopedia of Astrophysics (edited by I. Mandel, section editor J. Andrews) to be published by Elsevier as a Reference Module. 17 pages, 10 figures

Recent observations of core-collapse supernovae revealed that the existence of dense circumstellar matter (CSM) around their progenitors is ubiquitous. Interaction of supernova ejecta with such a dense CSM is a potential production sight of high-energy cosmic rays (CRs), gamma-rays, and neutrinos. We estimate the gamma-ray and neutrino signals from SN 2023ixf, a core-collapse supernova occurred in a nearby galaxy M101, which exhibits signatures of the interaction with the confined dense CSM. Using radiation-hydrodynamic simulation model calibrated by the optical and ultraviolet observations of SN 2023ixf, we find that the CRs cannot be accelerated in the early phase because the sharp velocity jump at the shock disappears due to strong radiation pressure. Roughly 4 days after the explosion, the collisionless sub-shock is formed in the CSM, which enables the CR production and leads to gamma-ray and neutrino emissions. The shock sweeps up the entire dense CSM roughly 9 days after the explosion, which ceases the high-energy radiation. Based on this scenario, we calculate the gamma-ray and neutrino signals, which have a peak around 9 days after the explosion. We can constrain the cosmic-ray production efficiency to be less than 30\% by comparing our prediction to the Fermi-LAT data. Future multi-messenger observations with an enlarged sample of nearby supernovae will provide a better constraint on the cosmic-ray production efficiency in the early phases of supernovae., Comment: 9 pages, 5 figures

SN 2023ixf is one of the most well-observed core-collapse supernova in recent decades, yet there is inconsistency in the inferred zero-age-main-sequence (ZAMS) mass $M_{\rm ZAMS}$ of its progenitor. Direct observations of the pre-SN red supergiant (RSG) estimate $M_{\rm ZAMS}$ spanning widely from 11 to 18 $M_{\rm \odot}$. Additional constraints, including host environment and the pulsation of its progenitor RSG, suggest a massive progenitor with $M_{\rm ZAMS}$ > 17 $M_{\rm \odot}$. However, the analysis of the properties of supernova, from light curve modeling to late phase spectroscopy, favor a relatively low mass scenario ($M_{\rm ZAMS}$ < 15 $M_{\rm \odot}$). In this work, we conduct systematic analysis of SN 2023ixf, from the RSG progenitor, plateau phase light curve to late phase spectroscopy. Using MESA+STELLA to simulate the RSG progenitor and their explosions, we find that, despite the ZAMS mass of the RSG models being varied from 12.0 to 17.5 $M_{\rm \odot}$, they can produce light curves that well match with SN 2023ixf if the envelope mass and the explosion energy are allowed to vary. Using late phase spectroscopy as independent measurement, the oxygen emission line [O I] suggests the ZAMS is intermediate massive (~16.0 $M_{\rm \odot}$), and the relatively weak H$\alpha$ emission line indicates the hydrogen envelope has been partially removed before the explosion. By incorporating the velocity structure derived from the light curve modeling into an axisymmetric model, we successfully generated [O I] line profiles that are consistent with the [O I] line observed in late phase spectroscopy of SN 2023ixf. Bringing these analyses together, we conclude that SN 2023ixf is the aspherical explosion of an intermediate massive star ($M_{\rm ZAMS}$ = 15-16 $M_{\rm \odot}$) with the hydrogen envelope being partially stripped to 4-5 $M_{\rm \odot}$ prior to its explosion., Comment: 12 pages, 4 figures. Submitted

We present the long-term photometric and spectroscopic analysis of a transitioning SN~IIn/Ibn from $-$10.8 d to 150.7 d post $V$-band maximum. SN~2021foa shows prominent He {\sc i} lines comparable in strength to the H$\alpha$ line around peak, placing SN~2021foa between the SN~IIn and SN~Ibn populations. The spectral comparison shows that it resembles the SN~IIn population at pre-maximum, becomes intermediate between SNe~IIn/Ibn and at post-maximum matches with SN~IIn 1996al. The photometric evolution shows a precursor at $-$50 d and a light curve shoulder around 17d. The peak luminosity and color evolution of SN 2021foa are consistent with most SNe~IIn and Ibn in our comparison sample. SN~2021foa shows the unique case of a SN~IIn where the narrow P-Cygni in H$\alpha$ appear at later stages. The H$\alpha$ profile consists of a narrow (500 -- 1200 km s$^{-1}$) component, intermediate width (3000 -- 8000 km s$^{-1}$) and broad component in absorption. Temporal evolution of the H$\alpha$ profile favours a disk-like CSM geometry. Hydrodynamical modelling of the lightcurve well reproduces a two-component CSM structure with different densities ($\rho$ $\propto$ r$^{-2}$ -- $\rho$ $\propto$ r$^{-5}$), mass-loss rates (10$^{-3}$ -- 10$^{-1}$ M$_{\odot}$ yr$^{-1}$) assuming a wind velocity of 1000 km s$^{-1}$ and having a CSM mass of 0.18 M$_{\odot}$. The overall evolution indicates that SN~2021foa most likely originated from a LBV star transitioning to a WR star with the mass-loss rate increasing in the period from 5 to 0.5 years before the explosion or it could be due to a binary interaction., Comment: Submitted to MNRAS; 20 pages, 16 figures, 4 tables

Pair-instability supernovae (PISNe) have long been predicted to be the final fates of near-zero-metallicity very massive stars ($Z < Z_\odot/3$, $\mathrm{M}_\mathrm{ZAMS} \gtrsim 140 \mathrm{M}_\odot$). However, no definite PISN has been observed to date, leaving theoretical modelling validation open. To investigate the observability of these explosive transients, we combine detailed stellar evolution models for PISNe formation, computed from the Binary Population and Spectral Synthesis code suite, BPASS, with the star formation history of all individual computational elements in the Illustris-TNG simulation. This allows us to compute comic PISN rates and predict their host galaxy properties. Of particular importance is that IllustrisTNG galaxies do not have uniform metallicities throughout, with metal-enriched galaxies often harbouring metal-poor pockets of gas where PISN progenitors may form. Accounting for the chemical inhomogeneities within these galaxies, we find that the peak redshift of PISNe formation is $z=3.5$ instead of the value of $z=6$ when ignoring chemical inhomogeneities within galaxies. Furthermore, the rate increases by an order of magnitude from 1.9 to 29 PISN Gpc$^{-3}$ yr$^{-1}$ at $z=0$, if the chemical inhomogeneities are considered. Using state-of-the-art theoretical PISN light curves, we find an observed rate of $13.8$ (1.2) visible PISNe per year for the Euclid-Deep survey, or $83$ (7.3) over the six-year lifetime of the mission when considering chemically inhomogeneous (homogenous) systems. Interestingly, only 12 per cent of helium PISN progenitors are sufficiently massive to power a super-luminous supernova event, which can potentially explain why PISN identification in time-domain surveys remains elusive and progress requires dedicated strategies., Comment: Accepted 2024 August 23 in Monthly Notices of the Royal Astronomical Society

Superluminous supernovae (SLSNe) are a population of supernovae (SNe) whose peak luminosities are much larger than those of canonical SNe. Although SLSNe were simply defined by their peak luminosity at first, it is currently recognized that they show rich spectroscopic diversities including hydrogen-poor (Type I) and hydrogen-rich (Type II) subtypes. The exact mechanisms making SLSNe luminous are still not fully understood, but there are mainly four major suggested luminosity sources (radioactive decay of 56Ni, circumstellar interaction, magnetar spin-down, and fallback accretion). We provide an overview of observational properties of SLSNe and major theoretical models for them. Future transient surveys are expected to discover SLSNe at high redshifts which will provide a critical information in revealing their nature., Comment: This is a pre-print of a chapter for the Encyclopedia of Astrophysics (edited by I. Mandel, section editor F.R.N. Schneider) to be published by Elsevier as a Reference Module

The current perspective about the explosions of massive hydrogen-rich blue supergiants is that they resemble SN 1987A. These so-called peculiar Type II supernovae, however, are one of the rarest types of supernovae and may not hence be the fate of all blue supergiants. In this work, we explore other explosion scenarios for blue supergiants. We create synthetic light curves from the explosions of blue supergiant models born from binary mergers, over a range of explosion energies and 56Ni masses. We find that blue supergiant explosions may also lead to intermediate-luminosity red transients. We thus identify two categories of supernovae possible from blue supergiant explosions: those with high 56Ni masses (> ~ 0.01 Msun) result in slow-rising, dome-shaped light curves like SN 1987A. Lower 56Ni masses result in low-luminosity, short-plateau light curves resembling some intermediate-luminosity red transients and Type II supernovae like SN 2008bp, which are possible from the explosions of compact blue supergiants and not from the far more extended red supergiants. Our results indicate that blue supergiant explosions are more diverse than SN 1987A-like events and may be hidden among different kinds of transients, explaining the possible discrepancies between the expected fraction of blue supergiants born from binary mergers and the observed fraction of SN 1987A-like supernovae., Comment: 8 pages, 5 figures, 1 table, accepted by Publications of the Astronomical Society of Japan

Interaction between supernova (SN) ejecta and dense circumstellar medium (CSM) with a flat density structure ($\rho \propto r^{-s}, s < 1.5$) was recently proposed as a possible mechanism behind interacting SNe that exhibit exceptionally long rise times exceeding 100 days. In such a configuration, the interaction luminosity keeps rising until the reverse shock propagates into the inner layers of the SN ejecta. We investigate the light curves of SNe interacting with a flatly distributed CSM in detail, incorporating the effects of photon diffusion inside the CSM into the model. We show that three physical processes - the shock breakout, the propagation of the reverse shock into the inner ejecta, and the departure of the shock from the dense CSM - predominantly determine the qualitative behaviour of the light curves. Based on the presence and precedence of these processes, the light curves of SNe interacting with flatly distributed CSM can be classified into five distinct morphological classes. We also show that our model can qualitatively reproduce doubly peaked SNe whose peaks are a few tens of days apart, such as SN 2005bf and SN 2022xxf. Our results show that the density distribution of the CSM is an important property of CSM that contributes to the diversity in light curves of interacting SNe., Comment: 17 pages, 7 figures, 2 tables, accepted for publication in the Astrophysical Journal

We present the discovery and early observations of the nearby Type II supernova (SN) 2024ggi in NGC 3621 at 6.64 +/- 0.3 Mpc. The SN was caught 5.8 (+1.9 -2.9) hours after its explosion by the ATLAS survey. Early-phase, high-cadence, and multi-band photometric follow-up was performed by the Kinder (Kilonova Finder) project, collecting over 1000 photometric data points within a week. The combined o- and r-band light curves show a rapid rise of 3.3 magnitudes in 13.7 hours, much faster than SN 2023ixf (another recent, nearby, and well-observed SN II). Between 13.8 and 18.8 hours after explosion SN 2024ggi became bluer, with u-g colour dropping from 0.53 to 0.15 mag. The rapid blueward evolution indicates a wind shock breakout (SBO) scenario. No hour-long brightening expected for the SBO from a bare stellar surface was detected during our observations. The classification spectrum, taken 17 hours after the SN explosion, shows flash features of high-ionization species such as Balmer lines, He I, C III, and N III. Detailed light curve modeling reveals critical insights into the properties of the circumstellar material (CSM). Our favoured model has an explosion energy of 2 x 10^51 erg, a mass-loss rate of 10^-3 solar_mass/yr (with an assumed 10 km/s wind), and a confined CSM radius of 6 x 10^14 cm. The corresponding CSM mass is 0.4 solar_mass. Comparisons with SN 2023ixf highlight that SN 2024ggi has a smaller CSM density, resulting in a faster rise and fainter UV flux. The extensive dataset and the involvement of citizen astronomers underscore that a collaborative network is essential for SBO searches, leading to more precise and comprehensive SN characterizations., Comment: 11 pages, 5 figures in manuscript, 6 pages in appendix, submitted to ApJL

The JWST Advanced Deep Extragalactic Survey (JADES) is a multi-cycle JWST program that has taken among the deepest near-/mid-infrared images to date (down to $\sim$30 ABmag) over $\sim$25 arcmin$^2$ in the GOODS-S field in two sets of observations with one year of separation. This presented the first opportunity to systematically search for transients, mostly supernovae (SNe), out to $z$$>$2. We found 79 SNe: 38 at $z$$<$2, 23 at 2$<$$z$$<$3, 8 at 3$<$$z$$<$4, 7 at 4$<$$z$$<$5, and 3 with undetermined redshifts, where the redshifts are predominantly based on spectroscopic or highly reliable JADES photometric redshifts of the host galaxies. At this depth, the detection rate is $\sim$1-2 per arcmin$^2$ per year, demonstrating the power of JWST as a supernova discovery machine. We also conducted multi-band follow-up NIRCam observations of a subset of the SNe to better constrain their light curves and classify their types. Here, we present the survey, sample, search parameters, spectral energy distributions (SEDs), light curves, and classifications. Even at $z$$\geq$2, the NIRCam data quality is high enough to allow SN classification via multi-epoch light-curve fitting with confidence. The multi-epoch SN sample includes a Type Ia SN at $z_{\mathrm{spec}}$$=$2.90, Type IIP SN at $z_{\mathrm{spec}}$$=$3.61, and a Type Ic-BL SN at $z_{\mathrm{spec}}$$=$2.845. We also found that two $z$$\sim$16 galaxy candidates from the first imaging epoch were actually transients that faded in the second epoch, illustrating the possibility that moderate/high-redshift SNe could mimic high-redshift dropout galaxies., Comment: 43 pages, 15 figures, 15 tables. Submitted to ApJ. Appendix A (64 MB) is available at https://drive.google.com/file/d/1xs5jXUVOvdDPgdghK72KR1FMGvPcK7dv/view?usp=sharing . Appendix B (81 MB) is available at https://drive.google.com/file/d/18ImLT80pQdPzXCZA-KEy21DaE2CQiGz1/view?usp=sharing . References updated, typos fixed, minor restructuring

We estimate the progenitor and explosion properties of the nearby Type II SN 2023ixf using a synthetic model grid of Type II supernova light curves. By comparing the light curves of SN 2023ixf with the pre-existing grid of Type II supernovae containing about 228,000 models with different combinations of the progenitor and explosion properties, we obtain the chi2 value for every model and evaluate the properties of the models providing small values of chi2. We found that the light-curve models with the progenitor zero-age main-sequence mass of 10 Msun, the explosion energy of (2-3)e51 erg, the 56Ni mass of 0.04-0.06 Msun, the mass-loss rate of 1e-3 - 1e-2 Msun/yr with a wind velocity of 10 km/s, and the dense, confined circumstellar matter radius of (6-10)e14 cm match well to the observed light curves of SN 2023ixf. The photospheric velocity evolution of these models is also consistent with the observed velocity evolution. We note that the progenitor mass estimate could be affected by the adopted progenitor models. Although our parameter estimation is based on a pre-existing model grid and we do not perform any additional computations, the estimated parameters are consistent with those obtained by the detailed modeling of SN 2023ixf previously reported. This result shows that comparing the pre-existing model grid is a reasonable way to obtain a rough estimate for the properties of Type II supernovae. This simple way to estimate the properties of Type II supernovae will be essential in the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) era when thousands of Type II supernovae are expected to be discovered yearly., Comment: 10 pages, 7 figures, accepted by Publications of the Astronomical Society of Japan

We report an analysis of a sample of 186 spectroscopically confirmed Type II supernova (SN) light curves (LCs) obtained from a combination of Zwicky Transient Facility (ZTF) and Asteroid Terrestrial-impact Last Alert System (ATLAS) observations. We implement a method to infer physical parameters from these LCs using hydrodynamic models that take into account the progenitor mass, the explosion energy, and the presence of circumstellar matter (CSM). The CSM is modelled via the mass loss rate, wind acceleration at the surface of the progenitor star with a $\beta$ velocity law, and the CSM radius. We also infer the time of explosion, attenuation (A$_V$), and the redshift for each SN. Our results favor low-mass progenitor stars (M$_{ZAMS}$\,$<$14\,$M_\odot$) with a dense CSM ($\dot{M}$ $>$ 10$^{-3}$ [M$_\odot$ yr$^{-1}$], a CSM radius of $\sim$ 10$^{15}$ cm, and $\beta$ $>$ 2). Additionally, we find that the redshift inferred from the supernova LCs is significantly more accurate than that inferred using the host galaxy photometric redshift, suggesting that this method could be used to infer more accurate host galaxy redshifts from large samples of SNe II in the LSST era. Lastly, we compare our results with similar works from the literature.

Core-collapse supernovae (CCSNe) are the terminal explosions of massive stars. While most massive stars explode as iron-core-collapse supernovae (FeCCSNe), slightly less massive stars explode as electron-capture supernovae (ECSNe), shaping the low-mass end of CCSNe. ECSNe was proposed $\sim 40$ years ago and first-principles simulations also predict their successful explosions. Observational identification and investigation of ECSNe are important for the completion of stellar evolution theory. To date, only one promising candidate has been proposed, SN 2018zd, other than the historical progenitor of the Crab Nebula, SN 1054. We present representative synthetic light curves of low-mass FeCCSNe and ECSNe exploding with energies in circumstellar media (CSM) estimated with theoretically or observationally plausible methods. The plateaus of the ECSNe are shorter, brighter, and bluer than those of the FeCCSNe. To investigate the robustness of their intrinsic differences, we adopted various explosion energies and CSM. Although they may have similar bolometric light-curve plateaus, ECSNe are bluer than FeCCSNe in the absence of strong CSM interaction, illustrating that multicolor observations are essential to identify ECSNe. This provides a robust indicator of ECSNe because the bluer plateaus stem from the low-density envelopes of their super-asymptotic-giant-branch progenitors. Furthermore, we propose a distance-independent method to identify ECSNe: $(g-r)_{t_{\rm PT}/2} < 0.008 \times t_{\rm PT} - 0.4$, i.e., blue $g-r$ at the middle of the plateau $(g-r)_{t_{\rm PT}/2}$, where $t_{\rm PT}$ is the transition epoch from plateau to tail. Using this method, we identified SN 2018zd as an ECSN, which we believe to be the first ECSN identified with modern observing techniques., Comment: 16 pages, 14 figures, 4 tables, submitted to Astrophysical Journal

Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and long (~100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew's spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak's photospheric phase, before H$\alpha$ resurfaces again during the nebular phase. The nebular lines ([O I], [Ca II], Mg I], H$\alpha$) exhibit a double-peaked structure which hints towards a clumpy or non-spherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than normal SESNe as well as requiring a very high $^{56}$Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak which also powers the second peak of the light curve through SN-CSM interaction., Comment: 22 pages, 11 figures, 5 tables

We present photometric, spectroscopic and polarimetric observations of the intermediate-luminosity Type IIP supernova (SN) 2021gmj from 1 to 386 days after the explosion. The peak absolute V-band magnitude of SN 2021gmj is -15.5 mag, which is fainter than that of normal Type IIP SNe. The spectral evolution of SN 2021gmj resembles that of other sub-luminous supernovae: the optical spectra show narrow P-Cygni profiles, indicating a low expansion velocity. We estimate the progenitor mass to be about 12 Msun from the nebular spectrum and the 56Ni mass to be about 0.02 Msun from the bolometric light curve. We also derive the explosion energy to be about 3 x 10^{50} erg by comparing numerical light curve models with the observed light curves. Polarization in the plateau phase is not very large, suggesting nearly spherical outer envelope. The early photometric observations capture the rapid rise of the light curve, which is likely due to the interaction with a circumstellar material (CSM). The broad emission feature formed by highly-ionized lines on top of a blue continuum in the earliest spectrum gives further indication of the CSM at the vicinity of the progenitor. Our work suggests that a relatively low-mass progenitor of an intermediate-luminosity Type IIP SN can also experience an enhanced mass loss just before the explosion, as suggested for normal Type IIP SNe., Comment: 18 pages, 16 figures, resubmitted to MNRAS after addressing referee comments

Thanks to the mid-infrared capability offered by Spitzer, systematic searches of dust in SNe have been carried out over the past decade. Studies have revealed the presence of a substantial amount of dust over a broad range of SN subtypes. How normal SNe present mid-IR excess at later time and turn out to be dusty SNe can be affected by several factors, such as mass-loss history and envelope structure of progenitors and their explosion environment. All these can be combined and related to their environmental properties. A systematic analysis of SNe that exploded under a dusty environment could be of critical importance to measure the properties of the dust-veiled exploding stars, and whether such an intense dust production process is associated with the local environment. In this work, we firstly use the IFS data to study the environmental properties of dusty SNe compared to those of normal ones, and analyze correlations between the environmental properties and their dust parameters. We find that dusty SNe have a larger proportion located at higher SFR regions compared to the normal types. The occurrence of dusty SNe is less dependent on metallicity, with the oxygen abundance spanning from subsolar to oversolar metallicity. We also find the host extinction of dusty SNe scatters a lot, with about 40% of dusty SN located at extremely low extinction environments, and another 30% of them with considerably high host extinction of E(B-V)>0.6 mag., Comment: 20 pages, 15 figures, submitted to ApJ

General relativistic instability supernovae at ~10 < z < ~15 are predicted to be observed as red faint point sources, and they can be detected only in the reddest filters in JWST/NIRCam (F444W and F356W). They should be observed as persistent point sources with little flux variations for a couple of decades because of time dilation. We search for static point sources detected only in the F444W filter or only in the F444W and F356W filters in the early JWST deep field data. No real point source of such kind is identified. Therefore, the general relativistic instability supernova rate at ~10 < z < ~15 is constrained to be less than ~ 8e-7 Mpc-3 yr-1 for the first time., Comment: 4 pages, 2 figures, 1 table, accepted by Monthly Notices of the Royal Astronomical Society

We present optical, near-infrared, and radio observations of supernova (SN) SN~IIb 2022crv. We show that it retained a very thin H envelope and transitioned from a SN~IIb to a SN~Ib; prominent H$\alpha$ seen in the pre-maximum phase diminishes toward the post-maximum phase, while He {\sc i} lines show increasing strength. \texttt{SYNAPPS} modeling of the early spectra of SN~2022crv suggests that the absorption feature at 6200\,\AA\ is explained by a substantial contribution of H$\alpha$ together with Si {\sc ii}, as is also supported by the velocity evolution of H$\alpha$. The light-curve evolution is consistent with the canonical stripped-envelope supernova subclass but among the slowest. The light curve lacks the initial cooling phase and shows a bright main peak (peak M$_{V}$=$-$17.82$\pm$0.17 mag), mostly driven by radioactive decay of $\rm^{56}$Ni. The light-curve analysis suggests a thin outer H envelope ($M_{\rm env} \sim$0.05 M$_{\odot}$) and a compact progenitor (R$_{\rm env}$ $\sim$3 R$_{\odot}$). An interaction-powered synchrotron self-absorption (SSA) model can reproduce the radio light curves with a mean shock velocity of 0.1c. The mass-loss rate is estimated to be in the range of (1.9$-$2.8) $\times$ 10$^{-5}$ M$_{\odot}$ yr$^{-1}$ for an assumed wind velocity of 1000 km s$^{-1}$, which is on the high end in comparison with other compact SNe~IIb/Ib. SN~2022crv fills a previously unoccupied parameter space of a very compact progenitor, representing a beautiful continuity between the compact and extended progenitor scenario of SNe~IIb/Ib., Comment: The paper contains 20 figures and 9 tables. Accepted for publication in The Astrophysical Journal (ApJ)

We present a comprehensive data set of supernova (SN) 2016adj located within the central dust lane of Centaurus A. SN 2016adj is significantly reddened and after correcting the peak apparent $B$-band magnitude ($m_B = 17.48\pm0.05$) for Milky Way reddening and our inferred host-galaxy reddening parameters (i.e., $R_{V}^{host} = 5.7\pm0.7$ and $A_{V}^{host} = 6.3\pm0.2$), we estimate it reached a peak absolute magnitude of $M_B \sim -18$. Detailed inspection of the optical/NIR spectroscopic time-series reveals a carbon-rich SN Ic and not a SN Ib/IIb as previously suggested in the literature. The NIR spectra shows prevalent carbon-monoxide formation occurring already by +41 days past $B$-band maximum, which is $\approx 11$ days earlier than previously reported in the literature for this object. Interestingly around two months past maximum, the NIR spectrum of SN~2016adj begins to exhibit H features, with a +97~d medium resolution spectrum revealing both Paschen and Bracket lines with absorption minima of $\sim 2000$ km/s, full-width-half-maximum emission velocities of $\sim 1000$ km/s, and emission line ratios consistent with a dense emission region. We speculate these attributes are due to circumstellar interaction (CSI) between the rapidly expanding SN ejecta and a H-rich shell of material formed during the pre-SN phase. A bolometric light curve is constructed and a semi-analytical model fit suggests the supernova synthesized 0.5 solar masses of $^{56}$Ni and ejected 4.2 solar masses of material, though these values should be approached with caution given the large uncertainties associated with the adopted reddening parameters, possible CSI contamination, and known light echo emission. Finally, inspection of Hubble Space Telescope archival data yielded no progenitor detection., Comment: Submitted to A&A, comments are welcome

We present optical photometric and spectroscopic observations of the 02es-like type Ia supernova (SN) 2022ywc. The transient occurred in the outskirts of an elliptical host galaxy and showed a striking double-peaked light curve with an early excess feature detected in the ATLAS orange and cyan bands. The early excess is remarkably luminous with an absolute magnitude $\sim -19$, comparable in luminosity to the subsequent radioactively-driven second peak. The spectra resemble the hybrid 02es-like SN 2016jhr, that is considered to be a helium shell detonation candidate. We investigate different physical mechanisms that could power such a prominent early excess and rule out massive helium shell detonation, surface $^{56}$Ni distribution and ejecta-companion interaction. We conclude that SN ejecta interacting with circumstellar material (CSM) is the most viable scenario. Semi-analytical modelling with MOSFiT indicates that SN ejecta interacting with $\sim 0.05\,$M$_{\odot}$ of CSM at a distance of $\sim 10^{14}$ cm can explain the extraordinary light curve. A double-degenerate scenario may explain the origin of the CSM, either by tidally-stripped material from the secondary white dwarf, or disk-originated matter launched along polar axes following the disruption and accretion of the secondary white dwarf. A non-spherical CSM configuration could suggest that a small fraction of 02es-like events viewed along a favourable line of sight may be expected to display a very conspicuous early excess like SN 2022ywc., Comment: Accepted to ApJL after minor revision

Some interaction-powered supernovae have long rise times of more than 100 days. We show that such long rise times are naturally expected if circumstellar matters (CSM) have a flat density structure (s <~ 1.5, where rho_CSM ~ r^{-s}). In such cases, bolometric luminosities from the CSM interaction keep increasing as long as the CSM interacts with the outer layers of the SN ejecta. Thus, the rise time is determined by the dynamical timescale in which the reverse shock propagates the outer layers of the SN ejecta, not by the timescales in which photons diffuse in the CSM as often considered. Interaction-powered supernovae with very long rise times can be an important probe of extensive non-steady mass loss in massive stars., Comment: 5 pages, 5 figures, 1 table, accepted by Monthly Notices of the Royal Astronomical Society

Binary neutron star mergers and neutron star-black hole mergers are multi-messenger sources that can be detected in gravitational waves and in electromagnetic radiation. The low electron fraction of neutron star merger ejecta favors the production of heavy elements such as lanthanides and actinides via rapid neutron capture (r-process). The decay of these unstable nuclei powers an infrared-bright transient called a "kilonova". The discovery of a population of kilonovae will allow us to determine if neutron star mergers are the dominant sites for r-process element nucleosynthesis, constrain the equation of state of nuclear matter, and make independent measurements of the Hubble constant. The Nancy Grace Roman Space Telescope (Roman) will have a unique combination of depth, near-infrared sensitivity, and wide field of view. These characteristics will enable Roman's discovery of GW counterparts that will be missed by optical telescopes, such as kilonova that are associated with large distances, high lanthanide fractions, high binary mass-ratios, large dust extinction in the line of sight, or that are observed from equatorial viewing angles. Our analysis suggests to (i) make available a rapid (about 1 week) Target of Opportunity mode for GW follow-up; (ii) include observations of the High Latitude Time-Domain survey footprint in at least two filters (preferably the F158 and F213 filters) with a cadence of < 8 days; (iii) operate in synergy with Rubin Observatory. Following these recommendations, we expect that 1-6 kilonovae can be identified by Roman via ToO observations of well localized (A < 10 sq. deg., 90% C.I.) neutron star mergers during 1.5 years of the LIGO-Virgo-KAGRA fifth (or about 4-21 in during the sixth) observing run. A sample of 5-40 serendipitously discovered kilonovae can be collected in a 5-year high latitude survey., Comment: accepted for publication in Astroparticle Physics

Type IIn supernovae occur when stellar explosions are surrounded by dense hydrogen-rich circumstellar matter. The dense circumstellar matter is likely formed by extreme mass loss from their progenitors shortly before they explode. The nature of Type IIn supernova progenitors and the mass-loss mechanism forming the dense circumstellar matter are still unknown. In this work, we investigate if there are any correlations between Type IIn supernova properties and their local environments. We use Type IIn supernovae with well-observed light-curves and host-galaxy integral field spectroscopic data so that we can estimate both supernova and environmental properties. We find that Type IIn supernovae with a higher peak luminosity tend to occur in environments with lower metallicity and/or younger stellar populations. The circumstellar matter density around Type IIn supernovae is not significantly correlated with metallicity, so the mass-loss mechanism forming the dense circumstellar matter around Type IIn supernovae might be insensitive to metallicity., Comment: 15 pages, 8 pages, 4 tables, accepted by Astronomy & Astrophysics

Nuclear reactions are key to our understanding of stellar evolution, particularly the $^{12}\rm{C}\left({\alpha},{\gamma}\right)^{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 $\texttt{MESA}$ stellar evolution code to explore the impact of uncertainty in the $^{12}\rm{C}\left({\alpha},{\gamma}\right)^{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 ($^{56}{\rm Ni}$) and the explosion energy increase with $^{12}\rm{C}\left({\alpha},{\gamma}\right)^{16}\!\rm{O}$ rate for the same initial mass, except in the high-mass edge region. With a high (about twice the $\texttt{STARLIB}$ standard value) rate, the maximum amount of nickel produced falls below 70 $M_\odot$, while with a low rate (about half of the standard value) it increases up to 83.9 $M_\odot$. 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}\left({\alpha},{\gamma}\right)^{16}\!\rm{O}$ reaction rate., Comment: 17 pages, 12 figures

The effects of the interaction between Type Ia supernova ejecta and their circumstellar wind on the photometric properties of Type Ia supernovae are investigated. We assume that a hydrogen-rich, dense, and extended circumstellar matter (CSM) is formed by the steady mass loss of their progenitor systems. The CSM density is assumed to be proportional to r^{-2}. When the mass-loss rate is above 1e-4 Msun/yr with a wind velocity of 100 km/s, CSM interaction results in an early flux excess in optical light-curves within 4 days of explosion. In these cases, the optical colour quickly evolves to the blue. The ultraviolet flux below 3000 A is found to have a persistent flux excess compared to Type Ia supernovae as long as CSM interaction continues. Type Ia supernovae with progenitor mass-loss rates between 1e-4 and 1e-3 Msun/yr may not have a CSM that is dense enough to affect spectra to make them Type Ia-CSM, but they may still result in Type Ia supernovae with an early optical flux excess. Because they have a persistent ultraviolet flux excess, ultraviolet light curves around the luminosity peak would be significantly different from those with a low-density CSM., Comment: 8 pages, 8 figures, accepted by Monthly Notices of the Royal Astronomical Society

We present ionizing spectra estimated at 13.6--100 eV for ten dwarf galaxies with strong high ionization lines of He {\sc {ii}}$\lambda$4686 and [Ne {\sc{v}}]$\lambda$3426 ([Ne {\sc{iv}}]$\lambda$2424) at $z=0$ ($z=8$) that are identified in our Keck/LRIS spectroscopy and the literature (the JWST ERO program). With the flux ratios of these high ionization lines and $>10$ low-ionization lines of hydrogen, helium, oxygen, neon, and sulfur, we determine ionizing spectra consisting of stellar and non-thermal power-law radiation by photoionization modeling with free parameters of nebular properties including metallicity and ionization parameter, cancelling out abundance ratio differences. We find that all of the observed flux ratios are well reproduced by the photoinization models with the power law index $\alpha_{\rm EUV}$ of $\alpha_{\rm EUV}\sim (-1)-0$ and the luminosity $L_{\rm EUV}$ of $L_{\rm EUV}\sim 10^{40}-10^{42}$ erg s$^{-1}$ at $\sim 55-100$ eV for six galaxies, while four galaxies include large systematics in $\alpha_{\rm EUV}$ caused by stellar radiation contamination. We then compare $\alpha_{\rm EUV}$ and $L_{\rm EUV}$ of these six galaxies with those predicted by the black hole (BH) accretion disk models, and find that $\alpha_{\rm EUV}$ and $L_{\rm EUV}$ are similar to those of the intermediate mass black holes (IMBHs) in BH accretion disk models {albeit with possibilities of the other scenarios.} Confirming these results with a known IMBH having a mass $M_{\rm BH}$ of $M_{\rm BH}=10^{5.75} \ M_\odot$, we find that four local galaxies and one $z=7.665$ galaxy have ionizing spectra consistent with those of IMBHs with $M_{\rm BH} \sim 10^3-10^5 \ M_\odot$., Comment: Accepted for publication in ApJ and 25 pages, 12 figures

We present gas-phase elemental abundance ratios of 7 local extremely metal-poor galaxies (EMPGs) including our new Keck/LRIS spectroscopy determinations together with 33 JWST $z\sim 4-10$ star-forming galaxies in the literature, and compare chemical evolution models. We develop chemical evolution models with the yields of core-collapse supernovae (CCSNe), Type Ia supernovae, hypernovae (HNe), and pair-instability supernovae (PISNe), and compare the EMPGs and high-$z$ galaxies in conjunction with dust depletion contributions. We find that high Fe/O values of EMPGs can (cannot) be explained by PISN metal enrichments (CCSN/HN enrichments even with the mixing-and-fallback mechanism enhancing iron abundance), while that the observed Ar/O and S/O values are much smaller than the predictions of the PISN models. The abundance ratios of the EMPGs can be explained by the combination of Type Ia SNe and CCSNe/HNe whose inner layers of argon and sulfur mostly fallback, which are comparable with Sculptor stellar chemical abundance distribution, suggesting that early chemical enrichment is taken place in the EMPGs. Comparing our chemical evolution models with the star-forming galaxies at $z\sim 4-10$, we find that the Ar/O and S/O ratios of the high-$z$ galaxies are comparable with those of the CCSNe/HNe models, while majority of the high-$z$ galaxies do not have constraints good enough to rule out contributions from PISNe. The high N/O ratio recently reported in GN-z11 cannot be explained even by rotating PISNe, but could be reproduced by the winds of rotating Wolf Rayet stars that end up as a direct collapse.

We present demography of the dynamics and gas-mass fraction of 33 extremely metal-poor galaxies (EMPGs) with metallicities of $0.015-0.195~Z_\odot$ and low stellar masses of $10^4-10^8~M_\odot$ in the local universe. We conduct deep optical integral-field spectroscopy (IFS) for the low-mass EMPGs with the medium high resolution ($R=7500$) grism of the 8m-Subaru FOCAS IFU instrument by the EMPRESS 3D survey, and investigate H$\alpha$ emission of the EMPGs. Exploiting the resolution high enough for the low-mass galaxies, we derive gas dynamics with the H$\alpha$ lines by the fitting of 3-dimensional disk models. We obtain an average maximum rotation velocity ($v_\mathrm{rot}$) of $15\pm3~\mathrm{km~s^{-1}}$ and an average intrinsic velocity dispersion ($\sigma_0$) of $27\pm10~\mathrm{km~s^{-1}}$ for 15 spatially resolved EMPGs out of the 33 EMPGs, and find that all of the 15 EMPGs have $v_\mathrm{rot}/\sigma_0<1$ suggesting dispersion dominated systems. There is a clear decreasing trend of $v_\mathrm{rot}/\sigma_0$ with the decreasing stellar mass and metallicity. We derive the gas mass fraction ($f_\mathrm{gas}$) for all of the 33 EMPGs, and find no clear dependence on stellar mass and metallicity. These $v_\mathrm{rot}/\sigma_0$ and $f_\mathrm{gas}$ trends should be compared with young high-$z$ galaxies observed by the forthcoming JWST IFS programs to understand the physical origins of the EMPGs in the local universe., Comment: 18 pages, 9 figures, accepted for publication in ApJ

A new model grid containing 228,016 synthetic red supergiant explosions (Type II supernovae) is introduced. Time evolution of spectral energy distributions from 1 A to 50,000 A (100 frequency bins in a log scale) is computed at each time step up to 500 days after explosion in each model. We provide light curves for the filters of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), Zwicky Transient Facility (ZTF), Sloan Digital Sky Servey (SDSS), and the Neil Gehrels Swift Observatory, but light curves for any photometric filters can be constructed by convolving any filter response functions to the synthetic spectral energy distributions. We also provide bolometric light curves and photosphere information such as photospheric velocity evolution. The parameter space covered by the model grid is five progenitor masses (10, 12, 14, 16, and 18 Msun at the zero-age main sequence, solar metallicity), ten explosion energies (0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0 x 10^51 erg), nine 56Ni masses (0.001, 0.01, 0.02, 0.04, 0.06, 0.08, 0.1, 0.2, and 0.3 Msun), nine mass-loss rates (1e-5.0, 1e-4.5, 1e-4.0, 1e-3.5, 1e-3.0, 1e-2.5, 1e-2.0, 1e-1.5, and 1e-1.0 Msun/yr with a wind velocity of 10 km/s), six circumstellar matter radii (1, 2, 4, 6, 8, and 10 x 10^14 cm), and ten circumstellar structures (beta = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0). 56Ni is assumed to be uniformly mixed up to the half mass of a hydrogen-rich envelope. This model grid can be a base for rapid characterizations of Type II supernovae with sparse photometric sampling expected in LSST through a Bayesian approach, for example. The model grid is available at https://doi.org/10.5061/dryad.pnvx0k6sj., Comment: 15 pages, 10 figures, 2 tables, data available at https://doi.org/10.5061/dryad.pnvx0k6sj

Early-time light curves/spectra of some hydrogen-rich supernovae (SNe) give firm evidence on the existence of confined, dense circumstellar matter (CSM) surrounding dying massive stars. We numerically and analytically study radiative acceleration of CSM in such systems, where the radiation is mainly powered by the interaction between the SN ejecta and the CSM. We find that the acceleration of the unshocked dense CSM ahead of the shock is larger for massive and compact CSM, with velocities reaching up to $\sim 10^3\ {\rm km\ s^{-1}}$ for a CSM of order $0.1\ M_\odot$ confined within $\sim 10^{15}$ cm. We show that the dependence of the acceleration on the CSM density helps us explain the diversity of the CSM velocity inferred from the early spectra of some Type II SNe. For explosions in even denser CSM, radiative acceleration can affect the dissipation of strong collisionless shocks formed after the shock breakout, which would affect early non-thermal emission expected from particle acceleration., Comment: 14 pages, 8 figures. Accepted to ApJ after minor revision

We analyze a sample of 45 Type II supernovae from the Zwicky Transient Facility (ZTF) public survey using a grid of hydrodynamical models in order to assess whether theoretically-driven forecasts can intelligently guide follow up observations supporting all-sky survey alert streams. We estimate several progenitor properties and explosion physics parameters including zero-age-main-sequence (ZAMS) mass, mass-loss rate, kinetic energy, 56Ni mass synthesized, host extinction, and the time of explosion. Using complete light curves we obtain confident characterizations for 34 events in our sample, with the inferences of the remaining 11 events limited either by poorly constraining data or the boundaries of our model grid. We also simulate real-time characterization of alert stream data by comparing our model grid to various stages of incomplete light curves (t less than 25 days, t less than 50 days, all data), and find that some parameters are more reliable indicators of true values at early epochs than others. Specifically, ZAMS mass, time of explosion, steepness parameter beta, and host extinction are reasonably constrained with incomplete light curve data, whereas mass-loss rate, kinetic energy and 56Ni mass estimates generally require complete light curves spanning greater than 100 days. We conclude that real-time modeling of transients, supported by multi-band synthetic light curves tailored to survey passbands, can be used as a powerful tool to identify critical epochs of follow up observations. Our findings are relevant to identify, prioritize, and coordinate efficient follow up of transients discovered by Vera C. Rubin Observatory., Comment: 27 pages, 14 figures, Accepted to The Astrophysical Journal

SN 2018ivc is an unusual type II supernova (SN II). It is a variant of SNe IIL, which might represent a transitional case between SNe IIP with a massive H-rich envelope, and IIb with only a small amount of the H-rich envelope. However, SN 2018ivc shows an optical light curve evolution more complicated than canonical SNe IIL. In this paper, we present the results of prompt follow-up observations of SN 2018ivc with the Atacama Large Millimeter/submillimeter Array (ALMA). Its synchrotron emission is similar to that of SN IIb 1993J, suggesting that it is intrinsically an SN IIb-like explosion of a He star with a modest (~0.5 - 1 Msun) extended H-rich envelope. Its radio, optical, and X-ray light curves are explained primarily by the interaction between the SN ejecta and the circumstellar material (CSM); we thus suggest that it is a rare example (and the first involving the `canonical' SN IIb ejecta) for which the multi-wavelength emission is powered mainly by the SN-CSM interaction. The inner CSM density, reflecting the progenitor activity in the final decade, is comparable to that of SN IIb 2013cu that showed a flash spectral feature. The outer CSM density, and therefore the mass-loss rate in the final ~200 years, is larger than that of SN 1993J by a factor of ~5. We suggest that SN 2018ivc represents a missing link between SNe IIP and IIb/Ib/Ic in the binary evolution scenario., Comment: 31 pages, 14 figures, 3 tables. Accepted for publication in ApJ

Type Ia supernovae (SNe Ia) are thermonuclear explosions of degenerate white dwarf (WD) stars destabilized by mass accretion from a companion star, but the nature of their progenitors remains poorly understood. A way to discriminate between progenitor systems is through radio observations; a non-degenerate companion star is expected to lose material through winds or binary interaction prior to explosion, and the SN ejecta crashing into this nearby circumstellar material (CSM) should result in radio synchrotron emission. However, despite extensive efforts, no SN Ia has ever been detected at radio wavelengths, which suggests a clean environment and a companion star that is itself a degenerate WD star. Here we report on the study of SN 2020eyj, a SN Ia showing helium-rich CSM, as revealed by its spectral features, infrared emission and, for the first time in a SN Ia, a radio counterpart. Based on our modeling, we conclude the CSM likely originates from a single-degenerate (SD) binary system where a WD accretes material from a helium donor star, an often hypothesized formation channel for SNe Ia. We describe how comprehensive radio follow-up of SN 2020eyj-like SNe Ia can improve the constraints on their progenitor systems., Comment: 62 pages, 10 figures, 4 tables

It was recently proposed that exotic particles can trigger a new stellar instability which is analogous to the e-e+ pair instability if they are produced and reach equilibrium in the stellar plasma. In this study, we construct axion instability supernova (AISN) models caused by the new instability to predict their observational signatures. We focus on heavy axion-like particles (ALPs) with masses of ~400 keV--2 MeV and coupling with photons of g_{ag}~10^{-5} GeV^{-1}. It is found that the 56Ni mass and the explosion energy are significantly increased by ALPs for a fixed stellar mass. As a result, the peak times of the light curves of AISNe occur earlier than those of standard pair-instability supernovae by 10--20 days when the ALP mass is equal to the electron mass. Also, the event rate of AISNe is 1.7--2.6 times higher than that of pair-instability supernovae, depending on the high mass cutoff of the initial mass function., Comment: 10 pages, 6 figures, 1 table, accepted for publication in ApJ

We present kinematics of 6 local extremely metal-poor galaxies (EMPGs) with low metallicities ($0.016-0.098\ Z_{\odot}$) and low stellar masses ($10^{4.7}-10^{7.6} M_{\odot}$). Taking deep medium-high resolution ($R\sim7500$) integral-field spectra with 8.2-m Subaru, we resolve the small inner velocity gradients and dispersions of the EMPGs with H$\alpha$ emission. Carefully masking out sub-structures originated by inflow and/or outflow, we fit 3-dimensional disk models to the observed H$\alpha$ flux, velocity, and velocity-dispersion maps. All the EMPGs show rotational velocities ($v_{\rm rot}$) of 5--23 km s$^{-1}$ smaller than the velocity dispersions ($\sigma_{0}$) of 17--31 km s$^{-1}$, indicating dispersion-dominated ($v_{\rm rot}/\sigma_{0}=0.29-0.80<1$) systems affected by inflow and/or outflow. Except for two EMPGs with large uncertainties, we find that the EMPGs have very large gas-mass fractions of $f_{\rm gas}\simeq 0.9-1.0$. Comparing our results with other H$\alpha$ kinematics studies, we find that $v_{\rm rot}/\sigma_{0}$ decreases and $f_{\rm gas}$ increases with decreasing metallicity, decreasing stellar mass, and increasing specific star-formation rate. We also find that simulated high-$z$ ($z\sim 7$) forming galaxies have gas fractions and dynamics similar to the observed EMPGs. Our EMPG observations and the simulations suggest that primordial galaxies are gas-rich dispersion-dominated systems, which would be identified by the forthcoming James Webb Space Telescope (JWST) observations at $z\sim 7$., Comment: ApJ in Press

In this Letter, we report the discovery of an ultraluminous fast-evolving transient in rest-frame UV wavelengths, MUSSES2020J, soon after its occurrence by using the Hyper Suprime-Cam (HSC) mounted on the 8.2 m Subaru telescope. The rise time of about 5 days with an extremely high UV peak luminosity shares similarities to a handful of fast blue optical transients whose peak luminosities are comparable with the most luminous supernovae while their timescales are significantly shorter (hereafter "fast blue ultraluminous transient," FBUT). In addition, MUSSES2020J is located near the center of a normal low-mass galaxy at a redshift of 1.063, suggesting a possible connection between the energy source of MUSSES2020J and the central part of the host galaxy. Possible physical mechanisms powering this extreme transient such as a wind-driven tidal disruption event and an interaction between supernova and circumstellar material are qualitatively discussed based on the first multiband early-phase light curve of FBUTs, although whether the scenarios can quantitatively explain the early photometric behavior of MUSSES2020J requires systematical theoretical investigations. Thanks to the ultrahigh luminosity in UV and blue optical wavelengths of these extreme transients, a promising number of FBUTs from the local to the high-z universe can be discovered through deep wide-field optical surveys in the near future., Comment: 14 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Mass loss of massive helium stars is not well understood even though it plays an essential role in determining their remnant neutron-star or black-hole masses as well as ejecta mass of Type Ibc supernovae. Radio emission from Type Ibc supernovae is strongly affected by circumstellar matter properties formed by mass loss of their massive helium star progenitors. In this study, we estimate the rise time and peak luminosity distributions of Type Ibc supernovae in radio based on a few massive helium star mass-loss prescriptions and compare them with the observed distribution to constrain the uncertain massive helium star mass-loss rates. We find that massive helium stars in the luminosity range expected for ordinary Type Ibc supernova progenitors (4.6 ~< log L/Lsun ~< 5.2) should generally have large mass-loss rates (> ~ 1e-6 Msun/yr) in order to account for the observed rise time and peak luminosity distribution. Therefore, mass-loss prescriptions that predict significantly low mass-loss rates for helium stars in this luminosity range is inconsistent with the supernova radio observations. It is also possible that massive helium stars shortly before their explosion generally undergo mass-loss enhancement in a different way from the standard radiation-driven wind mechanism., Comment: 5 pages, 3 figures, accepted by Monthly Notices of the Royal Astronomical Society Letters

Explosion mechanism of core-collapse supernovae is not fully understood yet. In this work, we give constraints on the explosion timescale based on $^{56}$Ni synthesized by supernova explosions. First, we systematically analyze multi-band light curves of 82 stripped-envelope supernovae (SESNe) to obtain bolometric light curves, which is among the largest samples of the bolometric light curves of SESNe derived from the multi-band spectral energy distribution. We measure the decline timescale and the peak luminosity of the light curves and estimate the ejecta mass ($M_{\rm ej}$) and $^{56}$Ni mass ($M_{\rm Ni}$) to connect the observed properties with the explosion physics. We then carry out one-dimensional hydrodynamics and nucleosynthesis calculations, varying the progenitor mass and the explosion timescale. From the calculations, we show that the maximum $^{56}$Ni mass that $^{56}$Ni-powered SNe can reach is expressed as $M_{\rm Ni} \lesssim 0.2 \ M_{\rm ej}$. Comparing the results from the observations and the calculations, we show that the explosion timescale shorter than 0.3 sec explains the synthesized $^{56}$Ni mass of the majority of the SESNe., Comment: 21 pages, 15 figures, accepted for publication in ApJ

We infer the expected detection number of pair instability supernovae (PISNe) during the operation of the Euclid space telescope based on binary population models. Our models reproduce the global maximum of the rate at the primary BH mass of $\sim 9-10$ $M_\odot$, and the overall gradient of the primary BH mass distribution in the binary BH merger rate consistent with recent observations. We consider different PISN conditions depending on the $^{12}$C$(\alpha, \gamma)^{16}$O reaction rate. The fiducial and $3\sigma$ models adopt the standard and $3\sigma$-smaller reaction rate, respectively. Our fiducial model predicts that Euclid detects several hydrogen-poor PISNe. For the $3\sigma$ model, detection of $\sim 1$ hydrogen-poor PISN by Euclid is expected if the stellar mass distribution extends to $M_{\max} = 600 M_\odot$, but the expected number becomes significantly smaller if $M_{\max} = 300 M_\odot$. We may be able to distinguish the fiducial and $3\sigma$ models by the observed PISN rate. This will help us to constrain the origin of binary BHs and the reaction rate, although there remains degeneracy between $M_{\max}$ and the reaction rate. PISN ejecta mass estimates from light curves and spectra obtained by follow-up observations would be important to disentangle the degeneracy., Comment: 7 pages, 4 figures, 1 table, accepted for MNRAS letters

An ultra-stripped supernova (USSN) is a type of core-collapse SN explosion proposed to be a candidate formation site of a double neutron star (DNS) binary. We investigate the dynamical evolution of an ultra-stripped supernova remnant (USSNR), which should host a DNS at its center. By accounting for the mass-loss history of the progenitor binary using a model developed by a previous study, we construct the large-scale structure of the {circumstellar medium (CSM)} up to a radius $\sim 100\,{\rm pc}$, and simulate the explosion and subsequent evolution of a USSN surrounded by such a CSM environment. We find that the CSM encompasses an extended region characterized by a hot plasma with a temperature $\sim 10^8\,$K located around the termination shock of the wind from the progenitor binary ($\sim 10\,$pc), and the USSNR blastwave is drastically weakened while penetrating through this hot plasma. Radio continuum emission from a young USSNR is sufficiently bright to be detectable if it inhabits our Galaxy but faint compared to the observed Galactic SNRs, and thereafter declines in luminosity through adiabatic cooling. Within our parameter space, USSNRs typically exhibit a low radio luminosity and surface brightness compared to the known Galactic SNRs. Due to the small event rate of USSNe and their relatively short observable lifespan, we calculate that USSNRs account for only $\sim 0.1$-$1$ % of the total SNR population. This is consistent with the fact that no SNR hosting a DNS binary has been discovered in the Milky Way so far., Comment: 23 pages, 10 figures, 3 tables, accepted for publication in The Astrophysical Journal

The primordial He abundance $Y_\mathrm{P}$ is a powerful probe of cosmology. Currently, $Y_\mathrm{P}$ is best determined by observations of metal-poor galaxies, while there are only a few known local extremely metal-poor ($<0.1 Z_\odot$) galaxies (EMPGs) having reliable He/H measurements with HeI$\lambda$10830 near-infrared (NIR) emission. Here we present deep Subaru NIR spectroscopy for 10 EMPGs. Combining the existing optical data, He/H values of 5 out of the 10 EMPGs are reliably derived by the Markov chain Monte Carlo algorithm. Adding the existing 3 EMPGs and 51 moderately metal-poor ($0.1-0.4 Z_\odot$) galaxies with reliable He/H estimates, we obtain $Y_\mathrm{P}=0.2370^{+0.0034}_{-0.0033}$ by linear regression in the $\mathrm{(He/H)}-\mathrm{(O/H)}$ plane, where we increase the number of EMPGs from 3 to 8 anchoring He/H of the most metal-poor gas in galaxies. Although our $Y_\mathrm{P}$ measurement and previous measurements are consistent, our result is slightly ($\sim 1\sigma$) smaller due to our EMPGs. With our $Y_\mathrm{P}$ and the existing primordial deuterium $D_\mathrm{P}$ measurement, we constrain the effective number of neutrino species $N_\mathrm{eff}$ and the baryon-to-photon ratio $\eta$ showing $\gtrsim 1-2\sigma$ tensions with the Standard Model and Planck Collaboration et al. (2020). Motivated by the tensions, we allow the degeneracy parameter of electron-neutrino $\xi_e$ to vary as well as $N_\mathrm{eff}$ and $\eta$. We obtain $\xi_e = 0.05^{+0.03}_{-0.02}$, $N_\mathrm{eff}=3.11^{+0.34}_{-0.31}$, and $\eta\times10^{10}=6.08^{+0.06}_{-0.06}$ from the $Y_\mathrm{P}$ and $D_\mathrm{P}$ measurements with a prior of $\eta$ taken from Planck Collaboration et al. (2020). Our constraints suggest a lepton asymmetry and allow for a high value of $N_\mathrm{eff}$ within the $1\sigma$ level, which could mitigate the Hubble tension., Comment: 22 pages, 9 figures, accepted for publication in ApJ

Luminosity evolution of some stripped-envelope supernovae such as Type I superluminous supernovae is difficult to be explained by the canonical 56Ni nuclear decay heating. A popular alternative heating source is rapid spin down of strongly-magnetized rapidly-rotating neutron stars (magnetars). Recent observations have indicated that Type I superluminous supernovae often have bumpy light curves with multiple luminosity peaks. The cause of bumpy light curves is unknown. In this study, we investigate the possibility that the light-curve bumps are caused by variations of the thermal energy injection from magnetar spin down. We find that a temporal increase in the thermal energy injection can lead to multiple luminosity peaks. The multiple luminosity peaks caused by the variable thermal energy injection is found to be accompanied by significant increase in photospheric temperature, and photospheric radii are not significantly changed. We show that the bumpy light curves of SN 2015bn and SN 2019stc can be reproduced by temporarily increasing magnetar spin-down energy input by a factor of 2-3 for 5-20 days. However, not all the light-curve bumps are accompanied by the clear photospheric temperature increase as predicted by our synthetic models. In particular, the secondary light-curve bump of SN 2019stc is accompanied by a temporal increase in photospheric radii rather than temperature, which is not seen in our synthetic models. We, therefore, conclude that not all the light-curve bumps observed in luminous supernovae are caused by the variable thermal energy injection from magnetar spin down and some bumps are likely caused by a different mechanism., Comment: 9 pages, 7 figures, accepted by Monthly Notices of the Royal Astronomical Society

We present an analysis of the nebular spectra of 103 stripped envelope (SE) supernovae (SNe) collected from the literature and observed with the Subaru Telescope from 2002 to 2012, focusing on [O I] 6300, 6363. The line profile and width of [O I] are employed to infer the ejecta geometry and the expansion velocity of the inner core. These two measurements are then compared with the SN sub types, and further with the [O I]/[Ca II] ratio, which is used as an indicator of the progenitor CO core mass. Based on the best fit results of the [O I] profile, the objects are classified into different morphological groups, and we conclude that the deviation from spherical symmetry is a common feature for all types of SESNe. There is a hint (at about 1 sigma level) that the distributions of the line profile fractions are different between canonical SESNe and broad-line SNe Ic. A correlation between [O I] width and [O I]/[Ca II] ratio is discerned, indicating that the oxygen-rich material tends to expand faster for objects with a more massive CO core. Such a correlation can be utilized to constrain the relation between the progenitor mass and the kinetic energy of the explosion. Further, when [O I]/[Ca II] ratio increases, the fraction of objects with Gaussian [O I] profile increases, while those with double-peaked profile decreases. This phenomenon connects ejecta geometry and the progenitor CO core mass., Comment: 40 pages, 23 figures, 10 tables. Accepted for publication in ApJ

Type Ibn Supernovae (SNe Ibn) show signatures of strong interaction between the SN ejecta and hydrogen-poor circumstellar matter (CSM). Deriving the ejecta and CSM properties of SNe Ibn provides a great opportunity to study the final evolution of massive stars. In the present work, we present a light curve (LC) model for the ejecta-CSM interaction, taking into account the processes in which the high-energy photons originally created at the forward and reverse shocks are converted to the observed emission in the optical. The model is applied to a sample of SNe Ibn and `SN Ibn' rapidly evolving transients. We show that the characteristic post-peak behavior commonly seen in the SN Ibn LCs, where a slow decay is followed by a rapid decay, is naturally explained by the transition of the forward-shock property from cooling to adiabatic regime without introducing a change in the CSM density distribution. The (commonly-found) slope in the rapid decay phase indicates a steep CSM density gradient (rho_CSM ~ r^{-3}), inferring a rapid increase in the mass-loss rate toward the SN as a generic properties of the SN Ibn progenitors. From the derived ejecta and CSM properties, we argue that massive Wolf-Rayet stars with the initial mass of >~ 18 Msun can be a potential class of the progenitors. The present work also indicates existence of currently missing population of UV-bright rapid transients for which the final mass-loss rate is lower than the optical SNe Ibn, which can be efficiently probed by future UV missions., Comment: 30 pages, 12 figures. Accepted for publication in ApJ. A few typos and references corrected/updated

Thorne-Zytkow objects are stars that have a neutron star core with an extended hydrogen-rich envelope. Massive Thorne-Zytkow objects are proposed to explode when the nuclear reactions sustaining their structure are terminated by the exhaustion of the seed elements. In this paper, we investigate the observational properties of the possible Thorne-Zytkow object explosions. We find that Thorne-Zytkow object explosions are observed as long-duration transients lasting for several years. If the accretion disk triggering the explosions does not last for a long time, Thorne-Zytkow object explosions have a luminosity plateau with about 1e39 erg/s lasting for a few years, and then they suddenly become faint. They would be observed as vanished stars after a bright phase lasting for a few years. If the accretion disk is sustained for long time, the Thorne-Zytkow object explosions become as bright as supernovae. They would be observed as supernovae with rise times of several hundred days. We found that their photospheric velocities are 2000 km/s at most, much smaller than those found in supernovae. Supernovae with extremely long rise times such as HSC16aayt and SN 2008iy may be related to the explosions of Thorne-Zytkow objects., Comment: 5 pages, 6 figures, accepted by Monthly Notices of the Royal Astronomical Society

Massive stars play critical roles for the reionization of the Universe. Individual massive stars at the reionization epoch (z > 6) are too faint to observe and quantify their contributions to reionization. Some massive stars, however, explode as superluminous supernovae (SLSNe) or pair-instability supernovae (PISNe) that are luminous enough to observe even at z > 6 and allow for the direct characterization of massive star properties at the reionization epoch. In addition, hypothetical long-sought-after PISNe are expected to be present preferentially at high redshifts, and their discovery will have a tremendous impact on our understanding of massive star evolution and the formation of stellar mass black holes. The near-infrared Wide Field Instrument on Nancy Grace Roman Space Telescope will excel at discovering such rare high-redshift supernovae. In this work, we investigate the best survey strategy to discover and identify SLSNe and PISNe at z > 6 with Roman. We show that the combination of the F158 and F213 filters can clearly separate both SLSNe and PISNe at z > 6 from nearby supernovae through their colors and magnitudes. The limiting magnitudes are required to be 27.0 mag and 26.5 mag in the F158 and F213 filters, respectively, to identify supernovae at z > 6. If we conduct a 10 deg2 transient survey with these limiting magnitudes for 5 years with a cadence of one year, we expect to discover 22.5 +- 2.8 PISNe and 3.1 +- 0.3 SLSNe at z > 6, depending on the cosmic star-formation history. The same survey is estimated to discover 76.1 +- 8.2 PISNe and 9.1 +- 0.9 SLSNe at 5 < z < 6. Such a supernova survey requires the total observational time of approximately 525 hours in 5 years. The legacy data acquired with the survey will also be beneficial for many different science cases including the study of high-redshift galaxies., Comment: 15 pages, 9 figures, 1 table, accepted by The Astrophysical Journal

We present the results of ALMA band 3 observations of a nearby type Ic supernova (SN) 2020oi. Under the standard assumptions on the SN-circumstellar medium (CSM) interaction and the synchrotron emission, the data indicate that the CSM structure deviates from a smooth distribution expected from the steady-state mass loss in the very vicinity of the SN (~10^{15} cm), which is then connected to the outer smooth distribution (~10^{16} cm). This structure is further confirmed through the light curve modeling of the whole radio data set as combined with data at lower frequency previously reported. Being an explosion of a bare carbon-oxygen (C+O) star having a fast wind, we can trace the mass-loss history of the progenitor of SN 2020oi in the final year. The inferred non-smooth CSM distribution corresponds to fluctuations on the sub-year time scale in the mass-loss history toward the SN explosion. Our finding suggests that the pre-SN activity is likely driven by the accelerated change in the nuclear burning stage in the last moments just before the massive star's demise. The structure of the CSM derived in this study is beyond the applicability of the other methods at optical wavelengths, highlighting an importance and uniqueness of quick follow-up observations of SNe by ALMA and other radio facilities., Comment: 20 pages, 8 figures, 1 table. Accepted for publication in ApJ