The Rapid ASKAP Continuum Survey V: Cataloguing the sky at 1 367.5 MHz and the second data release of RACS-mid

Abstract 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/).


Introduction
The Australian SKA Pathfinder (ASKAP; Hotan et al., 2021), located at Inyarrimanha Ilgari Bundara, the CSIRO a Murchison Radio-astronomy Observatory in Western Australia, is a survey radio telescope that operates with a 288-MHz bandwidth between 700 and 1800 MHz.ASKAP features a unique Phased Array Feed (PAF; Hotan et al., 2014;McConnell et al., 2016) which allows the formation of 36 primary beams covering up to ∼ 31 deg 2 .This instantaneous wide field of view coupled with 36 × 12 m antennas-with baselines ranging from 22 m to 6 km-provides the basis for an instrument optimised for survey speed, with the capability to rapidly survey ∼ 80% of the sky to heretofore unseen sensitivity.
The Rapid ASKAP Continuum Survey (RACS) is a multiwavelength survey covering the sky up to δ J2000 ≲ +50 • (Mc-Connell et al., 2020, hereinafter, Paper I).The survey is being conducted within the three operating bands used by ASKAP.For RACS, these three bands are centered on effective frequena Commonwealth Scientific and Industrial Research Organisation.cies of 887.5, 1367.5, and 1632.5 MHz.As the name suggests, RACS is a shallow survey with 15-min total integration per pointing reaching median root-mean-square (rms) noise between ∼ 200-260 µJy PSF -1 across the three bands.Due to its baseline distribution, ASKAP has excellent snapshot (u, v) coverage, providing a well-constrained point-spread function at zenith (δ J2000 ≈ -27 • ) with elongation of the major axis for lower-elevation pointings.The resultant median angular resolution of the whole survey across the three bands is in the range ∼ 9-18 arcsec.The goal of RACS is to provide a global sky model for calibration and validation of other ASKAP observations, and once images and source catalogues are available work can be done matching and spectrally modelling sources across the three bands.
The first pass of the sky in the lowest frequency band at 887.5 MHz (RACS-low) was released in 2020 and is described in Paper I. RACS-low features a median rms noise of ∼ 260 µJy PSF -1 and median angular resolution of 18 arcsec.Following the initial data release, a radio source catalogue was constructed and described by Hale et al. (2021, hereinafter Paper II).This radio source catalogue contains ∼ 2.1 million sources, excluding Galactic latitudes of |b| < 5 • .The catalogue was constructed after combining neighbouring images from the original data release to provide a more uniform sensitivity across the region covering -80 • ≤ δ J2000 ≤ +30 • .This catalogue is the first large-area catalogue produced as part of RACS.Alongside the Stokes I total intensity data releases, RACS-low data products are also being used to investigate the polarized sky.Using circular polarization data products from RACS-low, Pritchard et al. (2021) identified radio emission coincident with 33 stars, including 23 new radio star associations in this sample.In linear polarization, RACS data products are being used to generate Spectra and Polarisation In Cutouts of Extragalactic Sources (SPICE-RACS; Thomson et al., 2023, Paper III).The first data release from SPICE-RACS covers a 1 300 deg 2 pilot region and features 5 818 sources with Faraday rotation measures.
The second observing epoch of the sky at 1367.5 MHz (RACS-mid) was completed in 2021, with follow-up observations taken over 2021-2022, and the first data release is described by Duchesne et al. (2023, hereinafter Paper IV). RACS-mid data products, along with all available RACS data, are publicly available through the CSIRO ASKAP Science Data Archive (CASDA; Chapman et al., 2017;Huynh et al., 2020).RACS-mid features a number of changes from RACS-low that helped to improve the overall quality of the survey.These changes include: fully-autonomous scheduling with a limit to the hour angle of ±1 h , peeling and subtraction of bright offaxis sources, bespoke primary beam modelling, and correction of both on-axis and off-axis leakage of Stokes I into V.While RACS-mid has an effective bandwidth of 144-MHz due to RFI, the survey still boasts a median rms noise of ∼ 200 µJy PSF -1 , which varies smoothly across the sky as a function of declination, except around bright sources.Similarly, because of the close-to-meridian scheduling the angular resolution varies smoothly as a function of declination ranging from 8.1 arcsec near zenith and 47.5 arcsec in the lowest-elevation pointings at high declination.The first RACS-mid data release also included leakage-corrected Stokes V continuum images alongside the Stokes I images for each individual observation.
RACS-mid is both competitive with and complementary to existing completed large-area 1.4-GHz radio surveys in the Northern sky such as the NRAO b VLA c Sky Survey (NVSS; Condon et al., 1998) and the Faint Images of the Radio Sky at Twenty Centimeters (FIRST; Becker et al., 1995;White et al., 1997;Helfand et al., 2015).In the highest declination strip covered by RACS-mid (δ J2000 = +46 • ), the resolution and sensitivity approaches that of the NVSS.For the remainder of the survey the resolution and sensitivity can be expressed as an average of the survey properties of the NVSS and FIRST, with a median angular resolution of 10.1 arcsec and moderate sensitivity to extended sources up to a few arcmin.RACS-mid data products have so far been used to identify circularly polarized emission from 52 stars that had not previously been detected at b National Radio Astronomy Observatory.c Very Large Array.radio frequencies (Driessen et al. in prep), including the detection of circularly polarized emission from a T8 brown dwarf (Rose et al., 2023), and have been used in conjunction with FIRST and the VLA Sky Survey (VLASS; Gordon et al., 2021) in a search for radio stars via proper motion (Driessen et al., 2023).RACS-low and RACS-mid images are also used as an extra epoch for the Variability and Slow Transients (VAST) survey being performed by ASKAP (Murphy et al., 2013(Murphy et al., , 2021)).
At 1367.5 MHz, ASKAP is also undertaking the Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY), which is a neutral hydrogen survey being performed over the next five years with using the same frequency band as RACSmid (Koribalski et al., 2020).While WALLABY is primarily a spectral line survey, the 10-hour observations and subsequent imaging also result in deep 1367.5-MHzcontinuum images.Currently pre-pilot observations of the Eridanus supergroup (For et al., 2021) have been utilised to create an initial WAL-LABY continuum catalogue (hereinafter the WALLABY prepilot catalogue; Grundy et al., 2023) containing 9 416 sources covering ∼ 42 deg 2 .
In this fifth paper in the RACS series, we describe the cataloguing of the RACS-mid Stokes I data.Although Stokes V continuum images are available, cataloguing circularly polarized sources with RACS-mid will be described in a future work.

Catalogue creation
In cataloguing the RACS-mid Stokes I data, we follow the steps outlined for the construction of the RACS-low catalogue described by Paper II.We make some changes to this process to better suit the current data products.As a departure from Paper II, we opted to create a primary catalogue for general purpose use, as well as two auxiliary catalogues for more specific uses.The three catalogues cover the normal uses for Stokes I source catalogues, each with their own limitations and strengths.In the following sections, we describe the general steps taken to construct the catalogues, including selection of which observations to include (Section 2.1), making fullsensitivity images (Section 2.2), source-finding (Section 2.3), and merging source-lists and duplicate removal (Section 2.4).The primary catalogue and the two auxiliary catalogues are described in Section 3 and are created following the relevant subset of the following steps.

Image selection
As part of the RACS-mid observing process, some fields were observed and imaged several times.Fields have the naming convention RACS_HHMM±DD, and most duplicate field observations and their resultant Stokes I and V images are made available through CASDA.As described in section 2.1 in Paper IV, some of these re-observations were intended to provide a higher-resolution or to account for observations of < 12 min (see section 2.1 in Paper IV). Fields that were re-observed due to fatal errors with the instrument are not considered.We define a 'good' subset of observed fields to use in the cataloguing    process.Each RACS-mid observation has a dedicated scheduling block ID (SBID) which is used as a unique identifier.The selection of which SBID to use when multiple observations exist for a given field is based on the following criteria: 1 If either duplicate SBID had bright sources peeled (as described in section 2.3.3 of Paper IV), then use the SBID that has undergone peeling, 2 Else if either duplicate SBIDs have < 12 m total integration time, then Use the SBID with the longest integration time, 3 Else if the difference in PSF major axes is > 2 ′′ between duplicate SBIDs, then select the SBID with the smaller PSF major axes, 4 Else select the SBID with the lowest median rms noise.
The ordering is important, as peeling in a small number of cases was required to remove CLEAN divergence artefacts but was not necessarily done on the highest-resolution image.The SBIDs with < 12 m total integration time are significantly less sensitive than their re-observations so naturally should not be used when more sensitive observations exists.No short observations had better PSF or noise properties than their respective re-observations, and none of these short observations underwent peeling.Finally, resolution and noise properties are often closely related due to their direct relationship to (u, v) coverage and we find that generally the selection is based on the PSF size rather than the median rms noise if criteria (1) and ( 2) are not met.A single field was observed three times (RACS_1844+46), though only in the second observation was 3C 380 peeled and is chosen as the good SBID.The selection of good SBIDs is noted in the RACS database file d under the column SELECT, where SELECT=1 indicates the SBID is selected as part of the good SBID subset.We note that SBIDs not part of this subset can still be used for other science cases and may be particularly useful for variability or transient source studies.The selection is based on the Stokes I images only, and future work with Stokes V may wish to select different SBIDs if noise properties vary between them.

Full-sensitivity images
A 12% cut-off in the primary beam attenuation was used for primary beam correction and weighting the PAF beam mosaics for each observation.One of the unfortunate features of the linear mosaicking of individual PAF beam images is a tile mosaic with non-rectangular image boundaries and the usual primary beam roll-off.This results in an irregular tile image border that follows the 12% contour of the outer beams.From a science standpoint, the primary beam roll-off limits detection of sources towards image boundaries due to a drop in sensitivity.From a user-perspective, rectangular images with large, irregular areas of blanked pixels can make identifying which sources lie in which image difficult.
As RACS is observed with an overlap between adjacent observations that is at least the PAF beam separation (0.9 degrees for RACS-mid), we follow Paper II and use SWarp (Bertin   d https://bitbucket.csiro.au/projects/ASKAP_SURVEYS/repos/racs/browse.et al., 2002) to form a linear mosaic of adjacent tiles to create rectangular images with full sensitivity towards the image edges.We use the same weight maps used during PAF beam linear mosaicking to ensure the sensitivity profile is consistent across the full image.All adjacent tiles within a group are convolved to a common resolution e prior to mosaicking.We create one set of images convolved to the lowest common resolution of neighbouring images, and a second set convolved to a fixed 25-arcsec resolution.These images are provided as part of this data release.Convolving to the lowest common resolution results in mosaics that are close to the full resolution of the survey.These full-sensitivity images are created for Stokes I and V, though only Stokes I images are part of this data release.For the 25-arcsec images, we also increase the pixel size by a factor two in each celestial coordinate (i.e. from 2 ′′ × 2 ′′ to 4 ′′ × 4 ′′ ) to reduce the image file size.The 25-arcsec full-sensitivity images (and catalogue) are restricted to δ J2000 ≲ +30 • , similar to RACS-low, as the PSF major axis for some of the original images becomes larger than 25 arcsec beyond this declination.
Figure 1 shows example sources as they appear in the full-sensitivity images compared to the 25-arcsec image from RACS-low (Paper II), and comparison optical data from from the Dark Energy Survey data release 2 (DES DR2; Flaugher et al., 2015;Abbott et al., 2018Abbott et al., , 2021)).Figure 1(i) shows an example of radio sources projected onto the galaxy cluster Abell 2829, highlighting a range of angular scales and features.Figures 1(ii) and 1(iii) shows two radio galaxies, PKS 2357+00 and PKS 0229-208, and 1(iv) shows the spiral galaxy NGC 7582.Note these examples highlight situations where the angular resolution of the full-sensitivity image (9.8 ′′ ×8.5 ′′ and 11.4 ′′ ×9.2 ′′ ) is significantly reduced to form the 25-arcsec images.
Similarly, Figure 2 shows example sources in the highestdeclination strip of the survey (δ J2000 = +46 • ) with comparison to the NVSS, FIRST, and either the Digitized Sky Survey (DSS2) or the Sloan Digitized Sky Survey data release 7 (SDSS DR7; Abazajian et al., 2009).Figures 2(i)-2(iii) show the spiral galaxies M 51,M 106,and NGC 4096,and Figure 2(iv) shows the giant radio galaxy B3 1206+469 (Dabhade et al., 2020).Due to the declination, the 25-arcsec images from RACS-mid and RACS-low are not available, and the angular resolution of the RACS-mid data approaches that of the NVSS.

Source finding and measuring
Source-finding is performed on both the full-sensitivity images as well as the original images using PyBDSF (Mohan & Rafferty, 2015).To help limit the number of artefacts recovered during source-finding, we modify the box and grid size for the position-dependent rms calculations to account for the varying PSF (as a function of the minor axis, θ minor ).We find (15θ minor , 3θ minor ) provides a good balance in removing artefacts around bright sources though removes a similar number of faint (∼ 5σ rms ) sources from the resulting source-lists.This box and grid size is functionally the same (150, 30) pixel box e Using beamcon_2D from https://github.com/AlecThomson/RACS-tools,which is a generic radio image 2-D convolution tool.used by Paper II for their resolution and pixel size.The choice to scale with the minor axis of the PSF is to account for the smallest angular scale in the image.
For each of the individual source-lists, we identify sources with ab θ M θ m < median ab θ M θ m (1) for sources with major and minor axes a and b in images with PSF major and minor axes θ M and θ m that are comprised of more than one Gaussian component.Such sources are either point sources (or close to) and should not be decomposed into multiple components.Such decomposition is an artefact created during PyBDSF fitting where noise fluctuations lead to irregular boundaries of pixel islands.We find this to affect sources at all signal-to-noise ratios (SNRs).We group these components into a single source with the sum of integrated flux densities and a 2-D elliptical Gaussian to represent the size and shape.
For each source list the PSF of the image (major and minor axes FWHM and the position angle) is recorded for each source.For the individual image source lists these are single-valued over the individual source lists, but are retained when merging the individual source lists to ensure that each catalogued source has the correct PSF associated with it from the image within which it is found.

Merging source-lists and duplicate removal
As there is significant overlap between the full-sensitivity images (by construction of the survey) these overlapping regions contain similar information.For the primary catalogue, the overlap regions may differ as the PSF for each full-sensitivity image may be different.This difference is small, but can result in two effects: (1) a difference to detected sources due to a change in SNR, and (2) a difference in components fit to detected sources.For (1), this results in overlap regions containing different sources at the low-SNR end of the source counts.For (2), this becomes particularly noticeable for, e.g., double and triple radio sources where the components in one image are connected and catalogued as a single source in one image but decompose into multiple sources in an adjacent image.Additionally, a combination of (1) and (2) can occur where components from multi-component sources get pushed below SNR limits.
We construct the catalogues by merging individual image source-lists one-by-one, and cross-match sources in the incoming source-list to the partially constructed catalogue.Sources in the partially constructed catalogue that match to a source in the incoming source lists with an angular separation, s, satisfying are considered duplicates, where θ M,1 is the PSF major axis of the source in the partially constructed catalogue and θ M,2 is the PSF major axis of the source in the incoming source list.When a duplicate is found, the source that has the smallest separation from an image centre is kept.While constructing the catalogue, we also avoid sources within 2 arcmin of an image boundary to avoid cross-matching of extended sources that might be cut off at the edge of an image.The Gaussian components associated with the final set of sources are then collected into a separate component list.

The catalogues
The catalogues produced as part of this data release are summarised in Table 1.These catalogues include (1) the primary catalogue comprising our best characterisation of each source at the observed and position-dependent resolution; (2) an auxiliary fixed 25-arcsec resolution catalogue, prepared similarly to the primary catalogue but with images convolved to a resolution of 25-arcsec; (3) an auxiliary time-domain catalogue, which includes all sources including those detected multiple times in adjacent images or duplicate observations of a particular field.The catalogues have identical columns, which are described in Appendix A. Note that in Appendix A we also detail some changes to the columns and names used between the RACSmid catalogues and the RACS-low catalogue.Tables A3 and  A4 in Appendix A also show example rows from the primary source and component catalogues, respectively.The catalogues are detailed further in the following sections.a Uncertainties are reported from the 16 th and 84 th percentiles.b In parenthesis excluding the Galactic Plane (b ± 5 • ).

The primary catalogue
The first catalogue is considered the primary catalogue and will be suitable for most general users.This catalogue covers the full region observed as part of RACS-mid and remains close to the angular resolution of the original RACS-mid images described in Paper IV. Construction of this catalogue begins with identifying adjacent fields, convolving to the lowest common resolution of that subset of fields, then forming full-sensitivity images as described in Section 2.2.Examples of the image quality is shown in Figures 1 and 2 for a range of declinations, and Figure 3 shows the variation of the PSF major [3(i)] and minor [3(ii)] axes as a function of declination for sources in the catalogue.The PSF for the primary catalogue is elliptical, and the major axis reaches 48.4 arcsec (at high declination) while the minor axis has a maximum of 18.9 arcsec (at low declination).

Auxiliary fixed-resolution, 25-arcsec catalogue
The second catalogue is an auxiliary 'fixed-resolution' catalogue, with a position-independent resolution of 25×25 arcsec 2 .
Figure 3 shows the 25-arcsec PSF compared to other catalogues' PSF variation for reference.This catalogue is designed to match the RACS-low catalogue described by Paper II, including the same resolution and sky coverage up to δ J2000 ≤ +30 • f .

Auxiliary time-domain catalogue
As the previously mentioned catalogues are created after linearly mosaicking neighbouring observations, variable/transient radio sources may have their emission averaged between what may be two disjoint epochs.For many time-domain applications, it may be preferable to retain source detections and measurements at specific epochs.Searches for and characterisation of variable and transient sources have already been performed using individual images and sources lists from RACS-low (e.g.Leung et al., 2021) and RACS-mid (e.g.Driessen et al., 2023;Gulati et al., 2023).To help enable further time-domain science, we provide another auxiliary catalogue for Stokes I that is simply the concatenation of source-lists from the individual RACS-mid images.The individual images are not convolved f Note that the RACS-low catalogue described in Paper II features a hole in its coverage below δ J2000 < -80 • as some images did not meet the 25arcsec resolution requirement.For RACS-mid we provide full coverage below δ J2000 ≤ +30 • for the 25-arcsec catalogue as all images below this declination could be convolved to 25 ′′ × 25 ′′ .
prior to source-finding so they retain the original angular resolution.Figure 3 shows the variation of the PSF at the locations of sources in the catalogue compared to the primary and 25-arcsec catalogues.As with the primary catalogue the PSF is elliptical and has a similar range of values.Sources that are detected in multiple images in overlap regions (i.e.duplicates) are not removed.

Analysis of the catalogues and images
In the following sections, we provide an analysis of the overall quality of the available catalogues and full-sensitivity image mosaics made.As the catalogue are created from the same original images, there is some redundancy to the following analysis and validation work.For consistency we show results for all catalogues where appropriate.

Overall noise properties
The position-dependent rms noise properties of the images that are used for the three catalogues varies due to the difference in angular resolution and their overall construction.A local rms noise estimate is included for each source in the catalogue.Figure 4 shows the rms noise median-binned using Hierarchical Equal Area isoLatitude Pixelation (HEALPix; Górski et al., 2005) g with N side = 64 corresponding to ∼ 55 × 55 arcmin 2 bins.The rms values are reported at the location of the source, so will on average be marginally elevated compared to offsource regions.
The rms noise distribution of the primary catalogue [Figure 4(i)] follows closely the median Stokes I noise per tile shown in figure 24 of Paper IV, which was constructed by mosaicking position-dependent rms maps.The primary catalogue has a median rms noise of σ rms = 182 +41 -23 µJy PSF -1 (with uncertainties derived from the 16 th and 84 th percentiles of the rms noise distribution).This median rms noise is a decrease from the median rms noise reported over the original survey images reported in Paper IV due to the removal of the primary beam roll-off at the edge of convolved images.The larger spread in the distribution is due to the bias in local rms estimates taken from source locations, including the Galactic Plane and other areas of bright, extended sources.While not part of this data release, we note that the Stokes V images mosaicked and convolved in the same way achieve a median rms  noise of σ rms,V = 144 +20 -14 µJy PSF -1 at the catalogued source positions.
The 25-arcsec catalogue [Figure 4(ii)] has a median rms noise σ rms = 278 +68 -47 µJy PSF -1 , generally showing an increase in rms noise over most of the sky due to the decrease in angular resolution.Of particular note is an additional increase in noise near the celestial equator, corresponding to residuals from the poorer equatorial PSF becoming significant after convolution to the lower 25-arcsec resolution.Convolving the images to the lower 25 × 25 arcsec 2 PSF also results in a general loss of information in the images, translating to an overall increase in the image rms noise.The difference in the median Stokes I and Stokes V noise (σ rms,V = 167 +23  -17 µJy PSF -1 ) is comparatively larger than for the primary catalogue, suggesting another contribution to the increase in noise is the increase in classical source confusion at this lower resolution (Condon, 1974) though the 25-arcsec catalogue is not confusion-limited (see Section 4.2).
Finally, the time-domain catalogue [Figure 4(iii)] features increased noise at the boundaries between fields due to overlapping uncorrected primary beam roll-off.The remainder of the survey area has similar position-dependent noise characteristics to the primary catalogue due to the similarities towards the image centres.The median rms noise, σ rms = 203 +140 -33 µJy PSF -1 is closer to the tile median reported in Paper IV with an increase in spread due to bias in the sampled noise locations (i.e. at the locations of sources).

Classical source confusion
While we do not expect even the deeper 10-hr images from ASKAP to become limited by classical source confusion noise at 1.4 GHz (assuming a 10 × 10 arcsec 2 PSF; Condon et al., 2012), as this noise is a function of PSF solid angle [Ω PSF = πθ M θ m / (4 ln 2)] it is worth considering a possible limit given the variable PSFs of the RACS-mid catalogues.The number of PSF solid angles per source above some flux density threshold is often used as an estimate of the severity of confusion (e.g.    , 3, 4, 5, 10] mJy) for the primary catalogue (i), the 25-arcsec catalogue (ii), and the time-domain catalogue (iii).Dashed lines correspond to median source densities for the associated flux density limit.The grey, shaded region in (ii) is not covered in the 25-arcsec catalogue.Cohen, 2004;Condon et al., 2012;Heywood et al., 2016), with values of order ≲ 10 indicating that the images are limited in sensitivity by confusion.
Figure 5 shows the number of PSF solid angles per source above 5σ rms (N Ω /N sources ) binned as a function of declination, excluding the area covering Galactic latitudes of |b| < 10 • .We also show the equivalent product with the same threshold for some other surveys with largely contiguous sky coverage: the NVSS (Condon et al., 1998) with a 45 × 45 arcsec 2 PSF, the Sydney University Molonglo Sky Survey (SUMSS; Bock et al., 1999;Mauch et al., 2003) with a position-dependent 45 × 45/sin|δ J2000 | arcsec 2 PSF, and the Galactic and Extragalactic All-sky MWA h survey (GLEAM; Wayth et al., 2015) extragalactic catalogue (Hurley-Walker et al., 2017) with a position-dependent PSF major axis that varies between ∼ 128-240 arcsec and minor axis that varies between ∼ 122-173 arcsec at 200 MHz.Note we do not include the RACS-mid timedomain catalogue as the overlapping observations artificially increase the source density and are not reflective of the actual images used.The shaded regions in Figure 5 show the range of values within a declination bin for surveys with variable PSFs.
Despite the elongated PSF of the primary catalogue at high declination, the solid angle remains small due to the high ellipticity of the PSF (θ M /θ m ≲ 4.7).The number of PSF solid angles per source is thus sufficiently high that we do not consider classical source confusion to be the predominant noise in the RACS-mid images and subsequent catalogues.We note that the 25-arcsec catalogue generally has a larger PSF solid angle than the primary catalogue (except in seven of the full-sensitivity images used for the primary catalogue) and so is overall more affected, though still not limited by classical source confusion.

Source density
Figure 6 shows the HEALPix-binned (N side = 64, ∼ 55 × 55 arcmin 2 ) source and component density across the sky for the primary catalogue [6(i) and 6(ii)], the 25-arcsec catalogue [6(iii) and 6(iv)], and the time-domain catalogue [6(v) and  6 In the primary and 25-arcsec catalogues the Galactic Plane features a lower density of sources, largely due to higher noise on multiple scales.This is a combination of both the real extended sources and the artefacts on multiple angular scales from the poor modelling of the real extended sources.A similar decrease in source density is not as obvious in the time-domain catalogue due to the double-counting of sources in overlap regions.We also see a decrease in source density at the celestial equator corresponding to a reduction in image sensitivity.This is most prominent in the 25-arcsec catalogue.The component density of all three catalogues is higher than the source density, corresponding to the resolved sources being modelled with multiple Gaussian components by PyBDSF.Smaller-scale lowdensity regions can be seen around other bright sources, though generally restricted to the 1.2 degree radius left from peeling described in Paper IV.The regions in all three catalogues around δ J2000 ≈ -75 show significantly higher average source density caused by the large overlap in adjacent observations where the Celestial tiling changes to cover the South Celestial Pole.The time-domain catalogue features increased source density regions with duplicate observations as we do not remove the duplicated sources.
In Figure 7 we also show the source density, binned as a function of declination for six flux density limits (1, 2, 3, 4, 5, and 10 mJy) with the Galactic Plane excluded.The features described for the HEALPix maps are clear in the 1-and 2-mJy flux density populations but largely begin to disappear for the brighter source populations with the 10-mJy source density showing no variation as a function of declination.The exception to this is the time-domain catalogue, for which the 10-mJy population still shows variation in source density as a function of declination-a consequence of including duplicated source entries.

The fraction of unresolved sources
For certain comparisons to external surveys, it is generally preferable to consider isolated and unresolved sources to avoid biases introduced by differences in (u, v) coverage, and by extension angular resolution.While some comparison could be done for resolved/extended radio sources, comparison in that case becomes more heavily affected by differences in (u, v) coverage and frequency.
To create an unresolved subset of the primary catalogue and two auxiliary catalogues, we opt to follow the procedure outlined in section 5.2.1 of Paper II following similar methods employed by e.g., Bondi et al. (2008); Smolčić et al. (2017a); Shimwell et al. (2019).Firstly, this involves assuming the ratio of integrated flux density, S int , and peak flux density, S peak , is close to unity for an unresolved source.In practice, S int /S peak > 1 for a variety of reasons, including: source positions may move during the observation or due to mosaicking adjacent images with small astrometric offsets resulting in a blurred source in final images; self-calibration processes may pull sources in different directions in adjacent beam images (see section 3.4.3 in Paper I or section 3.7.1 in Paper IV), resulting in blurring of the source.In these situations integrated flux density is generally preserved and peak flux density decreases.Other factors such as as bandwidth or time smearing cause similar problems, which for RACS-mid can be up to ∼ 2-4 % reduction towards the edges of primary beam mainlobe, depending on pointing.The blurring, or smearing, of sources more heavily affects those in the low-SNR regime and low-SNR sources may also be pushed to S int /S peak < 1 due to image noise and uncertainty while fitting Gaussian components to sources.
As in Paper II, we begin by assuming that unresolved sources sit within an envelope defined by where C is an estimate of the median S int /S peak in the high-SNR regime unique to each catalogue as it will depend on the average PSF.A and B are found by fitting Equation 3 as in Paper II, after binning values below C and obtaining the 100 -95 th percentile in each bin.
Figure 8 shows the ratio S int /S peak as a function of SNR for single-component sources, along with the fitted envelopes highlighting the resolved sources for each of the three catalogues.The percentage of unresolved sources naturally increases as the catalogue resolution decreases, ranging from 48.4% in the primary catalogue to 58.7% in the 25-arcsec catalogue.Despite the lower rms noise in the RACS-mid data compared to RACS-low, we find a higher fraction of unresolved sources (cf.∼ 40% for RACS-low; Paper II).We suggest this is largely due to the loss in sensitivity to extended sources in the RACSmid images, and the lower number of resolved sources in the 25-arcsec catalogue is further evidence of this.
Identifying unresolved/resolved sources in this way also highlights a subset of the catalogued sources that fall both below S int /S peak < 1 and S int /S peak < C + A × SNR B (i.e.below the envelope).Such sources are unlikely to be real and may  be spurious source-finder detections or artefacts.We therefore add a column in each catalogue called Flag which can be one of the following: Flag = 0 An unresolved source, satisfying: The second condition covers high-SNR sources that have S int /S peak below the median S int /S peak but are not spurious detections.The requirement to include such sources suggests the fitted envelope is not a completely accurate in determining unresolved sources.
Flag = 1 A resolved source, satisfying: Flag = 2 A spurious source, not satisfying the above (i.e. S int /S peak below 1 and below the envelope), likely to be an artefact.

Reliability
Following Paper II, we investigate the reliability of our images and source-finding by quantifying the number of sources detected by PyBDSF below -5σ rms .We assume the noise is symmetric in the image and define the reliability as 1 -N negative /N, where N is the number of sources found above 5σ rms and N negative are sources detected below -5σ rms .We repeat sourcefinding on all images (used for the primary catalogue and the two auxiliary catalogues) after multiplying the images by -1 and setting the source-finding threshold to 5σ rms as in the normal source-finding procedure.Due to the choice of box size for rms calculations, for many images no sources are found in these inverted images.The source-lists generated from the inverted images are then merged following the same process Figures 10(i)-10(iii) show an estimate of the catalogues' reliability for sources outside of the Galactic Plane (|b| > 5 • ), binned as a function of declination for a range of flux density cuts (including 1, 5, and 10 mJy).We find that reliability is not constant in declination.In Table 2 we report median reliability as a percentage of sources without a negative counterpart for the given flux density limits and the full catalogues.For sources outside of the Galactic Plane, we find overall reliability of 99.28 (99.37) %, 99.88 (99.88) %, and 98.68 (98.86) %, for the primary, 25-arcsec, and time-domain catalogues, respectively, with (and without) the 'spurious' sources included.The time-domain catalogue shows the least reliability, which is a consequence of including the high-noise regions at the edges of  images.Conversely, the 25-arcsec catalogue has least number of negative sources and closely reflects the reliability of the RACS-low catalogue, also at 25-arcsec (cf.∼ 99.7% reliability reported in Paper II).We find that the primary catalogue suffers from its increased angular resolution, where smaller-scale, lowsignificance artefacts (not present in the 25-arcsec images) are harder to avoid during rms-thresholding and source-finding.As the reliability is lowest at declinations where the PSF is smallest, it is clear the angular resolution plays a significant role in affecting the reliability of the source-finding and catalogues.
We repeat this comparison for sources within the Galactic Plane (|b| < 5 • ), with medians reported in Table 2 and in Figures 10(iv)-10(vi).Generally there is an increase in negative sources detected in the Galactic Plane.This is most notable in the 25-arcsec catalogue, where reliability gets worse for brighter sources.This is due to artefacts caused by unmodelled (and poorly modelled) extended emission.When convolved to 25-arcsec, these artefacts become significant with respect to the background noise and the source-finder treats them as real sources.We suggest users of the 25-arcsec catalogue (and 25-arcsec images) be cautious when selecting sources in the Galactic Plane.

Completeness
To estimate the completeness of the RACS-mid catalogues, we follow Paper II and inject a realistic sky model into the individual images after subtracting the Gaussian components found and modelled by PyBDSF.For this purpose, we use the semi-empirical sky model from the SKA Simulated Skies project (Wilman et al., 2008;Levrier et al., 2009).We use a base sky model that covers a 100 deg 2 region including sources considered both resolved and unresolved at the angular resolution of the RACS-mid data.For each image, we generate five separate sub-sky models by randomising positions of sources from the 100-deg 2 sky model.We clip the sub-sky models at 3σ rms and at 2 × S max where S max is the maximum flux density of real sources in the image.We do not expect full recovery of sources at the brightest flux densities, as the brightest sources are typically extended at the RACS-mid resolution (e.g.Perley & Butler, 2017) and RACS-mid has poorer sensitivity to extended sources (see e.g.section 3.5 and figure 26 in Paper IV).The semi-empirical SKA sky model may not be accurate below ∼ 0.1 mJy (e.g.Smolčić et al., 2017a;Hale et al., 2023), which is generally below the rms noise in the RACS-mid images.The low flux density limit still allows sources below our detection thresholds to be modelled, but are typically above 0.1 mJy.This provides a total sub-sky model flux density of ∼ 150 Jy for each image.We do not distinguish between unresolved and resolved sources for these tests.
The sub-sky models are injected onto model images and then restored in a similar way to normal imaging: the model is first convolved with the image PSF then the convolved model      is added to the residual map.We use PyBDSF to source-find these model-injected images using the same settings as we used for the real images.The individual output source-lists are then cross-matched to the original sub-sky model catalogues for each image independently before concatenation.For this purpose we use the source list output by PyBDSF rather than the component list.This process is done separately for all image set used to generate the catalogues.
In Figures 11(i)-11(iii) we show the fraction of sources that are cross-matched to the input sky model after source-finding the three RACS-mid catalogues.As in Paper II, the detection fraction is calculated in logarithmically-sampled flux density bins, and we show both cumulative (orange line) and independent (blue line) detection fractions with 50% and 95% detection fractions highlighted.The detection fractions at 50% and 95% are higher than in RACS-low for the primary catalogue and conversely lower for for the 25-arcsec catalogue.The drop in detection fraction for the 25-arcsec catalogue is on account of increased blending of sources at the lower resolution of the 25-arcsec catalogue.
Figures 11(iv)-11(vi) show the completeness of each catalogue as a function of flux density.This is defined as the binned ratio of recovered flux density as a function of flux density.The primary catalogue is found to be 50% complete at 1.1 mJy and 95% complete at 1.6 mJy while we find 1.7 and 2.5 mJy completeness (at 50% and 95%, respectively) for the 25-arcsec catalogue and 1.4 and 5.3 mJy for the time-domain catalogue.Comparing to the 'resolved' source case from Paper II i this suggests the primary RACS-mid catalogue is complete to lower flux densities, accounting for frequency differences (assuming i Where 'resolved' in this case includes all sources, resolved or otherwise.a spectral index j of α = -0.8)whereas the 25-arcsec and timedomain catalogues are less complete due to lower resolution in the 25-arcsec catalogue and the high-noise image edges of the time-domain catalogue.
In addition to the overall completeness, we also repeat the estimate of completeness in 10-degree declination bins for each catalogue.Figures 11(vii)-11(ix) show the 50% and 95% completeness flux density limit (S complete ) as a function of declination.In all cases there are changes in the completeness as a function of declination, with the primary and time-domain catalogues being less complete beyond δ J2000 ≳ +20 • .The 25-arcsec catalogue shows a decrease in completeness around the equator, but with an increase beyond δ J2000 ≳ +20 • .These features correspond well to the variation in the noise properties seen in Figure 4.

Source counts
RACS provides catalogues at sensitivities, angular resolutions, and frequencies not previously available over such wide areas and flux density ranges.This makes the RACS releases valuable in exploring the radio source counts, which are typically restricted to high sensitivity and small areas (e.g.Vernstrom et al., 2016;Smolčić et al., 2017a,b;Mandal et al., 2021;Matthews et al., 2021;Hale et al., 2023) or lower sensitivity and wide areas (e.g.White et al., 1997;Franzen et al., 2021;Matthews et al., 2021;Hale et al., 2021).Paper II reported the normalised differential source counts for RACS-low, and here we do the same for the RACS-mid for both the primary and 25-arcsec catalogues for all sources outside of the Galactic Plane (|b| > 5 • ).
We summarise the source counts (and normalised source counts with and without completeness corrections) in Table 3.The normalised source counts are also shown in Figure 12 for the RACS-mid catalogues, scaled to 1.4 GHz assuming α = -0.8alongside similar source counts from a range of other surveys, including RACS-low, FIRST (reported by White et al. (1997)), the compilation of 1.4-GHz source counts by de Zotti et al. (2010, see references therein), and source counts from the DEEP2 catalogue (Matthews et al., 2021).We also show the 'corrected' RACS-mid source counts, after accounting for the flux density recovery fraction from Figures 11(iv) and 11(v).After applying corrections for incompleteness, we find generally good agreement with existing source counts, with some deviation at low and high flux densities.
At low flux densities (< 2 mJy), the source counts from the primary catalogue follow the low-flux density bins from the de Zotti et al. ( 2010) compilation and the DEEP2 source counts from Matthews et al. (2021) until flux density bins reach ∼ 0.6 mJy.This mirrors what is seen with the RACS-low catalogue.Conversely, the 25-arcsec RACS-mid catalogue shows a decrease in the corrected source counts at ∼ 1 mJy, though this is at the detection limit of the survey.In the highflux density bins, the primary catalogue fails to report some of the brighter sources in the sky.Some of the brightest extended sources such as Fornax A are not recovered well in the data j We define the spectral index as Sν ∝ ν α .as described in section 3.5 in Paper IV and may be split into smaller sources during source-finding.Such bright sources are therefore under-represented in the RACS-mid catalogues.We do not correct these high flux density bins for incompleteness as this requires a complex visibility-based process rather than our image-based sky model injection described in Section 4.6 and is beyond the scope of this brief source count assessment.

Matching to external catalogues
For the following sections, we investigate the brightness scale and astrometry of the primary catalogue and the 25-arcsec catalogue.Due to similarities with the source-lists used for validation in Paper IV, some comparisons with the time-domain catalogue are not repeated.We expect similarities with the primary catalogue and the 25-arcsec catalogue, though due to the additional mosaicking and image-based convolution it is worth confirming we have not introduced any significant systematic effects into the data.
For direct comparisons of flux density measurements and astrometry we look to the NVSS (Condon et al., 1998) at 1.4 GHz, FIRST (Becker et al., 1995;White et al., 1997;Helfand et al., 2015), and the 42-deg 2 region of the WALLABY pre-pilot catalogue (Grundy et al., 2023).FIRST is catalogued at both the same frequency as NVSS (for data prior to 2011) and at 1 335 MHz (for data from 2011).FIRST sources detected in data from 2011 onward comprise ∼ 6% of the full FIRST catalogue but are exclusively in the region that overlaps with RACS-mid.To account for this in the FIRST data, we scale the flux density measurements, S, of the 2011-onward data to match the 1 400-MHz data assuming a power law of the form S ∝ ν α with α = -0.8.For comparison to catalogues at other frequencies, we use the RACS-low catalogue (Paper II), SUMSS (Bock et al., 1999;Mauch et al., 2003), the TIFR k GMRT l Sky Survey alternate data release 1 (TGSS; Intema et al., 2017) catalogue, and the LOFAR m Two-metre Sky Survey data release 2 (LoTSS-DR2; Shimwell et al., 2022) catalogue.These catalogues provide a range of comparisons in the frequency range 144-1400 MHz.For additional astrometry assessment, we also compare to the VLA Sky Survey (VLASS; Lacy et al., 2020) first epoch 'quick-look' catalogue (hereinafter VLASS-QL; Gordon et al., 2021).
For the remainder of the validation work, unless otherwise specified, we are only considering unresolved (Flag = 0) sources in the RACS-mid catalogues.For the comparison catalogues, we generally only look at unresolved sources if the survey angular resolution is on average higher than RACSmid.In the case of RACS-low and the WALLABY pre-pilot catalogue, we use the same unresolved criteria described in Paper II and Grundy et al. (2023) (39.6% and 76.5% unresolved, respectively), but for the remaining comparison catalogues we employ the simpler criteria.For the FIRST, VLASS-QL, and LoTSS-DR2 catalogues we take sources that satisfy | log 10 (S int /S peak )| < log 10 (1.2) as unresolved.In the case of k Tata Insititute for Fundamental Research.l Giant Metrewave Radio Telescope.m LOw Frequency ARray.LoTSS-DR2, because of significantly blurring point sources show ratios closer to ∼ 1.25, and for that catalogue we shift the ratio criterion by an additional log 10 (0.25).With this criterion we find unresolved fractions for the FIRST, VLASS-QL, and LoTSS-DR2 of 44.9, 26.1, and 47.6%.We note for LoTSS-DR2 Shimwell et al. (2022) find a resolved fraction of 8.0%.
In this case we are more concerned with confidence in determining unresolved sources so we remove potentially resolved sources, whereas Shimwell et al. (2022) were focused on resolved sources and conversely remove potentially unresolved sources with their method.Finally, when cross-matching n the catalogues we consider only isolated sources in the respective catalogues.Isolated sources in this case are defined as those with no neighbours within 2θ major for the given catalogue.Cross-matching is then done with a match radius equal to half the PSF major FWHM for the survey with the lowest angular resolution.After crossmatching with these strict criteria we find that the 25-arcsec catalogue results in the largest number of matches in catalogues with angular resolutions ≳ 25 × 25 arcsec 2 .The time-domain catalogue results in the least matches simply because the duplication of sources makes the isolation criterion difficult to satisfy, and generally there will be few cross-matches in the overlap regions between tiles.For assessment of the overlap regions, we refer the reader to Paper IV or Section 4.12.
The numbers of sources cross-matched for each RACS-mid catalogue and external catalogue are reported in Table 4 .

Photometry
As part of the RACS-mid processing, we used cross-matches to NVSS and SUMSS to correct for time-dependent brightness scale variations (see section 2.6 in Paper IV and further discussion of this issue in section 3.3.3 in Paper I).By construction this ensures the bulk brightness scale between RACS-mid and the NVSS and SUMSS is consistent within the assumed spectral index of -0.82 o .Despite this it is still worth comparing the catalogues produced here due to a difference in the images and source-finding.

External brightness scale
For a direct brightness scale comparison, Figure 13 shows the ratio of integrated flux densities of sources cross-matched between NVSS, FIRST, and the WALLABY pre-pilot catalogue for the three RACS-mid catalogues independently as a function of the RACS-mid flux density.Sources in Figure 13 are coloured by S int /S peak as measured in the respective RACS-mid catalogues.Generally S int /S peak has a restricted range as these are intended to be unresolved sources as defined in Section 4.4.Figure 13 also displays a histogram of the distribution of the flux density ratios.
o We note users of the original images can 'unapply' the brightness scale correction factors as described in section 2.6 in Paper IV, with scaling factors retained in the FITS headers of the images available on CASDA.The time-domain catalogue retains the scaling factors which users interested in variability may wish to revert for e.g.comparison with NVSS and/or SUMSS.
The results for the three RACS-mid catalogues are generally similar, and we report in Table 4 the median flux density ratios for each cross-match alongside the median flux density ratio after applying a 10 mJy limit to the RACS-mid catalogue.For all three RACS-mid catalogues, the 10-mJy limit provides median ratios of ∼ 0.99-1.00for the NVSS and FIRST cross-matches broadly consistent with Paper IV.The RACS-mid data shows marginally higher flux densities (∼ 10%) than the WALLABY pre-pilot data, though the WALLABY pre-pilot data is already found to be underluminous with respect to the NVSS with additional positional-dependent variation (For et al., 2021;Grundy et al., 2023).The full sample shows more variation, with the flux density ratios showing systematic deviation, different for the low-resolution NVSS and high-resolution FIRST and WALLABY pre-pilot data.In the case of the NVSS, the low-brightness cross-matches tend towards < 1, which is a result of faint sources in the NVSS catalogue being pushed above the detection threshold due to scatter in the flux densities (i.e.Eddington bias).Conversely, the low-brightness end of the FIRST and WALLABY pre-pilot cross-matches shows a bump > 1 in the ratios before also dropping < 1.The increase in flux density as measured in RACS-mid can be explained by an increase in truly extended sources being cross-matched, which are not as well-detected in FIRST.This indicates our criteria for 'unresolved' sources defined in Section 4.4 has overestimated the true unresolved fraction for sources < 10 mJy.

Spectral indices
For comparisons with surveys at significantly different frequencies, we follow Paper II and opt to assess the spectral indices of sources cross-matched between the RACS-mid primary catalogue and 25-arcsec catalogue and other surveys.For this purpose, we use the RACS-low catalogue (Paper II), the SUMSS catalogue (Mauch et al., 2003), the TGSS ADR1 (Intema et al., 2017), and the LoTSS-DR2 catalogue (Shimwell et al., 2022).A two-point spectral index, α, is calculated for cross-matched sources assuming a power law, following α = log 10 S RACS-mid /S survey log 10 1367.5/νsurvey , where S RACS-mid and S survey are the integrated flux densities of the given source in the RACS-mid catalogue and comparison survey, respectively, and ν survey is the survey frequency.
Following Paper II we also create a subsample by applying a flux density limit to each survey to remove some sensitivity bias-this helps reduce artificial flattening or steepening of the median spectra.We assume limits in the range α = -0.8± 1.2 and set a flux density limit for both the RACS-mid catalogues and comparison catalogues to be either side of that range, depending on which survey is higher in frequency.The resultant limits are reported in Table 4. Table 4 also reports the median spectral indices for the full sample and flux-limited samples for each RACS-mid catalogue and comparison survey.We also plot the histogram of spectral indices (with and without the flux density limit) in Figure 14   The flux-limited subsample is included in the histograms to highlight the reduced number of sources, and the RACS-mid flux density limit is also shown in 2-D histogram for reference (though note the comparison survey limit is not shown).Where relevant the spectral indices of sources matched between catalogues are generally consistent with results from Paper IV.

Astrometry
As noted in Paper IV, there are two astrometric errors that exist in the RACS-mid data (see figures 36 and 38 in Paper IV). Firstly, the self-calibration process, performed independently for each PAF beam image prior to mosaicking, introduces up to pixel-scale offsets in both right ascension (R. A.) and declination.This effect is not dependent on the SNR of the sources in the field.The result of this effect is an uncertainty in the astrometric accuracy up to ±1-2 arcsec.The second effect is a systematic declination-dependent offset in declination which is currently undiagnosed.This effect results in up to ∼ 2 arcsec offset at the highest declination of the survey.While we cannot reasonably correct the beam-to-beam astrometric precision at this stage, we do attempt to remove the systematic declinationdependent offets.
For removal of the declination-dependent declination offset, we cross-match the unresolved time-domain RACS-mid catalogue to sources in the third realisation of the International Celestial Reference Frame (ICRF3; Charlot et al., 2020).The time-domain catalogue is used for this purpose to take advantage of sources in the overlap regions.Since the ICRF3 catalogue is sparse compared to RACS-mid, we do not use the         Figures 15(i)-15(iii) shows the RACS-mid --ICRF3 declination offsets as a function of declination in the top panels, with medians in declination bins shown as blue circles.The upper panels highlight both the lower astrometric precision at high declination (above δ J2000 ≈ 20 • ) and the general trend towards negative offsets at high-and low-declination.We show on the upper panels the selected 2 nd -order polynomial, derived from the time-domain catalogue, and the lower panels show the same data but with removal of the offset model (i.e. the 'corrected' declination).We report the median ∆α and ∆δ against the ICRF3 in Table 5 along with 1σ uncertainties.These remain similar to results from Paper IV.We also report the median declination offset after applying the model correction-this results in no substantial difference to the overall median and uncertainty.
Table 5 also reports similar median offsets for comparisons with a range of other surveys, both with and without the declination correction.Notably, the declination correction does not make a large difference in most cases as the correction is close to the level of the uncertainty (largely driven by the aforementioned beam-to-beam variation).The exception to this is for the comparison to LoTSS-DR2.For LoTSS-DR2, sky coverage begins at δ J2000 > 15 • so suffers from the most extreme declination offsets.A column containing the 'corrected' declination is added to all RACS-mid catalogues, though we retain the original declination column as well.Users may find crossmatching to high-declination regions more straightforward with the corrected declination measurements.
In Figure 16 we show coordinate offsets for the primary and 25-arcsec catalogues against the ICRF3, RACS-low, VLASS-QL, and LoTSS-DR2 as a function of source SNR in the RACSmid catalogues.In Figure 16 we also show binned median Right.Sources within the Galactic Plane (|b| < 5 • ).The distribution of sources that are not cross-matched between the two catalogues are also shown.
offsets as a function of SNR, with 1σ error ranges shown as shaded regions.Note that the uncertainty remains fairly constant as a function of SNR, as mentioned in Paper IV as the beam-to-beam astrometric offset introduced by self-calibration are direction-and source-independent.The lower panels of Figure 16 show the declination offsets after the correction described in Equation 9.The overall median offsets for each catalogue comparison are also reported in Table 5. Generally absolute median offsets are ≲ 0.3 arcsec with the except of the aforementioned declination offset in LoTSS-DR2 and right ascension offsets with RACS-low (∼ +0.4 arcsec) and SUMSS (∼ +0.6 arcsec).A right ascension offset for RACS-low was reported by Paper II and is intrinsic to that catalogue.It is likely this is an elevationdependent offset (the true form of the declination-dependent trend we see here) and with RACS-low being observed off the meridian the offset manifests in both declination and right ascension.

RACS-mid and other ASKAP surveys 4.11.1 RACS-low and RACS-mid at 25 arcsec
The goal of RACS is to create a global sky model for calibration and validation of the other ASKAP surveys and ASKAP observations in general.The global sky model is intended to be a combination of all three RACS bands from 887.5-1632.5 MHz (low, mid, and high).While we are still awaiting RACS-high images and catalogues, we begin with an initial assessment of a sky model built between 887.5 and 1367.5 MHz with the RACS-low and RACS-mid 25-arcsec catalogues.We perform simple cross-match between the source catalogues for each survey, taking sources that cross-match within (θ major θ minor )/4, where θ major and θ minor are the fitted major and minor axes.This variable match criterion is to account for extended sources, which may have different centres in each catalogue.More robust catalogue matching processes may be implemented in future combinations of RACS (e.g. the Positional Update and Matching Algorithm; Line et al., 2017), taking into account expected spectral information.We scale the flux densities of the RACS-low catalogue to 1367.5 MHz assuming the fluxlimit median α = -0.83(see Section 4.9.2).For this crossmatch, we are interested in the region that completely overlaps with both the RACS-low and RACS-mid 25-arcsec catalogues and we only consider sources within the declination range -80 • < δ J2000 < +29 • , including the Galactic Plane.Note that the RACS-low catalogue is missing a small region below δ J2000 = +30 • and has patchy coverage below δ J2000 = -80 • .
Figure 17 shows the distribution of sources of both the RACS-mid and RACS-low 25-arcsec catalogues, as a function of flux density (scaled to 1367.5 MHz) and a function of the deconvolved source size.Also shown are the distributions of sources in each catalogue which do not match to sources in the other catalogue with the aforementioned cross-matching limit.The distributions of sources in the RACS-mid and RACS-low 25-arcsec catalogues is similar both in flux density and source size.RACS-low sources below ∼ 2 mJy (at 1367.5 MHz) begin to fail to cross-match, which is due to the difference in sensitivity between the two catalogues.Conversely, as flux density increases there are more RACS-mid sources without matches.The Galactic Plane sources show a similar increase, particularly between ∼ 10-1000 mJy (and above 1 arcmin in size), which is a consequence of the increase in artefacts as described in Section 4.5.

RACS-mid and WALLABY at the full resolution
We also perform a similar variable-separation cross-match to the WALLABY pre-pilot catalogue (Grundy et al., 2023) using the primary RACS-mid catalogue.For this cross-match, we first extract a cut-out of the RACS-mid primary catalogue that covers the same image region of the WALLABY pre-pilot mosaic.The WALLABY pre-pilot data is entirely outside of the Galactic Plane, so no distinction between Galactic latitudes is done for this comparison.Figure 18 shows the distributions of sources in RACS-mid primary catalogue and the WALLABY pre-pilot catalogue along with non-matched sources as in Figure 17.Due to the increase in sensitivity of the WALLABY pre-pilot data (σ rms ≈ 50 µJy PSF -1 ), there are significantly more unmatched WALLABY pre-pilot sources at the low flux density end of the distribution, though no bias towards specific source sizes in the unmatched WALLABY pre-pilot sources.The RACS-mid primary catalogue, at least outside of the Galactic Plane and away from other particularly complex extended sources (e.g.Fornax A and other large angular extent, nearby radio sources), provides a good interpretation of the 1367.5 MHz sky.

The brightness scale of the full-sensitivity images
While the main data products released with this work are the 'all-sky' catalogues, we also provide the full-sensitivity, convolved maps generated prior to source-finding.With the exception of the time-domain catalogue (which uses the original images), these new images are mosaicked with nearby observations.This additional mosaicking provides a more uniform sensitivity pattern across the images, and some of the edge effects of each tile image discussed in Paper IV, particularly for the brightness scale, can be reduced during this process.
For assessing the brightness scale across the survey images, we apply the same unresolved flags described in Section 4.4 to the individual source lists then cross-match the individual sources lists to the various external catalogues as described in Section 4.8.Note that more sources will be cross-matched in this per-source-list assessment than when looking at the full catalogues, since we include sources in at the edge of images in overlap regions.For assessing the brightness scale, we look at the SUMSS and NVSS cross-matches to allow easy reference to the same comparisons made in Paper IV for the original images.As part of the additional metadata added prior to source-list concatenation, we add tile-reference coordinates, (l, m), for each source entry.In the case of the mosaicked images, these (l, m) values are with reference to the mosaiced image, and not the original image that was used in the mosaic.Following Paper IV we then median-bin the flux density ratios (S RACS-mid /S survey , for NVSS and SUMSS) in bins of angular sizes 8.6 × 8.6 arcmin 2 (NVSS) and 14.4 × 14.4 arcmin 2 (SUMSS).
Figures 19(i)-19(iv) shows (left) the flux density ratios of sources in the RACS-mid catalogues and comparison catalogues, median-binned in the image-based (l, m) coordinates.Figure 19 also shows the 2-D histogram of flux density ratios in m, where the majority of the edge-effects were seen in Paper IV.We normalise the ratios by the overall median flux density ratio to ensure bulk offsets introduced by a choice of spectral index do not detract from position-dependent effects.For both the primary and 25-arcsec catalogues, position-dependent brightness effects are almost completely removed due to the mosaicking process.Specifically, we refer to figures 30 and 31 in Paper IV that highlight these residual brightness scales features.Residual variation in the mosaiced images is ≲ a few %.
With the reduction of edge effects, we determine that the brightness scale uncertainty for the primary and 25-arcsec catalogues retain the uncertainty reported for the central tile region in Paper IV: 6%.The time-domain catalogue retains both uncertainties (6% and 14% for the tile edges, equation 5 in Paper IV), and the source entries will explicitly note the appropriate uncertainty based on its position in the image it is taken from.

Data availability
Paper IV describes the raw survey data and how to access the original images and calibrated visibility datasets.This data release comprises both the full-sensitivity mosaicked images used for the primary catalogue and the fixed-resolution 25-arcsec catalogue, as well as the three catalogues described in Section 3.For the catalogues and full-sensitivity images described in this paper, availability is similarly through the the CSIRO p ASKAP Science Data Archive (CASDA; Chapman et al., 2017;Huynh et al., 2020), and each product is summarised here for completeness: Primary catalogue.
• Source list.The final merged source lists.This contains all sources detected during source-finding on the full-sensitivity images that are not removed due to the de-duplication process when constructing the catalogue.This catalogue features the full sky up to δ J2000 = +49 • with sources detected in Stokes I images.This is the recommended catalogue for most use cases.• Component list.The Gaussian components associated with the source list.Generally only used for reconstruction of the source models.• Full-sensitivity images.Stokes I full-sensitivity images with associated weight maps.For most use-cases these are the preferred images to use over the original images.• Link to images and catalogues. https://doi.org/10.25919/p524-xb81.
• Source list.The merged 25-arcsec source lists.This is similar to the primary catalogue, except with the fixed 25arcsec resolution and with coverage up to δ J2000 = +30     Auxiliary time-domain catalogue.
• Source list.The concatenated source list from the original RACS-mid images.This covers the entire survey region, and retains duplicate source entries.• Component list.The Gaussian components associated with the time-domain source list.• Link to the catalogues. https://doi.org/10.25919/p8ns-da63.

Summary
The Australian SKA Pathfinder (ASKAP) has completed observing the sky in three bands centered at 887.5, 1367.5, and 1632.5 MHz as part of the Rapid ASKAP Continuum Survey (RACS).RACS-low (at 887.5 MHz) and RACS-mid (at 1367.5 MHz) have been processed and released (Paper I; Paper IV) with a catalogue for RACS-low also released (Paper II).While the 1632.5-MHzdata have been observed, these data are yet to be processed and will be released in the future as RACS-high.This work describes the release of an 'all-sky' general-purpose catalogue for RACS-mid, along with two auxiliary catalogues for more specific use-cases.The general purpose, primary catalogue was created by mosaicking nearby images to improve sensitivity in overlapping regions, convolving to the lowest common resolution of neighbouring images, source-finding with PyBDSF on these full-sensitivity mosaics, then merging the output source-lists while removing duplicate entries.The primary catalogue covers -90 • < δ J2000 ≤ +49 • (∼ 36 200 deg 2 ) and contains 3 107 143 sources (2 863 400 excluding the Galactic Plane) down to 5σ rms .We investigated the reliability and completeness of the catalogues, finding a significant decrease in reliability in the Galactic Plane due to higher incidence of artefacts from extended sources.The primary catalogue is found to be 95% complete at 2 mJy.We matched the catalogue to other external catalogues such as RACS-low, LoTSS-DR2, VLASS-QL, SUMSS, NVSS, and FIRST and find properties such as flux density and astrometry consistent with what is reported in Paper IV, suggesting additional mosaicking done here is not introducing extra systematic issues into the images and catalogue.
In addition to the primary catalogue, we also provide an auxiliary version of the catalogue with a fixed resolution of 25 ′′ × 25 ′′ .This catalogue is created in a similar way to the primary catalogue, except images are all convolved to 25 ′′ × 25 ′′ prior to mosaicking and source-finding.Because of this, the coverage is reduced to δ 2000 ≤ +30 • .This matches closely the existing RACS-low catalogue described in Paper II.This 25arcsec catalogue features lessened sensitivity and consequently fewer sources-2 156 393 sources above 5σ rms (1 992 406 excluding the Galactic Plane)-with a significant decrease in reliability in the Galactic Plane .Other properties and comparisons are similar to the primary catalogue.
A third auxiliary catalogue is provided for time-domain science.This version features no additional convolution or mosaicking of images prior to source-finding and retains all duplicate source entries in overlapping regions.The purpose of this auxiliary time-domain catalogue is to detect and characterise time-variable sources without averaging signal between epochs as will be the case for some sources in the primary catalogue.
All data products, including catalogues and full-sensitivity images are being made available through the CSIRO ASKAP Science Data Archive.This data release focuses on total intensity with Stokes I and work is currently underway to catalogue the sources of circular polarization using the Stokes V images produced as part of RACS-mid (e.g. Rose et al., 2023).

Acknowledgement
We would like to thank the anonymous referee for their feedback that has improved the quality of this work.This scientific work uses data obtained from Inyarrimanha Ilgari Bundara / the Murchison Radio-astronomy Observatory.We acknowledge the Wajarri Yamaji People as the Traditional Owners and native title holders of the Observatory site.CSIRO's ASKAP radio telescope is part of the Australia Telescope National Facility (https://ror.org/05qajvd42).Operation of ASKAP is funded by the Australian Government with support from the National Collaborative Research Infrastructure Strategy.ASKAP uses the resources of the Pawsey Supercomputing Research Centre.Establishment of ASKAP, Inyarrimanha Ilgari Bundara, the CSIRO Murchison Radio-astronomy Observatory and the Pawsey Supercomputing Research Centre are initiatives of the Australian Government, with support from the Government of Western Australia and the Science and Industry Endowment Fund.
This project used public archival data from the Dark Energy Survey (DES).Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Government grant NAG W-2166.The images of these surveys are based on photographic data obtained using the Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope.The plates were processed into the present compressed digital form with the permission of these institutions.
RACS images of radio galaxies near the galaxy cluster Abell 2829.i, r, g) (iv) RACS images of NGC 7582.

Figure 1 .
Figure 1.Example cutouts of sources in RACS images.The images are the RACS-mid full-sensitivity image (left), the RACS-mid 25-arcsec image (centre left), the RACS-low 25-arcsec image (centre right, Paper II), and a three-color DES DR2 image (i, r, g assigned to RGB, right).The radio colour scales follow a square-root stretch between [-2, 200]σrms.The ellipses in the lower left of the RACS images show the size of the PSF.5σrms contours from the left panel are drawn on the DES images for reference.

Figure 2 .
Figure 2. Examples cutouts of high-declination sources.The images are the RACS-mid full-sensitivity image (left), the NVSS (centre left), FIRST (centre right), and a three-color DSS2 or SDSS DR7 image right).The radio colour scales follow a square-root stretch between [-2, 100]σrms.The ellipses in the lower left of the RACS images show the size of the PSF.5σrms contours from the RACS-mid image are drawn on the DSS2 and SDSS DR7 images for reference.

Figure 3 .
Figure 3. Variation of the major (i) and minor (ii) axes of the PSF over the RACS-mid primary catalogue (red), the 25-arcsec catalogue (purple), and the time-domain catalogue (blue).The filled markers are medians within 4-deg bins, and the shaded regions show the range of PSF axes within the given bins.

Figure 4 .Figure 5 .
Figure 4. Median-binned HEALPix representation of the root-mean-square noise distributions of the primary catalogue (i), the 25-arcsec catalogue (ii), and the time-domain catalogue (iii).The black, dashed lines are drawn at Galactic latitudes b ± 5 • .

Figure 6 .
Figure 6.HEALPix representation of the source and component density of the primary catalogue [(i) and (ii)], the 25-arcsec catalogue [(iii) and (iv)], and the time-domain catalogue [(v) and (vi)].The black, dashed lines are drawn at Galactic latitudes b ± 5 • .Regions with no sources are coloured grey.The colourscales are kept the same to highlight the differences in source/component densities.

Figure 7 .
Figure7.Source density as a function of declination for six flux density limits([1, 2, 3, 4, 5, 10]  mJy) for the primary catalogue (i), the 25-arcsec catalogue (ii), and the time-domain catalogue (iii).Dashed lines correspond to median source densities for the associated flux density limit.The grey, shaded region in (ii) is not covered in the 25-arcsec catalogue.
(vi)].The average source density across the primary catalogue h Murchison Widefield Array. is ∼ 86 deg -2 , over the catalogued area of ∼ 36 200 deg 2 .The 25-arcsec catalogue has a source density of ∼ 70 deg -2 (coverage ∼ 30 900 deg 2 ) and the time-domain catalogue has a source density of ∼ 113 deg -2 (coverage ∼ 36 200 deg 2 , though note that by construction there are many duplicate sources in the time-domain catalogue).
Primary catalogue.
Time-domain catalogue.

Figure 8 .
Figure 8.The ratio of integrated flux density (S int ) to peak flux density (S peak ) as a function of signal-to-noise ratio (SNR = S peak /σrms).The left panels show binned 2-D histograms of the sources while the right panel shows a zoomed-in region with unresolved sources coloured blue.Red crosses in the right panel indicate the binned values used to fit the lower envelope.The lower and upper envelope are drawn as solid black lines, and the median S int /S peak value is shown as a black, dashed line.The single solid, black horizontal line indicates S int /S peak = 1.
Time-domain catalogue.

Figure 9 .
Figure9.Histograms of the negative (pink) and positive (blue) components found over the survey images used in the primary (i), 25-arcsec (ii), and time-domain (iii) catalogues.

Figure 10 .
Figure 10.Reliability for the three catalogues as a function of declination for different flux density cuts.Empty markers correspond to the subsample that excludes the Flag = 2 'spurious' sources.Horizontal lines correspond to medians for each flux density limit (dashed lines for the full samples, and dotted lines for the subsample).Top row.[(i)-(iii)] Sources outside of the Galactic Plane (|b| > 5 • ).Bottom row.[(iv)-(vi)] Sources within the Galactic Plane (|b| < 5 • ).Note all plots show a different y-axis scale to highlight different declination-dependent features.

Figure 11 .
Figure 11.Assessment of the completeness of the three RACS-mid catalogues.(i)-(iii): The fraction of detected model sources matched to the input SKA sky model.The red lines correspond to matched fraction in flux density bins, while the blue line corresponds to the integrated matched fraction.50 and 95 % matched fractions are indicated on the plots.(iv)-(vi): An estimate of completeness as the fraction of flux density recovered in the model images as a function of flux density.50 and 95 % completeness fractions are shown on the plots.(vii)-(ix): An estimate of the 50% and 95% completeness (S complete ) as a function of declination, with the mean completeness shown as horizontal lines.
Figure 12. corrected and uncorrected source counts from the primary and 25-arcsec RACS-mid catalogues (excluding |b| < 5 • ).Source counts from a selection of other surveys are shown for reference, including the 1.4-GHz compilation from de Zotti et al. (2010), the 1997 FIRST catalogue White et al. (1997), Matthews et al. (2021), and RACS-low (Paper II).
alongside a 2-D histogram showing RACS-mid flux density as a function of spectral index.

Figure 13 .
Figure 13.Flux density ratio of unresolved sources cross-matched in the RACS-mid catalogues and the NVSS (top panels), FIRST (centre panels), and the WALLABY pre-pilot catalogue (bottom panels) as a function of the RACS-mid flux density measurements.Background data points are coloured by their SNR in the RACS-mid data.The solid, black line in each panel indicates a ratio of 1. Pink points indicate median flux density ratios in bins, with the shaded pink region covering the 16 th and 84 th percentiles in each bin.

Figure 16 .
Figure 16.R. A. (∆α cos (δ)) and declination (∆δ) offsets as a function of SNR for the primary catalogue [(i)-(iv)] and the 25-arcsec catalogue [(v)-(viii)].Left.ICRF cross-match.Centre-left.RACS-low cross-match.Centre-right.VLASS-QL cross-match.Right.LoTSS-DR2 cross-match.In each panel, a dashed, grey line is drawn at ±2 arcsec corresponding to the pixel size of the original RACS-mid images.Median values in SNR bins are also shown as coloured markers and the shaded regions indicate ±2σ for the SNR bins.

Figure 17 .
Figure 17.The distribution of sources as a function of flux density (top) and fitted major axis (bottom) for the RACS-mid 25-arcsec catalogue and the RACS-low catalogue.Left.Sources outside of the Galactic Plane (|b| > 5 • ).Right.Sources within the Galactic Plane (|b| < 5 • ).The distribution of sources that are not cross-matched between the two catalogues are also shown.

Figure 18 .
Figure18.As in Figure17, comaparing the RACS-mid primary catalogue to the WALLABY pre-pilot catalogue (only considering sources within the WALLABY pre-pilot image region).No Galactic Plane cut is applied as the region is not within the Galactic Plane.

Figure 19 .
Figure 19.Brightness scale comparison as a function of tile reference coordinates for the primary catalogue and NVSS (i) and SUMSS (ii), and the 25-arcsec catalogue and NVSS (iii) and SUMSS (iv).Left.Median-binned flux density ratios.Right.2-D histogram of the flux density ratios and the tile coordinate m, with medians in m show as pink circles.16 th and 84 th percentiles are shown as pink, dashed lines.

Table 1 .
RACS-mid all-sky catalogues and image properties.

Table 2 .
Reliability estimates for the RACS-mid catalogue, for different flux density limits, S limit .

Table 3 .
Tabulated source counts for the primary and 25-arcsec RACS-mid catalogues.

Table 4 .
Median reported properties of sources cross-matched between the RACS-mid catalogues and other surveys, with and without flux density limits applied.
a Frequency of the comparison survey.b Sources detected in data from 2011 are scaled from 1 335 MHz assuming a power law with spectral index -0.7.
Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft, and the Collaborating Institutions in the Dark Energy Survey.The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciències de l'Espai (IEEC/CSIC), the Institut de Física d'Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, The Ohio State University, the OzDES Membership Consortium, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University.Based in part on observations at Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.

Table A1 .
Columns in the RACS-mid source catalogues.Uncertainty in the peak flux density from PyBDSF fitting.Quadrature sum of the brightness scale and PyBDSF uncertainties for peak flux density.

Table A2 .
Columns in the RACS-mid component catalogues.Gaussian_ID Unique identifier for the component based on the Field_ID.1Source_IDUniqueidentifierforthe source the component belongs to (see TableA1).Uncertainty on Dec_corr from fitting the component position and declination offset model.Uncertainty in the peak flux density from PyBDSF fitting.Quadrature sum of the brightness scale and PyBDSF uncertainties for peak flux density. •

Table A3 .
Example three rows from the primary catalogue source list.Rows with associated uncertainties are rounded to match the uncertainty.

Table A4 .
Example three rows from the primary catalogue component list.Rows with associated uncertainties are rounded to match the uncertainty.