Search

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.

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data., Comment: 39 pages, 27 figures. Submitted to ApJ. Comments welcome

The observed large-scale scatter in Lyman $\alpha$ opacity of the intergalactic medium at $z, Comment: Accepted to ApJ. 17 pages, 13 figures

We present deep Keck/MOSFIRE $H$-band spectroscopic observations covering the [CIII],CIII]$\lambda\lambda1907,1909$ doublet for three $z\sim8$ galaxy candidates in the AEGIS field. Along with non-detections in two galaxies, we obtain one of the highest-redshift detections to-date of [CIII]$\lambda 1907$ for the galaxy AEGIS-33376, finding $z_{\rm spec}=7.945\pm0.001$. We measure a [CIII]$\lambda$1907 flux of $2.24\pm0.71\times10^{-18} \mbox{ erg}\mbox{ s}^{-1} \mbox{ cm}^{-2}$, corresponding to a rest-frame equivalent width of $20.3\pm6.5 \unicode{x212B}$ for the single line. Given the not very constraining upper limit for CIII]$\lambda 1909$ based on strong sky-line contamination, we assume a [CIII]$\lambda$1907/CIII]$\lambda 1909$ doublet ratio of 1.5 and infer a total [CIII],CIII]$\lambda\lambda1907,1909$ equivalent width of $33.7\pm 10.8 \unicode{x212B}$. We repeat the same reductions and analysis on multiple subsets of our raw data divided on the basis of time and observing conditions, verifying that the [CIII]$\lambda 1907$ emission is present for AEGIS-33376 throughout our observations. We also confirm that the significance of the [CIII]$\lambda 1907$ detection in different subsets of our data tracks that of brighter emission features detected on the same multi-slit mask. These multiple tests suggest that the observed emission line is real and associated with the $z\sim 8$ target. The strong observed [CIII],CIII]$\lambda\lambda1907,1909$ in AEGIS-33376 likely indicates ISM conditions of low metallicity, high ionization parameter, and a hard ionizing spectrum, although AGN contributions are possible. This single detection represents a sizable increase in the current sample [CIII],CIII]$\lambda\lambda1907,1909$ detections at $z>7$, while $\textit{JWST}$ will provide the first statistical samples of such measurements at these redshifts., Comment: 6 pages, 3 figures, 1 table; Accepted to ApJL

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

We examine large-scale fluctuations in the He?II Ly? forest transmission during and after He?II reionization. We use a simple Monte Carlo model to distribute quasars throughout a large volume and compute the resulting radiation field along one-dimensional skewers. In agreement with previous studies, we find that the rarity of these sources induces order unity fluctuations in the mean optical depth after reionization, even when averaged over large segments (~10-100Mpc across). We compare our models to existing data along five He?II Ly? forest lines of sight spanning z ~ 2-3.2. The large cosmic variance contained in our model plausibly explains many of the observed fluctuations at z 2.5. But our models cannot accommodate the large fluctuations toward high optical depths on ~30Mpc scales observed at z ~ 2.7-2.9, and the measured optical depths (?eff 4) at z>2.9 are difficult to explain with a smoothly evolving mean radiation field. In order to better understand these data, we construct a toy model of He?II reionization, in which we assume that regions with the smallest radiation fields in a post-reionization universe (or farthest from strong ionizing sources) are completely dark during reionization. The observed fluctuations fit much more comfortably into this model, and we therefore argue that, according to present data, He?II reionization does not complete until z 2.9.

We use an analytic model to study how inhomogeneous hydrogen reionization affects the temperature distribution of the intergalactic medium (IGM). During this process, the residual energy of each ionizing photon is deposited in the IGM as heat, increasing its temperature to 20,000-30,000 K; subsequent expansion of the Universe then cools the gas. Because reionization most likely proceeds from high to low densities, underdense voids are ionized last, have less time to cool, and are (on average) warmer than mean-density gas immediately after reionization is complete (an "inverted" equation of state). From this initial configuration, the low-density gas cools quickly and eventually returns to a more normal equation of state. The rapidly evolving temperature introduces systematic uncertainties in measurements of the ionizing background at z~6. For example, late reionization implies rapid cooling, so that the ionizing background would have to evolve even more rapidly at z ~5-6 than typically claimed. This degeneracy is difficult to disentangle, because the Lyman-alpha forest probes only a narrow range in densities (over which the gas is nearly isothermal). However, higher Lyman-series transitions probe wider density ranges, sampling different effective temperatures, and offer a new way to measure the IGM equation of state that should work where nearly saturated absorption precludes other methods. This will help to separate evolution in temperature from that in the ionizing background. While more detailed study with hydrodynamic simulations is needed, we show that such measurements could potentially distinguish early and late reionization using only a handful of lines of sight.

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of co-moving side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, which evolves from a sub-Mpc scale at the beginning of reionization to R > 10 Mpc towards the end. We present a detailed statistical description of our results, and compare them with a numerical scheme based on the analytic model by Furlanetto, Zaldarriaga, and Hernquist. We find that the analytic calculation reproduces the size distribution of HII regions and the 21 cm power spectrum of the radiative transfer simulation remarkably well. The ionization field from the simulation, however, has more small scale structure than the analytic calculation, owing to Poisson scatter in the simulated abundance of galaxies on small scales. We propose and validate a simple scheme to incorporate this scatter into our calculations. Our results suggest that analytic calculations are sufficiently accurate to aid in predicting and interpreting the results of future 21 cm surveys. In particular, our fast numerical scheme is useful for forecasting constraints from future 21 cm surveys, and in constructing mock surveys to test data analysis procedures., 15 pages, 11 figures, submitted to ApJ

Before reionization, the intergalactic medium (IGM) may have been sufficiently cold for low-mass "minihalos" to condense out of the gas and subsequently affect reionization. Previous work has shown that minihalos generate reasonably large 21 cm fluctuations. Here we consider this signal in its proper cosmological context and show that isolating minihalos from the rest of the IGM is extremely difficult. Using the well-known halo model, we compute the power spectrum of 21 cm fluctuations from minihalos and show that the signal decreases rapidly as feedback increases the Jeans mass. We then show that even a small Lyman-alpha background increases the 21 cm fluctuations of the diffuse IGM well beyond those of the minihalos; because the mass fraction in the IGM is much larger, minihalos will lie buried within the IGM signal. The distinctive signatures of non-linear bias and minihalo structure emerge only at much smaller scales, well beyond the resolution of any upcoming instruments. Using simple, but representative, reionization histories, we then show that the required Lyman-alpha background level is most likely achieved at z>15, while minihalos are still rare, so that they are almost always degenerate with the diffuse IGM.

The presence of multiple luminous galaxies in clusters can be explained by the finite time over which a galaxy sinks to the center of the cluster and merges with the the central galaxy. The simplest measurable statistic to quantify the dynamical age of a system of galaxies is the luminosity (magnitude) gap, which is the difference in photometric magnitude between the two most luminous galaxies. We present a simple analytical estimate of the luminosity gap distribution in groups and clusters as a function of dark matter halo mass. The luminosity gap is used to define "fossil" groups; we expect the fraction of fossil systems to exhibit a strong and model-independent trend with mass: ~1-3% of massive clusters and ~5-40% of groups should be fossil systems. We also show that, on cluster scales, the observed intrinsic scatter in the central galaxy luminosity-halo mass relation can be ascribed to dispersion in the merger histories of satellites within the cluster. We compare our predictions to the luminosity gap distribution in a sample of 730 clusters in the Sloan Digital Sky Survey C4 Catalog and find good agreement. This suggests that theoretical excursion set merger probabilities and the standard theory of dynamical segregation are valid on cluster scales., Comment: 5 pages, 3 figures

During the epoch of reionization, local variations in the ionized fraction (patchiness) imprint arcminute-scale temperature anisotropies in the CMB through the kinetic Sunyaev-Zel'dovich (kSZ) effect. We employ an improved version of an analytic model of reionization devised in Furlanetto et al.(2004) to calculate the kSZ anisotropy from patchy reionization. This model uses extended Press-Schechter theory to determine the distribution and evolution of HII bubbles and produces qualitatively similar reionization histories to those seen in recent numerical simulations. We find that the angular power spectrum of the kSZ anisotropies depends strongly on the size distribution of the HII bubbles and on the duration of reionization. In addition, we show that upcoming measurements of the kSZ effect should be able to distinguish between several popular reionization scenarios. In particular, the amplitude of the patchy power spectrum for reionization scenarios in which the IGM is significantly ionized by Population III stars (or by mini-quasars/decaying particles) can be larger (or smaller) by over a factor of 3 than the amplitude in more traditional reionization histories (with temperature anisotropies that range between 0.5 and 3 micro-Kelvin at l = 5000). We highlight the differences in the kSZ signal between many possible reionization morphologies and discuss the constraints that future observations of the kSZ will place on this epoch., Comment: 14 pages, 10 figures. Accepted for publication in ApJ

The nature of the interaction between galaxies and the intergalactic medium (IGM) is one of the most fundamental problems in astrophysics. The accretion of gas onto galaxies provides fuel for star formation, while galactic winds transform the nearby IGM in a number of ways. One exciting technique to study this gas is through the imaging of hydrogen Lyman-alpha emission. We use cosmological simulations to study the Lyman-alpha signals expected from the growth of cosmic structure from z=0-5. We show that if dust absorption is negligible, recombinations following the absorption of stellar ionizing photons dominate the total Lyman-alpha photon production rate. However, galaxies are also surrounded by "Lyman-alpha coronae" of diffuse IGM gas. These coronae are composed of a combination of accreting gas and material ejected from the central galaxy by winds. The Lyman-alpha emission from this phase is powered by a combination of gravitational processes and the photoionizing background. While the former dominates at z~0, collisional excitation following photo-heating may well dominate the total emission at higher redshifts. The central regions of these systems are dense enough to shield themselves from the metagalactic ionizing background; unfortunately, in this regime our simulations are no longer reliable. We therefore consider several scenarios for the emission from the central cores, including one in which self-shielded gas does not emit at all. We show that the combination of star formation and cooling IGM gas can explain most of the observed "Lyman-alpha blobs" at z~3, with the important exception of the largest sources. On the other hand, except under the most optimistic assumptions, cooling IGM gas cannot explain the observations on its own.

We consider the degree to which "21 cm tomography" of the high-redshift Universe can distinguish different ionization histories. Using a new analytic model for the size distribution of HII regions that associates these ionized bubbles with large-scale galaxy overdensities, we compute the angular power spectrum and other statistical properties of the 21 cm brightness temperature during reionization. We show that the HII regions imprint features on the power spectrum that allow us to separate histories with discrete bubbles from those with partial uniform ionization (by, for example, X-rays). We also show that "double reionization" scenarios will modify the morphology of the bubbles in ways that depend on the mechanism through which the first generation of sources shuts off. If, for example, the transition occurs globally at a fixed redshift, the first generation imprints a persistent feature on the 21 cm power spectrum. Finally, we compare our model to one in which voids are ionized first. While the power spectra of these two models are qualitatively similar, we show that the underlying distributions of neutral hydrogen differ dramatically and suggest that other statistical tests can distinguish them. The next generation of low-frequency radio telescopes will have the sensitivity to distinguish all of these models and strongly constrain the history and morphology of reionization.

Clusters of galaxies are enriched with positrons from jets of active galactic nuclei (AGNs) or from the interaction of cosmic rays with the intracluster gas. We follow the cooling of these positrons and show that their eventual annihilation with cluster electrons yields a narrow annihilation line. Unlike annihilation in the interstellar medium of galaxies, the line produced in clusters is not smeared by three-photon decay of positronium, because positronium formation is suppressed at the high (>~ 1 keV) temperature of the cluster electrons. We show that if AGN jets are composed of e^+e^- pairs, then the annihilation line from rich clusters within a distance of 100 Mpc might be detectable with future space missions, such as INTEGRAL or EXIST., 39 pages, 11 figures, submitted to ApJ

We show that comparisons of HeII Lyman-alpha forest lines of sight to nearby quasar populations can strongly constrain the lifetimes and emission geometry of quasars. By comparing the HeII and HI Lyman-alpha forests along a particular line of sight, one can trace fluctuations in the hardness of the radiation field (which are driven by fluctuations in the HeII ionization rate). Because this high-energy background is highly variable - thanks to the rarity of the bright quasars that dominate it and the relatively short attenuation lengths of these photons - it is straightforward to associate features in the radiation field with their source quasars. Here we quantify how finite lifetimes and beamed emission geometries affect these expectations. Finite lifetimes induce a time delay that displaces the observed radiation peak relative to the quasar. For beamed emission, geometry dictates that sources invisible to the observer can still create a peak in the radiation field. We show that both these models produce substantial populations of "bare" peaks (without an associated quasar) for reasonable parameter values (lifetimes ~10^6-10^8 yr and beaming angles