We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the 5 deg 2 field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of ≈ 21.3 AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of arXiv:1906.06345v2 [astro-ph.HE]

We present photometric and spectroscopic observations of the nearby Type Ia SN 2019yvq, from its discovery ∼1 day after explosion to ∼100 days after its peak brightness. This SN exhibits several unusual features, most notably an extremely bright UV excess seen within ∼5 days of its explosion. As seen in Swift UV data, this early excess outshines its “peak” brightness, making this object more extreme than other supernovae (SNe) with early UV/blue excesses (e.g., iPTF14atg and SN 2017cbv). In addition, it was underluminous M B = −18.4, relatively quickly declining (Δm 15(B) = 1.37), and shows red colors past its early blue bump. Unusual (although not unprecedented) spectral features include extremely broad-lined and high-velocity Si absorption. Despite obvious differences in peak spectra, we classify SN 2019yvq as a transitional member of the 02es-like subclass due to its similarities in several respects (e.g., color, peak luminosity, peak Ti, and nebular [Ca ii]). We model this data set with a variety of published models, including SN ejecta–companion shock interaction and sub-Chandrasekhar-mass white dwarf (WD) double-detonation models. Radio constraints from the VLA place an upper limit of (4.5–20) × 10−8 M ⊙ yr−1 on the mass-loss rate from a symbiotic progenitor, which does not exclude a red giant or main-sequence companion. Ultimately, we find that no one model can accurately replicate all aspects of the data set, and further we find that the ubiquity of early excesses in 02es-like SNe Ia requires a progenitor system that is capable of producing isotropic UV flux, ruling out some models for this class of objects.

We present the results of our monitoring campaigns of the luminous red novae (LRNe) AT 2020hat in NGC 5068 and AT 2020kog in NGC 6106. The two objects were imaged (and detected) before their discovery by routine survey operations. They show a general trend of slow luminosity rise, lasting at least a few months. The subsequent major LRN outbursts were extensively followed in photometry and spectroscopy. The light curves present an initial short-duration peak, followed by a redder plateau phase. AT 2020kog is a moderately luminous event peaking at ∼7 × 1040 erg s−1, while AT 2020hat is almost one order of magnitude fainter than AT 2020kog, although it is still more luminous than V838 Mon. In analogy with other LRNe, the spectra of AT 2020kog change significantly with time. They resemble those of type IIn supernovae at early phases, then they become similar to those of K-type stars during the plateau, and to M-type stars at very late phases. In contrast, AT 2020hat already shows a redder continuum at early epochs, and its spectrum shows the late appearance of molecular bands. A moderate-resolution spectrum of AT 2020hat taken at +37 d after maximum shows a forest of narrow P Cygni lines of metals with velocities of 180 km s−1, along with an Hα emission with a full-width at half-maximum velocity of 250 km s−1. For AT 2020hat, a robust constraint on its quiescent progenitor is provided by archival images of the Hubble Space Telescope. The progenitor is clearly detected as a mid-K type star, with an absolute magnitude of MF606W = −3.33 ± 0.09 mag and a colour of F606W − F814W = 1.14 ± 0.05 mag, which are inconsistent with the expectations from a massive star that could later produce a core-collapse supernova. Although quite peculiar, the two objects nicely match the progenitor versus light curve absolute magnitude correlations discussed in the literature.

We present the discovery and high-cadence follow-up observations of SN2018ivc, an unusual SNe II that exploded in NGC1068 (D=10.1 Mpc). The light curve of SN2018ivc declines piecewise-linearly, changing slope frequently, with four clear slope changes in the first 30 days of evolution. This rapidly changing light curve indicates that interaction between the circumstellar material and ejecta plays a significant role in the evolution. Circumstellar interaction is further supported by a strong X-ray detection. The spectra are rapidly evolving and dominated by hydrogen, helium, and calcium emission lines. We identify a rare high-velocity emission-line feature blueshifted at ∼7800 km s 1 (in Hα, Hβ, Pβ, Pγ, He I, and Ca II), which is visible from day 18 until at least day 78 and could be evidence of an asymmetric progenitor or explosion. From the overall similarity between SN2018ivc and SN1996al, the Hα equivalent width of its parent H IIregion, and constraints from pre-explosion archival Hubble Space Telescope images, we find that the progenitor of SN2018ivc could be as massive as 52 M  but is more likely <12 M . SN2018ivc demonstrates the importance of the early discovery and rapid follow-up observations of nearby supernovae to study the physics and progenitors of these cosmic explosions.

With the conclusion of the third observing run for Advanced LIGO/Virgo (O3), we present a detailed analysis of both triggered and serendipitous observations of 17 gravitational-wave (GW) events (7 triggered and 10 purely serendipitous) from the Searches After Gravitational-waves Using ARizona Observatories (SAGUARO) program. We searched a total of 4935 deg2 down to a median 5σ transient detection depth of 21.1 AB mag using the Mt. Lemmon 1.5 m telescope, the discovery engine for SAGUARO. In addition to triggered events within 24 hr, our transient search encompassed a time interval following GW events of <120 hr, providing observations on ∼1/2 of the events accessible to the Mt. Lemmon 1.5 m telescope. We covered 2.1%–86% of the LVC total probability (P total) for individual events, with a median P total ≈ 8% within <120 hr. Following improvements to our pipeline and the addition of serendipitous observations, we find a total of seven new optical candidates across five GW events, which we are unable to rule out after searching for additional information and comparing to kilonova models. Using both publicly available and our own late-time data, we investigated a total of 252 optical candidates for these 17 events, finding that only 65% were followed up in some capacity by the community. Of the total 252 candidates, we are able to rule out an additional 12 previously reported counterpart candidates. In light of these results, we discuss lessons learned from the SAGUARO GW counterpart search. We discuss how community coordination of observations and candidate follow-up, as well as the role of archival data, are crucial to improving the efficiency of follow-up efforts and preventing unnecessary duplication of effort with limited electromagnetic resources.