Occurrence of metal-free galaxies in the early Universe

The character of the first galaxies at redshifts z ≳ 10 strongly depends on their level of pre-enrichment, which is in turn determined by the rate of primordial star formation prior to their assembly. In order for the first galaxies to remain metal-free, star formation in minihaloes must be highly suppressed, most likely by H₂-dissociating Lyman–Werner (LW) radiation. We show that the build-up of such a strong LW background is hindered by two effects. First, the level of the LW background is self-regulated, being produced by the Population III (Pop III) star formation which it, in turn, suppresses. Secondly, the high opacity to LW photons which is built up in the relic H  ii regions left by the first stars acts to diminish the global LW background. Accounting for a self-regulated LW background, we estimate a lower limit for the rate of Pop III star formation in minihaloes at z ≳ 15. Further, we simulate the formation of a 'first galaxy' with virial temperature   T _vir≳ 10⁴ K  and total mass  ≳10⁸ M_⊙  at z ≳ 10, and find that complete suppression of previous Pop III star formation is unlikely, with stars of  ≳100 M_⊙  (Pop III.1) and  ≳10 M_⊙  (Pop III.2) likely forming. Finally, we discuss the implications of these results for the nature of the first galaxies, which may be observed by future missions such as the James Webb Space Telescope .     

Two populations of metal-free stars in the early Universe

We construct star formation histories at redshifts z ≳ 5 for two physically distinct populations of primordial, metal-free stars, motivated by theoretical and observational arguments that have hinted towards the existence of an intermediate stellar generation between Population III and Population I/II. Taking into account the cosmological parameters as recently revised by the Wilkinson Microwave Anisotropy Probe after three years of operation, we determine self-consistent reionization histories and discuss the resulting chemical enrichment from these early stellar generations. We find that the bulk of ionizing photons and heavy elements produced at high redshifts must have originated in Population II.5 stars, which formed out of primordial gas in haloes with virial temperatures ≳10⁴ K, and had typical masses ≳10 M_⊙. Classical Population III stars, formed in minihaloes and having masses ≳100 M_⊙, on the other hand, had only a minor impact on reionization and early metal enrichment. Specifically, we conclude that only ≃10 per cent by mass of metal-free star formation went into Population III.     

Impact of cosmic rays on Population III star formation

We explore the implications of a possible cosmic-ray (CR) background generated during the first supernova explosions that end the brief lives of massive Population III stars. We show that such a CR background could have significantly influenced the cooling and collapse of primordial gas clouds in minihaloes around redshifts of   z ∼ 15–20  , provided the CR flux was sufficient to yield an ionization rate greater than about 10⁻¹⁹ s⁻¹ near the centre of the minihalo. The presence of CRs with energies  ≲10⁷  eV would indirectly enhance the molecular cooling in these regions, and we estimate that the resulting lower temperatures in these minihaloes would yield a characteristic stellar mass as low as  ∼10 M_⊙  . CRs have a less-pronounced effect on the cooling and collapse of primordial gas clouds inside more massive dark matter haloes with virial masses  ≳10⁸ M_⊙  at the later stages of cosmological structure formation around   z ∼ 10–15  . In these clouds, even without CR flux the molecular abundance is already sufficient to allow cooling to the floor set by the temperature of the cosmic microwave background.     

Hyperluminous infrared galaxies

39 galaxies are now known, from follow-up of faint IRAS sources and from submillimetre observations of high-redshift AGN, with far-infrared luminosities >10¹³ L_⊙. 13 of these, which have been found in 60- or 850-μm surveys, form an important unbiased subsample. 12 have been found by comparison of 60-μm surveys with quasar or radio galaxy catalogues, or from infrared surveys with colour selection biased towards AGN, while a further 14 have been found through submillimetre observations of known high-redshift AGN. In this paper I argue, on the basis of detailed modelling of the spectral energy distributions of hyperluminous galaxies with accurate radiative transfer models, and from evidence of high gas mass in several cases, that the bulk of the emission from these galaxies at rest frame wavelengths ≥50 μm is caused by star formation. Even after correction for the effects of lensing, hyperluminous galaxies with emission peaking at rest frame wavelengths ≥50 μm are therefore undergoing star formation at rates >10³ M_⊙ yr⁻¹ and are strong candidates for being primeval galaxies, in the process of a major episode of star formation.     

Simulating the formation of a protocluster at z∼ 2

We present results from two high-resolution hydrodynamical simulations of protocluster regions at   z ≃ 2.1  . The simulations have been compared to observational results for the so-called Spiderweb galaxy system, the core of a putative protocluster region at   z = 2.16  , found around a radio galaxy. The simulated regions have been chosen so as to form a poor cluster with   M ₂₀₀≃ 10¹⁴  h ⁻¹ M_⊙  (C1) and a rich cluster with   M ₂₀₀≃ 2 × 10¹⁵  h ⁻¹ M_⊙  (C2) at   z = 0  . The simulated protoclusters show evidence of ongoing assembly of a dominating central galaxy. The stellar mass of the brightest cluster galaxy of the C2 system is in excess with respect to observational estimates for the Spiderweb galaxy, with a total star formation rate which is also larger than indicated by observations. We find that the projected velocities of galaxies in the C2 cluster are consistent with observations, while those measured for the poorer cluster C1 are too low compared with the observed velocities. We argue that the Spiderweb complex resembles the high-redshift progenitor of a rich galaxy cluster. Our results indicate that the included supernovae feedback is not enough to suppress star formation in these systems, supporting the need of introducing active galactic nuclei feedback. According to our simulations, a diffuse atmosphere of hot gas in hydrostatic equilibrium should already be present at this redshift, and enriched at a level comparable to that of nearby galaxy clusters. The presence of this gas should be detectable with future deep X-ray observations.     

The changing interstellar medium of massive elliptical galaxies and cosmic evolution of radio galaxies and quasars

The recently discovered apparent dramatic expansion in the effective radii of massive elliptical galaxies from   z ≃ 2  to ≃0.1 has been interpreted in terms of either galaxy mergers or the rapid loss of cold gas due to active galactic nuclei (AGN) feedback. In examining the latter case, we have quantified the extent of the expansion, which is uncertain observationally, in terms of the star formation parameters and time of the expulsion of the cold gas. In either case, the large global decrease in stellar density should translate into a major drop in the interstellar medium density and pressure with cosmic epoch. These cosmological changes are expected to have a major influence on the gas accretion mode, which will shift from 'cold' thin disc accretion at high redshifts towards 'hot' Bondi fed Advection Dominated Accretion Flow (ADAF) accretion at low redshifts. The decline of angular momentum inflow would then lead to a spin down of the black hole, for which we have calculated more precise time-scales; a value of about 0.2 Gyr is typical for a  10⁹ M_⊙  central black hole. These results have implications for the different cosmological evolutionary patterns found for the luminosity functions of powerful and weak radio galaxies.     

Galaxy groups in the 2dF Galaxy Redshift Survey: effects of environment on star formation

We estimate the fraction of star-forming galaxies in a catalogue of groups, constructed from the 2dF Galaxy Redshift Survey by Merchán & Zandivarez. We use the η spectral type parameter of galaxies and subdivide the sample of galaxies in groups into four types depending on the values of the η parameter following Madgwick et al. We obtain a strong correlation between the relative fraction of galaxies with high star formation and the parent group virial mass. We find that even in the environment of groups with low virial mass   M ∼10¹³ M_⊙  the star formation of their member galaxies is significantly suppressed. The relation between the fraction of early-type galaxies and the group virial mass obeys a simple power law spanning over three orders of magnitude in virial mass. Our results show quantitatively the way that the presence of galaxies with high star formation rates is inhibited in massive galaxy systems.     

The cooling of shock-compressed primordial gas

We find that at redshifts   z ≳ 10, HD  line cooling allows strongly shocked primordial gas to cool to the temperature of the cosmic microwave background (CMB). This temperature is the minimum value attainable via radiative cooling. Provided that the abundance of HD, normalized to the total number density, exceeds a critical level of  ∼10⁻⁸  , the CMB temperature floor is reached in a time which is short in comparison to the Hubble time. We estimate the characteristic masses of stars formed out of shocked primordial gas in the wake of the first supernovae, and resulting from the virialization of dark matter haloes during hierarchical structure formation to be  ∼10 M_⊙  . In addition, we show that cooling by HD enables the primordial gas in relic H  ii regions to cool to temperatures considerably lower than those reached via H₂ cooling alone. We confirm that HD cooling is unimportant in cases where the primordial gas does not go through an ionized phase, as in the formation process of the very first stars in   z ≳ 20  minihaloes of mass  ∼10⁶ M_⊙  .     

The local star formation rate and radio luminosity density

We present a new determination of the local volume-averaged star formation rate from the 1.4-GHz luminosity function of star forming galaxies. Our sample, taken from the   B ≤12  Revised Shapley–Ames catalogue (231 normal spiral galaxies over an effective area of 7.1 sr) has ≃100 per cent complete radio detections and is insensitive to dust obscuration and cirrus contamination. After removal of known active galaxies, the best-fitting Schechter function has a faint-end slope of  −1.27±0.07  in agreement with the local H α luminosity function, characteristic luminosity   L ∗=(2.6±0.7)×10²² W Hz⁻¹  and density   φ ∗=(4.8±1.1)×10⁻⁴ Mpc⁻³.  The inferred local radio luminosity density of  (1.73±0.37±0.03)×10¹⁹ W Hz⁻¹ Mpc⁻³  (Poisson noise, large-scale structure fluctuations) implies a volume-averaged star formation rate ∼2 times larger than the Gallego et al. H α estimate, i.e.   ρ _1.4 GHz=(2.10±0.45±0.04)×10⁻² M_⊙ yr⁻¹ Mpc⁻³  for a Salpeter initial mass function from  0.1–125 M_⊙  and Hubble constant of 50 km s⁻¹ Mpc⁻¹. We demonstrate that the Balmer decrement is a highly unreliable extinction estimator, and argue that optical–ultraviolet (UV) star formation rates (SFRs) are easily underestimated, particularly at high redshift.     

Boundary layer on the surface of a neutron star

In an attempt to model the accretion on to a neutron star in low-mass X-ray binaries, we present 2D hydrodynamical models of the gas flow in close vicinity of the stellar surface. First, we consider a gas pressure-dominated case, assuming that the star is non-rotating. For the stellar mass we take   M _star= 1.4 × 10⁻² M_⊙  and for the gas temperature   T = 5 × 10⁶ K  . Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of T ≃ 10⁸ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.     

The growth of central and satellite galaxies in cosmological smoothed particle hydrodynamics simulations

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to  7 × 10⁹ M_⊙  . Most friends-of-friends haloes in the simulation have a distinct central galaxy, typically 2–5 times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive haloes form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5–1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Typically, satellites in our simulation are 0.1–0.2 mag bluer than in models that assume no gas accretion on to satellites after a halo merger. Since   z = 1  , 27 per cent of central galaxies (above  3 × 10¹⁰ M_⊙  ) and 22 per cent of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain 'central' objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.     

Baryonic conversion tree: the global assembly of stars and dark matter in galaxies from the Sloan Digital Sky Survey

Using the spectroscopic sample of the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed approximately 80 per cent of their stars at   z > 1  , while for the less massive ones the value is only approximately 20 per cent. By comparing the BCT with the dark matter merger tree, we find indications that star formation efficiency at   z > 1  had to be approximately a factor of two higher than today (∼10 per cent) in galaxies with present-day stellar mass larger than  2 × 10¹¹ M_⊙  , if this early star formation occurred in the main progenitor. Therefore, the λ cold dark matter (LCDM) paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than  10¹¹ M_⊙  , efficient star formation seems to have been triggered at   z ∼ 0.2  . We show that there is a characteristic mass  ( M _*∼ 10¹⁰ M_⊙)  for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical models of galaxy formation.     

Numerical simulations of hot halo gas in galaxy mergers

Galaxy merger simulations have explored the behaviour of gas within the galactic disc, yet the dynamics of hot gas within the galaxy halo have been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with metal-free hot halo gas. To isolate the effect of the halo gas, we simulate only the dark matter halo and the hot halo gas over a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that (i) a strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to   T ∼ 10^6.3 K  . (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is  ≳10³⁹ erg s⁻¹  for halo gas fractions larger than 10 per cent. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as  10⁴² erg s⁻¹  . This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) ∼10–20 per cent of the initial gas mass is unbound from the galaxies for equal-mass mergers, while 3–5 per cent of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the intergalactic medium with metals.     

Dry mergers: a crucial test for galaxy formation

We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a nearly constant dry merger rate of  ∼6 × 10⁻⁵ Mpc⁻³ Gyr⁻¹  at   z ≤ 1  and a steep decline at larger z . Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early- and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of   M _{*, c} ∼ 6.3 × 10¹⁰ M_⊙  by mergers only. As a consequence of this transition, we find that around   M _{*, c}   , the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggests that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z . We however find that only 10–20 per cent of galaxies more massive than   M _{*, c}   experience dry major mergers within their last Gyr at any given redshift   z ≤ 1  .     

The SCUBA HAlf Degree Extragalactic Survey – VI. 350-μm mapping of submillimetre galaxies

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 μm has been carried out to map the regions around several 850-μm-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared (IR) luminosities and hence star formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350-μm flux densities for 24 850-μm sources, seven of which are detected at ≥2.5σ within a 10 arcsec search radius of the 850-μm positions. When results from the literature are included the total number of 350-μm flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-IR (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich  ( M _dust≃ 9 × 10⁸ M_⊙)  , luminous  ( L _FIR≃ 2 × 10¹² L_⊙)  star-forming galaxies with intrinsic dust temperatures of ≃35 K and star formation rates of  ≃400 M_⊙ yr⁻¹  . We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350-μm source counts, N ₃₅₀(>25 mJy) ∼ 200–500 deg⁻².     

The census of nuclear activity of late-type galaxies in the Virgo cluster

The first spectroscopic census of active galactic nuclei (AGNs) associated with late-type galaxies in the Virgo cluster was carried out by observing 213 out of a complete set of 237 galaxies more massive than   M _dyn > 10^8.5 M_⊙  . Among them, 77 are classified as AGNs [including 21 transition objects, 47 low-ionization nuclear emission regions (LINERs) and nine Seyferts] and comprise 32 per cent of the late-type galaxies in Virgo. Due to spectroscopic incompleteness, at most 21 AGNs are missed in the survey, so that the fraction would increase up to 41 per cent. Using corollary near-infrared observations that enable us to estimate galaxy dynamical masses, it is found that AGNs are hosted exclusively in massive galaxies, i.e.   M _dyn≳ 10¹⁰ M_⊙  . Their frequency increases steeply with the dynamical mass from zero at   M _dyn≈ 10^9.5 M_⊙  to virtually 1 at   M _dyn > 10^11.5 M_⊙  . These frequencies are consistent with those of low-luminosity AGNs found in the general field by the Sloan Digital Sky Survey. Massive galaxies that harbour AGNs commonly show conspicuous r -band star-like nuclear enhancements. Conversely, they often, but not necessarily, contain massive bulges. A few well-known AGNs (e.g. M61, M100, NGC 4535) are found in massive Sc galaxies with little or no bulge. The AGN fraction seems to be only marginally sensitive to galaxy environment. We infer the black hole masses using the known scaling relations of quiescent black holes. No black holes lighter than  ∼10⁶ M_⊙  are found active in our sample.     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

On the galaxy stellar mass function, the mass–metallicity relation and the implied baryonic mass function

A comparison between published field galaxy stellar mass functions (GSMFs) shows that the cosmic stellar mass density is in the range 4–8 per cent of the baryon density (assuming  Ω_b= 0.045  ). There remain significant sources of uncertainty for the dust correction and underlying stellar mass-to-light ratio even assuming a reasonable universal stellar initial mass function. We determine the   z < 0.05  GSMF using the New York University Value-Added Galaxy Catalog sample of 49 968 galaxies derived from the Sloan Digital Sky Survey and various estimates of stellar mass. The GSMF shows clear evidence for a low-mass upturn and is fitted with a double Schechter function that has  α₂≃−1.6  . At masses below  ∼10^8.5 M_⊙  , the GSMF may be significantly incomplete because of missing low-surface-brightness galaxies. One interpretation of the stellar mass–metallicity relation is that it is primarily caused by a lower fraction of available baryons converted to stars in low-mass galaxies. Using this principle, we determine a simple relationship between baryonic mass and stellar mass and present an 'implied baryonic mass function'. This function has a faint-end slope,  α₂≃−1.9  . Thus, we find evidence that the slope of the low-mass end of the galaxy mass function could plausibly be as steep as the halo mass function. We illustrate the relationship between halo baryonic mass function → galaxy baryonic mass function → GSMF. This demonstrates the requirement for peak galaxy formation efficiency at baryonic masses  ∼10¹¹ M_⊙  corresponding to a minimum in feedback effects. The baryonic-infall efficiency may have levelled off at lower masses.     

The cooling function of HD molecule revisited

We report new calculations of the cooling rate of primordial gas by the HD molecule, taking into account its ro-vibrational structure. The HD cooling function is calculated including radiative and collisional transitions for   J ≤ 8  rotational levels, and for the vibrational levels v = 0, 1, 2 and 3. The ro-vibrational level population is calculated from the balance equation assuming steady state. The cooling function is evaluated in the ranges of the kinetic temperatures, T _k, from 10² to  2 × 10⁴ K  and the number densities, n _H, from 1 to  10⁸ cm⁻³  . We find that the inclusion of collisional ro-vibrational transitions increases significantly the HD cooling efficiency, in particular for high densities and temperatures. For   n _H≳ 10⁵  and   T _k∼ 10⁴ K  the cooling function becomes more than an order of magnitude higher than previously reported. We give also the HD cooling rate in the presence of the cosmic microwave radiation field for radiation temperatures of 30, 85 and 276 K (redshifts of 10, 30 and 100). The tabulated cooling functions are available at http://www.cifus.uson.mx/Personal_Pages/anton/DATA/HD_cooling/HD_cool.html . We discuss the relevance to explore the effects of including our results into models and simulations of galaxy formation, especially in the regime when gas cools down from temperatures above ∼3000 K.     

Dust and hydrogen molecules in the extremely metal-poor dwarf galaxy SBS 0335–052

During the early stages of galaxy evolution, the metallicity is generally low and nearby metal-poor star-forming galaxies may provide templates for primordial star formation. In particular, the dust content of such objects is of great importance, because early molecular formation can take place on grains. To gain insight into primeval galaxies at high redshift, we examine the dust content of the nearby extremely low-metallicity galaxy SBS  0335–052  which hosts a very young starburst (≲10⁷ yr). In young galaxies, the dust formation rate in Type II supernovae governs the amount of dust, and by incorporating recent results on dust production in Type II supernovae we model the evolution of dust content. If the star-forming region is compact (≲100 pc), as suggested by observations of SBS  0335–052  , our models consistently explain the quantity of dust, far-infrared luminosity, and dust temperature in this low-metallicity object. We also discuss the H₂ abundance. The compactness of the region is important to H₂ formation, because the optical depth of dust for UV photons becomes large and H₂ dissociation is suppressed. We finally focus on implications for damped Ly α systems.