M31-2015-LRN is a likely stellar merger discovered in the Andromeda Galaxy in 2015. We present new optical to mid-infrared photometry and optical spectroscopy for this event. Archival data shows that the source started to brighten $\sim$2 years before the nova event. During this precursor phase, the source brightened by $\sim$3 mag. The lightcurve at 6 and 1.5 months before the main outburst may show periodicity, with periods of 16$\pm$0.3 and 28.1$\pm$1.4 days respectively. This complex emission may be explained by runaway mass loss from the system after the binary undergoes Roche-lobe overflow, leading the system to coalesce in tens of orbital periods. While the progenitor spectral energy distribution shows no evidence of pre-existing warm dust in system, the remnant forms an optically thick dust shell at $\sim$4 months after the outburst peak. The optical depth of the shell increases dramatically after 1.5 years, suggesting the existence of shocks that enhance the dust formation process. We propose that the merger remnant is likely an inflated giant obscured by a cooling shell of gas with mass $\sim0.2$ M$_{\odot}$ ejected at the onset of the common envelope phase.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

We present a systematic study of the most luminous ($M_{\mathrm{IR}}$ [Vega magnitudes] brighter than $-14$) infrared (IR) transients discovered by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) between 2014 and 2018 in nearby galaxies ($D < 35$ Mpc). The sample consists of nine events that span peak IR luminosities of $M_{[4.5],\mathrm{peak}}$ between $-14$ and $-18.2$, show IR colors between $0.2 < ([3.6]{-}[4.5]) < 3.0$, and fade on timescales between $55$ days $< t_{\mathrm{fade}} < 480$ days. The two reddest events ($A_V > 12$) show multiple, luminous IR outbursts over several years and have directly detected, massive progenitors in archival imaging. With analyses of extensive, multiwavelength follow-up, we suggest the following possible classifications: five obscured core-collapse supernovae (CCSNe), two erupting massive stars, one luminous red nova, and one intermediate-luminosity red transient. We define a control sample of all optically discovered transients recovered in SPIRITS galaxies and satisfying the same selection criteria. The control sample consists of eight CCSNe and one Type Iax SN. We find that 7 of the 13 CCSNe in the SPIRITS sample have lower bounds on their extinction of $2 < A_V < 8$. We estimate a nominal fraction of CCSNe in nearby galaxies that are missed by optical surveys as high as $38.5^{+26.0}_{-21.9}$% (90% confidence). This study suggests that a significant fraction of CCSNe may be heavily obscured by dust and therefore undercounted in the census of nearby CCSNe from optical searches.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

We present a catalog of 417 luminous infrared variable stars with periods exceeding 250 days. These were identified in 20 nearby galaxies by the ongoing SPIRITS survey with the Spitzer Space Telescope. Of these, 359 variables have $M_{[4.5]}$ (phase-weighted mean magnitudes) fainter than $-12$ and periods and luminosities consistent with previously reported variables in the Large Magellanic Cloud. However, 58 variables are more luminous than $M_{[4.5]} = -12$, including 11 that are brighter than $M_{[4.5]} = -13$ with the brightest having $M_{[4.5]} = -15.51$. Most of these bright variable sources have quasi-periods longer than 1000 days, including four over 2000 days. We suggest that the fundamental period-luminosity relationship, previously measured for the Large Magellanic Cloud, extends to much higher luminosities and longer periods in this large galaxy sample. We posit that these variables include massive AGB stars (possibly super-AGB stars), red supergiants experiencing exceptionally high mass-loss rates, and interacting binaries. We also present 3.6, 4.5, 5.8 and 8.0 $\mu$m photometric catalogs for all sources in these 20 galaxies.