The properties of Lyα emitting galaxies in hierarchical galaxy formation models

We present detailed predictions for the properties of Lyα-emitting galaxies in the framework of the Λ cold dark matter cosmology, calculated using the semi-analytical galaxy formation model galform . We explore a model that assumes a top-heavy initial mass function in starbursts and that has previously been shown to explain the sub-millimetre number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Lyα photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Lyα emitters (LAEs) over the redshift range  3 < z < 6.6  . We also examine the distribution of Lyα equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of LAEs, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of LAEs at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of LAEs at   z > 7  , a redshift range that is starting to be be probed by near-infrared surveys and using new instruments such as the Dark Ages Z Lyman Explorer (DAzLE).     

Double starbursts triggered by mergers in hierarchical clustering scenarios

We use cosmological smooth particle hydrodynamical (SPH) simulations to study the effects of mergers in the star formation history of galactic objects in hierarchical clustering scenarios. We find that during some merger events, gaseous discs can experience two starbursts: the first one during the orbital decay phase, owing to gas inflows driven as the satellite approaches, and the second one when the two baryonic clumps collide. A trend for these first induced starbursts to be more efficient at transforming the gas into stars is also found. We detect that systems that do not experience early gas inflows have well-formed stellar bulges and more concentrated potential wells, which seem to be responsible for preventing further gas inward transport triggered by tidal forces. The potential wells concentrate owing to the accumulation of baryons in the central regions and of dark matter as the result of the pulling in by baryons. The coupled evolution of the dark matter and baryons would lead to an evolutionary sequence during which systems with shallower total potential wells suffer early gas inflows during the orbital decay phase that help to feed their central mass concentration, pulling in dark matter and contributing to build up more stable systems. Within this scenario, starbursts triggered by early gas inflows are more likely to occur at early stages of evolution of the systems and to be an important contributor to the formation of stellar bulges. Our results constitute the first proof that bulges can form as the product of collapse, collisions and secular evolution in a cosmological framework, and they are consistent with a rejuvenation of the stellar population in bulges at intermediate z with, at least, 50 per cent of the stars (in SCDM) being formed at high z .     

Clues for the origin of the fundamental metallicity relations – I. The hierarchical building up of the structure

We analyse the evolutionary history of galaxies formed in a hierarchical scenario consistent with the concordance Lambda cold dark matter (ΛCDM) model focusing on the study of the relation between their chemical and dynamical properties. Our simulations consistently describe the formation of the structure and its chemical enrichment within a cosmological context. Our results indicate that the luminosity–metallicity and the stellar mass–metallicity (LZR and MZR) relations are naturally generated in a hierarchical scenario. Both relations are found to evolve with redshift. In the case of the MZR, the estimated evolution is weaker than that deduced from observational works by approximately 0.10 dex. We also determine a characteristic stellar mass, M _c≈ 3 × 10¹⁰ M_⊙, which segregates the simulated galaxy population into two distinctive groups and which remains unchanged since z ∼ 3, with a very weak evolution of its metallicity content. The value and role played by M _c is consistent with the characteristic mass estimated from the SDSS galaxy survey by Kauffmann et al. Our findings suggest that systems with stellar masses smaller than M _c are responsible for the evolution of this relation at least from z ≈ 3. Larger systems are stellar dominated and have formed more than 50 per cent of their stars at   z ≥ 2  , showing very weak evolution since this epoch. We also found bimodal metallicity and age distributions from z ∼ 3, which reflects the existence of two different galaxy populations. Although SN feedback may affect the properties of galaxies and help to shape the MZR, it is unlikely that it will significantly modify M _c since, from   z = 3  this stellar mass is found in systems with circular velocities larger than 100 km s⁻¹.     

Star formation and chemical evolution of damped Lyman α systems

In this paper, we investigate the star formation and chemical evolution of damped Lyman α systems (DLAs) based on the disc galaxy formation model developed by Mo, Mao & White. We propose that the DLAs are the central galaxies of less-massive dark haloes present at redshifts z ∼3, and they should inhabit haloes of moderately low circular velocity. The empirical Schmidt law of star formation rates, and closed box model of chemical evolution that an approximation known as instantaneous recycling is assumed, are adopted. In our models, when the predicted distribution of metallicity for DLAs is calculated, two cases are considered. One is that, using the closed-box model, empirical Schmidt law and star formation time, the distribution of metallicity can be directly calculated. The other is that, when the simple gravitational instability of a thin isothermal gas disc as first discussed by Toomre is considered, the star formation occurs only in the region where the surface density of gas satisfies the critical value, not everywhere of a gas disc. In this case, we first obtain the region where the star formation can occur by assuming that the disc has a flat rotation curve and rotational velocity is equal to the circular velocity of the surrounding dark matter halo, and then calculate the metallicity distribution as in case one. We assume that star formation in each DLA lasts for a period of 1 Gyr from redshifts z =3. There is only one output parameter in our models, i.e. the stellar yield, which relates to the time of star formation history and is obtained by normalizing the predicted distribution of metallicity to the mean value of 1/13 Z_⊙ as presented by Pettini et al.. The predicted metallicity distribution is consistent with the current (rather limited) observational data. A random distribution of galactic discs is taken into account.     

Molecular hydrogen in damped Lyα systems: clues to interstellar physics at high redshift

In order to interpret H₂ quasar absorption-line observations of damped Lyα systems (DLAs) and subDLAs, we model their H₂ abundance as a function of dust-to-gas ratio, including H₂ self-shielding and dust extinction against dissociating photons. Then, we constrain the physical state of the gas by using H₂ data. Using H₂ excitation data for DLAs with H₂ detections, we derive a gas density  1.5 ≲ log n (cm⁻³) ≲ 2.5  , temperature  1.5 ≲ log T (K) ≲ 3  , and an internal ultraviolet (UV) radiation field (in units of the Galactic value)  0.5 ≲ log χ≲ 1.5  . We then find that the observed relation between the molecular fraction and the dust-to-gas ratio of the sample is naturally explained by the above conditions. However, it is still possible that H₂ deficient DLAs and subDLAs with H₂ fractions less than  ∼10⁻⁶  are in a more diffuse and warmer state. The efficient photodissociation by the internal UV radiation field explains the extremely small H₂ fraction  (≲10⁻⁶)  observed for  κ≲ 1/30  (κ is the dust-to-gas ratio in units of the Galactic value); H₂ self-shielding causes a rapid increase in, and large variations of, H₂ abundance for  κ≳ 1/30  . We finally propose an independent method to estimate the star formation rates of DLAs from H₂ abundances; such rates are then critically compared with those derived from other proposed methods. The implications for the contribution of DLAs to the cosmic star formation history are briefly discussed.     

NICMOS imaging search for high-redshift damped Lyα galaxies

We are engaged in a programme of imaging with the STIS and NICMOS (NIC2) instruments aboard the Hubble Space Telescope ( HST ), to search for the galaxy counterparts of 18 high-redshift z >1.75 damped Ly α absorption lines and five Lyman-limit systems seen in the spectra of 16 target quasars. This paper presents the results of the imaging campaign with the NIC2 camera. We describe the steps followed in reducing the data and combining in mosaics, and the methods used for subtracting the image of the quasar in each field, and for constructing error frames that include the systematic errors associated with the PSF subtraction. To identify candidate counterparts, which are either compact or diffuse, we convolved the image and variance frames with circular top-hat filters of diameter 0.45 and 0.90 arcsec respectively, to create frames of summed S /N within the aperture. For each target quasar we provide catalogues listing positions and aperture magnitudes of all sources within a square of side 7.5 arcsec centred on the quasar, detected at S / N >6 . We find a total of 41 candidates, of which three have already been confirmed spectroscopically as the counterparts. We provide the aperture magnitude detection limits as a function of impact parameter, for both detection filters, for each field. The average detection limit for compact (diffuse) sources is H _AB =25.0 (24.4) at an angular separation of 0.56 (0.79) arcsec from the quasar, improving to H _AB =25.5 (24.8) at large angular separations. For the brighter sources we have measured the half-light radius and the n parameter of the best-fitting deconvolved Sersic-law surface brightness profile, and the ellipticity and orientation.     

Constraining quasar host halo masses with the strength of nearby Lyα forest absorption

Using cosmological hydrodynamic simulations, we measure the mean transmitted flux in the Lyα forest for quasar sightlines that pass near a foreground quasar. We find that the trend of absorption with pixel quasar separation distance can be fitted using a simple power-law form including the usual correlation function parameters r ₀ and γ, so that     . From the simulations, we find the relation between r ₀ and quasar host mass, and formulate this as a way to estimate quasar host dark matter halo masses, quantifying uncertainties due to cosmological and IGM parameters, and redshift errors. With this method, we examine data for ∼9000 quasars from the Sloan Digital Sky Survey (SDSS) Data Release 5, assuming that the effect of ionizing radiation from quasars (the so-called transverse proximity effect) is unimportant (no evidence for it is seen in the data). We find that the best-fitting host halo mass for SDSS quasars with mean redshift z = 3 and absolute G -band magnitude −27.5 is  log  M /M_⊙= 12.68^+0.81_−0.67  . We also use the Lyman-Break Galaxy (LBG) and Lyα forest data of Adelberger et al. in a similar fashion to constrain the halo mass of LBGs to be  log₁₀  M /M_⊙= 11.41^+0.54_−0.59  , a factor of ∼20 lower than the bright quasars. In addition, we study the redshift distortions of the Lyα forest around quasars, using the simulations. We use the quadrupole to monopole ratio of the quasar Lyα forest correlation function as a measure of the squashing effect. We find its dependence on halo mass difficult to measure, but find that it may be useful for constraining cosmic geometry.