High-redshift Lyman Break Galaxies Research Articles

Lyman Break Analogs: Constraints on the Formation of Extreme Starbursts at Low and High Redshift

AbstractLyman Break Analogs (LBAs), characterized by high far-UV luminosities and surface brightnesses as detected by GALEX, are intensely star-forming galaxies in the low-redshift universe (z ~ 0.2), with star formation rates reaching up to 50 times that of the Milky Way. These objects present metallicities, morphologies and other physical properties similar to higher redshift Lyman Break Galaxies (LBGs), motivating the detailed study of LBAs as local laboratories of this high-redshift galaxy population. We present results from our recent integral-field spectroscopy survey of LBAs with Keck/OSIRIS, which shows that these galaxies have the same nebular gas kinematic properties as high-redshift LBGs. We argue that such kinematic studies are not an appropriate diagnostic to rule out merger events as the trigger for the observed starburst. Comparison between the kinematic analysis and morphological indices from HST imaging illustrates the difficulties of properly identifying (minor or major) merger events, with no clear correlation between the results using either of the two methods. Artificial redshifting of our data indicates that this problem becomes even worse at high redshift due to surface brightness dimming and resolution loss. Whether mergers could generate the observed kinematic properties is strongly dependent on gas fractions in these galaxies. We present preliminary results of a CARMA survey for LBAs and discuss the implications of the inferred molecular gas masses for formation models.

Constraints on star-forming galaxies atz≥ 6.5 from HAWK-IY-band imaging of GOODS-South

We present the results of our search for high-redshift Lyman-break galaxies over the GOODS-South field. We use Hubble Space Telescope (HST)-ACS data in B, V, i′ & z′, Very Large Telescope (VLT)-ISAAC J and Ks, Spitzer–Infrared Array Camera (IRAC) 3.6, 4.5, 5.8 and 8.0 μm data in conjunction with the new HAWK-I Y-band science verification data to search for dropout galaxies in the redshift range 6 < z < 9. We survey ≈119 arcmin2 to YAB= 25.7 (5σ), of which 37.5 arcmin2 reaches YAB= 25.9. Candidate z′ and Y dropouts were selected on the basis of a colour cut of (Y−J)AB > 0.75 mag and (z′−Y)AB > 1.0 mag, respectively. We find no robust Y-drops (z≈ 9) brighter than JAB < 25.4. In our search for z′-band dropouts (z≈ 6.5– 7.5), we identify four possible candidates, two with z′-drop colours and clear Spitzer–IRAC detections and two less likely candidates. We also identify two previously known Galactic T-dwarf stellar contaminants with these colours, and two likely transient objects seen in the Y-band data. The implications if all or none of our candidates is real on the ultraviolet galaxy luminosity functions at z > 6.5 are explored. We find our number of z′-drop candidates to be insufficient based on the expected number of z′ drops in a simple no-evolution scenario from the z= 3 Lyman-break galaxy luminosity function but we are consistent with the observed luminosity function at z≈ 6 (if all our candidates are real). However, if one or both of our best z′-drop candidates are not z > 6.5 galaxies, this would demand evolution of the luminosity function at early epochs, in the sense that the number density of ultraviolet luminous star-forming galaxies at z > 7 is less than at z∼ 6. We show that the future surveys to be conducted with the European Southern Observatory VISTA telescope over the next 5 yr will be able to measure the bulk of the luminosity function for both z′ and Y dropouts and thus provide the strongest constraints on the level of star-formation within the epoch of reionization.

Understanding the redshift evolution of the luminosity functions of Lymanαemitters

We present a semi-analytical model of star formation which explains simultaneously the observed UV luminosity function of high redshift Lyman break galaxies (LBGs) and luminosity functions of Lyman-alpha emitters. We consider both models that use the Press-Schechter (PS) and Sheth-Tormen (ST) halo mass functions to calculate the abundances of dark matter halos. The Lyman-alpha luminosity functions at z < 4 are well reproduced with only <10% of the LBGs emitting Lyman-alpha lines with rest equivalent width greater than the limiting equivalent width of the narrow band surveys. However, the observed luminosity function at z > 5 can be reproduced only when we assume that nearly all LBGs are Lyman-alpha emitters. Thus it appears that 4 < z < 5 marks the epoch when a clear change occurs in the physical properties of the high redshift galaxies. As Lyman-alpha escape depends on dust and gas kinematics of the inter stellar medium (ISM), this could mean that on an average the ISM at z > 5 could be less dusty, more clumpy and having more complex velocity field. All of these will enable easier escape of the Lyman-alpha photons. At z > 5 the observed Lyman-alpha luminosity function are well reproduced with the evolution in the halo mass function along with very minor evolution in the physical properties of high redshift galaxies. In particular, upto z = 6.5, we do not see the effect of evolving inter galactic medium (IGM) opacity on the Lyman-alpha escape from these galaxies.

BIASES AND UNCERTAINTIES IN PHYSICAL PARAMETER ESTIMATES OF LYMAN BREAK GALAXIES FROM BROADBAND PHOTOMETRY

We investigate the biases and uncertainties in estimates of physical parameters of high-redshift Lyman break galaxies (LBGs), such as stellar mass, mean stellar population age, and star formation rate (SFR), obtained from broadband photometry. These biases arise from the simplifying assumptions often used in fitting the spectral energy distributions (SEDs). By combining ?CDM hierarchical structure formation theory, semianalytic treatments of baryonic physics, and stellar population synthesis models, we construct model galaxy catalogs from which we select LBGs at redshifts z ~ 3.4, 4.0, and 5.0. The broadband photometric SEDs of these model LBGs are then analyzed by fitting galaxy template SEDs derived from stellar population synthesis models with smoothly declining SFRs. We compare the statistical properties of LBGs' physical parameters?such as stellar mass, SFR, and stellar population age?as derived from the best-fit galaxy templates with the intrinsic values from the semianalytic model. We find some trends in these distributions: first, when the redshift is known, SED-fitting methods reproduce the input distributions of LBGs' stellar masses relatively well, with a minor tendency to underestimate the masses overall, but with substantial scatter. Second, there are large systematic biases in the distributions of best-fit SFRs and mean ages, in the sense that single-component SED-fitting methods underestimate SFRs and overestimate ages. We attribute these trends to the different star formation histories predicted by the semianalytic models and assumed in the galaxy templates used in SED-fitting procedure, and to the fact that light from the current generation of star formation can hide older generations of stars. These biases, which arise from the SED-fitting procedure, can significantly affect inferences about galaxy evolution from broadband photometry.

STUDYING LARGE- AND SMALL-SCALE ENVIRONMENTS OF ULTRAVIOLET LUMINOUS GALAXIES

Studying the environments of 0.4 < z < 1.2 ultraviolet (UV)-selected galaxies, as examples of extreme star-forming galaxies (with star formation rates (SFRs) in the range of 3-30 M_☉ yr^(–1)), we explore the relationship between high rates of star formation, host halo mass, and pair fractions. We study the large- and small-scale environments of local ultraviolet luminous galaxies (UVLGs) by measuring angular correlation functions. We cross-correlate these systems with other galaxy samples: a volume-limited sample (ALL), a blue luminous galaxy sample, and a luminous red galaxy (LRG) sample. We determine the UVLG comoving correlation length to be r_0 = 4.8^(+11.6)_(–2.4) h^(–1) Mpc at z = 1.0, which is unable to constrain the halo mass for this sample. However, we find that UVLGs form close (separation <30 kpc) pairs with the ALL sample, but do not frequently form pairs with LRGs. A rare subset of UVLGs, those with the highest FUV surface brightnesses, are believed to be local analogs of high-redshift Lyman break galaxies (LBGs) and are called Lyman break analogs (LBAs). LBGs and LBAs share similar characteristics (i.e., color, size, surface brightness, specific SFRs, metallicities, and dust content). Recent Hubble Space Telescope images of z ~ 0.2 LBAs show disturbed morphologies, signs of mergers and interactions. UVLGs may be influenced by interactions with other galaxies and we discuss this result in terms of other high star-forming, merging systems.

AN OSIRIS STUDY OF THE GAS KINEMATICS IN A SAMPLE OF UV-SELECTED GALAXIES: EVIDENCE OF “HOT AND BOTHERED” STARBURSTS IN THE LOCAL UNIVERSE

We present data from Integral Field Spectroscopy for three supercompact UV-Luminous Galaxies (ScUVLGs). As nearby (z ~ 0.2) compact (R_(50) ~ 1-2 kpc) bright Paschen-α sources, with unusually high star formation rates (SFR = 3-100 M_☉ yr^(–1)), ScUVLGs are an ideal population for studying detailed kinematics and dynamics in actively star-forming galaxies. In addition, ScUVLGs appear to be excellent analogs to high-redshift Lyman Break Galaxies (LBGs), and our results may offer additional insight into the dynamics of LBGs. Previous work by our team has shown that the morphologies of these galaxies exhibit tidal features and companions, and in this study we find that the dynamics of ScUVLGs are dominated by disturbed kinematics of the emission line gas—suggesting that these galaxies have undergone recent feedback, interactions, or mergers. While two of the three galaxies do display rotation, v/σ<1—suggesting dispersion-dominated kinematics rather than smooth rotation. We also simulate how these observations would appear at z ~ 2. Lower resolution and loss of low surface brightness features cause some apparent discrepancies between the low-z (observed) and high-z (simulated) interpretations and quantitatively gives different values for v/σ, yet simulations of these low-z analogs manage to detect the brightest regions well and resemble actual high-z observations of LBGs.

IMPROVED PHOTOMETRIC REDSHIFTS WITH SURFACE LUMINOSITY PRIORS

We apply Bayesian statistics with prior probabilities of galaxy surface luminosity (SL) to improve photometric redshifts. We apply the method to a sample of 1266 galaxies with spectroscopic redshifts in the GOODS North and South fields at 0.1 < z < 2.0. We start with spectrophotometric redshifts (SPZs) based on Probing Evolution and Reionization Spectroscopically grism spectra, which cover a wavelength range of 6000-9000A, combined with (U)BViz(JHK) broadband photometry in the GOODS fields. The accuracy of SPZ redshifts is estimated to be \sigma (\Delta(z))=0.035 with an systematic offset of -0.026, where \Delta(z)=\Delta z / (1+z), for galaxies in redshift range of 0.5 < z < 1.25. The addition of the SL prior probability helps break the degeneracy of SPZ redshifts between low redshift 4000 A break galaxies and high-redshift Lyman break galaxies which are mostly catastrophic outliers. For the 1138 galaxies at z < 1.6, the fraction of galaxies with redshift deviation \Delta (z) > 0.2 is reduced from 15.0% to 10.4%, while the rms scatter of the fractional redshift error does not change much.

Hubble Space TelescopeMorphologies of Local Lyman Break Galaxy Analogs. I. Evidence for Starbursts Triggered by Merging

Heckman and coworkers used the GALEX UV imaging survey to show that there exists a rare population of nearby compact UV-luminous galaxies (UVLGs) that closely resemble high-redshift Lyman break galaxies (LBGs). We present HST images in the UV, optical, and Hα 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 eight 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 (M_* ~ 10^(10) M☉) 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 α ~ 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.

Extinction‐corrected Star Formation Rates Empirically Derived from Ultraviolet–Optical Colors

Using a sample of galaxies from the Sloan Digital Sky Survey spectroscopic catalog with measured star-formation rates (SFRs) and ultraviolet (UV) photometry from the GALEX Medium Imaging Survey, we derived empirical linear correlations between the SFR to UV luminosity ratio and the UV-optical colors of blue sequence galaxies. The relations provide a simple prescription to correct UV data for dust attenuation that best reconciles the SFRs derived from UV and emission line data. The method breaks down for the red sequence population as well as for very blue galaxies such as the local ``supercompact'' UV luminous galaxies and the majority of high redshift Lyman Break Galaxies which form a low attenuation sequence of their own.

The Young and the Dustless: Interpreting Radio Observations of Ultraviolet‐Luminous Galaxies

Ultraviolet-luminous galaxies (UVLGs) have been identified as intensely star-forming nearby galaxies. A subset of these, the supercompact UVLGs, are believed to be local analogs of high-redshift Lyman break galaxies. Here we investigate the radio continuum properties of this important population for the first time. We have observed 42 supercompact UVLGs with the VLA, all of which have extensive coverage in the UV/optical by GALEX and SDSS. Our analysis includes comparison samples of multiwavelength data from the Spitzer First Look Survey and from the SDSS-GALEX matched catalogs. In addition we have Spitzer MIPS data for 24 of our galaxies and find that they fall on the radio-FIR correlation of normal star-forming galaxies. We find that our galaxies have lower radio to UV ratios and lower Balmer decrements than other local galaxies with similar (high) star formation rates. Optical spectra show they have lower D_n(4000) and Hδ_A indices, higher Hβ emission-line equivalent widths, and higher [O III]5007/Hβ emission-line ratios than normal star-forming galaxies. Comparing these results to galaxy spectral evolution models we conclude that supercompact UVLGs are distinguished from normal star-forming galaxies firstly by their high specific star formation rates. Moreover, compared to other types of galaxies with similar star formation rates, they have significantly less dust attenuation. In both regards they are similar to Lyman break galaxies. This suggests that the process that causes star formation in the supercompact UVLGs differs from other local star-forming galaxies, but may be similar to Lyman break galaxies.

Do Local Analogs of Lyman Break Galaxies Exist?

The optical properties of a number of supercompact ultraviolet luminous galaxies (UVLGs), recently discovered in the local universe matching Galaxy Evolution Explorer (GALEX) and Sloan Digital Sky Survey (SDSS) data, are discussed. Detailed reanalysis of the SDSS data for these and other similar but nearer galaxies shows that their surface brightness radial profile in both R and u bands is in most cases well described by an extended disk plus a central unresolved component (possibly a bulge). Since the SDSS pipeline used a single disk component to derive the half-light radius of these UVLGs, their size was severely underestimated. Consequently, the average UV surface brightness is much lower that previously quoted, casting doubts on the claim that UVLGs are the local analogs of high-redshift Lyman break galaxies.

The Role of the Dust in Primeval Galaxies: A Simple Physical Model for Lyman Break Galaxies and Lyα Emitters

We explore the onset of star formation in the early Universe, exploiting the observations of high-redshift Lyman-break galaxies (LBGs) and Lyman alpha emitters (LAEs), in the framework of the galaxy formation scenario elaborated by Granato et al. (2004) already successfully tested against the wealth of data on later evolutionary stages. Complementing the model with a simple, physically plausible, recipe for the evolution of dust attenuation in metal poor galaxies we reproduce the luminosity functions (LFs) of LBGs and of LAEs at different redshifts. This recipe yields a much faster increase with galactic age of attenuation in more massive galaxies, endowed with higher star formation rates. These objects have therefore shorter lifetimes in the LAE and LBG phases, and are more easily detected in the dusty submillimeter bright phase. The short UV bright lifetimes of massive objects strongly mitigate the effect of the fast increase of the massive halo density with decreasing redshift, thus accounting for the weaker evolution of the LBG LF, compared to that of the halo mass function, and the even weaker evolution between z~6 and z~3 of the LAE LF. LAEs are on the average expected to be younger, with lower stellar masses, and associated to less massive halos than LBGs. Finally, we show that the intergalactic medium can be completely reionized at redshift z~6-7 by massive stars shining in protogalactic spheroids with halo masses from a few 10^10 to a few 10^11 M_sun, showing up as faint LBGs with magnitude in the range -17<M_1350<-20, without resorting to any special stellar initial mass function.

Feedback in the Local Lyman‐break Galaxy Analog Haro 11 as Probed by Far‐Ultraviolet and X‐Ray Observations

We have reanalyzed FUSE data and obtained new Chandra observations of Haro 11, a local (DL = 88 Mpc) UV luminous galaxy. Haro 11 has a far-UV luminosity (1010.3 L☉), UV surface brightness (109.4 L☉ kpc-2), star formation rate (SFR), and metallicity similar to those observed in Lyman-break galaxies (LBGs). We show that Haro 11 has extended soft thermal (kT ~ 0.68 keV) X-ray emission of a luminosity and size which scales with the physical properties (e.g., SFR, stellar mass) of the host galaxy. An enhanced α/Fe ratio of ~4 relative to solar abundance suggests significant supernova enrichment. These results are consistent with the X-ray emission being produced in a shock between a supernova-driven outflow and the ambient material. The far-ultraviolet (FUV) spectra show strong absorption lines similar to those observed in LBG spectra. A blueshifted absorption component is identified as a wind outflowing at ~200-280 km s-1. O VI λλ1032,1038 emission, the dominant cooling mechanism for coronal gas at T ~ 105.5 K, is also observed. If associated with the outflow, the luminosity of the O VI emission suggests that 20% of the total mechanical energy from the supernovae and solar winds is being radiated away. This implies that radiative cooling through O VI is not significantly inhibiting the growth of the outflowing gas. In contradiction to the findings of Bergvall and coworkers, we find no convincing evidence of Lyman continuum leakage in Haro 11. We conclude that the wind has not created a tunnel allowing the escape of a significant fraction of Lyman continuum photons, and place a limit on the escape fraction of fesc 2%. Overall, both Haro 11 and a previously observed LBG analog, VV 114, provide invaluable insight into the X-ray and FUV properties of high-redshift LBGs.

The Environments of Ultrastrong MgiiAbsorbers

We present r'- or i'-band WIYN images of the fields of 15 Sloan Digital Sky Survey quasars that have spectra exhibiting intervening MgII absorption-line systems with rest equivalent widths 2.7A \le REW \le 6.0A and redshifts 0.42 < z_{abs} < 0.84. Such systems are rare and exhibit projected absorption velocity spreads in excess of \approx 300-650 km/s. Approximately 60% are expected to be damped Ly\alpha systems. In each of our fields we detect at least one galaxy that, if at the absorption redshift, would have impact parameter b \lesssim 40 kpc and luminosity L \gtrsim 0.3 L*. We measure a significant excess of galaxies at low-b to the sightlines over a large range of luminosity. Many of the sightlines are found to pass either through or close to the optically-luminous extent of a galaxy. Considering the very large velocity spreads seen in absorption, this suggests that these absorbing regions are more kinematically complex than local spirals such as the Milky Way. Our data indicate that interactions and galaxy pairs may be a contributing factor to the production of such large velocity spreads. Finally, we also find evidence that a population of galaxies with luminosities in the range 4L* \lesssim L \lesssim 13 L* may contribute to the presence of ultra-strong MgII absorption. Thus, some of the absorbing galaxies may represent a population intermediate to the very luminous high-redshift Lyman break galaxies and the fainter local starburst population.

Far‐Ultraviolet and X‐Ray Observations of VV 114: Feedback in a Local Analog to Lyman Break Galaxies

We have analyzed Far Ultraviolet Spectroscopic Explorer (FUSE), XMM-Newton, and Chandra observations of VV 114, a local galaxy merger with strong similarities to typical high-redshift Lyman break galaxies (LBGs). Diffuse thermal X-ray emission encompassing VV 114 has been observed by Chandra and XMM-Newton. This region of hot (kT ~ 0.59 keV) gas has an enhanced α/Fe element ratio relative to solar abundances and follows the same relations as typical starbursts between its properties (luminosity, size, and temperature) and those of the starburst galaxy (star formation rate, dust temperature, and galaxy mass). These results are consistent with the X-ray gas having been produced by shocks driven by a galactic superwind. The FUSE observations of VV 114 show strong, broad interstellar absorption lines with a pronounced blueshifted component (similar to what is seen in LBGs). This implies an outflow of material moving at ~300-400 km s-1 relative to VV 114. The properties of the strong O VI absorption line are consistent with radiative cooling at the interface between the hot outrushing gas seen in X-rays and the cooler material seen in the other outflowing ions in the FUSE data. We show that the wind in VV 114 has not created a tunnel that enables more than a small fraction (<few percent) of the ionizing photons from VV 114 to escape into the intergalactic medium (IGM). Taken together, these data provide a more complete physical basis for understanding the outflows that seem to be generic in LBGs. This will lead to improved insight into the role that such outflows play in the evolution of galaxies and the IGM.

The Multiwavelength Survey by Yale‐Chile (MUSYC): Survey Design and Deep Public UBVRI z ′ Images and Catalogs of the Extended Hubble Deep Field–South

We present UBVRIz' optical images taken with MOSAIC on the CTIO 4 m telescope of the 0.32 deg2 Extended Hubble Deep Field-South. This is one of four fields comprising the MUSYC survey, which is optimized for the study of galaxies at z = 3, active galactic nucleus (AGN) demographics, and Galactic structure. Our methods used for astrometric calibration, weighted image combination, and photometric calibration in AB magnitudes are described. We calculate corrected aperture photometry and its uncertainties and find through tests that these provide a significant improvement upon standard techniques. Our photometric catalog of 62,968 objects is complete to a total magnitude of RAB = 25, with R-band counts consistent with results from the literature. We select z 3 Lyman break galaxy (LBG) candidates from their UVR colors and find a sky surface density of 1.4 arcmin-2 and an angular correlation function w(θ) = (2.3 ± 1.0)θ-0.8, consistent with previous findings that high-redshift Lyman break galaxies reside in massive dark matter halos. Our images and catalogs are available online.

An Analytic Model of Galactic Winds and Mass Outflows

Galactic winds and mass outflows are observed both in nearby starburst galaxies and in high-redshift star-forming galaxies. In this paper we develop a simple analytic model to understand the observed superwind phenomenon. Our model is built upon the model of McKee & Ostriker (1977) for the interstellar medium. It allows one to predict how properties of a superwind, such as wind velocity and mass outflow rate, are related to properties of its star-forming host galaxy, such as size, gas density and star formation rate. The model predicts a threshold of star formation rate density for the generation of observable galactic winds. Galaxies with more concentrated star formation produce superwinds with higher velocities. The predicted mass outflow rates are comparable to (or slightly larger than) the corresponding star formation rates. We apply our model to both local starburst galaxies and high-redshift Lyman break galaxies, and find its predictions to be in good agreement with current observations. Our model is simple, and so can be easily incorporated into numerical simulations and semi-analytical models of galaxy formation.

The first detection of [O III] emission from high-redshift damped Lyman- galaxies

We present the detection of [O iii] emission lines from the galaxies responsible for two high-redshift z > 1.75 damped Lyman-α (DLA) absorption lines. We find two sources of [O iii] emission corresponding to the z= 1.92 DLA absorber towards the quasar Q 2206−1958, and we also detect [O iii] emission from the galaxy responsible for the z= 3.10 DLA absorber towards the quasar 2233.9+1381. These are the first detections of rest-frame optical emission lines from high-redshift DLA galaxies. Unlike the Lyα line, the [O iii] line provides a measure of the systemic velocity of the galaxy. We compare the [O iii] redshifts with the velocity profile of the low-ionization metal lines in these two absorbers, with the goal of distinguishing between the model of Prochaska and Wolfe of DLA absorbers as large rapidly rotating cold thick discs, and the standard hierarchical cold dark matter model of structure formation, in which DLAs arise in protogalactic fragments. We find some discrepancies with the predictions of the former model. Furthermore, the image of the DLA galaxy towards Q 2206−1958 shows a complex disturbed morphology, which is more in accord with the hierarchical picture. We use the properties of the rest-frame optical emission lines to further explore the question posed by Moller et al.: are high-redshift DLA galaxies Lyman-break galaxies (LBGs) selected by gas cross-section? The measured velocity dispersions of the DLA galaxies are in agreement with this picture, while the data on the [O iii] luminosities and the velocity differences between the Lyα and [O iii] lines are inconclusive, as there are insufficient LBG measurements overlapping in luminosity. Finally, we estimate the star formation rates in these two DLA galaxies, using a variety of diagnostics, and include a discussion of the extent to which the [O iii] line is useful for this purpose.

A deeper view of extremely red galaxies: the redshift distribution in the GOODS/CDFS ISAAC field

We have analysed 5-epoch GOODS HST-ACS B, V, I_775 and z datasets (V1.0 release), in conjunction with existing VLT-ISAAC imaging in the J, H and Ks bands, to derive estimated redshifts for the sample of 198 Extremely Red Galaxies (ERGs) with Ks 3.92 selected by Roche et al.(2003) from 50.4 sq.arcmin of the GOODS/CDFS field. We find that, at this depth, the ERG population spans the redshift range 0.5 =1.0 to =2.5 and compare it with the global Ks-band LF at redshifts 1 =2.0 to =2.5 and we connect this fact with the presence of progenitors of the local L>L* population at redshifts z_phot>2. We determine lower limits of rho_c=(6.1 +/- 1.9) x 10^{-5} Mpc^{-3} and rho_c=(2.1 +/- 1.1) x 10^{-5}Mpc^{-3} on the comoving densities of progenitors of local massive galaxies already assembled at redshifts =2.5 and =3.5, respectively. We have investigated the existence of high-redshift Lyman break galaxies massive enough to be included in this ERG sample. Out of an initial list of 12 potential very high redshift candidates, we have identified 2 ERGs which have a high probability of lying at z_phot>4. We discuss the advantages of multi-colour to single-colour selection techniques in obtaining reliable lists of very high-redshift candidate sources, and present revised lower redshift estimates for sources previously claimed as potential z> 5 dropouts in recent studies.(abridged)

The nature of high-redshift galaxies

Using semi-analytic models of galaxy formation set within the Cold Dark Matter (CDM) merging hierarchy, we investigate several scenarios for the nature of the high-redshift ($z \ga 2$) Lyman-break galaxies (LBGs). We consider a ``collisional starburst'' model in which bursts of star formation are triggered by galaxy-galaxy mergers, and find that a significant fraction of LBGs are predicted to be starbursts. This model reproduces the observed comoving number density of bright LBGs as a function of redshift and the observed luminosity function at $z\sim 3$ and $z\sim 4$, with a reasonable amount of dust extinction. Model galaxies at $z \sim 3$ have star formation rates, half-light radii, I-K colours, and internal velocity dispersions that are in good agreement with the data. Global quantities such as the star formation rate density and cold gas and metal content of the Universe as a function of redshift also agree well. Two ``quiescent'' models without starbursts are also investigated. In one, the star formation efficiency in galaxies remains constant with redshift, while in the other, it scales inversely with disc dynamical time, and thus increases rapidly with redshift. The first quiescent model is strongly ruled out as it does not produce enough high redshift galaxies once realistic dust extinction is accounted for. The second quiescent model fits marginally, but underproduces cold gas and very bright galaxies at high redshift. A general conclusion is that star formation at high redshift must be more efficient than locally. The collisional starburst model appears to accomplish this naturally without violating other observational constraints.