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 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.

Over the last decade, the continuing decline in the cost of digital computing technology has brought about a dramatic transformation in how digital instrumentation for radio astronomy is developed and operated. In most cases, it is now possible to interface consumer computing hardware, e.g. inexpensive graphics processing units and storage devices, directly to the raw data streams produced by radio telescopes. Such systems bring with them myriad benefits: straightforward upgrade paths, cost savings through leveraging an economy of scale, and a lowered barrier to entry for scientists and engineers seeking to add new instrument capabilities. Additionally, the typical data-interconnect technology used with general-purpose computing hardware -- Ethernet -- naturally permits multiple subscribers to a single raw data stream. This allows multiple science programs to be conducted in parallel. When combined with broad bandwidths and wide primary fields of view, radio telescopes become capable of achieving many science goals simultaneously. Moreover, because many science programs are not strongly dependent on observing cadence and direction (e.g. searches for extraterrestrial intelligence and radio transient surveys), these so-called "commensal" observing programs can dramatically increase the scientific productivity and discovery potential of an observatory. In this whitepaper, we detail a project to add an Ethernet-based commensal observing mode to the Jansky Very Large Array (VLA), and discuss how this mode could be leveraged to conduct a powerful program to constrain the distribution of advanced life in the universe through a search for radio emission indicative of technology. We also discuss other potential science use-cases for the system, and how the system could be used for technology development towards next-generation processing systems for the Next Generation VLA.

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.

Here, we report on the detection and verification of Fast Radio Burst FRB 180301, which occurred on UTC 2018 March 1 during the Breakthrough Listen observations with the Parkes telescope. Full-polarization voltage data of the detection were captured--a first for non-repeating FRBs--allowing for coherent de-dispersion and additional verification tests. The coherently de-dispersed dynamic spectrum of FRB 180301 shows complex, polarized frequency structure over a small fractional bandwidth. As FRB 180301 was detected close to the geosynchronous satellite band during a time of known 1-2 GHz satellite transmissions, we consider whether the burst was due to radio interference emitted or reflected from an orbiting object. Based on the preponderance of our verification tests, we find that FRB 180301 is likely of astrophysical origin, but caution that anthropogenic sources cannot conclusively be ruled out.

Boyajian' s Star (KIC 8462852) has received significant attention due to its unusual periodic brightness fluctuations detected by the Kepler Spacecraft and subsequent ground based observations. Possible explanations for the dips in the photometric measurements include interstellar or circumstellar dust, and it has been speculated that an artificial megastructure could be responsible. We analyze 177 high-resolution spectra of Boyajian's Star in an effort to detect potential laser signals from extraterrestrial civilizations. The spectra were obtained by the Lick Observatory's Automated Planet Finder telescope as part of the Breakthrough Listen Project, and cover the wavelength range of visible light from 374 to 970 nm. We calculate that the APF would be capable of detecting lasers of power greater than approximately 24 MW at the distance of Boyajian's Star, d = 1470 ly. The top candidates from the analysis can all be explained as either cosmic ray hits, stellar emission lines or atmospheric air glow emission lines.

We report on high-resolution spectra obtained by the Automated Planet Finder and high resolution optical Levy Spectrometer and the search for periodic spectral modulations, such as those reported in Borra (2016). In their analysis of 2.5 million spectra from the Sloan Digital Sky Survey, Borra (2016) report periodic spectral modulations in 234 stars, and suggest that these signals may be evidence of extra-terrestrial civilizations. To further evaluate this claim, we observed a total of three of the 234 stars with the Automated Planet Finder Telescope and Levy Spectrometer including all stars brighter than a visual magnitude of 14. Fourier analysis of the resultant spectra of these three sources does not reveal any periodic spectral modulations at the reported period, nor at any other period.

We introduce Probing Radio Intensity at high-Z from Marion (PRIZM), a new experiment designed to measure the globally averaged sky brightness, including the expected redshifted 21 cm neutral hydrogen absorption feature arising from the formation of the first stars. PRIZM consists of two dual-polarization antennas operating at central frequencies of 70 and 100 MHz, and the experiment is located on Marion Island in the sub-Antarctic. We describe the initial design and configuration of the PRIZM instrument that was installed in 2017, and we present preliminary data that demonstrate that Marion Island offers an exceptionally clean observing environment, with essentially no visible contamination within the FM band.

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.

We report on a search for engineered signals from a sample of 692 nearby stars using the Robert C. Byrd Green Bank Telescope, undertaken as part of the $Breakthrough~Listen~Initiative$ search for extraterrestrial intelligence. Observations were made over 1.1$-$1.9 GHz (L band), with three sets of five-minute observations of the 692 primary targets, interspersed with five-minute observations of secondary targets. By comparing the "ON" and "OFF" observations we are able to identify terrestrial interference and place limits on the presence of engineered signals from putative extraterrestrial civilizations inhabiting the environs of the target stars. During the analysis, eleven events passed our thresholding algorithm, but a detailed analysis of their properties indicates they are consistent with known examples of anthropogenic radio frequency interference. We conclude that, at the time of our observations, none of the observed systems host high-duty-cycle radio transmitters emitting between 1.1 and 1.9 GHz with an Equivalent Isotropic Radiated Power of $\sim10^{13}$ W, which is readily achievable by our own civilization. Our results suggest that fewer than $\sim$ 0.1$\%$ of the stellar systems within 50 pc possess the type of transmitters searched in this survey.

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.

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.

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.

We describe in general terms the practical use in astronomy of a higher-order statistical quantity called Spectral Kurtosis (SK), and describe the first implementation of SK-enabled firmware in the F-engine (Fourier transform-engine) of a digital FX correlator for Expanded Owens Valley Solar Array (EOVSA). The development of the theory for SK is summarized, leading to an expression for generalized SK that is applicable to both SK spectrometers and those not specifically designed for SK. We also give the means for computing both the SK estimator and thresholds for its application as a discriminator of RFI contamination. Tests of the performance of EOVSA as an SK spectrometer are shown to agree precisely with theoretical expectations, and the methods for configuring the correlator for correct SK operation are described.

We present the target selection for the Breakthrough Listen search for extraterrestrial intelligence during the first year of observations at the Green Bank Telescope, Parkes Telescope and Automated Planet Finder. On the way to observing 1,000,000 nearby stars in search of technological signals, we present three main sets of objects we plan to observe in addition to a smaller sample of exotica. We choose the 60 nearest stars, all within 5.1 pc from the sun. Such nearby stars offer the potential to observe faint radio signals from transmitters having a power similar to those on Earth. We add a list of 1649 stars drawn from the Hipparcos catalog that span the Hertzprung-Russell diagram, including all spectral types along the main sequence, subgiants, and giant stars. This sample offers diversity and inclusion of all stellar types, but with thoughtful limits and due attention to main sequence stars. Our targets also include 123 nearby galaxies composed of a "morphological-type-complete" sample of the nearest spirals, ellipticals, dwarf spherioidals, and irregulars. While their great distances hamper the detection of technological electromagnetic radiation, galaxies offer the opportunity to observe billions of stars simultaneously and to sample the bright end of the technological luminosity function. We will also use the Green Bank and Parkes telescopes to survey the plane and central bulge of the Milky Way. Finally, the complete target list includes several classes of exotica, including white dwarfs, brown dwarfs, black holes, neutron stars, and asteroids in our Solar System.

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.

We present a new method for interferometric imaging that is ideal for the large fields of view and compact arrays common in 21 cm cosmology. We first demonstrate the method with simulations for two very different low frequency interferometers, the Murchison Widefield Array (MWA) and the MIT Epoch of Reionization (MITEoR) Experiment. We then apply the method to the MITEoR data set collected in July 2013 to obtain the first northern sky map from 128 MHz to 175 MHz at about 2 degree resolution, and find an overall spectral index of -2.73+/-0.11. The success of this imaging method bodes well for upcoming compact redundant low-frequency arrays such as HERA. Both the MITEoR interferometric data and the 150 MHz sky map are publicly available at http://space.mit.edu/home/tegmark/omniscope.html.

The Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) has been working for a decade to reduce the time and cost of designing, building and deploying new digital radio-astronomy instruments. Today, CASPER open-source technology powers over 45 scientific instruments worldwide, and is used by scientists and engineers at dozens of academic institutions. In this paper we catalog the current offerings of the CASPER collaboration, and instruments past and present built by CASPER users and developers. We describe the ongoing state of software development, as CASPER looks to support a broader range of programming environments and hardware and ensure compatibility with the latest vendor tools.

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.

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.

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.

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.