Hunting for Hidden Explosions: Exploring the Transient Infrared Sky with the Spitzer Space Telescope

The study of time-variable astronomical phenomena is undergoing an explosive renaissance ushered in by recent advancements in capabilities to monitor the sky from radio to gamma rays. The infrared region of the electromagnetic spectrum provides a unique window to uncover a vast array of stellar eruptions and explosions that are otherwise obscured; however, the dynamic infrared sky has remained largely unexplored. To uncover these hidden cosmic explosions, I undertook a systematic search in the infrared with the Spitzer Space Telescope called SPIRITS — the SPitzer InfraRed Intensive Transients Survey. Targeting a specially chosen sample of nearby galaxies, the search revealed a bounty of exceptionally red events, and may have discovered entirely new and diverse populations, including deeply embedded supernovae, catastrophic stellar mergers, and giant eruptions of massive, violently unstable stars. Providing a first census of stellar infrared transients, SPIRITS has paved the way for new and upcoming surveys to further expand our exploration of the dynamic infrared sky.</p

An Optical and Infrared Time-Domain Study of the Supergiant Fast X-ray Transient Candidate IC 10 X-2

We present an optical and infrared (IR) study of IC 10 X-2, a high-mass X-ray binary in the galaxy IC 10. Previous optical and X-ray studies suggest X-2 is a Supergiant Fast X-ray Transient: a large-amplitude (factor of $\sim$ 100), short-duration (hours to weeks) X-ray outburst on 2010 May 21. We analyze R- and g-band light curves of X-2 from the intermediate Palomar Transient Factory taken between 2013 July 15 and 2017 Feb 14 show high-amplitude ($\gtrsim$ 1 mag), short-duration ($\lesssim8$ d) flares and dips ($\gtrsim$ 0.5 mag). Near-IR spectroscopy of X-2 from Palomar/TripleSpec show He I, Paschen-$\gamma$, and Paschen-$\beta$ emission lines with similar shapes and amplitudes as those of luminous blue variables (LBVs) and LBV candidates (LBVc). Mid-IR colors and magnitudes from Spitzer/IRAC photometry of X-2 resemble those of known LBV/LBVcs. We suggest that the stellar companion in X-2 is an LBV/LBVc and discuss possible origins of the optical flares. Dips in the optical light curve are indicative of eclipses from optically thick clumps formed in the winds of the stellar counterpart. Given the constraints on the flare duration ($0.02 - 0.8$ d) and the time between flares ($15.1\pm7.8$ d), we estimate the clump volume filling factor in the stellar winds, $f_V$, to be $0.01 < f_V < 0.71$, which overlaps with values measured from massive star winds. In X-2, we interpret the origin of the optical flares as the accretion of clumps formed in the winds of an LBV/LBVc onto the compact object.Comment: 15 pages, 4 figures. Submitted to ApJ on Sep 26 201

The 21-SPONGE HI Absorption Survey I: Techniques and Initial Results

We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $\sigma_{\tau}=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases.Comment: Accepted for publication in ApJ; 24 pages, 14 figure

Mapping the Inner Structure of Quasars with Time-Domain Spectroscopy

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 deg^2 field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hβ to z=8

SPIRITS 16tn in NGC 3556: A heavily obscured and low-luminosity supernova at 8.8 Mpc

We present the discovery by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) of a likely supernova (SN) in NGC 3556 at only 8.8 Mpc, which, despite its proximity, was not detected by optical searches. A luminous infrared (IR) transient at $M_{[4.5]} = -16.7$ mag (Vega), SPIRITS 16tn is coincident with a dust lane in the inclined, star-forming disk of the host. Using IR, optical, and radio observations, we attempt to determine the nature of this event. We estimate $A_V \approx$ 8 - 9 mag of extinction, placing it among the three most highly obscured IR-discovered SNe to date. The [4.5] light curve declined at a rate of 0.013 mag day$^{-1}$, and the $[3.6] - [4.5]$ color grew redder from 0.7 to $\gtrsim$ 1.0 mag by 184.7 days post discovery. Optical/IR spectroscopy shows a red continuum, but no clearly discernible features, preventing a definitive spectroscopic classification. Deep radio observations constrain the radio luminosity of SPIRITS 16tn to $L_{\nu} \lesssim 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ between 3 - 15 GHz, excluding many varieties of radio core-collapse SNe. A type Ia SN is ruled out by the observed red IR color, and lack of features normally attributed to Fe-peak elements in the optical and IR spectra. SPIRITS 16tn was fainter at [4.5] than typical stripped-envelope SNe by $\approx$ 1 mag. Comparison of the spectral energy distribution to SNe II suggests SPIRITS 16tn was both highly obscured, and intrinsically dim, possibly akin to the low-luminosity SN 2005cs. We infer the presence of an IR dust echo powered by a peak luminosity of the transient of $5 \times 10^{40}$ erg s$^{-1} < L_{\mathrm{peak}} < 4\times10^{43}$ erg s$^{-1}$, consistent with the observed range for SNe II. This discovery illustrates the power of IR surveys to overcome the compounding effects of visible extinction and optically sub-luminous events in completing the inventory of nearby SNe.Comment: 25 pages, 10 figures, submitted to Ap

Recurring outbursts of the supernova impostor AT 2016blu in NGC 4559

We present the first photometric analysis of the supernova (SN) impostor AT 2016blu in NGC 4559. This transient was discovered by the Lick Observatory Supernova Search in 2012 and has continued its outbursts since then. Optical and infrared photometry of AT 2016blu reveals at least 19 outbursts in 2012-2022. Similar photometry from 1999-2009 shows no outbursts, indicating that the star was relatively stable in the decade before discovery. Archival {\it Hubble Space Telescope} observations suggest that the progenitor had a minimum initial mass of $M >= 33$ M$_{\odot}$ and a luminosity of $L >= 10^{5.7}$ L$_{\odot}$. AT 2016blu's outbursts show irregular variability with multiple closely spaced peaks having typical amplitudes of 1-2 mag and durations of 1-4 weeks. While individual outbursts have irregular light curves, concentrations of these peaks recur with a period of $\sim 113 \pm 2$ d. Based on this period, we predict times for upcoming outbursts in 2023 and 2024. AT 2016blu shares similarities with SN 2000ch in NGC 3432, where outbursts may arise from periastron encounters in an eccentric binary containing a luminous blue variable (LBV). We propose that AT 2016blu's outbursts are also driven by interactions that intensify around periastron in an eccentric system. Intrinsic variability of the LBV-like primary star may cause different intensity and duration of binary interaction at each periastron passage. AT 2016blu also resembles the periastron encounters of $\eta$ Carinae prior to its Great Eruption and the erratic pre-SN eruptions of SN 2009ip. This similarity and the onset of eruptions in the past decade hint that AT 2016blu may also be headed for a catastrophe, making it a target of great interest.Comment: 18 pages, 14 figures, 6 tables, MNRAS Accepte