Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

Future generations of radio interferometers targeting the 21\,cm with $N\gg 1000$ antennas could face a significant computational challenge in building correlators with the traditional architecture, whose computational resource requirement scales as $\mathcal{O}(N^2)$ with array size. The fundamental output of such correlators are the cross-correlation products of all antenna pairs in the array. The FFT-correlator architecture reduces the computational resources scaling to $\mathcal{O}(N\log_2N)$ by computing cross-correlation products through a spatial Fourier transform. However, the output of the FFT-correlator is meaningful only when the input antenna voltages are gain- and phase-calibrated. Traditionally, interferometric calibration has used the $\mathcal{O}(N^2)$ cross-correlations produced by a standard correlator. This paper proposes calibration schemes for a real-time calibrator that could work in parallel with an FFT-correlator as a self-contained $\mathcal{O}(N\log_2{N})$ correlator system that can be scaled to large-N redundant arrays. We discuss two calibration schemes that be employed by this calibrator, and compare their performance and scalability. We find that both the calibration schemes presented in this paper result in gain variance that decreases as $1/\log_2{(N)}$ with increase in the size of the array.

We introduce DAYENU, a linear, spectral filter for HI intensity mapping that achieves the desirable foreground mitigation and error minimization properties of inverse co-variance weighting with minimal modeling of the underlying data. Beyond 21 cm power-spectrum estimation, our filter is suitable for any analysis where high dynamic-range removal of spectrally smooth foregrounds in irregularly (or regularly) sampled data is required, something required by many other intensity mapping techniques. We show that DAYENU enables the access of large-scale line-of-sight modes that are inaccessible to tapered DFT estimators. Since these modes have the largest SNRs, DAYENU significantly increases the sensitivity of 21 cm analyses over tapered Fourier transforms. Slight modifications allow us to use DAYENU as a linear replacement for iterative delay CLEANing (DAYENUREST). An interactive jupyter tutorial on using DAYENU can be found at https://github.com/HERA-Team/uvtools/blob/master/examples/linear_clean_demo.ipynb DAYENU's source code can be found at https://github.com/HERA-Team/uvtools/blob/master/uvtools/dspec.py

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.

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 (EoR). HERA is a drift-scan array with a 10 degree 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 catalogues 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 EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

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 EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

We study the behavior of internal signal chain reflections and antenna cross coupling as systematics for 21 cm cosmological surveys. We outline the mathematics for how these systematics appear in interferometric visibilities and describe their phenomenology. We then describe techniques for modeling and removing these systematics without attenuating the 21 cm signal in the data. This has critical implications for low-frequency radio surveys aiming to characterize the 21cm signal from the Epoch of Reionization and Cosmic Dawn, as systematics can cause bright foreground emission to contaminate the EoR window and prohibit a robust detection. We also quantify the signal loss properties of the systematic modeling algorithms, and show that our techniques demonstrate resistance against EoR signal loss. In a companion paper, we demonstrate these methods on data from the Hydrogen Epoch of Reionization Array as a proof-of-concept.

We present limits on the 21cm power spectrum from the Epoch of Reionization (EoR) using data from the 64 antenna configuration of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) analyzed through a power spectrum pipeline independent from previous PAPER analyses. Previously reported results from PAPER have been found to contain significant signal loss (Cheng et al. 2018, arxiv:1810.05175). Several lossy steps from previous PAPER pipelines have not been included in this analysis, namely: delay-based foreground filtering, optimal fringe-rate filtering, and empirical covariance-based estimators. Steps which remain in common with previous analyses include redundant calibration and local sidereal time (LST) binning. The power spectra reported here are effectively the result of applying a linear Fourier transform analysis to the calibrated, LST binned data. This analysis also uses more data than previous publications, including the complete available redshift range of $z \sim 7.5$ to $11$. In previous PAPER analyses, many power spectrum measurements were found to be detections of noncosmological power at levels of significance ranging from two to hundreds of times the theoretical noise. Here, excess power is examined using redundancy between baselines and power spectrum jackknives. The upper limits we find on the 21cm power spectrum from reionization are ($1500$ mK)$^{2}$, ($1900$ mK)$^{2}$, ($280$ mK)$^{2}$, ($200$ mK)$^{2}$, ($380$ mK)$^{2}$, ($300$ mK)$^{2}$ at redshifts $z=10.87,\ 9.93,\ 8.68,\ 8.37,\ 8.13,$ and $7.48$, respectively. For reasons described in Cheng et al. 2018 (arxiv:1810.05175), these limits supersede all previous PAPER results (Ali et al. 2018, arxiv:1502.06016).

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Therefore, it is essential to understand and limit these chromatic effects. This paper describes a method to simulate the frequency and time responses of an antenna, by simultaneously taking into account the analogue RF receiver, the transmission cable, and the mutual coupling caused by adjacent antennas. Applied to the Hydrogen Epoch of Reionization Array (HERA), this study reveals the presence of significant reflections at high delays caused by the 150-m cable which links the antenna to the back-end. Besides, it shows that waves can propagate from one dish to another one through large sections of the array because of mutual coupling. In this more realistic approach, the simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h \rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity.

In this white paper, we lay out a US roadmap for high-redshift 21 cm cosmology (30 < z < 6) in the 2020s. Beginning with the currently-funded HERA and MWA Phase II projects and advancing through the decade with a coordinated program of small-scale instrumentation, software, and analysis projects targeting technology development, this roadmap incorporates our current best understanding of the systematics confronting 21 cm cosmology into a plan for overcoming them, enabling next-generation, mid-scale 21 cm arrays to be proposed late in the decade. Submitted for consideration by the Astro2020 Decadal Survey Program Panel for Radio, Millimeter, and Submillimeter Observations from the Ground as a Medium-Sized Project.

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI 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. 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. 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 the expected level. Further, 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.

The separation of cosmological signal from astrophysical foregrounds is a fundamental challenge for any effort to probe the evolution of neutral hydrogen during the Cosmic Dawn and epoch of reionization (EoR) using the 21 cm hyperfine transition. Foreground separation is made possible by their intrinsic spectral smoothness, making them distinguishable from spectrally complex cosmological signal even though they are ~5 orders of magnitude brighter. Precisely calibrated radio interferometers are essential to maintaining the smoothness and thus separability of the foregrounds. One powerful calibration strategy is to use redundant measurements between pairs of antennas with the same physical separation in order to solve for each antenna's spectral response without reference to a sky model. This strategy is being employed by the Hydrogen Epoch of Reionization Array (HERA), a large radio telescope in South Africa that is now observing while being built out to 350 14-m dishes. However, the deviations from perfect redundancy inherent in any real radio telescope complicate the calibration problem. Using simulations of HERA, we show how calibration with antenna-to-antenna variations in dish construction and placement generally lead to spectral structure in otherwise smooth foregrounds that significantly reduces the number of cosmological modes available to a 21 cm measurement. However, we also show that this effect can be largely eliminated by a modified redundant-baseline calibration strategy that relies predominantly on short baselines.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Observations of the 21cm Epoch of Reionization (EoR) signal are dominated by Galactic and extragalactic foregrounds. The need for foreground removal has led to the development of two main techniques, often referred to as "foreground avoidance" and "foreground subtraction." Avoidance is associated with filtering foregrounds in Fourier space, while subtraction uses an explicit foreground model that is removed. Using 1088 hours of data from the 64-element PAPER array, we demonstrate that subtraction of a foreground model prior to delay-space foreground filtering results in a modest but measurable improvement of the performance of the filter. This proof-of-concept result shows that improvement stems from the reduced dynamic range requirements needed for the foreground filter: subtraction of a foreground model reduces the total foreground power, so for a fixed dynamic range, the filter can push towards fainter limits. We also find that the choice of window function used in the foreground filter can have an appreciable affect on the performance near the edges of the observing band. We demonstrate these effects using a smaller 3 hour sampling of data from the MWA, and find that the hybrid filtering and subtraction removal approach provides similar improvements across the band as seen in the case with PAPER-64.

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 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" dataset 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$\times 10^{5}$ HERA time-ordered 1024 channeled 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 $F_{2}$ score of 0.75 as applied to our HERA-67 observations.

The Epoch of Reionization (EoR) is an uncharted era in our Universe's history during which the birth of the first stars and galaxies led to the ionization of neutral hydrogen in the intergalactic medium. There are many experiments investigating the EoR by tracing the 21cm line of neutral hydrogen. Because this signal is very faint and difficult to isolate, it is crucial to develop analysis techniques that maximize sensitivity and suppress contaminants in data. It is also imperative to understand the trade-offs between different analysis methods and their effects on power spectrum estimates. Specifically, with a statistical power spectrum detection in HERA's foreseeable future, it has become increasingly important to understand how certain analysis choices can lead to the loss of the EoR signal. In this paper, we focus on signal loss associated with power spectrum estimation. We describe the origin of this loss using both toy models and data taken by the 64-element configuration of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER). In particular, we highlight how detailed investigations of signal loss have led to a revised, higher 21cm power spectrum upper limit from PAPER-64. Additionally, we summarize errors associated with power spectrum error estimation that were previously unaccounted for. We focus on a subset of PAPER-64 data in this paper; revised power spectrum limits from the PAPER experiment are presented in a forthcoming paper by Kolopanis et al. (in prep.) and supersede results from previously published PAPER analyses.

Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes.

In this paper, we report new limits on 21cm emission from cosmic reionization based on a 135-day observing campaign with a 64-element deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. This work extends the work presented in Parsons et al. (2014) with more collecting area, a longer observing period, improved redundancy-based calibration, optimal fringe-rate filtering, and improved power-spectral analysis using optimal quadratic estimators. The result is a new $2\sigma$ upper limit on $\Delta^{2}(k)$ of (22.4 mK)$^2$ in the range $0.15 < k < 0.5h\ {\rm Mpc}^{-1}$ at $z = 8.4$. This represents a three-fold improvement over the previous best upper limit. As we discuss in more depth in a forthcoming paper (Pober et al. 2015, in prep), this upper limit supports and extends previous evidence against extremely cold reionization scenarios. We conclude with a discussion of implications for future 21cm reionization experiments, including the newly funded Hydrogen Epoch of Reionization Array (HERA). $\textbf{The limits presented in this paper have been retracted: The erratum can be found in Appendix A.}$

New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.

21 cm cosmology is a promising new probe of the evolution of visible matter in our universe, especially during the poorly-constrained 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 spectrally-smooth 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 redundant-baseline 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 dual-polarization instruments, derive the degeneracy structure, and explore the unique challenges to calibration and preserving spectral smoothness presented by a polarized measurement.

Radio interferometers designed to measure the cosmological 21 cm power spectrum require high sensitivity. Several modern low-frequency interferometers feature drift-scan antennas placed on a regular grid to maximize the number of instantaneously coherent (redundant) measurements. However, even for such maximum-redundancy arrays, significant sensitivity comes through partial coherence between baselines. Current visibility-based power-spectrum pipelines, though shown to ease control of systematics, lack the ability to make use of this partial redundancy. We introduce a method to leverage partial redundancy in such power-spectrum pipelines for drift-scan arrays. Our method cross-multiplies baseline pairs at a time lag and quantifies the sensitivity contributions of each pair of baselines. Using the configurations and beams of the 128-element Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER-128) and staged deployments of the Hydrogen Epoch of Reionization Array, we illustrate how our method applies to different arrays and predict the sensitivity improvements associated with pairing partially coherent baselines. As the number of antennas increases, we find partial redundancy to be of increasing importance in unlocking the full sensitivity of upcoming arrays.

Observations of the redshifted 21-cm line from the epoch of reionization have recently motivated the construction of low frequency radio arrays with highly redundant configurations. These configurations provide an alternative calibration strategy - "redundant calibration" - and boosts sensitivity on specific spatial scales. In this paper, we formulate calibration of redundant interferometric arrays as a complex optimization problem. We solve this optimization problem via the Levenberg-Marquardt algorithm. This calibration approach is more robust to initial conditions than current algorithms and, by leveraging an approximate matrix inversion, allows for further optimization and an efficient implementation ("redundant StEfCal"). We also investigated using the preconditioned conjugate gradient method as an alternative to the approximate matrix inverse, but found that its computational performance is not competitive with respect to "redundant StEfCal". The efficient implementation of this new algorithm is made publicly available.

Current and upcoming radio interferometric experiments are aiming to make a statistical characterization of the high-redshift 21cm fluctuation signal spanning the hydrogen reionization and X-ray heating epochs of the universe. However, connecting 21cm statistics to underlying physical parameters is complicated by the theoretical challenge of modeling the relevant physics at computational speeds quick enough to enable exploration of the high dimensional and weakly constrained parameter space. In this work, we use machine learning algorithms to build a fast emulator that mimics expensive simulations of the 21cm signal across a wide parameter space to high precision. We embed our emulator within a Markov-Chain Monte Carlo framework, enabling it to explore the posterior distribution over a large number of model parameters, including those that govern the Epoch of Reionization, the Epoch of X-ray Heating, and cosmology. As a worked example, we use our emulator to present an updated parameter constraint forecast for the Hydrogen Epoch of Reionization Array experiment, showing that its characterization of a fiducial 21cm power spectrum will considerably narrow the allowed parameter space of reionization and heating parameters, and could help strengthen Planck's constraints on $\sigma_8$. We provide both our generalized emulator code and its implementation specifically for 21cm parameter constraints as publicly available software.

A critical challenge in the observation of the redshifted 21-cm line is its separation from bright Galactic and extragalactic foregrounds. In particular, the instrumental leakage of polarized foregrounds, which undergo significant Faraday rotation as they propagate through the interstellar medium, may harmfully contaminate the 21-cm power spectrum. We develop a formalism to describe the leakage due to instrumental widefield effects in visibility-based power spectra measured with redundant arrays, extending the delay-spectrum approach presented in Parsons et al. (2012). We construct polarized sky models and propagate them through the instrument model to simulate realistic full-sky observations with the Precision Array to Probe the Epoch of Reionization. We find that the leakage due to a population of polarized point sources is expected to be higher than diffuse Galactic polarization at any $k$ mode for a 30~m reference baseline. For the same reference baseline, a foreground-free window at $k > 0.3 \, h$~Mpc$^{-1}$ can be defined in terms of leakage from diffuse Galactic polarization even under the most pessimistic assumptions. If measurements of polarized foreground power spectra or a model of polarized foregrounds are given, our method is able to predict the polarization leakage in actual 21-cm observations, potentially enabling its statistical subtraction from the measured 21-cm power spectrum.

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 time-domain 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 100-200 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 reflector-type antenna.

Polarized foreground emission is a potential contaminant of attempts to measure the fluctuation power spectrum of highly redshifted 21 cm HI emission from the epoch of reionization. Using the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), we present limits on the observed power spectra of all four Stokes parameters in two frequency bands, centered at 126 MHz ($z=10.3$) and 164 MHz ($z=7.66$) for a three-month observing campaign of a 32-antenna deployment, for which unpolarized power spectrum results have been reported at $z=7.7$ (Parsons et al 2014) and $7.5 < z < 10.5$ (Jacobs et al 2014). The power spectra in this paper are processed in the same way as in those works, and show no definitive detection of polarized power. This non-detection appears to be largely due to the suppression of polarized power by ionospheric rotation measure fluctuations, which strongly affect Stokes Q and U. We are able to show that the net effect of polarized leakage is a negligible contribution at the levels of of the limits reported in Parsons et al 2014 and Jacobs et al 2014.

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 HERA-320 should detect the EOR power spectrum at z~9 with a signal-to-noise 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 Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) has been working for a decade to reduce the time and cost of designing, building and deploying new digital radio-astronomy instruments. Today, CASPER open-source technology powers over 45 scientific instruments worldwide, and is used by scientists and engineers at dozens of academic institutions. In this paper we catalog the current offerings of the CASPER collaboration, and instruments past and present built by CASPER users and developers. We describe the ongoing state of software development, as CASPER looks to support a broader range of programming environments and hardware and ensure compatibility with the latest vendor tools.

Intensity mapping is a promising technique for surveying the large scale structure of our Universe from $z=0$ to $z \sim 150$, using the brightness temperature field of spectral lines to directly observe previously unexplored portions of out cosmic timeline. Examples of targeted lines include the $21\,\textrm{cm}$ hyperfine transition of neutral hydrogen, rotational lines of carbon monoxide, and fine structure lines of singly ionized carbon. Recent efforts have focused on detections of the power spectrum of spatial fluctuations, but have been hindered by systematics such as foreground contamination. This has motivated the decomposition of data into Fourier modes perpendicular and parallel to the line-of-sight, which has been shown to be a particularly powerful way to diagnose systematics. However, such a method is well-defined only in the limit of a narrow-field, flat-sky approximation. This limits the sensitivity of intensity mapping experiments, as it means that wide surveys must be separately analyzed as a patchwork of smaller fields. In this paper, we develop a framework for analyzing intensity mapping data in a spherical Fourier-Bessel basis, which incorporates curved sky effects without difficulty. We use our framework to generalize a number of techniques in intensity mapping data analysis from the flat sky to the curved sky. These include visibility-based estimators for the power spectrum, treatments of interloper lines, and the "foreground wedge" signature of spectrally smooth foregrounds.

To evaluate the feasibility of measuring differences in bulk composition among carbonaceous meteorite parent bodies from an asteroid or comet orbiter, we present the results of a performance simulation of an orbital gamma-ray spectroscopy ("GRS") experiment in a Dawn-like orbit around spherical model asteroids with a range of carbonaceous compositions. The orbital altitude was held equal to the asteroid radius for 4.5 months. Both the asteroid gamma-ray spectrum and the spacecraft background flux were calculated using the MCNPX Monte-Carlo code. GRS is sensitive to depths below the optical surface (to ~20--50 cm depth depending on material density). This technique can therefore measure underlying compositions beneath a sulfur-depleted (e.g., Nittler et al. 2001) or desiccated surface layer. We find that 3\sigma\ uncertainties of under 1 wt% are achievable for H, C, O, Si, S, Fe, and Cl for five carbonaceous meteorite compositions using the heritage Mars Odyssey GRS design in a spacecraft- deck-mounted configuration at the Odyssey end-of-mission energy resolution, FWHM = 5.7 keV at 1332 keV. The calculated compositional uncertainties are smaller than the compositional differences between carbonaceous chondrite subclasses.

We use time-domain electromagnetic simulations to determine the spectral characteristics of the Hydrogen Epoch of Reionization Arrays (HERA) antenna. These simulations are part of a multi-faceted 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 line-of-sight 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=6-12$), 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 large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m 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 delay-spectrum technique to motivate high-level 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, current-generation 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.

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. Next-generation telescopes are now being designed and built to meet these challenges, drawing lessons from first-generation experiments that showed the benefits of densely packed, highly redundant arrays--in which the same mode on the sky is sampled by many antenna pairs--for 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 sub-aperture sampling via "off-grid" antennas and increased angular resolution via far-flung "outrigger" antennas is possible with a redundantly calibratable array configuration.

When planet-hosting stars evolve off the main sequence and go through the red-giant branch, the stars become orders of magnitudes more luminous and, at the same time, lose mass at much higher rates than their main-sequence counterparts. Accordingly, if planetary companions exist around these stars at orbital distances of a few AU, they will be heated up to the level of canonical hot Jupiters and also be subjected to a dense stellar wind. Given that magnetized planets interacting with stellar winds emit radio waves, such "Red-Giant Hot Jupiters" (RGHJs) may also be candidate radio emitters. We estimate the spectral auroral radio intensity of RGHJs based on the empirical relation with the stellar wind as well as a proposed scaling for planetary magnetic fields. RGHJs might be intrinsically as bright as or brighter than canonical hot Jupiters and about 100 times brighter than equivalent objects around main-sequence stars. We examine the capabilities of low-frequency radio observatories to detect this emission and find that the signal from an RGHJ may be detectable at distances up to a few hundred parsecs with the Square Kilometer Array.

Current-generation low frequency interferometers constructed with the objective of detecting the high-redshift 21 cm background, aim to generate power spectra of the brightness-temperature contrast of neutral hydrogen in primordial intergalactic medium. Two-dimensional power spectra (power in Fourier modes parallel and perpendicular to the line of sight) formed from interferometric visibilities have been shown to delineate a boundary between spectrally-smooth foregrounds (known as the wedge) and spectrally-structured 21 cm background emission (the EoR-window). However, polarized foregrounds are known to possess spectral structure due to Faraday rotation, which can leak into the EoR window. In this work, we create and analyze 2D power spectra from the PAPER-32 imaging array in Stokes I, Q, U and V. These allow us to observe and diagnose systematic effects in our calibration at high signal-to-noise within the Fourier space most relevant to EoR experiments. We observe well-defined windows in the Stokes visibilities, with Stokes Q, U and V power spectra sharing a similar wedge shape to that seen in Stokes I. With modest polarization calibration, we see no evidence that polarization calibration errors move power outside the wedge in any Stokes visibility, to the noise levels attained. Deeper integrations will be required to confirm that this behavior persists to the depth required for EoR detection.

Unaccounted for systematics from foregrounds and instruments can severely limit the sensitivity of current experiments from detecting redshifted 21~cm signals from the Epoch of Reionization (EoR). Upcoming experiments are faced with a challenge to deliver more collecting area per antenna element without degrading the data with systematics. This paper and its companions show that dishes are viable for achieving this balance using the Hydrogen Epoch of Reionization Array (HERA) as an example. Here, we specifically identify spectral systematics associated with the antenna power pattern as a significant detriment to all EoR experiments which causes the already bright foreground power to leak well beyond ideal limits and contaminate the otherwise clean EoR signal modes. A primary source of this chromaticity is reflections in the antenna-feed assembly and between structures in neighboring antennas. Using precise foreground simulations taking wide-field effects into account, we provide a framework to set cosmologically-motivated design specifications on these reflections to prevent further EoR signal degradation. We show HERA will not be impeded by such spectral systematics and demonstrate that even in a conservative scenario that does not perform removal of foregrounds, HERA will detect EoR signal in line-of-sight $k$-modes, $k_\parallel \gtrsim 0.2\,h$~Mpc$^{-1}$, with high significance. All baselines in a 19-element HERA layout are capable of detecting EoR over a substantial observing window on the sky.

We generalize the technique of fringe-rate filtering, whereby visibilities measured by a radio interferometer are re-weighted according to their temporal variation. As the Earth rotates, radio sources traverse through an interferometer's fringe pattern at rates that depend on their position on the sky. Capitalizing on this geometric interpretation of fringe rates, we employ time-domain convolution kernels to enact fringe-rate filters that sculpt the effective primary beam of antennas in an interferometer. As we show, beam sculpting through fringe-rate filtering can be used to optimize measurements for a variety of applications, including mapmaking, minimizing polarization leakage, suppressing instrumental systematics, and enhancing the sensitivity of power-spectrum measurements. We show that fringe-rate filtering arises naturally in minimum variance treatments of many of these problems, enabling optimal visibility-based approaches to analyses of interferometric data that avoid systematics potentially introduced by traditional approaches such as imaging. Our techniques have recently been demonstrated in Ali et al. (2015), where new upper limits were placed on the 21 cm power spectrum from reionization, showcasing the ability of fringe-rate filtering to successfully boost sensitivity and reduce the impact of systematics in deep observations.

Amongst standard model parameters that are constrained by cosmic microwave background (CMB) observations, the optical depth $\tau$ stands out as a nuisance parameter. While $\tau$ provides some crude limits on reionization, it also degrades constraints on other cosmological parameters. Here we explore how 21 cm cosmology---as a direct probe of reionization---can be used to independently predict $\tau$ in an effort to improve CMB parameter constraints. We develop two complementary schemes for doing so. The first uses 21 cm power spectrum observations in conjunction with semi-analytic simulations to predict $\tau$. The other uses global 21 cm measurements to directly constrain low redshift (post-reheating) contributions to $\tau$ in a relatively model-independent way. Forecasting the performance of the upcoming Hydrogen Epoch of Reionization Array, we find that significant reductions in the errors on $\tau$ can be achieved. These results are particularly effective at breaking the CMB degeneracy between $\tau$ and the amplitude of the primordial fluctuation spectrum $A_s$, with errors on $\ln (10^{10} A_s)$ reduced by up to a factor of four. Stage 4 CMB constraints on the neutrino mass sum are also improved, with errors potentially reduced to $12\,\textrm{meV}$ regardless of whether CMB experiments can precisely measure the reionization bump in polarization power spectra. Observations of the 21 cm line are therefore capable of improving not only our understanding of reionization astrophysics, but also of cosmology in general.

Improvements in current instruments and the advent of next-generation instruments will soon push observational 21 cm cosmology into a new era, with high significance measurements of both the power spectrum and the mean ("global") signal of the 21 cm brightness temperature. In this paper we use the recently commenced Hydrogen Epoch of Reionization Array as a worked example to provide forecasts on astrophysical and cosmological parameter constraints. In doing so we improve upon previous forecasts in a number of ways. First, we provide updated forecasts using the latest best-fit cosmological parameters from the Planck satellite, exploring the impact of different Planck datasets on 21 cm experiments. We also show that despite the exquisite constraints that other probes have placed on cosmological parameters, the remaining uncertainties are still large enough to have a non-negligible impact on upcoming 21 cm data analyses. While this complicates high-precision constraints on reionization models, it provides an avenue for 21 cm reionization measurements to constrain cosmology. We additionally forecast HERA's ability to measure the ionization history using a combination of power spectrum measurements and semi-analytic simulations. Finally, we consider ways in which 21 cm global signal and power spectrum measurements can be combined, and propose a method by which power spectrum results can be used to train a compact parameterization of the global signal. This parameterization reduces the number of parameters needed to describe the global signal, increasing the likelihood of a high significance measurement.

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating.

A measurement of the cosmological 21 cm signal remains a promising but as-of-yet unattained ambition of radio astronomy. A positive detection would provide direct observations of key unexplored epochs of our cosmic history, including the cosmic dark ages and reionization. In this paper, we concentrate on measurements of the spatial monopole of the 21 cm brightness temperature as a function of redshift (the "global signal"). Most global experiments to date have been single-element experiments. In this paper, we show how an interferometer can be designed to be sensitive to the monopole mode of the sky, thus providing an alternate approach to accessing the global signature. We provide simple rules of thumb for designing a global signal interferometer and use numerical simulations to show that a modest array of tightly packed antenna elements with moderately sized primary beams (full-width-half-max of $\sim$40$^\circ$) can compete with typical single-element experiments in their ability to constrain phenomenological parameters pertaining to reionization and the pre-reionization era. We also provide a general data analysis framework for extracting the global signal from interferometric measurements (with analysis of single-element experiments arising as a special case) and discuss trade-offs with various data analysis choices. Given that interferometric measurements are able to avoid a number of systematics inherent in single-element experiments, our results suggest that interferometry ought to be explored as a complementary way to probe the global signal.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

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 foreground-free 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 well-controlled 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 trade-offs between accuracy and speed using one such telescope, the upcoming Hydrogen Epoch of Reionization Array, as a case study.

Strontium iodide doped with europium is a new scintillator material being developed as an alternative to lanthanum bromide doped with cerium for use in high energy astrophysical detectors. As with all scintillators, the issue of nonproportionality is important because it affects the energy resolution of the detector. In this study, we investigate how the nonproportionality of strontium iodide doped with europium changes as a function of temperature 16 deg. C to 60 deg. C by heating the strontium iodide doped with europium scintillator separate from the photomultiplier tube. In a separate experiment, we also investigate the nonproportionality at high energies (up to 6 MeV) of strontium iodide doped with europium at a testing facility located at NASA Goddard Space Flight Center. We find that the nonproportionality increases nearly monotonically as the temperature of the strontium iodide doped with europium scintillator is increased, although there is evidence of non-monotonic behavior near 40 deg. C, perhaps due to electric charge carriers trapping in the material. We also find that within the energy range of 662 keV to 6.1 MeV, the change in the nonproportionality of the strontium iodide doped with europium is about 1.5 to 2 percent.

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 next-generation 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 k-space. 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.

The epoch of reionization power spectrum is expected to evolve strongly with redshift, and it is this variation with cosmic history that will allow us to begin to place constraints on the physics of reionization. The primary obstacle to the measurement of the EoR power spectrum is bright foreground emission. We present an analysis of observations from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) telescope which place new limits on the HI power spectrum over the redshift range of $7.5<z<10.5$, extending previously published single redshift results to cover the full range accessible to the instrument. To suppress foregrounds, we use filtering techniques that take advantage of the large instrumental bandwidth to isolate and suppress foreground leakage into the interesting regions of $k$-space. Our 500 hour integration is the longest such yet recorded and demonstrates this method to a dynamic range of $10^4$. Power spectra at different points across the redshift range reveal the variable efficacy of the foreground isolation. Noise limited measurements of $\Delta^2$ at $k=$0.2hMpc$^{-1}$ and z$=7.55$ reach as low as (48mK)$^2$ ($1\sigma$). We demonstrate that the size of the error bars in our power spectrum measurement as generated by a bootstrap method is consistent with the fluctuations due to thermal noise. Relative to this thermal noise, most spectra exhibit an excess of power at a few sigma. The likely sources of this excess include residual foreground leakage, particularly at the highest redshift, and unflagged RFI. We conclude by discussing data reduction improvements that promise to remove much of this excess.

The 21 cm line provides a powerful probe of astrophysics and cosmology at high redshifts, but unlocking the potential of this probe requires the robust mitigation of foreground contaminants that are typically several orders of magnitude brighter than the cosmological signal. Recent simulations and observations have shown that the smooth spectral structure of foregrounds combines with instrument chromaticity to contaminate a "wedge"-shaped region in cylindrical Fourier space. While previous efforts have explored the suppression of foregrounds within this wedge, as well as the avoidance of this highly contaminated region, all such efforts have neglected a rigorous examination of the error statistics associated with the wedge. Using a quadratic estimator formalism applied to the interferometric measurement equation, we provide a framework for such a rigorous analysis (incorporating a fully covariant treatment of errors). Additionally, we find that there are strong error correlations at high spatial wavenumbers that have so far been neglected in sensitivity derivations. These error correlations substantially degrade the sensitivity of arrays relying on contributions from long baselines, compared to what one would estimate assuming uncorrelated errors.

For there to be a successful measurement of the 21 cm Epoch of Reionization (EoR) power spectrum, it is crucial that strong foreground contaminants be robustly suppressed. These foregrounds come from a variety of sources (such as Galactic synchrotron emission and extragalactic point sources), but almost all share the property of being spectrally smooth, and when viewed through the chromatic response of an interferometer, occupy a signature "wedge" region in cylindrical $k_\perp k_\parallel$ Fourier space. The complement of the foreground wedge is termed the "EoR window", and is expected to be mostly foreground-free, allowing clean measurements of the power spectrum. This paper is a sequel to a previous paper that established a rigorous mathematical framework for describing the foreground wedge and the EoR window. Here, we use our framework to explore statistical methods by which the EoR window can be enlarged, thereby increasing the sensitivity of a power spectrum measurement. We adapt the FKP approximation (commonly used in galaxy surveys) for 21 cm cosmology, and also compare the optimal quadratic estimator to simpler estimators that ignore covariances between different Fourier modes. The optimal quadratic estimator is found to suppress foregrounds by an extra factor of $\sim 10^5$ in power at the peripheries of the EoR window, boosting the detection of the cosmological signal from $12\sigma$ to $50\sigma$ at the midpoint of reionization in our fiducial models. If numerical issues can be finessed, decorrelation techniques allow the EoR window to be further enlarged, enabling measurements to be made deep within the foreground wedge. These techniques do not assume that foreground are Gaussian-distributed, and we additionally prove that a final round of foreground subtraction can be performed after decorrelation in a way that is guaranteed to have no cosmological signal loss.

We present new constraints on the 21cm Epoch of Reionization (EoR) power spectrum derived from 3 months of observing with a 32-antenna, dual-polarization deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. In this paper, we demonstrate the efficacy of the delay-spectrum approach to avoiding foregrounds, achieving over 8 orders of magnitude of foreground suppression (in $\textrm{mK}^2$). Combining this approach with a procedure for removing off-diagonal covariances arising from instrumental systematics, we achieve a best 2-sigma upper limit of $(41\,\textrm{mK})^2$ for $k=0.27 h\textrm{Mpc}^{-1}$ at $z=7.7$. This limit falls within an order of magnitude of the brighter predictions of the expected 21cm EoR signal level. Using the upper limits set by these measurements, we generate new constraints on the brightness temperature of 21cm emission in neutral regions for various reionization models. We show that for several ionization scenarios, our measurements are inconsistent with cold reionization. That is, heating of the neutral intergalactic medium (IGM) is necessary to remain consistent with the constraints we report. Hence, we have suggestive evidence that by $z=7.7$, the HI has been warmed from its cold primordial state, probably by X-rays from high-mass X-ray binaries or mini-quasars. The strength of this evidence depends on the ionization state of the IGM, which we are not yet able to constrain. This result is consistent with standard predictions for how reionization might have proceeded.

We present a catalog of spectral measurements covering a 100-200 MHz band for 32 sources, derived from observations with a 64-antenna deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. For transit telescopes such as PAPER, calibration of the primary beam is a difficult endeavor, and errors in this calibration are a major source of error in the determination of source spectra. In order to decrease reliance on accurate beam calibration, we focus on calibrating sources in a narrow declination range from -46d to -40d. Since sources at similar declinations follow nearly identical paths through the primary beam, this restriction greatly reduces errors associated with beam calibration, yielding a dramatic improvement in the accuracy of derived source spectra. Extrapolating from higher frequency catalogs, we derive the flux scale using a Monte-Carlo fit across multiple sources that includes uncertainty from both catalog and measurement errors. Fitting spectral models to catalog data and these new PAPER measurements, we derive new flux models for Pictor A and 31 other sources at nearby declinations. 90% of these confirm and refine a power-law model for flux density. Of note is the new Pictor A flux model, which is accurate to 1.4% and shows, in contrast to previous models, that between 100 MHz and 2 GHz, the spectrum of Pictor A is consistent with a single power law given by a flux at 150 MHz of 382+/-5.4 Jy, and a spectral index of -0.76+/-0.01. This accuracy represents an order of magnitude improvement over previous measurements in this band, and is limited by the uncertainty in the catalog measurements used to estimate the absolute flux scale. The simplicity and improved accuracy of Pictor A's spectrum make it an excellent calibrator for experiments seeking to measure 21cm emission from the Epoch of Reionization.

Experiments aimed at detecting highly-redshifted 21 centimeter emission from the Epoch of Reionization (EoR) are plagued by the contamination of foreground emission. A potentially important source of contaminating foregrounds may be Faraday-rotated, polarized emission, which leaks into the estimate of the intrinsically unpolarized EoR signal. While these foregrounds' intrinsic polarization may not be problematic, the spectral structure introduced by the Faraday rotation could be. To better understand and characterize these effects, we present a simulation of the polarized sky between 120 and 180 MHz. We compute a single visibility, and estimate the three-dimensional power spectrum from that visibility using the delay spectrum approach presented in Parsons et al. (2012b) . Using the Donald C. Backer Precision Array to Probe the Epoch of Reionization (PAPER) as an example instrument, we show the expected leakage into the unpolarized power spectrum to be several orders of magnitude above the expected 21cm EoR signal.

We present new observations with the Precision Array for Probing the Epoch of Reionization (PAPER) with the aim of measuring the properties of foreground emission for 21cm Epoch of Reionization experiments at 150 MHz. We focus on the footprint of the foregrounds in cosmological Fourier space to understand which modes of the 21cm power spectrum will most likely be compromised by foreground emission. These observations confirm predictions that foregrounds can be isolated to a "wedge"-like region of 2D (k-perpendicular, k-parallel)-space, creating a window for cosmological studies at higher k-parallel values. We also find that the emission extends past the nominal edge of this wedge due to spectral structure in the foregrounds, with this feature most prominent on the shortest baselines. Finally, we filter the data to retain only this "unsmooth" emission and image specific k-parallel modes of it. The resultant images show an excess of power at the lowest modes, but no emission can be clearly localized to any one region of the sky. This image is highly suggestive that the most problematic foregrounds for 21cm EoR studies will not be easily identifiable bright sources, but rather an aggregate of fainter emission.

We present observations taken with the Precision Array for Probing the Epoch of Reionization (PAPER) of the Centaurus A field in the frequency range 114 to 188 MHz. The resulting image has a 25' resolution, a dynamic range of 3500 and an r.m.s. of 0.5 Jy\beam (for a beam size of 25' x 23'). A spectral index map of Cen A is produced across the full band. The spectral index distribution is qualitatively consistent with electron reacceleration in regions of excess turbulence in the radio lobes, as previously identified morphologically. Hence, there appears to be an association of 'severe weather' in radio lobes with energy input into the relativistic electron population. We compare the PAPER large scale radio image with the X-ray image from the ROSAT All Sky Survey. There is a tentative correlation between radio and X-ray features at the end of the southern lobe, some 200 kpc from the nucleus, as might be expected from inverse Compton scattering of the CMB by the relativistic electrons also responsible for the radio synchrotron emission. The magnetic fields derived from the (possible) IC and radio emission are of similar magnitude to fields derived under the minimum pressure assumptions, ~ 1 {\mu}G. However, the X-ray field is complex, with large scale gradients and features possibly unrelated to Cen A. If these X-ray features are unrelated to Cen A, then these fields are lower limits.

This work describes a new instrument optimized for a detection of the neutral hydrogen 21cm power spectrum between redshifts of 0.5-1.5: the Baryon Acoustic Oscillation Broadband and Broad-beam (BAOBAB) Array. BAOBAB will build on the efforts of a first generation of 21cm experiments which are targeting a detection of the signal from the Epoch of Reionization at z ~ 10. At z ~ 1, the emission from neutral hydrogen in self-shielded overdense halos also presents an accessible signal, since the dominant, synchrotron foreground emission is considerably fainter than at redshift 10. The principle science driver for these observations are Baryon Acoustic Oscillations in the matter power spectrum which have the potential to act as a standard ruler and constrain the nature of dark energy. BAOBAB will fully correlate dual-polarization antenna tiles over the 600-900MHz band with a frequency resolution of 300 kHz and a system temperature of 50K. The number of antennas will grow in staged deployments, and reconfigurations of the array will allow for both traditional imaging and high power spectrum sensitivity operations. We present calculations of the power spectrum sensitivity for various array sizes, with a 35-element array measuring the cosmic neutral hydrogen fraction as a function of redshift, and a 132-element system detecting the BAO features in the power spectrum, yielding a 1.8% error on the z ~ 1 distance scale, and, in turn, significant improvements to constraints on the dark energy equation of state over an unprecedented range of redshifts from ~0.5-1.5.

A critical challenge in measuring the power spectrum of 21cm emission from cosmic reionization is compensating for the frequency dependence of an interferometer's sampling pattern, which can cause smooth-spectrum foregrounds to appear unsmooth and degrade the separation between foregrounds and the target signal. In this paper, we present an approach to foreground removal that explicitly accounts for this frequency dependence. We apply the delay transformation introduced in Parsons & Backer (2009) to each baseline of an interferometer to concentrate smooth-spectrum foregrounds within the bounds of the maximum geometric delays physically realizable on that baseline. By focusing on delay-modes that correspond to image-domain regions beyond the horizon, we show that it is possible to avoid the bulk of smooth-spectrum foregrounds. We map the point-spread function of delay-modes to k-space, showing that delay-modes that are uncorrupted by foregrounds also represent samples of the three-dimensional power spectrum, and can be used to constrain cosmic reionization. Because it uses only spectral smoothness to differentiate foregrounds from the targeted 21cm signature, this per-baseline analysis approach relies on spectrally- and spatially-smooth instrumental responses for foreground removal. For sufficient levels of instrumental smoothness relative to interfering foregrounds, this technique substantially reduces the level of calibration previously thought necessary to detect 21cm reionization. As a result, this approach places fewer constraints on antenna configuration within an array, facilitating the adoption of configurations optimized for power-spectrum sensitivity. Under these assumptions, we demonstrate the potential for PAPER to detect 21cm reionization at an amplitude of 10 mK^2 near k~0.2h Mpc^-1 with 132 dipoles in 7 months of observing.

Telescopes aiming to measure 21cm emission from the Epoch of Reionization must toe a careful line, balancing the need for raw sensitivity against the stringent calibration requirements for removing bright foregrounds. It is unclear what the optimal design is for achieving both of these goals. Via a pedagogical derivation of an interferometer's response to the power spectrum of 21cm reionization fluctuations, we show that even under optimistic scenarios, first-generation arrays will yield low-SNR detections, and that different compact array configurations can substantially alter sensitivity. We explore the sensitivity gains of array configurations that yield high redundancy in the uv-plane -- configurations that have been largely ignored since the advent of self-calibration for high-dynamic-range imaging. We first introduce a mathematical framework to generate optimal minimum-redundancy configurations for imaging. We contrast the sensitivity of such configurations with high-redundancy configurations, finding that high-redundancy configurations can improve power-spectrum sensitivity by more than an order of magnitude. We explore how high-redundancy array configurations can be tuned to various angular scales, enabling array sensitivity to be directed away from regions of the uv-plane (such as the origin) where foregrounds are brighter and where instrumental systematics are more problematic. We demonstrate that a 132-antenna deployment of the Precision Array for Probing the Epoch of Reionization (PAPER) observing for 120 days in a high-redundancy configuration will, under ideal conditions, have the requisite sensitivity to detect the power spectrum of the 21cm signal from reionization at a 3\sigma level at k<0.25h Mpc^{-1} in a bin of \Delta ln k=1. We discuss the tradeoffs of low- versus high-redundancy configurations.

We present a new technique for calibrating the primary beam of a wide-field, drift-scanning antenna element. Drift-scan observing is not compatible with standard beam calibration routines, and the situation is further complicated by difficult-to-parametrize beam shapes and, at low frequencies, the sparsity of accurate source spectra to use as calibrators. We overcome these challenges by building up an interrelated network of source "crossing points" -- locations where the primary beam is sampled by multiple sources. Using the single assumption that a beam has 180 degree rotational symmetry, we can achieve significant beam coverage with only a few tens of sources. The resulting network of crossing points allows us to solve for both a beam model and source flux densities referenced to a single calibrator source, circumventing the need for a large sample of well-characterized calibrators. We illustrate the method with actual and simulated observations from the Precision Array for Probing the Epoch of Reionization (PAPER).

Here we describe our ongoing efforts to develop high-performance and sensitive instrumentation for use in the search for extra-terrestrial intelligence (SETI). These efforts include our recently deployed Search for Extraterrestrial Emissions from Nearby Developed Intelligent Populations Spectrometer (SERENDIP V.v) and two instruments currently under development; the Heterogeneous Radio SETI Spectrometer (HRSS) for SETI observations in the radio spectrum and the Optical SETI Fast Photometer (OSFP) for SETI observations in the optical band. We will discuss the basic SERENDIP V.v instrument design and initial analysis methodology, along with instrument architectures and observation strategies for OSFP and HRSS. In addition, we will demonstrate how these instruments may be built using low-cost, modular components and programmed and operated by students using common languages, e.g. ANSI C.

We present observations from the Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa, observed in May and September 2010. Using two nights of drift scanning with PAPER's 60\arcdeg\ FWHM beam we have made a map covering the entire sky below +10 degrees declination with an effective center frequency of 145 MHz, a 70-MHz bandwidth, and a resolution of 26\arcmin. A 4800 square-degree region of this large map with the lowest Galactic emission reaches an RMS of 0.7 Jy. We establish an absolute flux scale using sources from the 160-MHz Culgoora catalog. Using the 408-MHz Molonglo Reference Catalog (MRC) as a finding survey, we identify counterparts to 480 sources in our maps, and compare our fluxes to the MRC and to 332 sources in the Culgoora catalog. For both catalogs, the ratio of PAPER to catalog flux averages to 1, with a standard deviation of 50%. This measured variation is consistent with comparisons between independent catalogs observed at different bands. The PAPER data represent new 145-MHz flux measurements for a large number of sources in the band expected to encompass cosmic reionization, and represents a significant step toward establishing a model for removing foregrounds to the reionization signal.

We present the second Swift Burst Alert Telescope (BAT) catalog of gamma-ray bursts (GRBs), which contains 476 bursts detected by the BAT between 2004 December 19 and 2009 December 21. This catalog (hereafter the BAT2 catalog) presents burst trigger time, location, 90% error radius, duration, fluence, peak flux, time-averaged spectral parameters and time-resolved spectral parameters measured by the BAT. In the correlation study of various observed parameters extracted from the BAT prompt emission data, we distinguish among long-duration GRBs (L-GRBs), short-duration GRBs (S-GRBs), and short-duration GRBs with extended emission (S-GRBs with E.E.) to investigate differences in the prompt emission properties. The fraction of L-GRBs, S-GRBs and S-GRBs with E.E. in the catalog are 89%, 8% and 2% respectively. We compare the BAT prompt emission properties with the BATSE, BeppoSAX and HETE-2 GRB samples. We also correlate the observed prompt emission properties with the redshifts for the GRBs with known redshift. The BAT T90 and T50 durations peak at 70 s and 30 s, respectively. We confirm that the spectra of the BAT S-GRBs are generally harder than those of the L-GRBs. The time-averaged spectra of the BAT S-GRBs with E.E. are similar to those of the L-GRBs. Whereas, the spectra of the initial short spikes of the S-GRBs with E.E. are similar to those of the S-GRBs. We show that the BAT GRB samples are significantly softer than the BATSE bright GRBs, and that the time-averaged Epeak of the BAT GRBs peaks at 80 keV which is significantly lower energy than those of the BATSE sample which peak at 320 keV. The time-averaged spectral properties of the BAT GRB sample are similar to those of the HETE-2 GRB samples. By time-resolved spectral analysis, we find that only 10% of the BAT observed photon indices are outside the allowed region of the synchrotron shock model.

We report on the spectral cross-calibration results of the Konus-Wind, the Suzaku/WAM, and the Swift/BAT instruments using simultaneously observed gamma-ray bursts (GRBs). This is the first attempt to use simultaneously observed GRBs as a spectral calibration source to understand systematic problems among the instruments. Based on these joint spectral fits, we find that 1) although a constant factor (a normalization factor) agrees within 20% among the instruments, the BAT constant factor shows a systematically smaller value by 10-20% compared to that of Konus-Wind, 2) there is a systematic trend that the low-energy photon index becomes steeper by 0.1-0.2 and Epeak becomes systematically higher by 10-20% when including the BAT data in the joint fits, and 3) the high-energy photon index agrees within 0.2 among the instruments. Our results show that cross-calibration based on joint spectral analysis is an important step to understanding the instrumental effects which could be affecting the scientific results from the GRB prompt emission data.

We present results from Swift, XMM-Newton, and deep INTEGRAL monitoring in the region of GRB 050925. This short Swift burst is a candidate for a newly discovered soft gamma-ray repeater (SGR) with the following observational burst properties: 1) galactic plane (b=-0.1 deg) localization, 2) 150 msec duration, and 3) a blackbody rather than a simple power-law spectral shape (with a significance level of 97%). We found two possible X-ray counterparts of GRB 050925 by comparing the X-ray images from Swift XRT and XMM-Newton. Both X-ray sources show the transient behavior with a power-law decay index shallower than -1. We found no hard X-ray emission nor any additional burst from the location of GRB 050925 in ~5 Ms of INTEGRAL data. We discuss about the three BATSE short bursts which might be associated with GRB 050925, based on their location and the duration. Assuming GRB 050925 is associated with the H II regions (W 58) at the galactic longitude of l=70 deg, we also discuss the source frame properties of GRB 050925.

We present a filtering technique that can be applied to individual baselines of wide-bandwidth, wide-field interferometric data to geometrically select regions on the celestial sphere that contain primary calibration sources. The technique relies on the Fourier transformation of wide-band frequency spectra from a given baseline to obtain one-dimensional "delay images", and then the transformation of a time-series of delay images to obtain two-dimensional "delay/delay-rate images." Source selection is possible in these images given appropriate combinations of baseline, bandwidth, integration time and source location. Strong and persistent radio frequency interference (RFI) limits the effectiveness of this source selection owing to the removal of data by RFI excision algorithms. A one-dimensional, complex CLEAN algorithm has been developed to compensate for RFI-excision effects. This approach allows CLEANed, source-isolated data to be used to isolate bandpass and primary beam gain functions. These techniques are applied to data from the Precision Array for Probing the Epoch of Reionization (PAPER) as a demonstration of their value in calibrating a new generation of low-frequency radio interferometers with wide relative bandwidths and large fields-of-view.

We present the catalog of sources detected in the first 22 months of data from the hard X-ray survey (14--195 keV) conducted with the BAT coded mask imager on the \swift satellite. The catalog contains 461 sources detected above the 4.8 sigma level with BAT. High angular resolution X-ray data for every source from Swift XRT or archival data have allowed associations to be made with known counterparts in other wavelength bands for over 97% of the detections, including the discovery of ~30 galaxies previously unknown as AGN and several new Galactic sources. A total of 266 of the sources are associated with Seyfert galaxies (median redshift z ~ 0.03) or blazars, with the majority of the remaining sources associated with X-ray binaries in our Galaxy. This ongoing survey is the first uniform all sky hard X-ray survey since HEAO-1 in 1977. Since the publication of the 9-month BAT survey we have increased the number of energy channels from 4 to 8 and have substantially increased the number of sources with accurate average spectra. The BAT 22-month catalog is the product of the most sensitive all-sky survey in the hard X-ray band, with a detection sensitivity (4.8 sigma) of 2.2e-11 erg/cm2/s (1 mCrab) over most of the sky in the 14--195 keV band.

We are developing the Precision Array for Probing the Epoch of Reionization (PAPER) to detect 21cm emission from the early Universe, when the first stars and galaxies were forming. We describe the overall experiment strategy and architecture and summarize two PAPER deployments: a 4-antenna array in the low-RFI environment of Western Australia and an 8-antenna array at our prototyping site in Green Bank, WV. From these activities we report on system performance, including primary beam model verification, dependence of system gain on ambient temperature, measurements of receiver and overall system temperatures, and characterization of the RFI environment at each deployment site. We present an all-sky map synthesized between 139 MHz and 174 MHz using data from both arrays that reaches down to 80 mJy (4.9 K, for a beam size of 2.15e-5 steradians at 154 MHz), with a 10 mJy (620 mK) thermal noise level that indicates what would be achievable with better foreground subtraction. We calculate angular power spectra ($C_\ell$) in a cold patch and determine them to be dominated by point sources, but with contributions from galactic synchrotron emission at lower radio frequencies and angular wavemodes. Although the cosmic variance of foregrounds dominates errors in these power spectra, we measure a thermal noise level of 310 mK at $\ell=100$ for a 1.46-MHz band centered at 164.5 MHz. This sensitivity level is approximately three orders of magnitude in temperature above the level of the fluctuations in 21cm emission associated with reionization.

Time-to-science is an important figure of merit for digital instrumentation serving the astronomical community. A digital signal processing (DSP) community is forming that uses shared hardware development, signal processing libraries, and instrument architectures to reduce development time of digital instrumentation and to improve time-to-science for a wide variety of projects. We suggest prioritizing technological development supporting the needs of this nascent DSP community. After outlining several instrument classes that are relying on digital instrumentation development to achieve new science objectives, we identify key areas where technologies pertaining to interoperability and processing flexibility will reduce the time, risk, and cost of developing the digital instrumentation for radio astronomy. These areas represent focus points where support of general-purpose, open-source development for a DSP community should be prioritized in the next decade. Contributors to such technological development may be centers of support for this DSP community, science groups that contribute general-purpose DSP solutions as part of their own instrumentation needs, or engineering groups engaging in research that may be applied to next-generation DSP instrumentation.

A new generation of radio telescopes is achieving unprecedented levels of sensitivity and resolution, as well as increased agility and field-of-view, by employing high-performance digital signal processing hardware to phase and correlate large numbers of antennas. The computational demands of these imaging systems scale in proportion to BMN^2, where B is the signal bandwidth, M is the number of independent beams, and N is the number of antennas. The specifications of many new arrays lead to demands in excess of tens of PetaOps per second. To meet this challenge, we have developed a general purpose correlator architecture using standard 10-Gbit Ethernet switches to pass data between flexible hardware modules containing Field Programmable Gate Array (FPGA) chips. These chips are programmed using open-source signal processing libraries we have developed to be flexible, scalable, and chip-independent. This work reduces the time and cost of implementing a wide range of signal processing systems, with correlators foremost among them,and facilitates upgrading to new generations of processing technology. We present several correlator deployments, including a 16-antenna, 200-MHz bandwidth, 4-bit, full Stokes parameter application deployed on the Precision Array for Probing the Epoch of Reionization.

We present the earliest ever ultraviolet spectrum of a gamma-ray burst (GRB) as observed with the Swift-UVOT. The GRB 081203A spectrum was observed for 50 seconds with the UV grism starting 251 seconds after the Swift-BAT trigger when the GRB was of u ~13.4 mag and still rising to its peak optical brightness. The UV grism spectrum shows a damped Ly-alpha line, Ly-beta, and the Lyman continuum break at a redshift z = 2.05 +/- 0.01. A model fit to the Lyman absorption implies log N(HI) = 22.0 +/- 0.2 cm-2, which is typical for GRB host galaxies with damped Ly-alpha absorbers. This observation of GRB 081203A demonstrates that for GRBs brighter than v ~14 mag and with 0.5 < z < 3.5 the UVOT will be able to provide redshifts, and probe for damped Ly-alpha absorbers within 4-6 minutes from the time of the Swift-BAT trigger.

We report a correlation based on a spectral simulation study of the prompt emission spectra of gamma-ray bursts (GRBs) detected by the Swift Burst Alert Telescope (BAT). The correlation is between the Epeak energy, which is the peak energy in the \nu F_\nu spectrum, and the photon index (\Gamma) derived from a simple power-law model. The Epeak - \Gamma relation, assuming the typical smoothly broken power-law spectrum of GRBs, is \log Epeak = 3.258 - 0.829\Gamma (1.3 < \Gamma < 2.3). We take into account not only a range of Epeak energies and fluences, but also distributions for both the low-energy photon index and the high-energy photon index in the smoothly broken power-law model. The distribution of burst durations in the BAT GRB sample is also included in the simulation. Our correlation is consistent with the index observed by BAT and Epeak measured by the BAT, and by other GRB instruments. Since about 85% of GRBs observed by the BAT are acceptably fit with the simple power-law model because of the relatively narrow energy range of the BAT, this relationship can be used to estimate Epeak when it is located within the BAT energy range.

We present the results from an ESO/VLT campaign aimed at studying the afterglow properties of the short/hard gamma ray burst GRB 070707. Observations were carried out at ten different epochs from ~0.5 to ~80 days after the event. The optical flux decayed steeply with a power-law decay index greater than 3, later levelling off at R~27.3 mag; this is likely the emission level of the host galaxy, the faintest yet detected for a short GRB. Spectroscopic observations did not reveal any line features/edges that could unambiguously pinpoint the GRB redshift, but set a limit z < 3.6. In the range of allowed redshifts, the host has a low luminosity, comparable to that of long-duration GRBs. The existence of such faint host galaxies suggests caution when associating short GRBs with bright, offset galaxies, where the true host might just be too dim for detection. The steepness of the decay of the optical afterglow of GRB 070707 challenges external shock models for the optical afterglow of short/hard GRBs. We argue that this behaviour might results from prolonged activity of the central engine or require alternative scenarios.

We describe and discuss the spectral and temporal characteristics of the prompt emission and X-ray afterglow emission of X-ray flashes (XRFs) and X-ray-rich gamma-ray bursts (XRRs) detected and observed by Swift between December 2004 and September 2006. We compare these characteristics to a sample of conventional classical gamma-ray bursts (C-GRBs) observed during the same period. We confirm the correlation between Epeak_obs and fluence noted by others and find further evidence that XRFs, XRRs and C-GRBs form a continuum. We also confirm that our known redshift sample is consistent with the correlation between the peak energy in the GRB rest frame (Epeak_src) and the isotropic radiated energy (Eiso), so called the Epeak_src-Eiso relation. The spectral properties of X-ray afterglows of XRFs and C-GRBs are similar, but the temporal properties of XRFs and C-GRBs are quite different. We found that the light curves of C-GRB afterglows show a break to steeper indices (shallow-to-steep break) at much earlier times than do XRF afterglows. Moreover, the overall luminosity of XRF X-ray afterglows is systematically smaller by a factor of two or more compared to that of C-GRBs. These distinct differences between the X-ray afterglows of XRFs and C-GRBs may be the key to understanding not only the mysterious shallow-to-steep break in X-ray afterglow light curves, but also the unique nature of XRFs.

We present the first Swift Burst Alert Telescope (BAT) catalog of gamma-ray bursts (GRBs), which contains bursts detected by the BAT between 2004 December 19 and 2007 June 16. This catalog (hereafter BAT1 catalog) contains burst trigger time, location, 90% error radius, duration, fluence, peak flux, and time averaged spectral parameters for each of 237 GRBs, as measured by the BAT. The BAT-determined position reported here is within 1.75' of the Swift X-ray Telescope (XRT)-determined position for 90% of these GRBs. The BAT T_90 and T_50 durations peak at 80 and 20 seconds, respectively. From the fluence-fluence correlation, we conclude that about 60% of the observed peak energies, Epeak, of BAT GRBs could be less than 100 keV. We confirm that GRB fluence to hardness and GRB peak flux to hardness are correlated for BAT bursts in analogous ways to previous missions' results. The correlation between the photon index in a simple power-law model and Epeak is also confirmed. We also report the current status for the on-orbit BAT calibrations based on observations of the Crab Nebula.

Almost the totality of the bright foreground sources in the WMAP CMB maps are blazars, a class of sources that show usually also X-ray emission. However, 23 objects in a flux-limited sample of 140 blazars of the WMAP catalog (first year) were never reported before as X-ray sources. We present here the results of 41 Swift observations which led to the detection of all these 23 blazars in the 0.3-10 keV band. We conclude that all micro-wave selected blazars are X-ray emitters and that the distribution of the micro-wave to X-ray spectral slope $\alpha_{mu x}$ of LBL blazars is very narrow, confirming that the X-ray flux of most blazars is a very good estimator of their micro-wave emission. The X-ray spectral shape of all the objects that were observed long enough to allow spectral analysis is flat and consistent with inverse Compton emission within the commonly accepted view where the radiation from blazars is emitted in a Sychrotron-Inverse-Compton scenario. We predict that all blazars and most radio galaxies above the sensitivity limit of the WMAP and of the Planck CMB missions are X-ray sources detectable by the present generation of X-ray satellites. An hypothetical all-sky soft X-ray survey with sensitivity of approximately $10^{-15}$ erg/s would be crucial to locate and remove over 100,000 blazars from CMB temperature and polarization maps and therefore accurately clean the primordial CMB signal from the largest population of extragalactic foreground contaminants.

We analyze Swift gamma-ray bursts (GRBs) and X-ray afterglows for three GRBs with spectroscopic redshift determinations -- GRB 050401, XRF 050416a, and GRB 050525a. We find that the relation between spectral peak energy and isotropic energy of prompt emissions (the Amati relation) is consistent with that for the bursts observed in pre-Swift era. However, we find that the X-ray afterglow lightcurves, which extend up to 10-70 days, show no sign of the jet break that is expected in the standard framework of collimated outflows. We do so by showing that none of the X-ray afterglow lightcurves in our sample satisfies the relation between the spectral and temporal indices that is predicted for the phase after jet break. The jet break time can be predicted by inverting the tight empirical relation between the peak energy of the spectrum and the collimation-corrected energy of the prompt emission (the Ghirlanda relation). We find that there are no temporal breaks within the predicted time intervals in X-ray band. This requires either that the Ghirlanda relation has a larger scatter than previously thought, that the temporal break in X-rays is masked by some additional source of X-ray emission, or that it does not happen because of some unknown reason.

Gamma ray bursts (GRBs) are known to come in two duration classes, separated at ~2 s. Long bursts originate from star forming regions in galaxies, have accompanying supernovae (SNe) when near enough to observe and are likely caused by massive-star collapsars. Recent observations show that short bursts originate in regions within their host galaxies with lower star formation rates consistent with binary neutron star (NS) or NS - black hole (BH) mergers. Moreover, although their hosts are predominantly nearby galaxies, no SNe have been so far associated with short GRBs. We report here on the bright, nearby GRB 060614 that does not fit in either class. Its ~102 s duration groups it with long GRBs, while its temporal lag and peak luminosity fall entirely within the short GRB subclass. Moreover, very deep optical observations exclude an accompanying supernova, similar to short GRBs. This combination of a long duration event without accompanying SN poses a challenge to both a collapsar and merging NS interpretation and opens the door on a new GRB classification scheme that straddles both long and short bursts.

The consortium for Galactic studies with the Arecibo L-band Feed Array (ALFA) is conducting a neutral hydrogen (HI) survey of the whole Arecibo sky (declination range from -1 to 38 degrees), with high angular (3.5') and velocity resolution (0.2 km/s). The precursor observations with ALFA of a region in the Galactic anti-center reveal numerous isolated, small (a few pc in size), and cold (T_k<400 K) HI clouds at low negative velocities, distinctly separated from the HI disk emission (`low-velocity clouds', LVCs). These clouds are most likely located in the transition region between the Galactic disk and halo (at scale heights of 60--900 pc), yet they have properties of typical cold neutral clouds. LVCs are colder and, most likely, smaller and less massive than Lockman's clouds in the disk/halo interface region of the inner Galaxy. Our observations demonstrate that the cloudy structure of the interface region is most likely a general phenomenon, not restricted to the inner Galaxy. LVCs have sizes and radial velocities in agreement with the expectations for clouds formed in low-temperature fountain flows, although we measure a factor of ten higher HI column densities. Alternatively, LVCs could represent the final stages of the infalling intergalactic material in the on-going construction of the Galaxy. In the same dataset at higher negative velocities, we have discovered a `companion' HI cloud located 50' southwest of CHVC186+19-114. CHVC186+19-114 is a typical compact high velocity cloud (HVC) with a well-defined core/envelope structure. The companion cloud has a diameter of only 7'x9', and is one of the smallest HVCs known, most likely stripped from the main cloud through the interactions with the halo medium.

Thanks to the rapid position notice and response by HETE-2 and Swift, the X-ray afterglow emissions have been found for four recent short gamma-ray bursts (GRBs; GRB 050509b, GRB 050709, GRB 050724, and GRB 050813). The positions of three out of four short GRBs are coincident with galaxies with no current or recent star formation. This discovery tightens the case for a different origin for short and long GRBs. On the other hand, from the prompt emission point of view, a short GRB shows a harder spectrum comparing to that of the long duration GRBs according to the BATSE observations. We investigate the prompt emission properties of four short GRBs observed by Swift/BAT. We found that the hardness of all four BAT short GRBs is in between the BATSE range for short and long GRBs. We will discuss the spectral properties of short GRBs including the short GRB sample of Konus-Wind and HETE-2 to understand the hard nature of the BATSE short GRBs.

The long burst GRB 050717 was observed simultaneously by the Burst Alert Telescope (BAT) on Swift and the Konus instrument on Wind. Significant hard to soft spectral evolution was seen. Early gamma-ray and X-ray emission was detected by both BAT and the X-Ray Telescope (XRT) on Swift. The XRT continued to observe the burst for 7.1 days and detect it for 1.4 days. The X-ray light curve showed a classic decay pattern; the afterglow was too faint for a jet break to be detected. No optical, infrared or ultraviolet counterpart was discovered despite deep searches within 14 hours of the burst. Two particular features of the prompt emission make GRB 050717 a very unusual burst. First, the peak of the nu Fnu spectrum was observed to be 2401(-568/+781) keV for the main peak, which is the highest value of Epeak ever observed. Secondly, the spectral lag for GRB 050717 was determined to be 2.5 +- 2.6 ms, consistent with zero and unusually short for a long burst. This lag measurement suggests that this burst has a high intrinsic luminosity and hence is at high redshift (z > 2.7). Despite these unusual features GRB 050717 exhibits the classic prompt and afterglow behaviour of a gamma-ray burst.

We present observations of XRF 050406, the first burst detected by Swift showing a flare in its X-ray light curve. During this flare, which peaks at t_peak ~210s after the BAT trigger, a flux variation of (delta F)/F~6 in a very short time (delta t)/t_peak<<1 was observed. Its measured fluence in the 0.2-10 keV band was ~1.4x10^-8 erg cm^-2, which corresponds to 1-15% of the prompt fluence. We present indications of spectral variations during the flare. We argue that the producing mechanism is late internal shocks, which implies that the central engine is still active at 210s, though with a reduced power with respect to the prompt emission. The X-ray light curve flattens to a very shallow slope with decay index of ~0.5 after ~4400s, which also supports continued central engine activity at late times. This burst is classified as an X-ray flash, with a relatively low fluence (~10^-7 erg cm^-2 in the 15-350 keV band, E_iso~10^51 erg), a soft spectrum (photon index 2.65), no significant flux above ~50 keV and a peak energy E_p<15 keV. XRF 050406 is one of the first examples of a well-studied X-ray light curve of an XRF. We show that the main afterglow characteristics are qualitatively similar to those of normal GRBs. In particular, X-ray flares superimposed on a power-law light curve have now been seen in both XRFs and GRBs. This indicates that a similar mechanism may be at work for both kinds of events.

We report Swift Burst Alert Telescope (BAT) observations of the X-ray Flash (XRF) XRF 050416A. The fluence ratio between the 15-25 keV and 25-50 keV energy bands of this event is 1.5, thus making it the softest gamma-ray burst (GRB) observed by BAT so far. The spectrum is well fitted by the Band function with E^{\rm obs}_{\rm peak} of 15.0_{-2.7}^{+2.3} keV. Assuming the redshift of the host galaxy (z = 0.6535), the isotropic-equivalent radiated energy E_{\rm iso} and the peak energy at the GRB rest frame (E^{\rm src}_{\rm peak}) of XRF 050416A are not only consistent with the correlation found by Amati et al. and extended to XRFs by Sakamoto et al., but also fill-in the gap of this relation around the 30 - 80 keV range of E^{\rm src}_{\rm peak}. This result tightens the validity of the E^{\rm src}_{\rm peak} - E_{\rm iso} relation from XRFs to GRBs. We also find that the jet break time estimated using the empirical relation between E^{\rm src}_{\rm peak} and the collimation corrected energy E_{\gamma} is inconsistent with the afterglow observation by Swift X-ray Telescope. This could be due to the extra external shock emission overlaid around the jet break time or to the non existence of a jet break feature for XRF, which might be a further challenging for GRB jet emission, models and XRF/GRB unification scenarios.

The Burst Alert Telescope (BAT) is one of 3 instruments on the Swift MIDEX spacecraft to study gamma-ray bursts (GRBs). The BAT first detects the GRB and localizes the burst direction to an accuracy of 1-4 arcmin within 20 sec after the start of the event. The GRB trigger initiates an autonomous spacecraft slew to point the two narrow field-of-view (FOV) instruments at the burst location within 20-70 sec so to make follow-up x-ray and optical observations. The BAT is a wide-FOV, coded-aperture instrument with a CdZnTe detector plane. The detector plane is composed of 32,768 pieces of CdZnTe (4x4x2mm), and the coded-aperture mask is composed of approximately 52,000 pieces of lead (5x5x1mm) with a 1-m separation between mask and detector plane. The BAT operates over the 15-150 keV energy range with approximately 7 keV resolution, a sensitivity of approximately 10E-8 erg*cm^-2*s^-1, and a 1.4 sr (half-coded) FOV. We expect to detect >100 GRBs/yr for a 2-year mission. The BAT also performs an all-sky hard x-ray survey with a sensitivity of approximately 2 mCrab (systematic limit) and it serves as a hard x-ray transient monitor.

The prompt optical emission that arrives with gamma-rays from a cosmic gamma-ray burst (GRB) is a signature of the engine powering the burst, the properties of the ultra-relativistic ejecta of the explosion, and the ejecta's interactions with the surroundings. Until now, only GRB 990123 had been detected at optical wavelengths during the burst phase. Its prompt optical emission was variable and uncorrelated with the prompt gamma-ray emission, suggesting that the optical emission was generated by a reverse shock arising from the ejecta's collision with the surrounding material. Here we report prompt optical emission from GRB 041219a. It is variable and correlated with the prompt gamma-rays, indicating a common origin for the optical light and the gamma-rays. Within the context of the standard fireball model of GRBs, we attribute this new optical component to internal shocks driven into the burst ejecta by variations of the inner engine. The correlated optical emission is a direct probe of the jet isolated from the medium. The timing of the uncorrelated optical emission is strongly dependent on the nature of the medium.

The explosion that results in a cosmic gamma-ray burst (GRB) is thought to produce emission from two physical processes -- the activity of the central engine gives rise to the high-energy emission of the burst through internal shocking and the subsequent interaction of the flow with the external environment produces long-wavelength afterglow. While afterglow observations continue to refine our understanding of GRB progenitors and relativistic shocks, gamma-ray observations alone have not yielded a clear picture of the origin of the prompt emission nor details of the central engine. Only one concurrent visible-light transient has been found and was associated with emission from an external shock. Here we report the discovery of infrared (IR) emission contemporaneous with a GRB, beginning 7.2 minutes after the onset of GRB 041219a. Our robotic telescope acquired 21 images during the active phase of the burst, yielding the earliest multi-colour observations of any long-wavelength emission associated with a GRB. Analysis of an initial IR pulse suggests an origin consistent with internal shocks. This opens a new possibility to study the central engine of GRBs with ground-based observations at long wavelengths.

Magnetars comprise two classes of rotating neutron stars (Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars), whose X-ray emission is powered by an ultrastrong magnetic field, B ~ 10^15 G. Occasionally SGRs enter into active episodes producing many short X-ray bursts; extremely rarely (about once per 50 years per source), SGRs emit a giant flare, an event with total energy at least 1000 times higher than their typical bursts. Here we report that, on 2004 December 27, SGR 1806-20 emitted the brightest extra-solar transient event ever recorded, even surpassing the full moon brightness for 0.2 seconds. The total (isotropic) flare energy is 2x10^46 erg, 100 times higher than the only two previous events, making this flare a once in a century event. This colossal energy release likely occurred during a catastrophic reconfiguration of the magnetar's magnetic field. Such an event would have resembled a short, hard Gamma Ray Burst (GRB) if it had occurred within 40 Mpc, suggesting that extragalactic SGR flares may indeed form a subclass of GRBs.

Swift has two main types of ``rate'' triggers: short and long. Short trigger time scales range from 4ms to 64ms, while long triggers are 64ms to ~ 16 seconds. While both short and long trigger have criteria with one background sample (traditional ``one-sided'' triggers), the long triggers can also have criteria with two background samples (``bracketed'' triggers) which remove trends in the background. Both long and short triggers can select energy ranges of 15-25, 15-50, 25-100 and 50-350 KeV. There are more than 180 short triggering criteria and approximately 500 long triggering criteria used to detect gamma ray bursts. To fully utilize these criteria, the thresholds must be set correctly. The optimum thresholds are determined by a tradeoff between avoiding false triggers and capturing as many bursts as possible. We use realistic simulated orbital variations, which are the prime cause of false triggers.

The BAT instrument tells the Swift satellite where to point to make immediate follow-up observations of GRBs. The science software on board must efficiently process gamma-ray events coming in at up to 34 kHz, identify rate increases that could be due to GRBs while disregarding those from known sources, and produce images to accurately and rapidly locate new Gamma-ray sources.

The Swift Burst Alert Telescope (BAT) is a huge (5200 cm2) coded aperture imager that will detect gamma-ray bursts in real time and provide a location that the Swift satellite will use to slew the optical and x-ray telescopes. The huge size of BAT is a challenge for the on-board triggering: a change as small as 1% is equivalent to a 1 sigma statistical variation in 1 second. There will be three types of triggers, two based on rates and one based on images. The first type of trigger is for short time scales (4 msec to 64 msec). These will be traditional triggers (single background) and we check about 25,000 combinations of time-energy-focal plane subregions per second. The second type of trigger will be similar to what is used on HETE: fits to multiple background regions to remove trends for time scales between 64 msec and 64 seconds. About 500 triggers will be checked per second. For these rate triggers, false triggers and variable non-GRB sources will be rejected by requiring a new source to be present in an image. The third type of trigger works on longer time scales (minutes), and will be based on routine images that are made of the field of view.

The Swift satellite will be a self-contained observatory that will bring new capabilities to the observing of the early afterglow emission of Gamma-ray Bursts. Swift is completely autonomous and will do all of the observations without help from the ground. There are three instruments on Swift. A large (5200 sq cm) coded aperture imager will locate the bursts within about 15 seconds. The satellite will be able to slew to point at the location within a minute or two. There are two narrow field of view instruments: an optical telescope and an x-ray telescope. Thus, Swift will provide simultaneous gamma-ray, x-ray, and optical observations of Gamma-ray bursts soon after the burst. A key to the success of Swift will be its ability to detect and locate a large number of gamma-ray bursts quick enough that the narrow field of view instruments can follow up. The results of simulations show that Swift will be able to detect about 300 bursts a year and locate about 150. The number that Swift will be able to slew to depends on constraints built into the satellite bus. Preliminary results indicate that we might be able to slew to 100 bursts per year, but that is heavily dependent on satellite operations.

The Gamma Ray Imaging Spectrometer (GRIS) recently observed the gamma-ray emission from the Galactic center region. We have detected the 1809 keV Galactic 26Al emission at a significance level of 6.8-sigma but have found no evidence for emission at 1173 keV and 1332 keV, expected from the decay chain of the nucleosynthetic 60Fe. The isotopic abundances and fluxes are derived for different source distribution models. The resulting abundances are between 2.6+-0.4 and 4.5+-0.7 Solar Masses for 26Al and a 2-sigma upper limit for 60Fe between 1.7 and 3.1 Solar Masses. The measured 26Al emission flux is significantly higher than that derived from the CGRO/COMPTEL 1.8 MeV sky map. This suggests that a fraction of the 26Al emission may come from extended sources with a low surface brightness that are invisible to COMPTEL. We obtain a 60Fe to 26Al flux ratio 2-sigma upper limit of 0.14, which is slightly lower than the 0.16 predicted from current nucleosynthesis models assuming that SNII are the major contributors to the galactic 26Al. Since the uncertainties in the predicted fluxes are large (up to a factor of 2), our measurement is still compatible with the theoretical expectations.