Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

We report the discovery and panchromatic followup observations of the young Type Ic supernova, SN2020oi, in M100, a grand design spiral galaxy at a mere distance of $14$ Mpc. We followed up with observations at radio, X-ray and optical wavelengths from only a few days to several months after explosion. The optical behaviour of the supernova is similar to those of other normal Type Ic supernovae. The event was not detected in the X-ray band but our radio observation revealed a bright mJy source ($L_{\nu} \approx 1.2 \times 10^{27} {\rm erg\,s}^{-1} {\rm Hz}^{-1}$). Given, the relatively small number of stripped envelope SNe for which radio emission is detectable, we used this opportunity to perform a detailed analysis of the comprehensive radio dataset we obtained. The radio emitting electrons initially experience a phase of inverse Compton cooling which leads to steepening of the spectral index of the radio emission. Our analysis of the cooling frequency points to a large deviation from equipartition at the level of $\epsilon_e/\epsilon_B \gtrsim 200$, similar to a few other cases of stripped envelope SNe. Our modeling of the radio data suggests that the shockwave driven by the SN ejecta into the circumstellar matter (CSM) is moving at $\sim 3\times 10^{4}\,{\rm km\,s}^{-1}$. Assuming a constant mass-loss from the stellar progenitor, we find that the mass-loss rate is $\dot{M} \approx 1.4\times 10^{-4}\,{M}_{\odot}\,{\rm yr}^{-1}$, for an assumed wind velocity of $1000\,{\rm km\,s}^{-1}$. The temporal evolution of the radio emission suggests a radial CSM density structure steeper than the standard $r^{-2}$.

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 the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

The jet opening angle and inclination of GW170817 -- the first detected binary neutron star merger -- were vital to understand its energetics, relation to short gamma-ray bursts, and refinement of the standard siren-based determination of the Hubble constant, $H_0$. These basic quantities were determined through a combination of the radio lightcurve and Very Long Baseline Interferometry (VLBI) measurements of proper motion. In this paper we discuss and quantify the prospects for the use of radio VLBI observations and observations of scintillation-induced variability to measure the source size and proper motion of merger afterglows, and thereby infer properties of the merger including inclination angle, opening angle and energetics. We show that these techniques are complementary as they probe different parts of the circum-merger density/inclination angle parameter space and different periods of the temporal evolution of the afterglow. We also find that while VLBI observations will be limited to the very closest events it will be possible to detect scintillation for a large fraction of events beyond the range of current gravitational wave detectors. Scintillation will also be detectable with next generation telescopes such as the Square Kilometre Array, 2000 antenna Deep Synoptic Array and the next generation Very Large Array, for a large fraction of events detected with third generation gravitational wave detectors. Finally, we discuss prospects for the measurement of the $H_0$ with VLBI observations of neutron star mergers and compare this technique to other standard siren methods.

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.

The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

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.

We present detailed submillimeter- through centimeter-wave observations of the extraordinary extragalactic transient AT2018cow. The apparent characteristics -- the high radio luminosity, the long-lived emission plateau at millimeter bands, and the sub-relativistic velocity -- have no precedent. A basic interpretation of the data suggests $E_k \gtrsim 10^{48}$ erg coupled to a fast but sub-relativistic ($v \approx 0.13c$) shock in a dense ($n_e \approx 3 \times 10^{5}\,$cm$^{-3}$) medium. We find that the X-ray emission is not naturally explained by an extension of the radio-submm synchrotron spectrum, nor by inverse Compton scattering of the dominant blackbody UVOIR photons by energetic electrons within the forward shock. By $\Delta t \approx20\,$days, the X-ray emission shows spectral softening and erratic inter-day variability. Taken together, we are led to invoke an additional source of X-ray emission: the central engine of the event. Regardless of the nature of this central engine, this source heralds a new class of energetic transients shocking a dense medium, which at early times are most readily observed at millimeter wavelengths.

We identified a highly-polarized, steep-spectrum radio source in a deep image with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope at 888 MHz. After considering and rejecting a stellar origin for this source, we discovered a new millisecond pulsar (MSP) using observations from the Parkes radio telescope. This pulsar has period 2.77 ms and dispersion measure 228.27 pc/cm**3. Although this pulsar does not yet appear to be particularly remarkable, the short spin period, wide profile and high dispersion measure do make it relatively hard to discover through traditional blind periodicity searches. Over the course of several weeks we see changes in the barycentric period of this pulsar that are consistent with orbital motion in a binary system, but the properties of any binary need to be confirmed by further observations. While even a deep ASKAP survey may not identify large numbers of new MSPs compared to the existing population, it would be competitive with existing all-sky surveys and could discover interesting new MSPs at high Galactic latitude without the need for computationally-expensive all-sky periodicity searches.

We revisit the K -band period-luminosity (P-L) relations of Galactic red supergiants using Gaia Data Release 2 parallaxes and up to 70 yr of photometry from AAVSO and ASAS campaigns. In addition, we examine 206 LMC red supergiants using 50 yr of photometric data from the Digitised Harvard Astronomical Plate Collection. We identified periods by computing power spectra and calculated the period-luminosity relations of our samples and compared them with the literature. Newly available data tighten the P-L relations substantially. Identified periods form two groups: one with periods of 300-1000 days, corresponding to pulsations, and another with Long Secondary Periods between 1000 and 8000 days. Among the 48 Galactic objects we find shorter periods in 25 stars and long secondary periods in 23 stars. In the LMC sample we identify 85 and 94 red supergiants with shorter and long secondary periods, respectively. The P-L relation of the Galactic red supergiants is in agreement with the red supergiants in both, the Large Magellanic Cloud and the Andromeda Galaxy. We and no clear continuity between the known red giant period-luminosity sequences, and the red supergiant sequences investigated here.

The detection of a neutron star merger by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) and Advanced Virgo gravitational wave detectors and the subsequent detection of an electromagnetic counterpart has opened a new era of transient astronomy. With upgrades to the Advanced LIGO and Advanced Virgo detectors and new detectors coming online in Japan and India, neutron star mergers will be detected at a higher rate in the future, starting with the O3 observing run which will begin in early 2019. The detection of electromagnetic emission from these mergers provides vital information about merger parameters and allows independent measurement of the Hubble constant. The Australian Square Kilometre Array Pathfinder (ASKAP) is expected to become fully operational early 2019 and its 30 deg$^2$ field of view will enable us to rapidly survey large areas of sky. In this work we explore prospects for detecting both prompt and long-term radio emission from neutron star mergers with ASKAP and determine an observing strategy that optimises the use of telescope time. We investigate different strategies to tile the sky with telescope pointings in order to detect radio counterparts with limited observing time, using 475 simulated gravitational wave events. Our results show a significant improvement in observing efficiency when compared with a na\"ive strategy of covering the entire localisation above some confidence threshold, even when achieving the same total probability covered.

We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817.

We present 2-9 GHz radio observations of GW170817 covering the period 125-200 days post-merger, taken with the Australia Telescope Compact Array and the Karl G. Jansky Very Large Array. Our observations demonstrate that the radio afterglow peaked at $149\pm2$ days post-merger and is now declining in flux density. We see no evidence for evolution in the radio-only spectral index, which remains consistent with optically-thin synchrotron emission connecting the radio, optical, and X-ray regimes. The peak implies a total energy in the synchrotron-emitting component of a ${\rm few}\times 10^{50}\,$erg. The temporal decay rate is most consistent with mildly- or non-relativistic material and we do not see evidence for a very energetic off-axis jet, but we cannot distinguish between a lower-energy jet and more isotropic emission.

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

Gravitational waves have been detected from a binary neutron star merger event, GW170817. The detection of electromagnetic radiation from the same source has shown that the merger occurred in the outskirts of the galaxy NGC 4993, at a distance of 40 megaparsecs from Earth. We report the detection of a counterpart radio source that appears 16 days after the event, allowing us to diagnose the energetics and environment of the merger. The observed radio emission can be explained by either a collimated ultra-relativistic jet viewed off-axis, or a cocoon of mildly relativistic ejecta. Within 100 days of the merger, the radio light curves will distinguish between these models and very long baseline interferometry will have the capability to directly measure the angular velocity and geometry of the debris.

We present low-frequency spectral energy distributions of 60 known radio pulsars observed with the Murchison Widefield Array (MWA) telescope. We searched the GaLactic and Extragalactic All-sky MWA (GLEAM) survey images for 200-MHz continuum radio emission at the position of all pulsars in the ATNF pulsar catalogue. For the 60 confirmed detections we have measured flux densities in 20 x 8 MHz bands between 72 and 231 MHz. We compare our results to existing measurements and show that the MWA flux densities are in good agreement.

We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < {\delta} < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar scintillation and variability on timescales of greater than 2 minutes. One further pulsar, PSR J0034-0721, showed significant variability, the physical origins of which are difficult to determine. The dynamic spectra for PSR J0953+0755 and PSR J0437-4715 show discrete time and frequency structure consistent with diffractive interstellar scintillation and we present the scintillation bandwidth and timescales from these observations. The remaining pulsars within our sample were statistically non-variable. We also explore the spectral properties of this sample and find spectral curvature in pulsars PSR J0835-4510, PSR J1752-2806 and PSR J0437-4715.