Line intensity mapping (LIM) is a promising tool to efficiently probe the three-dimensional large-scale structure by mapping the aggregate emission of a spectral line from all sources that trace the matter density field. Spectral lines from different redshifts can fall in the same observed frequency and be confused, however, which is a major challenge in LIM. In this work, we develop a line de-confusion technique in map space capable of reconstructing the three-dimensional spatial distribution of line-emitting sources. If multiple spectral lines of a source population are observable in multiple frequencies, using the sparse approximation our technique iteratively extracts sources along a given line-of-sight by fitting the LIM data to a set of spectral templates. We demonstrate that the technique successfully extracts sources with emission lines present at a few sigma above the noise level, taking into account uncertainties in the source modeling and presence of continuum foreground contamination and noise fluctuations. As an example, we consider a TIME/CONCERTO-like survey targeting [C II] at the Epoch of Reionization, and reliably reconstruct the 3D spatial distribution of the CO interlopers at $0.5\lesssim z\lesssim 1.5$. We also demonstrate a successful de-confusion for the SPHEREx mission in the near-infrared wavelengths. Potentially, the reconstructed maps can be further cross-correlated with a (galaxy) tracer population to estimate the total interloper power in the linear clustering regime. This technique is a general framework to extract the phase-space distribution of low-redshift interlopers, without the need of external information, for any line de-confusion problem.

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 present a formalism that incorporates hydrogen Lyman-alpha (Ly$\alpha$) polarization arising from the scattering of radiation in galaxy halos into the intensity mapping approach. Using the halo model, and emission profiles based on simulations and observations, we calculate auto and cross power spectra at redshifts $3\leq z \leq13$ for the Ly$\alpha$ total intensity, $I$, polarized intensity, $P$, degree of polarization, $\Pi=P/I$, and two new quantities, the astrophysical $E$ and $B$ modes of Ly$\alpha$ polarization. The one-halo terms of the $\Pi$ power spectra show a turnover that signals the average extent of the polarization signal, and thus the extent of the scattering medium. The position of this feature depends on redshift, and on the specific emission profile shape and extent, in our formalism. Therefore, the comparison of various Ly$\alpha$ polarization quantities and redshifts can break degeneracies between competing effects, and it can reveal the true shape of the emission profiles, and thus, the physical properties of the cool gas in halos. Furthermore, measurements of Ly$\alpha$ $E$ and $B$ modes may be used as probes of galaxy evolution, because they are related to the average degree of anisotropy in the emission and in the halo gas distribution across redshifts. The detection of the polarization signal at $z \sim 3-5$ requires improvements in the sensitivity of current ground-based experiments by a factor of $\sim 10$, and of $\sim 100$ for space-based instruments targeting the redshifts $z\sim 9-10$, the exact values depending on the specific redshift and experiment. Interloper contamination in polarization is expected to be small, because the interlopers need to also be polarized. Overall, Ly$\alpha$ polarization boosts the amount of physical information retrievable on galaxies and their surroundings, most of it not achievable with total emission alone. (abridged)

We use a semi-analytic model to explore the potential impact of a brief and violent period of radio-loud accretion onto black-holes (The Radio Scream) during the Cosmic Dawn on the HI hyperfine 21 cm signal. We find that radio emission from super-massive black hole seeds can impact the global 21 cm signal at the level of tens to hundreds of percent provided that they were as radio loud as $z\approx1$ black holes and obscured by gas with column depths of $N_\text{H}\gtrsim 10^{23}$ cm$^{-2}$. We determine plausible sets of parameters that reproduce some of the striking features of the EDGES absorption feature including its depth, timing, and side steepness while producing radio/X-ray backgrounds and source counts that are consistent with published limits. Scenarios yielding a dramatic 21 cm signature also predict large populations of $\sim \mu$Jy point sources that will be detectable in future deep surveys from the Square Kilometer Array (SKA). Thus, 21 cm measurements, complemented by deep point source surveys, have the potential to constrain optimistic scenarios where super-massive black-hole progenitors were radio-loud.

Line intensity mapping (LIM) is a promising approach to study star formation and the interstellar medium (ISM) in galaxies by measuring the aggregate line emission from the entire galaxy population. In this work, we develop a simple yet physically-motivated framework for modeling the line emission as would be observed in LIM experiments, which is done by building on analytic models of the Cosmic Infrared Background that connect total infrared luminosity of galaxies to their host dark matter halos. We present models of the HI 21cm, CO(1-0), [C II] 158$\mu$m, and [N II] 122 and 205$\mu$m lines consistent with current observational constraints. With four case studies of various combinations of these lines that probe different ISM phases, we demonstrate the potential for reliably extracting physical properties of the ISM, and the evolution of these properties with cosmic time, from auto and cross-correlation analysis of these lines as measured by future LIM experiments.

This is the final report submitted to NASA for a Probe-class concept study of the "Farside Array for Radio Science Investigations of the Dark ages and Exoplanets" (FARSIDE), a low radio frequency interferometric array on the farside of the Moon. The design study focused on the instrument, a deployment rover, the lander and base station, and delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dipole antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway, or an alternative relay satellite. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 150 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts from 50-100. It will also be a pathfinder for a larger 21-cm power spectrum instrument by carefully measuring the foreground with high dynamic range.

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.

FARSIDE (Farside Array for Radio Science Investigations of the Dark ages and Exoplanets) is a Probe-class concept to place a low radio frequency interferometric array on the farside of the Moon. A NASA-funded design study, focused on the instrument, a deployment rover, the lander and base station, delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dual polarization antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 100 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets, including the hypothetical Planet Nine. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts z=50-100. The unique observational window offered by FARSIDE would enable an abundance of additional science ranging from sounding of the lunar subsurface to characterization of the interstellar medium in the solar system neighborhood.

Intensity mapping has emerged as a promising tool to probe the three-dimensional structure of the universe. The traditional approach of galaxy redshift surveys is based on individual galaxy detection, typically performed by thresholding and digitizing large-scale intensity maps. By contrast, intensity mapping uses the integrated emission from all sources in a 3D pixel (or voxel) as an analog tracer of large-scale structure. In this work, we develop a formalism to quantify the performance of both approaches when measuring large-scale structures. We compute the Fisher information of an arbitrary observable, derive the optimal estimator, and study its performance as a function of source luminosity function, survey resolution, instrument sensitivity, and other survey parameters. We identify regimes where each approach is advantageous and discuss optimal strategies for different scenarios. To determine the best strategy for any given survey, we develop a metric that is easy to compute from the source luminosity function and the survey sensitivity, and we demonstrate the application with several planned intensity mapping surveys.

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 expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

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

The Cosmic Dawn Intensity Mapper (CDIM) will transform our understanding of the era of reionization when the Universe formed the first stars and galaxies, and UV photons ionized the neutral medium. CDIM goes beyond the capabilities of upcoming facilities by carrying out wide area spectro-imaging surveys, providing redshifts of galaxies and quasars during reionization as well as spectral lines that carry crucial information on their physical properties. CDIM will make use of unprecedented sensitivity to surface brightness to measure the intensity fluctuations of reionization on large-scales to provide a valuable and complementary dataset to 21-cm experiments. The baseline mission concept is an 83-cm infrared telescope equipped with a focal plane of 24 \times 20482 detectors capable of R = 300 spectro-imaging observations over the wavelength range of 0.75 to 7.5 {\mu}m using Linear Variable Filters (LVFs). CDIM provides a large field of view of 7.8 deg2 allowing efficient wide area surveys, and instead of moving instrumental components, spectroscopic mapping is obtained through a shift-and-stare strategy through spacecraft operations. CDIM design and capabilities focus on the needs of detecting faint galaxies and quasars during reionization and intensity fluctuation measurements of key spectral lines, including Lyman-{\alpha} and H{\alpha} radiation from the first stars and galaxies. The design is low risk, carries significant science and engineering margins, and makes use of technologies with high technical readiness level for space observations.

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.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Line-Intensity Mapping is an emerging technique which promises new insights into the evolution of the Universe, from star formation at low redshifts to the epoch of reionization and cosmic dawn. It measures the integrated emission of atomic and molecular spectral lines from galaxies and the intergalactic medium over a broad range of frequencies, using instruments with aperture requirements that are greatly relaxed relative to surveys for single objects. A coordinated, comprehensive, multi-line intensity-mapping experimental effort can efficiently probe over 80% of the volume of the observable Universe - a feat beyond the reach of other methods. Line-intensity mapping will uniquely address a wide array of pressing mysteries in galaxy evolution, cosmology, and fundamental physics. Among them are the cosmic history of star formation and galaxy evolution, the compositions of the interstellar and intergalactic media, the physical processes that take place during the epoch of reionization, cosmological inflation, the validity of Einstein's gravity theory on the largest scales, the nature of dark energy and the origin of dark matter.

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

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

Upcoming imaging surveys, such as LSST, will provide an unprecedented view of the Universe, but with limited resolution along the line-of-sight. Common ways to increase resolution in the third dimension, and reduce misclassifications, include observing a wider wavelength range and/or combining the broad-band imaging with higher spectral resolution data. The challenge with these approaches is matching the depth of these ancillary data with the original imaging survey. However, while a full 3D map is required for some science, there are many situations where only the statistical distribution of objects (dN/dz) in the line-of-sight direction is needed. In such situations, there is no need to measure the fluxes of individual objects in all of the surveys. Rather a stacking procedure can be used to perform an `ensemble photo-z'. We show how a shallow, higher spectral resolution survey can be used to measure dN/dz for stacks of galaxies which coincide in a deeper, lower resolution survey. The galaxies in the deeper survey do not even need to appear individually in the shallow survey. We give a toy model example to illustrate tradeoffs and considerations for applying this method. This approach will allow deep imaging surveys to leverage the high resolution of spectroscopic and narrow/medium band surveys underway, even when the latter do not have the same reach to high redshift.

We estimate the 21 cm Radio Background from accretion onto the first intermediate-mass Black Holes between $z\approx 30$ and $z\approx 16$. Combining potentially optimistic, but plausible, scenarios for black hole formation and growth with empirical correlations between luminosity and radio-emission observed in low-redshift active galactic nuclei, we find that a model of black holes forming in molecular cooling halos is able to produce a 21 cm background that exceeds the Cosmic Microwave Background (CMB) at $z \approx 17$ though models involving larger halo masses are not entirely excluded. Such a background could explain the surprisingly large amplitude of the 21 cm absorption feature recently reported by the EDGES collaboration. Such black holes would also produce significant X-ray emission and contribute to the $0.5-2$ keV soft X-ray background at the level of $\approx 10^{-13}-10^{-12}$ erg sec$^{-1}$ cm$^{-2}$ deg$^{-2}$, consistent with existing constraints. In order to avoid heating the IGM over the EDGES trough, these black holes would need to be obscured by Hydrogen column depths of $ N_\text{H} \sim 5 \times 10^{23} \text{cm}^{-2}$. Such black holes would avoid violating contraints on the CMB optical depth from Planck if their UV photon escape fractions were below $f_{\text{esc}} \lesssim 0.1$, which would be a natural result of $N_\text{H} \sim 5 \times 10^{23} \text{cm}^{-2}$ imposed by an unheated IGM.

The near-infrared background between 0.5 $\mu$m to 2 $\mu$m contains a wealth of information related to radiative processes in the universe. Infrared background anisotropies encode the redshift-weighted total emission over cosmic history, including any spatially diffuse and extended contributions. The anisotropy power spectrum is dominated by undetected galaxies at small angular scales and diffuse background of Galactic emission at large angular scales. In addition to these known sources, the infrared background also arises from intra-halo light (IHL) at $z < 3$ associated with tidally-stripped stars during galaxy mergers. Moreover, it contains information on the very first galaxies from the epoch of reionization (EoR). The EoR signal has a spectral energy distribution (SED) that goes to zero near optical wavelengths due to Lyman absorption, while other signals have spectra that vary smoothly with frequency. Due to differences in SEDs and spatial clustering, these components may be separated in a multi-wavelength-fluctuation experiment. To study the extent to which EoR fluctuations can be separated in the presence of IHL, extra-galactic and Galactic foregrounds, we develop a maximum likelihood technique that incorporates a full covariance matrix among all the frequencies at different angular scales. We apply this technique to simulated deep imaging data over a 2$\times$100 deg$^2$ sky area from 0.75 $\mu$m to 5 $\mu$m in 9 bands and find that such a "frequency tomography" can successfully reconstruct both the amplitude and spectral shape for representative EoR, IHL and the foreground signals.

We present a search for CII emission over cosmological scales at high-redshifts. The CII line is a prime candidate to be a tracer of star formation over large-scale structure since it is one of the brightest emission lines from galaxies. Redshifted CII emission appears in the submillimeter regime, meaning it could potentially be present in the higher frequency intensity data from the Planck satellite used to measure the cosmic infrared background (CIB). We search for CII emission over redshifts z=2-3.2 in the Planck 545 GHz intensity map by cross-correlating the 3 highest frequency Planck maps with spectroscopic quasars and CMASS galaxies from the Sloan Digital Sky Survey III (SDSS-III), which we then use to jointly fit for CII intensity, CIB parameters, and thermal Sunyaev-Zeldovich (SZ) emission. We report a measurement of an anomalous emission $\mathrm{I_\nu}=6.6^{+5.0}_{-4.8}\times10^4$ Jy/sr at 95% confidence, which could be explained by CII emission, favoring collisional excitation models of CII emission that tend to be more optimistic than models based on CII luminosity scaling relations from local measurements; however, a comparison of Bayesian information criteria reveal that this model and the CIB & SZ only model are equally plausible. Thus, more sensitive measurements will be needed to confirm the existence of large-scale CII emission at high redshifts. Finally, we forecast that intensity maps from Planck cross-correlated with quasars from the Dark Energy Spectroscopic Instrument (DESI) would increase our sensitivity to CII emission by a factor of 5, while the proposed Primordial Inflation Explorer (PIXIE) could increase the sensitivity further.

We review three distance measurement techniques beyond the local universe: (1) gravitational lens time delays, (2) baryon acoustic oscillation (BAO), and (3) HI intensity mapping. We describe the principles and theory behind each method, the ingredients needed for measuring such distances, the current observational results, and future prospects. Time delays from strongly lensed quasars currently provide constraints on $H_0$ with < 4% uncertainty, and with 1% within reach from ongoing surveys and efforts. Recent exciting discoveries of strongly lensed supernovae hold great promise for time-delay cosmography. BAO features have been detected in redshift surveys up to z <~ 0.8 with galaxies and z ~ 2 with Ly-$\alpha$ forest, providing precise distance measurements and $H_0$ with < 2% uncertainty in flat $\Lambda$CDM. Future BAO surveys will probe the distance scale with percent-level precision. HI intensity mapping has great potential to map BAO distances at z ~ 0.8 and beyond with precisions of a few percent. The next years ahead will be exciting as various cosmological probes reach 1% uncertainty in determining $H_0$, to assess the current tension in $H_0$ measurements that could indicate new physics.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

Intensity mapping provides a unique means to probe the epoch of reionization (EoR), when the neutral intergalactic medium was ionized by the energetic photons emitted from the first galaxies. The [CII] 158$\mu$m fine-structure line is typically one of the brightest emission lines of star-forming galaxies and thus a promising tracer of the global EoR star-formation activity. However, [CII] intensity maps at $6 \lesssim z \lesssim 8$ are contaminated by interloping CO rotational line emission ($3 \leq J_{\rm upp} \leq 6$) from lower-redshift galaxies. Here we present a strategy to remove the foreground contamination in upcoming [CII] intensity mapping experiments, guided by a model of CO emission from foreground galaxies. The model is based on empirical measurements of the mean and scatter of the total infrared luminosities of galaxies at $z < 3$ and with stellar masses $M_{*} > 10^{8}\,\rm M_{\rm \odot}$ selected in $K$-band from the COSMOS/UltraVISTA survey, which can be converted to CO line strengths. For a mock field of the Tomographic Ionized-carbon Mapping Experiment (TIME), we find that masking out the "voxels" (spectral-spatial elements) containing foreground galaxies identified using an optimized CO flux threshold results in a $z$-dependent criterion $m^{\rm AB}_{\rm K} \lesssim 22$ (or $M_{*} \gtrsim 10^{9} \,\rm M_{\rm \odot}$) at $z < 1$ and makes a [CII]/CO$_{\rm tot}$ power ratio of $\gtrsim 10$ at $k=0.1$ $h$/Mpc achievable, at the cost of a moderate $\lesssim 8\%$ loss of total survey volume.

We report results from 21-cm intensity maps acquired from the Parkes radio telescope and cross-correlated with galaxy maps from the 2dF galaxy survey. The data span the redshift range $0.057<z<0.098$ and cover approximately 1,300 square degrees over two long fields. Cross correlation is detected at a significance of $5.18\sigma$. The amplitude of the cross-power spectrum is low relative to the expected dark matter power spectrum, assuming a neutral hydrogen (HI) bias and mass density equal to measurements from the ALFALFA survey. The decrement is pronounced and statistically significant at small scales. At $k\sim1.5$ $ h \mathrm{Mpc^{-1}}$, the cross power spectrum is more than a factor of 6 lower than expected, with a significance of $14.8\,\sigma$. This decrement indicates either a lack of clustering of neutral hydrogen (HI), a small correlation coefficient between optical galaxies and HI, or some combination of the two. Separating 2dF into red and blue galaxies, we find that red galaxies are much more weakly correlated with HI on $k\sim1.5$ $h \mathrm{Mpc^{-1}}$ scales, suggesting that HI is more associated with blue star-forming galaxies and tends to avoid red galaxies.

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 present the first application of a new foreground removal pipeline to the current leading HI intensity mapping dataset, obtained by the Green Bank Telescope (GBT). We study the 15hr and 1hr field data of the GBT observations previously presented in Masui et al (2013) and Switzer et al (2013), covering about 41 square degrees at 0.6<z<1.0, for which cross-correlations may be measured with the galaxy distribution of the WiggleZ Dark Energy Survey. In the presented pipeline, we subtract the Galactic foreground continuum and the point source contamination using an independent component analysis technique (fastica), and develop a Fourier-based optimal estimator to compute the temperature power spectrum of the intensity maps and cross-correlation with the galaxy survey data. We show that fastica is a reliable tool to subtract diffuse and point-source emission through the non-Gaussian nature of their probability distributions. The temperature power spectra of the intensity maps is dominated by instrumental noise on small scales which fastica, as a conservative subtraction technique of non-Gaussian signals, can not mitigate. However, we determine similar GBT-WiggleZ cross-correlation measurements to those obtained by the Singular Value Decomposition (SVD) method, and confirm that foreground subtraction with fastica is robust against 21cm signal loss, as seen by the converged amplitude of these cross-correlation measurements. We conclude that SVD and fastica are complementary methods to investigate the foregrounds and noise systematics present in intensity mapping datasets.

Polarization leakage of foreground synchrotron emission is a critical issue in HI intensity mapping experiments. While the sought-after HI emission is unpolarized, polarized foregrounds such as Galactic and extragalactic synchrotron radiation, if coupled with instrumental impurity, can mimic or overwhelm the HI signals. In this paper we present the methodology for polarization calibration at 700-900 MHz, applied on data obtained from the Green Bank Telescope (GBT). We use astrophysical sources, both polarized and unpolarized sources including quasars and pulsars, as calibrators to characterize the polarization leakage and control systematic effects in our GBT HI intensity mapping project. The resulting fractional errors on polarization measurements on boresight are well controlled to within 0.6%-0.8% of their total intensity. The polarized beam patterns are measured by performing spider scans across both polarized quasars and pulsars. A dominant Stokes I to V leakage feature and secondary features of Stokes I to Q and I to U leakages in the 700-900 MHz frequency range are identified. These characterizations are important for separating foreground polarization leakage from the HI 21 cm signal.

Spectral line intensity mapping has been proposed as a promising tool to efficiently probe the cosmic reionization and the large-scale structure. Without detecting individual sources, line intensity mapping makes use of all available photons and measures the integrated light in the source confusion limit, to efficiently map the three-dimensional matter distribution on large scales as traced by a given emission line. One particular challenge is the separation of desired signals from astrophysical continuum foregrounds and line interlopers. Here we present a technique to extract large-scale structure information traced by emission lines from different redshifts, embedded in a three-dimensional intensity mapping data cube. The line redshifts are distinguished by the anisotropic shape of the power spectra when projected onto a common coordinate frame. We consider the case where high-redshift [CII] lines are confused with multiple low-redshift CO rotational lines. We present a semi-analytic model for [CII] and CO line estimates based on the cosmic infrared background measurements, and show that with a modest instrumental noise level and survey geometry, the large-scale [CII] and CO power spectrum amplitudes can be successfully extracted from a confusion-limited data set, without external information. We discuss the implications and limits of this technique for possible line intensity mapping experiments.

The sky-averaged, or global, background of redshifted $21$ cm radiation is expected to be a rich source of information on cosmological reheating and reionizaton. However, measuring the signal is technically challenging: one must extract a small, frequency-dependent signal from under much brighter spectrally smooth foregrounds. Traditional approaches to study the global signal have used single antennas, which require one to calibrate out the frequency-dependent structure in the overall system gain (due to internal reflections, for example) as well as remove the noise bias from auto-correlating a single amplifier output. This has motivated proposals to measure the signal using cross-correlations in interferometric setups, where additional calibration techniques are available. In this paper we focus on the general principles driving the sensitivity of the interferometric setups to the global signal. We prove that this sensitivity is directly related to two characteristics of the setup: the cross-talk between readout channels (i.e. the signal picked up at one antenna when the other one is driven) and the correlated noise due to thermal fluctuations of lossy elements (e.g. absorbers or the ground) radiating into both channels. Thus in an interferometric setup, one cannot suppress cross-talk and correlated thermal noise without reducing sensitivity to the global signal by the same factor -- instead, the challenge is to characterize these effects and their frequency dependence. We illustrate our general theorem by explicit calculations within toy setups consisting of two short dipole antennas in free space and above a perfectly reflecting ground surface, as well as two well-separated identical lossless antennas arranged to achieve zero cross-talk.

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a new 400-800MHz radio interferometer under development for deployment in South Africa. HIRAX will comprise 1024 six meter parabolic dishes on a compact grid and will map most of the southern sky over the course of four years. HIRAX has two primary science goals: to constrain Dark Energy and measure structure at high redshift, and to study radio transients and pulsars. HIRAX will observe unresolved sources of neutral hydrogen via their redshifted 21-cm emission line (`hydrogen intensity mapping'). The resulting maps of large-scale structure at redshifts 0.8-2.5 will be used to measure Baryon Acoustic Oscillations (BAO). HIRAX will improve upon current BAO measurements from galaxy surveys by observing a larger cosmological volume (larger in both survey area and redshift range) and by measuring BAO at higher redshift when the expansion of the universe transitioned to Dark Energy domination. HIRAX will complement CHIME, a hydrogen intensity mapping experiment in the Northern Hemisphere, by completing the sky coverage in the same redshift range. HIRAX's location in the Southern Hemisphere also allows a variety of cross-correlation measurements with large-scale structure surveys at many wavelengths. Daily maps of a few thousand square degrees of the Southern Hemisphere, encompassing much of the Milky Way galaxy, will also open new opportunities for discovering and monitoring radio transients. The HIRAX correlator will have the ability to rapidly and eXperimentciently detect transient events. This new data will shed light on the poorly understood nature of fast radio bursts (FRBs), enable pulsar monitoring to enhance long-wavelength gravitational wave searches, and provide a rich data set for new radio transient phenomena searches. This paper discusses the HIRAX instrument, science goals, and current status.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

Intensity mapping experiments survey the spectrum of diffuse line radiation rather than detect individual objects at high signal-to-noise. Spectral maps of unresolved atomic and molecular line radiation contain three-dimensional information about the density and environments of emitting gas, and efficiently probe cosmological volumes out to high redshift. Intensity mapping survey volumes also contain all other sources of radiation at the frequencies of interest. Continuum foregrounds are typically ~10^2-10^3 times brighter than the cosmological signal. The instrumental response to bright foregrounds will produce new spectral degrees of freedom that are not known in advance, nor necessarily spectrally smooth. The intrinsic spectra of foregrounds may also not be well-known in advance. We describe a general class of quadratic estimators to analyze data from single-dish intensity mapping experiments, and determine contaminated spectral modes from the data itself. The key attribute of foregrounds is not that they are spectrally smooth, but instead that they have fewer bright spectral degrees of freedom than the cosmological signal. Spurious correlations between the signal and foregrounds produce additional bias. Compensation for signal attenuation must estimate and correct this bias. A successful intensity mapping experiment will control instrumental systematics that spread variance into new modes, and it must observe a large enough volume that contaminant modes can be determined independently from the signal on scales of interest.

Cosmic Dawn Intensity Mapper is a "Probe Class" mission concept for reionization studies of the universe. It will be capable of spectroscopic imaging observations between 0.7 to 6-7 microns in the near-Infrared. The primary observational objective is pioneering observations of spectral emission lines of interest throughout the cosmic history, but especially from the first generation of distant, faint galaxies when the universe was less than 800 million years old. With spectro-imaging capabilities, using a set of linear variable filters (LVFs), CDIM will produce a three-dimensional tomographic view of the epoch of reionization (EoR). CDIM will also study galaxy formation over more than 90% of the cosmic history and will move the astronomical community from broad-band astronomical imaging to low-resolution (R=200-300) spectro-imaging of the universe.

Acoustic waves traveling through the early Universe imprint a characteristic scale in the clustering of galaxies, QSOs and inter-galactic gas. This scale can be used as a standard ruler to map the expansion history of the Universe, a technique known as Baryon Acoustic Oscillations (BAO). BAO offer a high-precision, low-systematics means of constraining our cosmological model. The statistical power of BAO measurements can be improved if the `smearing' of the acoustic feature by non-linear structure formation is undone in a process known as reconstruction. In this paper we use low-order Lagrangian perturbation theory to study the ability of $21\,$cm experiments to perform reconstruction and how augmenting these surveys with galaxy redshift surveys at relatively low number densities can improve performance. We find that the critical number density which must be achieved in order to benefit $21\,$cm surveys is set by the linear theory power spectrum near its peak, and corresponds to densities achievable by upcoming surveys of emission line galaxies such as eBOSS and DESI. As part of this work we analyze reconstruction within the framework of Lagrangian perturbation theory with local Lagrangian bias, redshift-space distortions, ${\bf k}$-dependent noise and anisotropic filtering schemes.

Fast Radio Bursts are bright, unresolved, non-repeating, broadband, millisecond flashes, found primarily at high Galactic latitudes, with dispersion measures much larger than expected for a Galactic source. The inferred all-sky burst rate is comparable to the core-collapse supernova rate out to redshift 0.5. If the observed dispersion measures are assumed to be dominated by the intergalactic medium, the sources are at cosmological distances with redshifts of 0.2 to 1. These parameters are consistent with a wide range of source models. One fast radio burst showed circular polarization [21(7)%] of the radio emission, but no linear polarization was detected, and hence no Faraday rotation measure could be determined. Here we report the examination of archival data revealing Faraday rotation in a newly detected burst - FRB 110523. It has radio flux at least 0.6 Jy and dispersion measure 623.30(5) pc cm$^{-3}$. Using Galactic contribution 45 pc cm$^{-3}$ and a model of intergalactic electron density, we place the source at a maximum redshift of 0.5. The burst has rotation measure -186.1(1.4) rad m$^{-2}$, much higher than expected for this line of sight through the Milky Way and the intergalactic medium, indicating magnetization in the vicinity of the source itself or within a host galaxy. The pulse was scattered by two distinct plasma screens during propagation, which requires either a dense nebula associated with the source or a location within the central region of its host galaxy. Keeping in mind that there may be more than one type of fast radio burst source, the detection in this instance of source-local magnetization and scattering favours models involving young stellar populations such as magnetars over models involving the mergers of older neutron stars, which are more likely to be located in low density regions of the host galaxy.

We present constraints on the abundance of carbon-monoxide in the early Universe from the CO Power Spectrum Survey (COPSS). We utilize a data set collected between 2005 and 2008 using the Sunyaev-Zel'dovich Array (SZA), which were previously used to measure arcminute-scale fluctuations of the CMB. This data set features observations of 44 fields, covering an effective area of 1.7 square degrees, over a frequency range of 27 to 35 GHz. Using the technique of intensity mapping, we are able to probe the CO(1-0) transition, with sensitivity to spatial modes between $k=0.5{-}2\ h\,\textrm{Mpc}^{-1}$ over a range in redshift of $z=2.3{-}3.3$, spanning a comoving volume of $3.6\times10^{6}\ h^{-3}\,\textrm{Mpc}^{3}$. We demonstrate our ability to mitigate foregrounds, and present estimates of the impact of continuum sources on our measurement. We constrain the CO power spectrum to $P_{\textrm{CO}}<2.6\times10^{4}\ \mu\textrm{K}^{2} (h^{-1}\,\textrm{Mpc})^{3}$, or $\Delta^{2}_{\textrm{CO}}(k\! = \! 1 \ h\,\textrm{Mpc}^{-1})<1.3 \times10^{3}\ \mu\textrm{K}^{2}$, at $95\%$ confidence. This limit resides near optimistic predictions for the CO power spectrum. Under the assumption that CO emission is proportional to halo mass during bursts of active star formation, this corresponds to a limit on the ratio of $\textrm{CO}(1{-}0)$ luminosity to host halo mass of $A_{\textrm{CO}}<1.2\times10^{-5}\ L_{\odot}\ M_{\odot}^{-1}$. Further assuming a Milky Way-like conversion factor between CO luminosity and molecular gas mass ($\alpha_{\textrm{CO}}=4.3\ M_{\odot}\ (\textrm{K}\ \textrm{km}\ \textrm{s}^{-1}\ \textrm{pc}^{-2})^{-1}$), we constrain the global density of molecular gas to $\rho_{z\sim3}(M_{\textrm{H}_{2}})\leq 2.8 \times10^{8}\ M_{\odot}\ \textrm{Mpc}^{-3}$.

We report here on key science topics for the Next Generation Very Large Array in the areas of time domain, fundamental physics, and cosmology. Key science cases considered are pulsars in orbit around the Galactic Center massive black hole, Sagittarius A*, electromagnetic counterparts to gravitational waves, and astrometric cosmology. These areas all have the potential for ground-breaking and transformative discovery. Numerous other topics were discussed during the preparation of this report and some of those discussions are summarized here, as well. There is no doubt that further investigation of the science case will reveal rich and compelling opportunities.

It has been suggested that antiscreening effects due to the running of the gravitational constant G might provide a partial solution to the dark matter mystery. It has also been hypothesized that renormalization group scaling transformations at large scales might supply the theoretical explanation. In this letter, we demonstrate that multifractal coarse-graining scaling effects due to classical fluctuations in the IR with consecutive symmetry breakings in gravitational evolution and induced running of the gravitational constant with fractal structures at larger scales may provide the plausible explanation of the observed results of weak lensing observations and beyond. The sporadic and localized antiscreening due to the running of the gravitational constant can also provide the backbone for the cosmic evolution and large scale structure formation. Our interpretation of this interesting finding is that such effects are the result of the complexity phenomenon involving the evolution of large-scale multifractal structures and accompanying fluctuations, not the conventional arguments suggesting quantum gravity being the primary cause. We also suggest that the running of the cosmological constant due to such stochastic complexity evolution may provide a key to the understanding of the observed cosmic acceleration.

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

We present the science enabled by cross-correlations of the SKA1-LOW 21-cm EoR surveys with other line mapping programs. In particular, we identify and investigate potential synergies with planned programs, such as the line intensity mapping of redshifted CO rotational lines, [CII] and Ly-$\alpha$ emissions during reionization. We briefly describe how these tracers of the star-formation rate at $z \sim 8$ can be modeled jointly before forecasting their auto- and cross-power spectra measurements with the nominal 21cm EoR survey. The use of multiple line tracers would be invaluable to validate and enrich our understanding of the EoR.

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 6 - 27 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2.

21cm intensity mapping experiments aim to observe the diffuse neutral hydrogen (HI) distribution on large scales which traces the Cosmic structure. The Square Kilometre Array (SKA) will have the capacity to measure the 21cm signal over a large fraction of the sky. However, the redshifted 21cm signal in the respective frequencies is faint compared to the Galactic foregrounds produced by synchrotron and free-free electron emission. In this article, we review selected foreground subtraction methods suggested to effectively separate the 21cm signal from the foregrounds with intensity mapping simulations or data. We simulate an intensity mapping experiment feasible with SKA phase 1 including extragalactic and Galactic foregrounds. We give an example of the residuals of the foreground subtraction with a independent component analysis and show that the angular power spectrum is recovered within the statistical errors on most scales. Additionally, the scale of the Baryon Acoustic Oscillations is shown to be unaffected by foreground subtraction.

One of the main goals of modern observational cosmology is to map the large scale structure of the Universe. A potentially powerful approach for doing this would be to exploit three-dimensional spectral maps, i.e. the specific intensity of extragalactic light as a function of wavelength and direction on the sky, to measure spatial variations in the total extragalactic light emission and use these as a tracer of the clustering of matter. A main challenge is that the observed intensity as a function of wavelength is a convolution of the source luminosity density with the rest-frame spectral energy distribution. In this paper, we introduce the method of spectral deconvolution as a way to invert this convolution and extract the clustering information. We show how one can use observations of the mean and angular fluctuations of extragalactic light as a function of wavelength, assuming statistical isotropy, to reconstruct jointly the rest-frame spectral energy distribution of the sources and the source spatial density fluctuations. This method is more general than the well known line mapping technique as it does not rely on spectral lines in the emitted spectra. After introducing the general formalism, we discuss its implementation and limitations. This formal paper sets the stage for future more practical studies.

In this report, we present a wide variety of ways in which information from multiple probes of dark energy may be combined to obtain additional information not accessible when they are considered separately. Fundamentally, because all major probes are affected by the underlying distribution of matter in the regions studied, there exist covariances between them that can provide information on cosmology. Combining multiple probes allows for more accurate (less contaminated by systematics) and more precise (since there is cosmological information encoded in cross-correlation statistics) measurements of dark energy. The potential of cross-correlation methods is only beginning to be realized. By bringing in information from other wavelengths, the capabilities of the existing probes of dark energy can be enhanced and systematic effects can be mitigated further. We present a mixture of work in progress and suggestions for future scientific efforts. Given the scope of future dark energy experiments, the greatest gains may only be realized with more coordination and cooperation between multiple project teams; we recommend that this interchange should begin sooner, rather than later, to maximize scientific gains.

Dynamic Complexity is a phenomenon exhibited by a nonlinearly interacting system within which multitudes of different sizes of large scale coherent structures emerge, resulting in a globally nonlinear stochastic behavior vastly different from that could be surmised from the underlying equations of interaction. The hallmark of such nonlinear, complex phenomena is the appearance of intermittent fluctuating events with the mixing and distributions of correlated structures at all scales. We briefly review here a relatively recent method, ROMA (rank-ordered multifractal analysis), explicitly constructed to analyze the intricate details of the distribution and scaling of such types of intermittent structures. This method is then applied to the analyses of selected examples related to the dynamical plasmas of the cusp region of the magnetosphere, velocity fluctuations of classical hydrodynamic turbulence, and the distribution of the structures of the cosmic gas obtained through large scale, moving mesh simulations. Differences and similarities of the analyzed results among these complex systems will be contrasted and highlighted. The first two examples have direct relevance to the geospace environment and are summaries of previously reported findings. The third example on the cosmic gas, though involving phenomena much larger in spatiotemporal scales, with its highly compressible turbulent behavior and the unique simulation technique employed in generating the data, provides direct motivations of applying such analysis to studies of similar multifractal processes in various extreme environments. These new results are both exciting and intriguing.

The large-scale distribution of neutral hydrogen in the Universe will be luminous through its 21 cm emission. Here, for the first time, we use the auto-power spectrum of 21 cm intensity fluctuations to constrain neutral hydrogen fluctuations at z~0.8. Our data were acquired with the Green Bank Telescope and span the redshift range 0.6 < z < 1 over two fields totalling ~41 deg. sq. and 190 h of radio integration time. The dominant synchrotron foregrounds exceed the signal by ~10^3, but have fewer degrees of freedom and can be removed efficiently. Even in the presence of residual foregrounds, the auto-power can still be interpreted as an upper bound on the 21 cm signal. Our previous measurements of the cross-correlation of 21 cm intensity and the WiggleZ galaxy survey provide a lower bound. Through a Bayesian treatment of signal and foregrounds, we can combine both fields in auto- and cross-power into a measurement of Omega_HI b_HI = [0.62^{+0.23}_{-0.15}] * 10^{-3} at 68% confidence with 9% systematic calibration uncertainty, where Omega_HI is the neutral hydrogen (HI) fraction and b_HI is the HI bias parameter. We describe observational challenges with the present data set and plans to overcome them.

The GMRT Epoch of Reionization (EoR) experiment is an ongoing effort to measure the power spectrum from neutral hydrogen at high redshift. We have previously reported an upper limit of (70 mK)^2 at wavenumbers of k=0.65 h/Mpc using a basic piecewise-linear foreground subtraction. In this paper we explore the use of a singular value decomposition to remove foregrounds with fewer assumptions about the foreground structure. Using this method we also quantify, for the first time, the signal loss due to the foreground filter and present new power spectra adjusted for this loss, providing a revised measurement of a 2-sigma upper limit at (248 mK)^2 for k=0.50 h/Mpc. While this revised limit is larger than previously reported, we believe it to be more robust and still represents the best current constraints on reionization at z=8.6.

We present a new procedure to measure the large-scale carbon monoxide (CO) emissions across cosmic history. As a tracer of large-scale structure (LSS), the CO gas content as a function of redshift can be quantified by its three-dimensional fluctuation power spectra. Furthermore, cross-correlating CO emission with other LSS tracers offers a way to measure the emission as a function of scale and redshift. Here we introduce the model relevant for such a cross-correlation measurement between CO and other LSS tracers, and between different CO rotational lines. We propose a novel use of Cosmic Microwave Background (CMB) data and attempt to extract redshifted CO emissions embedded in the Wilkinson Microwave Anisotropy Probe (WMAP) dataset. We cross-correlate the all-sky WMAP7 data with LSS data sets, namely, the photometric quasar sample and the luminous red galaxy sample from the Sloan Digital Sky Survey Data Release 6 and 7 respectively. We are unable to detect a cross-correlation signal with either CO(1-0) nor CO(2-1) lines, mainly due to the instrumental noise in the WMAP data. However, we are able to rule out models more than three times greater than our more optimistic model. We discuss the cross-correlation signal from the thermal Sunyaev-Zeldovich effect and dust as potential contaminants, and quantify their impact for our CO measurements. We discuss forecasts for current CMB experiments and a hypothetical future CO focused experiment, and propose to cross-correlate CO temperature data with the Hobby-Eberly Telescope Dark Energy Experiment Ly-alpha-emitter sample, for which a signal-to-noise ratio of 58 is possible.

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

In this letter, 21 cm intensity maps acquired at the Green Bank Telescope are cross-correlated with large-scale structure traced by galaxies in the WiggleZ Dark Energy Survey. The data span the redshift range 0.6 < z < 1 over two fields totaling ~41 deg. sq. and 190 hours of radio integration time. The cross-correlation constrains Omega_HI b_HI r = [0.43 \pm 0.07 (stat.) \pm 0.04(sys.)] x 10^-3, where Omega_HI is the neutral hydrogen HI fraction, r is the galaxy-hydrogen correlation coefficient, and b_HI is the HI bias parameter. This is the most precise constraint on neutral hydrogen density fluctuations in a challenging redshift range. Our measurement improves the previous 21 cm cross-correlation at z ~ 0.8 both in its precision and in the range of scales probed.

The upcoming new generation of spectroscopic galaxy redshift surveys will provide large samples of cosmic voids, the distinct, large underdense structures in the universe. Combining these with future galaxy imaging surveys, we study the prospects of probing the underlying matter distribution in and around cosmic voids via the weak gravitational lensing effects of stacked voids, utilizing both shear and magnification information. The statistical precision is greatly improved by stacking together a large number of voids along different lines of sight, even when taking into account the impact of inherent miscentering and projection effects. We show that Dark Energy Task Force Stage IV surveys, such as the Euclid satellite and the Large Synoptic Survey Telescope, should be able to detect the void lensing signal with sufficient precision from stacking abundant medium-sized voids, thus providing direct constraints on the matter density profile of voids independent of assumptions on galaxy bias.

We quantify the prospects for using emission lines from rotational transitions of the CO molecule to perform an `intensity mapping' observation at high redshift during the Epoch of Reionization (EoR). The aim of CO intensity mapping is to observe the combined CO emission from many unresolved galaxies, to measure the spatial fluctuations in this emission, and use this as a tracer of large scale structure at very early times in the history of our Universe. This measurement would help determine the properties of molecular clouds -- the sites of star formation -- in the very galaxies that reionize the Universe. We further consider the possibility of cross-correlating CO intensity maps with future observations of the redshifted 21 cm line. The cross spectrum is less sensitive to foreground contamination than the auto power spectra, and can therefore help confirm the high redshift origin of each signal. Furthermore, the cross spectrum measurement would help extract key information about the EoR, especially regarding the size distribution of ionized regions. We discuss uncertainties in predicting the CO signal at high redshift, and discuss strategies for improving these predictions. Under favorable assumptions, and feasible specifications for a CO survey mapping the CO(2-1) and CO(1-0) lines, the power spectrum of CO emission fluctuations and its cross spectrum with future 21 cm measurements from the MWA are detectable at high significance.

We present a new upper limit to the 21cm power spectrum during the Epoch of Reionization (EoR) which constrains reionization models with an unheated IGM. The GMRT-EoR experiment is an ongoing effort to make a statistical detection of the power spectrum of 21cm neutral hydrogen emission at redshift z~9. Data from this redshift constrain models of the EoR, the end of the Dark Ages arising from the formation of the first bright UV sources, probably stars or mini-quasars. We present results from approximately 50 hours of observations at the Giant Metrewave Radio Telescope in India from December 2007. We describe radio frequency interference (RFI) localisation schemes which allow bright sources on the ground to be identified and physically removed in addition to automated flagging. Singular-value decomposition is used to remove remaining broadband RFI by identifying ground sources with large eigenvalues. Foregrounds are modelled using a piecewise linear filter and the power spectrum is measured using cross-correlations of foreground subtracted images.

Observations of 21-cm radio emission by neutral hydrogen at redshifts z ~ 0.5 to ~ 2.5 are expected to provide a sensitive probe of cosmic dark energy. This is particularly true around the onset of acceleration at z ~ 1, where traditional optical cosmology becomes very difficult because of the infrared opacity of the atmosphere. Hitherto, 21-cm emission has been detected only to z=0.24. More distant galaxies generally are too faint for individual detections but it is possible to measure the aggregate emission from many unresolved galaxies in the 'cosmic web'. Here we report a three dimensional 21-cm intensity field at z=0.53 to 1.12. We then co-add HI emission from the volumes surrounding about ten thousand galaxies (from the DEEP2 optical galaxy redshift survey. We detect the aggregate 21-cm glow at a significance of ~ 4 sigma.

Before the universe was reionized, it was likely that the spin temperature of intergalactic hydrogen was decoupled from the CMB by UV radiation from the first stars through the Wouthuysen-Field effect. If the IGM had not yet been heated above the CMB temperature by that time, then the gas would appear in absorption relative to the CMB. Large, rare sources of X-rays could inject sufficient heat into the neutral IGM, so that the differential brightness temperature was greater than zero at comoving distances of tens to hundreds of Mpc, resulting in large 21-cm fluctuations with amplitudes of about 250 mK on arcminute to degree angular scales, an order of magnitude larger in amplitude than that caused by ionized bubbles during reionization, about 25 mK. This signal could therefore be easier to detect and probe higher redshifts than that due to patchy reionization. For the case in which the first objects to heat the IGM are QSOs hosting 10^7-solar mass black holes with an abundance exceeding about 1 per Gpc^3 at z~15, observations with either the Arecibo Observatory or the Five Hundred Meter Aperture Spherical Telescope (FAST) could detect and image their fluctuations at greater than 5-sigma significance in about a month of dedicated survey time. Additionally, existing facilities such as MWA and LOFAR could detect the statistical fluctuations arising from a population of 10^5-solar mass black holes with an abundance of about 10^4 per Gpc^3 at z~10-12.

Rank-Ordered Multifractal Analysis (ROMA), a recently developed technique that combines the ideas of parametric rank ordering and one parameter scaling of monofractals, has the capabilities of deciphering the multifractal characteristics of intermittent fluctuations. The method allows one to understand the multifractal properties through rank-ordered scaling or non-scaling parametric variables. The idea of the ROMA technique is applied to analyze the multifractal characteristics of the auroral zone electric field fluctuations observed by SIERRA. The observed fluctuations span across contiguous multiple regimes of scales with different multifractal characteristics. We extend the ROMA technique such that it can take into account the crossover behavior -- with the possibility of collapsing probability distributions functions (PDFs) -- over these contiguous regimes.

The GMRT reionization effort aims to map out the large scale structure of the Universe during the epoch of reionization (EoR). Removal of polarized Galactic emission is a difficult part of any 21 cm EoR program, and we present new upper limits to diffuse polarized foregrounds at 150 MHz. We find no high significance evidence of polarized emission in our observed field at mid galactic latitude (J2000 08h26m+26). We find an upper limit on the 2-dimensional angular power spectrum of diffuse polarized foregrounds of [l^2 C_l/(2 PI)]^{1/2}< 3K in frequency bins of width 1 MHz at 300<l<1000. The 3-dimensional power spectrum of polarized emission, which is most directly relevant to EoR observations, is [k^3 P_p(k)/(2 PI^2)]^{1/2}< 2K at k_perp > 0.03 h/Mpc, k < 0.1 h/Mpc. This can be compared to the expected EoR signal in total intensity of [k^3 P(k)/ (2 PI^2) ]^{1/2} ~ 10 mK. We find polarized structure is substantially weaker than suggested by extrapolation from higher frequency observations, so the new low upper limits reported here reduce the anticipated impact of these foregrounds on EoR experiments. We discuss Faraday beam and depth depolarization models and compare predictions of these models to our data. We report on a new technique for polarization calibration using pulsars, as well as a new technique to remove broadband radio frequency interference. Our data indicate that, on the edges of the main beam at GMRT, polarization squint creates ~ 3% leakage of unpolarized power into polarized maps at zero rotation measure. Ionospheric rotation was largely stable during these solar minimum night time observations.

Using the 21 cm line, observed all-sky and across the redshift range from 0 to 5, the large scale structure of the Universe can be mapped in three dimensions. This can be accomplished by studying specific intensity with resolution ~ 10 Mpc, rather than via the usual galaxy redshift survey. The data set can be analyzed to determine Baryon Acoustic Oscillation wavelengths, in order to address the question: 'What is the nature of Dark Energy?' In addition, the study of Large Scale Structure across this range addresses the questions: 'How does Gravity effect very large objects?' and 'What is the composition our Universe?' The same data set can be used to search for and catalog time variable and transient radio sources.

We present an overview for the reionization search at GMRT. The forecast sensitivities are promising for an early detection. RFI mitigation has been successful. Several hundred hours of telescope time have already been invested in this ongoing effort, and analysis of the data is in progress.

We present the first statistically significant detection of cosmic structure using broadly distributed hydrogen radio emission. This is accomplished using a cross correlation with optical galaxies. Statistical noise levels of $20 \mu $K are achieved, unprecedented in this frequency band. This lends support to the idea that large volumes of the universe can be rapidly mapped without the need to resolve individual faint galaxies, enabling precise constraints to dark energy models. We discuss strategies for improved intensity mapping.

The expansion of the universe appears to be accelerating, and the mysterious anti-gravity agent of this acceleration has been called ``dark energy''. To measure the dynamics of dark energy, Baryon Acoustic Oscillations (BAO) can be used. Previous discussions of the BAO dark energy test have focused on direct measurements of redshifts of as many as $10^9$ individual galaxies, by observing the 21cm line or by detecting optical emission. Here we show how the study of acoustic oscillation in the 21 cm brightness can be accomplished by economical three dimensional intensity mapping. If our estimates gain acceptance they may be the starting point for a new class of dark energy experiments dedicated to large angular scale mapping of the radio sky, shedding light on dark energy.

We describe a new method that is both physically explicable and quantitatively accurate in describing the multifractal characteristics of intermittent events based on groupings of rank-ordered fluctuations. The generic nature of such rank-ordered spectrum leads it to a natural connection with the concept of one-parameter scaling for monofractals. We demonstrate this technique using results obtained from a 2D MHD simulation. The calculated spectrum suggests a crossover from the near Gaussian characteristics of small amplitude fluctuations to the extreme intermittent state of large rare events.

The ideas of dynamical complexity induced intermittent turbulence by sporadic localized interactions of coherent structures are discussed. In particular, we address the phenomenon of magnetic reconfiguration due to coarse-grained dissipation as well as the interwoven connection between criticality and multifractal processes. Specific examples are provided.

The recently introduced models of reionization bubbles based on extended Press-Schechter theory (Furlanetto, Zaldarriaga & Hernquist 2004) are generalized to include mergers of ionization sources. Sources with a recent major merger are taken to have enhanced photon production due to star formation, and accretion onto a central black hole if a black hole is present. This produces a scatter in the number of ionized photons corresponding to a halo of a given mass and a change in photon production over time for any given halo mass. Photon production histories, bubble distributions, and ionization histories are computed for several different parameter and recombination assumptions; the resulting distributions interpolate between previously calculated limiting cases.

We present the first measurement of weak lensing by large-scale structure on 1-4 degree scales based on the FIRST Radio Survey, a quarter-sky, 20-cm survey produced with the NRAO Very Large Array (VLA). The large angular scales provide an unique measurement in the linear regime of the matter power spectrum. We use the shapelet method to estimate the shear of radio sources derived directly from the Fourier(uv) plane. With realistic simulations we verify that the method yields unbiased shear estimators. We quantify and correct for the systematic effects which can produce spurious shears, analytically and with simulations. We measure the shear correlation functions on angular scales of 0.5-40 degrees, and compute the corresponding aperture mass statistics. On 1-4 degrees scales, we find that the B-modes are consistent with zero, and detect a lensing E-mode signal significant at the 3.0 sigma level. After removing radio sources with an optical counterpart, the E-mode signal increases by 10-20%, as expected for a lensing signal derived from more distant sources. We use the E-mode measurement to constrain the mass power spectrum normalization sigma_8 and the median redshift z_m of the unidentified radio sources. We find sigma_8*(z_m/2)^(0.6) ~ 1.0 +/- 0.2, where the 1 sigma error bars include statistical errors, cosmic variance, and systematics. This is consistent with earlier determinations of sigma_8, and with our current knowledge of the redshift distribution of radio sources. Taking the prior sigma_8 = 0.9 +/- 0.1 (68%CL) from the WMAP experiment, this corresponds to z_m = 2.2 +/- 0.9 (68%CL) for radio sources without optical counterparts. Our results offer promising prospects for precision measurements of cosmic shear with future radio interferometers such as LOFAR and the SKA.

We present a new approach for image reconstruction and weak lensing measurements with interferometers. Based on the shapelet formalism presented in Refregier (2001), object images are decomposed into orthonormal Hermite basis functions. The shapelet coefficients of a collection of sources are simultaneously fit on the uv plane, the Fourier transform of the sky brightness distribution observed by interferometers. The resulting chi-square fit is linear in its parameters and can thus be performed efficiently by simple matrix multiplications. We show how the complex effects of bandwidth smearing, time averaging and non-coplanarity of the array can be easily and fully corrected for in our method. Optimal image reconstruction, co-addition, astrometry, and photometry can all be achieved using weighted sums of the derived coefficients. As an example we consider the observing conditions of the FIRST radio survey (Becker et al. 1995; White et al. 1997). We find that our method accurately recovers the shapes of simulated images even for the sparse uv sampling of this snapshot survey. Using one of the FIRST pointings, we find our method compares well with CLEAN, the commonly used method for interferometric imaging. Our method has the advantage of being linear in the fit parameters, of fitting all sources simultaneously, and of providing the full covariance matrix of the coefficients, which allows us to quantify the errors and cross-talk in image shapes. It is therefore well-suited for quantitative shape measurements which require high-precision. In particular, we show how our method can be combined with the results of Refregier & Bacon (2001) to provide an accurate measurement of weak lensing from interferometric data.

We present a new approach to measure the shapes of galaxies, a fundamental task in observational astronomy. This approach is based on the decomposition of a galaxy image into a series of orthogonal basis functions, or `shapelets'. Our choice of basis functions, namely the Gauss-Hermite series, has a number of remarkable properties under distortions, convolutions and noise, which makes them particularly well suited for astrophysical applications. In particular, we describe how they can be used to measure the shear induced by weak gravitational lensing, with the precision required for upcoming surveys. We also show how shapelets can be used to reconstruct images from interferometric observations. Other application of shapelets, such as image compression, PSF deconvolution, de-projection and the study of galaxy morphology, are also briefly discussed.

We present the results of HI line and radio continuum observations of five nearby E+A galaxies. These galaxies have spectra that are dominated by a young stellar component but lack the emission lines characteristic of significant, on-going star formation. They are selected from a unique sample of 21 E+A's identified by Zabludoff et al.(1996) in their spectroscopic search for E+A galaxies using the Las Campanas Redshift Survey, where over 11,000 nearby galaxies were examined. The five E+A galaxies span a range of environments: three are in the field and two are in clusters. Only one system was detected in HI emission, the field E+A galaxy EA1, with a total flux of 0.30 +/- 0.02 Jy km/s and an HI mass of (3.5 +/- 0.2) 10^9 h^(-2) M_sun. The HI morphology and kinematics of EA 1 suggest a galaxy-galaxy interaction, with a dynamical age of about 6 x 10^8 h^(-1) yr inferred from the HI tail lengths and velocities. This age estimate is consistent with the interpretation drawn from optical spectroscopy that starbursts in E+A galaxies began (and subsequently ended) within the last 10^9 yr. Our HI detection limits are such that if the other E+A's in our sample had the HI properties of EA 1, we would have detected (or marginally detected) them. We conclude that E+A galaxies have a range of HI properties. None of the galaxies were detected in radio continuum emission, with upper limits to the radio power of about 10^21 h^(-2) W/Hz. Our limits exclude the possibility that these E+A's are dust-enshrouded massive starburst galaxies, but are insufficient to exclude modest star formation rates of less than a few h^(-2) M_sun per yr.

The large H II complex N11 in the Large Magellanic Cloud contains OB associations at several different stages in their life histories. We have obtained ROSAT PSPC and HRI X-ray observations, Curtis Schmidt CCD images, echelle spectra in H-alpha and [N II] lines, and IUE interstellar absorption line observations of this region. The central bubble of N11 has an X-ray luminosity a factor of only 3-7 brighter than predicted for an energy-conserving superbubble, making this the first detection of X-ray emission from a superbubble without a strong X-ray excess. The region N11B contains an extremely young OB association analogous to the central association of the Carina nebula, apparently still embedded in its natal molecular cloud. We find that N11B emits diffuse X-ray emission, probably powered by stellar winds. Finally, we compare the tight cluster HD32228 in N11 to R136 in 30 Dor. The latter is a strong X-ray source, while the former is not detected, showing that strong X-ray emission from compact objects is not a universal property of such tight clusters.