We present the discovery of two extended, low surface-brightness radio continuum sources, each consisting of a near-circular body and an extended tail of emission, nicknamed Stingray 1 (ASKAP J0129–5350) and Stingray 2 (ASKAP J0245–5642). Both are found in the direction of the Magellanic Stream (MS) and were discovered in the Australian Square Kilometre Array Pathfinder (ASKAP) Evolutionary Map of the Universe (EMU) survey at 944 MHz. We combine the ASKAP data with low-frequency radio observations from the GaLactic and Extragalactic All-sky MWA Survey (GLEAM) to conduct a radio continuum analysis. galaxy pairs or groups, and Odd Radio Circles (ORCs). We explore both Galactic/near Galactic scenarios, including runaway or circumgalactic supernova remnants (SNRs) and parentless pulsar-wind nebulae (PWNe), and extragalactic scenarios including radio active galactic nuclei (AGNs), dying radio galaxies, galaxy clusters, galaxy pairs or groups, head-tail radio galaxies, and ORCs, as well as the possibility that the morphology is due to a chance alignment. The Stingrays exhibit non-thermal emission with spectral indices of $\alpha=-0.89\pm0.09$ for Stingray 1 and $\alpha=-1.77\pm0.06$ for Stingray 2. We find that none of the proposed scenarios can explain all of the observed properties, however we determine it most likely that their shape is caused by some kind of complex environmental interaction. The most likely scenario from the available data is that of a head-tail radio galaxy, but more data is required for a definitive classification.

We present the serendipitous radio-continuum discovery of a likely Galactic supernova remnant (SNR) G305.4–2.2. This object displays a remarkable circular symmetry in shape, making it one of the most circular Galactic SNRs known. Nicknamed Teleios due to its symmetry, it was detected in the new Australian Square Kilometre Array Pathfinder (ASKAP) Evolutionary Map of the Universe (EMU) radio–continuum images with an angular size of 1 320$^{\prime\prime}$$\times$1 260$^{\prime\prime}$ and PA = 0$^\circ$. While there is a hint of possible H$\alpha$ and gamma-ray emission, Teleios is exclusively seen at radio–continuum frequencies. Interestingly, Teleios is not only almost perfectly symmetric, but it also has one of the lowest surface brightnesses discovered among Galactic SNRs and a steep spectral index of $\alpha$=–0.6$\pm$0.3. Our best estimates from Hi studies and the $\Sigma$–D relation place Teleios as a type Ia SNR at a distance of either $\sim$2.2 kpc (near-side) or $\sim$7.7 kpc (far-side). This indicates two possible scenarios, either a young (under 1 000 yr) or a somewhat older SNR (over 10 000 yr). With a corresponding diameter of 14/48 pc, our evolutionary studies place Teleios at the either early or late Sedov phase, depending on the distance/diameter estimate. However, our modelling also predicts X-ray emission, which we do not see in the present generation of eROSITA images. We also explored a type Iax explosion scenario that would point to a much closer distance of $\lt$1 kpc and Teleios size of only $\sim$3.3 pc, which would be similar to the only known type Iax remnant SN1181. Unfortunately, all examined scenarios have their challenges, and no definitive Supernova (SN) origin type can be established at this stage. Remarkably, Teleios has retained its symmetrical shape as it aged even to such a diameter, suggesting expansion into a rarefied and isotropic ambient medium. The low radio surface brightness and the lack of pronounced polarisation can be explained by a high level of ambient rotation measure (RM), with the largest RM being observed at Teleios’s centre.

The Australian SKA Pathfinder (ASKAP) offers powerful new capabilities for studying the polarised and magnetised Universe at radio wavelengths. In this paper, we introduce the Polarisation Sky Survey of the Universe’s Magnetism (POSSUM), a groundbreaking survey with three primary objectives: (1) to create a comprehensive Faraday rotation measure (RM) grid of up to one million compact extragalactic sources across the southern $\sim50$% of the sky (20,630 deg$^2$); (2) to map the intrinsic polarisation and RM properties of a wide range of discrete extragalactic and Galactic objects over the same area; and (3) to contribute interferometric data with excellent surface brightness sensitivity, which can be combined with single-dish data to study the diffuse Galactic interstellar medium. Observations for the full POSSUM survey commenced in May 2023 and are expected to conclude by mid-2028. POSSUM will achieve an RM grid density of around 30–50 RMs per square degree with a median measurement uncertainty of $\sim$1 rad m$^{-2}$. The survey operates primarily over a frequency range of 800–1088 MHz, with an angular resolution of 20” and a typical RMS sensitivity in Stokes Q or U of 18 $\mu$Jy beam$^{-1}$. Additionally, the survey will be supplemented by similar observations covering 1296–1440 MHz over 38% of the sky. POSSUM will enable the discovery and detailed investigation of magnetised phenomena in a wide range of cosmic environments, including the intergalactic medium and cosmic web, galaxy clusters and groups, active galactic nuclei and radio galaxies, the Magellanic System and other nearby galaxies, galaxy halos and the circumgalactic medium, and the magnetic structure of the Milky Way across a very wide range of scales, as well as the interplay between these components. This paper reviews the current science case developed by the POSSUM Collaboration and provides an overview of POSSUM’s observations, data processing, outputs, and its complementarity with other radio and multi-wavelength surveys, including future work with the SKA.

We present a new radio continuum study of the Large Magellanic Cloud supernova remnant (SNR) MC SNR J0519–6902. With a diameter of $\sim$8 pc, this SNR shows a radio ring-like morphology with three bright regions towards the north, east, and south. Its linear polarisation is prominent with average values of $5\pm1$% and $6\pm1$% at 5 500 and 9 000 MHz, and we find a spectral index of ${-0.62\pm0.02}$, typical of a young SNR. The average rotation measure is estimated at ${-124\pm83}$ rad m$^{-2}$ and the magnetic field strength at ${\sim11}\,\mu$G. We also estimate an equipartition magnetic field of ${72\pm 5}\,\mu$G and minimum explosion energy of E$_\textrm{ min}$ = 2.6$\times10^{48}$ erg. Finally, we identified an H i cloud that may be associated with MC SNR J0519–6902, located in the southeastern part of the remnant, along with a potential wind-bubble cavity.

We present the Evolutionary Map of the Universe (EMU) survey conducted with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU aims to deliver the touchstone radio atlas of the southern hemisphere. We introduce EMU and review its science drivers and key science goals, updated and tailored to the current ASKAP five-year survey plan. The development of the survey strategy and planned sky coverage is presented, along with the operational aspects of the survey and associated data analysis, together with a selection of diagnostics demonstrating the imaging quality and data characteristics. We give a general description of the value-added data pipeline and data products before concluding with a discussion of links to other surveys and projects and an outline of EMU’s legacy value.

We present a new radio detection from the Australian Square Kilometre Array Pathfinder Evolutionary Map of the Universe (EMU) survey associated with the Reflection Nebula (RN) VdB-80. The radio detection is determined to be a previously unidentified Hii region, now named Lagotis. The RN is located towards Monoceros, centred in the molecular cloud feature known as the ‘Crossbones’. The 944 MHz EMU image shows a roughly semicircular Hii region with an integrated flux density of 30.2$\pm$0.3 mJy. The Hii region is also seen at 1.4 GHz by NRAO VLA Sky Survey (NVSS), yielding an estimated spectral index of 0.65$\pm$0.51, consistent with thermal radio emission. Gaia Data Release 3 (DR3) and Two Micron All Sky Survey (2MASS) data give a distance to the stars associated with the Hii region of $\sim$960 pc. This implies a size of 0.76$\times$0.68($\pm$0.09) pc for the Hii region. We derive an Hii region electron density of the bright radio feature to be 26 cm$^{-3}$, requiring a Lyman-alpha photon flux of $10^{45.6}$ s$^{-1}$, which is consistent with the expected Lyman flux of HD 46060, the B2 ii type star which is the likely ionising star of the region. The derived distance to this region implies that the Crossbones feature is a superposition of two filamentary clouds, with Lagotis embedded in the far cloud.

We present a re-discovery of G278.94+1.35a as possibly one of the largest known Galactic supernova remnants (SNRs) – that we name Diprotodon. While previously established as a Galactic SNR, Diprotodon is visible in our new Evolutionary Map of the Universe (EMU) and GaLactic and Extragalactic All-sky MWA (GLEAM) radio continuum images at an angular size of $3{{{{.\!^\circ}}}}33\times3{{{{.\!^\circ}}}}23$, much larger than previously measured. At the previously suggested distance of 2.7 kpc, this implies a diameter of 157$\times$152 pc. This size would qualify Diprotodon as the largest known SNR and pushes our estimates of SNR sizes to the upper limits. We investigate the environment in which the SNR is located and examine various scenarios that might explain such a large and relatively bright SNR appearance. We find that Diprotodon is most likely at a much closer distance of $\sim$1 kpc, implying its diameter is 58$\times$56 pc and it is in the radiative evolutionary phase. We also present a new Fermi-LAT data analysis that confirms the angular extent of the SNR in gamma rays. The origin of the high-energy emission remains somewhat puzzling, and the scenarios we explore reveal new puzzles, given this unexpected and unique observation of a seemingly evolved SNR having a hard GeV spectrum with no breaks. We explore both leptonic and hadronic scenarios, as well as the possibility that the high-energy emission arises from the leftover particle population of a historic pulsar wind nebula.

We report the discovery of a bow-shock pulsar wind nebula (PWN), named Potoroo, and the detection of a young pulsar J1638$-$4713 that powers the nebula. We present a radio continuum study of the PWN based on 20-cm observations obtained from the Australian Square Kilometre Array Pathfinder (ASKAP) and MeerKAT. PSR J1638$-$4713 was identified using Parkes radio telescope observations at frequencies above 3 GHz. The pulsar has the second-highest dispersion measure of all known radio pulsars (1 553 pc cm$^{-3}$), a spin period of 65.74 ms and a spin-down luminosity of $\dot{E}=6.1\times10^{36}$ erg s$^{-1}$. The PWN has a cometary morphology and one of the greatest projected lengths among all the observed pulsar radio tails, measuring over 21 pc for an assumed distance of 10 kpc. The remarkably long tail and atypically steep radio spectral index are attributed to the interplay of a supernova reverse shock and the PWN. The originating supernova remnant is not known so far. We estimated the pulsar kick velocity to be in the range of 1 000–2 000 km s$^{-1}$ for ages between 23 and 10 kyr. The X-ray counterpart found in Chandra data, CXOU J163802.6$-$471358, shows the same tail morphology as the radio source but is shorter by a factor of 10. The peak of the X-ray emission is offset from the peak of the radio total intensity (Stokes $\rm I$) emission by approximately 4.7$^{\prime\prime}$, but coincides well with circularly polarised (Stokes $\rm V$) emission. No infrared counterpart was found.

We demonstrate the importance of radio selection in probing heavily obscured galaxy populations. We combine Evolutionary Map of the Universe (EMU) Early Science data in the Galaxy and Mass Assembly (GAMA) G23 field with the GAMA data, providing optical photometry and spectral line measurements, together with Wide-field Infrared Survey Explorer (WISE) infrared (IR) photometry, providing IR luminosities and colours. We investigate the degree of obscuration in star-forming galaxies, based on the Balmer decrement (BD), and explore how this trend varies, over a redshift range of $0<z<0.345$. We demonstrate that the radio-detected population has on average higher levels of obscuration than the parent optical sample, arising through missing the lowest BD and lowest mass galaxies, which are also the lower star formation rate (SFR) and metallicity systems. We discuss possible explanations for this result, including speculation around whether it might arise from steeper stellar initial mass functions in low mass, low SFR galaxies.

We present a comparison between the performance of a selection of source finders (SFs) using a new software tool called Hydra. The companion paper, Paper I, introduced the Hydra tool and demonstrated its performance using simulated data. Here we apply Hydra to assess the performance of different source finders by analysing real observational data taken from the Evolutionary Map of the Universe (EMU) Pilot Survey. EMU is a wide-field radio continuum survey whose primary goal is to make a deep ($20\mu$Jy/beam RMS noise), intermediate angular resolution ($15^{\prime\prime}$), 1 GHz survey of the entire sky south of $+30^{\circ}$ declination, and expecting to detect and catalogue up to 40 million sources. With the main EMU survey it is highly desirable to understand the performance of radio image SF software and to identify an approach that optimises source detection capabilities. Hydra has been developed to refine this process, as well as to deliver a range of metrics and source finding data products from multiple SFs. We present the performance of the five SFs tested here in terms of their completeness and reliability statistics, their flux density and source size measurements, and an exploration of case studies to highlight finder-specific limitations.

We investigate the diversity in the sizes and average surface densities of the neutral atomic hydrogen (H i) gas discs in $\sim$280 nearby galaxies detected by the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY). We combine the uniformly observed, interferometric H i data from pilot observations of the Hydra cluster and NGC 4636 group fields with photometry measured from ultraviolet, optical, and near-infrared imaging surveys to investigate the interplay between stellar structure, star formation, and H i structural parameters. We quantify the H i structure by the size of the H i relative to the optical disc and the average H i surface density measured using effective and isodensity radii. For galaxies resolved by $>$$1.3$ beams, we find that galaxies with higher stellar masses and stellar surface densities tend to have less extended H i discs and lower H i surface densities: the isodensity H i structural parameters show a weak negative dependence on stellar mass and stellar mass surface density. These trends strengthen when we limit our sample to galaxies resolved by $>$2 beams. We find that galaxies with higher H i surface densities and more extended H i discs tend to be more star forming: the isodensity H i structural parameters have stronger correlations with star formation. Normalising the H i disc size by the optical effective radius (instead of the isophotal radius) produces positive correlations with stellar masses and stellar surface densities and removes the correlations with star formation. This is due to the effective and isodensity H i radii increasing with mass at similar rates while, in the optical, the effective radius increases slower than the isophotal radius. Our results are in qualitative agreement with previous studies and demonstrate that with WALLABY we can begin to bridge the gap between small galaxy samples with high spatial resolution H i data and large, statistical studies using spatially unresolved, single-dish data.

The latest generation of radio surveys are now producing sky survey images containing many millions of radio sources. In this context it is highly desirable to understand the performance of radio image source finder (SF) software and to identify an approach that optimises source detection capabilities. We have created Hydra to be an extensible multi-SF and cataloguing tool that can be used to compare and evaluate different SFs. Hydra, which currently includes the SFs Aegean, Caesar, ProFound, PyBDSF, and Selavy, provides for the addition of new SFs through containerisation and configuration files. The SF input RMS noise and island parameters are optimised to a 90% ‘percentage real detections’ threshold (calculated from the difference between detections in the real and inverted images), to enable comparison between SFs. Hydra provides completeness and reliability diagnostics through observed-deep ($\mathcal{D}$) and generated-shallow ($\mathcal{S}$) images, as well as other statistics. In addition, it has a visual inspection tool for comparing residual images through various selection filters, such as S/N bins in completeness or reliability. The tool allows the user to easily compare and evaluate different SFs in order to choose their desired SF, or a combination thereof. This paper is part one of a two part series. In this paper we introduce the Hydra software suite and validate its $\mathcal{D/S}$ metrics using simulated data. The companion paper demonstrates the utility of Hydra by comparing the performance of SFs using both simulated and real images.

Referrals to Child and Adolescent Mental Health Services (CAMHS) have increased in recent years. Services are already under-resourced and the adverse psychological impact of Covid-19 is likely to increase demand. Accordingly, an understanding of prevalence of mental health (MH) disorders among youth is imperative to help inform and plan services.

We present the first unbiased survey of neutral hydrogen absorption in the Small Magellanic Cloud. The survey utilises pilot neutral hydrogen observations with the Australian Square Kilometre Array Pathfinder telescope as part of the Galactic Australian Square Kilometre Array Pathfinder neutral hydrogen project whose dataset has been processed with the Galactic Australian Square Kilometre Array Pathfinder-HI absorption pipeline, also described here. This dataset provides absorption spectra towards 229 continuum sources, a 275% increase in the number of continuum sources previously published in the Small Magellanic Cloud region, as well as an improvement in the quality of absorption spectra over previous surveys of the Small Magellanic Cloud. Our unbiased view, combined with the closely matched beam size between emission and absorption, reveals a lower cold gas faction (11%) than the 2019 ATCA survey of the Small Magellanic Cloud and is more representative of the Small Magellanic Cloud as a whole. We also find that the optical depth varies greatly between the Small Magellanic Cloud’s bar and wing regions. In the bar we find that the optical depth is generally low (correction factor to the optically thin column density assumption of $\mathcal{R}_{\mathrm{HI}} \sim 1.04$) but increases linearly with column density. In the wing however, there is a wide scatter in optical depth despite a tighter range of column densities.

We present the most sensitive and detailed view of the neutral hydrogen (${\rm H\small I}$) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal ${\rm H\small I}$ in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K ($1.6\,\mathrm{mJy\ beam}^{-1}$) $\mathrm{per}\ 0.98\,\mathrm{km\ s}^{-1}$ spectral channel with an angular resolution of $30^{\prime\prime}$ (${\sim}10\,\mathrm{pc}$). We discuss the calibration scheme and the custom imaging pipeline that utilises a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire ${\sim}25\,\mathrm{deg}^2$ field-of-view. We provide an overview of the data products and characterise several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high-velocity cloud with previous ASKAP+Parkes ${\rm H\small I}$ test observations.