This report was submitted as part of the SKA Low Frequency Aperture Array Critical Design Review describing the design of the SKA1-LOW station that took place between 2013 and 2018. The SKA1 LOW field station is inscribed in a circular area having an effective station diameter (centre to centre) of 38 meters and has 256 SKALA4 elements. This document describes the electromagnetic design of the field station. In particular it describes the layout design and the electromagnetic modelling and characteristics of the station. This document describes the effects associated with the layout and array such as mutual coupling effects, side lobe pattern and beam shape (eg. smoothness, calibration models) and presents the state of the art of our ability to measure the array performance and validate the simulation work. The current LFAA field node requirements, derived from the SKA L1 requirements, have evolved over the last years since the LFAA PDR and the System Baseline Design. The SKA1 LOW field station has been designed to meet those requirements and has therefore tracked their evolution (eg. sensitivity requirements, array diameter, etc.). The aforementioned requirements represent a very tight space with a desire for very high sensitivity over a large frequency range (7 to 1) and wide field of view (90 degrees cone around zenith) while keeping the station diameter to a minimum, so as the filling factor but at the same time allowing for sufficient space between antennas to allow for easy maintenances, amongst many others. This results in a complex design.

The Square Kilometre Array (SKA) will be both the largest radio telescope ever constructed and the largest Big Data project in the known Universe. The first phase of the project will generate on the order of 5 zettabytes of data per year. A critical task for the SKA will be its ability to process data for science, which will need to be conducted by science pipelines. Together with polarization data from the LOFAR Multifrequency Snapshot Sky Survey (MSSS), we have been developing a realistic SKA-like science pipeline that can handle the large data volumes generated by LOFAR at 150 MHz. The pipeline uses task-based parallelism to image, detect sources, and perform Faraday Tomography across the entire LOFAR sky. The project thereby provides a unique opportunity to contribute to the technological development of the SKA telescope, while simultaneously enabling cutting-edge scientific results. In this paper, we provide an update on current efforts to develop a science pipeline that can enable tight constraints on the magnetised large-scale structure of the Universe.

The Square Kilometre Array's Low Frequency instrument (SKA-LOW) will be the most sensitive aperture array ever used for radio astronomy, and will operate in the under-sampled regime for most of the frequency band where grating-lobes pose particular challenges. To achieve the expected level of sensitivity for SKA-LOW, it is particularly important to understand how interfering sources in both near and far side-lobes of the station beam affect the imaging performance. We discuss options for station designs, and adopting a random element layout, we assess its effectiveness by investigating how sources far from the main lobe of the station beam degrade images of the target field. These sources have the effect of introducing a noise-like corruption to images, which we call the Far Side-lobe Source Noise (FSSN). Using OSKAR, a GPU-accelerated software simulator, we carried out end-to-end simulations using an all-sky model and telescope configuration representative of the SKA-LOW instrument. The FSSN is a function of both the station beam and the interferometric point spread function, and decreases with increasing observation time until the coverage of the aperture plane no longer improves. Using apodisation to reduce the level of near-in side-lobes of the station beam had a noticeable improvement on the level of FSSN at low frequencies. Our results indicate that the effects of picking up sources in the side-lobes are worse at low frequencies, where the array is less sparse.

The removal of the Galactic and extragalactic foregrounds remains a major challenge for those wishing to make a detection of the Epoch of Reionization 21-cm signal. Multiple methods of modelling these foregrounds with varying levels of assumption have been trialled and shown promising recoveries on simulated data. Recently however there has been increased discussion of using the expected shape of the foregrounds in Fourier space to define an EoR window free of foreground contamination. By carrying out analysis within this window only, one can avoid the foregrounds and any statistical bias they might introduce by instead removing these foregrounds. In this paper we discuss the advantages and disadvantages of both foreground removal and foreground avoidance. We create a series of simulations with noise levels in line with both current and future experiments and compare the recovered statistical cosmological signal from foreground avoidance and a simplified, frequency independent foreground removal model. We find that while, for current generation experiments, foreground avoidance enables a better recovery at $k_{perp} > 0.6 \mathrm{Mpc}^{-1}$, foreground removal is able to recover significantly more signal at small $k_{los}$ for both current and future experiments. We also relax the assumption that the foregrounds are smooth by introducing a Gaussian random factor along the line-of-sight and then also spatially. We find that both methods perform well for foreground models with line-of-sight and spatial variations around $0.1\%$ however at levels larger than this foregrounds removal shows a greater signal recovery.

Fast Radio Bursts (FRBs), are millisecond radio signals that exhibit dispersion larger than what the Galactic electron density can account for. We have conducted a 1446 hour survey for Fast Radio Bursts (FRBs) at 145~MHz, covering a total of 4193 sq. deg on the sky. We used the UK station of the LOFAR radio telescope -- the Rawlings Array -- , accompanied for a majority of the time by the LOFAR station at Nan\c{c}ay, observing the same fields at the same frequency. Our real-time search backend, ARTEMIS, utilizes graphics processing units to search for pulses with dispersion measures up to 320 cm$^{-3}$ pc. Previous derived FRB rates from surveys around 1.4~GHz, and favoured FRB interpretations, motivated this survey, despite all previous detections occurring at higher dispersion measures. We detected no new FRBs above a signal-to-noise threshold of 10, leading to the most stringent upper limit yet on the FRB event rate at these frequencies: 29 sky$^{-1}$ day$^{-1}$ for 5~ms-duration pulses above 62~Jy. The non-detection could be due to scatter-broadening, limitations on the volume and time searched, or the shape of FRB flux density spectra. Assuming the latter and that FRBs are standard candles, the non-detection is compatible with the published FRB sky rate, if their spectra follow a power law with frequency ($\propto \nu^{\alpha}$), with $\alpha\gtrsim+0.1$, demonstrating a marked difference from pulsar spectra. Our results suggest that surveys at higher frequencies, including the low frequency component of the Square Kilometre Array, will have better chances to detect, estimate rates and understand the origin and properties of FRBs.

The exceptional sensitivity of the SKA will allow observations of the Cosmic Dawn and Epoch of Reionization (CD/EoR) in unprecedented detail, both spectrally and spatially. This wealth of information is buried under Galactic and extragalactic foregrounds, which must be removed accurately and precisely in order to reveal the cosmological signal. This problem has been addressed already for the previous generation of radio telescopes, but the application to SKA is different in many aspects. In this chapter we summarise the contributions to the field of foreground removal in the context of high redshift and high sensitivity 21-cm measurements. We use a state-of-the-art simulation of the SKA Phase 1 observations complete with cosmological signal, foregrounds and frequency-dependent instrumental effects to test both parametric and non-parametric foreground removal methods. We compare the recovered cosmological signal using several different statistics and explore one of the most exciting possibilities with the SKA --- imaging of the ionized bubbles. We find that with current methods it is possible to remove the foregrounds with great accuracy and to get impressive power spectra and images of the cosmological signal. The frequency-dependent PSF of the instrument complicates this recovery, so we resort to splitting the observation bandwidth into smaller segments, each of a common resolution. If the foregrounds are allowed a random variation from the smooth power law along the line of sight, methods exploiting the smoothness of foregrounds or a parametrization of their behaviour are challenged much more than non-parametric ones. However, we show that correction techniques can be implemented to restore the performances of parametric approaches, as long as the first-order approximation of a power law stands.

Simulations of SKA1-low were performed to estimate the noise level in images produced by the telescope over a frequency range 50-600 MHz, which extends the 50-350 MHz range of the current baseline design. The root-mean-square (RMS) deviation between images produced by an ideal, error-free SKA1-low and those produced by SKA1-low with varying levels of uncorrelated gain and phase errors was simulated. The residual in-field and sidelobe noise levels were assessed. It was found that the RMS deviations decreased as the frequency increased. The residual sidelobe noise decreased by a factor of ~5 from 50 to 100 MHz, and continued to decrease at higher frequencies, attributable to wider strong sidelobes and brighter sources at lower frequencies. The thermal noise limit is found to range between ~10 - 0.3 $\mu$Jy and is reached after ~100-100 000 hrs integration, depending on observation frequency, with the shortest integration time required at ~100 MHz.

Astrophysical radio transients are excellent probes of extreme physical processes originating from compact sources within our Galaxy and beyond. Radio frequency signals emitted from these objects provide a means to study the intervening medium through which they travel. Next generation radio telescopes are designed to explore the vast unexplored parameter space of high time resolution astronomy, but require High Performance Computing (HPC) solutions to process the enormous volumes of data that are produced by these telescopes. We have developed a combined software /hardware solution (code named ARTEMIS) for real-time searches for millisecond radio transients, which uses GPU technology to remove interstellar dispersion and detect millisecond radio bursts from astronomical sources in real-time. Here we present an introduction to ARTEMIS. We give a brief overview of the software pipeline, then focus specifically on the intricacies of performing incoherent de-dispersion. We present results from two brute-force algorithms. The first is a GPU based algorithm, designed to exploit the L1 cache of the NVIDIA Fermi GPU. Our second algorithm is CPU based and exploits the new AVX units in Intel Sandy Bridge CPUs.

The identification, and subsequent discovery, of fast radio transients through blind-search surveys requires a large amount of processing power, in worst cases scaling as $\mathcal{O}(N^3)$. For this reason, survey data are generally processed offline, using high-performance computing architectures or hardware-based designs. In recent years, graphics processing units have been extensively used for numerical analysis and scientific simulations, especially after the introduction of new high-level application programming interfaces. Here we show how GPUs can be used for fast transient discovery in real-time. We present a solution to the problem of de-dispersion, providing performance comparisons with a typical computing machine and traditional pulsar processing software. We describe the architecture of a real-time, GPU-based transient search machine. In terms of performance, our GPU solution provides a speed-up factor of between 50 and 200, depending on the parameters of the search.

We present a comprehensive multiband spectral and polarimetric study of the jet of 3C 264 (NGC 3862). Included in this study are three HST optical and ultraviolet polarimetry data sets, along with new and archival VLA radio imaging and polarimetry, a re-analysis of numerous HST broadband data sets from the near infrared to the far ultraviolet, and a Chandra ACIS-S observation. We investigate similarities and differences between optical and radio polarimetry, in both degree of polarization and projected magnetic field direction. We also examine the broadband spectral energy distribution of both the nucleus and jet of 3C 264, from the radio through the X-rays. From this we place constraints on the physics of the 3C 264 system, the jet and its dynamics. We find significant curvature of the spectrum from the near-IR to ultraviolet, and synchrotron breaks steeper than 0.5, a situation also encountered in the jet of M87. This likely indicates velocity and/or magnetic field gradients and more efficient particle acceleration localized in the faster/higher magnetic field parts of the flow. The magnetic field structure of the 3C 264 jet is remarkably smooth; however, we do find complex magnetic field structure that is correlated with changes in the optical spectrum. We find that the X-ray emission is due to the synchrotron process; we model the jet spectrum and discuss mechanisms for accelerating particles to the needed energies, together with implications for the orientation of the jet under a possible spine-sheath model.

We investigate the brightest regions of the kpc-scale jet in the powerful radio galaxy 3C 346, using new optical HST ACS/F606W polarimetry together with Chandra X-ray data and 14.9 GHz and 22.5 GHz VLA radio polarimetry. The jet shows a close correspondence in optical and radio morphology, while the X-ray emission shows an 0.80 +/- 0.17 kpc offset from the optical and radio peak positions. Optical and radio polarimetry show the same apparent magnetic field position angle and fractional polarization at the brightest knot, where the jet undergoes a large kink of almost 70 degrees in the optical and radio images. The apparent field direction here is well-aligned with the new jet direction, as predicted by earlier work that suggested the kink was the result of an oblique shock. We have explored models of the polarization from oblique shocks to understand the geometry of the 3C 346 jet, and find that the upstream flow is likely to be highly relativistic (0.91 +0.05 / -0.07 c), where the plane of the shock front is inclined at an angle of 51 (+/- 11) degrees to the upstream flow which is at an angle 14 (+8 / -7) degrees to our line of sight. The actual deflection angle of the jet in this case is only 22 degrees.

We investigate the structure of the kpc-scale jet in the nearby (z = 0.073) radio galaxy 3C 15, using new optical Hubble Space Telescope (HST) ACS/F606W polarimetry together with archival multi-band HST imaging, Chandra X-ray data and 8.4 GHz VLA radio polarimetry. The new data confirm that synchrotron radiation dominates in the optical. With matched beams, the jet is generally narrower in the optical than in the radio, suggesting a stratified flow. We examine a simple two-component model comprising a highly relativistic spine and lower-velocity sheath. This configuration is broadly consistent with polarization angle differences seen in the optical and radio data. The base of the jet is relatively brighter in the ultraviolet and X-ray than at lower energies, and the radio and optical polarization angles vary significantly as the jet brightens downstream. Further out, the X-ray intensity rises again and the apparent magnetic field becomes simpler, indicating a strong shock. Modelling the synchrotron spectrum of this brightest X-ray knot provides an estimate of its minimum internal pressure, and a comparison with the thermal pressure from X-ray emitting gas shows that the knot is overpressured and likely to be a temporary, expanding feature.

We present an overview of new HST imaging polarimetry of six nearby radio galaxies with optical jets. These observations triple the number of extragalactic jets with subarcsecond-resolution optical polarimetry. We discuss the polarization characteristics and optical morphology of each jet. We find evidence of high optical polarization, averaging 20%, but reaching upwards of $\sim 50%$ in some objects, confirming that the optical emission is synchrotron, and that the components of the magnetic fields perpendicular to the line of sight are well ordered. We find a wide range of polarization morphologies, with each jet having a somewhat different relationship between total intensity and polarized flux and the polarization position angle. We find two trends in all of these jets. First, jet ``edges'' are very often associated with high fractional optical polarizations, as also found in earlier radio observations of these and other radio jets. In these regions, the magnetic field vectors appear to track the jet direction, even at bends, where we see particularly high fractional polarizations. This indicates a strong link between the local magnetic field and jet dynamics. Second, optical flux maximum regions are usually well separated from maxima in fractional polarization and often are associated with polarization minima. This trend is not found in radio data and was found in our optical polarimetry of M87 with HST. However, unlike in M87, we do not find a general trend for near-90$^\circ$ rotations in the optical polarization vectors near flux maxima. We discuss possibilities for interpreting these trends, as well as implications for jet dynamics, magnetic field structure and particle acceleration.