Nima Razavi-Ghods

Recent experiments in cosmology, particularly those aimed at detecting the faint, redshifted, global 21[ ]cm hydrogen line (depth [ ], [ ]), have imposed stringent new requirements on radiometer calibration. In this work, we present a framework for circuit modeling and parameter inference to strengthen these calibration pipelines. This new approach enables in situ characterization of otherwise immeasurable systematics using physically motivated models. A combination of frequentist and Bayesian techniques are employed in a pipeline that supports iterative modeling, robust parameter estimation, and detailed uncertainty quantification. The framework is applied to the REACH telescope, where the precise correction of variations in the radio signal paths arising from component aging or environmental effects is critical. Circuit models of REACH’s calibration sources are developed, with the goal of predicting source temperature corrections that are conventionally obtained from laboratory measurements. By fitting the models to measured data using a convolutional cost function, a strong agreement with RMS residuals no worse than [ ] is obtained. However, Bayesian inference reveals that the resulting temperature corrections can have uncertainties on the order of 1 to [ ], caused by reflection coefficient degeneracies, measurement noise, and errors in the models. To combat this, posteriors obtained from laboratory measurements are employed as updated priors, reducing correction uncertainties down to [ ]. Ultimately, the framework provides a means of dynamically accounting for drift in system non-idealities over time, addressing the increasing precision demands of global 21 cm radio astronomy.

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, we find at 95% confidence that Δ²(k = 0.34 h Mpc⁻¹) ≤ 457 mK² at z = 7.9 and that Δ²(k = 0.36 h Mpc⁻¹) ≤ 3496 mK² 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, 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.

The design and development process for the Square Kilometre Array (SKA) radio telescope’s Low Frequency Aperture Array component was progressed during the SKA pre-construction phase by an international consortium, with the goal of meeting requirements for a critical design review. As part of the development process a full-sized prototype SKA Low ‘station’ was deployed – the Aperture Array Verification System 1 (AAVS1). We provide a system overview and describe the commissioning results of AAVS1, which is a low frequency radio telescope with 256 dual-polarisation log-periodic dipole antennas working as a phased array. A detailed system description is provided, including an in-depth overview of relevant sub-systems, ranging from hardware, firmware, software, calibration, and control sub-systems. Early commissioning results cover initial bootstrapping, array calibration, stability testing, beam-forming, and on-sky sensitivity validation. Lessons learned are presented, along with future developments.

Detecting signals from the Cosmic Dawn using “global” experiments requires an unprecedented level of instrumental calibration, especially to remove possible systematics which arise from the interaction between the antenna and receiving systems. In this paper we describe the calibration methods which can be applied to a range of experiments including the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) which is an instrument being designed to detect such signals. We use the functions described by noise wave parameters in a fast Bayesian pipeline to determine physical calibration coefficients which can describe the systematics. This approach when used in the field and on the fly could calibrate the absolute temperature measured by the antenna down to sub-kelvin levels. In this paper we present the Bayesian framework and simulated results of our pipeline using models derived from early laboratory measurements. A major advantage of this approach is that we can properly model the uncertainties in our measurement system.

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. HERA is a drift-scan array with a 10° 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 catalogs 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 Epoch of Reionization (EoR) line-of-sight modes in the range 0.2 h ⁻¹ Mpc⁻¹ < < 0.5 h ⁻¹ Mpc⁻¹. In particular, we find evidence for nonnegligible amounts of spectral structure in the raw autocorrelations that overlaps with the EoR window and is suggestive of complex instrumental effects. Through systematic modeling on a single night of data, we find we can recover these modes in the power spectrum down to the integrated noise floor, achieving a dynamic range in the EoR window of 10⁶ in power (mK² units) with respect to the bright galactic foreground signal. Future work with deeper integrations will help determine whether these systematics can continue to be mitigated down to EoR levels. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted H i emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to k _∥) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the “wedge”), and the spectrally structured emission from H i 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. For instruments pursuing a “foreground avoidance” strategy of simply measuring in the EoR window, and not attempting to model and remove foregrounds, as is the plan for the first stage of the Hydrogen Epoch of Reionization Array (HERA), characterizing the intrinsic instrument polarization response is particularly important. 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, and discuss the achieved isolation of foreground power due to the intrinsic spectral smoothness of the foregrounds, the instrument chromaticity, and the calibration. 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 that expected from the interferometer array configuration and the modeled primary beam response. Furthermore, 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 (HERA) 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’ data set 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 × 105 HERA time-ordered 1024 channelled 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 F2 score of 0.75 as applied to our HERA-67 observations.