-
$\texttt{MEDEA}$: A New Model for Emulating Radio Antenna Beam Patterns for 21-cm Cosmology and Antenna Design Studies
Authors:
Joshua J. Hibbard,
Bang D. Nhan,
David Rapetti,
Jack O. Burns
Abstract:
In 21-cm experimental cosmology, accurate characterization of a radio telescope's antenna beam response is essential to measure the 21-cm signal. Computational electromagnetic (CEM) simulations estimate the antenna beam pattern and frequency response by subjecting the EM model to different dependencies, or beam hyper-parameters, such as soil dielectric constant or orientation with the environment.…
▽ More
In 21-cm experimental cosmology, accurate characterization of a radio telescope's antenna beam response is essential to measure the 21-cm signal. Computational electromagnetic (CEM) simulations estimate the antenna beam pattern and frequency response by subjecting the EM model to different dependencies, or beam hyper-parameters, such as soil dielectric constant or orientation with the environment. However, it is computationally expensive to search all possible parameter spaces to optimize the antenna design or accurately represent the beam to the level required for use as a systematic model in 21-cm cosmology. We therefore present $\texttt{MEDEA}$, an emulator which rapidly and accurately generates farfield radiation patterns over a large hyper-parameter space. $\texttt{MEDEA}$ takes a subset of beams simulated by CEM software, spatially decomposes them into coefficients in a complete, linear basis, and then interpolates them to form new beams at arbitrary hyper-parameters. We test $\texttt{MEDEA}$ on an analytical dipole and two numerical beams motivated by upcoming lunar lander missions, and then employ $\texttt{MEDEA}$ as a model to fit mock radio spectrometer data to extract covariances on the input beam hyper-parameters. We find that the interpolated beams have RMS relative errors of at most $10^{-2}$ using 20 input beams or less, and that fits to mock data are able to recover the input beam hyper-parameters when the model and mock derive from the same set of beams. When a systematic bias is introduced into the mock data, extracted beam hyper-parameters exhibit bias, as expected. We propose several future extensions to $\texttt{MEDEA}$ to potentially account for such bias.
△ Less
Submitted 28 August, 2024;
originally announced August 2024.
-
Toward Spectrum Coexistence: First Demonstration of the Effectiveness of Boresight Avoidance between the NRAO Green Bank Telescope and Starlink Satellites
Authors:
Bang D. Nhan,
Christopher G. De Pree,
Matt Iverson,
Brenne Gregory,
Daniel Dueri,
Anthony Beasley,
Brian Schepis
Abstract:
NRAO and SpaceX have been engaged in coordinated testing efforts since Fall 2021, including conducting experiments on different interference avoidance schemes for the Karl G. Jansky Very Large Array (VLA) in New Mexico, and the Green Bank Telescope (GBT) inside the National Radio Quiet Zone (NRQZ) in West Virginia. The Starlink system is capable of avoiding direct illumination of telescope sites w…
▽ More
NRAO and SpaceX have been engaged in coordinated testing efforts since Fall 2021, including conducting experiments on different interference avoidance schemes for the Karl G. Jansky Very Large Array (VLA) in New Mexico, and the Green Bank Telescope (GBT) inside the National Radio Quiet Zone (NRQZ) in West Virginia. The Starlink system is capable of avoiding direct illumination of telescope sites with their adaptive tasking to place downlink beams far away. Nevertheless, even satellites operating in this mode can potentially present strong signals into the telescope's receiver system if they pass close to the telescope's main beam at the boresight. For additional protection, Starlink satellites can either momentarily redirect or completely disable their downlink channels while they pass within some minimum angular separation threshold from the telescope's boresight, methods that are referred to as "telescope boresight avoidance". In two separate experiments conducted since Fall 2023, NRAO and SpaceX arranged to have the GBT observe a fixed RA/Dec position in the sky, chosen to have a large number of close-to-boresight Starlink passages. Preliminary analysis from these two experiments illustrate the feasibility of these avoidance methods to significantly reduce, if not eliminate, the negative impact of close-to-boresight satellite passages. Importantly, these experiments demonstrate the value of continuing cooperative efforts between NRAO and SpaceX, and expanding cooperation between the radio astronomy and satellite communities more generally.
△ Less
Submitted 31 July, 2024;
originally announced July 2024.
-
The Impact of Beam Variations on Power Spectrum Estimation for 21 cm Cosmology II: Mitigation of Foreground Systematics for HERA
Authors:
Honggeun Kim,
Nicholas S. Kern,
Jacqueline N. Hewitt,
Bang D. Nhan,
Joshua S. Dillon,
Eloy de Lera Acedo,
Scott B. C. Dynes,
Nivedita Mahesh,
Nicolas Fagnoni,
David R. DeBoer
Abstract:
One key challenge in detecting 21 cm cosmological signal at z > 6 is to separate the cosmological signal from foreground emission. This can be studied in a power spectrum space where the foreground is confined to low delay modes whereas the cosmological signal can spread out to high delay modes. When there is a calibration error, however, chromaticity of gain errors propagates to the power spectru…
▽ More
One key challenge in detecting 21 cm cosmological signal at z > 6 is to separate the cosmological signal from foreground emission. This can be studied in a power spectrum space where the foreground is confined to low delay modes whereas the cosmological signal can spread out to high delay modes. When there is a calibration error, however, chromaticity of gain errors propagates to the power spectrum estimate and contaminates the modes for cosmological detection. The Hydrogen Epoch of Reionization Array (HERA) employs a high-precision calibration scheme using redundancy in measurements. In this study, we focus on the gain errors induced by nonredundancies arising from feed offset relative to the HERA's 14 meter parabolic dish element, and investigate how to mitigate the chromatic gain errors using three different methods: restricting baseline lengths for calibration, smoothing the antenna gains, and applying a temporal filter prior to calibration. With 2 cm/2 degree perturbations for translation/tilting motions, a level achievable under normal HERA operating conditions, the combination of the baseline cut and temporal filtering indicates that the spurious gain feature due to nonredundancies is significantly reduced, and the power spectrum recovers the clean foreground-free region. We found that the mitigation technique works even for large feed motions but in order to keep a stable calibration process, the feed positions need to be constrained to 2 cm for translation motions and 2 degree for tilting offset relative to the dish's vertex.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
Constraining a Model of the Radio Sky Below 6 MHz Using the Parker Solar Probe/FIELDS Instrument in Preparation for Upcoming Lunar-based Experiments
Authors:
Neil Bassett,
David Rapetti,
Bang D. Nhan,
Brent Page,
Jack O. Burns,
Marc Pulupa,
Stuart D. Bale
Abstract:
We present a Bayesian analysis of data from the FIELDS instrument on board the Parker Solar Probe (PSP) spacecraft with the aim of constraining low frequency ($\lesssim$ 6 MHz) sky in preparation for several upcoming lunar-based experiments. We utilize data recorded during PSP's ``coning roll'' maneuvers, in which the axis of the spacecraft is pointed 45$^{\circ}$ off of the Sun. The spacecraft th…
▽ More
We present a Bayesian analysis of data from the FIELDS instrument on board the Parker Solar Probe (PSP) spacecraft with the aim of constraining low frequency ($\lesssim$ 6 MHz) sky in preparation for several upcoming lunar-based experiments. We utilize data recorded during PSP's ``coning roll'' maneuvers, in which the axis of the spacecraft is pointed 45$^{\circ}$ off of the Sun. The spacecraft then rotates about a line between the Sun and the spacecraft with a period of 24 minutes. We reduce the data into two formats: roll-averaged, in which the spectra are averaged over the roll, and phase-binned, in which the spectra are binned according to the phase of the roll. We construct a forward model of the FIELDS observations that includes numerical simulations of the antenna beam, an analytic emissivity function of the galaxy, and estimates of the absorption due to free electrons. Fitting 5 parameters, we find that the roll-averaged data can be fit well by this model and we obtain posterior parameter constraints that are in general agreement with previous estimates. The model is not, however, able to fit the phase-binned data well, likely due to limitations such as the lack of non-smooth emission structure at both small and large scales, enforced symmetry between the northern and southern galactic hemispheres, and large uncertainties in the free electron density. This suggests that significant improvement in the low frequency sky model is needed in order to fully and accurately represent the sky at frequencies below 6 MHz.
△ Less
Submitted 23 January, 2023;
originally announced January 2023.
-
The Impact of Beam Variations on Power Spectrum Estimation for 21-cm Cosmology I: Simulations of Foreground Contamination for HERA
Authors:
Honggeun Kim,
Bang D. Nhan,
Jacqueline N. Hewitt,
Nicholas S. Kern,
Joshua S. Dillon,
Eloy de Lera Acedo,
Scott Dynes,
Nivedita Mahesh,
Nicolas Fagnoni,
David R. DeBoer
Abstract:
Detecting cosmological signals from the Epoch of Reionization (EoR) requires high-precision calibration to isolate the cosmological signals from foreground emission. In radio interferometery, perturbed primary beams of antenna elements can disrupt the precise calibration, which results in contaminating the foreground-free region, or the EoR window, in the cylindrically averaged power spectrum. For…
▽ More
Detecting cosmological signals from the Epoch of Reionization (EoR) requires high-precision calibration to isolate the cosmological signals from foreground emission. In radio interferometery, perturbed primary beams of antenna elements can disrupt the precise calibration, which results in contaminating the foreground-free region, or the EoR window, in the cylindrically averaged power spectrum. For Hydrogen Epoch of Reionization Array (HERA), we simulate and characterize the perturbed primary beams induced by feed motions such as axial, lateral, and tilting motions, above the 14-meter dish. To understand the effect of the perturbed beams, visibility measurements are modeled with two different foreground components, point sources and diffuse sources, and we find different feed motions present a different reaction to each type of sky source. HERA's redundant-baseline calibration in the presence of non-redundant antenna beams due to feed motions introduces chromatic errors in gain solutions, which produces foreground power leakage into the EoR window. The observed leakage from vertical feed motions comes predominately from point sources around zenith. Furthermore, the observed leakage from horizontal and tilting feed motion comes predominately from the diffuse components near the horizon. Mitigation of chromatic gain errors will be necessary for robust detection of the EoR signals with minimal foreground bias, and this will be discussed in the subsequent paper.
△ Less
Submitted 28 October, 2022;
originally announced October 2022.
-
Lost Horizon: Quantifying the Effect of Local Topography on Global 21-cm Cosmology Data Analysis
Authors:
Neil Bassett,
David Rapetti,
Keith Tauscher,
Bang D. Nhan,
David D. Bordenave,
Joshua J. Hibbard,
Jack O. Burns
Abstract:
We present an investigation of the horizon and its effect on global 21-cm observations and analysis. We find that the horizon cannot be ignored when modeling low frequency observations. Even if the sky and antenna beam are known exactly, forward models cannot fully describe the beam-weighted foreground component without accurate knowledge of the horizon. When fitting data to extract the 21-cm sign…
▽ More
We present an investigation of the horizon and its effect on global 21-cm observations and analysis. We find that the horizon cannot be ignored when modeling low frequency observations. Even if the sky and antenna beam are known exactly, forward models cannot fully describe the beam-weighted foreground component without accurate knowledge of the horizon. When fitting data to extract the 21-cm signal, a single time-averaged spectrum or independent multi-spectrum fits may be able to compensate for the bias imposed by the horizon. However, these types of fits lack constraining power on the 21-cm signal, leading to large uncertainties on the signal extraction, in some cases larger in magnitude than the 21-cm signal itself. A significant decrease in signal uncertainty can be achieved by performing multi-spectrum fits in which the spectra are modeled simultaneously with common parameters. The cost of this greatly increased constraining power, however, is that the time dependence of the horizon's effect, which is more complex than its spectral dependence, must be precisely modeled to achieve a good fit. To aid in modeling the horizon, we present an algorithm and Python package for calculating the horizon profile from a given observation site using elevation data. We also address several practical concerns such as pixelization error, uncertainty in the horizon profile, and foreground obstructions such as surrounding buildings and vegetation. We demonstrate that our training set-based analysis pipeline can account for all of these factors to model the horizon well enough to precisely extract the 21-cm signal from simulated observations.
△ Less
Submitted 3 June, 2021;
originally announced June 2021.
-
Global 21-cm signal extraction from foreground and instrumental effects IV: Accounting for realistic instrument uncertainties and their overlap with foreground and signal models
Authors:
Keith Tauscher,
David Rapetti,
Bang D. Nhan,
Alec Handy,
Neil Bassett,
Joshua Hibbard,
David Bordenave,
Richard F. Bradley,
Jack O. Burns
Abstract:
All 21-cm signal experiments rely on electronic receivers that affect the data via both multiplicative and additive biases through the receiver's gain and noise temperature. While experiments attempt to remove these biases, the residuals of their imperfect calibration techniques can still confuse signal extraction algorithms. In this paper, the fourth and final installment of our pipeline series,…
▽ More
All 21-cm signal experiments rely on electronic receivers that affect the data via both multiplicative and additive biases through the receiver's gain and noise temperature. While experiments attempt to remove these biases, the residuals of their imperfect calibration techniques can still confuse signal extraction algorithms. In this paper, the fourth and final installment of our pipeline series, we present a technique for fitting out receiver effects as efficiently as possible. The fact that the gain and global signal, which are multiplied in the observation equation, must both be modeled implies that the model of the data is nonlinear in its parameters, making numerical sampling the only way to explore the parameter distribution rigorously. However, multi-spectra fits, which are necessary to extract the signal confidently as demonstrated in the third paper of the series, often require large numbers of foreground parameters, increasing the dimension of the posterior distribution that must be explored and therefore causing numerical sampling inefficiencies. Building upon techniques in the second paper of the series, we outline a method to explore the full parameter distribution by numerically sampling a small subset of the parameters and analytically marginalizing over the others. We test this method in simulation using a type-I Chebyshev band-pass filter gain model and a fast signal model based on a spline between local extrema. The method works efficiently, converging quickly to the posterior signal parameter distribution. The final signal uncertainties are of the same order as the noise in the data.
△ Less
Submitted 4 May, 2021;
originally announced May 2021.
-
Assessment of the Projection-induced Polarimetry Technique for Constraining the Foreground Spectrum in Global 21 cm Cosmology
Authors:
Bang D. Nhan,
David D. Bordenave,
Richard F. Bradley,
Jack O. Burns,
Keith Tauscher,
David Rapetti,
Patricia J. Klima
Abstract:
Detecting the cosmological sky-averaged (global) 21 cm signal as a function of observed frequency will provide a powerful tool to study the ionization and thermal history of the intergalactic medium (IGM) in the early Universe ($\sim$ 400 million years after the Big Bang). The greatest challenge in conventional total-power global 21 cm experiments is the removal of the foreground synchrotron emiss…
▽ More
Detecting the cosmological sky-averaged (global) 21 cm signal as a function of observed frequency will provide a powerful tool to study the ionization and thermal history of the intergalactic medium (IGM) in the early Universe ($\sim$ 400 million years after the Big Bang). The greatest challenge in conventional total-power global 21 cm experiments is the removal of the foreground synchrotron emission ($\sim 10^3$-$10^4$ K) to uncover the weak cosmological signal (tens to hundreds of mK), especially since the intrinsic smoothness of the foreground spectrum is corrupted by instrumental effects. Although the EDGES team has recently reported an absorption profile at 78 MHz in the sky-averaged spectrum, it is necessary to confirm this detection with an independent approach. The projection effect from observing anisotropic foreground source emission with a wide-view antenna pointing at the North Celestial Pole (NCP) can induce a net polarization, referred as the Projection-Induced Polarization Effect (PIPE). Due to Earth's rotation, observation centered at the circumpolar region will impose a dynamic sky modulation on the net polarization's waveforms which is unique to the foreground component. In this study, we review the implementation practicality and underlying instrumental effects of this new polarimetry-based technique with detailed numerical simulation and a testbed instrument, the Cosmic Twilight Polarimeter (CTP). In addition, we explore an SVD-based analysis approach for separating the foreground and instrumental effects from the background global 21 cm signal using the sky-modulated PIPE.
△ Less
Submitted 5 August, 2019; v1 submitted 12 November, 2018;
originally announced November 2018.
-
A polarimetric approach for constraining the dynamic foreground spectrum for cosmological global 21-cm measurements
Authors:
Bang D. Nhan,
Richard F. Bradley,
Jack O. Burns
Abstract:
The cosmological global (sky-averaged) 21-cm signal is a powerful tool to probe the evolution of the intergalactic medium (IGM) in high-redshift Universe ($z \leq 6$). One of the biggest observational challenges is to remove the foreground spectrum which is at least four orders of magnitude brighter than the cosmological 21-cm emission. Conventional global 21-cm experiments rely on the spectral sm…
▽ More
The cosmological global (sky-averaged) 21-cm signal is a powerful tool to probe the evolution of the intergalactic medium (IGM) in high-redshift Universe ($z \leq 6$). One of the biggest observational challenges is to remove the foreground spectrum which is at least four orders of magnitude brighter than the cosmological 21-cm emission. Conventional global 21-cm experiments rely on the spectral smoothness of the foreground synchrotron emission to separate it from the unique 21-cm spectral structures in a single total-power spectrum. However, frequency-dependent instrumental and observational effects are known to corrupt such smoothness and complicates the foreground subtraction. We introduce a polarimetric approach to measure the projection-induced polarization of the anisotropic foreground onto a stationary dual-polarized antenna. Due to Earth rotation, when pointing the antenna at a celestial pole, the revolving foreground will modulate this polarization with a unique frequency-dependent sinusoidal signature as a function of time. In our simulations, by harmonic decomposing this dynamic polarization, our technique produces two separate spectra in parallel from the same observation: (i) a total sky power consisting both the foreground and the 21-cm background, (ii) a model-independent measurement of the foreground spectrum at a harmonic consistent to twice the sky rotation rate. In the absence of any instrumental effects, by scaling and subtracting the latter from the former, we recover the injected global 21-cm model within assumed uncertainty. We further discuss several limiting factors and potential remedies for future implementation.
△ Less
Submitted 17 January, 2017; v1 submitted 18 November, 2016;
originally announced November 2016.