A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey.

We report characterization results for an engineering prototype of a next-generation low-frequency radio astronomy array. This prototype, which we refer to as the Aperture Array Verification System 0.5 (AAVS0.5), is a sparse pseudo-random array of 16 log-periodic antennas designed for 70-450 MHz. It is co-located with the Murchison Widefield Array (MWA) at the Murchison Radioastronomy Observatory (MRO) near the Australian Square Kilometre Array (SKA) core site. We characterize the AAVS0.5 using two methods: in-situ radio interferometry with astronomical sources and an engineering approach based on detailed full-wave simulation. In-situ measurement of the small prototype array is challenging due to the dominance of the Galactic noise and the relatively weaker calibration sources easily accessible in the southern sky. The MWA, with its 128 "tiles" and up to 3 km baselines, enabled in-situ measurement via radio interferometry. We present array sensitivity and beam pattern characterization results and compare to detailed full-wave simulation. We discuss areas where differences between the two methods exist and offer possibilities for improvement. Our work demonstrates the value of the dual astronomy-simulation approach in upcoming SKA design work.

The redshifted 21cm line of neutral hydrogen (Hi), potentially observable at low radio frequencies (~50-200 MHz), is a promising probe of the physical conditions of the inter-galactic medium during Cosmic Dawn and the Epoch of Reionisation (EoR). The sky-averaged Hi signal is expected to be extremely weak (~100 mK) in comparison to the Galactic foreground emission (~$10^4$ K). Moreover, the sky-averaged spectra measured by ground-based instruments are affected by chromatic propagation effects (of the order of tens of Kelvins) originating in the ionosphere. We analyze data collected with the upgraded BIGHORNS system deployed at the Murchison Radio-astronomy Observatory to assess the significance of ionospheric effects (absorption, emission and refraction) on the detection of the global EoR signal. We measure some properties of the ionosphere, such as the electron temperature ($T_e \approx$470 K at nighttime), magnitude, and variability of optical depth ($\tau_{100 MHz} \approx$0.01 and $\delta \tau \approx$0.005 at nighttime). According to the results of a statistical test applied on a large data sample, very long integrations lead to increased signal to noise even in the presence of ionospheric variability. This is further supported by the structure of the power spectrum of the sky temperature fluctuations, which has flicker noise characteristics at frequencies $\gtrsim 10^{-5}$ Hz, but becomes flat below $\approx 10^{-5}$ Hz. We conclude that the stochastic error introduced by the chromatic ionospheric effects tends to zero in an average. Therefore, the ionospheric effects and fluctuations are not fundamental impediments preventing ground-based instruments from integrating down to the precision required by global EoR experiments.

We present an in-situ antenna characterization method and results for a "low-frequency" radio astronomy engineering prototype array, characterized over the 75-300 MHz frequency range. The presence of multiple cosmic radio sources, particularly the dominant Galactic noise, makes in-situ characterization at these frequencies challenging; however, it will be shown that high quality measurement is possible via radio interferometry techniques. This method is well-known in the radio astronomy community but seems less so in antenna measurement and wireless communications communities, although the measurement challenges involving multiple undesired sources in the antenna field-of-view bear some similarities. We discuss this approach and our results with the expectation that this principle may find greater application in related fields.

The redshifted 21cm line of neutral hydrogen (HI), potentially observable at low radio frequencies (~50-200 MHz), should be a powerful probe of the physical conditions of the inter-galactic medium during Cosmic Dawn and the Epoch of Reionisation (EoR). The sky-averaged HI signal is expected to be extremely weak (~100 mK) in comparison to the foreground of up to 10000 K at the lowest frequencies of interest. The detection of such a weak signal requires an extremely stable, well characterised system and a good understanding of the foregrounds. Development of a nearly perfectly (~mK accuracy) calibrated total power radiometer system is essential for this type of experiment. We present the BIGHORNS (Broadband Instrument for Global HydrOgen ReioNisation Signal) experiment which was designed and built to detect the sky-averaged HI signal from the EoR at low radio frequencies. The BIGHORNS system is a mobile total power radiometer, which can be deployed in any remote location in order to collect radio-interference (RFI) free data. The system was deployed in remote, radio quiet locations in Western Australia and low RFI sky data have been collected. We present a description of the system, its characteristics, details of data analysis and calibration. We have identified multiple challenges to achieving the required measurement precision, which triggered two major improvements for the future system.

This paper offers an electromagnetic, more specifically array theory, perspective on understanding strong instrumental polarization effects for planar low-frequency "aperture arrays" with the Murchison Widefield Array (MWA) as an example. A long-standing issue that has been seen here is significant instrumental Stokes leakage after calibration, particularly in Stokes Q at high frequencies. A simple model that accounts for inter-element mutual coupling is presented which explains the prominence of Q leakage seen when the array is scanned away from zenith in the principal planes. On these planes, the model predicts current imbalance in the X (E-W) and Y (N-S) dipoles and hence the Q leakage. Although helpful in concept, we find that this model is inadequate to explain the full details of the observation data. This finding motivates further experimentation with more rigorous models that account for both mutual coupling and embedded element patterns. Two more rigorous models are discussed: the "full" and "average" embedded element patterns. The viability of the "full" model is demonstrated by simulating current MWA practice of using a Hertzian dipole model as a Jones matrix estimate. We find that these results replicate the observed Q leakage to approximately 2 to 5%. Finally, we offer more direct indication for the level of improvement expected from upgrading the Jones matrix estimate with more rigorous models. Using the "average" embedded pattern as an estimate for the "full" model, we find that Q leakage of a few percent is achievable.

The lowest frequency band (70 - 450 MHz) of the Square Kilometre Array will consist of sparse aperture arrays grouped into geographically-localised patches, or stations. Signals from thousands of antennas in each station will be beamformed to produce station beams which form the inputs for the central correlator. Two-stage beamforming within stations can reduce SKA-low signal processing load and costs, but has not been previously explored for the irregular station layouts now favoured in radio astronomy arrays. This paper illustrates the effects of two-stage beamforming on sidelobes and effective area, for two representative station layouts (regular and irregular gridded tile on an irregular station). The performance is compared with a single-stage, irregular station. The inner sidelobe levels do not change significantly between layouts, but the more distant sidelobes are affected by the tile layouts; regular tile creates diffuse, but regular, grating lobes. With very sparse arrays, the station effective area is similar between layouts. At lower frequencies, the regular tile significantly reduces effective area, hence sensitivity. The effective area is highest for a two-stage irregular station, but it requires a larger station extent than the other two layouts. Although there are cost benefits for stations with two-stage beamforming, we conclude that more accurate station modelling, and SKA-low configuration specifications, are required before design finalisation.