Search

Stellar mass black holes are formed from the terminal collapse of massive stars if the ensuing neutrino shock is unable to eject the stellar envelope. Direct observations of black hole formation remain inconclusive. We report observations of M31-2014-DS1, a massive, hydrogen-depleted supergiant in the Andromeda galaxy identified via a mid-infrared brightening in 2014. Its total luminosity remained nearly constant for the subsequent thousand days, before fading dramatically over the next thousand days by $\gtrsim 10\times$ and $\gtrsim 10^4\times$ in total and visible light, respectively. Together with the lack of a detected optical outburst, the observations are explained by the fallback of the stellar envelope into a newly formed black hole, moderated by the injection of a $\sim 10^{48}$ erg shock. Unifying these observations with a candidate in NGC 6946, we present a concordant picture for the birth of stellar mass black holes from stripped massive stars., Comment: Submitted for review

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!

The mechanism of X-ray outbursts in Be X-ray binaries remains a mystery, and understanding their circumstellar disks is crucial for a solution of the mass-transfer problem. In particular, it is important to identify the Be star activities (e.g., pulsations) that cause mass ejection and, hence, disk formation. Therefore, we investigated the relationship between optical flux oscillations and the infrared (IR) excess in a sample of five Be X-ray binaries. Applying the Lomb-Scargle technique to high-cadence optical light curves from the Transiting Exoplanet Survey Satellite (TESS), we detected several significant oscillation modes in the 3 to 24 hour period range for each source. We also measured the IR excess (a proxy for disk growth) of those five sources, using J-band light curves from Palomar Gattini-IR. In four of the five sources, we found anti-correlations between the IR excess and the amplitude of the main flux oscillation modes. This result is inconsistent with the conventional idea that non-radial pulsations drive mass ejections. We propose an alternative scenario where internal temperature variations in the Be star cause transitions between pulsation-active and mass-ejection-active states., Comment: 17 pages, 27 figures, 6 tables, accepted for publication in MNRAS

We study a magnitude-limited sample of 36 Broad-lined Type Ic Supernovae (SNe Ic-BL) from the Zwicky Transient Facility Bright Transient Survey (detected between March 2018 and August 2021), which is the largest systematic study of SNe Ic-BL done in literature thus far. We present the light curves (LCs) for each of the SNe, and analyze the shape of the LCs to derive empirical parameters, along with the explosion epochs for every event. The sample has an average absolute peak magnitude in the r band of $M_r^{max}$ = -18.51 $\pm$ 0.15 mag. Using spectra obtained around peak light, we compute expansion velocities from the Fe II 5169 Angstrom line for each event with high enough signal-to-noise ratio spectra, and find an average value of $v_{ph}$ = 16,100 $\pm$ 1,100 km $s^{-1}$. We also compute bolometric LCs, study the blackbody temperature and radii evolution over time, and derive the explosion properties of the SNe. The explosion properties of the sample have average values of $M_{Ni}$ = $0.37_{-0.06}^{+0.08}$ solar masses, $M_{ej}$ = $2.45_{-0.41}^{+0.47}$ solar masses, and $E_K$= $4.02_{-1.00}^{+1.37} \times 10^{51}$ erg. Thirteen events have radio observations from the Very Large Array, with 8 detections and 5 non-detections. We find that the populations that have radio detections and radio non-detections are indistinct from one another with respect to their optically-inferred explosion properties, and there are no statistically significant correlations present between the events' radio luminosities and optically-inferred explosion properties. This provides evidence that the explosion properties derived from optical data alone cannot give inferences about the radio properties of SNe Ic-BL, and likely their relativistic jet formation mechanisms., Comment: 52 pages, 34 Figures, 8 Tables; Accepted to ApJ, Revised Title from Proofs

The James Webb Space Telescope (JWST) has opened up a new window to study highly reddened explosive transients. We present results from late-time (1421 days post-explosion) JWST follow-up spectroscopic observations with NIRSpec and MIRI LRS of the intermediate luminosity red transient (ILRT) AT 2019abn located in the nearby Messier 51 galaxy (8.6 Mpc). ILRTs represent a mysterious class of transients which exhibit peak luminosities between those of classical novae and supernovae and which are known to be highly dust obscured. Similar to the prototypical examples of this class of objects, NGC 300 2008-OT and SN 2008S, AT 2019abn has an extremely red and dusty progenitor detected only in pre-explosion Spitzer/IRAC imaging at 3.6 and 4.5 micron and not in deep optical or near-infrared HST images. We find that late time observations of AT 2019abn from NEOWISE and JWST are consistent with the late time evolution of SN 2008S. In part because they are so obscured by dust, it is unknown what produces an ILRT with hypotheses ranging from high mass stellar merger events, non-terminal stellar outbursts, or terminal supernovae explosions through electron-capture in super-AGB stars. Our JWST observations show strong mid-IR Class C PAH features at 6.3 and 8.25 micron typical of carbon-rich post-AGB sources. These features suggest the dust around AT 2019abn, either pre-existing or newly formed in the ejecta, is composed of carbonaceous grains which are not typically observed around red supergiants. However, depending on the strength and temperature of hot bottom burning, SAGBs may be expected to exhibit a carbon-rich chemistry. Thus our JWST observations are consistent with AT 2019abn having an SAGB progenitor., Comment: 12 pages, 4 figures, submitted to ApJL

We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. R Coronae Borealis stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that surveys 18000 deg$^{2}$ of the northern sky ($\delta>-28^{o}$) at a cadence of 2 days. Using PGIR J-band lightcurves for $\sim$60 million stars together with mid-IR colors from WISE, we selected a sample of 530 candidate RCB stars. We obtained near-IR spectra for these candidates and identified 53 RCB stars in our sample. Accounting for our selection criteria, we find that there are a total of $\approx350^{+150}_{-100}$ RCB stars in the Milky Way. Assuming typical RCB lifetimes, this corresponds to an RCB formation rate of 0.8 - 5 $\times$ 10$^{-3}$ yr$^{-1}$, consistent with observational and theoretical estimates of the He-CO WD merger rate. We searched for quasi-periodic pulsations in the PGIR lightcurves of RCB stars and present pulsation periods for 16 RCB stars. We also examined high-cadenced TESS lightcurves for RCB and the chemically similar, but dustless hydrogen-deficient carbon (dLHdC) stars. We find that dLHdC stars show variations on timescales shorter than RCB stars, suggesting that they may have lower masses than RCB stars. Finally, we identified 3 new spectroscopically confirmed and 12 candidate Galactic DY Per type stars - believed to be colder cousins of RCB stars - doubling the sample of Galactic DY Per type stars., Comment: accepted for publication in PASP

Roughly 2% of white dwarfs harbor planetary debris disks detectable via infrared excesses, but only a few percent of these disks show a gaseous component, distinguished by their double-peaked emission at the near-infrared calcium triplet. Previous studies found most debris disks around white dwarfs are variable at 3.4 and 4.5 $\mu$m, but they analyzed only a few of the now 21 published disks showing calcium emission. To test if most published calcium emission disks exhibit large-amplitude stochastic variability in the near-infrared, we use light curves generated from the unWISE images at 3.4 $\mu$m that are corrected for proper motion to characterize the near-infrared variability of these disks against samples of disks without calcium emission, highly variable cataclysmic variables, and 3215 isolated white dwarfs. We find most calcium emission disks are extremely variable: 6/11 with sufficient signal-to-noise show high-amplitude variability in their 3.4-$\mu$m light curves. These results lend further credence to the notion that disks showing gaseous debris in emission are the most collisionally active. Under the assumption that 3.4-$\mu$m variability is characteristic of white dwarfs with dusty debris disks, we generate a catalog of 104 high-confidence near-infrared variable white dwarfs, 84 of which are published as variable for the first time. We do near-infrared spectroscopic follow-up of seven new candidate 3.4-$\mu$m variables, confirming at least one new remnant planetary system, and posit that empirical near-infrared variability can be a discovery engine for debris disks showing gaseous emission., Comment: 27 pages, 8 figures, accepted to The Astrophysical Journal. Supplemental catalogs and light curves available at https://zenodo.org/doi/10.5281/zenodo.12538705

The mid-infrared (MIR; $\lambda\simeq3 - 10\mu$m) bands offer a unique window into understanding accretion and its interplay with jet formation in Galactic black hole X-ray binaries (BHXRBs). Although extremely difficult to observe from the ground, the NEOWISE time domain survey offers an excellent data set to study MIR variability when combined with contemporaneous X-ray data from the MAXI all-sky survey over a $\approx15$ yr baseline. Using a new forced photometry pipeline for NEOWISE data, we present the first systematic study of BHXRB MIR variability in outburst. Analyzing a sample of 16 sources detected in NEOWISE, we show variability trends in the X-ray hardness and MIR spectral index wherein i) the MIR bands are typically dominated by jet emission during the hard states, constraining the electron power spectrum index to $p \approx 1-4$ in the optically thin regime and indicating emitting regions of a few tens of gravitational radii when evolving towards a flat spectrum, ii) the MIR luminosity ($L_{IR}$) scales as $L_{IR}\propto L_X^{0.82\pm0.12}$ with the $2-10$ keV X-ray luminosity ($L_X$) in the hard state, consistent with its origin in a jet, and iii) the thermal disk emission dominates the soft state as the jet switches off and dramatically suppresses ($\gtrsim 10\times$) the MIR emission into a inverted spectrum ($\alpha\approx -1$, where $F_\nu\propto\nu^{-\alpha}$). We highlight a population of `mini' BHXRB outbursts detected in NEOWISE (including two previously unreported episodes in MAXI J1828-249) but missed in MAXI due to their faint fluxes or source confusion, exhibiting MIR spectral indices suggestive of thermal emission from a large outer disk. We highlight that upcoming IR surveys and the Rubin observatory will be powerful discovery engines for the distinctively large amplitude and long-lived outbursts of BHXRBs, as an independent discovery route to X-ray monitors., Comment: 11 pages, 5 figures, 1 table, submitted to MNRAS. Comments welcome!

Palomar Gattini-IR is a wide-field, synoptic infrared time domain survey covering $\approx 15000$ sq. deg. of the sky at $\approx 1-3$ night cadence to a depth of $J\approx 13.0$ and $\approx 14.9$ Vega mag in and outside the Galactic plane, respectively. Here, we present the first data release of $J$-band light curves of 2MASS sources within the survey footprint covering approximately the first four years of operations. We describe the construction of the source catalog based on the 2MASS point sources, followed by exposure filtering criteria and forced PSF photometry. The catalog contains light curves of $\approx 286$ million unique sources with 2MASS magnitudes of $J < 15.5$ mag, with a total of $\approx 50$ billion photometric measurements and $\approx 20$ billion individual source detections at signal-to-noise-ratio $> 3$. We demonstrate the photometric fidelity of the catalog by i) quantifying the magnitude-dependent accuracy and uncertainty of the photometry with respect to 2MASS and ii) comparing against forced PGIR aperture photometry for known variable sources. We present simple filtering criteria for selecting reliable photometric measurements as well as example Python notebooks for users. This catalog is the largest compilation of nightly cadence, synoptic infrared light curves to date, comparable to those in the largest optical surveys, providing a stepping stone to upcoming infrared surveys in the coming decade., Comment: 10 Pages, 5 figures, submitted to PASP. Full catalog will be available upon publication, individual requests welcome!

While core-collapse supernovae (SNe) often show early and consistent signs of circumstellar (CSM) interaction, some exhibit delayed signatures due to interaction with distant material around the progenitor star. Here we present the discovery in NEOWISE data of WTP19aalnxx, a luminous mid-infrared (IR) transient in the outskirts of the galaxy KUG 0022-007 at $\approx 190$ Mpc. First detected in 2018, WTP19aalnxx reaches a peak absolute (Vega) magnitude of $\approx-22$ at $4.6 \, \mu$m in $\approx3$ yr, comparable to the most luminous interacting SNe. Archival data reveal a $\gtrsim 5\times$ fainter optical counterpart detected since 2015, while follow-up near-IR observations in 2022 reveal an extremely red ($Ks-W2 \approx 3.7$ mag) active transient. Deep optical spectroscopy confirm strong CSM interaction signatures via intermediate-width Balmer emission lines and coronal metal lines. Modeling the broadband spectral energy distribution, we estimate the presence of $\gtrsim 10^{-2}$ M$_\odot$ of warm dust, likely formed in the shock interaction region. Together with the lack of nebular Fe emission, we suggest that WTP19aalnxx is a missed, low (optical) luminosity SN in an emerging family of core-collapse SNe distinguished by their CSM-interaction-powered mid-IR emission that outshines the optical bands. Investigating the Zwicky Transient Facility sample of SNe in NEOWISE data, we find $17$ core-collapse SNe ($\gtrsim 3$% in a volume-limited sample) without early signs of CSM interaction that exhibit delayed IR brightening, suggestive of dense CSM shells at $\lesssim 10^{17}$cm. We suggest that synoptic IR surveys offer a new route to revealing late-time CSM interaction and the prevalence of intense terminal mass loss in massive stars., Comment: 15 pages, 5 figures, submitted to ApJ

We present SN 2023zaw $-$ a sub-luminous ($\mathrm{M_r} = -16.7$ mag) and rapidly-evolving supernova ($\mathrm{t_{1/2,r}} = 4.9$ days), with the lowest nickel mass ($\approx0.002$ $\mathrm{M_\odot}$) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He I and Ca NIR emission lines with velocities of $\sim10\ 000 - 12\ 000$ $\mathrm{km\ s^{-1}}$. The late-time spectra show prominent narrow He I emission lines at $\sim$1000$\ \mathrm{km\ s^{-1}}$, indicative of interaction with He-rich circumstellar material. SN 2023zaw is located in the spiral arm of a star-forming galaxy. We perform radiation-hydrodynamical and analytical modeling of the lightcurve by fitting with a combination of shock-cooling emission and nickel decay. The progenitor has a best-fit envelope mass of $\approx0.2$ $\mathrm{M_\odot}$ and an envelope radius of $\approx50$ $\mathrm{R_\odot}$. The extremely low nickel mass and low ejecta mass ($\approx0.5$ $\mathrm{M_\odot}$) suggest an ultra-stripped SN, which originates from a mass-losing low mass He-star (ZAMS mass $<$ 10 $\mathrm{M_\odot}$) in a close binary system. This is a channel to form double neutron star systems, whose merger is detectable with LIGO. SN 2023zaw underscores the existence of a previously undiscovered population of extremely low nickel mass ($< 0.005$ $\mathrm{M_\odot}$) stripped-envelope supernovae, which can be explored with deep and high-cadence transient surveys., Comment: Accepted for publication in ApJL

Stars in the Asymptotic Giant Branch (AGB) phase, dominated by low to intermediate-mass stars in the late stage of evolution, undergo periodic pulsations, with periods of several hundred days, earning them the name Long Period Variables (LPVs). These stars gradually shed their mass through stellar winds and mass ejections, enveloping themselves in dust. Infrared (IR) surveys can probe these dust-enshrouded phases and uncover populations of LPV stars in the Milky Way. In this paper, we present a catalog of 159,696 Long Period Variables using near-IR lightcurves from the Palomar Gattini - IR (PGIR) survey. PGIR has been surveying the entire accessible northern sky ($\delta > -28^{\circ}$) in the J-band at a cadence of 2-3 days since September 2018, and has produced J-band lightcurves for more than 60 million sources. We used a gradient-boosted decision tree classifier trained on a comprehensive feature set extracted from PGIR lightcurves to search for LPVs in this dataset. We developed a parallelized and optimized code to extract features at a rate of ~0.1 seconds per lightcurve. Our model can successfully distinguish LPVs from other stars with a true positive rate and weighted g-mean of 0.95. 73,346 (~46%) of the sources in our catalog are new, previously unknown LPVs., Comment: 15 pages, 15 figures

Most tidal disruption events (TDEs) are currently found in time-domain optical and soft X-ray surveys, both of which are prone to significant obscuration. The infrared (IR), however, is a powerful probe of dust-enshrouded environments, and hence, we recently performed a systematic search of NEOWISE mid-IR data for nearby, obscured TDEs within roughly 200 Mpc. We identified 18 TDE candidates in galactic nuclei, using difference imaging to uncover nuclear variability amongst significant host galaxy emission. These candidates were selected based on the following IR light curve properties: (1) $L_\mathrm{W2}\gtrsim10^{42}$ erg s$^{-1}$ at peak, (2) fast rise, followed by a slow, monotonic decline, (3) no significant prior variability, and (4) no evidence for AGN activity in WISE colors. The majority of these sources showed no variable optical counterpart, suggesting that optical surveys indeed miss numerous obscured TDEs. Using narrow line ionization levels and variability arguments, we identified 6 sources as possible underlying AGN, yielding a total of 12 TDEs in our gold sample. This gold sample yields a lower limit on the IR-selected TDE rate of $2.0\pm0.3\times10^{-5}$ galaxy$^{-1}$ year$^{-1}$ ($1.3\pm0.2\times10^{-7}$ Mpc$^{-3}$ year$^{-1}$), which is comparable to optical and X-ray TDE rates. The IR-selected TDE host galaxies do not show a green valley overdensity nor a preference for quiescent, Balmer strong galaxies, which are both overrepresented in optical and X-ray TDE samples. This IR-selected sample represents a new population of dusty TDEs that have historically been missed by optical and X-ray surveys and helps alleviate tensions between observed and theoretical TDE rates and the so-called missing energy problem., Comment: 28 pages, 13 figures (+ 7 pages, 4 figures in appendix), accepted for publication in ApJ

We present the results of a high-time resolution polarization study of single pulses from pulsars B0950+08 and B1642-03. Single pulses from pulsar B0950+08 sometimes show isolated micropulses without any significant associated subpulse emission. Assuming that the properties of such micropulses represent the intrinsic nature of micropulse emission, we characterize the width and polarization properties of these `intrinsic' microstructures. Most of the `intrinsic' micropulses (~90%) follow common characteristic polarization properties, while the average width of these micropulses is consistent with the general micropulse population from this pulsar. Single pulses from these pulsars show a diverse range of polarization properties, including depolarization and mixing of two orthogonal modes resulting in polarization position angle jumps. We present a superposition model of the two orthogonal modes which can explain depolarization, the observed position angle jumps, and associated changes in other polarization parameters., Comment: 12 pages, 18 figures, Accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS)

We present evidence in support of the hypothesis that the young stellar object RNO 54 is a mature-stage FU Ori type source. The star was first cataloged as a ``red nebulous object" in the 1980s but appears to have undergone its outburst prior to the 1890s. Present-day optical and near-infrared spectra are consistent with those of other FU Ori type stars, both in the details of spectral line presence and shape, and in the overall change in spectral type from an FGK-type in the optical, to the M-type presented in the near-infrared. In addition, the spectral energy distribution of RNO 54 is well-fit by a pure-accretion disk model with parameters: $\dot{M} = 10^{-3.45\pm0.06}$ $M_\odot$ yr$^{-1}$, $M_* = 0.23\pm0.06 \ M_\odot$, and $R_\mathrm{inner} = 3.68\pm0.76 \ R_\odot$, though we believe $R_\mathrm{inner}$ is likely close to its upper range of $4.5 R_\odot$ in order to produce a $T_\mathrm{max} = 7000$ K that is consistent with the optical to near-infrared spectra. The resulting $L_\mathrm{acc}$ is $\sim 265 \ L_\odot$. To find these values, we adopted a source distance $d=1400$ pc and extinction $A_V=3.9$ mag, along with disk inclination $i=50$ deg based on consideration of confidence intervals from our initial disk model, and in agreement with observational constraints. The new appreciation of a well-known source as an FU Ori type object suggests that other such examples may be lurking in extant samples., Comment: to appear in AAS Journals

Large-amplitude accretion outbursts in young stars are expected to play a central role in proto-stellar assembly. Outburst identification historically has taken place using optical techniques, but recent, systematic infrared searches are enabling their discovery in heavily dust-obscured regions of the Galactic plane. Here, we present the discovery of WTP$\,$10aaauow, a large-amplitude mid-infrared (MIR) outburst identified in a systematic search of NEOWISE data using new image subtraction techniques. The source is located towards the RCW$\,$49 star-forming region, and estimated to be at a distance of $\approx 4\,$kpc via Gaia parallax measurement. Concurrent with the MIR brightening, the source underwent a $\gtrsim5\,$mag outburst in the optical and near-infrared (NIR) bands, reaching a peak luminosity of $\approx260\,$L$_\odot$ in 2014-2015, followed by a slow decline over the next 7 years. Analysis of the pre- and post-outburst spectral energy distributions reveal a pre-outburst stellar photosphere at a temperature of $3600-4000\,$K, surrounded by a likely two-component dust structure similar to a flat-spectrum or Class I type YSO. We present optical and NIR spectroscopy that show a GK-type spectrum in the optical bands exhibiting complex line profiles in strong absorption features, and evidence for a wind reaching a terminal velocity of $\approx 400\,$km$\,$s$^{-1}$. The NIR bands are characterized by a cooler M-type spectrum exhibiting a forest of atomic and molecular features. All together, the spectra demonstrate that WTP$\,$10aaauow is an FU Ori type outburst. Ongoing systematic infrared searches will continue to reveal the extent of this population in the Galactic disk., Comment: 9 pages, 5 figures, 2 tables

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

High mass X-ray binaries hosting red supergiant (RSG) donors are a rare but crucial phase in massive stellar evolution, with only one source previously known in the Milky Way. In this letter, we present the identification of the second Galactic RSG X-ray binary SWIFT J0850.8-4219. We identify the source 2MASS 08504008-4211514 as the likely infrared counterpart with a chance coincidence probability $\approx 5 \times 10^{-6}$. We present a $1.0 - 2.5\,\mu$m spectrum of the counterpart, exhibiting features characteristic of late-type stars and an exceptionally strong He I emission line, corroborating the identification. Based on i) the strength of the $^{12}$CO(2,0) band, ii) strong CN bandheads and absent TiO bandheads at $\approx 1.1\,\mu$m and iii) equivalent width of the Mg I $1.71\,\mu$m line, we classify the counterpart to be a K3$-$K5 type RSG with an effective temperature of $3820 \pm 100$ K, located at a distance of $\approx 12$ kpc. We estimate the source X-ray luminosity to be $(4 \pm 1) \times 10^{35}$ erg s$^{-1}$, with a hard photon index ($\Gamma < 1$), arguing against a white dwarf accretor but consistent with a magnetized neutron star in the propeller phase. Our results highlight the potential of systematic NIR spectroscopy of Galactic hard X-ray sources in completing our census of the local X-ray binary population., Comment: 6 pages, 4 figures, 1 table. Submitted to MNRAS letters. Comments welcome!

We use the light curve and spectral synthesis code JEKYLL to calculate a set of macroscopically mixed Type IIb supernova (SN) models, which are compared to both previously published and new late-phase observations of SN 2020acat. The models differ in the initial mass, the radial mixing and expansion of the radioactive material, and the properties of the hydrogen envelope. The best match to the photospheric and nebular spectra and lightcurves of SN 2020acat is found for a model with an initial mass of 17 solar masses, strong radial mixing and expansion of the radioactive material, and a 0.1 solar mass hydrogen envelope with a low hydrogen mass-fraction of 0.27. The most interesting result is that strong expansion of the clumps containing radioactive material seems to be required to fit the observations of SN 2020acat both in the diffusion phase and the nebular phase. These "Ni bubbles" are expected to expand due to heating from radioactive decays, but the degree of expansion is poorly constrained. Without strong expansion there is a tension between the diffusion phase and the subsequent evolution, and models that fit the nebular phase produce a diffusion peak that is too broad. The diffusion phase lightcurve is sensitive to the expansion of the "Ni bubbles", as the resulting Swiss-cheese-like geometry decreases the effective opacity and therefore the diffusion time. This effect has not been taken into account in previous lightcurve modelling of stripped-envelope SNe, which may lead to a systematic underestimate of their ejecta masses. It should be emphasized, though, that JEKYLL is limited to a geometry that is spherically symmetric on average, and large-scale asymmetries may also play a role. The relatively high initial mass found for the progenitor of SN 2020acat places it at the upper end of the mass distribution of Type IIb SN progenitors, and a single star origin can not be excluded., Comment: Accepted for publication by Astronomy and Astrophysics

Kilonovae are a rare class of astrophysical transients powered by the radioactive decay of nuclei heavier than iron, synthesized in the merger of two compact objects. Over the first few days, the kilonova evolution is dominated by a large number of radioactive isotopes contributing to the heating rate. On timescales of weeks to months, its behavior is predicted to differ depending on the ejecta composition and merger remnant. However, late-time observations of known kilonovae are either missing or limited. Here we report observations of a luminous red transient with a quasi-thermal spectrum, following an unusual gamma-ray burst of long duration. We classify this thermal emission as a kilonova and track its evolution up to two months after the burst. At these late times, the recession of the photospheric radius and the rapidly-decaying bolometric luminosity ($L_{\rm bol}\propto t^{-2.7\pm 0.4}$) support the recombination of lanthanide-rich ejecta as they cool., Comment: 47 pages, 14 figures, 9 tables; submitted; a minor typo fixed

Gaia21bty, a pre-main sequence star that previously had shown aperiodic dips in its light curve, underwent a considerable $\Delta G\approx2.9$ mag brightening that occurred over a few months between 2020 October - 2021 February. The Gaia lightcurve shows that the star remained near maximum brightness for about $4-6$ months, and then started slowly fading over the next 2 years, with at least three superimposed $\sim$1 mag sudden rebrightening events. Whereas the amplitude and duration of the maximum is typical for EXors, optical and near-infrared spectra obtained at the maximum are dominated by features which are typical for FUors. Modelling of the accretion disc at the maximum indicates that the disc bolometric luminosity is 43 L$_{\odot}$ and the mass accretion rate is $2.5\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, which are typical values for FUors even considering the large uncertainty in the distance ($1.7_{-0.4}^{+0.8}$ kpc). Further monitoring is necessary to understand the cause of the quick brightness decline, the rebrightening, and the other post-outburst light changes, as our multi-colour photometric data suggest that they could be caused by a long and discontinuous obscuration event. We speculate that the outburst might have induced large-scale inhomogeneous dust condensations in the line of sight leading to such phenomena, whilst the FUor outburst continues behind the opaque screen., Comment: Accepted to MNRAS

We present the discovery of the Type II supernova SN 2023ixf in M101 and follow-up photometric and spectroscopic observations, respectively, in the first month and week of its evolution. Our discovery was made within a day of estimated first light, and the following light curve is characterized by a rapid rise ($\approx5$ days) to a luminous peak ($M_V\approx-18.2$ mag) and plateau ($M_V\approx-17.6$ mag) extending to $30$ days with a fast decline rate of $\approx0.03$ mag day$^{-1}$. During the rising phase, $U-V$ color shows blueward evolution, followed by redward evolution in the plateau phase. Prominent flash features of hydrogen, helium, carbon, and nitrogen dominate the spectra up to $\approx5$ days after first light, with a transition to a higher ionization state in the first $\approx2$ days. Both the $U-V$ color and flash ionization states suggest a rise in the temperature, indicative of a delayed shock breakout inside dense circumstellar material (CSM). From the timescales of CSM interaction, we estimate its compact radial extent of $\sim(3-7)\times10^{14}$ cm. We then construct numerical light-curve models based on both continuous and eruptive mass-loss scenarios shortly before explosion. For the continuous mass-loss scenario, we infer a range of mass-loss history with $0.1-1.0\,M_\odot\,{\rm yr}^{-1}$ in the final $2-1$ yr before explosion, with a potentially decreasing mass loss of $0.01-0.1\,M_\odot\,{\rm yr}^{-1}$ in $\sim0.7-0.4$ yr toward the explosion. For the eruptive mass-loss scenario, we favor eruptions releasing $0.3-1\,M_\odot$ of the envelope at about a year before explosion, which result in CSM with mass and extent similar to the continuous scenario. We discuss the implications of the available multiwavelength constraints obtained thus far on the progenitor candidate and SN 2023ixf to our variable CSM models., Comment: Updated to match the published letter in ApJL, 2023 September 19

Eruptive mass loss of massive stars prior to supernova (SN) explosion is key to understanding their evolution and end fate. An observational signature of pre-SN mass loss is the detection of an early, short-lived peak prior to the radioactive-powered peak in the lightcurve of the SN. This is usually attributed to the SN shock passing through an extended envelope or circumstellar medium (CSM). Such an early peak is common for double-peaked Type IIb SNe with an extended Hydrogen envelope but is uncommon for normal Type Ibc SNe with very compact progenitors. In this paper, we systematically study a sample of 14 double-peaked Type Ibc SNe out of 475 Type Ibc SNe detected by the Zwicky Transient Facility. The rate of these events is ~ 3-9 % of Type Ibc SNe. A strong correlation is seen between the peak brightness of the first and the second peak. We perform a holistic analysis of this sample's photometric and spectroscopic properties. We find that six SNe have ejecta mass less than 1.5 Msun. Based on the nebular spectra and lightcurve properties, we estimate that the progenitor masses for these are less than ~ 12 Msun. The rest have an ejecta mass > 2.4 Msun and a higher progenitor mass. This sample suggests that the SNe with low progenitor masses undergo late-time binary mass transfer. Meanwhile, the SNe with higher progenitor masses are consistent with wave-driven mass loss or pulsation-pair instability-driven mass loss simulations., Comment: Accepted for publication in ApJ

SN 1987A was an unusual hydrogen-rich core-collapse supernova originating from a blue supergiant star. Similar blue supergiant explosions remain a small family of events, and are broadly characterized by their long rises to peak. The Zwicky Transient Facility (ZTF) Census of the Local Universe (CLU) experiment aims to construct a spectroscopically complete sample of transients occurring in galaxies from the CLU galaxy catalog. We identify 13 long-rising (>40 days) Type II supernovae from the volume-limited CLU experiment during a 3.5 year period from June 2018 to December 2021, approximately doubling the previously known number of these events. We present photometric and spectroscopic data of these 13 events, finding peak r-band absolute magnitudes ranging from -15.6 to -17.5 mag and the tentative detection of Ba II lines in 9 events. Using our CLU sample of events, we derive a long-rising Type II supernova rate of $1.37^{+0.26}_{-0.30}\times10^{-6}$ Mpc$^{-3}$ yr$^{-1}$, $\approx$1.4% of the total core-collapse supernova rate. This is the first volumetric rate of these events estimated from a large, systematic, volume-limited experiment., Comment: 33 pages, 18 figures, 5 tables. Updated to ApJ accepted version

We present optical photometry and spectroscopy of the Type IIn supernova (SN) 2021qqp. Its unusual light curve is marked by a long precursor for $\approx300$ days, a rapid increase in brightness for $\approx60$ days, and then a sharp increase of $\approx1.6$ mag in only a few days to a first peak of $M_r \approx -19.5$ mag. The light curve then declines rapidly until it re-brightens to a second distinct peak of $M_r \approx -17.3$ mag centered at $\approx335$ days after the first peak. The spectra are dominated by Balmer lines with a complex morphology, including a narrow component with a width of $\approx 1300$ km s$^{-1}$ (first peak) and $\approx 2500$ km s$^{-1}$ (second peak) that we associate with the circumstellar medium (CSM) and a P Cygni component with an absorption velocity of $\approx 8500$ km s$^{-1}$ (first peak) and $\approx 5600$ km s$^{-1}$ (second peak) that we associate with the SN-CSM interaction shell. Using the luminosity and velocity evolution, we construct a flexible analytical model, finding two significant mass-loss episodes with peak mass loss rates of $\approx 10$ and $\approx 5\,M_{\odot}$ yr$^{-1}$ about $0.8$ and $2$ yr before explosion, respectively, with a total CSM mass of $\approx 2-4\,M_{\odot}$. We show that the most recent mass-loss episode could explain the precursor for the year preceding the explosion. The SN ejecta mass is constrained to be $\approx 5-30\,M_{\odot}$ for an explosion energy of $\approx (3-10)\times10^{51}$ erg. We discuss eruptive massive stars (luminous blue variable, pulsational pair instability) and an extreme stellar merger with a compact object as possible progenitor channels., Comment: Updated to match the published article in ApJ, 2024 March 29

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

Early observations of transient explosions can provide vital clues to their progenitor origins. In this paper we present the nearby Type Iax (02cx-like) supernova (SN), SN 2020udy that was discovered within hours ($\sim$7 hr) of estimated first light. An extensive dataset of ultra-violet, optical, and near-infrared observations was obtained, covering out to $\sim$150 d after explosion. SN 2020udy peaked at -17.86$\pm$0.43 mag in the r band and evolved similarly to other 'luminous' SNe Iax, such as SNe 2005hk and 2012Z. Its well-sampled early light curve allows strict limits on companion interaction to be placed. Main-sequence companion stars with masses of 2 and 6 M$_\odot$ are ruled out at all viewing angles, while a helium-star companion is allowed from a narrow range of angles (140-180$^\circ$ away from the companion). The spectra and light curves of SN2020udy are in good agreement with those of the 'N5def' deflagration model of a near Chandrasekhar-mass carbon-oxygen white dwarf. However, as has been seen in previous studies of similar luminosity events, SN 2020udy evolves slower than the model. Broad-band linear polarisation measurements taken at and after peak are consistent with no polarisation, in agreement with the predictions of the companion-star configuration from the early light curve measurements. The host galaxy environment is low metallicity and is consistent with a young stellar population. Overall, we find the most plausible explosion scenario to be the incomplete disruption of a CO white dwarf near the Chandrasekhar-mass limit, with a helium-star companion., Comment: 18 pages, 14 figures, submitted to MNRAS

The Galactic nova rate is intimately linked to our understanding of its chemical enrichment and progenitor channels of Type Ia supernovae. Yet past estimates have varied by more than an order of magnitude ($\approx10-300$ yr$^{-1}$) owing to limitations in both discovery methods as well as assumptions regarding the Galactic dust distribution and extragalactic stellar populations. Recent estimates utilizing synoptic near-infrared surveys have begun to provide a glimpse of a consensus ($\approx25-50$ yr$^{-1}$); however, a consistent estimate remains lacking. Here, we present the first all-sky search for Galactic novae using 8 years of data from the NEOWISE mid-infrared (MIR) survey. Operating at $3.4$ and $4.6$ $\mu$m where interstellar extinction is negligible, the 6-month cadence NEOWISE dataset offers unique sensitivity to discover slowly evolving novae across the entire Galaxy. Using a novel image subtraction pipeline together with systematic selection criteria, we identify a sample of 49 rapidly evolving MIR outbursts as candidate Galactic novae. While 27 of these sources are known novae, the remaining are previously missed nova candidates discovered in this work. The unknown events are spatially clustered along the densest and most heavily obscured regions of the Galaxy where previous novae are severely underrepresented. We use simulations of the NEOWISE survey strategy, the pipeline detection efficiency, and our criteria to derive a Galactic nova rate of $47.9^{+3.1}_{-8.3}$ yr$^{-1}$. The discovery of these exceptionally bright (yet overlooked) nova candidates confirm emerging suggestions that optical surveys have been highly incomplete in searches for Galactic novae, highlighting the potential for MIR searches in revealing the demographics of Galactic stellar outbursts., Comment: 13 pages, 4 figures, 1 table, 3 figures in Appendix. Submitted to MNRAS. Comments welcome

While the vast majority of Tidal Disruption Events (TDEs) has been identified by wide-field sky surveys in the optical and X-ray bands, recent studies indicate that a considerable fraction of TDEs may be dust obscured, and thus preferentially detected in the infrared (IR) wavebands. In this Letter, we present the discovery of a luminous mid-IR nuclear flare (termed WTP 14adbjsh) identified in a systematic transient search of archival images from the NEOWISE mid-IR survey. The source reached a peak luminosity of $L \simeq 10^{43} \text{erg s}^{-1}$ at 4.6 ${\mu}$m in 2015, before fading in the IR with a TDE-like $F \propto t^{-5/3}$ decline, radiating a total of more than $ 3\times 10^{51}$ erg in the last 7 years. The transient event took place in the nearby galaxy NGC 7392, at a distance of around 42 Mpc; yet, no optical or X-ray flare is detected. We interpret the transient as the nearest TDE candidate detected in the last decade, which was missed at other wavelengths due to dust obscuration, hinting at the existence of TDEs that have been historically overlooked. Unlike most previously detected TDEs, the transient was discovered in a star forming galaxy, corroborating earlier suggestions that dust obscuration suppresses significantly the detection of TDEs in these environments. Our results demonstrate that the study of IR-detected TDEs is critical in order to obtain a complete understanding of the physics of TDEs, and to conclude whether TDEs occur preferentially in a particular class of galaxies., Comment: 13 pages, 4 figures, accepted for publication by ApJL

One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing $r$-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of $r$-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of $r$-process nucleosynthesis in the binary neutron star merger GW170817 was its long-lasting near-infrared emission, thus motivating a systematic photometric study of the light curves of broadlined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL -- including 18 observed with the Zwicky Transient Facility and 7 from the literature -- in the optical/near-infrared bands to determine what quantity of $r$-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for $r$-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the $r$-process mass for these SNe. We also perform independent light curve fits to models without $r$-process. We find that the $r$-process-free models are a better fit to the light curves of the objects in our sample. Thus we find no compelling evidence of $r$-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of $r$-process ejecta mass or indicate whether all collapsars are completely devoid of $r$-process nucleosynthesis., Comment: 37 pages, 8 figures, 4 tables

This paper is one in a series reporting results from small telescope observations of variable young stars. Here, we study the repeating outbursts of three likely Be stars based on long-term optical, near-infrared, and mid-infrared photometry for all three objects, along with follow-up spectra for two of the three. The sources are characterised as rare, truly regularly outbursting Be stars. We interpret the photometric data within a framework for modelling light curve morphology, and find that the models correctly predict the burst shapes, including their larger amplitudes and later peaks towards longer wavelengths. We are thus able to infer the start and end times of mass loading into the circumstellar disks of these stars. The disk sizes are typically 3-6 times the areas of the central star. The disk temperatures are ~40%, and the disk luminosities are ~10% of those of the central Be star, respectively. The available spectroscopy is consistent with inside-out evolution of the disk. Higher excitation lines have larger velocity widths in their double-horned shaped emission profiles. Our observations and analysis support the decretion disk model for outbursting Be stars., Comment: Accepted for publication by MNRAS. 20 pages, 11 figures, 2 tables

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

Among the supernovae (SNe) that show strong interaction with the circumstellar medium, there is a rare subclass of Type Ia supernovae, SNe Ia-CSM, that show strong narrow hydrogen emission lines much like SNe IIn but on top of a diluted over-luminous Type Ia spectrum. In the only previous systematic study of this class (Silverman et al. 2013), 16 objects were identified, 8 historic and 8 from the Palomar Transient Factory (PTF). Now using the successor survey to PTF, the Zwicky Transient Facility (ZTF), we have classified 12 additional objects of this type through the systematic Bright Transient Survey (BTS). In this study, we present and analyze the optical and mid-IR light curves, optical spectra, and host galaxy properties of this sample. Consistent with previous studies, we find the objects to have slowly evolving light curves compared to normal SNe Ia with peak absolute magnitudes between -19.1 and -21, spectra having weak H$\beta$, large Balmer decrements of ~7 and strong Ca NIR emission. Out of 10 SNe from our sample observed by NEOWISE, 9 have $3\sigma$ detections, along with some showing a clear reduction in red-wing of H$\alpha$, indicative of newly formed dust. We do not find our SN Ia-CSM sample to have a significantly different distribution of equivalent width of He I $\lambda5876$ than SNe IIn as observed in Silverman et al. 2013. The hosts tend to be late-type galaxies with recent star formation. We also derive a rate estimate of 29$^{+27}_{-21}$ Gpc$^{-3}$ yr$^{-1}$ for SNe Ia-CSM which is ~0.02--0.2 % of the SN Ia rate. This work nearly doubles the sample of well-studied Ia-CSM objects in Silverman et al. 2013, increasing the total number to 28., Comment: 26 pages, 14 figures, 5 tables. Submitted to ApJ

Narrow transient emission lines (flash-ionization features) in early supernova (SN) spectra trace the presence of circumstellar material (CSM) around the massive progenitor stars of core-collapse SNe. The lines disappear within days after the SN explosion, suggesting that this material is spatially confined, and originates from enhanced mass loss shortly (months to a few years) prior to explosion. We performed a systematic survey of H-rich (Type II) SNe discovered within less than two days from explosion during the first phase of the Zwicky Transient Facility (ZTF) survey (2018-2020), finding thirty events for which a first spectrum was obtained within $< 2$ days from explosion. The measured fraction of events showing flash ionisation features ($>36\%$ at $95\%$ confidence level) confirms that elevated mass loss in massive stars prior to SN explosion is common. We find that SNe II showing flash ionisation features are not significantly brighter, nor bluer, nor more slowly rising than those without. This implies that CSM interaction does not contribute significantly to their early continuum emission, and that the CSM is likely optically thin. We measured the persistence duration of flash ionisation emission and find that most SNe show flash features for $\approx 5 $ days. Rarer events, with persistence timescales $>10$ days, are brighter and rise longer, suggesting these may be intermediate between regular SNe II and strongly-interacting SNe IIn.

We present observations of a peculiar hydrogen- and helium-poor stripped-envelope (SE) supernova (SN) 2020wnt, primarily in the optical and near-infrared (near-IR). Its peak absolute bolometric magnitude of -20.9 mag and a rise time of 69~days are reminiscent of hydrogen-poor superluminous SNe (SLSNe~I), luminous transients potentially powered by spinning-down magnetars. Before the main peak, there is a brief peak lasting <10 days post-explosion, likely caused by interaction with circumstellar medium (CSM) ejected ~years before the SN explosion. The optical spectra near peak lack a hot continuum and OII absorptions, which are signs of heating from a central engine; they quantitatively resemble those of radioactivity-powered H/He-poor Type Ic SESNe. At ~1 year after peak, nebular spectra reveal a blue pseudo-continuum and narrow OI recombination lines associated with magnetar heating. Radio observations rule out strong CSM interactions as the dominant energy source at +266 days post peak. Near-IR observations at +200-300 day reveal carbon monoxide and dust formation, which causes a dramatic optical light curve dip. Pair-instability explosion models predict slow light curve and spectral features incompatible with observations. SN 2020wnt is best explained as a magnetar-powered core-collapse explosion of a 28 Msun pre-SN star. The explosion kinetic energy is significantly larger than the magnetar energy at peak, effectively concealing the magnetar-heated inner ejecta until well after peak. SN 2020wnt falls into a continuum between normal SNe Ic and SLSNe I and demonstrates that optical spectra at peak alone cannot rule out the presence of a central engine., Comment: Submitted to ApJ

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

The dividing line between gamma-ray bursts (GRBs) and ordinary stripped-envelope core-collapse supernovae (SNe) is yet to be fully understood. Observationally mapping the variety of ejecta outcomes (ultra-relativistic, mildly-relativistic or non-relativistic) in SNe of Type Ic with broad lines (Ic-BL) can provide a key test to stellar explosion models. However, this requires large samples of the rare Ic-BL events with follow-up observations in the radio, where fast ejecta can be probed largely free of geometry and viewing angle effects. Here, we present the results of a radio (and X-ray) follow-up campaign of 16 SNe Ic-BL detected by the Zwicky Transient Facility (ZTF). Our radio campaign resulted in 4 counterpart detections and 12 deep upper limits. None of the events in our sample is as relativistic as SN 1998bw and we constrain the fraction of SN 1998bw-like explosions to $< 19\%$ (3$\sigma$ Gaussian equivalent), a factor of $\approx 2$ smaller than previously established. We exclude relativistic ejecta with radio luminosity densities in between $\approx 5\times10^{27}$ erg s$^{-1}$ Hz$^{-1}$ and $\approx 10^{29}$ erg s$^{-1}$ Hz$^{-1}$ at $t\gtrsim 20$ d since explosion for $\approx 60\%$ of the events in our sample. This shows that SNe Ic-BL similar to the GRB-associated SN 1998bw, SN 2003lw, SN 2010dh, or to the relativistic SN 2009bb and iPTF17cw, are rare. Our results also exclude an association of the SNe Ic-BL in our sample with largely off-axis GRBs with energies $E\gtrsim 10^{50}$ erg. The parameter space of SN2006aj-like events (faint and fast-peaking radio emission) is, on the other hand, left largely unconstrained and systematically exploring it represents a promising line of future research., Comment: 28 pages, 11 figures, submitted to AAS journal

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

The fate of stars in the zero-age main-sequence (ZAMS) range $\approx 8-12$ Msun is unclear. They could evolve to form white dwarfs or explode as electron-capture supernovae (SNe) or iron core-collapse SNe (CCSNe). Even though the initial mass function indicates that this mass range should account for over 40% of all CCSNe progenitors, few have been observationally confirmed, likely owing to the faintness and rapid evolution of these transients. In this paper, we present a sample of nine Ca-rich/O-poor Type IIb SNe detected by the Zwicky Transient Facility with progenitors likely in this mass range. We perform a holistic analysis of the spectroscopic and photometric properties of the sample. These sources have a flux ratio of [Ca II] $\lambda \lambda$7291, 7324 to [O I] $\lambda \lambda$6300, 6364 of $\gtrsim$ 2 in their nebular spectra. Comparing the measured [O I] luminosity ($\lesssim 10^{39} \mathrm{erg\ s^{-1}}$) and derived oxygen mass ($\lesssim 0.1$ Msun) with theoretical models, we infer that the progenitor ZAMS mass for these explosions is less than 12 Msun. These correspond to He-stars with core masses less than around 3 Msun. We find that the ejecta properties (Mej $\lesssim 1$ Msun) are also consistent with those expected for such low mass He-stars. The low ejecta mass of these sources indicates a class of strongly-stripped SNe that is a transition between the regular stripped-envelope SNe and ultra-stripped SNe. The progenitor could be stripped by a main sequence companion and result in the formation of a neutron star $-$ main sequence binary. Such binaries have been suggested to be progenitors of neutron star $-$ white dwarf systems that could merge within a Hubble time, and be detectable with LISA., Comment: Accepted for publication in ApJ

The detonation of a thin ($\lesssim$$0.03\,\mathrm{M_\odot}$) helium shell (He-shell) atop a $\sim$$1\,\mathrm{M_\odot}$ white dwarf (WD) is a promising mechanism to explain normal Type Ia supernovae (SNe Ia), while thicker He-shells and less massive WDs may explain some recently observed peculiar SNe Ia. We present observations of SN 2020jgb, a peculiar SN Ia discovered by the Zwicky Transient Facility (ZTF). Near maximum light, SN 2020jgb is slightly subluminous (ZTF $g$-band absolute magnitude $M_g$ between $-18.2$ and $-18.7$ mag depending on the amount of host galaxy extinction) and shows an unusually red color ($g_\mathrm{ZTF}-r_\mathrm{ZTF}$ between 0.4 and 0.2 mag) due to strong line-blanketing blueward of $\sim$5000 $\r{A}$. These properties resemble those of SN 2018byg, a peculiar SN Ia consistent with a thick He-shell double detonation (DDet) SN. Using detailed radiative transfer models, we show that the optical spectroscopic and photometric evolution of SN 2020jgb are broadly consistent with a $\sim$$0.95\,\mathrm{M_\odot}$ (C/O core + He-shell; up to $\sim$$1.00\,\mathrm{M_\odot}$ depending on the total host extinction) progenitor ignited by a thick ($\sim$$0.13\,\mathrm{M_\odot}$) He-shell. We detect a prominent absorption feature at $\sim$1 $\mu\mathrm{m}$ in the near-infrared (NIR) spectrum of SN 2020jgb, which could originate from unburnt helium in the outermost ejecta. While the sample size is limited, similar 1 $\mu\mathrm{m}$ features have been detected in all the thick He-shell DDet candidates with NIR spectra obtained to date. SN 2020jgb is also the first subluminous, thick He-shell DDet SN discovered in a star-forming galaxy, indisputably showing that He-shell DDet objects occur in both star-forming and passive galaxies, consistent with the normal SN Ia population., Comment: 23 pages, 10 figures. Updated to accepted version (ApJ)

S190426c / GW190426_152155 was the first probable neutron star - black hole merger candidate detected by the LIGO-Virgo Collaboration. We undertook a tiled search for optical counterparts of this event using the 0.7m GROWTH-India Telescope. Over a period of two weeks, we obtained multiple observations over a 22.1 deg^2 area, with a 17.5% probability of containing the source location. Initial efforts included obtaining photometry of sources reported by various groups, and a visual search for sources in all galaxies contained in the region. Subsequently, we have developed an image subtraction and candidate vetting pipeline with ~ 94% efficiency for transient detection. Processing the data with this pipeline, we find several transients, but none that are compatible with kilonova models. We present the details of our observations, working of our pipeline, results from the search, and our interpretations of the non-detections that will work as a pathfinder during the O4 run of LVK., Comment: 12 pages, 6 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

We present observations of an extreme radio flare, VT J024345.70-284040.08, hereafter VT J0243, from the nucleus of a galaxy with evidence for historic Seyfert activity at redshift $z=0.074$. Between NRAO VLA Sky Survey observations in 1993 to VLA Sky Survey observations in 2018, VT J0243 rose from a ${\sim}$GHz radio luminosity of $\nu L_\nu \lesssim 10^{38}$ erg s$^{-1}$ to $\nu L_\nu{\sim}10^{40}$ erg s$^{-1}$, and still continues to brighten. The radio spectral energy distribution (SED) evolution is consistent with a nascent jet that has slowed over ${\sim}3000$ days with an average $0.1 < \langle \beta \rangle < 0.6$. The jet is energetic (${\sim}10^{51-52}$ erg), and had a radius ${\sim}0.7$ pc in Dec. 2021. X-ray observations suggest a persistent or evolving corona, possibly associated with an accretion disk, and IR and optical observations constrain any high-energy counterpart to be sub-Eddington. VT J0243 may be an example of a young, off-axis radio jet from a slowly evolving tidal disruption event. Other more mysterious triggers for the accretion enhancement and jet launching are possible. In either case, VT J0243 is a unique example of a nascent jet, highlighting the unknown connection between supermassive black holes, the properties of their accretion flows, and jet launching., Comment: 20 pages, 5 figures, 3 tables. Submitted to ApJ

We present the design and performance of the GROWTH-India telescope, a 0.7 m robotic telescope dedicated to time-domain astronomy. The telescope is equipped with a 4k back-illuminated camera giving a 0.82-degree field of view and sensitivity of m_g ~20.5 in 5-min exposures. Custom software handles observatory operations: attaining high on-sky observing efficiencies (>~ 80%) and allowing rapid response to targets of opportunity. The data processing pipelines are capable of performing PSF photometry as well as image subtraction for transient searches. We also present an overview of the GROWTH-India telescope's contributions to the studies of Gamma-ray Bursts, the electromagnetic counterparts to gravitational wave sources, supernovae, novae and solar system objects., Comment: 17 pages, 8 figures, Accepted for publication in The Astronomical Journal

We present the pulse arrival times and high-precision dispersion measure estimates for 14 millisecond pulsars observed simultaneously in the 300-500 MHz and 1260-1460 MHz frequency bands using the upgraded Giant Metrewave Radio Telescope (uGMRT). The data spans over a baseline of 3.5 years (2018-2021), and is the first official data release made available by the Indian Pulsar Timing Array collaboration. This data release presents a unique opportunity for investigating the interstellar medium effects at low radio frequencies and their impact on the timing precision of pulsar timing array experiments. In addition to the dispersion measure time series and pulse arrival times obtained using both narrowband and wideband timing techniques, we also present the dispersion measure structure function analysis for selected pulsars. Our ongoing investigations regarding the frequency dependence of dispersion measures have been discussed. Based on the preliminary analysis for five millisecond pulsars, we do not find any conclusive evidence of chromaticity in dispersion measures. Data from regular simultaneous two-frequency observations are presented for the first time in this work. This distinctive feature leads us to the highest precision dispersion measure estimates obtained so far for a subset of our sample. Simultaneous multi-band uGMRT observations in Band 3 and Band 5 are crucial for high-precision dispersion measure estimation and for the prospect of expanding the overall frequency coverage upon the combination of data from the various Pulsar Timing Array consortia in the near future. Parts of the data presented in this work are expected to be incorporated into the upcoming third data release of the International Pulsar Timing Array., Comment: 23 pages, 21 figures, 3 tables. Published in PASA

We present photometric and spectroscopic observations of the Type IIn supernova SN 2019zrk (also known as ZTF20aacbyec). The SN shows a $\gtrsim$ 100 day precursor, with a slow rise, followed by a rapid rise to M $\sim -19.2$ in the $r$ and $g$ bands. The post-peak light-curve decline is well fit with an exponential decay with a timescale of $\sim 39$ days, but it shows prominent undulations, with an amplitude of $\sim 1$ mag. Both the light curve and spectra are dominated by an interaction with a dense circumstellar medium (CSM), probably from previous mass ejections. The spectra evolve from a scattering-dominated Type IIn spectrum to a spectrum with strong P-Cygni absorptions. The expansion velocity is high, $\sim 16,000$ km s$^{-1}$, even in the last spectra. The last spectrum $\sim 110$ days after the main eruption reveals no evidence for advanced nucleosynthesis. From analysis of the spectra and light curves, we estimate the mass-loss rate to be $\sim 4 \times 10^{-2}$ M$_\odot$ yr$^{-1}$ for a CSM velocity of 100 km s$^{-1}$, and a CSM mass of $\gtrsim 1$ M$_\odot$. We find strong similarities for both the precursor, general light curve, and spectral evolution with SN 2009ip and similar SNe, although SN 2019zrk displays a brighter peak magnitude. Different scenarios for the nature of the 09ip-class of SNe, based on pulsational pair instability eruptions, wave heating, and mergers, are discussed. }, Comment: 18 pages, 12 figures. Astronomy & Astrophysics, in press

We present the discovery and multi-wavelength characterization of SRGA J181414.6-225604, a Galactic hard X-ray transient discovered during the ongoing SRG/ART-XC sky survey. Using data from the Palomar Gattini-IR survey, we identify a spatially and temporally coincident variable infrared (IR) source, IRAS 18111-2257, and classify it as a very late-type (M7-M8), long period ($1502 \pm 24$ days) and luminous ($M_K\approx -9.9 \pm 0.2$) O-rich Mira donor star located at a distance of $\approx 14.6^{+2.9}_{-2.3}$ kpc. Combining multi-color photometric data over the last $\approx 25$ years, we show that the IR counterpart underwent a recent (starting $\approx 800$ days before the X-ray flare) enhanced mass loss (reaching $\approx 2.1 \times 10^{-5}$ M$_\odot$ yr$^{-1}$) episode resulting in an expanding dust shell obscuring the underlying star. Multi-epoch follow-up from Swift, NICER and NuSTAR reveal a $\approx 200$ day long X-ray outburst reaching a peak luminosity of $L_X \approx 2.5 \times 10^{36}$ erg s$^{-1}$, characterized by a heavily absorbed ($N_{\rm H} \approx 6\times 10^{22}$ cm$^{-2}$) X-ray spectrum consistent with an optically thick Comptonized plasma. The X-ray spectral and timing behavior suggest the presence of clumpy wind accretion together with a dense ionized nebula overabundant in silicate material surrounding the compact object. Together, we show that SRGA J181414.6-225604 is a new symbiotic X-ray binary in outburst, triggered by an intense dust formation episode of a highly evolved donor. Our results offer the first direct confirmation for the speculated connection between enhanced late-stage donor mass loss and active lifetimes of the symbiotic X-ray binaries., Comment: 32 pages, 15 figures, Accepted for publication in ApJ

We discuss the central role that dust condensation plays in shaping the observational appearance of outflows from coalescing binary systems. As binaries begin to coalesce, they shock-heat and expel material into their surroundings. Depending on the properties of the merging system, this material can expand to the point where molecules and dust form, dramatically increasing the gas opacity. We use the existing population of Luminous Red Novae (LRNe) to constrain the thermodynamics of these ejecta, then apply our findings to the progressive obscuration of merging systems in the lead in to their coalescence. Compact progenitor stars near the main sequence or in the Hertzsprung gap along with massive progenitor stars have sufficiently hot circumstellar material to remain unobscured by dust. By contrast, more extended, low-mass giants should become completely optically obscured by dust formation in the circumbinary environment. We predict that 30--50\% of stellar coalescence transients for solar-mass stars will be dusty, infrared-luminous sources. Of these, the optical transients may selectively trace complete merger outcomes while the infrared transients trace common envelope ejection outcomes., Comment: v2, published in AAS Journals

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution., Comment: 63 pages, 15 figures, 2 tables, published in Nature on May 4, 2022