We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

Astronomers have typically set out to solve supervised machine learning problems by creating their own representations from scratch. We show that deep learning models trained to answer every Galaxy Zoo DECaLS question learn meaningful semantic representations of galaxies that are useful for new tasks on which the models were never trained. We exploit these representations to outperform several recent approaches at practical tasks crucial for investigating large galaxy samples. The first task is identifying galaxies of similar morphology to a query galaxy. Given a single galaxy assigned a free text tag by humans (e.g. "#diffuse"), we can find galaxies matching that tag for most tags. The second task is identifying the most interesting anomalies to a particular researcher. Our approach is 100% accurate at identifying the most interesting 100 anomalies (as judged by Galaxy Zoo 2 volunteers). The third task is adapting a model to solve a new task using only a small number of newly-labelled galaxies. Models fine-tuned from our representation are better able to identify ring galaxies than models fine-tuned from terrestrial images (ImageNet) or trained from scratch. We solve each task with very few new labels; either one (for the similarity search) or several hundred (for anomaly detection or fine-tuning). This challenges the longstanding view that deep supervised methods require new large labelled datasets for practical use in astronomy. To help the community benefit from our pretrained models, we release our fine-tuning code Zoobot. Zoobot is accessible to researchers with no prior experience in deep learning.

Compact symmetric objects (CSOs) are thought to be short-lived radio sources with two lobes of emission that are separated by less than a kpc in projection. However, studies of such systems at high redshift is challenging due to the limited resolution of present-day telescopes, and can be biased to the most luminous objects. Here we report imaging of a gravitationally lensed CSO at a redshift of 2.059 using very long baseline interferometry at 1.7 GHz. The data are imaged using Bayesian forward modelling deconvolution, which reveals a spectacularly extended and thin gravitational arc, and several resolved features within the lensed images. The surface brightness of the lensing-corrected source shows two mini-lobes separated by 642 pc in projection, with evidence of multiple hotspots that have brightness temperatures of 10^8.6 to 10^9.2 K, and a total luminosity density of 10^26.3 W / Hz. By combining the well-resolved radio source morphology with previous multi-wavelength studies, we conclude that this object is likely a CSO of type 2, and that the properties are consistent with the bow-shock model for compact radio sources. Our analysis highlights the importance of combining high quality data sets with sophisticated imaging and modelling algorithms for studying the high redshift Universe.

Magnetic fields in large scale structure filaments beyond galaxy clusters remain poorly understood. Superclusters offer a unique setting to study these low density environments, where weak signals make detection challenging. The Faraday rotation measure (RM) of polarized sources along supercluster lines of sight helps constrain magnetic field properties in these regions. This study aims to determine magnetic field intensity in low density environments within superclusters using RM measurements at different frequencies. We analyzed three nearby (z<0.1) superclusters, Corona Borealis, Hercules, and Leo, where polarization observations were available at 1.4 GHz and 144 MHz. Our catalogue includes 4497 polarized background sources with RM values from literature and unpublished 144 MHz data. We constructed 3D density cubes for each supercluster to estimate density at RM measurement locations. By grouping RM values into three density bins (outskirts, filaments, and nodes) we examined RM variance linked to mean density. We found an RM variance excess of 2.5 \pm 0.5 rad^2 m^{-4} between the lowest-density regions outside the supercluster and the low-density filamentary regions within. This suggests an intervening magnetic field in the supercluster filaments. Modeling the RM variance with a single scale, randomly oriented magnetic field, we constrained the line of sight magnetic field to B_{//} = 19^{+50}_{-8} nG after marginalizing over reversal scale and path length. Our findings align with previous studies of large scale structure filaments, suggesting that adiabatic compression alone (B_{||} \sim 2 nG) cannot fully explain the observed field strengths. Other amplification mechanisms likely contribute to the evolution of magnetic fields in superclusters.

Context: Interferometric imaging is algorithmically and computationally challenging as there is no unique inversion from the measurement data back to the sky maps, and the datasets can be very large. Many imaging methods already exist, but most of them focus either on the accuracy or the computational aspect. Aims: This paper aims to reduce the computational complexity of the Bayesian imaging algorithm resolve, enabling the application of Bayesian imaging for larger datasets. Methods: By combining computational shortcuts of the CLEAN algorithm with the Bayesian imaging algorithm resolve we developed an accurate and fast imaging algorithm which we name fast-resolve. Results: We validate the accuracy of the presented fast-resolve algorithm by comparing it with results from resolve on VLA Cygnus A data. Furthermore, we demonstrate the computational advantages of fast-resolve on a large MeerKAT ESO 137-006 dataset which is computationally out of reach for resolve. Conclusions: The presented algorithm is significantly faster than previous Bayesian imaging algorithms, broadening the applicability of Bayesian interferometric imaging. Specifically for the single channel VLA Cygnus A datasets fast-resolve is about $144$ times faster than resolve. For the MeerKAT dataset with multiple channels the computational speedup of fast-resolve is even larger.

This research studies the impact of high-quality training datasets on the performance of Convolutional Neural Networks (CNNs) in detecting strong gravitational lenses. We stress the importance of data diversity and representativeness, demonstrating how variations in sample populations influence CNN performance. In addition to the quality of training data, our results highlight the effectiveness of various techniques, such as data augmentation and ensemble learning, in reducing false positives while maintaining model completeness at an acceptable level. This enhances the robustness of gravitational lens detection models and advancing capabilities in this field. Our experiments, employing variations of DenseNet and EfficientNet, achieved a best false positive rate (FP rate) of $10^{-4}$, while successfully identifying over 88 per cent of genuine gravitational lenses in the test dataset. This represents an 11-fold reduction in the FP rate compared to the original training dataset. Notably, this substantial enhancement in the FP rate is accompanied by only a 2.3 per cent decrease in the number of true positive samples. Validated on the KiDS dataset, our findings offer insights applicable to ongoing missions, like Euclid.

The popularity of the CLEAN algorithm in radio interferometric imaging stems from its maturity, speed, and robustness. While many alternatives have been proposed in the literature, none have achieved mainstream adoption by astronomers working with data from interferometric arrays operating in the big data regime. This lack of adoption is largely due to increased computational complexity, absence of mature implementations, and the need for astronomers to tune obscure algorithmic parameters. This work introduces pfb-imaging: a flexible library that implements the scaffolding required to develop and accelerate general radio interferometric imaging algorithms. We demonstrate how the framework can be used to implement a sparsity-based image reconstruction technique known as (unconstrained) SARA in a way that scales with image size rather than data volume and features interpretable algorithmic parameters. The implementation is validated on terabyte-sized data from the MeerKAT telescope, using both a single compute node and Amazon Web Services computing instances.

Recently several long-period radio transients have been discovered, with strongly polarised coherent radio pulses appearing on timescales between tens to thousands of seconds [1,2]. In some cases the radio pulses have been interpreted as coming from rotating neutron stars with extremely strong magnetic fields, known as magnetars; the origin of other, occasionally periodic and less well-sampled radio transients, is still debated [3]. Coherent periodic radio emission is usually explained by rotating dipolar magnetic fields and pair production mechanisms, but such models do not easily predict radio emission from such slowly-rotating neutron stars and maintain it for extended times. On the other hand, highly magnetic isolated white dwarfs would be expected to have long spin periodicities, but periodic coherent radio emission has not yet been directly detected from these sources. Here we report observations of a long-period (21 minutes) radio transient, which we have labeled GPMJ1839-10. The pulses vary in brightness by two orders of magnitude, last between 30 and 300 seconds, and have quasi-periodic substructure. The observations prompted a search of radio archives, and we found that the source has been repeating since at least 1988. The archival data enabled constraint of the period derivative to &lt;3.6\times10^{-13}s s$^{-1}$, which is at the very limit of any classical theoretical model that predicts dipolar radio emission from an isolated neutron star.

We present a detection of correlated clustering between MeerKAT radio intensity maps and galaxies from the WiggleZ Dark Energy Survey. We find a $7.7\sigma$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($\Omega_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constraint $\Omega_{\rm HI} b_{\rm HI} r\,{=}\,[0.86\,{\pm}\,0.10\,({\rm stat})\,{\pm}\,0.12\,({\rm sys})]\,{\times}\,10^{-3}$, at an effective scale of$k_{\rm eff}\,{\sim}\,0.13\,h\,{\rm Mpc}^{-1}$. The intensity maps were obtained from a pilot survey with the MeerKAT telescope, a 64-dish pathfinder array to the SKA Observatory (SKAO). The data were collected from 10.5 hours of observations using MeerKAT's L-band receivers over six nights covering the 11hr field of WiggleZ, in the frequency range $1015-973\,{\rm MHz}$ (0.400\,{&lt;}\,z\,{&lt;}\,0.459 in redshift). This detection is the first practical demonstration of the multi-dish auto-correlation intensity mapping technique for cosmology. This marks an important milestone in the roadmap for the cosmology science case with the full SKAO.

Repeating fast radio bursts (FRBs) present excellent opportunities to identify FRB progenitors and host environments, as well as decipher the underlying emission mechanism. Detailed studies of repeating FRBs might also hold clues to the origin of FRBs as a population. We aim to detect the first two repeating FRBs: FRB 121102 (R1) and FRB 180814.J0422+73 (R2), and characterise their repeat statistics. We also want to significantly improve the sky localisation of R2. We use the Westerbork Synthesis Radio Telescope to conduct extensive follow-up of these two repeating FRBs. The new phased-array feed system, Apertif, allows covering the entire sky position uncertainty of R2 with fine spatial resolution in one pointing. We characterise the energy distribution and the clustering of detected R1 bursts. We detected 30 bursts from R1. Our measurements indicate a dispersion measure of 563.5(2) pc cm$^{-3}$, suggesting a significant increase in DM over the past few years. We place an upper limit of 8% on the linear polarisation fraction of the brightest burst. We did not detect any bursts from R2. A single power-law might not fit the R1 burst energy distribution across the full energy range or widely separated detections. Our observations provide improved constraints on the clustering of R1 bursts. Our stringent upper limits on the linear polarisation fraction imply a significant depolarisation, either intrinsic to the emission mechanism or caused by the intervening medium, at 1400 MHz that is not observed at higher frequencies. The non-detection of any bursts from R2 implies either a highly clustered nature of the bursts, a steep spectral index, or a combination of both. Alternatively, R2 has turned off completely, either permanently or for an extended period of time.

Magnetic fields play an important role in the evolution of galaxies and in shaping the dynamics of their inter-stellar medium. However, the formation history of magnetic fields from initial seed-fields to well-ordered systems is not clear. Favoured scenarios include a turbulent dynamo that amplifies the field, and a mean-field dynamo that organizes it. Such a model can be tested through observing the magnetic-field structure of galaxies in the early Universe given the relative formation time-scales involved. Here, we combine the high angular resolution of the Atacama Large Milli-metre Array (ALMA) and gravitational lensing to resolve the magnetic field structure of a 4-kpc in extent grand-design spiral when the Universe was just 2.6 Gyr old. We find that the spiral arm structure, as traced by the heated dust emission, is coincident with the linearly polarized emission, which is consistent with a highly ordered magnetic field. The time-scale needed to produce such an ordered field is likely within at least several rotations of the disk. Our study highlights the importance of combining the long baselines of ALMA and gravitational lensing to resolve the structure of galaxies at cosmologically interesting epochs.

We present the first results of the HI intensity mapping power spectrum analysis with the MeerKAT International GigaHertz Tiered Extragalactic Exploration (MIGHTEE) survey. We use data covering $\sim$4 square degrees in the COSMOS field using a frequency range 962.5 MHz to 1008.42 MHz, equivalent to HI emission in $0.4

We present MeerKAT L-band (886-1682 MHz) observations of the extended radio structure of the peculiar galaxy pair PKS 2130-538 known as the "Dancing Ghosts". The complex of bending and possibly interacting jets and lobes originate from two Active Galactic Nuclei hosts in the Abell 3785 galaxy cluster, one of which is the brightest cluster galaxy. The radio properties of the PKS 2130-538 flux density, spectral index and polarization - are typical for large, bent-tail galaxies. We also investigate a number of thin extended low surface brightness filaments originating from the lobes. Southeast from the Dancing Ghosts, we detect a region of low surface brightness emission that has no clear origin. While it could originate from the Abell 3785 radio halo, we investigate the possibility that it is a radio relic or emission associated with the two PKS 2130-538 hosts. We find no evidence of interaction between the two PKS 2130-538 hosts.

We present the description and initial results of the TRAPUM (TRAnsients And PUlsars with MeerKAT) search for pulsars associated with supernova remnants (SNRs), pulsar wind nebulae and unidentified TeV emission. The list of sources to be targeted includes a large number of well-known candidate pulsar locations but also new candidate SNRs identified using a range of criteria. Using the 64-dish MeerKAT radio telescope, we use an interferometric beamforming technique to tile the potential pulsar locations with coherent beams which we search for radio pulsations, above a signal-to-noise of 9, down to an average flux density upper limit of 30 $\mu$Jy. This limit is target-dependent due to the contribution of the sky and nebula to the system temperature. Coherent beams are arranged to overlap at their 50 per cent power radius, so the sensitivity to pulsars is not degraded by more than this amount, though realistically averages around 65 per cent if every location in the beam is considered. We report the discovery of two new pulsars; PSR J1831$-$0941 is an adolescent pulsar likely to be the plerionic engine of the candidate PWN G20.0+0.0, and PSR J1818$-$1502 appears to be an old and faint pulsar that we serendipitously discovered near the centre of a SNR already hosting a compact central object. The survey holds importance for better understanding of neutron star birth rates and the energetics of young pulsars.

We conducted the MeerKAT Vela Supercluster survey, named Vela$-$HI, to bridge the gap between the Vela SARAO MeerKAT Galactic Plane Survey (Vela$-$SMGPS, $-2^{\circ} \leq b \leq 1^{\circ}$), and optical and near-infrared spectroscopic observations of the Vela Supercluster (hereafter Vela$-$OPT/NIR) at $|b| \gtrsim 7^{\circ}$. Covering coordinates from $263^{\circ} \leq \ell \leq 284^{\circ}$ and $1^{\circ} \leq b \leq 6.2^{\circ}$ above, and $-6.7^{\circ} \leq b \leq -2^{\circ}$ below the Galactic Plane (GP), we sampled 667 fields spread across an area of ${\sim} \rm 242 ~deg^2$. With a beam size of ${\sim} 38'' \times 31''$, Vela$-$HI achieved a sensitivity of $\langle \rm rms \rangle = 0.74$ mJy beam$^{-1}$ at 44.3 km s$^{-1}$ velocity resolution over ${\sim}$67 hours of observations. We cataloged 719 galaxies, with only 211 (29%) previously documented in the literature, primarily through the HIZOA, 2MASX, and WISE databases. Among these known galaxies, only 66 had optical spectroscopic redshift information. We found marginal differences of less than one channel resolution for all galaxies in common between HIZOA and Vela$-$SMGPS, and a mean difference of $70 \pm 15$ km s$^{-1}$ between optical and HI velocities. Combining data from Vela$-$SMGPS, Vela$-$HI, and Vela$-$OPT/NIR confirmed the connection of the Hydra/Antlia filament across the GP and revealed a previously unknown diagonal wall at a heliocentric velocity range of $6500-8000$ km s$^{-1}$. Vela$-$HI reinforces the connection between the first wall at $18500-20000$ km s$^{-1}$ and the inner ZOA. The second wall seems to traverse the GP at $270^{\circ} \leq \ell \leq 279^{\circ}$, where it appears that both walls intersect, jointly covering the velocity range $18500-21500$ km s$^{-1}$.

We present the serendipitous discovery of a new radio-continuum ring-like object nicknamed Kyklos (J1802-3353), with MeerKAT UHF and L-band observations. The radio ring, which resembles the recently discovered odd radio circles (ORCs), has a diameter of 80 arcsec and is located just 6 deg from the Galactic plane. However, Kyklos exhibits an atypical thermal radio-continuum spectrum ({\alpha} = -0.1 +/- 0.3), which led us to explore different possible formation scenarios. We concluded that a circumstellar shell around an evolved massive star, possibly a Wolf-Rayet, is the most convincing explanation with the present data.

The clusters Abell 2061 and Abell 2067 in the Corona Borealis supercluster have been studied at different radio frequencies and are both known to host diffuse radio emission. The aim of this work is to investigate the radio emission in between them, suggested by low resolution observations. We analyse deep LOFAR HBA observations at 144 MHz to follow up on the possible intercluster filament suggested by previous 1.4 GHz observations. We investigate the radial profiles and the point-to-point surface brightness correlation of the emission in A2061 with radio and Xray observations, to describe the nature of the diffuse emission. We report the detection of diffuse radio emission on 800 kpc scale, more extended than previously known, reaching beyond the radio halo in A2061 towards A2067 and over the separation outside the two clusters R500 radii. We confirm the presence of a radio halo in A2061, while do not find evidence of diffuse emission in A2067. The surface brightness profile from the centre of A2061 shows an excess of emission with respect to the azimuthally averaged radio halo profile and X-ray background. We explore three different dynamical scenario to explain the nature of the diffuse emission. We analyse a trail of emission of 760 kpc between the radio halo and radio relic in A2061. This pre merger system closely resembles the two other cluster pairs where radio bridges connecting the radio halos on Mpc scales have been detected. The diffuse emission extends beyond each cluster R500 radius but in this unique case, the absence of the radio halo in A2067 is likely the reason for the observed 'gap' between the two systems. However, the point-to-point correlation results are challenging to explain. The classification of the emission remains unclear, and detailed spectral analysis and further Xray observations are required to understand the origin of the diffuse emission.

Long-period radio transients are an emerging class of extreme astrophysical events of which only three are known. These objects emit highly polarised, coherent pulses of typically a few tens of seconds duration and minutes to hour-long periods. While magnetic white dwarfs and magnetars, either isolated or in binary systems, have been invoked to explain these objects, a consensus has not emerged. Here we report on the discovery of ASKAP J193505.1+214841.0 (henceforth ASKAPJ1935+2148) with a period of 53.8 minutes exhibiting three distinct emission states - a bright pulse state with highly linearly polarised pulses with widths of 10-50 seconds; a weak pulse state which is about 26 times fainter than the bright state with highly circularly polarised pulses of widths of approximately 370 milliseconds; and a quiescent or quenched state with no pulses. The first two states have been observed to progressively evolve over the course of 8 months with the quenched state interspersed between them suggesting physical changes in the region producing the emission. A constraint on the radius of the source for the observed period rules out a magnetic white dwarf origin. Unlike other long-period sources, ASKAPJ1935+2148 is the first to exhibit drastic variations in emission modes reminiscent of neutron stars. However, its radio properties challenge our current understanding of neutron star emission and evolution.

The Small and Large Magellanic Clouds are the only galaxies outside our own in which radio pulsars have been discovered to date. The sensitivity of the MeerKAT radio interferometer offers an opportunity to search for a population of more distant extragalactic pulsars. The TRAPUM (TRansients And PUlsars with MeerKAT) collaboration has performed a radio-domain search for pulsars and transients in the dwarf star-forming galaxies Sextans A and B, situated at the edge of the local group 1.4 Mpc away. We conducted three 2-hour multi-beam observations at L-band (856-1712 MHz) with the full array of MeerKAT. No pulsars were found down to a radio pseudo-luminosity upper limit of 7.9$\pm$0.4 Jy kpc$^{2}$ at 1400 MHz, which is 28 times more sensitive than the previous limit from the Murriyang telescope. This luminosity is 30 per cent greater than that of the brightest known radio pulsar and sets a cut-off on the luminosity distributions of the entire Sextans A and B galaxies for unobscured radio pulsars beamed in our direction. A Fast Radio Burst was detected in one of the Sextans A observations at a Dispersion Measure (DM) of 737 pc cm$^{-3}$. We believe this is a background event not associated with the dwarf galaxy due to its large DM and its S/N being strongest in the wide-field incoherent beam of MeerKAT.