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.

We investigate the prospects of detecting radio afterglows from long Gamma-Ray Bursts (GRBs) from Population III (Pop III) progenitors using the SKA precursor instruments WMA (Murchison Widefield Array) and ASKAP (Australian SKA Pathfinder). We derive a realistic model of GRB afterglows that encompasses the widest range of plausible physical parameters and observation angles. We define the best case scenario of Pop III GRB energy and redshift distributions. Using probability distribution functions fitted to the observed microphysical parameters of long GRBs, we simulate a large number of Pop III GRB afterglows to find the global probability of detection. We find that ASKAP may be able to detect 35% of Pop III GRB afterglows in the optimistic case, and 27% in the pessimistic case. A negligible number will be detectable by MWA in either case. Detections per image for ASKAP, found by incorporating intrinsic rates with detectable timescales, are as high as $\sim$ 6000 and as low as $\sim$ 11, which shows the optimistic case is unrealistic. We track how the afterglow flux density changes over various time intervals and find that, because of their very slow variability, the cadence for blind searches of these afterglows should be as long as possible. We also find Pop III GRBs at high redshift have radio afterglow lightcurves that are indistinguishable from those of regular long GRBs in the more local universe.

Gamma ray burst (GRB) optical and X-ray afterglow luminosity is expected to correlate with the GRB isotropic equivalent kinetic energy of the outflow in the standard synchrotron model for GRB afterglows. Previous studies, using prompt GRB isotropic equivalent energy ($E_{\rm iso}$) as a proxy for isotropic equivalent kinetic energy, have generally confirmed a correlation between X-ray and optical afterglow luminosities. Assuming that GRB afterglow luminosity does not evolve strongly with redshift, we identify a strong Malmquist bias in GRB optical and X-ray afterglow luminosity data. We show that selection effects dominate the observed E$_{\rm iso}$ -- L$_{\rm opt,X}$ correlations, and have likely been underestimated in other studies. The bias is strongest for a subset of optically faint bursts $m>24$ at 24 hr with $z>2$. After removing this optical selection bias, the E$_{\rm iso}$ -- L$_{\rm opt,X}$ correlation for long GRBs is not statistically significant, but combining both long and short GRB luminosity data the correlation is significant. Using the median of the $E_{\rm iso}$ and $L_{\rm opt,X}$ distributions, we apply the synchrotron model assuming the same power law index for short and long GRBs, but different microphysical parameter distributions. Comparing the ratio of optical and X-ray luminosities, we find tentative evidence that the fraction of post-shock energy in magnetic fields, $\epsilon_B$, could be systematically higher in SGRBs compared to LGRBs.

We investigate the detectability of a proposed population of Gamma-Ray Bursts (GRBs) from the collapse of Population III (Pop III) stars. The James Webb Space Telescope (JWST) and Space Infrared telescope for Cosmology and Astrophysics (SPICA) will be able to observe the late time infrared afterglows. We have developed a new method to calculate their detectability, which takes into account the fundamental initial mass function (IMF) and formation rates of Pop III stars, from which we find the temporal variability of the afterglows and ultimately the length of time JWST and SPICA can detect them. In the range of plausible Pop III GRB parameters, the afterglows are always detectable by these instruments during the isotropic emission, for a minimum of 55 days and a maximum of 3.7 years. The average number of detectable afterglows will be 2.96$\times 10^{-5}$ per SPICA field of view (FOV) and 2.78$\times 10^{-6}$ per JWST FOV. These are lower limits, using a pessimistic estimate of Pop III star formation. An optimal observing strategy with SPICA could identify a candidate orphan afterglow in $\sim 1.3$ years, with a 90 percent probability of confirmation with further detailed observations. A beamed GRB will align with the FOV of the planned GRB detector Energetic X-ray Imaging Survey Telescope once every 9 years. Pop III GRBs will be more easily detected by their isotropic emissions (i.e. orphan afterglows) rather than by their prompt emissions.

We employ realistic constraints on astrophysical and instrumental selection effects to model the Gamma-Ray Burst (GRB) redshift distribution using {\it Swift} triggered redshift samples acquired from optical afterglows (OA) and the TOUGH survey. Models for the Malmquist bias, redshift desert, and the fraction of afterglows missing because of host galaxy dust extinction, are used to show how the "true" GRB redshift distribution is distorted to its presently observed biased distribution. We also investigate another selection effect arising from a correlation between $E_{{\rm iso}}$ and $L_{{\rm opt}}$. The analysis, which accounts for the missing fraction of redshifts in the two data subsets, shows that a combination of selection effects (both instrumental and astrophysical) can describe the observed GRB redshift distribution. Furthermore, the observed distribution is compatible with a GRB rate evolution that tracks the global SFR, although the rate at high-$z$ cannot be constrained with confidence. Taking selection effects into account, it is not necessary to invoke high-energy GRB luminosity evolution with redshift to explain the observed GRB rate at high-$z$.