
21 cm cosmology is a promising new probe of the evolution of visible matter
in our universe, especially during the poorlyconstrained Cosmic Dawn and Epoch
of Reionization. However, in order to separate the 21 cm signal from bright
astrophysical foregrounds, we need an exquisite understanding of our telescopes
so as to avoid adding spectral structure to spectrallysmooth foregrounds. One
powerful calibration method relies on repeated simultaneous measurements of the
same interferometric baseline to solve for the sky signal and for instrumental
parameters simultaneously. However, certain degrees of freedom are not
constrained by asserting internal consistency between redundant measurements.
In this paper, we review the origin of these "degeneracies" of
redundantbaseline calibration and demonstrate how they can source unwanted
spectral structure in our measurement and show how to eliminate that
additional, artificial structure. We also generalize redundant calibration to
dualpolarization instruments, derive the degeneracy structure, and explore the
unique challenges to calibration and preserving spectral smoothness presented
by a polarized measurement.

Foreground power dominates the measurements of interferometers that seek a
statistical detection of highlyredshifted HI emission from the Epoch of
Reionization (EoR). The inherent spectral smoothness of synchrotron radiation,
the dominant foreground emission mechanism, and the chromaticity of the
instrument allows these experiments to delineate a boundary between spectrally
smooth and structured emission in Fourier space (the "wedge" or "pitchfork",
and the "EoR Window", respectively). Faraday rotation can inject spectral
structure into otherwise smooth polarized foreground emission, which through
instrument effects or miscalibration could possibly pollute the EoR Window.
Using data from the Hydrogen Epoch of Reionization Array (HERA) 19element
commissioning array, we investigate the polarization response of this new
instrument in the power spectrum domain. We confirm the expected structure of
foreground emission in Fourier space predicted by Thyagarajan et al. (2015a,
2016) for a HERAtype dish, and detect polarized power within the pitchfork.
Using simulations of the polarized response of HERA feeds, we find that almost
all of the power in Stokes Q, U and V can be attributed to instrumental leakage
effects. Power consistent with noise in the EoR window suggests a negligible
amount of spectrallystructured polarized power, to the noiselevels attained.
This lends confidence to deep integrations with HERA in the future, but with a
lower noise floor these future studies will also have to investigate their
polarized response.

We study the impact of skybased calibration errors from source mismodeling
on 21\,cm power spectrum measurements with an interferometer and propose a
method for suppressing their effects. While emission from faint sources that
are not accounted for in calibration catalogs is believed to be spectrally
smooth, deviations of true visibilities from model visibilities are not, due to
the inherent chromaticity of the interferometer's skyresponse (the "wedge").
Thus, unmodeled foregrounds, below the confusion limit of many instruments,
introduce frequency structure into gain solutions on the same lineofsight
scales on which we hope to observe the cosmological signal. We derive analytic
expressions describing these errors using linearized approximations of the
calibration equations and estimate the impact of this bias on measurements of
the 21\,cm power spectrum during the Epoch of Reionization (EoR). Given our
current precision in primary beam and foreground modeling, this noise will
significantly impact the sensitivity of existing experiments that rely on
skybased calibration. Our formalism describes the scaling of calibration with
array and skymodel parameters and can be used to guide future instrument
design and calibration strategy. We find that skybased calibration that
downweights long baselines can eliminate contamination in most of the region
outside of the wedge with only a modest increase in instrumental noise.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch
of Reionization (EoR) are limited predominantly by the chromatic instrumental
systematic effect. The delay spectrum methodology for 21 cm power spectrum
measurements brought new attention to the critical impact of an antenna's
chromaticity on the viability of making this measurement. This methodology
established a straightforward relationship between timedomain response of an
instrument and the power spectrum modes accessible to a 21 cm EoR experiment.
We examine the performance of a prototype of the Hydrogen Epoch of Reionization
Array (HERA) array element that is currently observing in Karoo desert, South
Africa. We present a mathematical framework to derive the beam integrated
frequency response of a HERA prototype element in reception from the return
loss measurements between 100200 MHz and determined the extent of additional
foreground contamination in the delay space. The measurement reveals excess
spectral structures in comparison to the simulation studies of the HERA
element. Combined with the HERA data analysis pipeline that incorporates
inverse covariance weighting in optimal quadratic estimation of power spectrum,
we find that in spite of its departure from the simulated response, HERA
prototype element satisfies the necessary criteria posed by the foreground
attenuation limits and potentially can measure the power spectrum at spatial
modes as low as $k_{\parallel} > 0.1h$~Mpc$^{1}$. The work highlights a
straightforward method for directly measuring an instrument response and
assessing its impact on 21 cm EoR power spectrum measurements for future
experiments that will use reflectortype antenna.

The Hydrogen Epoch of Reionization Array (HERA) is a radio interferometer
aiming to detect the power spectrum of 21 cm fluctuations from neutral hydrogen
from the Epoch of Reionization (EOR). Drawing on lessons from the Murchison
Widefield Array (MWA) and the Precision Array for Probing the Epoch of
Reionization (PAPER), HERA is a hexagonal array of large (14 m diameter) dishes
with suspended dipole feeds. Not only does the dish determine overall
sensitivity, it affects the observed frequency structure of foregrounds in the
interferometer. This is the first of a series of four papers characterizing the
frequency and angular response of the dish with simulations and measurements.
We focus in this paper on the angular response (i.e., power pattern), which
sets the relative weighting between sky regions of high and low delay, and
thus, apparent source frequency structure. We measure the angular response at
137 MHz using the ORBCOMM beam mapping system of Neben et al. We measure a
collecting area of 93 m^2 in the optimal dish/feed configuration, implying
HERA320 should detect the EOR power spectrum at z~9 with a signaltonoise
ratio of 12.7 using a foreground avoidance approach with a single season of
observations, and 74.3 using a foreground subtraction approach. Lastly we study
the impact of these beam measurements on the distribution of foregrounds in
Fourier space.

The Detection of redshifted 21 cm emission from the epoch of reionization
(EoR) is a challenging task owing to strong foregrounds that dominate the
signal. In this paper, we propose a general method, based on the delay spectrum
approach, to extract HI power spectra that is applicable to tracking
observations using an imaging radio interferometer (Delay Spectrum with Imaging
Arrays (DSIA)). Our method is based on modelling the HI signal taking into
account the impact of wide field effects such as the $w$term which are then
used as appropriate weights in crosscorrelating the measured visibilities. Our
method is applicable to any radio interferometer that tracks a phase center and
could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the
literature the delay spectrum approach has been implemented for nearredundant
baselines using drift scan observations. In this paper we explore the scheme
for nonredundant tracking arrays, and this is the first application of delay
spectrum methodology to such data to extract the HI signal. We analyze 3 hours
of MWA tracking data on the EoR1 field. We present both 2dimensional
($k_\parallel,k_\perp$) and 1dimensional (k) power spectra from the analysis.
Our results are in agreement with the findings of other pipelines developed to
analyse the MWA EoR data.

We present an improved Global Sky Model (GSM) of diffuse galactic radio
emission from 10 MHz to 5 THz, whose uses include foreground modeling for CMB
and 21 cm cosmology. Our model improves on past work both algorithmically and
by adding new data sets such as the Planck maps and the enhanced Haslam map.
Our method generalizes the Principal Component Analysis approach to handle
nonoverlapping regions, enabling the inclusion of 29 sky maps with no region
of the sky common to all. We also perform a blind separation of our GSM into
physical components with a method that makes no assumptions about physical
emission mechanisms (synchrotron, freefree, dust, etc). Remarkably, this blind
method automatically finds five components that have previously only been found
"by hand", which we identify with synchrotron, freefree, cold dust, warm dust,
and the CMB anisotropy, with maps and spectra agreeing with previous work but
in many cases with smaller error bars. The improved GSM is available online at
github.com/jeffzhen/gsm2016.

We use timedomain electromagnetic simulations to determine the spectral
characteristics of the Hydrogen Epoch of Reionization Arrays (HERA) antenna.
These simulations are part of a multifaceted campaign to determine the
effectiveness of the dish's design for obtaining a detection of redshifted 21
cm emission from the epoch of reionization. Our simulations show the existence
of reflections between HERA's suspended feed and its parabolic dish reflector
that fall below 40 dB at 150 ns and, for reasonable impedance matches, have a
negligible impact on HERA's ability to constrain EoR parameters. It follows
that despite the reflections they introduce, dishes are effective for
increasing the sensitivity of EoR experiments at relatively low cost. We find
that electromagnetic resonances in the HERA feed's cylindrical skirt, which is
intended to reduce cross coupling and beam ellipticity, introduces significant
power at large delays ($40$ dB at 200 ns) which can lead to some loss of
measurable Fourier modes and a modest reduction in sensitivity. Even in the
presence of this structure, we find that the spectral response of the antenna
is sufficiently smooth for delay filtering to contain foreground emission at
lineofsight wave numbers below $k_\parallel \lesssim 0.2$ $h$Mpc$^{1}$, in
the region where the current PAPER experiment operates. Incorporating these
results into a Fisher Matrix analysis, we find that the spectral structure
observed in our simulations has only a small effect on the tight constraints
HERA can achieve on parameters associated with the astrophysics of
reionization.

The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to
measure 21 cm emission from the primordial intergalactic medium (IGM)
throughout cosmic reionization ($z=612$), and to explore earlier epochs of our
Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes
heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is
designed to characterize the evolution of the 21 cm power spectrum to constrain
the timing and morphology of reionization, the properties of the first
galaxies, the evolution of largescale structure, and the early sources of
heating. The full HERA instrument will be a 350element interferometer in South
Africa consisting of 14m parabolic dishes observing from 50 to 250 MHz.
Currently, 19 dishes have been deployed on site and the next 18 are under
construction. HERA has been designated as an SKA Precursor instrument.
In this paper, we summarize HERA's scientific context and provide forecasts
for its key science results. After reviewing the current state of the art in
foreground mitigation, we use the delayspectrum technique to motivate
highlevel performance requirements for the HERA instrument. Next, we present
the HERA instrument design, along with the subsystem specifications that ensure
that HERA meets its performance requirements. Finally, we summarize the
schedule and status of the project. We conclude by suggesting that, given the
realities of foreground contamination, currentgeneration 21 cm instruments are
approaching their sensitivity limits. HERA is designed to bring both the
sensitivity and the precision to deliver its primary science on the basis of
proven foreground filtering techniques, while developing new subtraction
techniques to unlock new capabilities. The result will be a major step toward
realizing the widely recognized scientific potential of 21 cm cosmology.

The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch
of Reionization (EoR) data and now faces the challenge of overcoming foreground
and systematic contamination to reduce the data to a cosmological measurement.
We introduce several novel analysis techniques such as cable reflection
calibration, hyperresolution gridding kernels, diffuse foreground model
subtraction, and quality control methods. Each change to the analysis pipeline
is tested against a two dimensional power spectrum figure of merit to
demonstrate improvement. We incorporate the new techniques into a deep
integration of 32 hours of MWA data. This data set is used to place a
systematiclimited upper limit on the cosmological power spectrum of $\Delta^2
\leq 2.7 \times 10^4$ mK$^2$ at $k=0.27$ h~Mpc$^{1}$ and $z=7.1$, consistent
with other published limits, and a modest improvement (factor of 1.4) over
previous MWA results. From this deep analysis we have identified a list of
improvements to be made to our EoR data analysis strategies. These improvements
will be implemented in the future and detailed in upcoming publications.

Realizing the potential of 21 cm tomography to statistically probe the
intergalactic medium before and during the Epoch of Reionization requires large
telescopes and precise control of systematics. Nextgeneration telescopes are
now being designed and built to meet these challenges, drawing lessons from
firstgeneration experiments that showed the benefits of densely packed, highly
redundant arraysin which the same mode on the sky is sampled by many antenna
pairsfor achieving high sensitivity, precise calibration, and robust
foreground mitigation. In this work, we focus on the Hydrogen Epoch of
Reionization Array (HERA) as an interferometer with a dense, redundant core
designed following these lessons to be optimized for 21 cm cosmology. We show
how modestly supplementing or modifying a compact design like HERA's can still
deliver high sensitivity while enhancing strategies for calibration and
foreground mitigation. In particular, we compare the imaging capability of
several array configurations, both instantaneously (to address instrumental and
ionospheric effects) and with rotation synthesis (for foreground removal). We
also examine the effects that configuration has on calibratability using
instantaneous redundancy. We find that improved imaging with subaperture
sampling via "offgrid" antennas and increased angular resolution via farflung
"outrigger" antennas is possible with a redundantly calibratable array
configuration.

We present the 21 cm power spectrum analysis approach of the Murchison
Widefield Array Epoch of Reionization project. In this paper, we compare the
outputs of multiple pipelines for the purpose of validating statistical limits
cosmological hydrogen at redshifts between 6 and 12. Multiple, independent,
data calibration and reduction pipelines are used to make power spectrum limits
on a fiducial night of data. Comparing the outputs of imaging and power
spectrum stages highlights differences in calibration, foreground subtraction
and power spectrum calculation. The power spectra found using these different
methods span a space defined by the various tradeoffs between speed, accuracy,
and systematic control. Lessons learned from comparing the pipelines range from
the algorithmic to the prosaically mundane; all demonstrate the many pitfalls
of neglecting reproducibility. We briefly discuss the way these different
methods attempt to handle the question of evaluating a significant detection in
the presence of foregrounds.

We present first results from radio observations with the Murchison Widefield
Array seeking to constrain the power spectrum of 21 cm brightness temperature
fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three
hours of observations were conducted over two nights with significantly
different levels of ionospheric activity. We use these data to assess the
impact of systematic errors at low frequency, including the ionosphere and
radiofrequency interference, on a power spectrum measurement. We find that
after the 13 hours of integration presented here, our measurements at the
Murchison Radio Observatory are not limited by RFI, even within the FM band,
and that the ionosphere does not appear to affect the level of power in the
modes that we expect to be sensitive to cosmology. Power spectrum detections,
inconsistent with noise, due to fine spectral structure imprinted on the
foregrounds by reflections in the signalchain, occupy the spatial Fourier
modes where we would otherwise be most sensitive to the cosmological signal. We
are able to reduce this contamination using calibration solutions derived from
autocorrelations so that we achieve an sensitivity of $10^4$ mK on comoving
scales $k\lesssim 0.5 h$Mpc$^{1}$. This represents the first upper limits on
the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 18$
but is still limited by calibration systematics. While calibration improvements
may allow us to further remove this contamination, our results emphasize that
future experiments should consider carefully the existence of and their ability
to calibrate out any spectral structure within the EoR window.

We use the Fisher matrix formalism and seminumerical simulations to derive
quantitative predictions of the constraints that power spectrum measurements on
nextgeneration interferometers, such as the Hydrogen Epoch of Reionization
Array (HERA) and the Square Kilometre Array (SKA), will place on the
characteristics of the Xray sources that heated the high redshift
intergalactic medium. Incorporating observations between $z=5$ and $z=25$, we
find that the proposed 331 element HERA and SKA phase 1 will be capable of
placing $\lesssim 10\%$ constraints on the spectral properties of these first
Xray sources, even if one is unable to perform measurements within the
foreground contaminated "wedge" or the FM band. When accounting for the
enhancement in power spectrum amplitude from spin temperature fluctuations, we
find that the observable signatures of reionization extend well beyond the peak
in the power spectrum usually associated with it. We also find that lower
redshift degeneracies between the signatures of heating and reionization
physics lead to errors on reionization parameters that are significantly
greater than previously predicted. Observations over the heating epoch are able
to break these degeneracies and improve our constraints considerably. For these
two reasons, 21\,cm observations during the heating epoch significantly enhance
our understanding of reionization as well.

Detection of the cosmological neutral hydrogen signal from the Epoch of
Reionization, and estimation of its basic physical parameters, is the principal
scientific aim of many current lowfrequency radio telescopes. Here we describe
the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed
and implemented with data from the Murchison Widefield Array (MWA), to compute
the twodimensional and sphericallyaveraged power spectrum of brightness
temperature fluctuations. The principal motivations for CHIPS are the
application of realistic instrumental and foreground models to form the optimal
estimator, thereby maximising the likelihood of unbiased signal estimation, and
allowing a full covariant understanding of the outputs. CHIPS employs an
inversecovariance weighting of the data through the maximum likelihood
estimator, thereby allowing use of the full parameter space for signal
estimation ("foreground suppression"). We describe the motivation for the
algorithm, implementation, application to real and simulated data, and early
outputs. Upon application to a set of 3 hours of data, we set a 2$\sigma$ upper
limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{1}$ of
$\Delta_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.26.6]$,
consistent with previous estimates.

We confirm our recent prediction of the "pitchfork" foreground signature in
power spectra of highredshift 21 cm measurements where the interferometer is
sensitive to largescale structure on all baselines. This is due to the
inherent response of a widefield instrument and is characterized by enhanced
power from foreground emission in Fourier modes adjacent to those considered to
be the most sensitive to the cosmological H I signal. In our recent paper, many
signatures from the simulation that predicted this feature were validated
against Murchison Widefield Array (MWA) data, but this key pitchfork signature
was close to the noise level. In this paper, we improve the data sensitivity
through the coherent averaging of 12 independent snapshots with identical
instrument settings and provide the first confirmation of the prediction with a
signaltonoise ratio > 10. This widefield effect can be mitigated by careful
antenna designs that suppress sensitivity near the horizon. Simple models for
antenna apertures that have been proposed for future instruments such as the
Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate
they should suppress foreground leakage from the pitchfork by ~40 dB relative
to the MWA and significantly increase the likelihood of cosmological signal
detection in these critical Fourier modes in the threedimensional power
spectrum.

21 cm cosmology, the statistical observation of the high redshift universe
using the hyperfine transition of neutral hydrogen, has the potential to
revolutionize our understanding of cosmology and the astrophysical processes
that underlie the formation of the first stars, galaxies, and black holes
during the "Cosmic Dawn." By making tomographic maps with low frequency radio
interferometers, we can study the evolution of the 21 cm signal with time and
spatial scale and use it to understand the density, temperature, and ionization
evolution of the intergalactic medium over this dramatic period in the history
of the universe.
For my Ph.D. thesis, I explore a number of advancements toward detecting and
characterizing the 21 cm signal from the Cosmic Dawn, especially during its
final stage, the epoch of reionization. In seven different previously published
papers, I explore new techniques for the statistical analysis of
interferometric measurements, apply them to data from current generation
telescopes like the Murchison Widefield Array, and look forward to what we
might measure with the next generation of 21 cm observatories. I focus in
particular on estimating the power spectrum of 21 cm brightness temperature
fluctuations in the presence enormous astrophysical foregrounds and how those
measurements may constrain the physics of the Cosmic Dawn.
Thesis Supervisor: Max Tegmark

The separation of the faint cosmological background signal from bright
astrophysical foregrounds remains one of the most daunting challenges of
mapping the highredshift intergalactic medium with the redshifted 21 cm line
of neutral hydrogen. Advances in mapping and modeling of diffuse and point
source foregrounds have improved subtraction accuracy, but no subtraction
scheme is perfect. Precisely quantifying the errors and error correlations due
to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm
power spectrum and for the maximal isolation of the "EoR window" from
foreground contamination. We present a method to infer the covariance of
foreground residuals from the data itself in contrast to previous attempts at a
priori modeling. We demonstrate our method by setting limits on the power
spectrum using a 3 h integration from the 128tile Murchison Widefield Array.
Observing between 167 and 198 MHz, we find at 95% confidence a best limit of
Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^1 and at z = 6.8,
consistent with existing limits.

Detection of 21~cm emission of HI from the epoch of reionization, at
redshifts z>6, is limited primarily by foreground emission. We investigate the
signatures of widefield measurements and an allsky foreground model using the
delay spectrum technique that maps the measurements to foreground object
locations through signal delays between antenna pairs. We demonstrate
interferometric measurements are inherently sensitive to all scales, including
the largest angular scales, owing to the nature of widefield measurements.
These widefield effects are generic to all observations but antenna shapes
impact their amplitudes substantially. A dishshaped antenna yields the most
desirable features from a foreground contamination viewpoint, relative to a
dipole or a phased array. Comparing data from recent Murchison Widefield Array
observations, we demonstrate that the foreground signatures that have the
largest impact on the HI signal arise from power received far away from the
primary field of view. We identify diffuse emission near the horizon as a
significant contributing factor, even on wide antenna spacings that usually
represent structures on small scales. For signals entering through the primary
field of view, compact emission dominates the foreground contamination. These
two mechanisms imprint a characteristic "pitchfork" signature on the
"foreground wedge" in Fourier delay space. Based on these results, we propose
that selective downweighting of data based on antenna spacing and time can
mitigate foreground contamination substantially by a factor ~100 with
negligible loss of sensitivity.

In order to study the "Cosmic Dawn" and the Epoch of Reionization with 21 cm
tomography, we need to statistically separate the cosmological signal from
foregrounds known to be orders of magnitude brighter. Over the last few years,
we have learned much about the role our telescopes play in creating a
putatively foregroundfree region called the "EoR window." In this work, we
examine how an interferometer's effects can be taken into account in a way that
allows for the rigorous estimation of 21 cm power spectra from interferometric
maps while mitigating foreground contamination and thus increasing sensitivity.
This requires a precise understanding of the statistical relationship between
the maps we make and the underlying true sky. While some of these calculations
would be computationally infeasible if performed exactly, we explore several
wellcontrolled approximations that make mapmaking and the calculation of map
statistics much faster, especially for compact and highly redundant
interferometers designed specifically for 21 cm cosmology. We demonstrate the
utility of these methods and the parametrized tradeoffs between accuracy and
speed using one such telescope, the upcoming Hydrogen Epoch of Reionization
Array, as a case study.

A number of experiments are currently working towards a measurement of the 21
cm signal from the Epoch of Reionization. Whether or not these experiments
deliver a detection of cosmological emission, their limited sensitivity will
prevent them from providing detailed information about the astrophysics of
reionization. In this work, we consider what types of measurements will be
enabled by a nextgeneration of larger 21 cm EoR telescopes. To calculate the
type of constraints that will be possible with such arrays, we use simple
models for the instrument, foreground emission, and the reionization history.
We focus primarily on an instrument modeled after the $\sim 0.1~\rm{km}^2$
collecting area Hydrogen Epoch of Reionization Array (HERA) concept design, and
parameterize the uncertainties with regard to foreground emission by
considering different limits to the recently described "wedge" footprint in
kspace. Uncertainties in the reionization history are accounted for using a
series of simulations which vary the ionizing efficiency and minimum virial
temperature of the galaxies responsible for reionization, as well as the mean
free path of ionizing photons through the IGM. Given various combinations of
models, we consider the significance of the possible power spectrum detections,
the ability to trace the power spectrum evolution versus redshift, the
detectability of salient power spectrum features, and the achievable level of
quantitative constraints on astrophysical parameters. Ultimately, we find that
$0.1~\rm{km}^2$ of collecting area is enough to ensure a very high significance
($\gtrsim30\sigma$) detection of the reionization power spectrum in even the
most pessimistic scenarios. This sensitivity should allow for meaningful
constraints on the reionization history and astrophysical parameters,
especially if foreground subtraction techniques can be improved and
successfully implemented.

We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder
lowfrequency radio interferometer whose goal is to test technologies that
improve the calibration precision and reduce the cost of the highsensitivity
3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using
massive baseline redundancy, which enables both automated precision calibration
and correlator cost reduction. We demonstrate and quantify the power and
robustness of redundancy for scalability and precision. We find that the
calibration parameters precisely describe the effect of the instrument upon our
measurements, allowing us to form a model that is consistent with $\chi^2$ per
degree of freedom < 1.2 for as much as 80% of the observations. We use these
results to develop an optimal estimator of calibration parameters using Wiener
filtering, and explore the question of how often and how finely in frequency
visibilities must be reliably measured to solve for calibration coefficients.
The success of MITEoR with its 64 dualpolarization elements bodes well for the
more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other
nextgeneration instruments, which would incorporate many identical or similar
technologies.

We describe a new technique for constraining the radio loud population of
active galactic nuclei at high redshift by measuring the imprint of 21 cm
spectral absorption features (the 21 cm forest) on the 21 cm power spectrum.
Using seminumerical simulations of the intergalactic medium and a
semiempirical source population we show that the 21 cm forest dominates a
distinctive region of $k$space, $k \gtrsim 0.5 \text{Mpc}^{1}$. By simulating
foregrounds and noise for current and potential radio arrays, we find that a
next generation instrument with a collecting area on the order of $\sim
0.1\text{km}^2$ (such as the Hydrogen Epoch of Reionization Array) may
separately constrain the Xray heating history at large spatial scales and
radio loud active galactic nuclei of the model we study at small ones. We
extrapolate our detectability predictions for a single radio loud active
galactic nuclei population to arbitrary source scenarios by analytically
relating the 21 cm forest power spectrum to the optical depth power spectrum
and an integral over the radio luminosity function.

We present techniques for bridging the gap between idealized inverse
covariance weighted quadratic estimation of 21 cm power spectra and the
realworld challenges presented universally by interferometric observation. By
carefully evaluating various estimators and adapting our techniques for large
but incomplete data sets, we develop a robust power spectrum estimation
framework that preserves the socalled "EoR window" and keeps track of
estimator errors and covariances. We apply our method to observations from the
32tile prototype of the Murchinson Widefield Array to demonstrate the
importance of a judicious analysis technique. Lastly, we apply our method to
investigate the dependence of the clean EoR window on frequencyespecially the
frequency dependence of the socalled "wedge" featureand establish upper
limits on the power spectrum from z = 6.2 to z = 11.7. Our lowest limit is
Delta(k) < 0.3 Kelvin at 95% confidence at a comoving scale k = 0.046 Mpc^1
and z = 9.5.

Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral
hydrogen has the potential to overtake the cosmic microwave background as our
most powerful cosmological probe, because it can map a much larger volume of
our Universe, shedding new light on the epoch of reionization, inflation, dark
matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder
lowfrequency radio interferometer whose goal is to test technologies that
greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR
accomplishes this by using massive baseline redundancy both to enable automated
precision calibration and to cut the correlator cost scaling from N^2 to NlogN,
where N is the number of antennas. The success of MITEoR with its 64
dualpolarization elements bodes well for the more ambitious HERA project,
which would incorporate many identical or similar technologies using an order
of magnitude more antennas, each with dramatically larger collecting area.