Calibration and wide field imaging with PAPER: a catalogue of compact sources

Observations of the redshifted 21 cm HI line promise to be a formidable tool for cosmology, allowing the investigation of the end of the so-called dark ages, when the first galaxies formed, and the subsequent Epoch of Reionization when the intergalactic medium transitioned from neutral to ionized. Such observations are plagued by foreground emission which is a few orders of magnitude brighter than the 21 cm line. In this thesis I analyzed data from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in order to improve the characterization of the extragalactic foreground component. I derived a catalogue of unresolved radio sources down to a 5 Jy flux density limit at 150 MHz and derived their spectral index distribution using literature data at 408 MHz. I implemented advanced techniques to calibrate radio interferometric data that led to a few percent accuracy on the flux density scale of the derived catalogue. This work, therefore, represents a further step towards creating an accurate, global sky model that is crucial to improve calibration of Epoch of Reionization observations

The design, construction and deployment of PRIᶻM.

Doctoral Degree. University of KwaZulu-Natal, Durban.Abstract available in PDF.Author's publication list can be found on page ii of the thesis. Publications that form part and/or include research presented in this thesis: Probing Radio Intensity at high-Z from Marion: 2017 Instrument (Philip et al., 2018

Captivity-induced behaviour and spatial learning abilities in an enigmatic, aquifer-dwelling blind eel, Rakthamichthys digressus

We investigated the impact of captive life on behaviour and learning abilities in an enigmatic, aquifer-dwelling blind eel, Rakthamichthys digressus. Of eight major behavioural traits related to exploration and activity in a novel arena, four were significantly altered by life in captivity. While the startle response upon introduction into the arena and overall swimming away from the walls increased after captivity, inactivity exhibited immediately after the startle and the reaction to an external disturbance decreased. We also observed behavioural syndromes between ‘startle responses’ and ‘horizontal wall following’, and between ‘overall activity’ and ‘vertical wall following’; however, these behavioural syndromes were not altered by maintenance in captivity. Interestingly, this blind-eel failed to learn a simple spatial task in a Y-maze apparatus. Captive-associated behavioural changes in R. digressus may influence their survival after reintroduction into natural habitats, and such changes must be taken into account while developing protocols for ex-situ conservation and subsequent release

COMAP Early Science: VIII. A Joint Stacking Analysis with eBOSS Quasars

We present a new upper limit on the cosmic molecular gas density at z = 2.4 − 3.4 obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 243 quasars selected from the Extended Baryon Oscillation SpectroscopicSurvey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.129 Jykm/s. Depending on the balance of the emission between the quasar host and its environment, this value can be interpreted as an average CO line luminosity L′CO of eBOSS quasars of ≤ 1.26 × 1011 K km pc2s−1, or an average molecular gas density ρH2 in regions of the universe containing a quasar of ≤ 1.52 × 108 M⊙ cMpc−3. The L′ CO upper limit falls among CO line luminosities obtained fromindividually-targeted quasars in the COMAP redshift range, and the ρH2 value is comparable to upper limits obtained from other Line Intensity Mapping (LIM) surveys and their joint analyses. Further, we forecast the values obtainable with the COMAP/eBOSS stack after the full 5-year COMAP Pathfinder survey. We predict that a detection is probable with this method, depending on the CO properties of the quasar sample. Based on the achieved sensitivity, we believe that this technique of stacking LIM data on the positions of traditional galaxy or quasar catalogs is extremely promising, both asa technique for investigating large galaxy catalogs efficiently at high redshift and as a technique for bolstering the sensitivity of LIM experiments, even with a fraction of their total expected survey data

COMAP Early Science: VIII. A Joint Stacking Analysis with eBOSS Quasars

We present a new upper limit on the cosmic molecular gas density at $z=2.4-3.4$ obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 282 quasars selected from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.210 Jy km/s. Depending on the assumptions made, this value can be interpreted as either an average CO line luminosity $L'_\mathrm{CO}$ of eBOSS quasars of $\leq 7.30\times10^{10}$ K km pc$^2$ s$^{-1}$, or an average molecular gas density $\rho_\mathrm{H_2}$ in regions of the universe containing a quasar of $\leq 2.02\times10^8$ M$_\odot$ cMpc$^{-3}$. The $L'_\mathrm{CO}$ upper limit falls among CO line luminosities obtained from individually-targeted quasars in the COMAP redshift range, and the $\rho_\mathrm{H_2}$ value is comparable to upper limits obtained from other Line Intensity Mapping (LIM) surveys and their joint analyses. Further, we forecast the values obtainable with the COMAP/eBOSS stack after the full 5-year COMAP Pathfinder survey. We predict that a detection is probable with this method, depending on the CO properties of the quasar sample. Based on these achieved sensitivities, we believe that this technique of stacking LIM data on the positions of traditional galaxy or quasar catalogs is extremely promising, both as a technique for investigating large galaxy catalogs efficiently at high redshift and as a technique for bolstering the sensitivity of LIM experiments, even with a fraction of their total expected survey data.Comment: 15 pages, 8 figures. To be submitted to Ap

COMAP Early Science: III. CO Data Processing

We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, including $\chi^2$ and multi-scale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a dataset with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers.Comment: Paper 3 of 7 in series. 26 pages, 23 figures, submitted to Ap

COMAP Early Science: IV. Power Spectrum Methodology and Results

We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scales $k=0.051-0.62 \,\mathrm{Mpc}^{-1}$ we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4\mu\textrm{K}^2\mathrm{Mpc}^3$, the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum in the literature.Comment: Paper 4 of 7 in series. 18 pages, 11 figures, as accepted in Ap

COMAP Early Science: VII. Prospects for CO Intensity Mapping at Reionization

We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts ($z\sim5-8$) in addition to providing a significant boost to the $z\sim3$ sensitivity of the Pathfinder. We examine a set of existing models of the EoR CO signal, and find power spectra spanning several orders of magnitude, highlighting our extreme ignorance about this period of cosmic history and the value of the COMAP-EoR measurement. We carry out the most detailed forecast to date of an intensity mapping cross-correlation, and find that five out of the six models we consider yield signal to noise ratios (S/N) $\gtrsim20$ for COMAP-EoR, with the brightest reaching a S/N above 400. We show that, for these models, COMAP-EoR can make a detailed measurement of the cosmic molecular gas history from $z\sim2-8$, as well as probe the population of faint, star-forming galaxies predicted by these models to be undetectable by traditional surveys. We show that, for the single model that does not predict numerous faint emitters, a COMAP-EoR-type measurement is required to rule out their existence. We briefly explore prospects for a third-generation Expanded Reionization Array (COMAP-ERA) capable of detecting the faintest models and characterizing the brightest signals in extreme detail.Comment: Paper 7 of 7 in series. 19 pages, 10 figures, to be submitted to Ap

COMAP Early Science: VI. A First Look at the COMAP Galactic Plane Survey

We present early results from the COMAP Galactic Plane Survey conducted between June 2019 and April 2021, spanning $20^\circ<\ell<40^\circ$ in Galactic longitude and |b|<1.\!\!^{\circ}5 in Galactic latitude with an angular resolution of $4.5^{\prime}$. The full survey will span $\ell \sim 20^{\circ}$- $220^{\circ}$ and will be the first large-scale radio continuum survey at $30$ GHz with sub-degree resolution. We present initial results from the first part of the survey, including diffuse emission and spectral energy distributions (SEDs) of HII regions and supernova remnants. Using low and high frequency surveys to constrain free-free and thermal dust emission contributions, we find evidence of excess flux density at $30\,$GHz in six regions that we interpret as anomalous microwave emission. Furthermore we model UCHII contributions using data from the $5\,$GHz CORNISH catalogue and reject this as the cause of the $30\,$GHz excess. Six known supernova remnants (SNR) are detected at $30\,$GHz, and we measure spectral indices consistent with the literature or show evidence of steepening. The flux density of the SNR W44 at $30\,$GHz is consistent with a power-law extrapolation from lower frequencies with no indication of spectral steepening in contrast with recent results from the Sardinia Radio Telescope. We also extract five hydrogen radio recombination lines to map the warm ionized gas, which can be used to estimate electron temperatures or to constrain continuum free-free emission. The full COMAP Galactic plane survey, to be released in 2023/2024, will be an invaluable resource for Galactic astrophysics.Comment: Paper 6 of 7 in series. 28 pages, 10 figures, submitted to Ap

COMAP Early Science: II. Pathfinder Instrument

Line intensity mapping (LIM) is a new technique for tracing the global properties of galaxies over cosmic time. Detection of the very faint signals from redshifted carbon monoxide (CO), a tracer of star formation, pushes the limits of what is feasible with a total-power instrument. The CO Mapping Project (COMAP) Pathfinder is a first-generation instrument aiming to prove the concept and develop the technology for future experiments, as well as delivering early science products. With 19 receiver channels in a hexagonal focal plane arrangement on a 10.4 m antenna, and an instantaneous 26-34 GHz frequency range with 2 MHz resolution, it is ideally suited to measuring CO($J$=1-0) from $z\sim3$. In this paper we discuss strategies for designing and building the Pathfinder and the challenges that were encountered. The design of the instrument prioritized LIM requirements over those of ancillary science. After a couple of years of operation, the instrument is well understood, and the first year of data is already yielding useful science results. Experience with this Pathfinder will drive the design of the next generations of experiments.Comment: Paper 2 of 7 in series. 27 pages, 28 figures, submitted to Ap