While unobscured and radio-quiet active galactic nuclei are regularly being found at redshifts $z > 6$ , their obscured and radio-loud counterparts remain elusive. We build upon our successful pilot study, presenting a new sample of low-frequency-selected candidate high-redshift radio galaxies (HzRGs) over a sky area 20 times larger. We have refined our selection technique, in which we select sources with curved radio spectra between 72–231 MHz from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. In combination with the requirements that our GLEAM-selected HzRG candidates have compact radio morphologies and be undetected in near-infrared $K_{\rm s}$ -band imaging from the Visible and Infrared Survey Telescope for Astronomy Kilo-degree Infrared Galaxy (VIKING) survey, we find 51 new candidate HzRGs over a sky area of approximately $1200\ \mathrm{deg}^2$ . Our sample also includes two sources from the pilot study: the second-most distant radio galaxy currently known, at $z=5.55$ , with another source potentially at $z \sim 8$ . We present our refined selection technique and analyse the properties of the sample. We model the broadband radio spectra between 74 MHz and 9 GHz by supplementing the GLEAM data with both publicly available data and new observations from the Australia Telescope Compact Array at 5.5 and 9 GHz. In addition, deep $K_{\rm s}$ -band imaging from the High-Acuity Widefield K-band Imager (HAWK-I) on the Very Large Telescope and from the Southern Herschel Astrophysical Terahertz Large Area Survey Regions $K_{\rm s}$ -band Survey (SHARKS) is presented for five sources. We discuss the prospects of finding very distant radio galaxies in our sample, potentially within the epoch of reionisation at $z \gtrsim 6.5$ .

We present Hubble Space Telescope Wide Field Camera 3 photometric and grism observations of the candidate ultra-high-redshift ( $z>7$ ) radio galaxy, GLEAM J0917–0012. This radio source was selected due to the curvature in its 70–230 MHz, low-frequency Murchison Widefield Array radio spectrum and its faintness in K-band. Follow-up spectroscopic observations of this source with the Jansky Very Large Array and Atacama Large Millimetre Array were inconclusive as to its redshift. Our F105W and F0986M imaging observations detect the host of GLEAM J0917–0012 and a companion galaxy, $\sim$ one arcsec away. The G102 grism observations reveal a single weak line in each of the spectra of the host and the companion. To help identify these lines we utilised several photometric redshift techniques including template fitting to the grism spectra, fitting the ultraviolet (UV)-to-radio photometry with galaxy templates plus a synchrotron model, fitting of the UV-to-near-infrared photometry with EAZY, and fitting the radio data alone with RAiSERed. For the host of GLEAM J0917–0012 we find a line at $1.12\,\mu$ m and the UV-to-radio spectral energy distribution (SED) fitting favours solutions at $z\sim 2$ or $z\sim 8$ . While this fitting shows a weak preference for the lower redshift solution, the models from the higher redshift solution are more consistent with the strength of the spectral line. The redshift constraint by RAiSERed of $>6.5$ also supports the interpretation that this line could be Lyman $-\alpha$ at $z=8.21$ ; however EAZY favours the $z\sim 2$ solution. We discuss the implications of both solutions. For the companion galaxy we find a line at $0.98\,\mu$ m and the SED fitting favours solutions at $z<3$ implying that the line could be the [OII] $\lambda3727$ doublet at $z=1.63$ (although the EAZY solution is $z\sim 2.6\pm 0.5$ ). Further observations are still required to unambiguously determine the redshift of this intriguing candidate ultra-high-redshift radio galaxy.