Direct Optimal Mapping Image Power Spectrum and its Window Functions

Zhilei Xu,Honggeun Kim,Jacqueline N. Hewitt,Kai-Feng Chen,Nicholas S. Kern,Eleanor Rath,Ruby Byrne,Adélie Gorce,Robert Pascua,Zachary E. Martinot,Joshua S. Dillon,Bryna J. Hazelton,Adrian Liu,Miguel F. Morales,Zara Abdurashidova,Tyrone Adams,James E. Aguirre,Paul Alexander,Zaki S. Ali,Rushelle Baartman,Yanga Balfour,Adam P. Beardsley,Gianni Bernardi,Tashalee S. Billings,Judd D. Bowman,Richard F. Bradley,Philip Bull,Jacob Burba,Steven Carey,Chris L. Carilli,Carina Cheng,David R. DeBoer,Eloy de Lera Acedo,Matt Dexter,Nico Eksteen,John Ely,Aaron Ewall-Wice,Nicolas Fagnoni,Randall Fritz,Steven R. Furlanetto,Kingsley Gale-Sides,Brian Glendenning,Deepthi Gorthi,Bradley Greig,Jasper Grobbelaar,Ziyaad Halday,Jack Hickish,Daniel C. Jacobs,Austin Julius,MacCalvin Kariseb,Joshua Kerrigan,Piyanat Kittiwisit,Saul A. Kohn,Matthew Kolopanis,Adam Lanman,Paul La Plante,Anita Loots,David Harold Edward MacMahon,Lourence Malan,Cresshim Malgas,Keith Malgas,Bradley Marero,Andrei Mesinger,Mathakane Molewa,Tshegofalang Mosiane,Steven G. Murray,Abraham R. Neben,Bojan Nikolic,Hans Nuwegeld,Aaron R. Parsons,Nipanjana Patra,Samantha Pieterse,Nima Razavi-Ghods,James Robnett,Kathryn Rosie,Peter Sims,Craig Smith,Hilton Swarts,Nithyanandan Thyagarajan,Pieter van Wyngaarden,Peter K. G. Williams,Haoxuan Zheng

2024-07-06

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

Instrumentation and Methods for Astrophysics,Cosmology and Nongalactic Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to separate the cosmological signal from the dominant foreground radiation when detecting the neutral hydrogen signal during the Epoch of Reionization (EoR). Specifically, the paper proposes the Direct Optimal Mapping (DOM) method for mapping the interferometric visibility data to an image and calculating the power spectrum and its window function based on these images.

### Problem Background

During the Epoch of Reionization, neutral hydrogen in the universe emits radiation in the 21 - centimeter band. After redshift, this radiation can be observed in different frequency channels, thereby tracking the matter distribution during the Epoch of Reionization. However, the foreground radiation (such as the continuous - spectrum radiation of the Milky Way and other galaxies) is several orders of magnitude stronger than the cosmological signal. Therefore, how to effectively separate the foreground and cosmological signals is a key challenge.

### Solution

The DOM method proposed in the paper solves this problem in the following ways:

1. **Linear - operation mapping**: The DOM algorithm maps the visibility data to an image using only linear operations, ensuring the accurate calculation of the point - spread function (PSF) and the pixel covariance matrix.
2. **Fast Fourier Transform (FFT)**: Based on the DOM - mapped image, use FFT to calculate the image power spectrum and its window function.
3. **Simulation verification**: The effectiveness of the DOM method is verified through noiseless simulation data. In particular, it is shown that the window function can control the power leakage in the foreground area at an extremely low level (< 10^-11), thereby effectively separating the foreground and EoR signals.

### Formula Summary

The core formulas of the DOM method include:

- **Point - spread function (PSF) matrix**:
\[
P_f = D_f A^\dagger_f N^{-1}_f A_f
\]
where \(D_f\) is the normalization matrix, \(A_f\) is the measurement matrix, and \(N_f\) is the noise covariance matrix.

- **Power spectrum estimator**:
\[
\hat{q}_\alpha = V \frac{\hat{m}^\dagger R^\dagger E_\alpha R \hat{m}}{N^2}
\]
where \(R\) is the 3D smoothing function of the image cube, and \(E_\alpha\) is the discrete Fourier transform operator.

- **Window function matrix**:
\[
W_{\alpha\beta} = \frac{H_{\alpha\beta}}{\sum_\delta H_{\alpha\delta}}
\]
where \(H\) is the power spectrum response function matrix.

Through these formulas, the DOM method not only provides an effective means of separating cosmological signals but also provides new tools and verification methods for future 21 - centimeter cosmological observations.