Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3
aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UVsim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of zsimeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to zsimeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.
Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate
that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.
Our heuristic understanding of the abundance of dark matter halos centers around the concept of a density threshold, or barrier, for gravitational collapse. If one adopts the ansatz that regions of the linearly evolved density field smoothed on mass
scale M with an overdensity that exceeds the barrier will undergo gravitational collapse into halos of mass M, the corresponding abundance of such halos can be estimated simply as a fraction of the mass density satisfying the collapse criterion divided by the mass M. The key ingredient of this ansatz is therefore the functional form of the collapse barrier as a function of mass M or, equivalently, of the variance sigma^2(M). Several such barriers based on the spherical, Zeldovich, and ellipsoidal collapse models have been extensively discussed. Using large scale cosmological simulations, we show that the relation between the linear overdensity and the mass variance for regions that collapse to form halos by the present epoch resembles expectations from dynamical models of ellipsoidal collapse. However, we also show that using such a collapse barrier with the excursion set ansatz predicts a halo mass function inconsistent with that measured directly in cosmological simulations. This inconsistency demonstrates a failure of the excursion set ansatz as a physical model for halo collapse. We discuss implications of our results for understanding the collapse epoch for halos as a function of mass, and avenues for improving consistency between analytical models for the collapse epoch and the results of cosmological simulations.
The efficiency of star formation governs many observable properties of the cosmological galaxy population, yet many current models of galaxy formation largely ignore the important physics of star formation and the interstellar medium (ISM). Using hyd
rodynamical simulations of disk galaxies that include a treatment of the molecular ISM and star formation in molecular clouds (Robertson & Kravtsov 2008), we study the influence of star formation efficiency and molecular hydrogen abundance on the properties of high-redshift galaxy populations. In this work, we focus on a model of low-mass, star forming galaxies at 1<~z<~2 that may host long duration gamma-ray bursts (GRBs). Observations of GRB hosts have revealed a population of faint systems with star formation properties that often differ from Lyman-break galaxies (LBGs) and more luminous high-redshift field galaxies. Observed GRB sightlines are deficient in molecular hydrogen, but it is unclear to what degree this deficiency owes to intrinsic properties of the galaxy or the impact the GRB has on its environment. We find that hydrodynamical simulations of low-stellar mass systems at high-redshifts can reproduce the observed star formation rates and efficiencies of GRB host galaxies at redshifts 1<~z<~2. We show that the compact structure of low-mass high-redshift GRB hosts may lead to a molecular ISM fraction of a few tenths, well above that observed in individual GRB sightlines. However, the star formation rates of observed GRB host galaxies imply molecular gas masses of 10^8 - 10^9 M_sun similar to those produced in the simulations, and may therefore imply fairly large average H_2 fractions in their ISM.
Integral field spectroscopy of galaxies at redshift z~2 has revealed a population of early-forming, rotationally-supported disks. These high-redshift systems provide a potentially important clue to the formation processes that build disk galaxies in
the universe. A particularly well-studied example is the z=2.38 galaxy BzK-15504, which was shown by Genzel et al. (2006) to be a rotationally supported disk despite the fact that its high star formation rate and short gas consumption timescale require a very rapid acquisition of mass. Previous kinematical analyses have suggested that z~2 disk galaxies like BzK-15504 did not form through mergers because their line-of-sight velocity fields display low levels of asymmetry. We perform the same kinematical analysis on a set of simulated disk galaxies formed in gas-rich mergers of the type that may be common at high redshift, and show that the remnant disks display low velocity field asymmetry and satisfy the criteria that have been used to classify high-redshift galaxies as disks observationally. Further, we compare one of our remnants to the bulk properties of BzK-15504 and show that it has a star formation rate, gas surface density, and a circular velocity-to-velocity dispersion ratio that matches BzK-15504 remarkably well. We suggest that observations of high-redshift disk galaxies like BzK-15504 are consistent with the hypothesis that gas-rich mergers play an important role in disk formation at high redshift.
