Star Formation in a Cosmological Simulation of Reionization

We study the luminosity functions of high-redshift galaxies in detailed hydrodynamic simulations of cosmic reionization, which are designed to reproduce the evolution of the Lyman-alpha forest between z=5 and z=6. We find that the luminosity functions and total stellar mass densities are in agreement with observations when plausible assumptions about reddenning at z=6 are made. Our simulations support the conclusion that stars alone reionized the universe.Comment: Accepted for publication in Ap

Distribution of Damped Lyman-alpha Absorbers in a Lambda Cold Dark Matter Universe

We present the results of a numerical study of a galactic wind model and its implications on the properties of damped Lyman-alpha absorbers (DLAs) using cosmological hydrodynamic simulations. We vary both the wind strength and the internal parameters of the the wind model in a series of cosmological SPH simulations that include radiative cooling and heating by a UV background, star formation, and feedback from supernovae and galactic winds. To test our simulations, we examine the DLA `rate-of-incidence' as a function of halo mass, galaxy apparent magnitude, and impact parameter. We find that the statistical distribution of DLAs does not depend on the exact values of internal numerical parameters that control the decoupling of hydrodynamic forces when the gas is ejected from starforming regions. The DLA rate-of-incidence in our simulations at z=3 is dominated by the faint galaxies with apparent magnitude R_AB < 25.5. However, interestingly in a `strong wind' run, the differential distribution of DLA sight-lines is peaked at Mhalo = 10^{12} Msun/h (R_AB~27), and the mean DLA halo mass is Mmean=10^{12.4} Msun/h (R_AB ~ 26). These mass-scales are much larger than those if we ignore winds, because galactic wind feedback suppresses the DLA cross section in low-mass halos and increases the relative contribution to the DLA incidence from more massive halos. The DLAs in our simulations are more compact than the present-day disk galaxies, and the impact parameter distribution is very narrow unless we limit the search for the host galaxy to only bright LBGs. The comoving number density of DLAs is higher than that of LBGs down to R_AB=30 mag if the physical radius of each DLA is smaller than 5 kpc/h_70. We discuss conflicts between current simulations and observations, and potential problems with simulations based on the CDM model.Comment: 37 pages, 11 figures. Accepted to ApJ. Additional numerical tests of the internal parameters of the galactic wind model are presente

Detectability of [CII] 158 micron emission from high-redshift galaxies: predictions for ALMA and SPICA

We discuss the detectability of high-redshift galaxies via [CII] 158 micron line emission by coupling an analytic model with cosmological Smoothed Particle Hydrodynamics (SPH) simulations that are based on the concordance Lambda cold dark matter (CDM) model. Our analytic model describes a multiphase interstellar medium irradiated by the far ultra-violet radiation from local star-forming regions, and it calculates thermal and ionization equilibrium between cooling and heating. The model allows us to predict the mass fraction of a cold neutral medium (CNM) embedded in a warm neutral medium (WNM). Our cosmological SPH simulations include a treatment of radiative cooling/heating, star formation, and feedback effects from supernovae and galactic winds. Using our method, we make predictions for the [CII] luminosity from high-redshift galaxies which can be directly compared with upcoming observations by the Atacama Large Millimeter Array (ALMA) and the Space Infrared Telescope for Cosmology and Astrophysics (SPICA). We find that the number density of high-redshift galaxies detectable by ALMA and SPICA via [CII] emission depends significantly on the amount of neutral gas which is highly uncertain. Our calculations suggest that, in a CDM universe, most [CII] sources at z=3 are faint objects with \Snu < 0.01 mJy. Lyman-break galaxies (LBGs) brighter than R_AB=23.5 mag are expected to have flux densities \Snu = 1-3 mJy depending on the strength of galactic wind feedback. The recommended observing strategy for ALMA and SPICA is to aim at very bright LBGs or star-forming DRG/BzK galaxies.Comment: 39 pages, 13 figures, accepted to ApJ. Matched to the accepted versio

HST morphologies of local Lyman break galaxy analogs I: Evidence for starbursts triggered by merging

Heckman et al. (2005) used the Galaxy Evolution Explorer (GALEX) UV imaging survey to show that there exists a rare population of nearby compact UV-luminous galaxies (UVLGs) that closely resembles high redshift Lyman break galaxies (LBGs). We present HST images in the UV, optical, and Ha, and resimulate them at the depth and resolution of the GOODS/UDF fields to show that the morphologies of UVLGs are also similar to those of LBGs. Our sample of 8 LBG analogs thus provides detailed insight into the connection between star formation and LBG morphology. Faint tidal features or companions can be seen in all of the rest-frame optical images, suggesting that the starbursts are the result of a merger or interaction. The UV/optical light is dominated by unresolved (~100-300 pc) super starburst regions (SSBs). A detailed comparison with the galaxies Haro 11 and VV 114 at z=0.02 indicates that the SSBs themselves consist of diffuse stars and (super) star clusters. The structural features revealed by the new HST images occur on very small physical scales and are thus not detectable in images of high redshift LBGs, except in a few cases where they are magnified by gravitational lensing. We propose, therefore, that LBGs are mergers of gas-rich, relatively low-mass (~10^10 Msun) systems, and that the mergers trigger the formation of SSBs. If galaxies at high redshifts are dominated by SSBs, then the faint end slope of the luminosity function is predicted to have slope alpha~2. Our results are the most direct confirmation to date of models that predict that the main mode of star formation in the early universe was highly collisional.Comment: 32 pages, 15 figures. ApJ In pres

Ice emission and the redshifts of submillimeter sources

Observations at submillimeter wavelengths have revealed a population of sources thought to be at relatively large redshifts. The position of the 850 $\mu$m passband on the Rayleigh-Jeans portion of the Planck function leads to a maximum redshift estimate of $z\sim$4.5 since sources will not retain their redshift independent brightness close to the peak of the Planck function and thus drop out of surveys. Here we review evidence that ice absorption is present in the spectra of local ultraluminous infrared galaxies which are often taken as analogs for the 850 $\mu$m source population. We consider the implication of this absorption for ice induced spectral structure at far infrared wavelengths and present marginal astronomical evidence that amorphous ice may have a feature similar to crystalline ice near 150 $\mu$m. Recent corroborative laboratory evidence is supportive of this conclusion. It is argued that early metal enrichment by pair instability SN may lead to a high ice content relative to refractory dust at high redshift and a fairly robust detection of ice emission in a $z=6.42$ quasar is presented. It is further shown that ice emission is needed to understand the 450 $\mu$m sources observed in the GOODS-N field. We are thus encouraged to apply far infrared ice emission models to the available observations of HDF 850.1, the brightest submillimeter source in the {\it Hubble Deep Field}. We suggest that a redshift as large as 13 may need to be considered for this source, nearly a factor of three above the usual top estimate. Inclusion of the possibility of far infrared ice emission in the spectral energy distributions of model sources generally broadens the range of redshifts to be considered for submillimeter sources compared to models without ice emission.Comment: 37 pages, 8 figures, accepted for publication in the Astrophysical Journa

Feedback-Driven Evolution of the Far-Infrared Spectral Energy Distributions of Luminous and Ultraluminous Infrared Galaxies

We calculate infrared spectral energy distributions (SEDs) from simulations of major galaxy mergers and study the effect of AGN and starburst driven feedback on the evolution of the SED as a function of time. We use a self-consistent three-dimensional radiative equilibrium code to calculate the emergent SEDs and to make images. To facilitate a simple description of our findings, we describe our results in reference to an approximate analytic solution for the far-IR SED. We focus mainly on the luminous infrared galaxy (LIRG) and ultraluminous infrared galaxy (ULIRG) phases of evolution. We contrast the SEDs of simulations performed with AGN feedback to simulations performed with starburst driven wind feedback. We find that the feedback processes critically determine the evolution of the SED. Changing the source of illumination (whether stellar or AGN) has virtually no impact on the reprocessed far-infrared SED. We find that AGN feedback is particularly effective at dispersing gas and rapidly injecting energy into the ISM. The observational signature of such powerful feedback is a warm SED. In general, simulations performed with starburst driven winds have colder spectra and reprocess more of their emission into the infrared, resulting in higher infrared to bolometric luminosities compared to (otherwise equivalent) simulations performed with AGN feedback. We depict our results in IRAS bands, as well as in Spitzer's MIPS bands, and in Herschel's PACS bands.Comment: Accepted to ApJ, 11 figures, 26 pages, added 2 figure

Adverse impact of banana Xanthomonas Wilt on farmers’ livelihoods in Eastern and Central Africa

Banana is a key crop in the livelihoods of many people in the Great Lakes region of East and Central Africa. For more than a decade now, the crop has been threatened by Banana Xanthomonas Wilt (BXW) which has spread throughout the region but at different rates. The disease attacks all banana cultivars and can cause up to 100% yield losses at farm level if effective control measures are not put in place. However, limited information on impact of BXW at regional level is available to guide interventions. Thus, this study assessed the impact of BXW on farmers’ livelihoods in Kagera basin of Tanzania, Burundi and Rwanda. A total of 436 households (Tanzania 120, Burundi 208 and Rwanda 108) mostly from major banana-producing and BXW-affected districts were sampled and interviewed in a household survey. Thirty-three to seventy-five of the total banana mats per farm in the three countries were infected with BXW. Banana production losses caused by BXW were valued at US$10.2 million and US$ 2.95 million in Tanzania and Rwanda, respectively, banana sales by farmers dropped by 35% while bunch prices unpredictably doubled. Since banana is a key component of these farming communities, the banana production losses resulted in significant reduction in household food security and incomes. To cope with these challenges, most households are diversifying into other food crops such as maize, cassava and sweet potatoes. This poses a number of socio-economic and biological implications that require further investigation

Searching for massive galaxies at z>=3.5 in GOODS-North

We constrain the space density and properties of massive galaxy candidates at redshifts of z>=3.5 in the GOODS-N field. By selecting sources in the Spitzer+IRAC bands, a highly stellar mass-complete sample is assembled,including massive galaxies which are very faint in the optical/near-IR bands that would be missed by samples selected at shorter wavelengths. The z>=3.5 sample was selected down to 23 mag at 4.5 micron using photometric redshifts that have been obtained by fitting the galaxies SEDs at optical, near-IR and IRAC bands. We also require that the brightest band in which candidates are detected is the IRAC 8 micron band in order to ensure that the near-IR 1.6 micron (rest-frame) peak is falling in or beyond this band. We found 53 z>=3.5 candidates, with masses in the range of M~10^10-10^11M_sun. At least ~81% of these galaxies are missed by traditional Lyman Break selection methods based on UV light. Spitzer+MIPS emission is detected for 60% of the sample of z>=3.5 galaxy candidates. Although in some cases this might suggest a residual contamination from lower redshift star-forming galaxies or AGN, 37% of these objects are also detected in the sub-mm/mm bands in recent SCUBA,AzTEC and MAMBO surveys, and have properties fully consistent with vigorous starburst galaxies at z>=3.5. The comoving number density of galaxies with stellar masses >= 5x10^10M_sun(a reasonable stellar mass completeness limit for our sample) is 2.6x10^-5Mpc^-3 (using the volume within 3.5<z<5), and the corresponding stellar mass density is ~2.9x10^6M_sunMpc^-3, or~3% of the local density above the same stellar mass limit.For the sub-sample of MIPS-undetected galaxies,we find a number density of ~0.97x10^-5Mpc^-3 and a stellar mass density of ~1.15x10^6M_sun Mpc^-3.[abridged]Comment: Accepted by A&A; 35 pages, 15 figures, references update

Photometric Properties of the Most Massive High-Redshift Galaxies

We calculate the observable properties of the most massive high-redshift galaxies in the hierarchical formation scenario where stellar spheroid and supermassive black hole growth are fueled by gas-rich mergers. Combining high-resolution hydrodynamical simulations of the hierarchical formation of a z~6 quasar, stellar population synthesis models, template AGN spectra, prescriptions for interstellar and intergalactic absorption, and the response of modern telescopes, the photometric evolution of galaxies destined to host z~6 quasars are modeled at redshifts z~4-14. These massive galaxies, with enormous stellar masses of M_star ~10^11.5-10^12 M_sun. and star formation rates of SFR~10^3-10^4 M_sun yr^-1 at z>~7, satisfy a variety of photometric selection criteria based on Lyman-break techniques including V-band dropouts at z>~5, i-band dropouts at z>~6, and z-band dropouts at z>~7. The observability of the most massive high-redshift galaxies is assessed and compared with a wide range of existing and future photometric surveys including SDSS, GOODS/HUDF, NOAO WDFS, UKIDSS, the IRAC Shallow Survey, Pan-STARRS, LSST, and SNAP. Massive stellar spheroids descended from z~6 quasars will likely be detected at z~4 by existing surveys, but owing to their low number densities the discovery of quasar progenitor galaxies at z>7 will likely require future surveys of large portions of the sky (>~0.5%) at wavelengths lambda>1 micron. The detection of rare, star-bursting, massive galaxies at redshifts z>~6 would provide support for the hierarchical formation of the earliest quasars and characterize the primitive star-formation histories of the most luminous elliptical galaxies.Comment: Version accepted by ApJ (5/21/2007), for high-resolution figures see http://kicp.uchicago.edu/~brant/astro-ph/massive_high_redshift_galaxies

The dusty, albeit ultraviolet bright infancy of galaxies

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMAComment: 5 pages, 3 figures, submitted to MNRAS-le