Search

The search for radio technosignatures is an anomaly detection problem: candidate signals represent needles of interest in the proverbial haystack of radio-frequency interference (RFI). Current search frameworks find an enormity of false-positive signals, especially in large surveys, requiring manual follow-up to a sometimes prohibitive degree. Unsupervised learning provides an algorithmic way to winnow the most anomalous signals from the chaff, as well as group together RFI signals that bear morphological similarities. We present GLOBULAR (Grouping Low-frequency Observations By Unsupervised Learning After Reduction) clustering, a signal processing method that uses HDBSCAN to reduce the false-positive rate and isolate outlier signals for further analysis. When combined with a standard narrowband signal detection and spatial filtering pipeline, such as turboSETI, GLOBULAR clustering offers significant improvements in the false-positive rate over the standard pipeline alone, suggesting dramatic potential for the amelioration of manual follow-up requirements for future large surveys. By removing RFI signals in regions of high spectral occupancy, GLOBULAR clustering may also enable the detection of signals missed by the standard pipeline. We benchmark our method against the Choza et al. (2024) turboSETI-only search of 97 nearby galaxies at L-band, demonstrating a false-positive hit reduction rate of 93.1% and a false-positive event reduction rate of 99.3%., Comment: 20 pages, 14 figures, submitted to AJ

The precise characterization and mitigation of systematic effects is one of the biggest roadblocks impeding the detection of the fluctuations of cosmological 21cm signals. Missing data in radio cosmological experiments, often due to radio frequency interference (RFI), poses a particular challenge to power spectrum analysis as it could lead to the ringing of bright foreground modes in Fourier space, heavily contaminating the cosmological signals. Here we show that the problem of missing data becomes even more arduous in the presence of systematic effects. Using a realistic numerical simulation, we demonstrate that partially flagged data combined with systematic effects can introduce significant foreground ringing. We show that such an effect can be mitigated through inpainting the missing data. We present a rigorous statistical framework that incorporates the process of inpainting missing data into a quadratic estimator of the 21cm power spectrum. Under this framework, the uncertainties associated with our inpainting method and its impact on power spectrum statistics can be understood. These results are applied to the latest Phase II observations taken by the Hydrogen Epoch of Reionization Array, forming a crucial component in power spectrum analyses as we move toward detecting 21cm signals in the ever more noisy RFI environment., Comment: 25 pages, 11 figures; Replaced to match accepted ApJ version. New version contains small editorial changes throughout in response to referee comments, no changes to results

The quest for radio signals from technologically-advanced extraterrestrial intelligence has traditionally concentrated on the vicinity of 1.4 GHz. In this paper, we extend the search to unprecedented territories, detailing our extensive observations at 6 GHz and initiating the first-ever survey at 18 GHz with the Sardinia Radio Telescope (SRT). Our strategy entailed rigorous observation sessions, totaling 36 hours, directed towards the Galactic Center and 72 selected TESS targets-making this the most comprehensive high-frequency technosignature search to date. Our narrowband signal search found no definitive evidence of drifting signals that could suggest an extraterrestrial origin from the surveyed regions. Nevertheless, our efforts have enabled us to set new constraints on the presence of radio emissions from approximately $5\times 10^{5}$ stars, establishing an isotropic radiated power limit of $1.8\times 10^{19} W$. We also provide a comparative analysis of the 'hits' recorded across both frequencies to highlight the significance of pursuing technosignature searches at higher frequencies, where the spectral landscape is less congested and more conducive to detection., Comment: Submitted to Acta Astronautica

Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correlated with the rotating sky vs. those relative to the ground, down-weighting emission in the primary beam sidelobes, and suppressing noise. FR filtering causes the noise contributions to the visibility data to become correlated in time however, making interpretation of subsequent averaging and error estimation steps more subtle. In this paper, we describe fringe rate filters that are implemented using discrete prolate spheroidal sequences, and designed for two different purposes -- beam sidelobe/horizon suppression (the `mainlobe' filter), and ground-locked systematics removal (the `notch' filter). We apply these to simulated data, and study how their properties affect visibilities and power spectra generated from the simulations. Included is an introduction to fringe-rate filtering and a demonstration of fringe-rate filters applied to simple situations to aid understanding., Comment: 21 pages, 18 figures, submitted to Monthly Notices of the Royal Astronomical Society

This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system and discuss progress on commissioning and future upgrades. As HERA is a designated Square Kilometer Array (SKA) pathfinder instrument, we also show a number of "case studies" that investigate systematics seen while commissioning the phase II system, which may be of use in the design and operation of future arrays. Common pathologies are likely to manifest in similar ways across instruments, and many of these sources of contamination can be mitigated once the source is identified.

Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability to perform high-fidelity simulations of the kinds of data that are produced by the large, many-element, radio interferometric arrays that have been purpose-built for these studies. The large scale of these arrays presents a computational challenge, as one must simulate a detailed sky and instrumental model across many hundreds of frequency channels, thousands of time samples, and tens of thousands of baselines for arrays with hundreds of antennas. In this paper, we present a fast matrix-based method for simulating radio interferometric measurements (visibilities) at the necessary scale. We achieve this through judicious use of primary beam interpolation, fast approximations for coordinate transforms, and a vectorised outer product to expand per-antenna quantities to per-baseline visibilities, coupled with standard parallelisation techniques. We validate the results of this method, implemented in the publicly-available matvis code, against a high-precision reference simulator, and explore its computational scaling on a variety of problems., Comment: 24 pages, 10 figures, accepted to RAS Techniques and Instruments, matvis is publicly available at https://github.com/HERA-Team/matvis

The Breakthrough Listen search for intelligent life is, to date, the most extensive technosignature search of nearby celestial objects. We present a radio technosignature search of the centers of 97 nearby galaxies, observed by Breakthrough Listen at the Robert C. Byrd Green Bank Telescope. We performed a narrowband Doppler drift search using the turboSETI pipeline with a minimum signal-to-noise parameter threshold of 10, across a drift rate range of $\pm$ 4 Hz\ $s^{-1}$, with a spectral resolution of 3 Hz and a time resolution of $\sim$ 18.25 s. We removed radio frequency interference by using an on-source/off-source cadence pattern of six observations and discarding signals with Doppler drift rates of 0. We assess factors affecting the sensitivity of the Breakthrough Listen data reduction and search pipeline using signal injection and recovery techniques and apply new methods for the investigation of the RFI environment. We present results in four frequency bands covering 1 -- 11 GHz, and place constraints on the presence of transmitters with equivalent isotropic radiated power on the order of $10^{26}$ W, corresponding to the theoretical power consumption of Kardashev Type II civilizations., Comment: 25 pages, 22 figures

Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method allows us to form statistical uncertainty estimates for both our models and the recovered visibilities, which is an important ingredient in establishing robust upper limits on the Epoch of Reionisation (EoR) power spectrum. In cases where the noise is large compared to the EoR signal, this approach can constrain the systematics well enough to mitigate them down to the noise level for both systematics studied. Where the noise is smaller than the EoR, our modelling can mitigate the majority of the reflections with there being only a minor level of residual systematics, while cross-coupling sees essentially complete mitigation. Our approach performs similarly to existing filtering/fitting techniques used in the HERA pipeline, but with the added benefit of rigorously propagating uncertainties. In all cases it does not significantly attenuate the underlying signal., Comment: 19 pages, 14 figures, submitted to MNRAS

The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based image power spectrum and its window functions computed from the DOM images. We use noiseless simulation, based on the Hydrogen Epoch of Reionization Array Phase I configuration, to study the image power spectrum properties. The window functions show $<10^{-11}$ of the integrated power leaks from the foreground-dominated region into the EoR window; the 2D and 1D power spectra also verify the separation between the foregrounds and the EoR., Comment: Published in ApJ

The primary goal of the search for extraterrestrial intelligence (SETI) is to gain an understanding of the prevalence of technologically advanced beings (organic or inorganic) in the Galaxy. One way to approach this is to look for technosignatures: remotely detectable indicators of technology, such as temporal or spectral electromagnetic emissions consistent with an artificial source. With the new Commensal Open-Source Multimode Interferometer Cluster (COSMIC) digital backend on the Karl G. Jansky Very Large Array (VLA), we aim to conduct a search for technosignatures that is significantly larger, more sensitive, and more efficient than previously attempted. The COSMIC system is currently operational on the VLA, recording data, and designed with the flexibility to provide user-requested modes. This paper describes the hardware system design, the current software pipeline, and plans for future development., Comment: 30 pages, 17 figures. Submitted to AJ

Radio searches for extraterrestrial intelligence have mainly targeted the discovery of narrowband continuous-wave beacons and artificially dispersed broadband bursts. Periodic pulse trains, in comparison to the above technosignature morphologies, offer an energetically efficient means of interstellar transmission. A rotating beacon at the Galactic Center (GC), in particular, would be highly advantageous for galaxy-wide communications. Here, we present blipss, a CPU-based open-source software that uses a fast folding algorithm (FFA) to uncover channel-wide periodic signals in radio dynamic spectra. Running blipss on 4.5 hours of 4-8 GHz data gathered with the Robert C. Byrd Green Bank Telescope, we searched the central 6' of our Galaxy for kHz-wide signals with periods between 11-100 s and duty cycles ($\delta$) between 10-50%. Our searches, to our knowledge, constitute the first FFA exploration for periodic alien technosignatures. We report a non-detection of channel-wide periodic signals in our data. Thus, we constrain the abundance of 4-8 GHz extraterrestrial transmitters of kHz-wide periodic pulsed signals to fewer than one in about 600,000 stars at the GC above a 7$\sigma$ equivalent isotropic radiated power of $\approx 2 \times 10^{18}$ W at $\delta \simeq 10\%$. From an astrophysics standpoint, blipss, with its utilization of a per-channel FFA, can enable the discovery of signals with exotic radio frequency sweeps departing from the standard cold plasma dispersion law., Comment: 20 pages, 11 figures, published in AJ, in press (http://seti.berkeley.edu/blipss/)

The application of Very Long Baseline Interferometry (VLBI) to the Search for Extraterrestrial Intelligence (SETI) has been limited to date, despite the technique offering many advantages over traditional single-dish SETI observations. In order to further develop interferometry for SETI, we used the European VLBI Network (EVN) at $21$~cm to observe potential secondary phase calibrators in the Kepler field. Unfortunately, no secondary calibrators were detected. However, a VLBA primary calibrator in the field, J1926+4441, offset only $\sim1.88'$ from a nearby exoplanet Kepler-111~b, was correlated with high temporal $\left(0.25 \ \rm{s}\right)$ and spectral $\left(16384 \times 488\ \rm{Hz \ channels}\right)$ resolution. During the analysis of the high-resolution data, we identified a spectral feature that was present in both the auto and cross-correlation data with a central frequency of $1420.424\pm0.0002$ MHz and a width of 0.25 MHz. We demonstrate that the feature in the cross-correlations is an artefact in the data, associated with a significant increase in each telescope's noise figure due to the presence of \ion{H}{i} in the beam. This would typically go unnoticed in data correlated with standard spectral resolution. We flag (excluded from the subsequent analysis) these channels and phase rotate the data to the location of Kepler-111~b aided by the GAIA catalogue and search for signals with $\rm{SNR}>7$. At the time of our observations, we detect no transmitters with an Equivalent Isotropically Radiated Power (EIRP) > $\sim4\times10^{15}$ W., Comment: 11 pages, 12 figures, accepted MNRAS

Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standard analysis techniques makes use of the closure phase, which allows one to bypass antenna-based direction-independent calibration. Similarly to standard approaches, we use a delay spectrum technique to search for the EoR signal. Using 94 nights of data observed with Phase I of the Hydrogen Epoch of Reionization Array (HERA), we place approximate constraints on the 21 cm power spectrum at $z=7.7$. We find at 95% confidence that the 21 cm EoR brightness temperature is $\le$(372)$^2$ "pseudo" mK$^2$ at 1.14 "pseudo" $h$ Mpc$^{-1}$, where the "pseudo" emphasises that these limits are to be interpreted as approximations to the actual distance scales and brightness temperatures. Using a fiducial EoR model, we demonstrate the feasibility of detecting the EoR with the full array. Compared to standard methods, the closure phase processing is relatively simple, thereby providing an important independent check on results derived using visibility intensities, or related., Comment: 16 pages, 14 figures, accepted for publication by MNRAS

The goal of the Search for Extraterrestrial Intelligence (SETI) is to quantify the prevalence of technological life beyond Earth via their "technosignatures". One theorized technosignature is narrowband Doppler drifting radio signals. The principal challenge in conducting SETI in the radio domain is developing a generalized technique to reject human radio frequency interference (RFI). Here, we present the most comprehensive deep-learning based technosignature search to date, returning 8 promising ETI signals of interest for re-observation as part of the Breakthrough Listen initiative. The search comprises 820 unique targets observed with the Robert C. Byrd Green Bank Telescope, totaling over 480, hr of on-sky data. We implement a novel beta-Convolutional Variational Autoencoder to identify technosignature candidates in a semi-unsupervised manner while keeping the false positive rate manageably low. This new approach presents itself as a leading solution in accelerating SETI and other transient research into the age of data-driven astronomy., Comment: 10 pages of main paper followed by 16 pages of methods; 17 figures total and 7 tables; published in Nature Astronomy

A growing avenue for determining the prevalence of life beyond Earth is to search for "technosignatures" from extraterrestrial intelligences/agents. Technosignatures require significant energy to be visible across interstellar space and thus intentional signals might be concentrated in frequency, in time, or in space, to be found in mutually obvious places. Therefore, it could be advantageous to search for technosignatures in parts of parameter space that are mutually-derivable to an observer on Earth and a distant transmitter. In this work, we used the L-band (1.1-1.9 GHz) receiver on the Robert C. Byrd Green Bank Telescope (GBT) to perform the first technosignature search pre-synchronized with exoplanet transits, covering 12 Kepler systems. We used the Breakthrough Listen turboSETI pipeline to flag narrowband hits ($\sim$3 Hz) using a maximum drift rate of $\pm$614.4 Hz/s and a signal-to-noise threshold of 5 - the pipeline returned $\sim 3.4 \times 10^5$ apparently-localized features. Visual inspection by a team of citizen scientists ruled out 99.6% of them. Further analysis found 2 signals-of-interest that warrant follow-up, but no technosignatures. If the signals-of-interest are not re-detected in future work, it will imply that the 12 targets in the search are not producing transit-aligned signals from 1.1-1.9 GHz with transmitter powers $>$60 times that of the former Arecibo radar. This search debuts a range of innovative technosignature techniques: citizen science vetting of potential signals-of-interest, a sensitivity-aware search out to extremely high drift rates, a more flexible method of analyzing on-off cadences, and an extremely low signal-to-noise threshold., Comment: 18 pages, 11 figures

Radio Frequency Interference (RFI) is one of the systematic challenges preventing 21cm interferometric instruments from detecting the Epoch of Reionization. To mitigate the effects of RFI on data analysis pipelines, numerous inpaint techniques have been developed to restore RFI corrupted data. We examine the qualitative and quantitative errors introduced into the visibilities and power spectrum due to inpainting. We perform our analysis on simulated data as well as real data from the Hydrogen Epoch of Reionization Array (HERA) Phase 1 upper limits. We also introduce a convolutional neural network that capable of inpainting RFI corrupted data in interferometric instruments. We train our network on simulated data and show that our network is capable at inpainting real data without requiring to be retrained. We find that techniques that incorporate high wavenumbers in delay space in their modeling are best suited for inpainting over narrowband RFI. We also show that with our fiducial parameters Discrete Prolate Spheroidal Sequences (DPSS) and CLEAN provide the best performance for intermittent ``narrowband'' RFI while Gaussian Progress Regression (GPR) and Least Squares Spectral Analysis (LSSA) provide the best performance for larger RFI gaps. However we caution that these qualitative conclusions are sensitive to the chosen hyperparameters of each inpainting technique. We find these results to be consistent in both simulated and real visibilities. We show that all inpainting techniques reliably reproduce foreground dominated modes in the power spectrum. Since the inpainting techniques should not be capable of reproducing noise realizations, we find that the largest errors occur in the noise dominated delay modes. We show that in the future, as the noise level of the data comes down, CLEAN and DPSS are most capable of reproducing the fine frequency structure in the visibilities of HERA data., Comment: 21 pages, 13 figures

We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $\Delta^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $\Delta^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars., Comment: 57 pages, 37 figures. Updated to match the accepted ApJ version. Corresponding author: Joshua S. Dillon

Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process. In a delay-based analysis, the mapping between instrumental and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand the power measured by an interferometer, we assess the impact of instrument characteristics and analysis choices on these window functions. Focusing on the Hydrogen Epoch of Reionization Array (HERA) as a case study, we find that long-baseline observations correspond to enhanced low-k tails of the window functions, which facilitate foreground leakage, whilst an informed choice of bandwidth and frequency taper can reduce said tails. With simple test cases and realistic simulations, we show that, apart from tracing mode mixing, the window functions help accurately reconstruct the power spectrum estimator of simulated visibilities. The window functions depend strongly on the beam chromaticity, and less on its spatial structure - a Gaussian approximation, ignoring side lobes, is sufficient. Finally, we investigate the potential of asymmetric window functions, down-weighting the contribution of low-k power to avoid foreground leakage. The window functions presented here correspond to the latest HERA upper limits for the full Phase I data. They allow an accurate reconstruction of the power spectrum measured by the instrument and will be used in future analyses to confront theoretical models and data directly in cylindrical space., Comment: 18 pages, 19 figures, accepted for publication in MNRAS

The search for extraterrestrial intelligence at radio frequencies has largely been focused on continuous-wave narrowband signals. We demonstrate that broadband pulsed beacons are energetically efficient compared to narrowband beacons over longer operational timescales. Here, we report the first extensive survey searching for such broadband pulsed beacons towards 1883 stars as a part of the Breakthrough Listen's search for advanced intelligent life. We conducted 233 hours of deep observations across 4 to 8 GHz using the Robert C. Byrd Green Bank Telescope and searched for three different classes of signals with artificial (or negative) dispersion. We report a detailed search -- leveraging a convolutional neural network classifier on high-performance GPUs -- deployed for the very first time in a large-scale search for signals from extraterrestrial intelligence. Due to the absence of any signal-of-interest from our survey, we place a constraint on the existence of broadband pulsed beacons in our solar neighborhood: $\lesssim$1 in 1000 stars have transmitter power-densities $\gtrsim$10$^5$ W/Hz repeating $\leq$500 seconds at these frequencies., Comment: Accepted for publication in ApJ

Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipeline with simulated data, we develop a maximum likelihood figure-of-merit for comparing four sky models at 166MHz with a bandwidth of 100kHz. The HERA data agree with the GLEAM catalogs to <10%. After subtracting the GLEAM point sources, the HERA data discriminate between the different continuum sky models, providing most support for the model of Byrne et al. 2021. We report the computation cost for mapping the HERA Phase I data and project the computation for the HERA 320-antenna data; both are feasible with a modern server. The algorithm is broadly applicable to other interferometers and is valid for wide-field and non-coplanar arrays., Comment: 16 pages, 10 figures, 2 tables, published on ApJ

The Galactic Center (GC), with its high density of massive stars, is a promising target for radio transient searches. In particular, the discovery and timing of a pulsar orbiting the central supermassive black hole (SMBH) of our Galaxy will enable stringent strong-field tests of gravity and accurate measurements of SMBH properties. We performed multi-epoch 4-8 GHz observations of the inner $\approx$ 15 pc of our Galaxy using the Robert C. Byrd Green Bank Telescope in 2019 August-September. Our investigations constitute the most sensitive 4-8 GHz GC pulsar survey conducted to date, reaching down to a 6.1 GHz pseudo-luminosity threshold of $\approx$ 1 mJy kpc$^2$ for a pulse duty cycle of 2.5$\%$. We searched our data in the Fourier domain for periodic signals incorporating a constant or linearly changing line-of-sight pulsar acceleration. We report the successful detection of the GC magnetar PSR J1745$-$2900 in our data. Our pulsar searches yielded a non-detection of novel periodic astrophysical emissions above a 6$\sigma$ detection threshold in harmonic-summed power spectra. We reconcile our non-detection of GC pulsars with inadequate sensitivity to a likely GC pulsar population dominated by millisecond pulsars. Alternatively, close encounters with compact objects in the dense GC environment may scatter pulsars away from the GC. The dense central interstellar medium may also favorably produce magnetars over pulsars., Comment: 20 pages, 10 figures, accepted by ApJ

The Breakthrough Listen Initiative, as part of its larger mission, is performing the most thorough technosignature search of nearby stars. Additionally, Breakthrough Listen is collaborating with scientists working on NASAs Transiting Exoplanet Survey Satellite (TESS), to examine TESS Targets of Interest (TOIs) for technosignatures. Here, we present a $1-11$ $\textrm{GHz}$ radio technosignature search of $61$ TESS TOIs that were in transit during their Breakthrough Listen observation at the Robert C. Byrd Green Bank Telescope. We performed a narrowband Doppler drift search with a minimum S/N threshold of $10$, across a drift rate range of $\pm 4$ $\textrm{Hz}$ $\textrm{s}$ $^{-1}$, with a resolution of $3$ $\textrm{Hz}$. We removed radio frequency interference by comparing signals across cadences of target sources. After interference removal, there are no remaining events in our survey, and therefore no technosignature signals-of-interest detected in this work. This null result implies that at L, S, C, and X bands, fewer than $52\textrm{%}$, $20\textrm{%}$, $16\textrm{%}$, and $15\textrm{%}$, respectively, of TESS TOIs possess a transmitter with an equivalent isotropic radiated power greater than a few times $10^{14}$ $\textrm{W}$., Comment: 12 pages; 7 figures; Accepted to The Astronomical Journal

Observation of redshifted 21-cm signal from the Epoch of Reionization (EoR) is challenging due to contamination from the bright foreground sources that exceed the signal by several orders of magnitude. The removal of this very high foreground relies on accurate calibration to keep the intrinsic property of the foreground with frequency. Commonly employed calibration techniques for these experiments are the sky model-based and the redundant baseline-based calibration approaches. However, the sky model-based and redundant baseline-based calibration methods could suffer from sky-modeling error and array redundancy imperfection issues, respectively. In this work, we introduce the hybrid correlation calibration ("CorrCal") scheme, which aims to bridge the gap between redundant and sky-based calibration by relaxing redundancy of the array and including sky information into the calibration formalisms. We demonstrate the slight improvement of power spectra, about $-6\%$ deviation at the bin right on the horizon limit of the foreground wedge-like structure, relative to the power spectra before the implementation of "CorrCal" to the data from the Precision Array for Probing the Epoch of Reionization (PAPER) experiment, which was otherwise calibrated using redundant baseline calibration. This small improvement of the foreground power spectra around the wedge limit could be suggestive of reduced spectral structure in the data after "CorrCal" calibration, which lays the foundation for future improvement of the calibration algorithm and implementation method., Comment: 18 pages, 12 figures, MNRAS in press

The detection of life beyond Earth is an ongoing scientific endeavour, with profound implications. One approach, known as the search for extraterrestrial intelligence (SETI), seeks to find engineered signals (`technosignatures') that indicate the existence technologically-capable life beyond Earth. Here, we report on the detection of a narrowband signal-of-interest at ~982 MHz, recorded during observations toward Proxima Centauri with the Parkes Murriyang radio telescope. This signal, `BLC1', has characteristics broadly consistent with hypothesized technosignatures and is one of the most compelling candidates to date. Analysis of BLC1 -- which we ultimately attribute to being an unusual but locally-generated form of interference -- is provided in a companion paper (Sheikh et al., 2021). Nevertheless, our observations of Proxima Centauri are the most sensitive search for radio technosignatures ever undertaken on a star target., Comment: 14 pages, 4 figures (+3 supplementary figures). Published open-access in Nature Astronomy

The aim of the search for extraterrestrial intelligence (SETI) is to find technologically-capable life beyond Earth through their technosignatures. On 2019 April 29, the Breakthrough Listen SETI project observed Proxima Centauri with the Parkes 'Murriyang' radio telescope. These data contained a narrowband signal with characteristics broadly consistent with a technosignature near 982 MHz ('blc1'). Here we present a procedure for the analysis of potential technosignatures, in the context of the ubiquity of human-generated radio interference, which we apply to blc1. Using this procedure, we find that blc1 is not an extraterrestrial technosignature, but rather an electronically-drifting intermodulation product of local, time-varying interferers aligned with the observing cadence. We find dozens of instances of radio interference with similar morphologies to blc1 at frequencies harmonically related to common clock oscillators. These complex intermodulation products highlight the necessity for detailed follow-up of any signal-of-interest using a procedure such as the one outlined in this work., Comment: 15 pages of main paper followed by 19 pages of supplementary material; 22 figures total; published open-access in Nature Astronomy

We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets., Comment: 31 pages, 17 figures

Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,\nu}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results., Comment: 40 pages, 19 figures, accepted to ApJ

We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $\Delta^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $\Delta^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment., Comment: Accepted to ApJ. https://reionization.org/science/public-data-release-1/

A line-of-sight towards the Galactic Center (GC) offers the largest number of potentially habitable systems of any direction in the sky. The Breakthrough Listen program is undertaking the most sensitive and deepest targeted SETI surveys towards the GC. Here, we outline our observing strategies with Robert C. Byrd Green Bank Telescope (GBT) and Parkes telescope to conduct 600 hours of deep observations across 0.7--93 GHz. We report preliminary results from our survey for ETI beacons across 1--8 GHz with 7.0 and 11.2 hours of observations with Parkes and GBT, respectively. With our narrowband drifting signal search, we were able to place meaningful constraints on ETI transmitters across 1--4 GHz and 3.9--8 GHz with EIRP limits of $\geq$4$\times$10$^{18}$ W among 60 million stars and $\geq$5$\times$10$^{17}$ W among half a million stars, respectively. For the first time, we were able to constrain the existence of artificially dispersed transient signals across 3.9--8 GHz with EIRP $\geq$1$\times$10$^{14}$ W/Hz with a repetition period $\leq$4.3 hours. We also searched our 11.2 hours of deep observations of the GC and its surrounding region for Fast Radio Burst-like magnetars with the DM up to 5000 pc cm$^{-3}$ with maximum pulse widths up to 90 ms at 6 GHz. We detected several hundred transient bursts from SGR J1745$-$2900, but did not detect any new transient burst with the peak luminosity limit across our observed band of $\geq$10$^{31}$ erg s$^{-1}$ and burst-rate of $\geq$0.23 burst-hr$^{-1}$. These limits are comparable to bright transient emission seen from other Galactic radio-loud magnetars, constraining their presence at the GC., Comment: Accepted for publication in AJ

Precision calibration poses challenges to experiments probing the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization (z~30-6). In both interferometric and global signal experiments, systematic calibration is the leading source of error. Though many aspects of calibration have been studied, the overlap between the two types of instruments has received less attention. We investigate the sky based calibration of total power measurements with a HERA dish and an EDGES style antenna to understand the role of auto-correlations in the calibration of an interferometer and the role of sky in calibrating a total power instrument. Using simulations we study various scenarios such as time variable gain, incomplete sky calibration model, and primary beam model. We find that temporal gain drifts, sky model incompleteness, and beam inaccuracies cause biases in the receiver gain amplitude and the receiver temperature estimates. In some cases, these biases mix spectral structure between beam and sky resulting in spectrally variable gain errors. Applying the calibration method to the HERA and EDGES data, we find good agreement with calibration via the more standard methods. Although instrumental gains are consistent with beam and sky errors similar in scale to those simulated, the receiver temperatures show significant deviations from expected values. While we show that it is possible to partially mitigate biases due to model inaccuracies by incorporating a time-dependent gain model in calibration, the resulting errors on calibration products are larger and more correlated. Completely addressing these biases will require more accurate sky and primary beam models., Comment: 16 pages, 13 figures, 1 table; accepted for publication in MNRAS main journal

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis., Comment: 32 pages, 20 figures. Submitted to the Astrophysical Journal

As current- and next-generation astronomical instruments come online, they will generate an unprecedented deluge of data. Analyzing these data in real time presents unique conceptual and computational challenges, and their long-term storage and archiving is scientifically essential for generating reliable, reproducible results. We present here the real-time processing (RTP) system for the Hydrogen Epoch of Reionization Array (HERA), a radio interferometer endeavoring to provide the first detection of the highly redshifted 21 cm signal from Cosmic Dawn and the Epoch of Reionization by an interferometer. The RTP system consists of analysis routines run on raw data shortly after they are acquired, such as calibration and detection of radio-frequency interference (RFI) events. RTP works closely with the Librarian, the HERA data storage and transfer manager which automatically ingests data and transfers copies to other clusters for post-processing analysis. Both the RTP system and the Librarian are public and open source software, which allows for them to be modified for use in other scientific collaborations. When fully constructed, HERA is projected to generate over 50 terabytes (TB) of data each night, and the RTP system enables the successful scientific analysis of these data., Comment: 15 pages, 3 figures, published in Astronomy and Computing

New radio telescope arrays offer unique opportunities for large-scale commensal SETI surveys. Ethernet-based architectures are allowing multiple users to access telescope data simultaneously by means of multicast Ethernet subscriptions. Breakthrough Listen will take advantage of this by conducting a commensal SETI survey on the MeerKAT radio telescope in South Africa. By subscribing to raw voltage data streams, Breakthrough Listen will be able to beamform commensally anywhere within the field of view during primary science observations. The survey will be conducted with unprecedented speed by forming and processing 64 coherent beams simultaneously, allowing the observation of several million objects within a few years. Both coherent and incoherent observing modes are planned. We present the list of desired sources for observation and explain how these sources were selected from the Gaia DR2 catalog. Given observations planned by MeerKAT's primary telescope users, we discuss their effects on the commensal survey and propose a commensal observing strategy in response. Finally, we outline our proposed approach towards observing one million nearby stars and analyse expected observing progress in the coming years.

Precise measurements of the 21 cm power spectrum are crucial for understanding the physical processes of hydrogen reionization. Currently, this probe is being pursued by low-frequency radio interferometer arrays. As these experiments come closer to making a first detection of the signal, error estimation will play an increasingly important role in setting robust measurements. Using the delay power spectrum approach, we have produced a critical examination of different ways that one can estimate error bars on the power spectrum. We do this through a synthesis of analytic work, simulations of toy models, and tests on small amounts of real data. We find that, although computed independently, the different error bar methodologies are in good agreement with each other in the noise-dominated regime of the power spectrum. For our preferred methodology, the predicted probability distribution function is consistent with the empirical noise power distributions from both simulated and real data. This diagnosis is mainly in support of the forthcoming HERA upper limit, and also is expected to be more generally applicable., Comment: 35 Pages, 9 Figures, 4 Tables. Replaced with accepted ApJ version; some clarifying text added in response to referee comments with no changes to results

This work gives an update to existing reconstructions of the Galactic Faraday rotation sky by processing almost all Faraday rotation data sets available at the end of the year 2020. Observations of extra-Galactic sources in recent years have, among other regions, further illuminated the previously under-constrained southern celestial sky, as well as parts of the inner disc of the Milky Way. This has culminated in an all-sky data set of 55,190 data points, which is a significant expansion on the 41,330 used in previous works, hence making an updated separation of the Galactic component a promising venture. The increased source density allows us to present our results in a resolution of about $1.3\cdot 10^{-2}\, \mathrm{deg}^2$ ($46.8\,\mathrm{arcmin}^2$), which is a twofold increase compared to previous works. As for previous Faraday rotation sky reconstructions, this work is based on information field theory, a Bayesian inference scheme for field-like quantities which handles noisy and incomplete data. In contrast to previous reconstructions, we find a significantly thinner and pronounced Galactic disc with small-scale structures exceeding values of several thousand $\mathrm{rad}\,\mathrm{m}^{-2}$. The improvements can mainly be attributed to the new catalog of Faraday data, but are also supported by advances in correlation structure modeling within numerical information field theory. We furthermore give a detailed discussion on statistical properties of the Faraday rotation sky and investigate correlations to other data sets., Comment: accepted in A&A; 15 pages, 12 Figures; results at https://wwwmpa.mpa-garching.mpg.de/~ensslin/research/data/faraday2020.html and http://cutouts.cirada.ca/rmcutout

Exoplanetary systems are prime targets for the Search for Extraterrestrial Intelligence (SETI). With the recent uptick in the identification of candidate and confirmed exoplanets through the work of missions like the Transiting Exoplanet Survey Satellite (TESS), we are beginning to understand that Earth-like planets are common. In this work, we extend the Breakthrough Listen (BL) search for extraterrestrial intelligence to include targeted searches of stars identified by TESS as potential exoplanet hosts. We report on 113 30-min cadence observations collected for 28 targets selected from the TESS Input Catalog (TIC) from among those identified as containing signatures of transiting planets. The targets were searched for narrowband signals from 1-11 GHz using the turboSETI pipeline architecture modified for compatibility with the Google Cloud environment. Data were searched for drift rates of +/-4 Hz/s above a minimum signal-to-noise threshold of 10, following the parameters of previous searches conducted by Price et al. (2020) and Enriquez et al. (2017). The observations presented in this work establish some of the deepest limits to date over such a wide band (1-11 GHz) for life beyond Earth. We determine that fewer than 12.72% of the observed targets possess transmitters operating at these frequencies with an Equivalent Isotropic Radiated Power (EIRP) greater than our derived threshold of 4.9*10^(14) W., Comment: Submitted to AAS Journals

We report polarization properties for eight narrowband bursts from FRB 121102 that have been re-detected in a high-frequency (4-8 GHz) Breakthrough Listen observation with the Green Bank Telescope, originally taken on 2017 August 26. The bursts were found to exhibit nearly 100% linear polarization, Faraday rotation measures (RM) bordering 9.3$\times$10$^4$ rad-m$^{-2}$, and stable polarization position angles (PA), all of which agree with burst properties previously reported for FRB 121102 at the same epoch. We confirm that these detections are indeed physical bursts with limited spectral occupancies and further support the use of sub-banded search techniques in FRB detection., Comment: 4 pages, 1 table. Accepted to RNAAS

The Event Horizon Telescope (EHT) has recently delivered the first resolved images of M87*, the supermassive black hole in the center of the M87 galaxy. These images were produced using 230 GHz observations performed in 2017 April. Additional observations are required to investigate the persistence of the primary image feature - a ring with azimuthal brightness asymmetry - and to quantify the image variability on event horizon scales. To address this need, we analyze M87* data collected with prototype EHT arrays in 2009, 2011, 2012, and 2013. While these observations do not contain enough information to produce images, they are sufficient to constrain simple geometric models. We develop a modeling approach based on the framework utilized for the 2017 EHT data analysis and validate our procedures using synthetic data. Applying the same approach to the observational data sets, we find the M87* morphology in 2009-2017 to be consistent with a persistent asymmetric ring of ~40 uas diameter. The position angle of the peak intensity varies in time. In particular, we find a significant difference between the position angle measured in 2013 and 2017. These variations are in broad agreement with predictions of a subset of general relativistic magnetohydrodynamic simulations. We show that quantifying the variability across multiple observational epochs has the potential to constrain the physical properties of the source, such as the accretion state or the black hole spin.

We have conducted a search for artificial radio emission associated with the Kepler-160 system following the report of the discovery of the Earth-like planet candidate KOI-456.04 on 2020 June 4 (arXiv:1905.09038v2). Our search targeted both narrowband (2.97 Hz) drifting ($\pm 4$ Hz s$^{-1})$ and wideband pulsed (5 ms at all bandwidths) artificially-dispersed technosignatures using the turboSETI (arXiv:1709.03491v2) and SPANDAK pipelines, respectively, from 1-8 GHz. No candidates were identified above an upper limit Equivalent Isotropic Radiated Power (EIRP) of $5.9 \times 10^{14}$ W for narrowband emission and $7.3 \times 10^{12}$ W for wideband emission. Here we briefly describe our observations and data reduction procedure.

We present Breakthrough Listen's "Exotica" Catalog as the centerpiece of our efforts to expand the diversity of targets surveyed in the Search for Extraterrestrial Intelligence (SETI). As motivation, we introduce the concept of survey breadth, the diversity of objects observed during a program. Several reasons for pursuing a broad program are given, including increasing the chance of a positive result in SETI, commensal astrophysics, and characterizing systematics. The Exotica Catalog is a 963 entry collection of 816 distinct targets intended to include "one of everything" in astronomy. It contains four samples: the Prototype sample, with an archetype of every known major type of non-transient celestial object; the Superlative sample of objects with the most extreme properties; the Anomaly sample of enigmatic targets that are in some way unexplained; and the Control sample with sources not expected to produce positive results. As far as we are aware, this is the first object list in recent times with the purpose of spanning the breadth of astrophysics. We share it with the community in hopes that it can guide treasury surveys and as a general reference work. Accompanying the catalog is extensive discussion of classification of objects and a new classification system for anomalies. Extensive notes on the objects in the catalog are available online. We discuss how we intend to proceed with observations in the catalog, contrast it with our extant Exotica efforts, and suggest similar tactics may be applied to other programs., Comment: Corresponds to version published in ApJS. 122 pages (29 pages + full appendices and references), 13 tables, 6 figures

As it stands today, the search for extraterrestrial intelligence (SETI) is highly dependent on our ability to detect interesting candidate signals, or technosignatures, in radio telescope observations and distinguish these from human radio frequency interference (RFI). Current signal search pipelines look for signals in spectrograms of intensity as a function of time and frequency (which can be thought of as images), but tend to do poorly in identifying multiple signals in a single data frame. This is especially apparent when there are dim signals in the same frame as bright, high signal-to-noise ratio (SNR) signals. In this work, we approach this problem using convolutional neural networks (CNN) as a computationally efficient method for localizing signals in synthetic observations resembling data collected by Breakthrough Listen using the Green Bank Telescope. We generate two synthetic datasets, the first with exactly one signal at various SNR levels and the second with exactly two signals, one of which represents RFI. We find that a residual CNN with strided convolutions and using multiple image normalizations as input outperforms a more basic CNN with max pooling trained on inputs with only one normalization. Training each model on a smaller subset of the training data at higher SNR levels results in a significant increase in model performance, reducing root mean square errors by at least a factor of 3 at an SNR of 25 dB. Although each model produces outliers with significant error, these results demonstrate that using CNNs to analyze signal location is promising, especially in image frames that are crowded with multiple signals., Comment: 12 pages, 10 figures, 1 table, submitted to PASP

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., Comment: 22 pages, 22 figures, accepted for publication in ApJ

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., Comment: 22 pages, 12 figures (including sub-figures). Published in PhRvD. Abstract may be slightly abridged compared to the actual manuscript due to length limitations on arXiv

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., Comment: 15 pages, 15 figures, 1 table, Accepted to MNRAS

The discovery of ubiquitous habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, has ushered in a renaissance in the search for extra-terrestrial intelligence (SETI). Large scale SETI activities are now underway at numerous international facilities. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is the largest single-aperture radio telescope in the world, well positioned to conduct sensitive searches for radio emission indicative of exo-intelligence. SETI is one of the five key science goals specified in the original FAST project plan. A collaboration with the Breakthrough Listen Initiative has been initiated in 2016 with a joint statement signed both by Dr. Jun Yan, the then director of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), and Dr. Peter Worden, the Chairman of the Breakthrough Prize Foundation. In this paper, we highlight some of the unique features of FAST that will allow for novel SETI observations. We identify and describe three different signal types indicative of a technological source, namely, narrow-band, wide-band artificially dispersed, and modulated signals. We here propose observations with FAST to achieve sensitivities never before explored., Comment: 14 pages, 8 figures, Accepted by Research in Astron. Astrophys. (FAST special issue)

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally-smooth astrophysical foregrounds. 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., Comment: 24 Pages, 19 Figures. Updated to match the accepted MNRAS version

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., Comment: Accepted to ApJ

The star EPIC 249706694 (HD 139139) was found to exhibit 28 transit-like events over an 87 day period during the Kepler mission's K2 Campaign 15 (Rappaport et al. 2019). These events did not fall into an identifiable pattern and could not be explained by a multitude of transit scenarios explored by the authors. We conduct follow-up observations at C-band frequencies with the Green Bank Telescope as part of the ongoing Breakthrough Listen search for technosignatures. We search for narrow band signals above a signal-to-noise threshold of 10 and with Doppler drift rates within +-5 Hz/s. We detect no evidence of technosignatures from EPIC 249706694 and derive an upper limit for the EIRP (Equivalent Isotropic Radiated Power) of putative transmissions to be 10 TW., Comment: 3 pages, 1 table. Accepted for publication in RNAAS

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., Comment: Submitted to ApJ

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. Applied to the Hydrogen Epoch of Reionization Array (HERA), this paper combines detailed electromagnetic and electrical simulations in order to model the chromatic effects of the instrument, and quantify its frequency and time responses. In particular, the effects of the analogue receiver, transmission cables, and mutual coupling are included. These simulations are able to accurately predict the intensity of the reflections occurring in the 150-m cable which links the antenna to the back-end. They also show that electromagnetic waves can propagate from one dish to another one through large sections of the array due to mutual coupling. 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., Comment: 13 pages, 17 figures - Submitted to MNRAS - 2nd revision