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