We describe the latest iteration of upgrades (designated Phase III) to the Murchison Widefield Array (MWA), in the fourth paper in a series that covers the evolution of the telescope from design concept to initial operational facility, and through two major upgrades. As part of the Phase III upgrade of the MWA, we report the completion of work to design, build, and deploy a new fleet of digital receivers that further optimise the MWA for Epoch of Reionisation observations. These receivers complement existing receivers, such that the MWA now supports the full correlation of all 256 antenna tiles currently in the array. This step releases the MWA from the prior constraint of having to correlate only 128 of the 256 tiles at any given time, which means that the maximum instantaneous sensitivity of the MWA is doubled and the maximum number of interferometric baselines is approximately quadrupled. The upgrade is fundamentally enabled by the new MWAX correlator and various other improvements to the MWA sub-systems. In this paper we describe the new digital receivers and the other improvements that result in the Phase III system. A range of operational benefits arise from the upgrade and scientific flexibility is increased. We also comment on the transition from the MWA to the SKA-Low facility near the end of the decade, including a description of some unique science opportunities utilising joint MWA/SKA-Low data during the Science Verification phase of the SKA-Low Array Assembly 2 (AA2) period.

Understanding the evolution of the cosmic star formation rate density (SFRD) is key to uncovering how the Universe arrived at its present state. This paper presents a novel and efficient method to estimate the local SFRD, which uses supervised machine learning to first identify a population of star-forming galaxies (SFGs). Next, star-formation rates (SFRs) are determined using the 1.4-GHz radio-continuum emission detected by the Australian Square Kilometre Array Pathfinder (ASKAP). Specifically, a gradient-boosted decision tree model was imple-mented to classify extragalactic sources from the Beck et al. (2022) catalogue as either galaxies or quasars using RACS-mid and WISE photometry. The full sample, consisting of 389,392 sources, was partitioned into a 70%-15%-15% split for training, validating, and testing. The optimised model achieved a weighted F1 score of 0.93 and an accuracy of 0.94 on the test dataset, ultimately classifying 336,674 sources as galaxies and 52,718 sources as quasars. Using the resulting z <0.1 depth-matched galaxy sample and the photometric redshift predictions from Beck et al. (2022), a modified 1.4-GHz SFR calibration was determined, yielding a local, completeness-corrected, z < 0.1 SFRD of (1.4 ± 0.5) × 10−2 M⊙ yr−1 Mpc−3 using 11,293 sources. This value is consistent with previous results. Thus, this study demonstrates the feasibility of using supervised learning to identify large populations of SFGs in order to investigate the SFRD evolution. This presents an exciting prospect for future, deeper surveys such as EMU, which will enable the cosmic SFRD to be probed out to higher redshifts.

The Southern-sky MWA Rapid Two-metre (SMART) survey, which capitalises on the MWA’s large field of view and voltage recording capability, is an ambitious effort to conduct sensitive searches for pulsars and fast transients in the 140–170MHz band. The novelty of voltage recording, long dwell times (4800 s) and the high-time and -frequency resolutions (100 μs/10-kHz) exchange a large survey speed (∼ 450 deg2 h–1) for high computational cost. The survey covers the entire sky south of +30° in declination through a series of dedicated observing campaigns, accumulating nearly four petabytes of data. The large volumes of data and the processing challenges at low frequencies necessitate data processing to be approached in multiple phases, and the initial searches focused on a first-pass (shallow) survey of parts of the skies, as reported in earlier papers in this series. These data are also processed for re-detections of hundreds of known pulsars in the southern sky, many of which are also the first detections at frequencies below 400 MHz. This paper is motivated by the need to address the inherent difficulties (for the wider community) in handling large amounts of voltage data and software/processing challenges for routine pulsar detections, and also by the fast-evolving landscape of the SKA Observatory (SKAO).With the construction and commissioning ramping up towards the full-scale SKA-Low, a low-frequency catalogue of detectable pulsars in the southern sky will prove to be a valuable reference for the science verification exercise. A growing sample of low-frequency pulsar detections and measurements will also prove invaluable in a variety of science applications including population studies, survey simulations and emission beam models, refining interstellar medium models for electron densities and the spatial distribution of turbulence, and also for forecasting the detection prospects and survey yield from pulsar surveys planned with SKA-Low. We also present an electronic catalogue of various data products, including pulse profiles, time series and multi-channel folded archives, along with the measurements of dispersion and rotation measures, and mean flux densities for the detected pulsars, and this will be periodically updated as more detections flow on from the ongoing data processing.

Measurements of peculiar velocities in the local Universe are a powerful tool to study the nature of dark energy at low (z < 0.1) redshifts. Here we present the largest single set of z < 0.1 peculiar velocity measurements to date, obtained using the Fundamental Plane (FP) of galaxies in the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI).We describe the photometric and spectroscopic selection criteria used to define the sample, as well as extensive quality control checks on the photometry and velocity dispersion measurements. Additionally, we perform detailed systematics checks for the many analysis parameters in our pipeline. Our DESI DR1 catalogue contains FP-based distances and peculiar velocities for 98, 292 unique early-type galaxies, increasing the total number of z < 0.1 FP distances ever measured by a factor of ∼ 2. We achieve a precision of 26% random error in our distance measurements which is comparable to previous surveys. A series of companion DESI papers use the distances and peculiar velocities presented in this paper to measure cosmological parameters.

We present a new data release from the Rapid ASKAP Continuum Survey (RACS), a widefield snapshot radio survey conducted with the Australian SKA Pathfinder (ASKAP). This data release contains the second RACS epoch to make use of ASKAP’s low-frequency band, centred on 887.5MHz with a bandwidth of 288 MHz, referred to as RACS-low2. This RACS-low2 data release includes both Stokes I and V imaging and catalogue data products covering the whole sky up to a declination of ≈ +48°. RACS-low2 largely follows the observation footprint of the first low-band epoch, though includes additional coverage in the Northern Hemisphere. The observation scheduling and data processing follow previous RACS epochs, making use of autonomous scheduling and holography-derived primary beam models. The Stokes I and V catalogues are derived from images with a median point-spread function (PSF) of 14.2′′ × 11.8′′, and have median root-mean-square noise properties of ≈ 195 and ≈ 163 μJy PSF−1, for the Stokes I and V images, respectively. The consolidated Stokes I catalogue contains 3 922 151 sources. We also constructed a catalogue of Stokes V sources by measuring Stokes V images at the Stokes I source positions. For the 221 Stokes V measurements above the estimated leakage and detection thresholds we provide likely identifications, including detections of 61 radio stars, 85 pulsars, 43 AGN (many likely to be residual leakage), and one source that is not associated with a known astronomical object. These data products, including calibrated visibilities, images, source lists, and consolidated catalogues, are made publicly available through the CSIRO ASKAP Science Data Archive (CASDA).

Fast radio bursts (FRBs) emit broadband radio emission that may, in rare cases, encode atomic hydrogen (Hi) absorption signals as they traverse the interstellar medium of their host galaxies. Though considered in the early FRB literature, the demanding observational prerequisites and the rarity of suitable events have meant that no thorough search for Hi absorption in FRB spectra has yet been undertaken. Here, we present an updated systematic analysis assessing the likelihood of modern facilities to detect such absorption features. As a proof of concept, we search for absorption in the spectrum of the bright ASKAP-localised FRB 20211127I, finding a 3σ opacity upper limit of 0.51. While this test case offers little constraining power, we find that narrow FRBs with fluences exceeding 20/70/150 Jy ms observed with MeerKAT/ASKAP/DSA can probe opacities below 0.1 — a regime in which absorption detections become physically meaningful. We further highlight that stacking thousands of bursts from hyperactive repeaters with FAST offers a very powerful avenue toward detection. Finally, we discuss the broad scientific potential of such detections, including constraints on extragalactic Hi spin temperatures, a means to physically probe the environment surrounding the progenitor, and a path towards disentangling host galaxy contributions to dispersion and scattering.

We present the Galactic Plane Monitor (GPM), a dedicated low-frequency radio transient survey conducted with the Murchison Widefield Array at 200 MHz. The GPM targets an underexplored region of transient parameter space spanning timescales from seconds to weeks, where recent discoveries of unusual Galactic radio sources—including long-period transients and other coherent emitters—suggest a rich and largely uncharted population. Leveraging the MWA’s wide field of view, we carried out 72 observing epochs between June 2022 and March 2025, searching ∼ 5,000 sq. deg. of the Southern Galactic plane. We developed a multi-stage processing pipeline to search for variability across a broad range of timescales, from intra-snapshot (4 s to minutes), to inter-snapshot (minutes to hours), and inter-epoch (days to years). This approach combines image differencing, statistical variability metrics, and visual inspection to robustly identify transient candidates while mitigating the effects of ionospheric variability and imaging artefacts. In parallel, we produced deep mosaics using a “lucky imaging” strategy to maximise sensitivity to faint persistent sources. We present the detections of five long-period radio transients, seven eclipsing binary millisecond pulsars, seven intrinsically variable pulsars, and the detection of refractive interstellar scintillation and night-time interplanetary scintillation. The GPM represents a key step toward fully mapping the dynamic low-frequency radio sky and informs the design of next-generation transient surveys with the SKA Observatory.

Searches of unassociated gamma-ray sources in the Fermi-LAT catalogues have led to the discoveries of around a fifth of all known millisecond pulsars (MSPs). These searches have almost exclusively been performed at radio frequencies above 300MHz, where dispersion and scattering in the interstellar medium are less significant. We report on a shallow survey for pulsars targeting 308 unassociated Fermi-LAT sources in archival MurchisonWidefield Array (MWA) observations from the Southern-sky MWA Rapid Two-metre (SMART) pulsar survey at 154MHz. This is the largest radio survey of unassociated Fermi-LAT sources to date, and only the second to be conducted below 300MHz after a survey with the Low Frequency Array (LOFAR) that discovered three MSPs. Each source was observed for 20 min by digitally beamforming the MWA tile voltages. Searches were then performed using a new pipeline that implements a semi-coherent dispersion removal scheme for MWA data, enabling greater sensitivities to MSPs than is possible with fully-incoherent dispersion removal (e.g. 2–3 times better sensitivity for dispersion measures between 20–40 cm−3 pc). The pipeline was tested by blindly detecting five known MSPs, four of which are in short-orbit binaries. No new pulsars were identified in the survey, which we attribute to insufficient sensitivity. We estimate flux density limits of approximately 30–220 mJy at 154MHz (or 0.7–5.2 mJy at 1.4GHz) for a spin period of 2 ms and a duty cycle of 28%, with a dependence on the sky temperature and the offset from the phase-centre of the primary beam. We discuss how the improved instantaneous sensitivity from the Phase III upgrade of the MWA will increase the number of detectable gamma-ray pulsars by ∼ 30% for the same integration time. Additionally, the real-time beamformer (under development) will enable longer observations with sensitivities that are more competitive with previous surveys of Fermi-LAT sources. The semi-coherent search pipeline we have developed will also be useful for searches of supernova remnants, globular clusters, and pulsar candidates identified in imaging surveys, all of which will help to inform the significance of future surveys with SKA-Low.

The stellar initial mass function (IMF) is a fundamental ingredient in galaxy evolution, linking observed integrated light to stellar masses, star-formation rates, and chemical enrichment histories. Constraining the full IMF shape beyond the Milky Way remains challenging, as most studies focus either on the low-mass end of quiescent galaxies or the high-mass end of star-forming galaxies. Here we present the first simultaneous analysis of both ends of the stellar initial mass function (IMF) in 214 star-forming galaxies from the Hector survey (z ∼ 0.01–0.07). We estimate the low-mass end slope (αlow) using a stellar population approach that fits IMF-sensitive absorption features with extended star formation histories, while the high-mass end slope (αhigh) is derived via the Kennicutt diagnostic, which compares the observed Hα equivalent width and g – r colour with stellar population synthesis models. We find substantial diversity in IMF shapes, with galaxies spanning combinations of bottom-heavy/light and top-heavy/light slopes. A weak but statistically robust correlation between the low- and high-mass IMF slopes is observed, but partial correlation analysis indicates that this apparent link is largely driven by their mutual dependence on stellar mass and metallicity. Both IMF slopes show significant correlations with stellar mass, star formation activity (traced by Hα luminosity and surface density), and stellar metallicity ([M/H]). In general, higher stellar mass, stronger star formation activity, and higher metallicity are associated with both bottom-heavy and top-heavy IMFs. We find that the full IMF shape seems to be modulated by total stellar mass. Partial correlation analysis reveals that αlow is primarily driven by [M/H], whereas αhigh is mainly linked to stellar mass and recent star formation. Because αlow traces the IMF over long-term averages and αhigh captures only recent (≲ 10 Myr) star formation, the processes shaping each end likely occur over different and possibly decoupled timescales. Our findings challenge the universality of the IMF and emphasise the need for galaxy evolution and stellar population models to incorporate a flexible IMF prescription. Accounting for these variations is essential to build an IMF-consistent picture of galaxy evolution across cosmic time.

We present 29 successfully recovered Civ time lags in Active Galactic Nuclei from the complete Dark Energy Survey Reverberation Mapping campaign. The AGN in this sample span a redshift range of 1.9 < z <3.5. We successfully measure the velocity dispersion from the Civ spectral linewidth for 25 of these 29 sources, and use these to calculate new high-redshift black hole mass estimates, finding masses between 0.8 and 1.3 billion solar masses. We also identify a selection effect due to the duration of the survey that can impact the radius-luminosity relation derived from this and other (high-redshift) data. This paper represents the culmination of the OzDES Civ campaign.

Over a hundred fast radio burst (FRB) host galaxies have now been identified, enabling both comparisons of host redshift with FRB dispersion measure to study the cosmological distribution of ionised gas, and analyses of host properties in order to identify FRB progenitors. The standard method for determining the most likely FRB host galaxy in an optical image is the Bayesian framework Probabilistic Association of Transients to their Hosts (PATH), which accounts for uncertainties in the radio localisation, and simplified prior distributions on the host being observable. In this work we extend PATH, incorporating physically-motivated priors that are based on expectations about FRB host galaxy magnitudes. We develop three different models for the apparent r-band magnitude distribution based on an FRB’s expected host galaxy redshift, P(mr|z) and combine these with expectations for redshift based on an FRB’s dispersion measure, P(z|DM). We fit the parameters of these prior models using host galaxy candidates for 32 FRBs detected by the Australian SKA Pathfinder (ASKAP) in incoherent sum (ICS) mode by the Commensal Real-time ASKAP Fast Transients (CRAFT) survey.

Employing PATH with the new priors on the host magnitudes, we find increased confidence in the most probable hosts of all ASKAP ICS FRB host galaxies. All three models predict similar distributions of FRB host magnitudes at low redshift (z ∼ 0.1), and we confirm previous results that the true FRB host galaxy distribution is fainter than expected for a star-formation-weighted distribution (p-value of 0.12%). However, a mass-weighted distribution provides an even worse fit (p-value of 10−9). Tests against more FRBs in the z > 0.5 range, where the models differ, and extensions of the models to account for e.g. host metallicity, may help to resolve these uncertainties in the FRB host distribution.

We investigate how galaxy evolution varies with environment in the nearby Universe by comparing an “average” reference volume in the Southern Galactic Pole (SGP) dataset from Van Kempen et al. (2024) to the Nexus region, a dynamically assembling superstructure centred on the Abell 4038 galaxy cluster. We quantify environmental effects using the quenched fraction (fQ) and the mean and scatter of the specific star formation rate (sSFR) for the star-forming population, measured as functions of stellar mass and improved group-scale halo mass estimates from Van Kempen et al. (2026). One-dimensional binned trends are summarised with logistic fits for fQ and a power-law describes the binned mean log sSFR trends. We also de-couple the stellar–halo mass dependence in both the fQ and mean log sSFR analyses, demonstrating a joint dependence: fQ increases with stellar mass in both field and group environments, while group galaxies show an additional dependence on halo mass. The Nexus exhibits systematic differences relative to the SGP baseline, consistent with increased heterogeneity in accretion histories and pre-processing within a forming superstructure. For star-forming galaxies, the mean log sSFR declines strongly with stellar mass and shows additional environment-linked suppression in group-scale halos, while the scatter in log sSFR varies primarily with stellar mass and shows comparatively weaker dependence on halo mass. However, these differences are generally within the uncertainties, and larger samples of dynamically evolving Nexus-like structures are required to determine whether they reflect genuine environmental effects or cosmic variance. Within the Nexus, splitting the sample into three projected radial zones around Abell 4038 shows that both quenching and the properties of the star-forming population vary systematically with distance from the node and forming supercluster, largely driven by differences in the sampled halo mass function, indicating that environmental regulation is not spatially uniform across the structure. Finally, a projected phase-space (PPS) analysis of Abell 4038 shows higher fQ in regions associated with earlier infall, linking quenching to orbital history within the cluster. However, when grouping infall PPS zones and splitting by stellar mass, this trend is strongly mass dependent, with low-mass galaxies (logMstellar < 10) showing no significant change in fQ. These results demonstrate that the drivers of galaxy evolution depend jointly on stellar mass, the halo mass of the local group environment, and location within the surrounding large-scale structure. This motivates future, larger-statistics, multi-wavelength studies that combine wide-area spectroscopy with tracers of gas, dust, and hot halos to test quenching and star formation regulation mechanisms across the cosmic web.

Over the last decade, the Australian Dark Energy (OzDES) collaboration has used Reverberation Mapping to measure the masses of high redshift supermassive black holes. Here we present the final review and analysis of this OzDES reverberation mapping campaign. These observations use $6-7$ years of photometric and spectroscopic observations of 735 Active Galactic Nuclei (AGN) in the redshift range $z\in [0.13, 3.85]$ and bolometric luminosity range $\log_{10}(L_{\mathrm{bol}})\in [44.3, 47.5] \; \mathrm{erg/s}$. Both photometry and spectra are observed in visible wavelengths, allowing for the physical scale of the AGN broad line region to be estimated from reverberations of the H$\beta$, MgII and CIV emission lines. We successfully use reverberation mapping to constrain the masses of 62 super-massive black holes, and combine with existing data to fit a power law to the lag-luminosity relation for the H$\beta$ and MgII lines with a scatter of $\sim0.25$ dex, the tightest yet identified, fit specifically for consistency with high redshift AGN. We fit a similarly constrained relation for CIV, resolving a tension with the low luminosity literature AGN by accounting for selection effects arising from finite survey length. We also examine the impact of emission line width and luminosity (related to accretion rate) in reducing the scatter of these scaling relationships and find no significant improvement over the lag-only approach for any of the three lines. Using these relations, we further estimate the masses and accretion rates of 246 AGN with single epoch methods. We also use these relations to estimate the relative sizes of the H$\beta$, MgII and CIV emitting regions, and find evidence that the MgII emission may occur further out than H$\beta$. In short, we provide a comprehensive benchmark of high redshift AGN reverberation mapping at the close of this most recent generation of surveys, including light curves, time-delays, and a set of significantly improved radius-luminosity relations for use with high-redshift populations.

The primordial helium mass fraction, $Y_{p}$, is a key observational pillar of Big Bang nucleosynthesis and a sensitive probe of early-Universe physics. Over the past several decades, numerous observational $Y_{p}$ determinations have been published using a wide range of astrophysical tracers and cosmological techniques. Although recent measurements exhibit striking convergence and increasingly small uncertainties, the statistical and historical context of this consensus has not been examined systematically. Here, we compile and analyse a comprehensive dataset of observational $Y_{p}$ determinations published between the late-1960s and 2022. The final sample comprises 143 reported values spanning multiple tracers. We find clear evidence for long-term convergence in published $Y_{p}$ values, punctuated by statistically significant change points in the mid-2000s and early 2010s. Careful examination reveals that many extragalactic Hii-region determinations are not fully independent, relying on re-analyses or partial reuse of a limited number of observational datasets. This reduces the effective number of independent constraints and provides important context for interpreting the precision of recent results. Our findings do not challenge the overall consistency of modern $Y_{p}$ determinations with standard cosmology, but they underscore the importance of accounting for data dependence, methodological homogeneity and historical evolution when synthesising measurements.

Techniques using multiple tracers of the large scale structure of the universe show great promise for examining the fundamentals of our Universe’s cosmology. Such techniques rely on the different relationship between the overdensity of tracers and the broader matter overdensity, enabling cosmic-variance-free tests of primordial non-Gaussianity in the initial curvature perturbations. There is a great opportunity for current and future all-sky extra-galactic radio surveys to make use of this technique to test for non-Gaussianity at a precision greater than existing all-sky constraints from the cosmic microwave background. To realise this goal, there is a need for accurate simulations. Previous radio galaxy simulations have either been realistic but covering only a small area (and so unhelpful for cosmological forecasts), or all-sky dark matter only cosmological simulations but having no connection to a real radio galaxy population. In this study, we use the FLAMINGO suite of cosmological surveys, as well as the matching of dark matter halos to radio galaxy population, to create an accurate sky simulation in order to examine the feasibility of multi-tracer techniques. We present an analysis of the clustering (with a bias model for the simulation), as well as redshift distributions, source counts and radio luminosity functions, and discuss future work on non-Gaussianity detection.

Comparing the redshifts of fast radio burst (FRB) host galaxies to FRB dispersion measures has unlocked a new probe of the cosmological distribution of ionised gas. However, the necessary optical observations to identify FRB hosts and measure their redshifts are becoming increasingly onerous. Here, we analyse the ability of the Legacy Survey of Space and Time (LSST), being conducted by the Vera C. Rubin Observatory, to identify FRB host galaxies, and the utility of LSST photometric redshifts for FRB cosmology. By combining a model of FRB host galaxy r-band magnitudes, $m_r$, with predictions for the FRB z–DM distribution, we create a method to predict the $m_r(z)$ distribution for the host galaxies of FRBs detected by radio surveys. We then predict these distributions for the coherent modes of the Australian Square Kilometre Array Pathfinder (ASKAP) and MeerKAT. We find that even a single visit with Rubin will be able to identify 65% of FRB host galaxies detected by ASKAP’s coherent upgrade, ‘CRACO’; while the final 10 year co-added images will identify 81% of those from MeerKAT’s tied array beams. We find that estimated photo-z errors result in a decreased precision of only 7% on $H_0$ for ASKAP’s CRACO system. The impact of missing faint FRB hosts, however, will degrade sensitivity to $H_0$ by 47% or 62% when combined with photo-z errors. All told, Rubin’s LSST will be an incredibly powerful survey for facilitating FRB cosmology, although supplemental observations may be useful for particularly faint and distant host galaxies.

Polarimetric observations at sub-GHz frequencies offer unique access to the magnetised universe through Faraday rotation and depolarization studies, but achieving reliable polarization calibration at these frequencies remains challenging. We report the identification and resolution of a systematic polarization calibration instability in the upgraded Giant Metrewave Radio Telescope (uGMRT) Band 4 (550–750 MHz). Through diagnostic observations of multiple calibrators, we discovered that the cross-hand phase response varies with the fractional polarization of the observed source, violating the fundamental assumption of calibration transferability in radio interferometry. Systematic engineering tests traced this behaviour to the quadrature hybrid (QH) polariser in the frontend signal chain. We conducted a controlled experiment in which the QH was bypassed in seven antennas, converting them to linear polarization feeds. The bypassed system shows dramatically improved performance: instrumental leakage reduced from 10–15% to 2–5%, residual leakage after calibration reduced from $\sim$0.5% to less than $0.2\%$, and stable cross-hand phases independent of source polarization. For the polarised source DA 240 (RM $=$ 3.3 rad m$^{-2}$), the QH-bypassed system accurately recovers the expected $25^\circ$ polarization angle rotation across the band, which the with QH system fails to reproduce. These results establish that the QH polariser is the dominant source of polarimetric instability in uGMRT Band 4 and demonstrate that its removal enables reliable sub-GHz polarimetry. We recommend the linear feed configuration for science cases requiring accurate polarization angle and rotation measure measurements.

We present the second data release (DR2) of Spectra and Polarisation in Cutouts of Extragalactic sources from RACS (SPICE-RACS). SPICERACS DR2 is derived from the third low-band epoch of the Rapid ASKAP Continuum Survey (RACS-low3) and covers the entire sky from the South celestial pole up to a declination of +49°; approximately 87.5% of the celestial sphere. We produce ‘cutout’ spectral cubes in Stokes I, Q, U around 4 million radio sources and extract spectra towards 5 million radio components. Across our observed band of 799.5 to 1087.5MHz we find an rms noise of ∼ 200 μJy/PSF, an angular resolution of ∼ 15″, and residual wide-field instrumental polarisation on the order of 0.1 %. After de-duplication, our polarisation catalogue contains the detection of 2.5×105 (3.4×105) Faraday rotation measures (RM) for components with a linearly polarised signal above 8σ (6σ). This places SPICE-RACS DR2 as the largest single RM catalogue ever produced by nearly an order of magnitude; the number of RMs in our catalogue alone is ∼ 5 times larger than every previous RM catalogue combined. Our resulting RM grid has an areal density of , providing an effective ‘resolution’ of ∼ 23′, and reveals striking features across the sky. The broad-band RMs have a median uncertainty of ∼2 radm−2, and include complexity metrics and information from the time domain. The breadth and quality of the SPICE-RACS DR2 dataset will enable a new generation of RM science. Further, SPICE-RACS will provide an ideal reference for forthcoming deep polarisation surveys such as the ASKAP POSSUM survey. All of our data products are publicly available on the CSIRO Data Access Portal (DAP) and the CSIRO ASKAP Science Data Archive (CASDA).