Extreme Galaxies During Reionization: Testing ISM and Disk Models

We test the ability of equilibrium galactic disk and one-zone interstellar medium models to describe the physical and emission properties of quasar hosts, submillimeter galaxies, and Lyman-alpha emitters at z>~6. The size, line widths, star formation rates, black hole accretion rates, gas masses and temperatures, and the relationships between these properties are all well-described by our model, and we provide approximate fitting formulae for comparison with future observations. However, comparing our carbon line predictions to observations reveals differences between the ISM at low and high redshifts. Our underestimate of the [CII] line emission indicates either higher star formation efficiencies in high-redshift molecular clouds or less depletion of metals into dust at fixed metallicity. Further, our over-prediction of the CO(6-5)/CO(1-0) ratio suggests that molecular clouds in real high-redshift galaxies have a lower turbulent Mach number and more subthermal CO(6-5) emission than expected owing either to sizes smaller than the local Jeans mass or to a pressure support mechanism other than turbulence.Comment: Accepted in MNRAS; 19 pages; 10 figures; 4 table

Fast Large Volume Simulations of the 21 cm Signal from the Reionization and pre-Reionization Epochs

While limited to low spatial resolution, the next generation low-frequency radio interferometers that target 21 cm observations during the era of reionization and prior will have instantaneous fields-of-view that are many tens of square degrees on the sky. Predictions related to various statistical measurements of the 21 cm brightness temperature must then be pursued with numerical simulations of reionization with correspondingly large volume box sizes, of order 1000 Mpc on one side. We pursue a semi-numerical scheme to simulate the 21 cm signal during and prior to Reionization by extending a hybrid approach where simulations are performed by first laying down the linear dark matter density field, accounting for the non-linear evolution of the density field based on second-order linear perturbation theory as specified by the Zel'dovich approximation, and then specifying the location and mass of collapsed dark matter halos using the excursion-set formalism. The location of ionizing sources and the time evolving distribution of ionization field is also specified using an excursion-set algorithm. We account for the brightness temperature evolution through the coupling between spin and gas temperature due to collisions, radiative coupling in the presence of Lyman-alpha photons and heating of the intergalactic medium, such as due to a background of X-ray photons. The hybrid simulation method we present is capable of producing the required large volume simulations with adequate resolution in a reasonable time so a large number of realizations can be obtained with variations in assumptions related to astrophysics and background cosmology that govern the 21 cm signal.Comment: 14 pages and 15 figures. New version to match accepted version for MNRAS. Code available in: http://www.SimFast21.or

Lyman-alpha Damping Wing Constraints on Inhomogeneous Reionization

One well-known way to constrain the hydrogen neutral fraction, x_H, of the high-redshift intergalactic medium (IGM) is through the shape of the red damping wing of the Lya absorption line. We examine this method's effectiveness in light of recent models showing that the IGM neutral fraction is highly inhomogeneous on large scales during reionization. Using both analytic models and "semi-numeric" simulations, we show that the "picket-fence" absorption typical in reionization models introduces both scatter and a systematic bias to the measurement of x_H. In particular, we show that simple fits to the damping wing tend to overestimate the true neutral fraction in a partially ionized universe, with a fractional error of ~ 30% near the middle of reionization. This bias is generic to any inhomogeneous model. However, the bias is reduced and can even underestimate x_H if the observational sample only probes a subset of the entire halo population, such as quasars with large HII regions. We also find that the damping wing absorption profile is generally steeper than one would naively expect in a homogeneously ionized universe. The profile steepens and the sightline-to-sightline scatter increases as reionization progresses. Of course, the bias and scatter also depend on x_H and so can, at least in principle, be used to constrain it. Damping wing constraints must therefore be interpreted by comparison to theoretical models of inhomogeneous reionization.Comment: 11 pages, 10 figures; submitted to MNRA

A Flexible Analytic Model of Cosmic Variance in the First Billion Years

Cosmic variance is the intrinsic scatter in the number density of galaxies due to fluctuations in the large-scale dark matter density field. In this work, we present a simple analytic model of cosmic variance in the high redshift Universe ($z\sim5-15$). We assume that galaxies grow according to the evolution of the halo mass function, which we allow to vary with large-scale environment. Our model produces a reasonable match to the observed ultraviolet luminosity functions in this era by regulating star formation through stellar feedback and assuming that the UV luminosity function is dominated by recent star formation. We find that cosmic variance in the UVLF is dominated by the variance in the underlying dark matter halo population, and not by differences in halo accretion or the specifics of our stellar feedback model. We also find that cosmic variance dominates over Poisson noise for future high-$z$ surveys except for the brightest sources or at very high redshifts ($z \gtrsim 12$). We provide a linear approximation of cosmic variance for a variety of redshifts, magnitudes, and survey areas through the public Python package galcv. Finally, we introduce a new method for incorporating priors on cosmic variance into estimates of the galaxy luminosity function and demonstrate that it significantly improves constraints on that important observable

Reionization History from Coupled CMB/21cm Line Data

We study CMB secondary anisotropies produced by inhomogeneous reionization by means of cosmological simulations coupled with the radiative transfer code CRASH. The reionization history is consistent with the WMAP Thomson optical depth determination. We find that the signal arising from this process dominates over the primary CMB component for l > 4000 and reaches a maximum amplitude of l(l+1)C_l/2\pi ~ 1.6 x 10^{-13} on arcmin scale, i.e. l as large as several thousands. We then cross-correlate secondary CMB anisotropy maps with neutral hydrogen 21cm line emission fluctuations obtained from the same simulations. The two signals are highly anti-correlated on angular scales corresponding to the typical size of HII regions (including overlapping) at the 21cm map redshift. We show how the CMB/21cm cross-correlation can be used to: (a) study the nature of the reionization sources, (b) reconstruct the cosmic reionization history, (c) infer the mean cosmic ionization level at any redshift. We discuss the feasibility of the proposed experiment with forthcoming facilities.Comment: 7 pages, 5 figures, MNRAS in pres

Constraining the Collisional Nature of the Dark Matter Through Observations of Gravitational Wakes

We propose to use gravitational wakes as a direct observational probe of the collisional nature of the dark matter. We calculate analytically the structure of a wake generated by the motion of a galaxy in the core of an X-ray cluster for dark matter in the highly-collisional and collisionless limits. We show that the difference between these limits can be recovered from detailed X-ray or weak lensing observations. We also discuss the sizes of sub-halos in these limits. Preliminary X-ray data on the motion of NGC 1404 through the Fornax group disfavors fluid-like dark matter but does not exclude scenarios in which the dark matter is weakly collisional.Comment: 29 pages, 3 figures, submitted to Ap

The new galaxy evolution paradigm revealed by the Herschel surveys

The Herschel Space Observatory has revealed a very different galaxyscape from that shown by optical surveys which presents a challenge for galaxy-evolution models. The Herschel surveys reveal (1) that there was rapid galaxy evolution in the very recent past and (2) that galaxies lie on a single Galaxy Sequence (GS) rather than a star-forming ‘main sequence’ and a separate region of ‘passive’ or ‘red-and-dead’ galaxies. The form of the GS is now clearer because far-infrared surveys such as the Herschel ATLAS pick up a population of optically red starforming galaxies that would have been classified as passive using most optical criteria. The space-density of this population is at least as high as the traditional star-forming population. By stacking spectra of H-ATLAS galaxies over the redshift range 0.001 < z < 0.4, we show that the galaxies responsible for the rapid low-redshift evolution have high stellar masses, high star-formation rates but, even several billion years in the past, old stellar populations – they are thus likely to be relatively recent ancestors of early-type galaxies in the Universe today. The form of the GS is inconsistent with rapid quenching models and neither the analytic bathtub model nor the hydrodynamical EAGLE simulation can reproduce the rapid cosmic evolution. We propose a new gentler model of galaxy evolution that can explain the new Herschel results and other key properties of the galaxy population

Cosmological Parameter Estimation Using 21 cm Radiation from the Epoch of Reionization

A number of radio interferometers are currently being planned or constructed to observe 21 cm emission from reionization. Not only will such measurements provide a detailed view of that epoch, but, since the 21 cm emission also traces the distribution of matter in the Universe, this signal can be used to constrain cosmological parameters at 6 < z < 20. The sensitivity of an interferometer to the cosmological information in the signal may depend on how precisely the angular dependence of the 21 cm 3-D power spectrum can be measured. Utilizing an analytic model for reionization, we quantify all the effects that break the spherical symmetry of the 3-D 21 cm power spectrum and produce physically motivated predictions for this power spectrum. We find that upcoming observatories will be sensitive to the 21 cm signal over a wide range of scales, from larger than 100 to as small as 1 comoving Mpc. We consider three methods to measure cosmological parameters from the signal: (1) direct fitting of the density power spectrum to the signal, (2) using only the velocity field fluctuations in the signal, (3) looking at the signal at large enough scales such that all fluctuations trace the density field. With the foremost method, the first generation of 21 cm observations should moderately improve existing constraints on cosmological parameters for certain low-redshift reionization scenarios, and a two year observation with the second generation interferometer MWA5000 can improve constraints on Omega_w, Omega_m h^2, Omega_b h^2, Omega_nu, n_s, and alpha_s. If the Universe is substantially ionized by z = 12 or if spin temperature fluctuations are important, we show that it will be difficult to place competitive constraints on cosmological parameters with any of the considered methods.Comment: 20 pages, 12 figures, accepted by Ap