Search

The cores of stars are the cosmic furnaces where light elements are fused into heavier nuclei. The fusion of hydrogen to helium initially powers all stars. The ashes of the fusion reactions are then predicted to serve as fuel in a series of stages, eventually transforming massive stars into a structure of concentric shells. These are composed of natal hydrogen on the outside, and consecutively heavier compositions inside, predicted to be dominated by helium, carbon/oxygen, oxygen/neon/magnesium, and oxygen/silicon/sulphur. Silicon and sulphur are fused into inert iron, leading to the collapse of the core and either a supernova explosion or the direct formation of a black hole. Stripped stars, where the outer hydrogen layer has been removed and the internal He-rich layer (in Wolf-Rayet WN stars) or even the C/O layer below it (in Wolf-Rayet WC/WO stars) are exposed, provide evidence for this shell structure, and the cosmic element production mechanism it reflects. The types of supernova explosions that arise from stripped stars embedded in shells of circumstellar material (most notably Type Ibn supernovae from stars with outer He layers, and Type Icn supernovae from stars with outer C/O layers) confirm this scenario. However, direct evidence for the most interior shells, which are responsible for the production of elements heavier than oxygen, is lacking. Here, we report the discovery of the first-of-its-kind supernova arising from a star peculiarly stripped all the way to the silicon and sulphur-rich internal layer. Whereas the concentric shell structure of massive stars is not under debate, it is the first time that such a thick, massive silicon and sulphur-rich shell, expelled by the progenitor shortly before the SN explosion, has been directly revealed., Comment: 48 pages, 12 figures and 10 tables. Submitted to a high-impact journal. The reduced spectra and photometry will be made available via the journal webpage and the WISeREP archive after the acceptance of the paper

The progenitor system(s) as well as the explosion mechanism(s) of thermonuclear (Type Ia) supernovae are long-standing issues in astrophysics. Here we present ejecta masses and other physical parameters for 28 recent Type Ia supernovae inferred from multiband photometric and optical spectroscopic data. Our results confirm that the majority of SNe Ia show {\it observable} ejecta masses below the Chandrasekhar-limit (having a mean $M_{\rm ej} \approx 1.1 \pm 0.3$ M$_\odot$), consistent with the predictions of recent sub-M$_{\rm Ch}$ explosion models. They are compatible with models assuming either single- or double-degenerate progenitor configurations. We also recover a sub-sample of supernovae within $1.2 $ M$_\odot$ $< M_{\rm {ej}} < 1.5$ M$_\odot$ that are consistent with near-Chandrasekhar explosions. Taking into account the uncertainties of the inferred ejecta masses, about half of our SNe are compatible with both explosion models. We compare our results with those in previous studies, and discuss the caveats and concerns regarding the applied methodology.

SN 2023ixf was discovered in M101 within a day of explosion and rapidly classified as a Type II supernova with flash features. Here we present ultraviolet (UV) spectra obtained with the Hubble Space Telescope 14, 19, 24, and 66 days after explosion. Interaction between the supernova ejecta and circumstellar material (CSM) is seen in the UV throughout our observations in the flux of the first three epochs and asymmetric Mg II emission on day 66. We compare our observations to CMFGEN supernova models that include CSM interaction ($\dot{M}<10^{-3}$ M$_{\odot}$ yr$^{-1}$) and find that the power from CSM interaction is decreasing with time, from $L_{\rm sh}\approx5\times10^{42}$ erg s$^{-1}$ to $L_{\rm sh}\approx1\times10^{40}$ erg s$^{-1}$ between days 14 and 66. We examine the contribution of individual atomic species to the spectra on days 14 and 19, showing that the majority of the features are dominated by iron, nickel, magnesium, and chromium absorption in the ejecta. The UV spectral energy distribution of SN 2023ixf sits between that of supernovae which show no definitive signs of CSM interaction and those with persistent signatures assuming the same progenitor radius and metallicity. Finally, we show that the evolution and asymmetric shape of the Mg II $\lambda\lambda$ 2796, 2802 emission are not unique to SN 2023ixf. These observations add to the early measurements of dense, confined CSM interaction, tracing the mass-loss history of SN 2023ixf to $\sim33$ yr prior to the explosion and the density profile to a radius of $\sim5.7\times10^{15}$ cm. They show the relatively short evolution from a quiescent red supergiant wind to high mass loss., Comment: Accepted ApJL

We present extensive photometric and spectroscopic observations of the nearby Type Ia supernova (SN) 2023wrk at a distance of about 40 Mpc. The earliest detection of this SN can be traced back to a few hours after the explosion. Within the first few days the light curve shows a bump feature, while the B - V color is blue and remains nearly constant. The overall spectral evolution is similar to that of an SN 1991T/SN 1999aa-like SN Ia, while the C II $\lambda6580$ absorption line appears to be unusually strong in the first spectrum taken at $t \approx -$15.4 days after the maximum light. This carbon feature disappears quickly in subsequent evolution but it reappears at around the time of peak brightness. The complex evolution of the carbon line and the possible detection of Ni III absorption around 4700 {\AA} and 5300 {\AA} in the earliest spectra indicate macroscopic mixing of fuel and ash. The strong carbon lines is likely related to collision of SN ejecta with unbound carbon, consistent with the predictions of pulsational delayed-detonation or carbon-rich circumstellar-matter interaction models. Among those carbon-rich SNe Ia with strong C II $\lambda6580$ absorption at very early times, the line-strength ratio of C II to Si II and the B-V color evolution are found to exhibit large diversity, which may be attributed to different properties of unbound carbon and outward-mixing $^{56}$Ni., Comment: Accepted for publication in the Astrophysical Journal (27 pages, 14 figures, 7 tables)

We present multiwavelength observations of the Swift short $\gamma$-ray burst GRB 231117A, localized to an underlying galaxy at redshift $z = 0.257$ at a small projected offset ($\sim 2~$kpc). We uncover long-lived X-ray (Chandra) and radio/millimeter (VLA, MeerKAT, and ALMA) afterglow emission, detected to $\sim 37~$days and $\sim 20~$days (rest frame), respectively. We measure a wide jet ($\sim 10.4^\circ$) and relatively high circumburst density ($\sim 0.07~{\rm cm}^{-3}$) compared to the short GRB population. Our data cannot be easily fit with a standard forward shock model, but they are generally well fit with the incorporation of a refreshed forward shock and a reverse shock at $< 1~$day. We incorporate GRB 231117A into a larger sample of 132 X-ray detected events, 71 of which were radio-observed (17 cm-band detections), for a systematic study of the distributions of redshifts, jet and afterglow properties, galactocentric offsets, and local environments of events with and without detected radio afterglows. Compared to the entire short GRB population, the majority of radio-detected GRBs are at relatively low redshifts ($z < 0.6$) and have high circumburst densities ($> 10^{-2}~{\rm cm}^{-3}$), consistent with their smaller ($< 8~$kpc) projected galactocentric offsets. We additionally find that 70% of short GRBs with opening angle measurements were radio-detected, indicating the importance of radio afterglows in jet measurements, especially in the cases of wide ($> 10^\circ$) jets where observational evidence of collimation may only be detectable at radio wavelengths. Owing to improved observing strategies and the emergence of sensitive radio facilities, the number of radio-detected short GRBs has quadrupled in the past decade., Comment: 30 pages, 11 figures, submitted to ApJ

We present extensive observations and analysis of supernova (SN) 2021dbg, utilizing optical photometry and spectroscopy. For approximately 385 days following the explosion, SN 2021dbg exhibited remarkable luminosity, surpassing most SNe II. This initial high luminosity is potentially attributed to the interaction between the ejected material and the surrounding circumstellar material (CSM), as evidenced by the pronounced interaction signatures observed in its spectra. The subsequent high luminosity is primarily due to the significant $^{56}$Ni ($0.17 \pm 0.05$ M$_{\odot}$) produced in the explosion. Based on the flux of flash emission lines detected in the initial spectra, we estimate that the CSM mass near the progenitor amounted to $\sim$(1.0--2.0) $\times 10^{-3}$ M$_{\odot}$, likely resulting from intense stellar wind activity 2--3 yr preceding the explosion. Considering the bolometric light curve, nebular spectrum modeling, and mass-loss rate, we suggest that the progenitor of SN 2021dbg was a red supergiant (RSG) with a mass of $\sim 20$ M$_{\odot}$ and a radius of 1200 R$_{\odot}$. This RSG featured a thick hydrogen shell, which may have contained a region with a sharp decrease in material density, electron density, and temperature, contributing to its layered structure. This object demonstrates mixed features of SNe IIP and SNe IIL, making it as a transitional event linking the above two subclasses of SNe II.

We present optical photometry for the afterglow of GRB 221009A, in some respects the most extraordinary gamma-ray burst (GRB) ever observed. Good quality in the R-band light curve is obtained, covering 0.32-19.57 days since the Fermi-GBM trigger. We find that a weak bump emerges fromthe declining afterglow at $t \approx 11$ days; a supernova (SN) may be responsible. We use a smooth broken power-law and $^{56}\mathrm{Ni}$ model to fit the light curve. The best-fitting results reveal that the SN ejected a total mass of $M_\mathrm{ej} = 3.70 M_\odot$, a $^{56}\mathrm{Ni}$ mass of $M_\mathrm{Ni} = 0.23 M_\odot$, and a kinetic energy of $E_\mathrm{SN,K} = 2.35 \times 10^{52} \mathrm{erg}$. We also compare GRB 221009A with other GRB-SN events based on a GRB-associated SN sample, and find that only SN 2003lw and SN 2011kl can be obviously revealed in the afterglow of GRB 221009A by setting these objects at its distance. This suggests that a supernova (SN 2022xiw) is possibly obscured by the brighter afterglow emission from GRB 221009A.

The reactivated Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE-R) serendipitously caught the Type II supernova SN 2023ixf in Messier 101 on the rise, starting day 3.6 through day 10.9, and on the late-time decline from days 211 through 213 and days 370 through 372. We have considered these mid-infrared (mid-IR) data together with observations from the ultraviolet (UV) through the near-IR, when possible. At day 3.6 we approximated the optical emission with a hot, ~26,630 K blackbody, with a notable UV excess likely from strong SN shock interaction with circumstellar matter (CSM). In the IR, however, a clear excess is also obvious, and we fit it with a cooler, ~1,620 K blackbody with radius of ~2.6 x 10^{15} cm, consistent with dust in the progenitor's circumstellar shell likely heated by the UV emission from the CSM interaction. On day 10.8, the light detected was consistent with SN ejecta-dominated emission. At late times we also observed a clear NEOWISE-R excess, which could arise either from newly formed dust in the inner ejecta or in the contact discontinuity between the forward and reverse shocks, or from more distant pre-existing dust grains in the SN environment. Furthermore, the large 4.6 micron excess at late times can also be explained by the emergence of the carbon monoxide 1--0 vibrational band. SN 2023ixf is the best-observed SN IIP in the mid-IR during the first several days after explosion and one of the most luminous such SNe ever seen., Comment: 22 pages, 8 figures, submitted to AAS Journals

Mrk 501 is a prototypical high-synchrotron-peaked blazar (HBL) and serves as one of the primary targets for the {\it Imaging X-ray Polarimetry Explorer} ({\it IXPE}). In this study, we report X-ray polarization measurements of Mrk 501 based on six {\it IXPE} observations. The detection of X-ray polarization at a confidence level exceeding 99\% is achieved in four out of the six observations conducted across the entire energy range (2--8 keV) of {\it IXPE}. The maximum polarization degree ($\Pi_{\rm X}$) is measured to be $15.8\%\pm2.8\%$, accompanied by a polarization angle ($\psi_{\rm X}$) of $98.0\deg\pm5.1\deg$ at a confidence level of $5.6 \sigma$. During the remaining two observations, only an upper limit of $\Pi_{\rm X}<$12\% could be derived at the 99\% confidence level. No temporal variability in polarization is observed throughout all six {\it IXPE} observations for Mrk 501. A discernible trend of energy-dependent variation in the polarization degree is detected in optical spectropolarimetry; however, no analogous indication is observed in $\Pi_{\rm X}$. The chromatic behavior of $\Pi$ and the consistent values of $\psi$ across different frequencies from X-rays to radio waves, along with the agreement between $\psi$ and jet position angle, strongly support the interpretation of the energy-stratified model with shock-accelerated particles in the jet of Mrk 501. Additionally, the possibility of the presence of a global helical magnetic field in the jet of Mrk 501 is discussed., Comment: 20 pages, 8 figures, 4 tables, accepted for publication in ApJL

Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secure example of this phenomenon to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising (t ~ 2 hr), luminous (M_UV,peak ~ -24.4 mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing strong limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of Gamma_0 ~ 30-50, significantly lower than in any well-observed GRB and consistent with the theoretically-predicted "dirty fireball" scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-Gamma, off-axis from a classical high-Gamma jet core. This event represents a milestone in orphan afterglow searches, demonstrating that luminous afterglows with weak or no detectable gamma-ray radiation exist in nature and can be discovered by high-cadence optical surveys., Comment: Submitted to MNRAS

Over recent decades, robotic (or highly automated) searches for supernovae (SNe) have discovered several thousand events, many of them in quite nearby galaxies (distances < 30 Mpc). Most of these SNe, including some of the best-studied events to date, were found before maximum brightness and have associated with them extensive follow-up photometry and spectroscopy. Some of these discoveries are so-called SN impostors, thought to be superoutbursts of luminous blue variable stars, although possibly a new, weak class of massive-star explosions. We conducted a Snapshot program with the Hubble Space Telescope(HST) and obtained images of the sites of 31 SNe and four impostors, to acquire late-time photometry through two filters. The primary aim of this project was to reveal the origin of any lingering energy for each event, whether it is the result of radioactive decay or, in some cases, ongoing late-time interaction of the SN shock with pre-existing circumstellar matter, or the presence of a light echo. Alternatively, lingering faint light at the SN position may arise from an underlying stellar population (e.g., a host star cluster, companion star, or a chance alignment). The results from this study complement and extend those from Snapshot programs by various investigators in previous HST cycles., Comment: Added details on paper: 37 pages, 42 figures, accepted for publication in ApJ

Binary evolution theory predicts that the second common envelope (CE) ejection can produce low-mass (0.32-0.36 Msun) subdwarf B (sdB) stars inside ultrashort-orbital-period binary systems, as their helium cores are ignited under nondegenerate conditions. With the orbital decay driven by gravitational-wave (GW) radiation, the minimum orbital periods of detached sdB binaries could be as short as ~20 minutes. However, only four sdB binaries with orbital periods below an hour have been reported so far, while none of them has an orbital period approaching the above theoretical limit. Here we report the discovery of a 20.5-minute-orbital-period ellipsoidal binary, TMTS J052610.43+593445.1, in which the visible star is being tidally deformed by an invisible carbon-oxygen white dwarf (WD) companion. The visible component is inferred to be an sdB star with a mass of ~0.33 Msun, approaching that of helium-ignition limit, although a He-core WD cannot be completely ruled out. In particular, the radius of this low-mass sdB star is only 0.066 Rsun, about seven Earth radii, possibly representing the most compact nondegenerate star ever known. Such a system provides a key clue to map the binary evolution scheme from the second CE ejection to the formation of AM CVn stars having a helium-star donor, and it will also serve as a crucial verification binary of space-borne GW detectors in the future., Comment: 24 pages, 11 figures, 1 table, published on Nature Astronomy, URL: https://www.nature.com/articles/s41550-023-02188-2

AT 2019azh is a H+He tidal disruption event (TDE) with one of the most extensive ultraviolet and optical data sets available to date. We present our photometric and spectroscopic observations of this event starting several weeks before and out to approximately two years after the g-band peak brightness and combine them with public photometric data. This extensive data set robustly reveals a change in the light-curve slope and a possible bump in the rising light curve of a TDE for the first time, which may indicate more than one dominant emission mechanism contributing to the pre-peak light curve. Indeed, we find that the MOSFiT-derived parameters of AT 2019azh, which assume reprocessed accretion as the sole source of emission, are not entirely self-consistent. We further confirm the relation seen in previous TDEs whereby the redder emission peaks later than the bluer emission. The post-peak bolometric light curve of AT 2019azh is better described by an exponential decline than by the canonical t^{-5/3} (and in fact any) power-law decline. We find a possible mid-infrared excess around the peak optical luminosity, but cannot determine its origin. In addition, we provide the earliest measurements of the Halpha emission-line evolution and find no significant time delay between the peak of the V-band light curve and that of the Halpha luminosity. These results can be used to constrain future models of TDE line formation and emission mechanisms in general. More pre-peak 1-2 days cadence observations of TDEs are required to determine whether the characteristics observed here are common among TDEs. More importantly, detailed emission models are needed to fully exploit such observations for understanding the emission physics of TDEs., Comment: Submitted to ApJ

We present photometry and spectroscopy of the slowly evolving superluminous Type IIn SN2015da. SN2015da is extraordinary for its very high peak luminosity, and also for sustaining a high luminosity for several years. Even at 8\,yr after explosion, SN2015da remains as luminous as the peak of a normal SNII-P. The total radiated energy integrated over this time period (with no bolometric correction) is at least 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 Msun of H-rich CSM, which in turn implies a massive progenitor system above 30 Msun. Narrow P Cyg features show steady CSM expansion at 90 km/s, requiring a high average mass-loss rate of roughly 0.1 Msun/yr sustained for 2 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, disfavor a pulsational pair instability model for the pre-SN mass loss. Instead, violent pre-SN binary interaction is a likely cuprit. Finally, SN2015da exhibits the characteristic asymmetric blueshift in its emission lines from shortly after peak until the present epoch, adding another well-studied superluminous SNeIIn with unambiguous evidence of post-shock dust formation., Comment: 18 pages, 11 figs. submitted

In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole., Comment: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

The origins of quasi-periodic eruptions (QPEs) are poorly understood, although most theoretical explanations invoke an accretion disk around a supermassive black hole. The gas and stellar environments in the galactic nuclei of these sources are also poorly constrained. In this paper, we present an analysis of archival Hubble Space Telescope (HST) images to study the narrow-line [O III] emission in the QPE source GSN 069. We find strong evidence for a compact nuclear [O III] emission region of size $\lesssim 35$ pc, overlaid on top of extended [O III] emission up to 2 kpc away from the nucleus. The age of the accretion system is estimated to be between 10 and 100 yr. The [O III] luminosity of the compact region was measured to be $(2.1 \pm 0.3) \times 10^{40}\,\rm erg\,s^{-1}$. Based on CLOUDY simulations, we constrain that the [O III] emitting gas has a hydrogen number density in the range $5 \times 10^{3} < n_{\rm H} \lesssim 10^{8}\,\rm cm^{-3}$ and volume filling factor $f_{\rm V} < 2 \times 10^{-3}$. We suggest that the dense gas in the nuclear region of GSN 069 originates from molecular clouds (with total mass $\gtrsim 3 \times 10^{3}\,M_{\odot}$), which are freshly ionised by the soft X-ray photons from the accretion disk. We predict possible evolution of the compact narrow-line region on emission-line diagnostic diagrams, and hence future HST or integral-field unit observations can be used to further pin down the age of this puzzling system., Comment: 10 pages, 7 figures. Comments welcome

We present the discovery and studies of the helium-rich, fast-evolving supernova (SN) 2021agco at a distance of $\sim$ 40 Mpc. Its early-time flux is found to rise from half peak to the peak of $-16.06\pm0.42$ mag in the $r$ band within $2.4^{+1.5}_{-1.0}$ days, and the post-peak light curves also decline at a much faster pace relative to normal stripped-envelope SNe of Type Ib/Ic. The early-time spectrum of SN~2021agco ($t \approx 1.0$ days after the peak) is characterized by a featureless blue continuum superimposed with a weak emission line of ionized C III, and the subsequent spectra show prominent He I lines. Both the photometric and spectroscopic evolution shows close resemblances to SN 2019dge, which is believed to have an extremely stripped progenitor. We reproduce the multicolor light curves of SN 2021agco with a model combining shock-cooling emission with \Ni decay. The best-fit results give an ejecta mass of $\approx 0.3$~M$_\odot$ and a synthesized nickel mass of $\approx 2.2\times10^{-2}$~M$_\odot$. The progenitor is estimated to have an envelope radius $R_{\rm env} \approx 80$~R$_\odot$ and a mass $M_{\rm env} \approx 0.10$~M$_\odot$. All these suggest that SN~2021agco can be categorized as an ultrastripped SN~Ib, representing the closest object of this rare subtype. This SN is found to explode in the disk of an Sab-type galaxy with an age of $\sim 10.0$~Gyr and low star-forming activity. Compared to normal SNe Ib/c, the host galaxies of SN 2021agco and other ultrastripped SNe tend to have relatively lower metallicity, which complicates the properties of their progenitor populations., Comment: 20 figures, 10 figures

We present mid-infrared (mid-IR) imaging of the Type IIL supernova (SN) 1980K with the James Webb Space Telescope (JWST) more than 40 yr post-explosion. SN 1980K, located in the nearby ($D\approx7$ Mpc) "SN factory" galaxy NGC 6946, was serendipitously captured in JWST/MIRI images taken of the field of SN 2004et in the same galaxy. SN 1980K serves as a promising candidate for studying the transitional phase between young SNe and older SN remnants and also provides a great opportunity to investigate its the close environment. SN 1980K can be identified as a clear and bright point source in all eight MIRI filters from F560W up to F2550W. We fit analytical dust models to the mid-IR spectral energy distribution that reveal a large amount ($M_d \approx 0.002 {M}_{\odot}$) of Si-dominated dust at $T_{dust}\approx 150$ K (accompanied by a hotter dust/gas component), and also computed numerical SED dust models. Radiative transfer modeling of a late-time optical spectrum obtained recently with Keck discloses that an even larger ($\sim 0.24-0.58~{M}_{\odot}$) amount of dust is needed in order for selective extinction to explain the asymmetric line profile shapes observed in SN 1980K. As a conclusion, with JWST, we may see i) pre-existing circumstellar dust heated collisionally (or, partly radiatively), analogous to the equatorial ring of SN 1987A, or ii) the mid-IR component of the presumed newly-formed dust, accompanied by much more colder dust present in the ejecta (as suggested by the late-time the optical spectra)., Comment: 14 pages, 9 figures

Light echoes give us a unique perspective on the nature of supernovae and non-terminal stellar explosions. Spectroscopy of light echoes can reveal details on the kinematics of the ejecta, probe asymmetry, and reveal details on its interaction with circumstellar matter, thus expanding our understanding of these transient events. However, the spectral features arise from a complex interplay between the source photons, the reflecting dust geometry, and the instrumental setup and observing conditions. In this work we present an improved method for modeling these effects in light echo spectra, one that relaxes the simplifying assumption of a light curve weighted sum, and instead estimates the true relative contribution of each phase. We discuss our logic, the gains we obtain over light echo analysis method(s) used in the past, and prospects for further improvements. Lastly, we show how the new method improves our analysis of echoes from Tycho's supernova (SN 1572) as an example.

We present optical and near-infrared observations of SN~2022crv, a stripped envelope supernova in NGC~3054, discovered within 12 hrs of explosion by the Distance Less Than 40 Mpc Survey. We suggest SN~2022crv is a transitional object on the continuum between SNe Ib and SNe IIb. A high-velocity hydrogen feature ($\sim$$-$20,000 -- $-$16,000 $\rm km\,s^{-1}$) was conspicuous in SN~2022crv at early phases, and then quickly disappeared around maximum light. By comparing with hydrodynamic modeling, we find that a hydrogen envelope of $\sim 10^{-3}$ \msun{} can reproduce the behaviour of the hydrogen feature observed in SN~2022crv. The early light curve of SN~2022crv did not show envelope cooling emission, implying that SN~2022crv had a compact progenitor with extremely low amount of hydrogen. The analysis of the nebular spectra shows that SN~2022crv is consistent with the explosion of a He star with a final mass of $\sim$4.5 -- 5.6 \msun{} that has evolved from a $\sim$16 -- 22 \msun{} zero-age main sequence star in a binary system with about 1.0 -- 1.7 \msun{} of oxygen finally synthesized in the core. The high metallicity at the supernova site indicates that the progenitor experienced a strong stellar wind mass loss. In order to retain a small amount of residual hydrogen at such a high metallicity, the initial orbital separation of the binary system is likely larger than $\sim$1000~$\rm R_{\odot}$. The near-infrared spectra of SN~2022crv show a unique absorption feature on the blue side of He I line at $\sim$1.005~$\mu$m. This is the first time that such a feature has been observed in a Type Ib/IIb, and could be due to \ion{Sr}{2}. Further detailed modelling on SN~2022crv can shed light on the progenitor and the origin of the mysterious absorption feature in the near infrared., Comment: 33 pages, 23 figures, submitted to ApJ

We present extensive photometric and spectroscopic observations of a peculiar type Ia supernova (SN Ia) 2022vqz. It shares many similarities with the SN 2002es-like SNe Ia, such as low luminosity (i.e., $M_{B,\rm max}=-18.11\pm0.16$ mag) and moderate post-peak decline rate (i.e., $\Delta m_{15,B}=1.33\pm0.11$ mag). The nickel mass synthesized in the explosion is estimated as $0.20\pm0.04~{\rm M}_\odot$ from the bolometric light curve, which is obviously lower than normal SNe Ia. SN 2022vqz is also characterized by a slow expanding ejecta, with Si II velocities persisting around 7000 km s$^{-1}$ since 16 days before the peak, which is unique among all known SNe Ia. While all these properties imply a less energetic thermonuclear explosion that should leave considerable amount of unburnt materials, however, absent signature of unburnt carbon in the spectra of SN 2022vqz is puzzling. A prominent early peak is clearly detected in the $c$- and $o$-band light curves of ATLAS and in the $gr$-band data of ZTF within days after the explosion. Possible mechanisms for the early peak are discussed, including sub-Chandrasekhar mass double detonation model and interaction of SN ejecta with circumstellar material (CSM). We found both models face some difficulties in replicating all aspects of the observed data. As an alternative, we propose a hybrid CONe white dwarf as progenitor of SN 2022vqz which can simultaneously reconcile the tension between low ejecta velocity and absence of carbon. We further discuss the diversity of 02es-like objects and possible origins of different scenarios., Comment: 24 pages, 12 figures, submitted to MNRAS

Late-time spectra of Type Ia supernovae (SNe Ia) are important in clarifying the physics of their explosions, as they provide key clues to the inner structure of the exploding white dwarfs. We examined late-time optical spectra of 36 SNe Ia, including five from our own project (SNe 2019np, 2019ein, 2021hpr, 2021wuf, and 2022hrs), with phase coverage of $\sim 200$ to $\sim 400$ days after maximum light. At this late phase, the outer ejecta have become transparent and the features of inner iron-group elements emerge in the spectra. Based on multicomponent Gaussian fits and reasonable choices for the pseudocontinuum around Ni and Fe emission features, we get reliable estimates of the Ni to Fe ratio, which is sensitive to the explosion models of SNe Ia. Our results show that the majority (about 67%) of our SNe Ia are more consistent with the sub-Chandrasekhar-mass (i.e., double-detonation) model, although they could be affected by evolutionary or ionisation effects. Moreover, we find that the Si II $\lambda$6355 velocity measured around the time of maximum light tends to increase with the Ni to Fe ratio for the subsample with either redshifted or blueshifted nebular velocities, suggesting that progenitor metallicity might play an important role in accounting for the observed velocity diversity of SNe Ia., Comment: 20 pages, 10 figures, accepted for publication in MNRAS

We follow our first paper with an analysis of the ensemble of the extensive pre-explosion ground- and space-based infrared observations of the red supergiant (RSG) progenitor candidate for the nearby core-collapse supernova SN 2023ixf in Messier 101, together with optical data prior to explosion obtained with the Hubble Space Telescope (HST). We have confirmed the association of the progenitor candidate with the SN, as well as constrained the metallicity at the SN site, based on SN observations with instruments at Gemini-North. The internal host extinction to the SN has also been confirmed from a high-resolution Keck spectrum. We fit the observed spectral energy distribution (SED) for the star, accounting for its intrinsic variability, with dust radiative-transfer modeling, which assume a silicate-rich dust shell ahead of the underlying stellar photosphere. The star is heavily dust-obscured, likely the dustiest progenitor candidate yet encountered. We found median estimates of the star's effective temperature and luminosity of 2770 K and 9.0e4 L_Sun, with 68% credible intervals of 2340--3150 K and (7.5--10.9)e4 L_sun. The candidate may have a Galactic RSG analog, IRC -10414, with a strikingly similar SED and luminosity. Via comparison with single-star evolutionary models we have constrained the initial mass of the progenitor candidate from 12 M_sun to as high as 14 M_sun. We have had available to us an extraordinary view of the SN 2023ixf progenitor candidate, which should be further followed up in future years with HST and the James Webb Space Telescope., Comment: 40 pages, substantive modifications relative to the previous, although the overall conclusions remain the same; to appear in AAS Journals

Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon-oxygen rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity ($M_{B}=-18.9$ mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak $B$-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O I] $\lambda\lambda 6300,\ 6364$ (${\rm FWHM} \approx 2{,}000$ km s$^{-1}$), strong, broad emission from [Ca II] $\lambda\lambda 7291,\ 7323$ (${\rm FWHM} \approx 7{,}300$ km s$^{-1}$), and a rapid Fe III to Fe II ionization change. Finally, we present the first-ever optical-to-mid-infrared (MIR) nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with $T \approx 500$ K), combined with prominent [O I] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within carbon/oxygen-rich CSM., Comment: 23 pages, 11 figures, submitted to ApJ

We present an analysis of ground-based and JWST observations of SN~2022pul, a peculiar "03fg-like" (or "super-Chandrasekhar") Type Ia supernova (SN Ia), in the nebular phase at 338d post explosion. Our combined spectrum continuously covers 0.4--14 $\mu$m and includes the first mid-infrared spectrum of an 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. Strong, broad, centrally peaked [Ne II] line at 12.81 $\mu$m was previously predicted as a hallmark of "violent merger'' SN Ia models, where dynamical interaction between two sub-$M_{ch}$ white dwarfs (WDs) causes disruption of the lower mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to better reproduce the optical iron emission, and add mass in the innermost region ($< 2000$ km s$^{-1}$) to account for the observed narrow [O I]~$\lambda\lambda6300$, 6364 emission. A violent WD-WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SN Ia., Comment: 20 pages, 10 figures, published in ApJ

We present observations of SN 2022joj, a peculiar Type Ia supernova (SN Ia) discovered by the Zwicky Transient Facility (ZTF). SN 2022joj exhibits an unusually red $g_\mathrm{ZTF}-r_\mathrm{ZTF}$ color at early times and a rapid blueward evolution afterwards. Around maximum brightness, SN 2022joj shows a high luminosity ($M_{g_\mathrm{ZTF},\mathrm{max}}\simeq-19.7$ mag), a blue broadband color ($g_\mathrm{ZTF}-r_\mathrm{ZTF}\simeq-0.2$ mag), and shallow Si II absorption lines, consistent with those of overluminous, SN 1991T-like events. The maximum-light spectrum also shows prominent absorption around 4200 \r{A}, which resembles the Ti II features in subluminous, SN 1991bg-like events. Despite the blue optical-band colors, SN 2022joj exhibits extremely red ultraviolet $-$ optical colors at maximum luminosity ($u-v\simeq1.6$ mag and $uvw1 - v\simeq4.0$ mag), suggesting a suppression of flux between $\sim$2500--4000 \r{A}. Strong C II lines are also detected at peak. We show that these unusual spectroscopic properties are broadly consistent with the helium-shell double detonation of a sub-Chandrasekhar mass ($M\simeq1\mathrm{M_\odot}$) carbon/oxygen (C/O) white dwarf (WD) from a relatively massive helium shell ($M_s\simeq0.04$--$0.1\mathrm{M_\odot}$), if observed along a line of sight roughly opposite to where the shell initially detonates. None of the existing models could quantitatively explain all the peculiarities observed in SN 2022joj. The low flux ratio of [Ni II] $\lambda$7378 to [Fe II] $\lambda$7155 emission in the late-time nebular spectra indicates a low yield of stable Ni isotopes, favoring a sub-Chandrasekhar mass progenitor. The significant blueshift measured in the [Fe II] $\lambda$7155 line is also consistent with an asymmetric chemical distribution in the ejecta, as is predicted in double-detonation models., Comment: 24 pages, 11 figures, 6 tables. Submitted to ApJ

We present highly sampled photometry of the supernova (SN) 2023ixf, a Type II SN in M101, beginning 2 days before its first known detection. To gather these data, we enlisted the global Unistellar Network of citizen scientists. These 252 observations from 115 telescopes show the SN's rising brightness associated with shock emergence followed by gradual decay. We measure a peak $M_{V}$ = -18.18 $\pm$ 0.09 mag at 2023-05-25 21:37 UTC in agreement with previously published analyses., Comment: 4 pages, 1 figure

We present six epochs of optical spectropolarimetry of the Type II supernova (SN) 2023ixf ranging from $\sim$ 2 to 15 days after the explosion. Polarimetry was obtained with the Kast double spectrograph on the Shane 3 m telescope at Lick Observatory, representing the earliest such observations ever captured for an SN. We observe a high continuum polarization $p_{\text{cont}} \approx 1$ % on days +1.4 and +2.5 before dropping to 0.5 % on day +3.5, persisting at that level up to day +14.5. Remarkably, this change coincides temporally with the disappearance of highly ionized "flash" features. The decrease of the continuum polarization is accompanied by a $\sim 70^\circ$ rotation of the polarization position angle ($PA$) as seen across the continuum. The early evolution of the polarization may indicate different geometric configurations of the electron-scattering atmosphere as seen before and after the disappearance of the emission lines associated with highly-ionized species (e.g., He II, C IV, N III), which are likely produced by elevated mass loss shortly prior to the SN explosion. We interpret the rapid change of polarization and $PA$ from days +2.5 to +4.5 as the time when the SN ejecta emerge from the dense asymmetric circumstellar material (CSM). The temporal evolution of the continuum polarization and the $PA$ is consistent with an aspherical SN explosion that exhibits a distinct geometry compared to the CSM. The rapid follow-up spectropolarimetry of SN 2023ixf during the shock ionization phase reveals an exceptionally asymmetric mass-loss process leading up to the explosion., Comment: Accepted to ApJL

The AGN STORM 2 campaign is a large, multiwavelength reverberation mapping project designed to trace out the structure of Mrk 817 from the inner accretion disk to the broad emission line region and out to the dusty torus. As part of this campaign, Swift performed daily monitoring of Mrk 817 for approximately 15 months, obtaining observations in X-rays and six UV/optical filters. The X-ray monitoring shows that Mrk 817 was in a significantly fainter state than in previous observations, with only a brief flare where it reached prior flux levels. The X-ray spectrum is heavily obscured. The UV/optical light curves show significant variability throughout the campaign and are well correlated with one another, but uncorrelated with the X-rays. Combining the Swift UV/optical light curves with Hubble UV continuum light curves, we measure interband continuum lags, $\tau(\lambda)$, that increase with increasing wavelength roughly following $\tau(\lambda) \propto \lambda^{4/3}$, the dependence expected for a geometrically thin, optically thick, centrally illuminated disk. Modeling of the light curves reveals a period at the beginning of the campaign where the response of the continuum is suppressed compared to later in the light curve - the light curves are not simple shifted and scaled versions of each other. The interval of suppressed response corresponds to a period of high UV line and X-ray absorption, and reduced emission line variability amplitudes. We suggest that this indicates a significant contribution to the continuum from the broad line region gas that sees an absorbed ionizing continuum., Comment: 20 pages, 13 figures, 3 tables, accepted for publication in ApJ

We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities $\gtrsim 1300\,\rm km\,s^{-1}$. These are most likely the surviving donors from double-degenerate binaries in which the other WD exploded. The other two objects have lower minimum velocities, $\gtrsim 600\,\rm km\,s^{-1}$, and may have formed through a different mechanism, such as pure deflagration of a WD in a Type Iax supernova. The four fastest stars are hotter and smaller than the previously known "D$^6$ stars," with effective temperatures ranging from $\sim$20,000 to $\sim$130,000 K and radii of $\sim 0.02-0.10\,R_{\odot}$. Three of these have carbon-dominated atmospheres, and one has a helium-dominated atmosphere. Two stars have RVs of $-1694$ and $-2285\rm \,km\,s^{-1}$ -- the fastest systemic stellar RVs ever measured. Their inferred birth velocities, $\sim 2200-2500\,\rm km\,s^{-1}$, imply that both WDs in the progenitor binary had masses $>1.0\,M_{\odot}$. The high observed velocities suggest that a dominant fraction of the observed hypervelocity WD population comes from double-degenerate binaries whose total mass significantly exceeds the Chandrasekhar limit. However, the two nearest and faintest D$^6$ stars have the lowest velocities and masses, suggesting that observational selection effects favor rarer, higher-mass stars. A significant population of fainter low-mass runaways may still await discovery. We infer a birth rate of D$^6$ stars that is consistent with the SN Ia rate. The birth rate is poorly constrained, however, because the luminosities and lifetimes of $\rm D^6$ stars are uncertain., Comment: 26 pages, 17 figures. Accepted to OJA

We search the archival Zwicky Transient Facility public survey for rapidly evolving transient (RET) candidates based on well-defined criteria between 2018 May and 2021 December. The search yielded 19 bona-fide RET candidates, corresponding to a discovery rate of $\sim 5.2$ events per year. Even with a Galactic latitude cut of $20^\circ$, 8 of the 19 events ($\sim 42$%) are Galactic, including one with a light-curve shape closely resembling that of the GW170817 kilonova (KN). An additional event is a nova in M31. Four out of the 19 events ($\sim 21$%) are confirmed extragalactic RETs (one confirmed here for the first time) and the origin of 6 additional events cannot be determined. We did not find any extragalactic events resembling the GW170817 KN, from which we obtain an upper limit on the volumetric rate of GW170817-like KNe of $R \le$ 2400 Gpc$^{-3}$ yr$^{-1}$ (95% confidence). These results can be used for quantifying contaminants to RET searches in transient alert streams, specifically when searching for kilonovae independently of gravitational-wave and gamma-ray-burst triggers., Comment: 16 pages, 7 figures, submitted to AAS Journals

We have conducted photometric and spectroscopic observations of the peculiar Type Ia supernova (SN Ia) 2016ije that was discovered through the Tsinghua-NAOC Transient Survey. This peculiar object exploded in the outskirts of a metal-poor, low-surface brightness galaxy (i.e., $M_{g}$ = $-$14.5 mag). Our photometric analysis reveals that SN 2016ije is subluminous ($M_{B,\rm{max}}$ = $-$17.65$\pm$0.06 mag) but exhibits relatively broad light curves (${\Delta}m_{15}(B)$ = 1.35$\pm$0.14 mag), similar to the behavior of SN 2002es. Our analysis of the bolometric light curve indicates that only 0.14$\pm$0.04 $M_{\odot}$ of $^{56}$Ni was synthesized in the explosion of SN 2016ije, which suggests a less energetic thermonuclear explosion when compared to normal SNe Ia, and this left a considerable amount of unburned materials in the ejecta. Spectroscopically, SN 2016ije resembles other SN 2002es-like SNe Ia, except that the ejecta velocity inferred from its carbon absorption line ($\sim$ 4500 km s$^{-1}$) is much lower than that from silicon lines ($\sim$ 8300 km s$^{-1}$) at around the maximum light. Additionally, most of the absorption lines are broader than other 02es-like SNe Ia. These peculiarities suggest the presence of significant unburned carbon in the inner region and a wide line-forming region along the line of sight. These characteristics suggest that SN 2016ije might originate from the violent merger of a white dwarf binary system, when viewed near an orientation along the iron-group-element cavity caused by the companion star., Comment: 25 pages, 13 figures

Abridged - Stars with ZAMS masses between 140 and $260 M_\odot$ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN2018ibb is a H-poor SLSN at $z=0.166$ that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the NIR with 2-10m class telescopes. SN2018ibb radiated $>3\times10^{51} \rm erg$ during its evolution, and its bolometric light curve reached $>2\times10^{44} \rm erg\,s^{-1}$ at peak. The long-lasting rise of $>93$ rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source ($^{56}$Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions for their photometric and spectroscopic properties. SN2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25-44 $M_\odot$ of freshly nucleosynthesised $^{56}$Ni, pointing to the explosion of a metal-poor star with a He-core mass of 120-130 $M_\odot$ at the time of death. This interpretation is also supported by the tentative detection of [Co II]$\lambda$1.025$\mu$m, which has never been observed in any other PISN candidate or SLSN before. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN2018ibb by far the best candidate for being a PISN, to date., Comment: Accepted in A&A, the revised version includes a PISN rate estimate and an additional test with PISN models. 47 pages, main text 41 pages, 38 figures, 16 Tables

A growing number of supernovae (SNe) are now known to exhibit evidence for significant interaction with a dense, pre-existing, circumstellar medium (CSM). SNe Ibn comprise one such class that can be characterised by both rapidly evolving light curves and persistent narrow He I lines. The origin of such a dense CSM in these systems remains a pressing question, specifically concerning the progenitor system and mass-loss mechanism. In this paper, we present multi-wavelength data of the Type Ibn SN 2020nxt, including $HST$/STIS ultraviolet spectra. We fit the data with recently updated CMFGEN models designed to handle configurations for SNe Ibn. The UV coverage yields strong constraints on the energetics and, when combined with the CMFGEN models, offer new insight on potential progenitor systems. We find the most successful model is a $\lesssim4 {\rm M}_\odot$ helium star that lost its $\sim 1\,{\rm M}_\odot$ He-rich envelope in the years preceding core collapse. We also consider viable alternatives, such as a He white dwarf merger. Ultimately, we conclude at least some SNe Ibn do not arise from single, massive ($>30 {\rm M}_\odot$) Wolf-Rayet-like stars., Comment: 17 pages, 13 figures, 1 table, submitted to MNRAS

We present multiepoch spectropolarimetry of the superluminous interacting Type IIn supernova SN2017hcc, covering 16 to 391 days after explosion. In our first epoch we measure continuum polarization as high as 6%, making SN 2017hcc the most intrinsically polarized SN ever reported. During the first 29 days of coverage, when the polarization is strongest, the continuum polarization has a wavelength dependence that rises toward blue wavelengths, and becomes wavelength independent by day 45. The polarization strength drops rapidly during the first month, even as the SN flux is still climbing to peak brightness. Nonetheless, record-high polarization is maintained until day 68, at which point the source polarization declines to 1.9%, comparable to peak levels in previous well-studied SNe IIn. Thereafter the SN continues in polarization decline, while exhibiting only minor changes in position angle on the sky. The blue slope of the polarized continuum during the first month, accompanied by short-lived polarized flux for Balmer emission, suggests that an aspherical distribution of dust grains in pre-shock circumstellar material (CSM) is echoing the SN IIn spectrum and strongly influencing the polarization, while the subsequent decline during the wavelength-independent phase appears broadly consistent with electron scattering near the SN/CSM interface. The persistence of the polarization position angle between these two phases suggests that the pre-existing CSM responsible for the dust scattering at early times is part of the same geometric structure as the electron-scattering region that dominates the polarization at later times. SN2017hcc appears to be yet another, but much more extreme, case of aspherical yet well-ordered CSM in Type IIn SNe, possibly resulting from pre-SN mass loss shaped by a binary progenitor system., Comment: Submitted to MNRAS April 21, 2023

Superluminous supernovae are among the most energetic stellar explosions in the Universe, but their energy sources remain an open question. Here we present long-term observations of one of the closest examples of the hydrogen-poor subclass (SLSNe-I), SN~2017egm, revealing the most complicated known luminosity evolution of SLSNe-I. Three distinct post-peak bumps were recorded in its light curve collected at about $100$--350\,days after maximum brightness, challenging current popular power models such as magnetar, fallback accretion, and interaction between ejecta and a circumstellar shell. However, the complex light curve can be well modelled by successive interactions with multiple circumstellar shells with a total mass of about $6.8$--7.7\,M$_\odot$. In this scenario, large energy deposition from interaction-induced reverse shocks results in ionization of neutral oxygen in the supernova ejecta and hence a much lower nebular-phase line ratio of [O\,\textsc{i}] $\lambda6300$/([Ca\,\textsc{ii}] + [O\,\textsc{ii}]) $\lambda7300$ ($\sim 0.2$) compared with that derived for other superluminous and normal stripped-envelope SNe. The pre-existing multiple shells indicate that the progenitor of SN~2017egm experienced pulsational mass ejections triggered by pair instability within 2 years before explosion, in robust agreement with theoretical predictions for a pre-pulsation helium-core mass of 48--51\,M$_{\odot}$. Finally, this work shows that the final explosion product may be a black hole with about 40\,M$_{\odot}$, and has significant implication for the formation of such heavy black holes that have been recently observed by LIGO-Virgo gravitational wave detectors., Comment: 34 pages, 13 figures, 6 tables

We report early-time ultraviolet (UV) and optical spectroscopy of the young, nearby Type II supernova (SN) 2022wsp obtained by the Hubble Space Telescope (HST)/STIS at about 10 and 20 days after the explosion. The SN 2022wsp UV spectra are compared to those of other well-observed Type II/IIP SNe, including the recently studied Type IIP SN 2021yja. Both SNe exhibit rapid cooling and similar evolution during early phases, indicating a common behavior among SNe II. Radiative-transfer modeling of the spectra of SN 2022wsp with the TARDIS code indicates a steep radial density profile in the outer layer of the ejecta, a supersolar metallicity, and a relatively high total extinction of E(B-V) = 0.35 mag. The early-time evolution of the photospheric velocity and temperature derived from the modeling agree with the behavior observed from other previously studied cases. The strong suppression of hydrogen Balmer lines in the spectra suggests interaction with a pre-existing circumstellar environment could be occurring at early times. In the SN 2022wsp spectra, the absorption component of the Mg II P Cygni profile displays a double-trough feature on day +10 that disappears by day +20. The shape is well reproduced by the model without fine-tuning the parameters, suggesting that the secondary blueward dip is a metal transition that originates in the SN ejecta., Comment: Submitted to ApJ Letters on 4/11/2023

Over the past few years, wide-field time-domain surveys like ZTF and OGLE have led to discoveries of various types of interesting short-period stellar variables, such as ultracompact eclipsing binary white dwarfs, rapidly rotating magnetised white dwarfs (WDs), transitional cataclysmic variables between hydrogen-rich and helium accretion, and blue large-amplitude pulsators (BLAPs), which greatly enrich our understandings of stellar physics under some extreme conditions. In this paper, we report the first-two-year discoveries of short-period variables (i.e., P<2 hr) by the Tsinghua University-Ma Huateng Telescopes for Survey (TMTS). TMTS is a multi-tube telescope system with a field of view up to 18 deg^2, which started to monitor the LAMOST sky areas since 2020 and generated uninterrupted minute-cadence light curves for about ten million sources within 2 years. Adopting the Lomb-Scargle periodogram with period-dependent thresholds for the maximum powers, we identify over 1 100 sources that exhibit a variation period shorter than 2 hr. Compiling the light curves with the Gaia magnitudes and colours, LAMOST spectral parameters, VSX classifications, and archived observations from other prevailing time-domain survey missions, we identified 1 076 as delta Scuti stars, which allows us study their populations and physical properties in the short-period regime. The other 31 sources include BLAPs, subdwarf B variables (sdBVs), pulsating WDs, ultracompact/short-period eclipsing/ellipsoidal binaries, cataclysmic variables below the period gap, etc., which are highly interesting and worthy of follow-up investigations., Comment: 22 pages, 15 figures, 5 tables, accepted by MNRAS

We present six epochs of optical spectropolarimetry of the Type IIP supernova (SN) 2021yja ranging from $\sim$ 25 to 95 days after the explosion. An unusually high continuum linear polarization of $p \sim 0.9\%$ is measured during the early photospheric phase, followed by a steady decrease well before the onset of the nebular phase. This behavior has not been observed before in Type IIP supernovae (SNe IIP). The observed continuum polarization angle does not change significantly during the photospheric phase. We find a pronounced axis of symmetry in the global ejecta that is shared in common with the H$\alpha$ and Ca II near-infrared triplet lines. These observations are consistent with an ellipsoidal geometry. The temporal evolution of the continuum polarization is also compatible with the SN ejecta interacting with aspherical circumstellar matter, although no spectroscopic features that may be associated with strong interaction can be identified. Alternatively, we consider the source of the high polarization to be an extended hydrogen envelope that is indistinguishable from low-density circumstellar matter., Comment: Submitted to MNRAS

When discovered, SN~2017egm was the closest (redshift $z=0.03$) hydrogen-poor superluminous supernova (SLSN-I) and a rare case that exploded in a massive and metal-rich galaxy. Thus, it has since been extensively observed and studied. We report spectroscopic data showing strong emission at around He~I $\lambda$10,830 and four He~I absorption lines in the optical. Consequently, we classify SN~2017egm as a member of an emerging population of helium-rich SLSNe-I (i.e., SLSNe-Ib). We also present our late-time photometric observations. By combining them with archival data, we analyze high-cadence ultra-violet, optical, and near-infrared light curves spanning from early pre-peak ($\sim -20\,d$) to late phases ($\sim +300\,d$). We obtain its most complete bolometric light curve, in which multiple bumps are identified. None of the previously proposed models can satisfactorily explain all main light-curve features, while multiple interactions between the ejecta and circumstellar material (CSM) may explain the undulating features. The prominent infrared excess with a blackbody luminosity of $10^7$--$10^8\,L_{sun}$ detected in SN~2017egm could originate from the emission of either an echo of a pre-existing dust shell, or newly-formed dust, offering an additional piece of evidence supporting the ejecta-CSM interaction model. Moreover, our analysis of deep $Chandra$ observations yields the tightest-ever constraint on the X-ray emission of an SLSN-I, amounting to an X-ray-to-optical luminosity ratio $\lesssim 10^{-3}$ at late phases ($\sim100-200\,d$), which could help explore its close environment and central engine., Comment: 25 pages, 14 Figures, 4 Tables; accepted for publication in ApJ (Mar. 2023)

The AGN STORM 2 collaboration targeted the Seyfert 1 galaxy Mrk 817 for a year-long multiwavelength, coordinated reverberation mapping campaign including HST, Swift, XMM-Newton, NICER, and ground-based observatories. Early observations with NICER and XMM revealed an X-ray state ten times fainter than historical observations, consistent with the presence of a new dust-free, ionized obscurer. The following analysis of NICER spectra attributes variability in the observed X-ray flux to changes in both the column density of the obscurer by at least one order of magnitude ($N_\mathrm{H}$ ranges from $2.85\substack{+0.48\\ -0.33} \times 10^{22}\text{ cm}^{-2}$ to $25.6\substack{+3.0\\ -3.5} \times 10^{22} \text{ cm}^{-2}$) and the intrinsic continuum brightness (the unobscured flux ranges from $10^{-11.8}$ to $10^{-10.5}$ erg s$^{-1}$ cm$^{-2}$ ). While the X-ray flux generally remains in a faint state, there is one large flare during which Mrk 817 returns to its historical mean flux. The obscuring gas is still present at lower column density during the flare but it also becomes highly ionized, increasing its transparency. Correlation between the column density of the X-ray obscurer and the strength of UV broad absorption lines suggests that the X-ray and UV continua are both affected by the same obscuration, consistent with a clumpy disk wind launched from the inner broad line region., Comment: 19 pages, 9 figures, Accepted to ApJ

Type Iax supernovae (SN Iax) are the largest known class of peculiar white dwarf supernovae, distinct from normal Type Ia supernovae (SN Ia). The unique properties of SN Iax, especially their strong photospheric lines out to extremely late times, allow us to model their optical spectra and derive physical parameters for the long-lasting photosphere. We present an extensive spectral timeseries, including 21 new spectra, of SN Iax 2014dt from +11 to +562 days after maximum light. We are able to reproduce the entire timeseries with a self-consistent, nearly unaltered deflagration explosion model from Fink et al. (2014) using TARDIS, an open-source radiative transfer code (Kerzendorf & Sim 2014; Kerzendorf et al. 2023). We find that the photospheric velocity of SN 2014dt slows its evolution between +64 and +148 days, which closely overlaps the phase when we see SN 2014dt diverge from the normal spectral evolution of SN Ia (+90 to +150 days). The photospheric velocity at these epochs, ~400$-$1000 km s$^{-1}$, may demarcate a boundary within the ejecta below which the physics of SN Iax and normal SN Ia differ. Our results suggest that SN 2014dt is consistent with a weak deflagration explosion model that leaves behind a bound remnant and drives an optically thick, quasi-steady-state wind creating the photospheric lines at late times. The data also suggest that this wind may weaken at epochs past +450 days, perhaps indicating a radioactive power source that has decayed away., Comment: Accepted to ApJ, 22 pages, 8 figures, 3 tables

We report on analysis using the James Webb Space Telescope (JWST) to identify a candidate progenitor star of the Type II-plateau supernova SN 2022acko in the nearby, barred spiral galaxy NGC 1300. To our knowledge, our discovery represents the first time JWST has been used to localize a progenitor system in pre-explosion archival Hubble Space Telescope (HST) images. We astrometrically registered a JWST NIRCam image from 2023 January, in which the SN was serendipitously captured, to pre-SN HST F160W and F814W images from 2017 and 2004, respectively. An object corresponding precisely to the SN position has been isolated with reasonable confidence. That object has a spectral energy distribution (SED) and overall luminosity consistent with a single-star model having an initial mass possibly somewhat less than the canonical 8 Msun theoretical threshold for core collapse (although masses as high as 9 Msun for the star are also possible); however, the star's SED and luminosity are inconsistent with that of a super-asymptotic giant branch star which might be a forerunner of an electron-capture SN. The properties of the progenitor alone imply that SN 2022acko is a relatively normal SN II-P, albeit most likely a low-luminosity one. The progenitor candidate should be confirmed with follow-up HST imaging at late times, when the SN has sufficiently faded. This potential use of JWST opens a new era of identifying SN progenitor candidates at high spatial resolution., Comment: 8 pages, substantial changes from v1, to appear in MNRAS

Supernova (SN) explosions have been sought for decades as a possible source of dust in the Universe, providing the seeds of galaxies, stars, and planetary systems. SN 1987A offers one of the most promising examples of significant SN dust formation, but until the James Webb Space Telescope (JWST), instruments have traditionally lacked the sensitivity at both late times (>1 yr post-explosion) and longer wavelengths (i.e., >10 um) to detect analogous dust reservoirs. Here we present JWST/MIRI observations of two historic Type IIP SNe, 2004et and SN 2017eaw, at nearly 18 and 5 yr post-explosion, respectively. We fit the spectral energy distributions as functions of dust mass and temperature, from which we are able to constrain the dust geometry, origin, and heating mechanism. We place a 90% confidence lower limit on the dust masses for SNe 2004et and 2017eaw of >0.014 and >4e-4 M_sun, respectively. More dust may exist at even colder temperatures or may be obscured by high optical depths. We conclude dust formation in the ejecta to be the most plausible and consistent scenario. The observed dust is radiatively heated to ~100-150 K by ongoing shock interaction with the circumstellar medium. Regardless of the best fit or heating mechanism adopted, the inferred dust mass for SN 2004et is the second highest (next to SN 1987A) inferred dust mass in extragalactic SNe thus far, promoting the prospect of SNe as potential significant sources of dust in the Universe., Comment: 12 pages, 7 figures, submitting to MNRAS

We present extensive optical/ultraviolet observations and modelling analysis for the nearby SN 1987A-like peculiar Type II supernova (SN) 2018hna. Both photometry and spectroscopy covered phases extending to $>$500 days after the explosion, making it one of the best-observed SN II of this subtype. SN 2018hna is obviously bluer than SN 1987A during the photospheric phase, suggesting higher photospheric temperature, which may account for weaker BaII $\mathrm{\lambda}$6142 lines in its spectra. Analysis of early-time temperature evolution suggests a radius of $\sim$45 $\mathrm{R_{\odot}}$ for the progenitor of SN 2018hna, consistent with a blue supergiant (BSG). By fitting the bolometric light curve with hydrodynamical models, we find that SN 2018hna has an ejecta mass of $\sim$(13.7--17.7) $\mathrm{M_{\odot}}$, a kinetic energy of $\sim$ (1.0--1.2) $\times 10^{51}$ erg, and a $^{56}$Ni mass of about 0.05 $\mathrm{M_{\odot}}$. Moreover, based on standard stellar evolution and the oxygen mass (0.44--0.73 $\mathrm{M_{\odot}}$) deduced from nebular [OI] lines, the progenitor of SN 2018hna is expected to have an initial main-sequence mass $<$16 $\mathrm{M_{\odot}}$. In principle, such a relatively low-mass star cannot end as a BSG just before core-collapse, except some unique mechanisms are involved, such as rapid rotation, restricted semiconvection, etc. On the other hand, binary scenario may be more favourable, like in the case of SN 1987A. While the much lower oxygen mass inferred for SN~2018hna may imply that its progenitor system also had much lower initial masses than that of SN 1987A., Comment: 18 pages; accepted for publication in MNRAS

Much excitement surrounds the intense mass loss that seems to take place in some massive stars immediately before core collapse. However, occurring too late, it has a negligible impact on the star's evolution or the final yields, which are influenced instead by the longer-term, quasi-steady mass loss taking place during H and He burning. Late-time observations of core-collapse supernovae interacting with the progenitor wind are one means to constrain this secular mass loss. Here, we present radiative transfer calculations for a Type II SN with and without this interaction power, focusing on the phase between 350 and 1000d after explosion. Without interaction power, the ejecta are powered through radioactive decay whose exponential decline produces an ever-fading SN. Instead, with a constant interaction power of $10^{40}$ erg s$^{-1}$ (representative of an SN II ramming into a steady-state $10^{-6} M_\odot$ yr$^{-1}$ wind), the spectrum morphs from decay powered at 350d, with narrow lines forming in the inner metal-rich ejecta, to interaction powered at 1000d, with broad boxy lines forming in the outer H-rich ejecta. Intermediate times are characterized by a hybrid and complex spectrum made of overlapping narrow and broad lines. While interaction boosts primarily the flux in the ultraviolet, which remains largely unobserved today, a knee in the $R$-band light curve or a $U$-band boost are clear signatures of interaction at late times. The model predictions compare favorably with a number of Type II supernovae including SN 2004et or SN 2017eaw at 500-1000d after explosion., Comment: submitted to A&A

Optical spectropolarimetry of the normal thermonuclear supernova SN2019np from -14.5 to +14.5 days relative to B-band maximum detected an intrinsic continuum polarization, p(cont), of 0.21+-0.09% at the first epoch. Between days -11.5 to +05, p(cont) remained about 0 and by day +14.5 was again significant at 0.19+-0.10%. Not considering the first epoch, the dominant axis of SiII(6355A) was roughly constant, staying close to the continuum until both rotated in opposite directions on day +14.5. Detailed radiation-hydrodynamical simulations produce a very steep density slope in the outermost ejecta so that the low first-epoch p(const) of about 0.2% nevertheless suggests a separate structure with an axis ratio of about 2 in the outer carbon-rich 3...5E-3 Mo. Large-amplitude fluctuations in the polarization profiles and a flocculent appearance of the polar diagram for the CaII near-infrared triplet (NIR3) may be related by a common origin. The temporal evolution of the polarization spectra agrees with an off-center delayed detonation. The late-time increase in polarization and the possible change in position angle are also consistent with an aspherical 56Ni core. The p(cont) and the absorptions due to Si II(6355A) and the CaII NIR3 form in the same region of the extended photosphere, with an interplay between line occultation and thermalization producing p. Small-scale polarization features may be due to small-scale structures, but many could be related to atomic patterns of the quasi-continuum; the hardly have an equivalent in the total-flux spectrum. We compare SN2019np to other SNe and develop future objectives and strategies for SNIa spectropolarimetry., Comment: 23 pages, 18 figures, 2 tables, MNRAS, submitted 11/9/22, revised 12/9/22

As part of a larger completed Hubble Space Telescope (HST) Snapshot program, we observed the sites of six nearby core-collapse supernovae (SNe) at high spatial resolution: SN 2012A, SN 2013ej, SN 2016gkg, SN 2017eaw, SN 2018zd, and SN 2018aoq. These observations were all conducted at sufficiently late times in each SN's evolution to demonstrate that the massive-star progenitor candidate identified in each case in pre-explosion imaging data had indeed vanished and was therefore most likely the actual progenitor. However, we have determined for SN 2016gkg that the progenitor candidate was most likely a blend of two objects: the progenitor, which itself has likely vanished, and another closely-neighbouring star. We thus provide a revised estimate of that progenitor's properties: a binary system with a hydrogen-stripped primary star at explosion with effective temperature ~6300--7900 K, bolometric luminosity ~10^{4.65} L_sun, radius ~118--154 R_sun, and initial mass 9.5--11 M_sun. Utilising late-time additional archival HST data nearly contemporaneous with our Snapshots, we also show that SN 2017eaw had a luminous ultraviolet excess, which is best explained as a result of ongoing interaction of the SN shock with pre-existing circumstellar matter. We offer the caveat, particularly in the case of SN 2013ej, that obscuration from SN dust may be compromising our conclusions. This sample adds to the growing list of confirmed or likely core-collapse SN progenitors., Comment: 12 pages, to appear in MNRAS

We analyse photometric observations of the supernova (SN) impostor SN 2000ch in NGC 3432 covering the time since its discovery. This source was previously observed to have four outbursts in 2000-2010. Observations now reveal at least three additional outbursts in 2004-2007, and sixteen outbursts in 2010-2022. Outburst light curves are irregular and multipeaked, exhibiting a wide variety of peak magnitude, duration, and shape. The outbursts after 2008 repeat with a period of 200.7$\pm{2}$ d, while the outburst in 2000 seems to match with a shorter period. The next outburst should occur around January/February 2023. We propose that these periodic eruptions arise from violent interaction around times of periastron in an eccentric binary system, similar to the periastron encounters of $\eta$ Carinae leading up to its Great Eruption, and resembling the erratic pre-SN eruptions of SN 2009ip. We attribute the irregularity of the eruptions to the interplay between the orbit and the variability of the luminous blue variable (LBV) primary star, wherein each successive periastron pass may have a different intensity or duration due to the changing radius and mass-loss rate of the LBV-like primary. Such outbursts may occasionally be weak or undetectable if the LBV is relatively quiescent at periastron, but can be much more extreme when the LBV is active. The observed change in orbital period may be a consequence of mass lost in outbursts. Given the similarity to the progenitor of SN 2009ip, SN 2000ch deserves continued attention in the event it is headed for a stellar merger or a SN-like explosion., Comment: 18 pages, 13 figures, 4 tables, MNRAS accepted

Luminous red novae (LRNe) are transients characterized by low luminosities and expansion velocities, and are associated with mergers or common envelope ejections in stellar binaries. Intermediate-luminosity red transients (ILRTs) are an observationally similar class with unknown origins, but generally believed to either be electron capture supernovae (ECSN) in super-AGB stars, or outbursts in dusty luminous blue variables (LBVs). In this paper, we present a systematic sample of 8 LRNe and 8 ILRTs detected as part of the Census of the Local Universe (CLU) experiment on the Zwicky Transient Facility (ZTF). The CLU experiment spectroscopically classifies ZTF transients associated with nearby ($<150$ Mpc) galaxies, achieving 80% completeness for m$_{r}<20$\,mag. Using the ZTF-CLU sample, we derive the first systematic LRNe volumetric-rate of 7.8$^{+6.5}_{-3.7}\times10^{-5}$ Mpc$^{-3}$ yr$^{-1}$ in the luminosity range $-16\leq$M$_{\rm{r}}$$\leq -11$ mag. We find that in this luminosity range, the LRN rate scales as dN/dL $\propto L^{-2.5\pm0.3}$ - significantly steeper than the previously derived scaling of $L^{-1.4\pm0.3}$ for lower luminosity LRNe (M$_{V}\geq-10$). The steeper power law for LRNe at high luminosities is consistent with the massive merger rates predicted by binary population synthesis models. We find that the rates of the brightest LRNe (M$_{r}\leq-13$ mag) are consistent with a significant fraction of them being progenitors of double compact objects (DCOs) that merge within a Hubble time. For ILRTs, we derive a volumetric rate of $2.6^{+1.8}_{-1.4}\times10^{-6}$ Mpc$^{-3}$yr$^{-1}$ for M$_{\rm{r}}\leq-13.5$, that scales as dN/dL $\propto L^{-2.5\pm0.5}$. This rate is $\approx1-5\%$ of the local core-collapse supernova rate, and is consistent with theoretical ECSN rate estimates., Comment: 32 pages, 16 figures. Submitted to ApJ, comments welcome