# 1 Introduction

## 1 Introduction

In the last few years, the remarkable progress in our understanding of faint galaxy data made possible by the combination of HST deep imaging  and ground-based spectroscopy , ,  has permitted to shed some light on the evolution of the stellar birthrate in the universe, to tentatively identify the epoch $`1\mathrm{}<z\mathrm{}<2`$ where most of the optical extragalactic background light was produced, and to set important contraints on galaxy evolution scenarios , , , . The explosion in the quantity of information available on the high-redshift universe at optical wavelengths has been complemented by the detection of the far-IR/sub-mm background by DIRBE and FIRAS , , that has revelead the optically ‘hidden’ side of galaxy formation, and shown that a significant fraction of the energy released by stellar nucleosynthesis is re-emitted as thermal radiation by dust. The underlying goal of all these efforts is to understand the growth of cosmic structures and the mechanisms that shaped the Hubble sequence, and ultimately to map the transition from the cosmic ‘dark age’  to a ionized universe populated with luminous sources. While one of the important questions recently emerged is the nature (starbursts or AGNs?) and redshift distribution of the ultraluminous sub-mm sources discovered by SCUBA , , , of perhaps equal interest is the possible existence of a large population of faint galaxies still undetected at high-$`z`$, as the color-selected ground-based and Hubble Deep Field (HDF) samples include only the brightest and bluest star-forming objects. In hierarchical clustering cosmogonies, high-$`z`$ dwarfs and/or mini-quasars (i.e. an early generation of stars and accreting black holes in dark matter halos with circular velocities $`v_c50\mathrm{km}\mathrm{s}^1`$) may actually be one of the main source of UV photons and heavy elements at early epochs , , .

In this talk I will focus on some of the open issues and controversies surrounding our present understanding of the history of the conversion of cold gas into stars within galaxies, and of the evolution with cosmic time of the space density of luminous sources. An Einstein-de Sitter universe ($`q_0=0.5`$) with $`H_0=50h_{50}\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$ will be adopted in the following.

## 2 Extragalactic background light

The extragalactic background light (EBL) is an indicator of the total luminosity of the universe. It provides unique information on the evolution of cosmic structures at all epochs, as the cumulative emission from galactic systems and AGNs is expected to be recorded in this background. Figure 1 shows the optical EBL from known galaxies together with the recent COBE results. The value derived by integrating the galaxy counts  down to very faint magnitude levels \[because of the flattening at faint magnitudes of the $`N(m)`$ differential counts most of the contribution to the optical EBL comes from relatively bright galaxies\] implies a lower limit to the EBL intensity in the 0.3–2.2 $`\mu `$m interval of $`I_{\mathrm{opt}}12\mathrm{nW}\mathrm{m}^2\mathrm{sr}^1`$. When combined with the FIRAS and DIRBE measurements ($`I_{\mathrm{FIR}}16\mathrm{nW}\mathrm{m}^2\mathrm{sr}^1`$ in the 125–5000 $`\mu `$m range), this gives an observed EBL intensity in excess of $`28\mathrm{nW}\mathrm{m}^2\mathrm{sr}^1`$. The correction factor needed to account for the residual emission in the 2.2 to 125 $`\mu `$m region is probably $`\mathrm{}<2`$ . We shall see below how a population of dusty AGNs could make a significant contribution to the FIR background. In this talk I shall adopt a conservative reference value for the total EBL intensity associated with star formation activity over the entire history of the universe of $`I_{\mathrm{EBL}}=40I_{40}\mathrm{nW}\mathrm{m}^2\mathrm{sr}^1`$.

## 3 Star formation history

It has become familiar to interpret recent observations of high-redshift sources via the comoving volume-averaged history of star formation. This is the mean over cosmic time of the stochastic, possibly short-lived star formation episodes of individual galaxies, and follows a relatively simple dependence on redshift. Its latest version, uncorrected for dust extinction, is plotted in Figure 2 (left panel). The measurements are based upon the rest-frame UV luminosity function (at 1500 and 2800 Å), assumed to be from young stellar populations . The prescription for a ‘correct’ de-reddening of these values has been the subject of an ongoing debate. Dust may play a role in obscuring the UV continuum of Canada-France Reshift Survey (CFRS, $`0.3<z<1`$) and Lyman-break ($`z3`$) galaxies, as their colors are too red to be fit with an evolving stellar population and a Salpeter initial mass function (IMF) . The fiducial model of  had an upward correction factor of 1.4 at 2800 Å, and 2.1 at 1500 Å. Much larger corrections have been argued for by  ($`\times 10`$ at $`z=1`$),  ($`\times 15`$ at $`z=3`$), and  ($`\times 16`$ at $`z>2`$). As noted already by  and , a consequence of such large extinction values is the possible overproduction of metals and red light at low redshifts. Most recently, the evidence for more moderate extinction corrections has included measurements of star-formation rates (SFR) from Balmer lines by  ($`\times 2`$ at $`z=0.2`$),  ($`\times 3.1\pm 0.4`$ at $`z=1`$), and  ($`\times 26`$ at $`z=3`$). ISO follow-up of CFRS fields  has shown that the star-formation density derived by FIR fluxes ($`\times 2.3\pm 0.7`$ at $`0z1`$) is about 3.5 times lower than in . Figure 2 (right panel) depicts an extinction-corrected (with $`A_{1500}=1.2`$ mag, 0.4 mag higher than in ) version of the same plot. The best-fit cosmic star formation history (shown by the dashed-line) produces a total EBL of $`37\mathrm{nW}\mathrm{m}^2\mathrm{sr}^1`$. About 60% of this is radiated in the UV$`+`$optical$`+`$near-IR between 0.1 and 5 $`\mu `$m; the total amount of starlight that is absorbed by dust and reprocessed in the far-IR is $`13\mathrm{nW}\mathrm{m}^2\mathrm{sr}^1`$. Because of the uncertainties associated with the incompleteness of the data sets, photometric redshift technique, dust reddening, and UV-to-SFR conversion, these numbers are only meant to be indicative. On the other hand, the model is not in obvious disagreement with any of the observations, and is able, in particular, to provide a reasonable estimate of the near-IR luminosity density in the range $`0\mathrm{}<z\mathrm{}<1`$ (Fig. 3).

## 4 The stellar mass density today

With the help of some simple stellar population synthesis tools it is possible at this stage to make an estimate of the integrated stellar mass density today. The total bolometric luminosity of a simple stellar population (a single generation of coeval stars) having mass $`M`$ can be well approximated by a power-law with time for all ages $`t\mathrm{}>100`$ Myr,

$$L(t)=1.3L_{}\frac{M}{M_{}}\left(\frac{t}{1\mathrm{Gyr}}\right)^{0.8}$$
(1)   
(cf. ), where we have assumed solar metallicity and a Salpeter IMF truncated at 0.1 and 125 $`M_{}`$. In a stellar system with arbitrary star-formation rate per unit cosmological volume, $`\dot{\rho }_{}`$, the comoving bolometric emissivity at time $`t`$ is given by the convolution integral

$$\rho _{\mathrm{bol}}(t)=_0^tL(\tau )\dot{\rho }_{}(t\tau )𝑑\tau .$$
(2)   
The total background light observed at Earth ($`t=t_H`$) is

$$I_{\mathrm{EBL}}=\frac{c}{4\pi }_0^{t_H}\frac{\rho _{\mathrm{bol}}(t)}{1+z}𝑑t,$$
(3)   
where the factor $`(1+z)`$ at the denominator is lost to cosmic expansion when converting from observed to radiated luminosity density. From the above equations it is easy to derive

$$I_{\mathrm{EBL}}=740\mathrm{nW}\mathrm{m}^2\mathrm{sr}^1\frac{\dot{\rho }_{}}{\mathrm{M}_{}\mathrm{yr}^1\mathrm{Mpc}^3}\left(\frac{t_H}{13\mathrm{Gyr}}\right)^{1.87}.$$
(4)   
The observations shown in Figure 1 therefore imply a “fiducial” mean star formation density of $`\dot{\rho }_{}=0.054I_{40}\mathrm{M}_{}\mathrm{yr}^1\mathrm{Mpc}^3`$. In the instantaneous recycling approximation, the total stellar mass density observed today is

$$\rho _{}(t_H)=(1R)_0^{t_H}\dot{\rho }_{}(t)𝑑t5\times 10^8I_{40}\mathrm{M}_{}\mathrm{Mpc}^3,$$
(5)   
(corresponding to $`\mathrm{\Omega }_{}=0.007I_{40}`$) where $`R`$ is the mass fraction of a generation of stars that is returned to the interstellar medium, $`R0.3`$ for a Salpeter IMF. The optical/COBE background therefore requires that about 10% of the nucleosynthetic baryons ($`\mathrm{\Omega }_bh_{50}=0.08`$ ) are in the forms of stars and their remnants. The predicted stellar mass-to-blue light ratio is $`M/L_B5`$. Note that these values are quite sensitive to the lower-mass cutoff of the IMF, as very-low mass stars can contribute significantly to the mass but not to the integrated light of the whole stellar population. A lower cutoff of 0.5$`\mathrm{M}_{}`$ instead of the 0.1$`\mathrm{M}_{}`$ adopted would decrease the mass-to-light ratio (and $`\mathrm{\Omega }_{}`$) by a factor of 1.9 for a Salpeter function.

## 5 A constant star-formation density?

Based on the agreement between the $`z3`$ and $`z4`$ luminosity functions at the bright end, it has been recently argued by  that the decline in the luminosity density of faint HDF Lyman-break galaxies observed in the same redshift interval  may not be real, but simply due to sample variance in the HDF. When extinction corrections are applied, the emissivity per unit comoving volume due to star formation may then remain essentially flat for all redshift $`z\mathrm{}>1`$ (see Fig. 2). While this has obvious implications for hierarchical models of structure formation, the epoch of first light, and the reionization of the intergalactic medium (IGM), it is also interesting to speculate on the possibility of a constant star-formation density at all epochs $`0z5`$, as recently advocated by . Figure 3 (left panel) shows the time evolution of the near-IR rest-frame luminosity density of a stellar population characterized by a Salpeter IMF, solar metallicity, and a (constant) star-formation rate of $`\dot{\rho _{}}=0.054\mathrm{M}_{}\mathrm{yr}^1\mathrm{Mpc}^3`$ (needed to produce the observed EBL). <sup>1</sup><sup>1</sup>1The near-IR light is dominated by near-solar mass evolved stars, the progenitors of which make up the bulk of a galaxy’s stellar mass, and is sensitive to the past star-formation history. The predicted evolution appears to be a poor match to the observations: it overpredicts the local $`K`$-band luminosity density  and undepredicts the 1$`\mu `$m emissivity at $`z1`$ from the CFRS survey .

## 6  A population of hidden AGNs?

Recent dynamical evidence indicates that supermassive black holes reside at the center of most nearby galaxies. The available data (about 30 objects) show a strong correlation (but with a large scatter) between bulge and black hole mass , with $`M_{\mathrm{bh}}=0.006M_{\mathrm{bulge}}`$ as a best-fit. The total mass density in spheroids today is $`\mathrm{\Omega }_{\mathrm{bulge}}=0.0036_{0.0017}^{+0.0024}`$ , implying a mean mass density of dead quasars

$$\rho _{\mathrm{bh}}=1.34_{0.6}^{+0.9}\times 10^6\mathrm{M}_{}\mathrm{Mpc}^3.$$
(6)   
Noting that the observed energy density from all quasars is equal to the emitted energy divided by the average quasar redshift , the total contribution to the EBL from accretion onto black holes is

$$I_{\mathrm{bh}}=\frac{c^3}{4\pi }\frac{\eta \rho _{\mathrm{bh}}}{1+z}18\mathrm{nW}\mathrm{m}^2\mathrm{sr}^1\eta _{0.1}1+z^1,$$
(7)   
where $`\eta _{0.1}`$ is the efficiency for transforming accreted rest-mass energy into radiation (in units of 10%). A population of AGNs at (say) $`z1`$ could then make a significant contribution to the FIR background if dust-obscured accretion onto supermassive black holes is an efficient process .

## 7 Reionization of the IGM

The history of the transition from a neutral universe to one that is almost fully ionized can reveal the character of cosmological ionizing sources and constrain the star formation activity at high redshifts. The existence of a filamentary, low-density intergalactic medium (IGM), which contains the bulk of the hydrogen and helium in the universe, is predicted as a product of primordial nucleosynthesis  and of hierarchical models of gravitational instability with “cold dark matter” (CDM) , . The application of the Gunn-Peterson constraint on the amount of smoothly distributed neutral material along the line of sight to distant objects requires the hydrogen component of the diffuse IGM to have been highly ionized by $`z5`$ , and the helium component by $`z2.5`$ . From QSO absorption studies we also know that neutral hydrogen accounts for only a small fraction, $`10\%`$, of the nucleosynthetic baryons at early epochs . It thus appears that substantial sources of ultraviolet photons were present at $`z\mathrm{}>5`$, perhaps low-luminosity quasars  or a first generation of stars in virialized dark matter halos with $`T_{\mathrm{vir}}10^410^5`$K , , .

The existence of a decline in the space density of bright quasars at redshifts beyond $`3`$ was first suggested by , and has been since then the subject of a long-standing debate. In recent years, several optical surveys have consistently provided new evidence for a turnover in the QSO counts , , , . The interpretation of the drop-off observed in optically selected samples is equivocal, however, because of the possible bias introduced by dust obscuration arising from intervening systems. Radio emission, on the other hand, is unaffected by dust, and it has recently been shown  that the space density of radio-loud quasars also decreases strongly for $`z>3`$. This argues that the turnover is indeed real and that dust along the line of sight has a minimal effect on optically-selected QSOs (Figure 4, left panel). The QSO emission rate of hydrogen ionizing photons per unit comoving volume is shown in Figure 4 (right panel) . It is important to notice that the procedure adopted to derive this quantity implies a large correction for incompleteness at high-$`z`$. With a fit to the quasar luminosity function (LF) which goes as $`\varphi (L)L^{1.64}`$ at the faint end , the contribution to the emissivity converges rather slowly, as $`L^{0.36}`$. At $`z=4`$, for example, the blue magnitude at the break of the LF is $`M_{}25.4`$, comparable or slightly fainter than the limits of current high-$`z`$ QSO surveys. A large fraction, about 90% at $`z=4`$ and even higher at earlier epochs, of the ionizing quasar emissivity is therefore produced by sources that have not been actually observed, and are assumed to be present based on an extrapolation from lower redshifts.

Galaxies with ongoing star-formation are another obvious source of Lyman continuum photons. Since the rest-frame UV continuum at 1500 Å (redshifted into the visible band for a source at $`z3`$) is dominated by the same short-lived, massive stars which are responsible for the emission of photons shortward of the Lyman edge, the needed conversion factor, about one ionizing photon every 10 photons at 1500 Å, is fairly insensitive to the assumed IMF and is independent of the galaxy history for $`t10^7`$ yr. Figure 4 shows the estimated Lyman-continuum luminosity density of galaxies at $`z3`$.<sup>2</sup><sup>2</sup>2At all ages $`\mathrm{}>0.1`$ Gyr one has $`L(1500)/L(912)6`$ for a Salpeter mass function and constant SFR . This number neglects any correction for intrinsic $`\mathrm{I}`$absorption. The data point assumes a value of $`f_{\mathrm{esc}}=0.5`$ for the unknown fraction of ionizing photons which escapes the galaxy $`\mathrm{I}`$layers into the intergalactic medium. A substantial population of dwarf galaxies below the detection threshold, i.e. having star-formation rates $`<0.3\mathrm{M}_{}\mathrm{yr}^1`$, and with a space density in excess of that predicted by extrapolating to faint magnitudes the $`\alpha =1.38`$ best-fit Schechter function, may be expected to form at early times in hierarchical clustering models, and has been recently proposed by  and  as a possible candidate for photoionizing the IGM at these epochs. One should note that, while highly reddened galaxies at high redshifts would be missed by the dropout color technique (which isolates sources that have blue colors in the optical and a sharp drop in the rest-frame UV), it seems unlikely that very dusty objects (with $`f_{\mathrm{esc}}1`$) would contribute in any significant manner to the ionizing metagalactic flux.

As the hydrogen mean recombination timescale, $`\overline{t}_{\mathrm{rec}}`$, at high redshifts is much smaller than the then Hubble time , it is possible to compute at any given epoch a critical value for the photon emission rate per unit cosmological comoving volume,

$$\dot{𝒩}_{\mathrm{ion}}(z)=\frac{\overline{n}_\mathrm{H}(0)}{\overline{t}_{\mathrm{rec}}(z)}=(10^{51.2}\mathrm{s}^1\mathrm{Mpc}^3)C_{30}\left(\frac{1+z}{6}\right)^3\left(\frac{\mathrm{\Omega }_bh_{50}^2}{0.08}\right)^2,$$
(8)   
independently of the (unknown) previous emission history of the universe: only rates above this value will provide enough UV photons to ionize the IGM by that epoch. Here $`\overline{n}_\mathrm{H}(0)`$ is the mean hudrogen density of the expanding IGM at the present-epoch, and $`C`$ is the ionized hydrogen clumping factor. One can then compare our determinations of $`\dot{𝒩}_{\mathrm{ion}}`$ to the estimated contribution from QSOs and star-forming galaxies. The uncertainty on this critical rate is difficult to estimate, as it depends on the clumpiness of the IGM (scaled in the expression above to the value inferred at $`z=5`$ from numerical simulations ) and the nucleosynthesis constrained baryon density. The evolution of the critical rate as a function of redshift is plotted in Figure 4. While $`\dot{𝒩}_{\mathrm{ion}}`$ is comparable to the quasar contribution at $`z\mathrm{}>3`$, there is some indication of a deficit of Lyman continuum photons at $`z=5`$. For bright, massive galaxies to produce enough UV radiation at $`z=5`$, their space density would have to be comparable to the one observed at $`z3`$, with most ionizing photons being able to escape freely from the regions of star formation into the IGM. This scenario may be in conflict with direct observations of local starbursts below the Lyman limit showing that at most a few percent of the stellar ionizing radiation produced by these luminous sources actually escapes into the IGM .<sup>3</sup><sup>3</sup>3Note that, at $`z=3`$, Lyman-break galaxies would radiate more ionizing photons than QSOs for $`f_{\mathrm{esc}}\mathrm{}>30\%`$. If, on the other hand, faint QSOs with (say) $`M_{\mathrm{AB}}=19`$ at rest-frame ultraviolet frequencies were to provide all the required ionizing flux, their comoving space density would be such ($`0.0015`$Mpc<sup>-3</sup>) that about 50 of them would expected in the HDF down to $`I_{\mathrm{AB}}=27.2`$. At $`z\mathrm{}>5`$, they would appear very red in $`VI`$ as the Ly$`\alpha `$forest is shifted into the visible. This simple model can be ruled out, however, as there is only a handful (7) of sources in the HDF with $`(VI)_{\mathrm{AB}}>1.5`$ mag down to this magnitude limit.

It is interesting to convert the derived value of $`\dot{𝒩}_{\mathrm{ion}}`$ into a “minimum” SFR per unit (comoving) volume, $`\dot{\rho }_{}`$ (hereafter we assume $`\mathrm{\Omega }_bh_{50}^2=0.08`$ and $`C=30`$):

$$\dot{\rho }_{}(z)=\dot{𝒩}_{\mathrm{ion}}(z)\times 10^{53.1}f_{\mathrm{esc}}^10.013f_{\mathrm{esc}}^1\left(\frac{1+z}{6}\right)^3\mathrm{M}_{}\mathrm{yr}^1\mathrm{Mpc}^3.$$
(9)   
(The conversion factor assumes a Salpeter IMF with solar metallicity). The star-formation density given in equation (9) is comparable with the value directly “observed” (i.e., uncorrected for dust reddening) at $`z3`$ .