No Arabic abstract
SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $approx 20%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/ddot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $dot{M}$, but reaches a maximum value at a characteristic mass $sim 2 times 10^{11} M_odot$ and a characteristic accretion rate $sim 6 times 10^2 M_odot/{rm yr}$ at $z approx 6$. We find a universal EoR luminosity-accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $phi_star$ decreases, $M_star$ is more positive (fainter), and $alpha$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{rm AB} approx 31 (32)$, it can observe $gtrsim 11 (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_star$ at $z lesssim 13 (14)$.
According to the current understanding of cosmic structure formation, the precursors of the most massive structures in the Universe began to form shortly after the Big Bang, in regions corresponding to the largest fluctuations in the cosmic density field. Observing these structures during their period of active growth and assembly - the first few hundred million years of the Universe - is challenging because it requires surveys that are sensitive enough to detect the distant galaxies that act as signposts for these structures and wide enough to capture the rarest objects. As a result, very few such objects have been detected so far. Here we report observations of a far-infrared-luminous object at redshift 6.900 (less than 800 Myr after the Big Bang) that was discovered in a wide-field survey. High-resolution imaging reveals this source to be a pair of extremely massive star-forming galaxies. The larger of these galaxies is forming stars at a rate of 2900 solar masses per year, contains 270 billion solar masses of gas and 2.5 billion solar masses of dust, and is more massive than any other known object at a redshift of more than 6. Its rapid star formation is probably triggered by its companion galaxy at a projected separation of just 8 kiloparsecs. This merging companion hosts 35 billion solar masses of stars and has a star-formation rate of 540 solar masses per year, but has an order of magnitude less gas and dust than its neighbor and physical conditions akin to those observed in lower-metallicity galaxies in the nearby Universe. These objects suggest the presence of a dark-matter halo with a mass of more than 400 billion solar masses, making it among the rarest dark-matter haloes that should exist in the Universe at this epoch.
We present measurements of carbon, oxygen, silicon, and iron in quasar absorption systems existing when the universe was roughly one billion years old. We measure column densities in nine low-ionization systems at 4.7 < z < 6.3 using Keck, Magellan, and VLT optical and near-infrared spectra with moderate to high resolution. The column density ratios among C II, O I, Si II, and Fe II are nearly identical to sub-DLAs and metal-poor ([M/H] < -1) DLAs at lower redshifts, with no significant evolution over 2 < z < 6. The estimated intrinsic scatter in the ratio of any two elements is also small, with a typical r.m.s. deviation of <0.1 dex. These facts suggest that dust depletion and ionization effects are minimal in our z > 4.7 systems, as in the lower-redshift DLAs, and that the column density ratios are close to the intrinsic relative element abundances. The abundances in our z > 4.7 systems are therefore likely to represent the typical integrated yields from stellar populations within the first gigayear of cosmic history. Due to the time limit imposed by the age of the universe at these redshifts, our measurements thus place direct constraints on the metal production of massive stars, including iron yields of prompt supernovae. The lack of redshift evolution further suggests that the metal inventories of most metal-poor absorption systems at z > 2 are also dominated by massive stars, with minimal contributions from delayed Type Ia supernovae or AGB winds. The relative abundances in our systems broadly agree with those in very metal-poor, non-carbon-enhanced Galactic halo stars. This is consistent with the picture in which present-day metal-poor stars were potentially formed as early as one billion years after the Big Bang.
We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]
The intergalactic medium (IGM) is the dominant reservoir of baryons at all cosmic epochs. We investigate the evolution of the IGM from z=2-0 in 48 Mpc/h, 110-million particle cosmological hydrodynamic simulations using three prescriptions for galactic outflows. We focus on the evolution of IGM physical properties, and how such properties are traced by Ly-alpha absorption as detectable using HST/COS. Our results broadly confirm the canonical picture that most Ly-alpha absorbers arise from highly ionized gas tracing filamentary large-scale structure. Growth of structure causes gas to move from the diffuse photoionized IGM into other cosmic phases, namely stars, cold and hot gas within galaxy halos, and the unbound and shock-heated warm-hot intergalactic medium (WHIM). By today, baryons are roughly equally divided between bound phases (35%), the diffuse IGM (41%), and the WHIM (24%). Here we (re)define the WHIM as gas with overdensities lower than that in halos and temperatures >10^5 K, in order to more closely align it with missing baryons. When we tune our photoionizing background to match the observed evolution of the Ly-alpha mean flux decrement, we obtain a line count evolution that broadly agrees with available data. We predict a column density distribution slope of -1.70 for our favored momentum-driven wind model, in agreement with recent observations, and it becomes shallower with redshift. With improved statistics, the frequency of strong lines can be a valuable diagnostic of outflows, and our favored wind model matches existing data best among our models. The relationship between column density and physical density is fairly tight from z=2-0, and evolves as rho N_HI^0.74 10^(-0.37z) for diffuse absorbers. Linewidths only loosely reflect the temperature of the absorbing gas, which will hamper attempts to quantify the WHIM using broad Ly-alpha absorbers. [Abridged]
Empirical models of galaxy formation require assumptions about the correlations between galaxy and halo properties. These may be calibrated against observations or inferred from physical models such as hydrodynamical simulations. In this Letter, we use the EAGLE simulation to investigate the correlation of galaxy size with halo properties. We motivate this analysis by noting that the common assumption of angular momentum partition between baryons and dark matter in rotationally supported galaxies overpredicts both the spread in the stellar mass-size relation and the anticorrelation of size and velocity residuals, indicating a problem with the galaxy-halo connection it implies. We find the EAGLE galaxy population to perform significantly better on both statistics, and trace this success to the weakness of the correlations of galaxy size with halo mass, concentration and spin at fixed stellar mass. Using these correlations in empirical models will enable fine-grained aspects of galaxy scalings to be matched.