This paper deals with the calibration of the analogue chains of a Square Kilometre Array (SKA) phased aperture array station, using embedded element patterns (one per antenna in the array, thus accounting for the full effects of mutual coupling) or average element patterns to generate model visibilities. The array is composed of 256 log-periodic dipole array antennas. A simulator capable of generating such per-baseline model visibility correlation matrices was implemented, which allowed for a direct comparison of calibration results using StEFCal (Statistically Efficient and Fast Calibration) with both pattern types. Calibrating the array with StEFCal using simulator-generated model visibilities was successful and thus constitutes a possible routine for calibration of an SKA phase aperture array station. In addition, results indicate that there was no significant advantage in calibrating with embedded element patterns, with StEFCal successfully retrieving similar per-element coefficients with model visibilities generated with either pattern type. This can be of significant importance for mitigating computational costs for calibration, particularly for the consideration of real-time calibration strategies. Data from the AAVS-1 (Aperture Array Verification System 1) prototype station in Western Australia was used for demonstration purposes.

A violation of the distance-duality relation is directly linked with a temporal variation of the electromagnetic fine-structure constant. We consider a number of well-studied $f(T)$ gravity models and we revise the theoretical prediction of their corresponding induced violation of the distance-duality relationship. We further extract constraints on the involved model parameters through fine-structure constant variation data, alongside with supernovae data, and Hubble parameter measurements. Moreover, we constrain the evolution of the effective $f(T)$ gravitational constant. Finally, we compare with revised constraints on the phenomenological parametrisations of the violation of the equivalence principle in the electromagnetic sector.

The signal processing firmware that has been developed for the Low Frequency Aperture Array component of the Square Kilometre Array is described. The firmware is implemented on a dual FPGA board, that is capable of processing the streams from 16 dual polarization antennas. Data processing includes channelization of the sampled data for each antenna, correction for instrumental response and for geometric delays and formation of one or more beams by combining the aligned streams. The channelizer uses an oversampling polyphase filterbank architecture, allowing a frequency continuous processing of the input signal without discontinuities between spectral channels. Each board processes the streams from 16 antennas, as part of larger beamforming system, linked by standard Ethernet interconnections. There are envisaged to be 8192 of these signal processing platforms in the first phase of the Square Kilometre array so particular attention has been devoted to ensure the design is low cost and low power.

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

Galaxy clusters are known to host a variety of extended radio sources: tailed radio galaxies whose shape is modelled by the interaction with the intra-cluster medium (ICM); radio bubbles filling cavities in the ICM distribution and rising buoyantly through the thermal gas; diffuse giant radio sources ("halos" and "relics") revealing the presence of relativistic electrons and magnetic fields in the intra-cluster volume. It is currently the subject of an active debate how the non-thermal components that we observe at radio wavelengths affect the physical properties of the ICM and depend on the dynamical state of galaxy clusters. In this work we start our SKA1 feasibility study of the "radio cluster zoo" through simulations of a typical radio-loud cluster, hosting several bright tailed radio galaxies and a diffuse radio halo. Realistic simulations of SKA1 observations are obtained through the MeqTrees software. A new deconvolution algorithm, based on sparse representations and optimised for the detection of faint diffuse astronomical sources, is tested and compared to the classical CLEAN method.

GPU-based beamforming is a relatively unexplored area in radio astronomy, possibly due to the assumption that any such system will be severely limited by the PCIe bandwidth required to transfer data to the GPU. We have developed a CUDA-based GPU implementation of a coherent beamformer, specifically designed and optimised for deployment at the BEST-2 array which can generate an arbitrary number of synthesized beams for a wide range of parameters. It achieves $\sim$1.3 TFLOPs on an NVIDIA Tesla K20, approximately 10x faster than an optimised, multithreaded CPU implementation. This kernel has been integrated into two real-time, GPU-based time-domain software pipelines deployed at the BEST-2 array in Medicina: a standalone beamforming pipeline and a transient detection pipeline. We present performance benchmarks for the beamforming kernel as well as the transient detection pipeline with beamforming capabilities as well as results of test observation.

We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder low-frequency radio interferometer whose goal is to test technologies that improve the calibration precision and reduce the cost of the high-sensitivity 3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using massive baseline redundancy, which enables both automated precision calibration and correlator cost reduction. We demonstrate and quantify the power and robustness of redundancy for scalability and precision. We find that the calibration parameters precisely describe the effect of the instrument upon our measurements, allowing us to form a model that is consistent with $\chi^2$ per degree of freedom < 1.2 for as much as 80% of the observations. We use these results to develop an optimal estimator of calibration parameters using Wiener filtering, and explore the question of how often and how finely in frequency visibilities must be reliably measured to solve for calibration coefficients. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other next-generation instruments, which would incorporate many identical or similar technologies.

In this memo we investigate the applicability of NVIDIA Graphics Processing Units (GPUs) for SKA1-Low station and Central Signal Processing (CSP)-level processing. Station-level processing primarily involves generating a single station beam which will then be correlated with other beams in CSP. Fine channelisation can be performed either at the station of CSP-level, while coarse channelisation is assumed to be performed on FPGA-based Tile Processors, together with A/D conversion, equilisation and other processes. Rough estimates for number of GPUs required and power requirements will also be provided.

A new digital backend has been developed for the BEST-2 array at Radiotelescopi di Medicina, INAF-IRA, Italy which allows concurrent operation of an FX correlator, and a direct-imaging correlator and beamformer. This backend serves as a platform for testing some of the spatial Fourier transform concepts which have been proposed for use in computing correlations on regularly gridded arrays. While spatial Fourier transform-based beamformers have been implemented previously, this is to our knowledge, the first time a direct-imaging correlator has been deployed on a radio astronomy array. Concurrent observations with the FX and direct-imaging correlator allows for direct comparison between the two architectures. Additionally, we show the potential of the direct-imaging correlator for time-domain astronomy, by passing a subset of beams though a pulsar and transient detection pipeline. These results provide a timely verification for spatial Fourier transform-based instruments that are currently in commissioning. These instruments aim to detect highly-redshifted hydrogen from the Epoch of Reionization and/or to perform widefield surveys for time-domain studies of the radio sky. We experimentally show the direct-imaging correlator architecture to be a viable solution for correlation and beamforming.

Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral hydrogen has the potential to overtake the cosmic microwave background as our most powerful cosmological probe, because it can map a much larger volume of our Universe, shedding new light on the epoch of reionization, inflation, dark matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder low-frequency radio interferometer whose goal is to test technologies that greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR accomplishes this by using massive baseline redundancy both to enable automated precision calibration and to cut the correlator cost scaling from N^2 to NlogN, where N is the number of antennas. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious HERA project, which would incorporate many identical or similar technologies using an order of magnitude more antennas, each with dramatically larger collecting area.

Radio transient discovery using next generation radio telescopes will pose several digital signal processing and data transfer challenges, requiring specialized high-performance backends. Several accelerator technologies are being considered as prototyping platforms, including Graphics Processing Units (GPUs). In this paper we present a real-time pipeline prototype capable of processing multiple beams concurrently, performing Radio Frequency Interference (RFI) rejection through thresholding, correcting for the delay in signal arrival times across the frequency band using brute-force dedispersion, event detection and clustering, and finally candidate filtering, with the capability of persisting data buffers containing interesting signals to disk. This setup was deployed at the BEST-2 SKA pathfinder in Medicina, Italy, where several benchmarks and test observations of astrophysical transients were conducted. These tests show that on the deployed hardware eight 20MHz beams can be processed simultaneously for $\sim$640 Dispersion Measure (DM) values. Furthermore, the clustering and candidate filtering algorithms employed prove to be good candidates for online event detection techniques. The number of beams which can be processed increases proportionally to the number of servers deployed and number of GPUs, making it a viable architecture for current and future radio telescopes.

The Square Kilometre Array (SKA) is a radio telescope designed to operate between 70MHz and 10GHz. Due to this large bandwidth, the SKA will be built out of different collectors, namely antennas and dishes to cover the frequency range adequately. In order to deal with this bandwidth, innovative feeds and detectors must be designed and introduced in the initial phases of development. Moreover, the required level of resolution may only be achieved through a groundbreaking configuration of dishes and antennas. Due to the large collecting area and the specifications required for the SKA to deliver the promised science, the configuration of the dishes and the antennas within stations is an important question. This research builds on the work done before by Cohanim et al. (2004), Hassan et al. (2005) and Grigorescu et al. (2009) to further investigate the applicability of machine learning techniques to determine the optimum configurations for the collecting elements within the SKA. This work primarily uses genetic algorithms to search a large space of optimum layouts. Every genetic step provides a population with candidate individuals each of which encodes a possible solution. These are randomly generated or created through the combination of previous encodings. In this study, a number of fitness functions that rank individuals within a population of dish configurations are investigated. The UV density, connecting wire length and power spectra are considered to determine a good dish layout.

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.

The authors apologise for the preliminary submission of this paper.

The aim of this paper is to present a new fast-convergent numerically stable space-time adaptive processing (STAP) algorithm derived using a novel technique of feedback orthogonalization. The main advantages of this approach lie in its perfected stability to computational errors and faults which makes its real-time implementation on substantially faster and cheaper regular fixed-point processors possible.

Beamforming is central to the processing function of all phased arrays and becomes particularly challenging with a large number of antenna element (e.g. >100,000). The ability to beamform efficiently with reasonable power requirements is discussed in this paper. Whilst the most appropriate beamforming technology will change over time due to semiconductor and processing developments, we present a hierarchical structure which is technology agnostic and describe both Radio-Frequency (RF) and digital hierarchical beamforming approaches. We present implementations of both RF and digital beamforming systems on two antenna array demonstrators, namely the Electronic Multi Beam Radio Astronomy ConcEpt (EMBRACE) and the dualpolarisation all-digital array (2-PAD). This paper will compare and contrast both digital and analogue implementations without considering the deep system design of these arrays.

We describe an hierarchical, frequency-domain beamforming architecture for synthesising a sky beam from the wideband antenna feeds of digital aperture arrays. The development of densely-packed, all-digital aperture arrays is an important area of research required for the Square Kilometre Array (SKA) radio telescope. The design of real-time signal processing systems for digital aperture arrays is currently a central challenge in pathfinder projects worldwide. In particular, this work describes a specific implementation of the beamforming architecture to the 2-Polarisation All-Digital (2-PAD) aperture array demonstrator.

In this work, decision tree learning algorithms and fuzzy inferencing systems are applied for galaxy morphology classification. In particular, the CART, the C4.5, the Random Forest and fuzzy logic algorithms are studied and reliable classifiers are developed to distinguish between spiral galaxies, elliptical galaxies or star/unknown galactic objects. Morphology information for the training and testing datasets is obtained from the Galaxy Zoo project while the corresponding photometric and spectra parameters are downloaded from the SDSS DR7 catalogue.

Densely-packed, all-digital aperture arrays form a key area of technology development required for the Square Kilometre Array (SKA) radio telescope. The design of real-time signal processing systems for digital aperture arrays is currently a central challenge in pathfinder projects worldwide. We describe interim results of such work; an heirarchical, frequency-domain beamforming architecture for synthesising a sky beam from the wideband antenna feeds of digital aperture arrays.