The redshifted 21-cm background is expected to be a powerful probe of the early Universe, carrying both cosmological and astrophysical information from a wide range of redshifts. In particular, the power spectrum of fluctuations in the 21-cm brightne
ss temperature is anisotropic due to the line-of-sight velocity gradient, which in principle allows for a simple extraction of this information in the limit of linear fluctuations. However, recent numerical studies suggest that the 21-cm signal is actually rather complex, and its analysis likely depends on detailed model fitting. We present the first realistic simulation of the anisotropic 21-cm power spectrum over a wide period of early cosmic history. We show that on observable scales, the anisotropy is large and thus measurable at most redshifts, and its form tracks the evolution of 21-cm fluctuations as they are produced early on by Lyman-a radiation from stars, then switch to X-ray radiation from early heating sources, and finally to ionizing radiation from stars. In particular, we predict a redshift window during cosmic heating (at z ~ 15), when the anisotropy is small, during which the shape of the 21-cm power spectrum on large scales is determined directly by the average radial distribution of the flux from X-ray sources. This makes possible a model-independent reconstruction of the X-ray spectrum of the earliest sources of cosmic heating.
Models and simulations of the epoch of reionization predict that spectra of the 21-cm transition of atomic hydrogen will show a clear fluctuation peak, at a redshift and scale, respectively, that mark the central stage of reionization and the charact
eristic size of ionized bubbles. This is based on the assumption that the cosmic gas was heated by stellar remnants - particularly X-ray binaries - to temperatures well above the cosmic microwave background at that time (~ 30 K). Here we show instead that the hard spectra (that is, spectra with more high-energy photons than low-energy photons) of X-ray binaries make such heating ineffective, resulting in a delayed and spatially uniform heating that modifies the 21-cm signature of reionization. Rather than looking for a simple rise and fall of the large-scale fluctuations (peaking at several millikelvin), we must expect a more complex signal also featuring a distinct minimum (at less than a millikelvin) that marks the rise of the cosmic mean gas temperature above the microwave background. Observing this signal, possibly with radio telescopes in operation today, will demonstrate the presence of a cosmic background of hard X-rays at that early time.
Understanding the formation and evolution of the first stars and galaxies represents one of the most exciting frontiers in astronomy. Since the universe was filled with neutral hydrogen at early times, the most promising method for observing the epoc
h of the first stars is using the prominent 21-cm spectral line of the hydrogen atom. Current observational efforts are focused on the reionization era (cosmic age t~500 Myr), with earlier times considered much more challenging. However, the next frontier of even earlier galaxy formation (t~200 Myr) is emerging as a promising observational target. This is made possible by a recently noticed effect of a significant relative velocity between the baryons and dark matter at early times. The velocity difference significantly suppresses star formation. The spatial variation of this suppression enhances large-scale clustering and produces a prominent cosmic web on 100 comoving Mpc scales in the 21-cm intensity distribution. This structure makes it much more feasible for radio astronomers to detect these early stars, and should drive a new focus on this era, which is rich with little-explored astrophysics.
We present the first complete calculation of the history of the inhomogeneous 21-cm signal from neutral hydrogen during the era of the first stars. We use hybrid computational methods to capture the large-scale distribution of the first stars, whose
radiation couples to the neutral hydrogen emission, and to evaluate the 21-cm signal from z ~ 15-35. In our realistic picture large-scale fluctuations in the 21-cm signal are sourced by the inhomogeneous density field and by the Ly-alpha and X-ray radiative backgrounds. The star formation is suppressed by two spatially varying effects: negative feedback provided by the Lyman-Werner radiative background, and supersonic relative velocities between the gas and dark matter. Our conclusions are quite promising: we find that the fluctuations imprinted by the inhomogeneous Ly-alpha background in the 21-cm signal at z ~ 25 should be detectable with the Square Kilometer Array.
Understanding the formation and evolution of the first stars and galaxies represents one of the most exciting frontiers in astronomy. Since the universe was filled with neutral hydrogen at early times, the most promising method for observing the epoc
h of the first stars is using the prominent 21-cm spectral line of the hydrogen atom. Current observational efforts are focused on the reionization era (cosmic age t ~ 500 Myr), with earlier times considered much more challenging. However, the next frontier of even earlier galaxy formation (t ~ 200 Myr) is emerging as a promising observational target. This is made possible by a recently noticed effect of a significant relative velocity between the baryons and dark matter at early times. The velocity difference suppresses star formation, causing a unique form of early luminosity bias. The spatial variation of this suppression enhances large-scale clustering and produces a prominent cosmic web on 100 comoving Mpc scales in the 21-cm intensity distribution. This structure makes it much more feasible for radio astronomers to detect these early stars, and should drive a new focus on this era, which is rich with little-explored astrophysics.
Dark and baryonic matter moved at different velocities in the early Universe, which strongly suppressed star formation in some regions. This was estimated to imprint a large-scale fluctuation signal of about 2 mK in the 21-cm spectral line of atomic
hydrogen associated with stars at a redshift of 20, although this estimate ignored the critical contribution of gas heating due to X-rays and major enhancements of the suppression. A large velocity difference reduces the abundance of halos and requires the first stars to form in halos of about a million solar masses, substantially greater than previously expected. Here we report a simulation of the distribution of the first stars at z=20 (cosmic age of ~180 Myr), incorporating all these ingredients within a 400 Mpc box. We find that the 21-cm signature of these stars is an enhanced (10 mK) fluctuation signal on the 100-Mpc scale, characterized by a flat power spectrum with prominent baryon acoustic oscillations. The required sensitivity to see this signal is achievable with an integration time of a thousand hours with an instrument like the Murchison Wide-field Array or the Low Frequency Array but designed to operate in the range of 50-100 MHz.
Baryonic acoustic oscillations (BAOs) modulate the density ratio of baryons to dark matter across large regions of the Universe. We show that the associated variation in the mass-to-light ratio of galaxies should generate an oscillatory, scale-depend
ent bias of galaxies relative to the underlying distribution of dark matter. A measurement of this effect would calibrate the dependence of the characteristic mass-to-light ratio of galaxies on the baryon mass fraction in their large scale environment. This bias, though, is unlikely to significantly affect measurements of BAO peak positions.
Gravitational lensing observations of massive X-ray clusters imply a steep characteristic density profile marked by a central concentration of dark matter. The observed mass fraction within a projected radius of 150 kpc is twice that found in state-o
f-the-art dark matter simulations of the standard Lambda-CDM cosmology. A central baryon enhancement that could explain this discrepancy is not observed, leaving a major puzzle. We propose a solution based on the merger histories of clusters. A significant fraction of the final dark matter content of a cluster halo originates within galaxy-sized halos, in which gas can cool and compress the dark matter core to high densities. The subsequent tidal stripping of this compressed dark matter occurs in denser regions that are closer to the center of the cluster halo. Eventually, the originally cooled gas must be dispersed into the intracluster medium through feedback, for consistency with observations that do not find central baryon enhancements in clusters. Still, the early adiabatic compression of the galactic dark matter leaves a net effect on the cluster. Using a simple model for this process, we show that the central cluster profile is substantially modified, potentially explaining the observed discrepancy.
A new generation of radio telescopes are currently being built with the goal of tracing the cosmic distribution of atomic hydrogen at redshifts 6-15 through its 21-cm line. The observations will probe the large-scale brightness fluctuations sourced b
y ionization fluctuations during cosmic reionization. Since detailed maps will be difficult to extract due to noise and foreground emission, efforts have focused on a statistical detection of the 21-cm fluctuations. During cosmic reionization, these fluctuations are highly non-Gaussian and thus more information can be extracted than just the one-dimensional function that is usually considered, i.e., the correlation function. We calculate a two-dimensional function that if measured observationally would allow a more thorough investigation of the properties of the underlying ionizing sources. This function is the probability distribution function (PDF) of the difference in the 21-cm brightness temperature between two points, as a function of the separation between the points. While the standard correlation function is determined by a complicated mixture of contributions from density and ionization fluctuations, we show that the difference PDF holds the key to separately measuring the statistical properties of the ionized regions.
Analytical approaches to galaxy formation and reionization are based on the mathematical problem of random walks with barriers. The statistics of a single random walk can be used to calculate one-point distributions ranging from the mass function of
virialized halos to the distribution of ionized bubble sizes during reionization. However, an analytical calculation of two-point correlation functions or of spatially-dependent feedback processes requires the joint statistics of random walks at two different points. An accurate analytical expression for the statistics of two correlated random walks has been previously found only for the case of a constant barrier height. However, calculating bubble sizes or accurate statistics for halo formation involves more general barriers that can often be approximated as linear barriers. We generalize the two-point solution with constant barriers to linear barriers, and apply it as an illustration to calculate the correlation function of cosmological 21-cm fluctuations during reionization.
