Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration.

Radio Frequency Interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array grow larger in number of receivers. To address this, we present a Deep Fully Convolutional Neural Network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known "ground truth" dataset for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6$\times 10^{5}$ HERA time-ordered 1024 channeled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time-frequency context which increases discrimination between RFI and Non-RFI. The inclusion of phase when predicting achieves a Recall of 0.81, Precision of 0.58, and $F_{2}$ score of 0.75 as applied to our HERA-67 observations.

New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.

Low-frequency radio continuum observations of edge-on galaxies are ideal to study cosmic-ray electrons (CREs) in halos via radio synchrotron emission and to measure magnetic field strengths. We obtained new observations of the edge-on spiral galaxy NGC 891 at 129-163 MHz with the LOw Frequency ARray (LOFAR) and at 13-18 GHz with the Arcminute Microkelvin Imager (AMI) and combine them with recent high-resolution Very Large Array (VLA) observations at 1-2 GHz, enabling us to study the radio continuum emission over two orders of magnitude in frequency. The spectrum of the integrated nonthermal flux density can be fitted by a power law with a spectral steepening towards higher frequencies or by a curved polynomial. Spectral flattening at low frequencies due to free-free absorption is detected in star-forming regions of the disk. The mean magnetic field strength in the halo is 7 +- 2 $\mu$G. The scale heights of the nonthermal halo emission at 146 MHz are larger than those at 1.5 GHz everywhere, with a mean ratio of 1.7 +- 0.3, indicating that spectral ageing of CREs is important and that diffusive propagation dominates. The halo scale heights at 146 MHz decrease with increasing magnetic field strengths which is a signature of dominating synchrotron losses of CREs. On the other hand, the spectral index between 146 MHz and 1.5 GHz linearly steepens from the disk to the halo, indicating that advection rather than diffusion is the dominating CRE transport process. This issue calls for refined modelling of CRE propagation.

The Arcminute Microkelvin Imager (AMI) telescopes located at the Mullard Radio Astronomy Observatory near Cambridge have been significantly enhanced by the implementation of a new digital correlator with 1.2 MHz spectral resolution. This system has replaced a 750-MHz resolution analogue lag-based correlator, and was designed to mitigate the effects of radio frequency interference, particularly from geostationary satellites that contaminate observations at low declinations. The upgraded instrument consists of 18 ROACH2 Field Programmable Gate Array platforms used to implement a pair of real-time FX correlators -- one for each of AMI's two arrays. The new system separates the down-converted RF baseband signal from each AMI receiver into two 2.3 GHz-wide sub-bands which are each digitized at 5-Gsps with 8 bits of precision. These digital data streams are filtered into 2048 frequency channels and cross-correlated using FPGA hardware, with a commercial 10 Gb Ethernet switch providing high-speed data interconnect. Images formed using data from the new digital correlator show over an order of magnitude improvement in dynamic range over the previous system. The ability to observe at low declinations has also been significantly improved.

QUIJOTE (Q-U-I JOint TEnerife) is an experiment designed to achieve CMB B-mode polarization detection and sensitive enough to detect a primordial gravitational-wave component if the B-mode amplitude is larger than r = 0.05. It consists in two telescopes and three instruments observing in the frequency range 10-42 GHz installed at the Teide Observatory in the Canary Islands, Spain. The observing strategy includes three raster scan deep integration fields for cosmology, a nominal wide survey covering the Northen Sky and specific raster scan deep integration observations in regions of specific interest. The main goals of the project are presented and the first scientific results obtained with the first instrument are reviewed.

We report on late time radio and X-ray observations of the tidal disruption event candidate ASASSN-14li, covering the first 1000 days of the decay phase. For the first $\sim200$ days the radio and X-ray emission fade in concert. This phase is better fit by an exponential decay at X-ray wavelengths, while the radio emission is well described by either an exponential or the canonical $t^{-5/3}$ decay assumed for tidal disruption events. The correlation between radio and X-ray emission during this period can be fit as $L_{R}\propto L_{X}^{1.9\pm0.2}$. After 400 days the radio emission at $15.5\,\textrm{GHz}$ has reached a plateau level of $244\pm8\,\mu\textrm{Jy}$ which it maintains for at least the next 600 days, while the X-ray emission continues to fade exponentially. This steady level of radio emission is likely due to relic radio lobes from the weak AGN-like activity implied by historical radio observations. We note that while most existing models are based upon the evolution of ejecta which are decoupled from the central black hole, the radio : X-ray correlation during the declining phase is also consistent with core jet emission coupled to a radiatively efficient accretion flow.

We present a deep survey of the SuperCLASS super-cluster - a region of sky known to contain five Abell clusters at redshift $z\sim0.2$ - performed using the Arcminute Microkelvin Imager (AMI) Large Array (LA) at 15.5$~$GHz. Our survey covers an area of approximately 0.9 square degrees. We achieve a nominal sensitivity of $32.0~\mu$Jy beam$^{-1}$ toward the field centre, finding 80 sources above a $5\sigma$ threshold. We derive the radio colour-colour distribution for sources common to three surveys that cover the field and identify three sources with strongly curved spectra - a high-frequency-peaked source and two GHz-peaked-spectrum sources. The differential source count (i) agrees well with previous deep radio source count, (ii) exhibits no evidence of an emerging population of star-forming galaxies, down to a limit of 0.24$~$mJy, and (iii) disagrees with some models of the 15$~$GHz source population. However, our source count is in agreement with recent work that provides an analytical correction to the source count from the SKADS Simulated Sky, supporting the suggestion that this discrepancy is caused by an abundance of flat-spectrum galaxy cores as-yet not included in source population models.

The candidate black hole X-ray binary Swift J1753.5-0127 faded to quiescence in 2016 November, after a prolonged outburst that was discovered in 2005. Nearly three months later the system displayed renewed activity that lasted through 2017 July. Here, we present radio and X-ray monitoring over ~3 months of the renewed activity to study the coupling between the jet and the inner regions of the disk/jet system. Our observations cover low X-ray luminosities that have not historically been well-sampled (Lx~2e33 - 1e36 erg/s; 1-10 keV), including time periods when the system was both brightening and fading. At these low luminosities Swift J1753.5-0127 occupies a parameter space in the radio/X-ray luminosity plane that is comparable to "canonical" systems (e.g., GX 339-4), regardless of whether the system was brightening or fading, even though during its >11-year outburst Swift J1753.5-0127 emitted less radio emission from its jet than expected. We discuss implications for the existence of a single radio/X-ray luminosity correlation for black hole X-ray binaries at the lowest luminosities (Lx < 1e35 erg/s), and we compare to supermassive black holes. Our campaign includes the lowest luminosity quasi-simultaneous radio/X-ray detection to date for a black hole X-ray binary during its rise out of quiescence, thanks to early notification from optical monitoring combined with fast responses from sensitive multiwavelength facilities.

The Arcminute Microkelvin Imager (AMI) carried out a blind survey for galaxy clusters via their Sunyaev-Zel'dovich effect decrements between 2008 and 2011. The first detection, known as AMI-CL J0300+2613, has been reobserved with AMI equipped with a new digital correlator with high dynamic range. The combination of the new AMI data and more recent high-resolution sub-mm and infra-red maps now shows the feature in fact to be a ring of positive dust-correlated Galactic emission, which is likely to be anomalous microwave emission (AME). If so, this is the first completely blind detection of AME at arcminute scales.

In this paper we present the electromagnetic modeling and beam pattern measurements of a 16-elements ultra wideband sparse random test array for the low frequency instrument of the Square Kilometer Array telescope. We discuss the importance of a small array test platform for the development of technologies and techniques towards the final telescope, highlighting the most relevant aspects of its design. We also describe the electromagnetic simulations and modeling work as well as the embedded-element and array pattern measurements using an Unmanned Aerial Vehicle system. The latter are helpful both for the validation of the models and the design as well as for the future instrumental calibration of the telescope thanks to the stable, accurate and strong radio frequency signal transmitted by the UAV. At this stage of the design, these measurements have shown a general agreement between experimental results and numerical data and have revealed the localized effect of un-calibrated cable lengths in the inner side-lobes of the array pattern.

We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna.

The connection between accretion and jet production in accreting white dwarf binary systems, especially dwarf novae, is not well understood. Radio wavelengths provide key insights into the mechanisms responsible for accelerating electrons, including jets and outflows. Here we present densely-sampled radio coverage, obtained with the Arcminute MicroKelvin Imager Large Array, of the dwarf nova SS Cyg during its February 2016 anomalous outburst. The outburst displayed a slower rise (3 days per mag) in the optical than typical ones, and lasted for more than 3 weeks. Rapid radio flaring on timescales <1 hour was seen throughout the outburst. The most intriguing behavior in the radio was towards the end of the outburst where a fast, luminous ("giant"), flare peaking at ~20 mJy and lasting for 15 minutes was observed. This is the first time that such a flare has been observed in SS Cyg, and insufficient coverage could explain its non-detection in previous outbursts. These data, together with past radio observations, are consistent with synchrotron emission from plasma ejection events as being the origin of the radio flares. However, the production of the giant flare during the declining accretion rate phase remains unexplained within the standard accretion-jet framework and appears to be markedly different to similar patterns of behavior in X-ray binaries.

After more than 26 years in quiescence, the black hole transient V404 Cyg went into a luminous outburst in June 2015, and additional activity was detected in late December of the same year. Here, we present an optical spectroscopic follow-up of the December mini-outburst, together with X-ray, optical and radio monitoring that spanned more than a month. Strong flares with gradually increasing intensity are detected in the three spectral ranges during the ~10 days following the Swift trigger. Our optical spectra reveal the presence of a fast outflowing wind, as implied by the detection of a P-Cyg profile (He I - 5876 A) with a terminal velocity of ~2500 km/s. Nebular-like spectra -- with an H_alpha equivalent width of ~500 A -- are also observed. All these features are similar to those seen during the main June 2015 outburst. Thus, the fast optical wind simultaneous with the radio jet is most likely present in every V404 Cyg outburst. Finally, we report on the detection of a strong radio flare in late January 2016, when X-ray and optical monitoring had stopped due to Sun constraints.

The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to measure 21 cm emission from the primordial intergalactic medium (IGM) throughout cosmic reionization ($z=6-12$), and to explore earlier epochs of our Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is designed to characterize the evolution of the 21 cm power spectrum to constrain the timing and morphology of reionization, the properties of the first galaxies, the evolution of large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m parabolic dishes observing from 50 to 250 MHz. Currently, 19 dishes have been deployed on site and the next 18 are under construction. HERA has been designated as an SKA Precursor instrument. In this paper, we summarize HERA's scientific context and provide forecasts for its key science results. After reviewing the current state of the art in foreground mitigation, we use the delay-spectrum technique to motivate high-level performance requirements for the HERA instrument. Next, we present the HERA instrument design, along with the subsystem specifications that ensure that HERA meets its performance requirements. Finally, we summarize the schedule and status of the project. We conclude by suggesting that, given the realities of foreground contamination, current-generation 21 cm instruments are approaching their sensitivity limits. HERA is designed to bring both the sensitivity and the precision to deliver its primary science on the basis of proven foreground filtering techniques, while developing new subtraction techniques to unlock new capabilities. The result will be a major step toward realizing the widely recognized scientific potential of 21 cm cosmology.

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 very demanding requirements of the SKA-low instrument call for a challenging antenna design capable of delivering excellence performance in radiation patterns, impedance matching, polarization purity, cost, longevity, etc. This paper is devoted to the development (design and test of first prototypes) of an active ultra-wideband antenna element for the low-frequency instrument of the SKA radio telescope. The antenna element and differential low noise amplifier described here were originally designed to cover the former SKA-low band (70-450MHz) but it is now aimed to cover the re-defined SKA-low band (50-350MHz) and furthermore the antenna is capable of performing up to 650MHz with the current design. The design is focused on maximum sensitivity in a wide field of view (+/- 45deg from zenith) and low cross-polarization ratios. Furthermore, the size and cost of the element has to be kept to a minimum as millions of these antennas will need to be deployed for the full SKA in very compact configurations. The primary focus of this paper is therefore to discuss various design implications for the SKA-low telescope.

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.

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating.

QUIJOTE (Q-U-I JOint TEnerife) is a new polarimeter aimed to characterize the polarization of the Cosmic Microwave Background and other Galactic and extragalactic signals at medium and large angular scales in the frequency range 10-40 GHz. The multi-frequency (10-20~GHz) instrument, mounted on the first QUIJOTE telescope, saw first light on November 2012 from the Teide Observatory (2400~m a.s.l). During 2014 the second telescope has been installed at this observatory. A second instrument at 30~GHz will be ready for commissioning at this telescope during summer 2015, and a third additional instrument at 40~GHz is now being developed. These instruments will have nominal sensitivities to detect the B-mode polarization due to the primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

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.

Anomalous microwave emission (AME) has been observed in numerous sky regions, in the frequency range ~10-60 GHz. One of the most scrutinized regions is G159.6-18.5, located within the Perseus molecular complex. In this paper we present further observations of this region (194 hours in total over ~250 deg^2), both in intensity and in polarization. They span four frequency channels between 10 and 20 GHz, and were gathered with QUIJOTE, a new CMB experiment with the goal of measuring the polarization of the CMB and Galactic foregrounds. When combined with other publicly-available intensity data, we achieve the most precise spectrum of the AME measured to date, with 13 independent data points being dominated by this emission. The four QUIJOTE data points provide the first independent confirmation of the downturn of the AME spectrum at low frequencies, initially unveiled by the COSMOSOMAS experiment in this region. We accomplish an accurate fit of these data using models based on electric dipole emission from spinning dust grains, and also fit some of the parameters on which these models depend. We also present polarization maps with an angular resolution of ~1 deg and a sensitivity of ~25 muK/beam. From these maps, which are consistent with zero polarization, we obtain upper limits of Pi<6.3% and <2.8% (95% C.L.) respectively at 12 and 18 GHz, a frequency range where no AME polarization observations have been reported to date. These constraints are compatible with theoretical predictions of the polarization fraction from electric dipole emission originating from spinning dust grains. At the same time, they rule out several models based on magnetic dipole emission from dust grains ordered in a single magnetic domain, which predict higher polarization levels. Future QUIJOTE data in this region may allow more stringent constraints on the polarization level of the AME.

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.