We demonstrate that the Bayesian evidence can be used to find a good approximation of the ground truth likelihood function of a dataset, a goal of the likelihood-free inference (LFI) paradigm. As a concrete example, we use forward modelled sky-averaged 21-cm signal antenna temperature datasets where we artificially inject noise structures of various physically motivated forms. We find that the Gaussian likelihood performs poorly when the noise distribution deviates from the Gaussian case, for example, heteroscedastic radiometric or heavy-tailed noise. For these non-Gaussian noise structures, we show that the generalised normal likelihood is on a similar Bayesian evidence scale with comparable sky-averaged 21-cm signal recovery as the ground truth likelihood function of our injected noise. We therefore propose the generalised normal likelihood function as a good approximation of the true likelihood function if the noise structure is a priori unknown.

We present observations of the four $^2 \Pi _{3/2}\,J=3/2$ ground-rotational state transitions of the hydroxyl molecule (OH) along 107 lines of sight both in and out of the Galactic plane: 92 sets of observations from the Arecibo telescope and 15 sets of observations from the Australia Telescope Compact Array (ATCA). Our Arecibo observations included off-source pointings, allowing us to measure excitation temperature ($T_{\rm ex}$) and optical depth, while our ATCA observations give optical depth only. We perform Gaussian decomposition using the Automated Molecular Excitation Bayesian line-fitting Algorithm ‘Amoeba’ (Petzler, Dawson, & Wardle 2021, ApJ, 923, 261) fitting all four transitions simultaneously with shared centroid velocity and width. We identify 109 features across 38 sightlines (including 58 detections along 27 sightlines with excitation temperature measurements). While the main lines at 1665 and 1667 MHz tend to have similar excitation temperatures (median $|\Delta T_{\rm ex}({\rm main})|=0.6\,$K, 84% show $|\Delta T_{\rm ex}({\rm main})|<2\,$K), large differences in the 1612 and 1720 MHz satellite line excitation temperatures show that the gas is generally not in LTE. For a selection of sightlines, we compare our OH features to associated (on-sky and in velocity) Hi cold gas components (CNM) identified by Nguyen et al. (2019, ApJ, 880, 141) and find no strong correlations. We speculate that this may indicate an effective decoupling of the molecular gas from the CNM once it accumulates.

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

We describe the first results from the All-sky BRIght, Complete Quasar Survey (AllBRICQS), which aims to discover the last remaining optically bright quasars. We present 156 spectroscopically confirmed quasars (140 newly identified) having $|b|>10^{\circ}$. 152 of the quasars have Gaia DR3 magnitudes brighter than $B_{P}=16.5$ or $R_{P}=16$ mag, while four are slightly fainter. The quasars span a redshift range of $z=0.07-3.93$. In particular, we highlight the properties of J0529-4351 at $z=3.93$, which, if unlensed, is one of the most intrinsically luminous quasars in the Universe. The AllBRICQS sources have been selected by combining data from the Gaia and WISE all-sky satellite missions, and we successfully identify quasars not flagged as candidates by Gaia Data Release 3. We expect the completeness to be $\approx$96% within our magnitude and latitude limits, while the preliminary results indicate a selection purity of $\approx$96%. The optical spectroscopy used for source classification will also enable detailed quasar characterisation, including black hole mass measurements and identification of foreground absorption systems. The AllBRICQS sources will greatly enhance the number of quasars available for high-signal-to-noise follow-up with present and future facilities.

Active galactic nuclei (AGN) have been observed as far as redshift $z \sim 7$. They are crucial in investigating the early Universe as well as the growth of supermassive black holes at their centres. Radio-loud AGN with their jets seen at a small viewing angle are called blazars and show relativistic boosting of their emission. Thus, their apparently brighter jets are easier to detect in the high-redshift Universe. DES J014132.4–542749.9 is a radio-luminous but X-ray weak blazar candidate at $z = 5$. We conducted high-resolution radio interferometric observations of this source with the Australian Long Baseline Array at $1.7$ and $8.5$ GHz. A single, compact radio-emitting feature was detected at both frequencies with a flat radio spectrum. We derived the milliarcsecond-level accurate position of the object. The frequency dependence of its brightness temperature is similar to that of blazar sources observed at lower redshifts. Based on our observations, we can confirm its blazar nature. We compared its radio properties with those of two other similarly X-ray-weak and radio-bright AGN, and found that they show very different relativistic boosting characteristics.

Multi-messenger observations of the transient sky to detect cosmic explosions and counterparts of gravitational wave mergers critically rely on orbiting wide-FoV telescopes to cover the wide range of wavelengths where atmospheric absorption and emission limit the use of ground facilities. Thanks to continuing technological improvements, miniaturised space instruments operating as distributed-aperture constellations are offering new capabilities for the study of high-energy transients to complement ageing existing satellites. In this paper we characterise the performance of the upcoming joint SpIRIT and HERMES-TP/SP constellation for the localisation of high-energy transients through triangulation of signal arrival times. SpIRIT is an Australian technology and science demonstrator satellite designed to operate in a low-Earth Sun-synchronous Polar orbit that will augment the science operations for the equatorial HERMES-TP/SP constellation. In this work we simulate the improvement to the localisation capabilities of the HERMES-TP/SP constellation when SpIRIT is included in an orbital plane nearly perpendicular (inclination = 97.6°) to the HERMES-TP/SP orbits. For the fraction of GRBs detected by three of the HERMES satellites plus SpIRIT, we find that the combined constellation is capable of localising 60% of long GRBs to within ${\sim}30\,\textrm{deg}^{2}$ on the sky, and 60% of short GRBs within ${\sim}1850\,\textrm{deg}^{2}$ ($1\sigma$ confidence regions), though it is beyond the scope of this work to characterise or rule out systematic uncertainty of the same order of magnitude. Based purely on statistical GRB localisation capabilities (i.e., excluding systematic uncertainties and sky coverage), these figures for long GRBs are comparable to those reported by the Fermi Gamma Burst Monitor instrument. These localisation statistics represents a reduction of the uncertainty for the burst localisation region for both long and short GRBs by a factor of ${\sim}5$ compared to the HERMES-TP/SP alone. Further improvements by an additional factor of 2 (or 4) can be achieved by launching an additional 4 (or 6) SpIRIT-like satellites into a Polar orbit, respectively, which would both increase the fraction of sky covered by multiple satellite elements, and also enable localisation of ${\geq} 60\%$ of long GRBs to within a radius of ${\sim}1.5^{\circ}$ (statistical uncertainty) on the sky, clearly demonstrating the value of a distributed all-sky high-energy transient monitor composed of nano-satellites.

Pulsars have been studied extensively over the last few decades and have proven instrumental in exploring a wide variety of physics. Discovering more pulsars emitting at low radio frequencies is crucial to further our understanding of spectral properties and emission mechanisms. The Murchison Widefield Array Voltage Capture System (MWA VCS) has been routinely used to study pulsars at low frequencies and discover new pulsars. The MWA VCS offers the unique opportunity of recording complex voltages from all individual antennas (tiles), which can be off-line beamformed or correlated/imaged at millisecond time resolution. Devising imaged-based methods for finding pulsar candidates, which can be verified in beamformed data, can accelerate the complete process and lead to more pulsar detections. Image-based searches for pulsar candidates can reduce the number of tied-array beams required, increasing compute resource efficiency. Despite a factor of $\sim$4 loss in sensitivity, searching for pulsar candidates in images from the MWA VCS, we can explore a larger parameter space, potentially leading to discoveries of pulsars missed by high-frequency surveys such as steep spectrum pulsars, exotic binary systems, or pulsars obscured in high-time resolution time series data by propagation effects. Image-based searches are also essential to probing parts of parameter space inaccessible to traditional beamformed searches with the MWA (e.g. at high dispersion measures). In this paper we describe the innovative approach and capability of dual-processing MWA VCS data, that is forming 1-s visibilities and sky images, finding pulsar candidates in these images, and verifying by forming tied-array beam. We developed and tested image-based methods of finding pulsar candidates, which are based on pulsar properties such as steep spectral index, polarisation and variability. The efficiency of these methodologies has been verified on known pulsars, and the main limitations explained in terms of sensitivity and low-frequency spectral turnover of some pulsars. No candidates were confirmed to be a new pulsar, but this new capability will now be applied to a larger subset of observations to accelerate pulsar discoveries with the MWA and potentially speed up future searches with the SKA-Low.

We use the MaNGA integral field spectroscopic survey of low-redshift galaxies to compare the stellar populations of the bulge and disc components, identified from their Sérsic profiles, for various samples of galaxies. Bulge-dominated regions tend to be more metal-rich and have slightly older stellar ages than their associated disc-dominated regions. The metallicity difference is consistent with the deeper gravitational potential in bulges relative to discs, which allows bulges to retain more of the metals produced by stars. The age difference is due to star formation persisting longer in discs relative to bulges. Relative to galaxies with lower stellar masses, galaxies with higher stellar masses tend to have bulge-dominated regions that are more metal-rich and older (in light-weighted measurements) than their disc-dominated regions. This suggests high-mass galaxies quench from the inside out, while lower-mass galaxies quench across the whole galaxy simultaneously. Early-type galaxies tend to have bulge-dominated regions the same age as their disc-dominated regions, while late-type galaxies tend to have disc-dominated regions significantly younger than their bulge-dominated regions. Central galaxies tend to have a greater metallicity difference between their bulge-dominated regions and disc-dominated regions than satellite galaxies at similar stellar mass. This difference may be explained by central galaxies being subject to mergers or extended gas accretion bringing new, lower-metallicity gas to the disc, thereby reducing the average metallicity and age of the stars; quenching of satellite discs may also play a role.

The amount and complexity of data delivered by modern galaxy surveys has been steadily increasing over the past years. New facilities will soon provide imaging and spectra of hundreds of millions of galaxies. Extracting coherent scientific information from these large and multi-modal data sets remains an open issue for the community and data-driven approaches such as deep learning have rapidly emerged as a potentially powerful solution to some long lasting challenges. This enthusiasm is reflected in an unprecedented exponential growth of publications using neural networks, which have gone from a handful of works in 2015 to an average of one paper per week in 2021 in the area of galaxy surveys. Half a decade after the first published work in astronomy mentioning deep learning, and shortly before new big data sets such as Euclid and LSST start becoming available, we believe it is timely to review what has been the real impact of this new technology in the field and its potential to solve key challenges raised by the size and complexity of the new datasets. The purpose of this review is thus two-fold. We first aim at summarising, in a common document, the main applications of deep learning for galaxy surveys that have emerged so far. We then extract the major achievements and lessons learned and highlight key open questions and limitations, which in our opinion, will require particular attention in the coming years. Overall, state-of-the-art deep learning methods are rapidly adopted by the astronomical community, reflecting a democratisation of these methods. This review shows that the majority of works using deep learning up to date are oriented to computer vision tasks (e.g. classification, segmentation). This is also the domain of application where deep learning has brought the most important breakthroughs so far. However, we also report that the applications are becoming more diverse and deep learning is used for estimating galaxy properties, identifying outliers or constraining the cosmological model. Most of these works remain at the exploratory level though which could partially explain the limited impact in terms of citations. Some common challenges will most likely need to be addressed before moving to the next phase of massive deployment of deep learning in the processing of future surveys; for example, uncertainty quantification, interpretability, data labelling and domain shift issues from training with simulations, which constitutes a common practice in astronomy.

We present a catalogue of over 7000 sources from the GLEAM survey which have significant structure on sub-arcsecond scales at 162 MHz. The compact nature of these sources was detected and quantified via their Interplanetary Scintillation (IPS) signature, measured in interferometric images from the Murchison Widefield Array. The advantage of this approach is that all sufficiently compact sources across the survey area are included down to a well-defined flux density limit. The survey is based on ${\sim}250\times 10\hbox{-}\mathrm{min}$ observations, and the area covered is somewhat irregular, but the area within $1\,\mathrm{h}<\mathrm{RA}<11\,\mathrm{h}$ ; $-10^\circ<\mathrm{Decl.}<+20^\circ$ is covered entirely, and over 85% of this area has a detection limit for compact structure below 0.2 Jy. 7839 sources clearly showing IPS were detected ( ${>}5\sigma$ confidence), with a further 5550 tentative ( ${>}2\sigma$ confidence) detections. Normalised Scintillation Indices (NSI; a measure of the fraction of flux density coming from a compact component) are reported for these sources. Robust and informative upper limits on the NSI are reported for a further 31081 sources. This represents the largest survey of compact sources at radio frequencies ever undertaken.

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

The three-body problem is famously chaotic, with no closed-form analytical solutions. However, hierarchical systems of three or more bodies can be stable over indefinite timescales. A system is considered hierarchical if the bodies can be divided into separate two-body orbits with distinct time and length scales, such that one orbit is only mildly affected by the gravitation of the other bodies. Previous work has mapped the stability of such systems at varying resolutions over a limited range of parameters, and attempts have been made to derive analytic and semi-analytic stability boundary fits to explain the observed phenomena. Certain regimes are understood relatively well. However, there are large regions of the parameter space which remain unmapped, and for which the stability boundary is poorly understood. We present a comprehensive numerical study of the stability boundary of hierarchical triples over a range of initial parameters. Specifically, we consider the mass ratio of the inner binary to the outer third body ( $q_\mathrm{out}$ ), mutual inclination (i), initial mean anomaly and eccentricity of both the inner and outer binaries ( $e_\mathrm{in}$ and $e_\mathrm{out}$ respectively). We fit the dependence of the stability boundary on $q_\mathrm{ out}$ as a threshold on the ratio of the inner binary’s semi-major axis to the outer binary’s pericentre separation $a_\mathrm{in}/R_\mathrm{p, out} \leq 10^{-0.6 + 0.04q_\mathrm{out}}q_\mathrm{out}^{0.32+0.1q_\mathrm{out}}$ for coplanar prograde systems. We develop an additional factor to account for mutual inclination. The resulting fit predicts the stability of $10^4$ orbits randomly initialised close to the stability boundary with $87.7\%$ accuracy.

We present WALLABY pilot data release 1, the first public release of H i pilot survey data from the Wide-field ASKAP L-band Legacy All-sky Blind Survey (WALLABY) on the Australian Square Kilometre Array Pathfinder. Phase 1 of the WALLABY pilot survey targeted three $60\,\mathrm{deg}^{2}$ regions on the sky in the direction of the Hydra and Norma galaxy clusters and the NGC 4636 galaxy group, covering the redshift range of $z \lesssim 0.08$ . The source catalogue, images and spectra of nearly 600 extragalactic H i detections and kinematic models for 109 spatially resolved galaxies are available. As the pilot survey targeted regions containing nearby group and cluster environments, the median redshift of the sample of $z \approx 0.014$ is relatively low compared to the full WALLABY survey. The median galaxy H i mass is $2.3 \times 10^{9}\,{\rm M}_{{\odot}}$ . The target noise level of $1.6\,\mathrm{mJy}$ per 30′′ beam and $18.5\,\mathrm{kHz}$ channel translates into a $5 \sigma$ H i mass sensitivity for point sources of about $5.2 \times 10^{8} \, (D_{\rm L} / \mathrm{100\,Mpc})^{2} \, {\rm M}_{{\odot}}$ across 50 spectral channels ( ${\approx} 200\,\mathrm{km \, s}^{-1}$ ) and a $5 \sigma$ H i column density sensitivity of about $8.6 \times 10^{19} \, (1 + z)^{4}\,\mathrm{cm}^{-2}$ across 5 channels ( ${\approx} 20\,\mathrm{km \, s}^{-1}$ ) for emission filling the 30′′ beam. As expected for a pilot survey, several technical issues and artefacts are still affecting the data quality. Most notably, there are systematic flux errors of up to several 10% caused by uncertainties about the exact size and shape of each of the primary beams as well as the presence of sidelobes due to the finite deconvolution threshold. In addition, artefacts such as residual continuum emission and bandpass ripples have affected some of the data. The pilot survey has been highly successful in uncovering such technical problems, most of which are expected to be addressed and rectified before the start of the full WALLABY survey.

We present the Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) Pilot Phase I Hi kinematic models. This first data release consists of Hi observations of three fields in the direction of the Hydra and Norma clusters, and the NGC 4636 galaxy group. In this paper, we describe how we generate and publicly release flat-disk tilted-ring kinematic models for 109/592 unique Hi detections in these fields. The modelling method adopted here—which we call the WALLABY Kinematic Analysis Proto-Pipeline (WKAPP) and for which the corresponding scripts are also publicly available—consists of combining results from the homogeneous application of the FAT and 3DBarolo algorithms to the subset of 209 detections with sufficient resolution and $S/N$ in order to generate optimised model parameters and uncertainties. The 109 models presented here tend to be gas rich detections resolved by at least 3–4 synthesised beams across their major axes, but there is no obvious environmental bias in the modelling. The data release described here is the first step towards the derivation of similar products for thousands of spatially resolved WALLABY detections via a dedicated kinematic pipeline. Such a large publicly available and homogeneously analysed dataset will be a powerful legacy product that that will enable a wide range of scientific studies.

The formation of the first stars and galaxies during ‘Cosmic Dawn’ is thought to have imparted a faint signal onto the 21-cm spin temperature from atomic Hydrogen gas in the early Universe. Observationally, an absorption feature should be measurable as a frequency dependence in the sky-averaged (i.e. global) temperature at meter wavelengths. This signal should be separable from the smooth—but orders of magnitude brighter—foregrounds by jointly fitting a log-polynomial and absorption trough to radiometer spectra. A majority of approaches to measure the global 21-cm signal use radiometer systems on dipole-like antennas. Here, we argue that beamforming-based methods may allow radio arrays to measure the global 21-cm signal. We simulate an end-to-end drift-scan observation of the radio sky at 50–100 MHz using a zenith-phased array, and find that the complex sidelobe structure introduces a significant frequency-dependent systematic. However, the $\lambda/D$ evolution of the beam width with frequency does not confound detection. We conclude that a beamformed array with a median sidelobe level ${\sim}-50$ dB may offer an alternative method to measure the global 21-cm signal. This level is achievable by arrays with $O(10^5)$ antennas.

We present the pulsar_spectra software repository, an open-source pulsar flux density catalogue and automated spectral fitting software that finds the best spectral model and produces publication-quality plots. The python-based software includes features that enable users in the astronomical community to add newly published spectral measurements to the catalogue as they become available. The spectral fitting software is an implementation of the method described in Jankowski et al. (2018, MNRAS, 473, 4436) which uses robust statistical methods to decide on the best-fitting model for individual pulsar spectra. pulsar_spectra is motivated by the need for a centralised repository for pulsar flux density measurements to make published measurements more accessible to the astronomical community and provide a suite of tools for measuring spectra.

The shape of emission lines in the optical spectra of star-forming galaxies reveals the kinematics of the diffuse gaseous component. We analyse the shape of prominent emission lines in a sample of $\sim$ 53000 star-forming galaxies from the Sloan Digital Sky Survey, focusing on departures from gaussianity. Departures from a single gaussian profile allow us to probe the motion of gas and to assess the role of outflows. The sample is divided into groups according to their stellar velocity dispersion and star formation rate (SFR). The spectra within each group are stacked to improve the signal-to-noise ratio of the emission lines, to remove individual signatures, and to enhance the effect of SFR on the shapes of the emission lines. The moments of the emission lines, including kurtosis and skewness, are determined. We find that most of the emission lines in strong star-forming systems unequivocally feature negative kurtosis. This signature is present in $\mathrm{H}\unicode{x03B2}$ , $\mathrm{H}\unicode{x03B1}$ , [N ii], and [S ii] in massive galaxies with high SFRs. We attribute it as evidence of radial outflows of ionised gas driven by the star formation of the galaxies. Also, most of the emission lines in low-mass systems with high SFRs feature negative skewness, and we interpret it as evidence of dust obscuration in the galactic disk. These signatures are however absent in the [O iii] line, which is believed to trace a different gas component. The observed trend is significantly stronger in face-on galaxies, indicating that star formation drives the outflows along the galactic rotation axis, presumably the path of least resistance. The data suggest that outflows driven by star formation exert accumulated impacts on the interstellar medium, and the outflow signature is more evident in older galaxies as they have experienced a longer total duration of star formation.

We demonstrate the effectiveness of a Bayesian evidence -based analysis for diagnosing and disentangling the sky-averaged 21-cm signal from instrumental systematic effects. As a case study, we consider a simulated REACH pipeline with an injected systematic. We demonstrate that very poor performance or erroneous signal recovery is achieved if the systematic remains unmodelled. These effects include sky-averaged 21-cm posterior estimates resembling a very deep or wide signal. However, when including parameterised models of the systematic, the signal recovery is dramatically improved in performance. Most importantly, a Bayesian evidence-based model comparison is capable of determining whether or not such a systematic model is needed as the true underlying generative model of an experimental dataset is in principle unknown. We, therefore, advocate a pipeline capable of testing a variety of potential systematic errors with the Bayesian evidence acting as the mechanism for detecting their presence.

The redshifted cosmological 21-cm signal emitted by neutral hydrogen during the first billion years of the universe is much fainter relative to other galactic and extragalactic radio emissions, posing a great challenge towards detection of the signal. Therefore, precise instrumental calibration is a vital prerequisite for the success of radio interferometers such as the Murchison Widefield Array (MWA), which aim for a 21-cm detection. Over the previous years, novel calibration techniques targeting the power spectrum paradigm of EoR science have been actively researched and where possible implemented. Some of these improvements, for the MWA, include the accuracy of sky models used in calibration and the treatment of ionospheric effects, both of which introduce unwanted contamination to the EoR window. Despite sophisticated non-traditional calibration algorithms being continuously developed over the years to incorporate these methods, the large datasets needed for EoR measurements require high computational costs, leading to trade-offs that impede making use of these new tools to maximum benefit. Using recently acquired computation resources for the MWA, we test the full capabilities of the state-of-the-art calibration techniques available for the MWA EoR project, with a focus on both direction-dependent and direction-independent calibration. Specifically, we investigate improvements that can be made in the vital calibration stages of sky modelling, ionospheric correction, and compact source foreground subtraction as applied in the hybrid foreground mitigation approach (one that combines both foreground subtraction and avoidance). Additionally, we investigate a method of ionospheric correction using interpolated ionospheric phase screens and assess its performance in the power spectrum space. Overall, we identify a refined RTS calibration configuration that leads to an at least 2 factor reduction of the EoR window power contamination at the $0.1 \; \textrm{hMpc}^{-1}$ scale. The improvement marks a step further towards detecting the 21-cm signal using the MWA and the forthcoming SKA low telescope.

We examine the redshifts of a comprehensive set of published Type Ia supernovae, and provide a combined, improved catalogue with updated redshifts. We improve on the original catalogues by using the most up-to-date heliocentric redshift data available; ensuring all redshifts have uncertainty estimates; using the exact formulae to convert heliocentric redshifts into the Cosmic Microwave Background (CMB) frame; and utilising an improved peculiar velocity model that calculates local motions in redshift-space and more realistically accounts for the external bulk flow at high-redshifts. We review 2607 supernova redshifts; 2285 are from unique supernovae and 322 are from repeat-observations of the same supernova. In total, we updated 990 unique heliocentric redshifts, and found 5 cases of missing or incorrect heliocentric corrections, 44 incorrect or missing supernova coordinates, 230 missing heliocentric or CMB frame redshifts, and 1200 missing redshift uncertainties. The absolute corrections range between $10^{-8} \leq \Delta z \leq 0.038$ , and RMS $(\Delta z) \sim 3{\times 10^{-3}}$ . The sign of the correction was essentially random, so the mean and median corrections are small: $4{\times 10^{-4}}$ and $4{\times 10^{-6}}$ respectively. We examine the impact of these improvements for $H_0$ and the dark energy equation of state w and find that the cosmological results change by $\Delta H_0 = -0.12\,\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$ and $\Delta w = 0.003$ , both significantly smaller than previously reported uncertainties for $H_0$ of 1.0 $\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$ and w of 0.04 respectively.