We have conducted a widefield, wideband, snapshot survey using the Australian SKA Pathfinder (ASKAP) referred to as the Rapid ASKAP Continuum Survey (RACS). RACS covers ≈ 90% of the sky, with multiple observing epochs in three frequency bands sampling the ASKAP frequency range of 700–1 800 MHz. This paper describes the third major epoch at 1 655.5 MHz, RACS-high, and the subsequent imaging and catalogue data release. The RACS-high observations at 1 655.5MHz are otherwise similar to the previously released RACS-mid (at 1367.5 MHz), and were calibrated and imaged with minimal changes. From the 1 493 images covering the sky up to declination ≈ +48°, we present a catalogue of 2 677 509 radio sources. The catalogue is constructed from images with a median root-mean-square noise of ≈ 195 μJy PSF−1 (point-spread function) and a median angular resolution of 11″. 8 × 8″. 1. The overall reliability of the catalogue is estimated to be 99.18 %, and we find a decrease in reliability as angular resolution improves. We estimate the brightness scale to be accurate to 10 %, and the astrometric accuracy to be within ≈ 0″. 6 in right ascension and ≈ 0″. 7 in declination after correction of a systematic declination-dependent offset. All data products from RACS-high, including calibrated visibility datasets, images from individual observations, full-sensitivity mosaics, and the all-sky catalogue are available at the CSIRO ASKAP Science Data Archive.

Next generations of radio surveys are expected to identify tens of millions of new sources and identifying and classifying their morphologies will require novel and more efficient methods. Self-organising maps (SOMs), a type of unsupervised machine learning, can be used to address this problem. We map 251 259 multi-Gaussian sources from Rapid ASKAP Continuum Survey (RACS) onto a SOM with discrete neurons. Similarity metrics, such as Euclidean distances, can be used to identify the best-matching neuron or unit (BMU) for each input image. We establish a reliability threshold by visually inspecting a subset of input images and their corresponding BMU. We label the individual neurons based on observed morphologies, and these labels are included in our value-added catalogue of RACS sources. Sources for which the Euclidean distance to their BMU is $\lesssim$5 (accounting for approximately 79$\%$ of sources) have an estimated $ \gt $90% reliability for their SOM-derived morphological labels. This reliability falls to less than 70$\%$ at Euclidean distances $\gtrsim$7. Beyond this threshold it is unlikely that the morphological label will accurately describe a given source. Our catalogue of complex radio sources from RACS with their SOM-derived morphological labels from this work will be made publicly available.

There is growing evidence that the broadband radio spectral energy distributions (SEDs) of star-forming galaxies (SFGs) contain a wealth of complex physics. In this paper we aim to determine the physical emission and loss processes causing radio SED curvature and steepening to see what observed global astrophysical properties, if any, are correlated with radio SED complexity. To do this, we have acquired radio continuum data between 70 MHz and 17 GHz for a sample of 19 southern local ($z \lt 0.04$) SFGs. Of this sample 11 are selected to contain low-frequency ($ \lt $300 MHz) turnovers (LFTOs) in their SEDs and eight are control galaxies with similar global properties. We model the radio SEDs for our sample using a Bayesian framework whereby radio emission (synchrotron and free-free) and absorption or loss processes are included modularly. We find that without the inclusion of higher frequency data ($ \gt $17 GHz) single synchrotron power-law based models are always preferred for our sample; however, additional processes including free-free absorption (FFA) and synchrotron losses are often required to accurately model radio SED complexity in SFGs. The fitted synchrotron spectral indices range from $-0.45$ to $-1.07$ and are strongly anticorrelated with stellar mass suggesting that synchrotron losses are the dominant mechanism acting to steepen the spectral index in larger/more massive nearby SFGs. We find that LFTOs in the radio SED are independent from the inclination of SFGs; however, higher inclination galaxies tend to have steeper fitted spectral indices indicating losses to diffusion of cosmic ray electrons into the galactic halo. Four of five of the merging systems in our SFG sample have elevated specific star formation rates and flatter fitted spectral indices with unconstrained LFTOs. Lastly, we find no significant separation in global properties between SFGs with or without modelled LFTOs. Overall these results suggest that LFTOs are likely caused by a combination of FFA and ionisation losses in individual recent starburst regions with specific orientations and interstellar medium properties that, when averaged over the entire galaxy, do not correlate with global astrophysical properties.

Fast radio bursts (FRBs) are millisecond-duration radio waves from the Universe. Even though more than 50 physical models have been proposed, the origin and physical mechanism of FRB emissions are still unknown. The classification of FRBs is one of the primary approaches to understanding their mechanisms, but previous studies classified conventionally using only a few observational parameters, such as fluence and duration, which might be incomplete. To overcome this problem, we use an unsupervised machine-learning model, the Uniform Manifold Approximation and Projection to handle seven parameters simultaneously, including amplitude, linear temporal drift, time duration, central frequency, bandwidth, scaled energy, and fluence. We test the method for homogeneous 977 sub-bursts of FRB 20121102A detected by the Arecibo telescope. Our machine-learning analysis identified five distinct clusters, suggesting the possible existence of multiple different physical mechanisms responsible for the observed FRBs from the FRB 20121102A source. The geometry of the emission region and the propagation effect of FRB signals could also make such distinct clusters. This research will be a benchmark for future FRB classifications when dedicated radio telescopes such as the square kilometer array or Bustling Universe Radio Survey Telescope in Taiwan discover more FRBs than before.

We present deep near-infrared $K_\textrm{s}$-band imaging for 35 of the 53 sources from the high-redshift ($z \gt 2$) radio galaxy candidate sample defined in Broderick et al. (2022, PASA, 39, e061). These images were obtained using the High-Acuity Widefield K-band Imager (HAWK-I) on the Very Large Telescope. Host galaxies are detected for 27 of the sources, with $K_\textrm{s} \approx 21.6$–23.0 mag (2$^{\prime\prime}$ diameter apertures; AB). The remaining eight targets are not detected to a median $3\unicode{x03C3}$ depth of $K_\textrm{s} \approx 23.3$ mag (2$^{\prime\prime}$ diameter apertures). We examine the radio and near-infrared flux densities of the 35 sources, comparing them to the known $z \gt 3$ powerful radio galaxies with 500-MHz radio luminosities $L_{500\,\textrm{MHz}} \gt 10^{27}$ W Hz$^{-1}$. By plotting 150-MHz flux density versus $K_\textrm{s}$-band flux density, we find that, similar to the sources from the literature, these new targets have large radio to near-infrared flux density ratios, but extending the distribution to fainter flux densities. Five of the eight HAWK-I deep non-detections have a median $3\unicode{x03C3}$ lower limit of $K_\textrm{s} \gtrsim 23.8$ mag (1$.\!^{\prime\prime}$5 diameter apertures); these five targets, along with a further source from Broderick et al. (2022, PASA, 39, e061) with a deep non-detection ($K_\textrm{s} \gtrsim 23.7$ mag; $3\unicode{x03C3}$; 2$^{\prime\prime}$ diameter aperture) in the Southern H-ATLAS Regions $K_\textrm{s}$-band Survey, are considered candidates to be ultra-high-redshift ($z \gt 5$) radio galaxies. The extreme radio to near-infrared flux density ratios ($\gt 10^5$) for these six sources are comparable to TN J0924$-$2201, GLEAM J0856$+$0223 and TGSS J1530$+$1049, the three known powerful radio galaxies at $z \gt 5$. For a selection of galaxy templates with different stellar masses, we show that $z \gtrsim 4.2$ is a plausible scenario for our ultra-high-redshift candidates if the stellar mass $M_\textrm{*} \gtrsim 10^{10.5}$ M$_\odot$. In general, the 35 targets studied have properties consistent with the previously known class of infrared-faint radio sources. We also discuss the prospects for finding more UHzRG candidates from wide and deep near-infrared surveys.

We present radio observations of the galaxy cluster Abell S1136 at 888 MHz, using the Australian Square Kilometre Array Pathfinder radio telescope, as part of the Evolutionary Map of the Universe Early Science program. We compare these findings with data from the Murchison Widefield Array, XMM-Newton, the Wide-field Infrared Survey Explorer, the Digitised Sky Survey, and the Australia Telescope Compact Array. Our analysis shows the X-ray and radio emission in Abell S1136 are closely aligned and centered on the Brightest Cluster Galaxy, while the X-ray temperature profile shows a relaxed cluster with no evidence of a cool core. We find that the diffuse radio emission in the centre of the cluster shows more structure than seen in previous low-resolution observations of this source, which appeared formerly as an amorphous radio blob, similar in appearance to a radio halo; our observations show the diffuse emission in the Abell S1136 galaxy cluster contains three narrow filamentary structures visible at 888 MHz, between $\sim$80 and 140 kpc in length; however, the properties of the diffuse emission do not fully match that of a radio (mini-)halo or (fossil) tailed radio source.

We present source detection and catalogue construction pipelines to build the first catalogue of radio galaxies from the 270 $\rm deg^2$ pilot survey of the Evolutionary Map of the Universe (EMU-PS) conducted with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The detection pipeline uses Gal-DINO computer vision networks (Gupta et al. 2024, PASA, 41, e001) to predict the categories of radio morphology and bounding boxes for radio sources, as well as their potential infrared host positions. The Gal-DINO network is trained and evaluated on approximately 5 000 visually inspected radio galaxies and their infrared hosts, encompassing both compact and extended radio morphologies. We find that the Intersection over Union (IoU) for the predicted and ground-truth bounding boxes is larger than 0.5 for 99% of the radio sources, and 98% of predicted host positions are within $3^{\prime \prime}$ of the ground-truth infrared host in the evaluation set. The catalogue construction pipeline uses the predictions of the trained network on the radio and infrared image cutouts based on the catalogue of radio components identified using the Selavy source finder algorithm. Confidence scores of the predictions are then used to prioritise Selavy components with higher scores and incorporate them first into the catalogue. This results in identifications for a total of 211 625 radio sources, with 201 211 classified as compact and unresolved. The remaining 10 414 are categorised as extended radio morphologies, including 582 FR-I, 5 602 FR-II, 1 494 FR-x (uncertain whether FR-I or FR-II), 2 375 R (single-peak resolved) radio galaxies, and 361 with peculiar and other rare morphologies. Each source in the catalogue includes a confidence score. We cross-match the radio sources in the catalogue with the infrared and optical catalogues, finding infrared cross-matches for 73% and photometric redshifts for 36% of the radio galaxies. The EMU-PS catalogue and the detection pipelines presented here will be used towards constructing catalogues for the main EMU survey covering the full southern sky.

Odd Radio Circles (ORCs) are a class of low surface brightness, circular objects approximately one arcminute in diameter. ORCs were recently discovered in the Australian Square Kilometre Array Pathfinder (ASKAP) data and subsequently confirmed with follow-up observations on other instruments, yet their origins remain uncertain. In this paper, we suggest that ORCs could be remnant lobes of powerful radio galaxies, re-energised by the passage of a shock. Using relativistic hydrodynamic simulations with synchrotron emission calculated in post-processing, we show that buoyant evolution of remnant radio lobes is alone too slow to produce the observed ORC morphology. However, the passage of a shock can produce both filled and edge-brightnened ORC-like morphologies for a wide variety of shock and observing orientations. Circular ORCs are predicted to have host galaxies near the geometric centre of the radio emission, consistent with observations of these objects. Significantly offset hosts are possible for elliptical ORCs, potentially causing challenges for accurate host galaxy identification. Observed ORC number counts are broadly consistent with a paradigm in which moderately powerful radio galaxies are their progenitors.

All very massive early-type galaxies contain supermassive blackholes, but are these blackholes all sufficiently active to produce detectable radio continuum sources? We have used the 887.5 MHz Rapid ASKAP Continuum Survey DR1 to measure the radio emission from morphological early-type galaxies brighter than $K_S=9.5$ selected from the 2MASS Redshift Survey, HyperLEDA, and RC3. In line with previous studies, we find median radio power increases with infrared luminosity, with $P_{1.4} \propto L_K^{2.2}$, although the scatter about this relation spans several orders of magnitude. All 40 of the $M_K<-25.7$ early-type galaxies in our sample have measured radio flux densities that are more than $2\sigma$ above the background noise, with $1.4\,{\rm GHz}$ radio powers spanning ${\sim} 3 \times 10^{20}$ to ${\sim} 3\times 10^{25}\,{\rm W/Hz^{-1}}$. Cross-matching our sample with integral field spectroscopy of early-type galaxies reveals that the most powerful radio sources preferentially reside in galaxies with relatively low angular momentum (i.e. slow rotators). While the infrared colours of most galaxies in our early-type sample are consistent with passive galaxies with negligible star formation and the radio emission produced by active galactic nuclei or AGN remnants, very low levels of star formation could power the weakest radio sources with little effect on many other star formation rate tracers.

The Australian SKA Pathfinder (ASKAP) has surveyed the sky at multiple frequencies as part of the Rapid ASKAP Continuum Survey (RACS). The first two RACS observing epochs, at 887.5 (RACS-low) and 1 367.5 (RACS-mid) MHz, have been released (McConnell, et al. 2020, PASA, 37, e048; Duchesne, et al. 2023, PASA, 40, e034). A catalogue of radio sources from RACS-low has also been released, covering the sky south of declination $+30^{\circ}$ (Hale, et al., 2021, PASA, 38, e058). With this paper, we describe and release the first set of catalogues from RACS-mid, covering the sky below declination $+49^{\circ}$. The catalogues are created in a similar manner to the RACS-low catalogue, and we discuss this process and highlight additional changes. The general purpose primary catalogue covering 36 200 deg$^2$ features a variable angular resolution to maximise sensitivity and sky coverage across the catalogued area, with a median angular resolution of $11.2^{\prime\prime} \times 9.3^{\prime\prime}$. The primary catalogue comprises 3 105 668 radio sources, including those in the Galactic Plane (2 861 923 excluding Galactic latitudes of $|b|<5^{\circ}$), and we estimate the catalogue to be 95% complete for sources above 2 mJy. With the primary catalogue, we also provide two auxiliary catalogues. The first is a fixed-resolution, 25-arcsec catalogue approximately matching the sky coverage of the RACS-low catalogue. This 25-arcsec catalogue is constructed identically to the primary catalogue, except images are convolved to a less-sensitive 25-arcsec angular resolution. The second auxiliary catalogue is designed for time-domain science and is the concatenation of source lists from the original RACS-mid images with no additional convolution, mosaicking, or de-duplication of source entries to avoid losing time-variable signals. All three RACS-mid catalogues, and all RACS data products, are available through the CSIRO ASKAP Science Data Archive (https://research.csiro.au/casda/).

Creating radio galaxy catalogues from next-generation deep surveys requires automated identification of associated components of extended sources and their corresponding infrared hosts. In this paper, we introduce RadioGalaxyNET, a multimodal dataset, and a suite of novel computer vision algorithms designed to automate the detection and localization of multi-component extended radio galaxies and their corresponding infrared hosts. The dataset comprises 4 155 instances of galaxies in 2 800 images with both radio and infrared channels. Each instance provides information about the extended radio galaxy class, its corresponding bounding box encompassing all components, the pixel-level segmentation mask, and the keypoint position of its corresponding infrared host galaxy. RadioGalaxyNET is the first dataset to include images from the highly sensitive Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope, corresponding infrared images, and instance-level annotations for galaxy detection. We benchmark several object detection algorithms on the dataset and propose a novel multimodal approach to simultaneously detect radio galaxies and the positions of infrared hosts.

We here present 0.02–0.04″ resolution ALMA observation of the compact obscured nucleus (CON) of IRAS 17578-0400. A dusty torus within the nucleus, approximately 4 pc in radius, has been uncovered, exhibiting a usually flat spectral index at ALMA band 3, likely due to the millimeter corona emission from the central supermassive black hole (SMBH). The dense gas disk, traced by 13CO(1-0), spans 7 pc in radius and suggests an outflow driven by a disk wind due to its asymmetrical structure along the minor axis. Collimated molecular outflows (CMO), traced by the low-velocity components of the HCN(3-2) and HCO+(3-2) lines, align with the minor axis gas disk. Examination of position-velocity plots of HCN(3-2) and HCO+(3-2) reveals a flared dense gas disk extended a radius of ∼60 pc, infalling and rotating at speeds of about 200 km s−1 and 300 km s−1 respectively. A centrifugal barrier, located around 4 pc from the dynamical center, implies an SMBH mass of approximately 108 Mȯ, consistent with millimeter corona emission estimates. The CMO maintains a steady rotation speed of 200 km s−1 over the 100 pc scale along the minor axis. The projected speed of the CMO is about 80 km s−1, corresponding to around ∼500 km s−1, assuming an inclination angle of 80°. Such a kinematics structure of disk-driven collimated rotating molecular outflow with gas supplies from a falling rotating disk indicates that the feedback of the compact obscured nucleus is likely regulated by the momentum transfer of the molecular gas that connects to both the feeding of the nuclear starburst and supermassive black hole.

The present work discusses the use of a weakly-supervised deep learning algorithm that reduces the cost of labelling pixel-level masks for complex radio galaxies with multiple components. The algorithm is trained on weak class-level labels of radio galaxies to get class activation maps (CAMs). The CAMs are further refined using an inter-pixel relations network (IRNet) to get instance segmentation masks over radio galaxies and the positions of their infrared hosts. We use data from the Australian Square Kilometre Array Pathfinder (ASKAP) telescope, specifically the Evolutionary Map of the Universe (EMU) Pilot Survey, which covered a sky area of 270 square degrees with an RMS sensitivity of 25–35 $\mu$Jy beam$^{-1}$. We demonstrate that weakly-supervised deep learning algorithms can achieve high accuracy in predicting pixel-level information, including masks for the extended radio emission encapsulating all galaxy components and the positions of the infrared host galaxies. We evaluate the performance of our method using mean Average Precision (mAP) across multiple classes at a standard intersection over union (IoU) threshold of 0.5. We show that the model achieves a mAP$_{50}$ of 67.5% and 76.8% for radio masks and infrared host positions, respectively. The network architecture can be found at the following link: https://github.com/Nikhel1/Gal-CAM

The Australian SKA Pathfinder (ASKAP) radio telescope has carried out a survey of the entire Southern Sky at 887.5 MHz. The wide area, high angular resolution, and broad bandwidth provided by the low-band Rapid ASKAP Continuum Survey (RACS-low) allow the production of a next-generation rotation measure (RM) grid across the entire Southern Sky. Here we introduce this project as Spectral and Polarisation in Cutouts of Extragalactic sources from RACS (SPICE-RACS). In our first data release, we image 30 RACS-low fields in Stokes I, Q, U at 25$^{\prime\prime}$ angular resolution, across 744–1032 MHz with 1 MHz spectral resolution. Using a bespoke, highly parallelised, software pipeline we are able to rapidly process wide-area spectro-polarimetric ASKAP observations. Notably, we use ‘postage stamp’ cutouts to assess the polarisation properties of 105912 radio components detected in total intensity. We find that our Stokes Q and U images have an rms noise of $\sim$80 $\unicode{x03BC}$Jy PSF$^{-1}$, and our correction for instrumental polarisation leakage allows us to characterise components with $\gtrsim$1% polarisation fraction over most of the field of view. We produce a broadband polarised radio component catalogue that contains 5818 RM measurements over an area of $\sim$1300 deg$^{2}$ with an average error in RM of $1.6^{+1.1}_{-1.0}$ rad m$^{-2}$, and an average linear polarisation fraction $3.4^{+3.0}_{-1.6}$ %. We determine this subset of components using the conditions that the polarised signal-to-noise ratio is $>$8, the polarisation fraction is above our estimated polarised leakage, and the Stokes I spectrum has a reliable model. Our catalogue provides an areal density of $4\pm2$ RMs deg$^{-2}$; an increase of $\sim$4 times over the previous state-of-the-art (Taylor, Stil, Sunstrum 2009, ApJ, 702, 1230). Meaning that, having used just 3% of the RACS-low sky area, we have produced the 3rd largest RM catalogue to date. This catalogue has broad applications for studying astrophysical magnetic fields; notably revealing remarkable structure in the Galactic RM sky. We will explore this Galactic structure in a follow-up paper. We will also apply the techniques described here to produce an all-Southern-sky RM catalogue from RACS observations. Finally, we make our catalogue, spectra, images, and processing pipeline publicly available.

We present the Cosmological Double Radio Active Galactic Nuclei (CosmoDRAGoN) project: a large suite of simulated AGN jets in cosmological environments. These environments sample the intra-cluster media of galaxy clusters that form in cosmological smooth particle hydrodynamics (SPH) simulations, which we then use as inputs for grid-based hydrodynamic simulations of radio jets. Initially conical jets are injected with a range of jet powers, speeds (both relativistic and non-relativistic), and opening angles; we follow their collimation and propagation on scales of tens to hundreds of kiloparsecs, and calculate spatially resolved synthetic radio spectra in post-processing. In this paper, we present a technical overview of the project, and key early science results from six representative simulations which produce radio sources with both core- (Fanaroff-Riley Type I) and edge-brightened (Fanaroff-Riley Type II) radio morphologies. Our simulations highlight the importance of accurate representation of both jets and environments for radio morphology, radio spectra, and feedback the jets provide to their surroundings.

We present the highest resolution and sensitivity $\sim$$1.4\,$GHz continuum observations of the Eridanus supergroup obtained as a part of the Widefield Australian Square Kilometer Array Pathfinder (ASKAP) L-band Legacy All-sky Blind surveY (WALLABY) pre-pilot observations using the ASKAP. We detect 9461 sources at 1.37 GHz down to a flux density limit of $\sim$$0.1$ mJy at $6.1''\times 7.9''$ resolution with a median root mean square of 0.05 mJy beam$^{-1}$. We find that the flux scale is accurate to within 5 % (compared to NVSS at 1.4 GHz). We then determine the global properties of eight Eridanus supergroup members, which are detected in both radio continuum and neutral hydrogen (HI) emission, and find that the radio-derived star formation rates (SFRs) agree well with previous literature. Using our global and resolved radio continuum properties of the nearby Eridanus galaxies, we measure and extend the infrared-radio correlation (IRRC) to lower stellar masses and inferred SFRs than before. We find the resolved IRRC to be useful for: (1) discriminating between active galactic nuclei and star-forming galaxies; (2) identifying background radio sources; and (3) tracing the effects of group environment pre-processing in NGC 1385. We find evidence for tidal interactions and ram-pressure stripping in the HI, resolved spectral index and IRRC morphologies of NGC 1385. There appears to be a spatial coincidence (in projection) of double-lobed radio jets with the central HI hole of NGC 1367. The destruction of polycyclic aromatic hydrocarbons by merger-induced shocks may be driving the observed WISE W3 deficit observed in NGC 1359. Our results suggest that resolved radio continuum and IRRC studies are excellent tracers of the physical processes that drive galaxy evolution and will be possible on larger sample of sources with upcoming ASKAP radio continuum surveys.

Spectral variability offers a new technique to identify small scale structures from scintillation, as well as determining the absorption mechanism for peaked-spectrum (PS) radio sources. In this paper, we present very long baseline interferometry (VLBI) imaging using the long baseline array (LBA) of two PS sources, MRC 0225–065 and PMN J0322–4820, identified as spectrally variable from observations with the Murchison Widefield Array (MWA). We compare expected milliarcsecond structures based on the detected spectral variability with direct LBA imaging. We find MRC 0225–065 is resolved into three components, a bright core and two fainter lobes, roughly 430 pc projected separation. A comprehensive analysis of the magnetic field, host galaxy properties, and spectral analysis implies that MRC 0225–065 is a young radio source with recent jet activity over the last $10^2$–$10^3$ yr. We find PMN J0322–4820 is unresolved on milliarcsecond scales. We conclude PMN J0322–4820 is a blazar with flaring activity detected in 2014 with the MWA. We use spectral variability to predict morphology and find these predictions consistent with the structures revealed by our LBA images.

Investigating the thermal and non-thermal processes in galaxies is vital to understand their evolution over cosmic time. This can best be studied by combining the radio and optical/near-infrared observations of galaxies. The JWST can resolve the evolution of the thermal processes by mapping ionized gas and dust in distant galaxies. This information combined with the upcoming surveys with the Square Kilometer Array (SKA) will make a major breakthrough in mapping the non-thermal processes and understanding their role in the evolution of galaxies. Our simulations show that SKA surveys will be able to trace the evolution history of spiral galaxies such as M 51 and NGC 6946 back to a redshift of 3 already in its first phase of construction. This study indicates the important role of the non-thermal pressure inserted by cosmic rays and magnetic fields in deriving winds and outflows at cosmic noon as deduced by a flat synchrotron spectrum in star forming galaxies.

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 a set of peculiar radio sources detected using an unsupervised machine learning method. We use data from the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to train a self-organizing map (SOM). The radio maps from three ASKAP surveys, Evolutionary Map of Universe pilot survey (EMU-PS), Deep Investigation of Neutral Gas Origins pilot survey (DINGO), and Survey With ASKAP of GAMA-09 + X-ray (SWAG-X), are used to search for the rarest or unknown radio morphologies. We use an extension of the SOM algorithm that implements rotation and flipping invariance on astronomical sources. The SOM is trained using the images of all ‘complex’ radio sources in the EMU-PS which we define as all sources catalogued as ‘multi-component’. The trained SOM is then used to estimate a similarity score for complex sources in all surveys. We select 0.5% of the sources that are most complex according to the similarity metric and visually examine them to find the rarest radio morphologies. Among these, we find two new odd radio circle (ORC) candidates and five other peculiar morphologies. We discuss multiwavelength properties and the optical/infrared counterparts of selected peculiar sources. In addition, we present examples of conventional radio morphologies including: diffuse emission from galaxy clusters, and resolved, bent-tailed, and FR-I and FR-II type radio galaxies. We discuss the overdense environment that may be the reason behind the circular shape of ORC candidates.