No Arabic abstract
[Abridged] The detection of the rotational lines of CO in proto-galaxies in the early Universe provides one of the most promising ways of probing the fundamental physical properties of a galaxy, such as its size, dynamical mass, gas density, and temperature. While such observations are currently limited to the most luminous galaxies, the advent of ALMA will change the situation dramatically, resulting in the detection of numerous normal galaxies at high redshifts. Maps and spectra of rotational CO line emission were calculated from a cosmological N-body/hydrodynamical TreeSPH simulation of a z ~ 3 Lyman break galaxy of UV luminosity about one order of magnitude below L*. To simulate a typical observation of our system with ALMA, we imposed characteristic noise, angular, and spectral resolution constraints. The CO line properties predicted by our simulation are in good agreement with the two Lyman break systems detected in CO to date. We find that while supernovae explosions from the ongoing star formation carve out large cavities in the molecular ISM, they do not generate large enough gas outflows to make a substantial imprint on the CO line profile. This implies that for most proto-galaxies - except possibly the most extreme cases - stellar feedback effects do not affect CO as a dynamical mass tracer. Detecting CO in sub-L* galaxies at z ~3 will push ALMA to the limits of its cababilities, and whether a source is detected or not may depend critically on its inclination angle. Both these effects (sensitivity and inclination) will severely impair the ability of ALMA to infer the gas kinematics and dynamical masses using line observations.
What type of objects are being detected as $zsim 3$ Lyman break galaxies? Are they predominantly the most massive galaxies at that epoch, or are many of them smaller galaxies undergoing a short-lived burst of merger-induced star formation? We attempt to address this question using high-resolution cosmological hydrodynamic simulations including star formation and feedback. Our $Lambda$CDM simulation, together with Bruzual-Charlot population synthesis models, reproduces the observed number density and luminosity function of Lyman break galaxies when dust is incorporated. The inclusion of dust is crucial for this agreement. In our simulation, these galaxies are predominantly the most massive objects at this epoch, and have a significant population of older stars. Nevertheless, it is possible that our simulations lack the resolution and requisite physics to produce starbursts, despite having a physical resolution of $la 700$ pc at z=3. Thus we cannot rule out merger-induced starburst galaxies also contributing to the observed population of high-redshift objects.
We perform a spectrophotometric analysis of galaxies at redshifts z = 4 - 6 in cosmological SPH simulations of a Lambda CDM universe. Our models include radiative cooling and heating by a uniform UV background, star formation, supernova feedback, and a phenomenological model for galactic winds. Analysing a series of simulations of varying boxsize and particle number allows us to isolate the impact of numerical resolution on our results. Specifically, we determine the luminosity functions in B, V, R, i, and z filters, and compare the results with observed galaxy surveys done with the Subaru telescope and the Hubble Space Telescope. We find that the simulated galaxies have UV colours consistent with observations and fall in the expected region of the colour-colour diagrams used by the Subaru group. Assuming a uniform extinction of E(B-V) = 0.15, we also find reasonable agreement between simulations and observations in the space density of UV bright galaxies at z = 3 - 6, down to the magnitude limit of each survey. For the same moderate extinction level of E(B-V) ~ 0.15, the simulated luminosity functions match observational data, but have a steep faint-end slope with alpha ~ -2.0. We discuss the implications of the steep faint-end slope found in the simulations.
We study the properties of Lyman-alpha emitters (LAEs) and Lyman-break galaxies (LBGs) at z=3-6 using cosmological SPH simulations. We investigate two simple scenarios for explaining the observed Ly-a and rest-frame UV luminosity functions (LFs) of LAEs: (i) the escape fraction scenario, in which the effective escape fraction (including the IGM attenuation) of Ly-a photons is f_Lya ~0.1 (0.15) at z=3 (6), and (ii) the stochastic scenario, in which the fraction of LAEs that are turned on at z=3 (6) is Cstoc ~0.07 (0.2) after correcting for the IGM attenuation. Our comparisons with a number of different observations suggest that the stochastic scenario is preferred over the escape fraction scenario. We find that the mean values of stellar mass, metallicity and black hole mass hosted by LAEs are all smaller in the stochastic scenario than in the escape fraction scenario. In our simulations, the galaxy stellar mass function evolves rapidly, as expected in hierarchical structure formation. However, its evolution is largely compensated by a beginning decline in the specific star formation rate, resulting in little evolution of the rest-frame UV LF from z=6 to 3. The rest-frame UV LF of both LAEs and LBGs at z=3 & 6 can be described well by the stochastic scenario provided the extinction is moderate, E(B-V) ~0.15, for both populations, although our simulation might be overpredicting the number of bright LBGs at z=6. We also discuss the correlation function and bias of LAEs. The Ly-a LFs at z=6 in a field-of-view of 0.2 deg^2 show a significantly larger scatter owing to cosmic variance relative to that in a 1 deg^2 field, and the scatter seen in the current observational estimates of the Ly-a LF can be accounted for by cosmic variance.
We investigate several fundamental properties of z ~ 4 Lyman-break galaxies by comparing observations with the predictions of a semi-analytic model based on the Cold Dark Matter theory of hierarchical structure formation. We use a sample of B_{435}-dropouts from the Great Observatories Origins Deep Survey, and complement the ACS optical B_{435}, V_{606}, i_{775}, and z_{850} data with the VLT ISAAC J, H, and K_{s} observations. We extract B_{435}-dropouts from our semi-analytic mock catalog using the same color criteria and magnitude limits that were applied to the observed sample. We find that the i_{775} - K_{s} colors of the model-derived and observed B_{435}-dropouts are in good agreement. However, we find that the i_{775}-z_{850} colors differ significantly, indicating perhaps that either too little dust or an incorrect extinction curve have been used. Motivated by the reasonably good agreement between the model and observed data we present predictions for the stellar masses, star formation rates, and ages for the z ~ 4 Lyman-break sample. We find that according to our model the color selection criteria used to select our z ~ 4 sample surveys 67% of all galaxies at this epoch down to z_{850} < 26.5. We find that our model predicts a roughly 40% mass build-up between the z ~ 4 and z ~ 3 epochs for the UV rest-frame L* galaxies. Furthermore, according to our model, at least 50% of the total stellar mass resides in relatively massive UV-faint objects that fall below our observational detection limit.
We compute the molecular line emission of massive protostellar disks by solving the equation of radiative transfer through the cores and disks produced by the recent radiation-hydrodynamic simulations of Krumholz, Klein, & McKee. We find that in several representative lines the disks show brightness temperatures of hundreds of Kelvin over velocity channels ~10 km s^-1 wide, extending over regions hundreds of AU in size. We process the computed intensities to model the performance of next-generation radio and submillimeter telescopes. Our calculations show that observations using facilities such as the EVLA and ALMA should be able to detect massive protostellar disks and measure their rotation curves, at least in the nearest massive star-forming regions. They should also detect significant sub-structure and non-axisymmetry in the disks, and in some cases may be able to detect star-disk velocity offsets of a few km s^-1, both of which are the result of strong gravitational instability in massive disks. We use our simulations to explore the strengths and weaknesses of different observational techniques, and we also discuss how observations of massive protostellar disks may be used to distinguish between alternative models of massive star formation.