Massive, rapidly-spinning magnetar remnants produced as a result of binary neutron star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a synchrotron radio signal from the interaction of the ejecta with the circumburst medium. Here, we present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (GRBs; $z<0.5$) on rest-frame timescales of $\approx 2.4-13.9$ yr following the bursts. We place $3\sigma$ limits on radio continuum emission of $F_{\nu}\lesssim 6-20\,\mu$Jy at the burst positions, or $L_{\nu} \lesssim (0.6-8.3) \times10^{28}$ erg s$^{-1}$ Hz$^{-1}$. Comparing these limits with new light curve modeling which properly incorporates relativistic effects, we obtain limits on the energy deposited into the ejecta of $E_{\rm ej}\lesssim (0.6-6.7)\times 10^{52}$ erg ($E_{\rm ej}\lesssim (1.7-14.1) \times 10^{52}$ erg) for an ejecta mass of $0.03\,M_{\odot}$ ($0.1\,M_{\odot}$). We present a uniform re-analysis of 27 short GRBs with 5.5-6.0 GHz observations, and find that $\gtrsim 50\%$ of short GRBs did not form stable magnetar remnants in their mergers. Assuming short GRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction $f_{\rm GW190425}$ of cases, we place constraints on the maximum mass of a non-rotating neutron star (NS) ($M_{\rm TOV}$), finding $M_{\rm TOV} \lesssim 2.4\,M_{\odot}$ for $f_{\rm GW190425}=0.4$. The detection (or lack thereof) of radio remnants in untargeted surveys such as the VLA Sky Survey (VLASS) could provide more stringent constraints on the fraction of mergers that produce stable remnants. If $\gtrsim 100-1000$ radio remnants are discovered in VLASS, this suggests that short GRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce.

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 $\approx 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 optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

The orientation of the magnetic field (B-field) in the filamentary dark cloud GF 9 was traced from the periphery of the cloud into the L1082C dense core that contains the low-mass, low-luminosity Class 0 young stellar object (YSO) GF 9-2 (IRAS 20503+6006). This was done using SOFIA HAWC+ dust thermal emission polarimetry (TEP) at 216 um in combination with Mimir near-infrared background starlight polarimetry (BSP) conducted at H-band (1.6 um) and K-band (2.2 um). These observations were augmented with published I-band (0.77 um) BSP and Planck 850 um TEP to probe B-field orientations with offset from the YSO in a range spanning 6000 AU to 3 pc. No strong B-field orientation change with offset was found, indicating remarkable uniformity of the B-field from the cloud edge to the YSO environs. This finding disagrees with weak-field models of cloud core and YSO formation. The continuity of inferred B-field orientations for both TEP and BSP probes is strong evidence that both are sampling a common B-field that uniformly threads the cloud, core, and YSO region. Bayesian analysis of Gaia DR2 stars matched to the Mimir BSP stars finds a distance to GF 9 of 270 +/- 10 pc. No strong wavelength dependence of B-field orientation angle was found, contrary to previous claims.