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 Epoch of Reionization (EoR)---when neutral hydrogen in the intergalactic medium was systematically ionized---is a period in our Universe's history that is currently poorly understood. However, a key prediction of most models is a correlation between the density and ionization field during the EoR. This has consequences for the 21cm power spectrum. Here, we propose a parametrization for the density-ionization correlation and study the dependence of the 21cm power spectrum on this parameterization. We use this formalism to forecast the ability of current and future observations to constrain these correlations. We find that upper limits on the dimensionless power spectrum at redshifts $7.5 < z < 8.5$ using $k$ bins between $0.1\,\textrm{Mpc}^{-1} < k < 0.75\,\textrm{Mpc}^{-1}$ with error bars at the level of $\sim\!\!20\,\textrm{mK}^2$ about our fiducial model would rule out uncorrelated reionization at $99\%$ credibility. Beyond upper limits, we find that at its full sensitivity, the Hydrogen Epoch of Reionization Array (HERA) will be able to place strong constraints on the sign and magnitude of density-ionization correlations.

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.

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

We study the chemical evolution of galaxy clusters by measuring the iron mass in the ICM after dissecting the abundance profiles into different components. We use Chandra archival observations of 186 morphologically regular clusters in the redshift range [0.04, 1.07]. For each cluster we compute the iron abundance and gas density profiles. We aim at identifying in the iron distribution a central peak associated with the BCG, and an approximately constant plateau associated with early enrichment. We are able to firmly identify the two components in a significant fraction of the sample, simply relying on the fit of the abundance profile. We compute the iron mass included in the iron peak and plateau, and the gas mass-weighted iron abundance out to $r_{500}$. While the iron plateau shows no evolution, we find marginal decrease with redshift in the iron peak. We measure that the fraction of iron peak mass is typically a few percent (~1%) of the total iron mass within $r_{500}$. Therefore, since the total iron mass budget is dominated by the plateau, we find consistently that the global gas mass-weighted iron abundance does not evolve significantly. We are also able to reproduce past claims of evolution in the global iron abundance, which turn out to be due to the use of cluster samples with different selection methods combined to the use of emission-weighted instead of gas mass-weighted abundance values. Finally, while the intrinsic scatter in the iron plateau mass is consistent with zero, the iron peak mass exhibits a large scatter, in line with the fact that the peak is produced after the virialization of the halo and depends on the formation of the hosting cool core and the associated feedback processes. We conclude that only a spatially-resolved approach can resolve the issue of the ICM iron evolution, reconciling the contradictory results obtained in the last ten years.

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.

The Packed Ultra-wideband Mapping Array (PUMA) is a proposed low-resolution transit interferometric radio telescope operating over the frequency range 200 - 1100MHz. Its rich science portfolio will include measuring structure in the universe from redshift z = 0.3 to 6 using 21cm intensity mapping, detecting one million fast radio bursts, and monitoring thousands of pulsars. It will allow PUMA to advance science in three different areas of physics (the physics of dark energy, the physics of cosmic inflation and time-domain astrophysics). This document is a response to a request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Decadal Survey on Astronomy and Astrophysics 2020. We present the science case of PUMA, the development path and major risks to the project.

We use semi-analytic galaxy catalogs based on two high-resolution cosmological $N$-body simulations, Millennium-WMAP7 and Millennium-II, to investigate the formation of the Local Group (LG) analogs. Unlike previous studies, we use the observed stellar masses to select the LG member (Milky Way (MW) and M31) analogs, and then impose constrains using the observed separation, isolation, and kinematics of the two main member galaxies. By comparing radial and low-ellipticity orbits between the MW and M31, we find higher tangential velocity results in higher total mass, which are 4.4$^{+2.4}_{-1.5}\times$10$^{12}\rm M_{\odot}$ and 6.6$^{+2.7}_{-1.5}\times$10$^{12}\rm M_{\odot}$ for radial and low-ellipticity orbits. The orbits also influence the individual mass distribution of MW and M31 analogs. For radial orbits, the typical host halo masses of the MW and M31 are 1.5$^{+1.4}_{-0.7}\times$10$^{12}\rm M_{\odot}$ and 2.5$^{+1.3}_{-1.1}\times$10$^{12}\rm M_{\odot}$; for low-ellipticity orbits, the masses are 2.5$^{+2.2}_{-1.4}\times$10$^{12}\rm M_{\odot}$ and 3.8$^{+2.8}_{-1.8}\times$10$^{12} \rm M_{\odot}$. The LG is located primarily in filaments with tails extending toward higher densities up to $\delta\sim4.5$. The dark matter velocity anisotropy parameters $\beta$ of both the MW and M31 analogs are close to zero in the center, increasing to 0.2--0.3 at 50--80 kpc and decreasing slowly outward. The slope is much flatter than computed from the MW satellites, and the amplitude is smaller than traced by halo stars. Values of $\beta$ from different tracers agree at $\sim$120 kpc where $\beta \sim$ 0.2. We also find that model predictions agree broadly with observations in the radial distribution and luminosity function of satellites around the MW and M31.

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.

The redshifted 21 cm line is an emerging tool in cosmology, in principle permitting three-dimensional surveys of our Universe that reach unprecedentedly large volumes, previously inaccessible length scales, and hitherto unexplored epochs of our cosmic timeline. Large radio telescopes have been constructed for this purpose, and in recent years there has been considerable progress in transforming 21 cm cosmology from a field of considerable theoretical promise to one of observational reality. Increasingly, practitioners in the field are coming to the realization that the success of observational 21cm cosmology will hinge on software algorithms and analysis pipelines just as much as it does on careful hardware design and telescope construction. This review provides a pedagogical introduction to state-of-the-art ideas in 21 cm data analysis, covering a wide variety of steps in a typical analysis pipeline, from calibration to foreground subtraction to mapmaking to power spectrum estimation to parameter estimation.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

The Canadian Hydrogen Observatory and Radio-transient Detector (CHORD) is a next-generation radio telescope, proposed for construction to start immediately. CHORD is a pan-Canadian project, designed to work with and build on the success of the Canadian Hydrogen Intensity Mapping Experiment (CHIME). It is an ultra-wideband, "large-N, small-D" telescope, consisting of a central array of 512x6-m dishes, supported by a pair of distant outrigger stations, each equipped with CHIME-like cylinders and a 64-dish array. CHORD will measure the distribution of matter over a huge swath of the Universe, detect and localize tens of thousands of Fast RadioBursts (FRBs), and undertake cutting-edge measurements of fundamental physics.

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.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

The application of machine learning (ML) methods to the analysis of astrophysical datasets is on the rise, particularly as the computing power and complex algorithms become more powerful and accessible. As the field of ML enjoys a continuous stream of breakthroughs, its applications demonstrate the great potential of ML, ranging from achieving tens of millions of times increase in analysis speed (e.g., modeling of gravitational lenses or analysing spectroscopic surveys) to solutions of previously unsolved problems (e.g., foreground subtraction or efficient telescope operations). The number of astronomical publications that include ML has been steadily increasing since 2010. With the advent of extremely large datasets from a new generation of surveys in the 2020s, ML methods will become an indispensable tool in astrophysics. Canada is an unambiguous world leader in the development of the field of machine learning, attracting large investments and skilled researchers to its prestigious AI Research Institutions. This provides a unique opportunity for Canada to also be a world leader in the application of machine learning in the field of astrophysics, and foster the training of a new generation of highly skilled researchers.

Line-intensity mapping of the 21cm line is a powerful probe of large scale structure at z<6, tracing large-scale structure via neutral hydrogen content that is found within galaxies. In principle, it enables cost-efficient surveys of the matter distribution up to z~6, unlocking orders of magnitude more modes for observational cosmology. Canada has been a traditional leader in this field, having led the first detections of the cosmological 21cm signal via cross-correlations with optical galaxy surveys and having constructed the Canadian Hydrogen Intensity Mapping Experiment (CHIME). The field is now entering a new era where data is abundant, allowing studies in how to overcome systematics to be tackled in an empirical, head-on fashion. In the next few years, this will produce the first detection of the 21cm auto power spectrum, which will pave the way towards a large suite of scientific possibilities. These potentially include precision measurements on the dark energy equation of state and other LCDM parameters, constraints on how HI mass traces dark matter, a detection of neutrino effects on large-scale structure, and the use of 21cm lensing to further constrain cosmology. To turn these promising directions into reality, we recommend a sustained program of investment in 21cm cosmology, starting with funding for the Canadian Hydrogen Observatory and Radio transient Detector (CHORD), followed by small-scale development efforts targeting next-generation hardware and sustained support for theory and technical staff support. Additionally, Canada should invest in complementary line-intensity mapping efforts (such as with CO or [CII] lines) and maintain participation in next-generation international efforts such as the Packed Ultra-wideband Mapping Array (PUMA) and the Square Kilometre Array (SKA).

In the next few years, the 21cm line will enable direct observations of the Dark Ages, Cosmic Dawn, and Reionization, which represent previously unexplored periods in our cosmic history. With a combination of sky-averaged global signal measurements and spatial mapping surveys, the possible science return is enormous. This potentially includes (but is not limited to) constraints on first-generation galaxies (such as their typical masses and luminosities), constraints on cosmological parameters, a measurement of the Hubble parameter at z~15 to 20, the elimination of the optical depth nuisance parameter in the CMB, and searches for exotic phenomena such as baryon-dark matter couplings. To enable continued Canadian leadership in these science opportunities, we recommend 1) continued investments in 21cm experiments at all redshifts, 2) detailed analysis efforts of current data to overcome systematic effects, 3) new investments in preliminary experiments to explore the truly low-frequency sky as a stepping stone towards the Dark Ages, 4) new investments in line-intensity mapping experiments beyond the 21cm line, 5) continued theory support for 21cm cosmology, and 6) continued participation and knowledge transfer to next-generation international efforts such as the Square Kilometre Array.

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.

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

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.

Based on HST and MUSE data, we probe the stellar and gas properties (i.e. kinematics, stellar mass, star formation rate) of the radio-loud brightest cluster galaxy (BCG) located at the centre of the X-ray luminous cool core cluster Abell 2667 (z = 0.2343). The bi-dimensional modelling of the BCG surface brightness profile reveals the presence of a complex system of substructures extending all around the galaxy. Clumps of different size and shape plunged into a more diffuse component constitute these substructures, whose intense 'blue' optical colour hints to the presence of a young stellar population. Our results depict the BCG as a massive (M_star ~ 1.38 x 10^11 M_sun) dispersion-supported spheroid (v_star < 150 km/s, sigma_0 ~ 216 km/s) hosting an active supermassive black hole (M_SMBH ~ 3.8 x 10^9 M_sun) whose optical features are typical of low ionisation nuclear emission line regions. Although the velocity pattern of the stars in the BCG is irregular, the stellar kinematics in the regions of the clumps show a positive velocity of ~ 100 km/s, similarly to the gas component. An analysis of the mechanism giving rise to the observed lines in the clumps through empirical diagnostic diagrams points out that the emission is composite, suggesting the contribution from both star formation and AGN. We conclude our analysis describing how scenarios of both chaotic cold accretion and merging with a gas-rich disc galaxy can efficaciously explain the phenomena the BCG is undergoing.

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.

A central drop of ICM Fe abundance has been observed in several cool-core clusters. It has been proposed that this abundance drop may be due, at least partially, to the depletion of Fe into dust grains in the central, high-density regions. According to this scenario, noble gas elements such as Ar and Ne are not expected to be depleted into dust, and therefore should not show any drop, but follow the general increase of metal abundance toward the center. In this work, we test this scenario by measuring with {\sl Chandra} data the radial profiles of Ar and Ne in a sample of 12 groups and clusters where a central drop in Fe abundance has been detected. We confirm the presence of the Fe drop in 10 out of 12 clusters at more than 2$\sigma$ c.l., and 4 Ar drops with similar significance. We also find 4 Ne drops, with the caveat that Ne abundance measurement from CCD spectra suffers from systematics not completely under control. Our results are consistent with an abundance drop common to the three elements. When comparing the profiles, we find that, on average, the abundance profiles of Ar and Ne are significantly higher than Fe and steeper toward the center, while they all gradually become consistent with solar composition at $r\geq 0.05r_{500}$. We also check that Si and S profiles are mostly consistent with Fe. This result confirms a scenario in which some fraction of Fe is depleted into dust grains in the inner regions, although the global central abundance drop is mostly due to mechanical processes, like the displacement of metal-rich ICM from the very center to larger radii by AGN-driven feedback. Finally, we report the detection of an Fe drop in the cluster MACSJ1423.8+2404 at $z=0.543$, showing that this feature appears early on in cool-core clusters.

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 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 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 "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.

We search for global rotation of the intracluster medium (ICM) in the galaxy cluster Abell 2107, where previous studies have detected rotational motion in the member galaxies with a high significance level. By fitting the centroid of the iron $K_{\alpha}$ line complex at 6.7--6.9 keV rest frame in {\sl Chandra} ACIS-I spectra, we identify the possible rotation axis with the line that maximizes the difference between the emission-weighted spectroscopic redshift measured in the two halves defined by the line itself. Then, we measure the emission-weighted redshift in linear regions parallel to the preferred rotation axis, and find a significant gradient as a function of the projected distance from the rotation axis, compatible with a rotation pattern with maximum tangential velocity ${\tt v}_{\rm max}=1380\pm 600$ km/s at a radius $\lambda_0\sim 160$ kpc. This result, if interpreted in the framework of hydrostatic equilibrium, as suggested by the regular morphology of Abell 2107, would imply a large mass correction of the order of $\Delta M = (6 \pm 4)\times 10^{13} M_\odot$ at $\sim 160$ kpc, which is incompatible with the cluster morphology itself. A more conservative interpretation may be provided by an unnoticed off-center, head-on collision between two comparable halos. Our analysis confirms the peculiar dynamical nature of the otherwise regular cluster Abell 2107, but is not able to resolve the rotation vs merger scenario, a science case that can be addressed by the next-generation X-ray facilities carrying X-ray bolometers onboard.

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 21-cm power spectrum (PS) has been shown to be a powerful discriminant of reionization and cosmic dawn astrophysical parameters. However, the 21-cm tomographic signal is highly non-Gaussian. Therefore there is additional information which is wasted if only the PS is used for parameter recovery. Here we showcase astrophysical parameter recovery directly from 21-cm images, using deep learning with convolutional neural networks (CNN). Using a database of 2D images taken from 10,000 21-cm lightcones (each generated from different cosmological initial conditions), we show that a CNN is able to recover parameters describing the first galaxies: (i) Tvir , their minimum host halo virial temperatures (or masses) capable of hosting efficient star formation; (ii) {\zeta} , their typical ionizing efficiencies; (iii) LX/SFR , their typical soft-band X-ray luminosity to star formation rate; and (iv) E0 , the minimum X-ray energy capable of escaping the galaxy into the IGM. For most of their allowed ranges, log Tvir and log LX/SFR are recovered with < 1% uncertainty, while {\zeta} and E0 are recovered within 10% uncertainty. Our results are roughly comparable to the accuracy obtained from Monte Carlo Markov Chain sampling of the PS with 21CMMC for the two mock observations analyzed previously, although we caution that we do not yet include noise and foreground contaminants in this proof-of-concept study.

44325 measured and assigned transitions of H$_2^{32}$S, the parent isotopologue of the hydrogen sulfide molecule, are collated from 33 publications into a single database and reviewed critically. Based on this information, rotation-vibration energy levels are determined for the ground electronic state using the Measured Active Rotational-Vibrational Energy Levels (MARVEL) technique. The ortho and para principal components of the measured spectroscopic network of H$_2^{32}$S are considered separately. The verified set of 25293 ortho- and 18778 para- H$_2^{32}$S transitions determine 3969 ortho and 3467 para energy levels. The Marvel results are compared with alternative data compilations, including a theoretical variational linelist.

Active Galactic Nuclei (AGN) and star-forming galaxies are leading candidates for being the luminous sources that reionized our Universe. Next-generation 21cm surveys are promising to break degeneracies between a broad range of reionization models, hence revealing the nature of the source population. While many current efforts are focused on a measurement of the 21cm power spectrum, some surveys will also image the 21cm field during reionization. This provides further information with which to determine the nature of reionizing sources. We create a Convolutional Neural Network (CNN) that is efficiently able to distinguish between 21cm maps that are produced by AGN versus galaxies scenarios with an accuracy of 92-100%, depending on redshift and neutral fraction range. An exception to this is when our Universe is highly ionized, since the source models give near-identical 21cm maps in that case. When adding thermal noise from typical 21cm experiments, the classification accuracy depends strongly on the effectiveness of foreground removal. Our results show that if foregrounds can be removed reasonably well, SKA, HERA and LOFAR should be able to discriminate between source models with greater accuracy at a fixed redshift. Only future SKA 21cm surveys are promising to break the degeneracies in the power spectral analysis.

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.

We investigate the spatial distribution of iron in the intra-cluster medium in a selected sample of 41 relaxed clusters in the redshift range $0.05 < z < 1.03$ using {\sl Chandra} archival data. We compute the azimuthally-averaged, deprojected $Z_{\rm Fe}$ profile of each cluster out to $\sim 0.4r_{500}$, and identify a peak in the distribution of iron followed by a flatter distribution at larger radii. Due to the steep gradient both in gas density and abundance, we find that the emission-weighted iron abundance within $0.2r_{500}$, which entirely includes the iron peak in most of the cases, is on average $\sim$25\% higher than the mass-weighted value, showing that spatially resolved analysis and accurate deprojection are key to study the evolution of iron enrichment in the very central regions of cool core clusters. We quantify the extent of the iron distribution in each cluster with a normalized scale parameter $r_{\rm Fe}$, defined as the radius where the iron abundance excess is half of its peak value. We find that $r_{\rm Fe}$ increases by a factor of $\sim 3$ from $z\sim 1$ to $z\sim 0.1$, suggesting that the spatial distribution of iron in the ICM extends with time, possibly due to the mixing with the mechanical-mode feedback from the central galaxy. We also find that the iron mass excess within $0.3r_{500}$, when normalized to the total baryonic mass within the same region, does not evolve significantly, showing that this iron mass component is already established at $z\sim 1$.

We study anisotropic inflationary solutions in massive Gauge-flation. We work with the theory in both the Stueckelberg and dynamical symmetry-breaking limits and demonstrate that extended periods of accelerated anisotropic expansion are possible. In the case of dynamical symmetry breaking, we show that spacetime can transition from isotropic quasi-de Sitter space to an accelerating Bianchi spacetime due to a rolling Higgs field - the spacetime can develop hair. Similarly, symmetry restoring transitions are possible from accelerating Bianchi spacetime to quasi-de Sitter space - the spacetime can lose its hair. These transitions can be arranged to occur quickly, within an $e$-folding or so, or over tens of $e$-folds.

To investigate the dynamics of the galaxy cluster A2142, we compile an extended catalog of 2239 spectroscopic redshifts of sources, including newly measured 237 redshifts, within 30 arcmin from the cluster center. With the $\sigma$-plateau algorithm from the caustic method, we identify 868 members and a number of substructures in the galaxy distribution both in the outskirts, out to $\sim$3.5 Mpc from the cluster center, and in the central region. In the outskirts, one substructure overlaps a falling clump of gas previously identified in the X-ray band. These substructures suggests the presence of multiple minor mergers, which are responsible for the complex dynamics of A2142, and the absence of recent or ongoing major mergers. We show that the distribution of the galaxies in the cluster core and in several substructures are consistent with the mass distribution inferred from the weak lensing signal. Moreover, we use spatially-resolved X-ray spectroscopy to measure the redshift of different regions of the intracluster medium within $\sim$3 arcmin from the cluster center. We find a ring of gas near the two X-ray cold fronts identified in previous analyses and measure a velocity of this ring of $810 \pm 330 {\rm km s^{-1}}$ larger than the cluster mean velocity. Our analysis suggests the presence of another ring surrounding the core, whose velocity is $660 \pm 300 {\rm km s^{-1}}$ larger than the cluster velocity. These X-ray features are not associated to any optical substructures, and support the core-sloshing scenario suggested in previous work.

Next-generation 21cm observations will enable imaging of reionization on very large scales. These images will contain more astrophysical and cosmological information than the power spectrum, and hence providing an alternative way to constrain the contribution of different reionizing sources populations to cosmic reionization. Using Convolutional Neural Networks, we present a simple network architecture that is sufficient to discriminate between Galaxy-dominated versus AGN-dominated models, even in the presence of simulated noise from different experiments such as the HERA and SKA.

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.

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.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

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

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between 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.

We present an X-ray photometric analysis of six gravitationally lensed quasars, with observation campaigns spanning from 5 to 14 years, measuring the total (0.83 - 21.8 keV restframe), soft (0.83 - 3.6 keV), and hard (3.6 - 21.8 keV) band image flux ratios or each epoch. Using the ratios of the model-predicted macro-magnifications as baselines, we build differential microlensing light curves and obtain joint likelihood functions for the average X-ray emission region sizes. Our analysis yields a Probability Distribution Function for the average half-light radius of the X-Ray emission region in the sample that peaks slightly above 1 gravitational radius and with nearly indistinguishable 68% confidence (one-sided) upper limits of 17.8 and 18.9 gravitational radii for the soft and hard X-ray emitting regions, assuming a mean stellar mass of 0.3 solar masses. We see hints of energy dependent microlensing between the soft and hard bands in two of the objects. In a separate analysis on the root-mean-square of the microlensing variability, we find significant differences between the soft and hard bands but the sign of the difference is not consistent across the sample. This suggests the existence of some kind of spatial structure to the X-ray emission in an otherwise extremely compact source. We also discover a correlation between the RMS microlensing variability and the average microlensing amplitude.

The recent discovery of gravitational waves from merging black holes has generated interest in primordial black holes as a possible component of the dark matter. In this paper, we show that pulsar timing may soon have sufficient data to constrain $1$-$1000\,M_{\odot}$ primordial black holes via the non-detection of a third-order Shapiro time delay as the black holes move around the Galactic halo. We present the results of a Monte Carlo simulation which suggests that future data from known pulsars may be capable of constraining the PBH density more stringently than other existing methods in the mass range ~1-30$\,M_{\odot}$. We find that timing new pulsars discovered using the proposed Square Kilometre Array may constrain primordial black holes in this mass range to comprise less than ~1-10% of the dark matter.

We present a new method for interferometric imaging that is ideal for the large fields of view and compact arrays common in 21 cm cosmology. We first demonstrate the method with simulations for two very different low frequency interferometers, the Murchison Widefield Array (MWA) and the MIT Epoch of Reionization (MITEoR) Experiment. We then apply the method to the MITEoR data set collected in July 2013 to obtain the first northern sky map from 128 MHz to 175 MHz at about 2 degree resolution, and find an overall spectral index of -2.73+/-0.11. The success of this imaging method bodes well for upcoming compact redundant low-frequency arrays such as HERA. Both the MITEoR interferometric data and the 150 MHz sky map are publicly available at http://space.mit.edu/home/tegmark/omniscope.html.

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.

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

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.

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

We search for bulk motions in the intracluster medium (ICM) of massive clusters showing evidence of an ongoing or recent major merger with spatially resolved spectroscopy in {\sl Chandra} CCD data. We identify a sample of 6 merging clusters with $>$150 ks {\sl Chandra} exposure in the redshift range $0.1 < z < 0.3$. By performing X-ray spectral analysis of projected ICM regions selected according to their surface brightness, we obtain the projected redshift maps for all of these clusters. After performing a robust analysis of the statistical and systematic uncertainties in the measured X-ray redshift $z_{\rm X}$, we check whether or not the global $z_{\rm X}$ distribution differs from that expected when the ICM is at rest. We find evidence of significant bulk motions at more than 3$\sigma$ in A2142 and A115, and less than 2$\sigma$ in A2034 and A520. Focusing on single regions, we identify significant localized velocity differences in all of the merging clusters. We also perform the same analysis on two relaxed clusters with no signatures of recent mergers, finding no signs of bulk motions, as expected. Our results indicate that deep {\sl Chandra} CCD data enable us to identify the presence of bulk motions at the level of $v_{\rm BM} >$ 1000\ ${\rm km\ s^{-1}}$ in the ICM of massive merging clusters at $0.1<z<0.3$. Although the CCD spectral resolution is not sufficient for a detailed analysis of the ICM dynamics, {\sl Chandra} CCD data constitute a key diagnostic tool complementing X-ray bolometers on board future X-ray missions.

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.

We propose a strategy to search for bulk motions in the intracluster medium (ICM) of merging clusters based on {\sl Chandra} CCD data. Our goal is to derive robust measurements of the average redshift of projected ICM regions obtained from the centroid of the $K_\alpha$ line emission. We thoroughly explore the effect of the unknown temperature structure along the line of sight to accurately evaluate the systematic uncertainties on the ICM redshift. We apply our method to the "Bullet cluster" (1E~0657-56). We directly identify 23 independent regions on the basis of the surface brightness contours, and measure the redshift of the ICM averaged along the line of sight in each. We find that the redshift distribution across these regions is marginally inconsistent with the null hypothesis of a constant redshift or no bulk motion in the ICM, at a confidence level of about $2\, \sigma$. We tentatively identify the regions most likely affected by bulk motions and find a maximum velocity gradient of about $(46\pm 13)$ $\rm km~s^{-1}~kpc^{-1}$ along the line of sight on a scale of $\sim 260 $ kpc along the path of the "bullet." We interpret this as the possible signature of a significant mass of ICM pushed away along a direction perpendicular to the merging. This preliminary result is promising for a systematic search for bulk motions in bright, moderate-redshift clusters based on spatially resolved spectral analysis of {\sl Chandra} CCD data. This preliminary result is promising for a systematic search for bulk motions in bright, moderate-redshift clusters based on spatially resolved spectral analysis of {\sl Chandra} CCD data.

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.

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.

Spectral measurements of the 21 cm monopole background have the promise of revealing the bulk energetic properties and ionization state of our universe from z ~ 6-30. Synchrotron foregrounds are orders of magnitude larger than the cosmological signal, and are the principal challenge faced by these experiments. While synchrotron radiation is thought to be spectrally smooth and described by relatively few degrees of freedom, the instrumental response to bright foregrounds may be much more complex. To deal with such complexities, we develop an approach that discovers contaminated spectral modes using spatial fluctuations of the measured data. This approach exploits the fact that foregrounds vary across the sky while the signal does not. The discovered modes are projected out of each line-of-sight of a data cube. An angular weighting then optimizes the cosmological signal amplitude estimate by giving preference to lower-noise regions. Using this method, we show that it is essential for the passband to be stable to at least ~10^{-4}. In contrast, the constraints on the spectral smoothness of the absolute calibration are mainly aesthetic if one is able to take advantage of spatial information. To the extent it is understood, controlling polarization to intensity leakage at the ~10^{-2} level will also be essential to rejecting Faraday rotation of the polarized synchrotron emission.

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.

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

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

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

[Abridged] The spatially averaged global spectrum of the redshifted 21cm line has generated much experimental interest, for it is potentially a direct probe of the Epoch of Reionization and the Dark Ages. Since the cosmological signal here has a purely spectral signature, most proposed experiments have little angular sensitivity. This is worrisome because with only spectra, the global 21cm signal can be difficult to distinguish from foregrounds such as Galactic synchrotron radiation, as both are spectrally smooth and the latter is orders of magnitude brighter. We establish a mathematical framework for global signal data analysis in a way that removes foregrounds optimally, complementing spectra with angular information. We explore various experimental design trade-offs, and find that 1) with spectral-only methods, it is impossible to mitigate errors that arise from uncertainties in foreground modeling; 2) foreground contamination can be significantly reduced for experiments with fine angular resolution; 3) most of the statistical significance in a positive detection during the Dark Ages comes from a characteristic high-redshift trough in the 21cm brightness temperature; and 4) Measurement errors decrease more rapidly with integration time for instruments with fine angular resolution. We show that if observations and algorithms are optimized based on these findings, an instrument with a 5 degree beam can achieve highly significant detections (greater than 5-sigma) of even extended (high Delta-z) reionization scenarios after integrating for 500 hrs. This is in contrast to instruments without angular resolution, which cannot detect gradual reionization. Abrupt ionization histories can be detected at the level of 10-100's of sigma. The expected errors are also low during the Dark Ages, with a 25-sigma detection of the expected cosmological signal after only 100 hrs of integration.

We develop and demonstrate an acceleration of the Liu & Tegmark quadratic estimator formalism for inverse variance foreground subtraction and power spectrum estimation in 21 cm tomography from O(N^3) to O(N log N), where N is the number of voxels of data. This technique makes feasible the megavoxel scale analysis necessary for current and upcoming radio interferometers by making only moderately restrictive assumptions about foreground models and survey geometry. We exploit iterative and Monte Carlo techniques and the symmetries of the foreground covariance matrices to quickly estimate the 21 cm brightness temperature power spectrum, P(k_parallel, k_perpendicular), the Fisher information matrix, the error bars, the window functions, and the bias. We also extend the Liu & Tegmark foreground model to include bright point sources with known positions in a way that scales as O[(N log N)(N point sources)] < O(N^5/3). As a first application of our method, we forecast error bars and window functions for the upcoming 128-tile deployment of the Murchinson Widefield Array, showing that 1000 hours of observation should prove sufficiently sensitive to detect the power spectrum signal from the Epoch of Reionization.

Before it becomes a sensitive probe of the Epoch of Reionization, the Dark Ages, and fundamental physics, 21 cm tomography must successfully contend with the issue of foreground contamination. Broadband foreground sources are expected to be roughly four orders of magnitude larger than any cosmological signals, so precise foreground models will be necessary. Such foreground models often contain a large number of parameters, reflecting the complicated physics that governs foreground sources. In this paper, we concentrate on spectral modeling (neglecting, for instance, bright point source removal from spatial maps) and show that 21 cm tomography experiments will likely not be able to measure these parameters without large degeneracies, simply because the foreground spectra are so featureless and generic. However, we show that this is also an advantage, because it means that the foregrounds can be characterized to a high degree of accuracy once a small number of parameters (likely three or four, depending on one's instrumental specifications) are measured. This provides a simple understanding for why 21 cm foreground subtraction schemes are able to remove most of the contaminants by suppressing just a small handful of simple spectral forms. In addition, this suggests that the foreground modeling process should be relatively simple and will likely not be an impediment to the foreground subtraction schemes that are necessary for a successful 21 cm tomography experiment.

We report the discovery of rapid periodic signals in the light curves of two cataclysmic variables with prominent white-dwarf components in their spectra, SDSS J1457+51 and BW Sculptoris. These stars therefore appear to be new members of the GW Lib class of variable star, in which the fast periodic (and non-commensurate with the orbital period) signals are believed to arise from non-radial pulsations in the underlying white dwarf. The power spectra of both stars show complex signals with primary periods near 10 and 20 minutes. These signals change in frequency by a few percent on a timescale of weeks or less, and probably contain an internal fine structure unresolved by our observations. We also detect double-humped waves signifying the underlying orbital periods, near 78 minutes for both stars. In addition, BW Scl shows a transient but powerful signal with a period near 87 minutes, a quiescent superhump. The 11 % excess over the orbital period is difficult to understand, and may arise from an eccentric instability near the 2:1 resonance in the accretion disc.

21cm tomography promises to be a powerful tool for estimating cosmological parameters, constraining the epoch of reionization, and probing the so-called dark ages. However, realizing this promise will require the extraction of a cosmological power spectrum from beneath overwhelmingly large sources of foreground contamination. In this paper, we develop a unified matrix-based framework for foreground subtraction and power spectrum estimation, which allows us to quantify the errors and biases that arise in the power spectrum as a result of foreground subtraction. We find that existing line-of-sight foreground subtraction proposals can lead to substantial mode-mixing as well as residual noise and foreground biases, whereas our proposed inverse variance foreground subtraction eliminates noise and foreground biases, gives smaller error bars, and produces less correlated measurements of the power spectrum. We also numerically confirm the intuitive belief in the literature that 21cm foreground subtraction is best done using frequency rather than angular information.

Growing interest in 21 cm tomography has led to the design and construction of broadband radio interferometers with low noise, moderate angular resolution, high spectral resolution, and wide fields of view. With characteristics somewhat different from traditional radio instruments, these interferometers may require new calibration techniques in order to reach their design sensitivities. Self-calibration or redundant calibration techniques that allow an instrument to be calibrated off complicated sky emission structures are ideal. In particular, the large number of redundant baselines possessed by these new instruments makes redundant calibration an especially attractive option. In this paper, we explore the errors and biases in existing redundant calibration schemes through simulations, and show how statistical biases can be eliminated. We also develop a general calibration formalism that includes both redundant baseline methods and basic point source calibration methods as special cases, and show how slight deviations from perfect redundancy and coplanarity can be taken into account.

The Terrestrial Planet Finder Coronagraph (TPF-C) mission presented here is an existence proof for a flagship-class internal coronagraph space mission capable of detecting and characterizing Earth-like planets and planetary systems at visible wavelengths around nearby stars, using an existing launch vehicle. TPF-C will use spectroscopy to measure key properties of exoplanets including the presence of atmospheric water or oxygen, powerful signatures in the search for habitable worlds.

21 cm tomography is expected to be difficult in part because of serious foreground contamination. Previous studies have found that line-of-sight approaches are capable of cleaning foregrounds to an acceptable level on large spatial scales, but not on small spatial scales. In this paper, we introduce a Fourier-space formalism for describing the line-of-sight methods, and use it to introduce an improved new method for 21 cm foreground cleaning. Heuristically, this method involves fitting foregrounds in Fourier space using weighted polynomial fits, with each pixel weighted according to its information content. We show that the new method reproduces the old one on large angular scales, and gives marked improvements on small scales at essentially no extra computational cost.

21 cm tomography is emerging as a promising probe of the cosmological dark ages and the epoch of reionization, as well as a tool for observational cosmology in general. However, serious sources of foreground contamination must be subtracted for experimental efforts to be viable. In this paper, we focus on the removal of unresolved extragalactic point sources with smooth spectra, and evaluate how the residual foreground contamination after cleaning depends on instrumental and algorithmic parameters. A crucial but often ignored complication is that the synthesized beam of an interferometer array shrinks towards higher frequency, causing complicated frequency structure in each sky pixel as "frizz" far from the beam center contracts across unresolved radio sources. We find that current-generation experiments should nonetheless be able to clean out this points source contamination adequately, and quantify the instrumental and algorithmic design specifications required to meet this foreground challenge.

We present the period analysis of unfiltered photometric observations of V5116 Sgr (Nova Sgr 2005 #2) and we search for superhump candidates in novae remnants. The PDM method for period analysis is used. The masses of the novae componets are estimated from the secondary mass -- orbital period and primary mass -- decline time relations. We found that 13 nights of V5116 Sgr observations in the year 2006 are modulated with a period of $0.1238 \pm 0.0001$ d ($2.9712 \pm 0.0024$ h). Following the shape of the phased light curves and no apparent change in the value of the periodicity in different subsamples of the data, we interpret the period as orbital in nature. The binary system then falls within the period gap of the orbital period distribution of cataclysmic variables. From the maximum magnitude -- rate of decline relation, we estimate the maximum absolute visual magnitude of $M_{\rm Vmax} = -8.85 \pm 0.04$ mag using the measured value of decline $t_{\rm 2} = 6.5 \pm 1.0$ d. The mass-period relation for cataclysmic variables yields a secondary mass estimate of about $0.26 \pm 0.05 {\rm M}_{\rm \odot}$. We propose that V5116 Sgr is a high inclination system showing an irradiation effect of the secondary star. No fully developed accretion disc up to the tidal radius with the value lower than $3.5 10^{10}$ cm is probable. The mass ratio was estimated in a few novae and the presence or absence of superhumps in these systems was compared with the mass ratio limit for superhumps of about 0.35. We found that in the majority of novae with expected superhumps, this variability has not been found yet. Therefore, more observations of these systems is encouraged.

A period analysis of CCD unfiltered photometry of V4745 Sgr (Nova Sgr 2003 #1) performed during 23 nights in the years 2003 - 2005 is presented. The photometric data are modulated with a period of 0.20782 +- 0.00003 d ($4.98768 +- 0.00072 h). Following the shape of the phased light curve and the presence of the periodicity in all data sets with no apparent change in its value, we interpret this periodicity as orbital in nature and this is consistent with a cataclysmic variable above the period gap. We found a probable short-term periodicity of 0.017238 +- 0.000037 d (24.82272 +- 0.05328 m) which we interpret as the probable spin period of the white dwarf or the beat period between the orbital and spin period. Therefore, we propose that nova V4745 Sgr should be classified as an intermediate polar candidate, supporting the proposed link between transition-oscillation novae and intermediate polars. The mass-period relation for cataclysmic variables yields a secondary mass of about 0.52 +- 0.05 M_solar.

We present the results of 11 nights of CCD unfiltered photometry of V4743 Sgr (Nova Sgr 2002 # 3) from 2003 and 2005. We find two periods of 0.2799 d ~ 6.7 h and 0.01642 d ~ 24 min in the 2005 data. The long period is also present in the 2003 data, but only weak evidence of the shorter period is found in this year. The 24-min period is somewhat longer than the 22-min period, which was detected from X-ray observations. We suggest that the 6.7-h periodicity represents the orbital period of the underlying binary system and that the 24-min period is the beat periodicity between the orbital period and the X-ray period, which is presumably the spin period of the white dwarf. Thus, V4743 Sgr should be classified as an intermediate polar (DQ Her star). About six months after the nova outburst, the optical light curve of V4743 Sgr seemed to show quasi-periodic oscillations, which are typical of the transient phase in classical nova. Therefore, our results support the previous suggestion that the trans ition phase in novae may be related to intermediate polars.

We report on time-resolved photometry of the 2005 July superoutburst of the dwarf nova, 2QZ J021927.9-304545. The resultant light curves showed conspicuous superhumps with a period of 0.081113(19) days, confirming the SU UMa nature of the object. Although we missed the maximum phase of the outburst, the amplitude of the superoutburst well exceeded 5 mag. This value is slightly larger than that of typical SU UMa-type dwarf novae. The superhump period decreased as time elapsed, as can be seen in most SU UMa-type dwarf novae. Based on the archive of ASAS-3, the recurrence time of a superoutburst of the variable turned out to be about 400 days. This value is typical of well known SU UMa stars. The distance to this system was roughly estimated as 370(+20, -60) pc using an empirical relation.

We present the results of 37 nights of CCD unfiltered photometry of nova V2574 Oph (2004) from 2004 and 2005. We find two periods of 0.14164 d (~3.40 h) and 0.14773 d (~3.55 h) in the 2005 data. The 2004 data show variability on a similar timescale, but no coherent periodicity was found. We suggest that the longer periodicity is the orbital period of the underlying binary system and that the shorter period represents a negative superhump. The 3.40 h period is about 4% shorter than the orbital period and obeys the relation between superhump period deficit and binary period. The detection of superhumps in the light curve is evidence of the presence of a precessing accretion disk in this binary system shortly after the nova outburst. From the maximum magnitude - rate of decline relation, we estimate the decay rate t_2 = 17+/-4 d and a maximum absolute visual magnitude of M_Vmax = -7.7+/-1.7 mag.

We present the results of 26 nights of CCD photometry of the nova V2540 Oph (2002) from 2003 and 2004. We find a period of 0.284781 +/- 0.000006 d (6.8347 +/- 0.0001 h) in the data. Since this period was present in the light curves taken in both years, with no apparent change in its value or amplitude, we interpret it as the orbital period of the nova binary system. The mass-period relation for cataclysmic variables yields a secondary mass of about 0.75 +/- 0.04 Msun. From maximum magnitude - rate of decline relation, we estimate a maximum absolute visual magnitude of M(V)=-6.2 +/- 0.4 mag. This value leads to an uncorrected distance modulus of (m-M) = 14.7 +/- 0.7. By using the interstellar reddening for the location of V2540 Oph, we find a rough estimate for the distance of 5.2 +/- 0.8 kpc. We propose that V2540 Oph is either: 1) a high inclination cataclysmic variable showing a reflection effect of the secondary star, or having a spiral structure in the accretion disc, 2) a high inclination intermediate polar system, or less likely 3) a polar.

We present the results obtained from unfiltered photometric CCD observations of V1193 Ori made during 24 nights between November 2002 and January 2003. We found periods of 0.1430 d, 0.1362 d and possibly 2.98 d in the data. The 0.1430 d period is consistent with the 1d^{-1} alias of the proposed orbital period of P_{orb}=0.165 d (Ringwald et al., 1994, Papadaki et al., 2004). Thus and using the known relation between the orbital and superhump periods, we interpret these periods as the orbital period of 0.1430 d, the negative superhump period of 0.1362 d and the precession period of 2.98 d. V1193 Ori can then be classified as a permanent superhump system.

Our previous analyses of radio Doppler and ranging data from distant spacecraft in the solar system indicated that an apparent anomalous acceleration is acting on Pioneer 10 and 11, with a magnitude a_P ~ 8 x 10^{-8} cm/s^2, directed towards the Sun (anderson,moriond). Much effort has been expended looking for possible systematic origins of the residuals, but none has been found. A detailed investigation of effects both external to and internal to the spacecraft, as well as those due to modeling and computational techniques, is provided. We also discuss the methods, theoretical models, and experimental techniques used to detect and study small forces acting on interplanetary spacecraft. These include the methods of radio Doppler data collection, data editing, and data reduction. There is now further data for the Pioneer 10 orbit determination. The extended Pioneer 10 data set spans 3 January 1987 to 22 July 1998. [For Pioneer 11 the shorter span goes from 5 January 1987 to the time of loss of coherent data on 1 October 1990.] With these data sets and more detailed studies of all the systematics, we now give a result, of a_P = (8.74 +/- 1.33) x 10^{-8} cm/s^2. (Annual/diurnal variations on top of a_P, that leave a_P unchanged, are also reported and discussed.)

We present the results obtained from unfiltered photometric CCD observations of the newly discovered cataclysmic variable SDSS J040714.78-064425.1 made during 7 nights in November 2003. We establish the dwarf nova nature of the object as it was in outburst during our observations. We also confirm the presence of deep eclipses with a period of 0.17017d+/-0.00003 in the optical light curve of the star. In addition, we found periods of 0.166d+/-0.001 and possibly also 5.3d+/-0.7 in the data. The 0.17017d periodicity is consistent within the errors with the proposed orbital period of 0.165d (Szkody et. al. 2003) and 0.1700d (Monard 2004). Using the known relation between the orbital and superhump periods, we interpret the 0.166d and 5.3d periods as the negative superhump and the nodal precession period respectively. SDSS J040714.78-064425.1 is then classified as a negative superhump system with one of the largest orbital periods.

We conducted a time-resolved photometric campaign of the helium dwarf nova, 2003aw in 2004 May--June. 2003aw stayed at 14.7--15.7 mag for the first several days during this campaign, which is the plateau phase of this superoutburst. This variable then rapidly decayed to $\sim$18.0 mag, still brighter by about 2 mag than its quiescence magnitude, and maintained this brightness for about 20 days, having short flares of $\sim$2 mag. A long fading tail followed it. We detected superhumps with a period of 0.02357(4) d [= 2036(3) s] during the plateau phase. The whole light curve of the superoutburst in 2003aw, taking into account the present data and those in the literature, perfectly resembles that of the 1996-1997 superoutburst of the peculiar WZ Sge-type hydrogen-rich dwarf nova, EG Cnc.

Since the discovery of the largest positive superhump period in TV Col, we have started a program to search for superhumps in CVs with large orbital periods. Here, we summarize preliminary results of TX Col and V4742 Sgr. TX Col is an intermediate polar with a 5.7-h orbital period. V4742 Sgr is a recent nova with no known periods. CCD unfiltered continuous photometry of these 2 objects was carried out during 56 nights in 2002-3. In TX Col, in addition to the orbital period of 5.7 h, we found peaks at 7.1 h and 5.0 h. These are interpreted as positive and negative superhumps correspondingly, although the effects of the quasi-periodic oscillations at about 2 h were not taken into consideration. In the light curve of V4742 Sgr 2 long periods are detected -- 6.1 and 5.4 h as well as a short-term period at 1.6 h. This result suggests that V4742 Sgr is an intermediate polar candidate and a permanent superhump system with a large orbital period (5.4 h) and a superhump period excess of 13 percent. If these results are confirmed, TX Col, V4742 Sgr and TV Col form a group of intermediate polars with extremely large superhump periods. There seems to be now growing evidence that superhumps can occur in intermediate polars with long orbital periods, which is very likely inconsistent with the theoretical prediction that superhumps can only occur in systems with mass ratios below 0.33. Alternatively, if the mass ratio in these systems is nevertheless below the theoretical limit, they should harbour undermassive secondaries and massive white dwarfs, near the Chandrasekhar limit, which would make them excellent candidates for progenitors of supernovae type Ia.

Rapid, efficient deactivation of the photoexcited PAH cations accounts for their remarkable photostability and have important implications for astrochemistry, as these cations are the leading candidates for the species responsible for the diffuse interstellar bands (DIB) observed throughout the Galaxy.Ultrafast relaxation dynamics for photoexcited PAH cations isolated in boric acid glass have been studied using femtosecond and picosecond transient grating spectroscopy. With the exception of perylene+, the recovery kinetics for the ground doublet (D0) states of these radical cations are biexponential, containing a fast (< 200 fs) and a slow (3-20 ps) components. No temperature dependence or isotope effect was observed for the fast component, whereas the slow component exhibits both the H/D isotope effect (1.1-1.3) and strong temperature dependence (15 to 300 K). We suggest that the fast component is due to internal Dn to D0 conversion and the slow component is due to vibrational energy transfer (VET) from a hot D0 state to the glass matrix.

We conclude that Katz's proposal (anisotropic heat reflection off of the back of the spacecraft high-gain antennae, the heat coming from the RTGs) does not provide enough power and so can not explain the Pioneer anomaly.

We conclude that Murphy's proposal (radiation of the power of the main-bus electrical systems from the rear of the craft) can not explain the anomalous Pioneer acceleration.

Recently we reported that radio Doppler data generated by NASA's Deep Space Network (DSN) from the Pioneer 10 and 11 spacecraft indicate an apparent anomalous, constant, spacecraft acceleration with a magnitude $\sim 8.5\times 10^{-8}$ cm s$^{-2}$, directed towards the Sun (gr-qc/9808081). Analysis of similar Doppler and ranging data from the Galileo and Ulysses spacecraft yielded ambiguous results for the anomalous acceleration, but it was useful in that it ruled out the possibility of a systematic error in the DSN Doppler system that could easily have been mistaken as a spacecraft acceleration. Here we present some new results, including a critique suggestions that the anomalous acceleration could be caused by collimated thermal emission. Based partially on a further data for the Pioneer 10 orbit determination, the data now spans January 1987 to July 1998, our best estimate of the average Pioneer 10 acceleration directed towards the Sun is $\sim 7.5 \times 10^{-8}$ cm s$^{-2}$.

Radio metric data from the Pioneer 10/11, Galileo, and Ulysses spacecraft indicate an apparent anomalous, constant, acceleration acting on the spacecraft with a magnitude $\sim 8.5\times 10^{-8}$ cm/s$^2$, directed towards the Sun. Two independent codes and physical strategies have been used to analyze the data. A number of potential causes have been ruled out. We discuss future kinematic tests and possible origins of the signal.