Rates of tidal disruption of stars by massive central black holes

There is strong evidence for some kind of massive dark object in the centres of many galaxy bulges. The detection of flares from tidally disrupted stars could confirm that these objects are black holes (BHs). Here we present calculations of the stellar disruption rates in detailed dynamical models of real galaxies, taking into account the refilling of the loss cone of stars on disruptable orbits by two-body relaxation and tidal forces in non-spherical galaxies. The highest disruption rates (one star per 10⁴ yr) occur in faint ( L ≲10¹⁰ L_⊙) galaxies, which have steep central density cusps. More luminous galaxies are less dense and have much longer relaxation times and more massive BHs. Dwarf stars in such galaxies are swallowed whole by the BH and hence do not emit flares; giant stars could produce flares as often as every 10⁵ yr, although the rate depends sensitively on the shape of the stellar distribution function. We discuss the possibility of detecting disruption flares in current supernova searches. The total mass of stars consumed over the lifetime of the galaxy is of the order of 10⁶ M_⊙, independent of galaxy luminosity; thus, disrupted stars may contribute significantly to the present BH mass in galaxies fainter than ∼10⁹ L_⊙.     

The dynamical evolution of stellar superclusters

Recent images taken with the Hubble Space Telescope ( HST ) of the interacting disc galaxies NGC 4038/4039 (the Antennae) reveal clusters of many dozens and possibly hundreds of young compact massive star clusters within projected regions spanning about 100 to 500 pc. It is shown here that a large fraction of the individual star clusters merge within a few tens to a hundred Myr. Bound stellar systems with radii of a few hundred parsecs, masses ≲ 10⁹ M⊙ and relaxation times of 10¹¹ − 10¹² yr may form from these. These spheroidal dwarf galaxies contain old stars from the pre-merger galaxy and much younger stars formed in the massive star clusters, and possibly from later gas accretion events. The possibility that star formation in the outer regions of gas-rich tidal tails may also lead to superclusters is raised. The mass-to-light ratio of these objects is small, because they contain an insignificant amount of dark matter. After many hundred Myr such systems may resemble dwarf spheroidal satellite galaxies with large apparent mass-to-light ratios, if tidal shaping is important.     

Galactic winds and circulation of the interstellar medium in dwarf galaxies

We study, through 2D hydrodynamical simulations, the feedback of a starburst on the ISM of typical gas-rich dwarf galaxies. The main goal is to address the circulation of the ISM and metals following the starburst. We assume a single-phase rotating ISM in equilibrium in the galactic potential generated by a stellar disc and a spherical dark halo. The starburst is assumed to occur in a small volume in the centre of the galaxy, and it generates a mechanical power of 3.8×10³⁹ or 3.8×10⁴⁰ erg s⁻¹ for 30 Myr. We find, in accordance with previous investigations, that the galactic wind is not very effective in removing the ISM. The metal-rich stellar ejecta, however, can be efficiently expelled from the galaxy and dispersed in the intergalactic medium.
Moreover, we find that the central region of the galaxy is always replenished with cold and dense gas a few 100 million years after the starburst, achieving the requisite for a new star formation event in ≈0.5–1 Gyr. The hydrodynamical evolution of galactic winds is thus consistent with the episodic star formation regime suggested by many chemical evolution studies.
We also discuss the X-ray emission of these galaxies and find that the observable (emission-averaged) abundance of the hot gas underestimates the real one if thermal conduction is effective. This could explain the very low hot-gas metallicities estimated in starburst galaxies.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

On the fate of processed matter in dwarf galaxies

Two-dimensional calculations of the evolution of remnants generated by the strong mechanical energy deposited by stellar clusters in dwarf galaxies ( M  ∼ 10⁹–10¹⁰ M⊙) are presented. The evolution is followed for times longer than both the blowout time and the presumed span of energy injection generated by a coeval massive stellar cluster. The remnants are shown to end up wrapping around the central region of the host galaxy, while growing to kiloparsec-scale dimensions. Properties of the remnants such as luminosity, size, swept-up mass and expansion speed are given as functions of time for all calculated cases.   The final fate of the swept-up galactic gas and of the matter processed by the central starburst is shown to be highly dependent on the properties of the low-density galactic halo. Superbubbles powered by star clusters, with properties similar to those inferred from the observations, slow down in the presence of an extended halo to expansion speeds smaller than the host galaxy escape velocity. Values of the critical luminosity required for the superbubbles to reach the edge of the galaxies with a speed comparable to the escape speed are derived analytically and numerically. The critical luminosities are larger than those in the detected sources, and thus the superbubbles in amorphous dwarf galaxies must have already undergone blowout and are presently evolving into an extended low-density halo. This will inhibit the loss of the swept-up and processed matter from the galaxy.     

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.     

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.     

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.     

Evolution of massive binary black holes

Since many or most galaxies have central massive black holes (BHs), mergers of galaxies can form massive binary black holes (BBHs). In this paper we study the evolution of massive BBHs in realistic galaxy models, using a generalization of techniques used to study tidal disruption rates around massive BHs. The evolution of BBHs depends on BH mass ratio and host galaxy type. BBHs with very low mass ratios (say, ≲0.001) are hardly ever formed by mergers of galaxies, because the dynamical friction time-scale is too long for the smaller BH to sink into the galactic centre within a Hubble time. BBHs with moderate mass ratios are most likely to form and survive in spherical or nearly spherical galaxies and in high-luminosity or high-dispersion galaxies; they are most likely to have merged in low-dispersion galaxies (line-of-sight velocity dispersion ≲90 km s⁻¹) or in highly flattened or triaxial galaxies.
The semimajor axes and orbital periods of surviving BBHs are generally in the range  10^-3–10 pc  and  10–10⁵ yr;  they are also larger in high-dispersion galaxies than in low-dispersion galaxies, larger in nearly spherical galaxies than in highly flattened or triaxial galaxies, and larger for BBHs with equal masses than for BBHs with unequal masses. The orbital velocities of surviving BBHs are generally in the range  10²–10⁴ km s^-1  . The methods of detecting surviving BBHs are also discussed.
If no evidence of BBHs is found in AGNs, this may be either because gas plays a major role in BBH orbital decay or because nuclear activity switches on soon after a galaxy merger, and ends before the smaller BH has had time to spiral to the centre of the galaxy.     

K-band imaging of strong Ca ii-absorber host galaxies at z∼1

We present K -band imaging of fields around 30 strong Ca  ii absorption-line systems, at  0.7 < z < 1.2  , three of which are confirmed damped Lyman α systems. A significant excess of galaxies is found within 6.0 arcsec (≃50 kpc) from the absorber line of sight. The excess galaxies are preferentially luminous compared to the population of field galaxies. A model in which field galaxies possess a luminosity-dependent cross-section for Ca  ii absorption of the form  ( L / L *)^0.7  reproduces the observations well. The luminosity-dependent cross-section for the Ca  ii absorbers appears to be significantly stronger than the established  ( L / L *)^0.4  dependence for Mg  ii absorbers. The associated galaxies lie at large physical distances from the Ca  ii -absorbing gas; we find a mean impact parameter of 24 kpc  ( H ₀= 70 km s⁻¹ Mpc⁻¹)  . Combined with the observed number density of Ca  ii absorbers the large physical separations result in an inferred filling factor of only ∼10 per cent. The physical origin of the strong Ca  ii absorption remains unclear, possible explanations vary from very extended discs of the luminous galaxies to associated dwarf galaxy neighbours, remnants of outflows from the luminous galaxies, or tidal debris from cannibalism of smaller galaxies.     

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.     

Density–potential pairs for spherical stellar systems with Sérsic light profiles and (optional) power-law cores

Popular models for describing the luminosity-density profiles of dynamically hot stellar systems (e.g. Jaffe, Hernquist, Dehnen) were constructed with the desire to match the deprojected form of an   R ^1/4  light profile. Real galaxies, however, are now known to have a range of different light-profile shapes that scale with mass. Consequently, although highly useful, the above models have implicit limitations, and this is illustrated here through their application to a number of real galaxy density profiles. On the other hand, the analytical density profile given by Prugniel & Simien closely matches the deprojected form of Sérsic   R ^{1/ n}   light profiles – including deprojected exponential light profiles. It is thus applicable for describing bulges in spiral galaxies, dwarf elliptical galaxies, and both ordinary and giant elliptical galaxies. Moreover, the observed Sérsic quantities define the parameters of the density model. Here we provide simple equations, in terms of elementary and special functions, for the gravitational potential and force associated with this density profile. Furthermore, to match galaxies with partially depleted cores, and better explore the supermassive black hole/galaxy connection, we have added a power-law core to this density profile and derived similar expressions for the potential and force of this hybrid profile. Expressions for the mass and velocity dispersion, assuming isotropy, are also given. These spherical models may also prove appropriate for describing the dark matter distribution in haloes formed from ΛCDM cosmological simulations.     

The CO content of the Local Group dwarf irregular galaxies IC 5152, UGCA 438 and the Phoenix dwarf

We present a search for  CO(1 → 0)  emission in three Local Group dwarf irregular galaxies: IC 5152, the Phoenix dwarf and UGCA 438, using the ATNF Mopra radio telescope. Our scans largely cover the optical extent of the galaxies and the stripped H  i cloud west of the Phoenix dwarf. Apart from a tentative but non-significant emission peak at one position in the Phoenix dwarf, no significant emission was detected in the CO spectra of these galaxies. For a velocity width of 6 km s⁻¹, we derive 4σ upper limits of 0.03, 0.04 and 0.06 K km s⁻¹ for IC 5152, the Phoenix dwarf and UGCA 438, respectively. This is an improvement of over a factor of 10 compared with previous observations of IC 5152; the other two galaxies had not yet been observed at millimetre wavelengths. Assuming a Galactic CO-to-H₂ conversion factor, we derive upper limits on the molecular gas mass of  6.2 × 10⁴, 3.7 × 10³  and  1.4 × 10⁵ M_⊙  for IC 5152, the Phoenix dwarf and UGCA 438, respectively. We investigate two possible causes for the lack of CO emission in these galaxies. On the one hand, there may be a genuine lack of molecular gas in these systems, in spite of the presence of large amounts of neutral gas. However, in the case of IC 5152 which is actively forming stars, molecular gas is at least expected to be present in the star-forming regions. On the other hand, there may be a large increase in the CO-to-H₂ conversion factor in very low-metallicity dwarfs  (−2 ≤[Fe/H]≤−1)  , making CO a poor tracer of the molecular gas content in dwarf galaxies.     

Diffuse X-ray emission from late-type galaxy haloes

Current theories of galaxy formation predict that spiral galaxies are embedded in a reservoir of hot gas. This gas is able to cool on to the galaxy, replenishing cold gas that is consumed by star formation. Estimates of the X-ray luminosity emitted in the cooling region suggest a bolometric luminosity of the order of 10×10⁴¹ erg s⁻¹ in massive systems. We have used ROSAT PSPC data to search for extended X-ray emission from the haloes of three nearby, massive, late-type galaxies: NGC 2841, 4594 and 5529. We infer 95 per cent upper limits on the bolometric X-ray luminosities of the haloes of NGC 2841, 4594 and 5529 of 0.4, 1.2 and 3.8×10⁴¹ erg s⁻¹ respectively. Thus, the true luminosity lies well below the straightforward theoretical prediction. We discuss this discrepancy and suggest a number of ways in which the theoretical model might be brought into agreement with the observational results. A possible solution is that the gravitational potentials of the dark matter haloes of these galaxies are weaker than assumed in the current model. Alternatively, the present-day accretion may be substantially less than is required on average to build the disc over the Hubble time. Our results are, however, based on only three galaxies, none of which is ideal for this kind of study. A larger data set is required to explore this important problem further.     

The clustering of Lyman-break galaxies

We calculate the statistical clustering of Lyman-break galaxies predicted in a selection of currently fashionable structure formation scenarios. These models are all based on the cold dark matter model, but vary in the amount of dark matter, the initial perturbation spectrum, the background cosmology and the presence or absence of a cosmological constant term. If Lyman-break galaxies form as a result of hierarchical merging, the amplitude of clustering depends quite sensitively on the minimum halo mass that can host such a galaxy. Interpretation of the recent observations by Giavalisco et al. would therefore be considerably clarified by a direct determination of the relevant halo properties. For a typical halo mass around 10¹¹  h ⁻¹ M⊙ the observations do not discriminate strongly between cosmological models, but if the appropriate mass is larger, say 10¹²  h ⁻¹ M⊙ (which seems likely on theoretical grounds), then the data strongly favour models with a low matter density.     

Unresolved emission and ionized gas in the bulge of M31

We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 10²⁷ erg s⁻¹ L⁻¹_⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 10⁶ M_⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L _X/SFR  (star formation rate) ratio equals  ∼9 × 10³⁸ (erg s⁻¹)(M_⊙ yr⁻¹)⁻¹  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.     

The anomalous properties of Markarian 1460

We present and discuss optical, near-infrared and H  i measurements of the galaxy Markarian 1460 at a distance of 19 Mpc in the Ursa Major Cluster. This low-luminosity ( M _B =−14) galaxy is unusual because (i) it is blue ( B − R =0.8) and has the spectrum of an H  ii galaxy, (ii) it has a light profile that is smooth and well fitted by an r ^1/4 and not an exponential function at all radii larger than the seeing, and (iii) it has an observed central brightness of about μ _B =20 mag arcsec⁻² , intermediate between those of elliptical galaxies (on the bright μ _B side) and normal low-luminosity dwarf irregular (on the low μ _B side) galaxies. No other known galaxy exhibits all these properties in conjunction. On morphological grounds this galaxy looks like a normal distant luminous elliptical galaxy, since the Fundamental Plane tells us that higher luminosity normal elliptical galaxies tend to have lower surface-brightnesses. Markarian 1460 has 2×10⁷ M_⊙ of H  i and a ratio M (H  i )/ L _B of 0.2, which is low compared to the typical values for star-forming dwarf galaxies. From the high surface-brightness and r ^1/4 profile, we infer that the baryonic component of Markarian 1460 has become self-gravitating through dissipative processes. From the colours, radio continuum, H  i and optical emission line properties, and yet smooth texture, we infer that Markarian 1460 has had significant star formation as recently as ∼1 Gyr ago but not today.     

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.     

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.