Populating a cluster of galaxies – I. Results at

We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star formation history and the orbital motion of all galaxies more luminous than the Fornax dwarf spheroidal, while dark halo structure is tracked consistently throughout the cluster for all galaxies more luminous than the SMC. Within its virial radius this model contains about     dark matter particles and almost 5000 distinct dynamically resolved galaxies. Simulations of this same cluster at a variety of resolutions allow us to check explicitly for numerical convergence both of the dark matter structures produced by our new parallel N -body and substructure identification codes, and of the galaxy populations produced by the phenomenological models we use to follow cooling, star formation, feedback and stellar aging. This baryonic modelling is tuned so that our simulations reproduce the observed properties of isolated spirals outside clusters. Without further parameter adjustment our simulations then produce a luminosity function, a mass-to-light ratio, luminosity, number and velocity dispersion profiles, and a morphology–radius relation which are similar to those observed in real clusters. In particular, since our simulations follow galaxy merging explicitly, we can demonstrate that it accounts quantitatively for the observed cluster population of bulges and elliptical galaxies.     

Globular cluster systems in nearby dwarf galaxies – I. HST/ACS observations and dynamical properties of globular clusters at low environmental density

We investigate the old globular cluster (GC) population of 68 faint  ( M _V > −16 mag)  dwarf galaxies located in the halo regions of nearby (≲12 Mpc) loose galaxy groups and in the field environment based on archival Hubble Space Telescope ( HST )/Advanced Camera for Surveys (ACS) images in F606W and F814W filters. The combined colour distribution of 175 GC candidates peaks at  ( V − I ) = 0.96 ± 0.07 mag  and the GC luminosity function turnover for the entire sample is found at   M _{V ,TO}=−7.6 ± 0.11 mag  , similar to the old metal-poor Large Magellanic Cloud (LMC) GC population. Our data reveal a tentative trend of   M _{V ,TO}  becoming fainter from late- to early-type galaxies. The luminosity and colour distributions of GCs in dIrrs show a lack of faint blue GCs (bGCs). Our analysis reveals that this might reflect a relatively younger GC system than typically found in luminous early-type galaxies. If verified by spectroscopy, this would suggest a later formation epoch of the first metal-poor star clusters in dwarf galaxies. We find several bright (massive) GCs which reside in the nuclear regions of their host galaxies. These nuclear clusters have similar luminosities and structural parameters as the peculiar Galactic clusters suspected of being the remnant nuclei of accreted dwarf galaxies, such as M54 and ωCen. Except for these nuclear clusters, the distribution of GCs in dIrrs in the half-light radius versus cluster mass plane is very similar to that of Galactic young halo clusters, which suggests comparable formation and dynamical evolution histories. A comparison with theoretical models of cluster disruption indicates that GCs in low-mass galaxies evolve dynamically as self-gravitating systems in a benign tidal environment.     

Faint Lyman-break galaxies as a crucial test for galaxy formation models

It has recently been shown that galaxy formation models within the Λ cold dark matter cosmology predict that, compared to the observed population, small galaxies (with stellar masses  <10¹¹ M_⊙  ) form too early, are too passive since   z ∼ 3  and host too old stellar populations at   z = 0  . We then expect an overproduction of small galaxies at   z ≳ 4  that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the morgana galaxy formation model and grasil spectrophotometric  +  radiative transfer code to generate mock catalogues of deep fields observed with Hubble Space Telescope Advanced Camera for Surveys. We add observational noise and the effect of Lyman α emission, and perform colour–colour selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i -dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of   z ₈₅₀∼ 27 V   -dropouts at   z ∼ 5  . We quantify the properties of these 'excess' galaxies and discuss the implications: these galaxies are hosted in dark matter haloes with circular velocities in excess of 100 km s⁻¹, and their suppression may require a deep rethinking of stellar feedback processes taking place in galaxy formation.     

Measurement of dwarf galaxies in the rich clusters Abell 665 and 963 at z=0.2

We show that the luminosity functions of the distant rich clusters Abell 665 ( z =0.182) and Abell 963 ( z =0.206) are flat or gradually rising down to M_R =−14, with α≈−1.2±0.4 [here α is the logarithmic slope of the luminosity function: φ( L )∝ L ^α at the faint end]. We do not confirm the steep luminosity functions (α≤−1.8) that have been recently proposed for these two clusters.
Several technical points are discussed in detail. In particular, we compute the corrections to the background contamination caused by gravitational lensing from the cluster dark matter, and show that the corrections are small unless we wish to determine variations in the luminosity function on small scales.
Recent observations have also shown that the field galaxy luminosity function at z ≈0.2 is also shallow between M_B =−19 and M_B =−13. Abell 665 and 963 are two of the richest clusters known at that redshift. We therefore propose that the galaxy luminosity function might be universal in this magnitude range at z =0.2.
The dwarf galaxies that we see in Abell 665 have a colour distribution that is strongly peaked at B − R =1.9. We compute K -corrections based on the spectral energy distributions of local galaxies, and show that these are probably dwarf spheroidal galaxies. This might suggest that the dwarf spheroidal population observed in Virgo already existed at z =0.2.     

Recent arrival of faint cluster galaxies on the red sequence: luminosity functions from 119 deg2 of CFHTLS

The global star formation rate has decreased significantly since   z ∼ 1  , for reasons that are not well understood. Red-sequence galaxies, dominating in galaxy clusters, represent the population that have had their star formation shut off, and may therefore be the key to this problem. In this work, we select 127 rich galaxy clusters at  0.17 ≤ z ≤ 0.36  , from 119 deg² of the Canada–France–Hawaii Telescope Legacy Survey (CFHTLS) optical imaging data, and construct the r '-band red-sequence luminosity functions (LFs). We show that the faint end of the LF is very sensitive to how red-sequence galaxies are selected, and an optimal way to minimize the contamination from the blue cloud is to mirror galaxies on the redder side of the colour–magnitude relation. The LFs of our sample have a significant inflexion centred at     , suggesting a mixture of two populations. Combining our survey with low-redshift samples constructed from the Sloan Digital Sky Survey, we show that there is no strong evolution of the faint end of the LF (or the red-sequence dwarf-to-giant ratio) over the redshift range  0.2 ≲ z ≲ 0.4  , but from   z ∼ 0.2  to ∼0 the relative number of red-sequence dwarf galaxies has increased by a factor of ∼3, implying a significant build-up of the faint end of the cluster red sequence over the last 2.5 Gyr.     

Luminosity distributions within rich clusters — III. A comparative study of seven Abell/ACO clusters

We recover the luminosity distributions over a wide range of absolute magnitude (−24.5 <  M _R  < −16.5) for a sample of seven rich southern galaxy clusters. We find a large variation in the ratio of dwarf to giant galaxies, DGR: 0.8 ≤ DGR ≤ 3.1. This variation is shown to be inconsistent with a ubiquitous cluster luminosity function. The DGR shows a smaller variation from cluster to cluster in the inner regions ( r  ≲ 0.56 Mpc). Outside these regions we find the DGR to be strongly anticorrelated with the mean local projected galaxy density, with the DGR increasing towards lower densities. In addition, the DGR in the outer regions shows some correlation with Bautz–Morgan type. Radial analysis of the clusters indicates that the dwarf galaxies are less centrally clustered than the giants, and they form a significant halo around clusters. We conclude that measurements of the total cluster luminosity distribution based on the inner core alone are likely to be severe underestimates of the dwarf component, the integrated cluster luminosity and the contribution of galaxy masses to the cluster's total mass. Further work is required to quantify this. The observational evidence that the unrelaxed, lower density outer regions of clusters are dwarf-rich adds credence to the recent evidence and conjecture that the field is a predominantly dwarf-rich environment, and that the dwarf galaxies are under-represented in measures of the local field luminosity function.     

Supersonic motions of galaxies in clusters

We study motions of galaxies in galaxy clusters formed in the concordance Λ cold dark matter cosmology. We use high-resolution cosmological simulations that follow the dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing the motions of galaxies and the properties of intracluster gas in a sample of eight simulated clusters at z = 0, we study the velocity dispersion profiles of the dark matter, gas and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ≈1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move slightly faster than the dark matter particles. The magnitude of the velocity bias,   b _v ≈ 1.1  , is, however, smaller than the bias estimated for subhaloes in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ≈20–30 per cent of the gas sound speed. These random motions provide about 10 per cent of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.     

Galaxy–dark matter correlations applied to galaxy–galaxy lensing: predictions from the semi-analytic galaxy formation models

We use semi-analytic models of galaxy formation combined with high-resolution N -body simulations to make predictions for galaxy–dark matter correlations and apply them to galaxy–galaxy lensing. We analyse cross-power spectra between the dark matter and different galaxy samples selected by luminosity, colour or star formation rate. We compare the predictions with the recent detection by the Sloan Digital Sky Survey (SDSS). We show that the correlation amplitude and the mean tangential shear depend strongly on the luminosity of the sample on scales below 1  h ⁻¹ Mpc, reflecting the correlation between the galaxy luminosity and the halo mass. The cross-correlation cannot, however, be used to infer the halo profile directly because different halo masses dominate on different scales and because not all galaxies are at the centres of the corresponding haloes. We compute the redshift evolution of the cross-correlation amplitude and compare it with those of galaxies and dark matter. We also compute the galaxy–dark matter correlation coefficient and show that it is close to unity on scales above 1  h ⁻¹ Mpc for all considered galaxy types. This would allow one to extract the bias and the dark matter power spectrum on large scales from the galaxy and galaxy–dark matter correlations.     

The effects of photoionization on galaxy formation – I. Model and results at z=0

We develop a coupled model for the evolution of the global properties of the intergalactic medium (IGM) and the formation of galaxies, in the presence of a photoionizing background due to stars and quasars. We use this model to predict the thermodynamic history of the IGM when photoionized by galaxies forming in a cold dark matter (CDM) universe. The evolution of the galaxies is calculated using a semi-analytical model, including a detailed treatment of the effects of tidal stripping and dynamical friction on satellite galaxies orbiting inside larger dark matter haloes. We include in the model the negative feedback on galaxy formation from the photoionizing background. Photoionization inhibits galaxy formation in low-mass dark matter haloes in two ways: (i) heating of the IGM and inhibition of the collapse of gas into dark haloes by the IGM pressure, and (ii) reduction in the rate of radiative cooling of gas within haloes. The result of our method is a self-consistent model of galaxy formation and the IGM. The IGM is reheated twice (during reionization of H  i and He  ii ), and we find that the star formation rate per unit volume is slightly suppressed after each episode of reheating. We find that galaxies brighter than L _★ are mostly unaffected by reionization, while the abundance of faint galaxies is significantly reduced, leading to present-day galaxy luminosity functions with shallow faint-end slopes, in good agreement with recent observational data. Reionization also affects other properties of these faint galaxies, in a readily understandable way.     

Near-infrared bulge–disc correlations of lenticular galaxies

We consider the luminosity and environmental dependence of structural parameters of lenticular galaxies in the near-infrared K band. Using a 2D galaxy image decomposition technique, we extract bulge and disc structural parameters for a sample of 36 lenticular galaxies observed by us in the K band. By combining data from the literature for field and cluster lenticulars with our data, we study correlations between parameters that characterize the bulge and the disc as a function of luminosity and environment. We find that scaling relations such as the Kormendy relation, photometric plane and other correlations involving bulge and disc parameters show a luminosity dependence. This dependence can be explained in terms of galaxy formation models in which faint lenticulars  ( M _T > −24.5)  formed via secular formation processes that likely formed the pseudo-bulges of late-type disc galaxies, while brighter lenticulars  ( M _T < −24.5)  formed through a different formation mechanism most likely involving major mergers. On probing variations in lenticular properties as a function of environment, we find that faint cluster lenticulars show systematic differences with respect to faint field lenticulars. These differences support the idea that the bulge and disc components fade after the galaxy falls into a cluster, while simultaneously undergoing a transformation from spiral to lenticular morphologies.     

Galaxy cluster masses without non-baryonic dark matter

We apply the modified acceleration law obtained from Einstein gravity coupled to a massive skew symmetric field,   F _μνλ  , to the problem of explaining X-ray galaxy cluster masses without exotic dark matter. Utilizing X-ray observations to fit the gas mass profile and temperature profile of the hot intracluster medium (ICM) with King 'β-models', we show that the dynamical masses of the galaxy clusters resulting from our modified acceleration law fit the cluster gas masses for our sample of 106 clusters without the need of introducing a non-baryonic dark matter component. We are further able to show for our sample of 106 clusters that the distribution of gas in the ICM as a function of radial distance is well fitted by the dynamical mass distribution arising from our modified acceleration law without any additional dark matter component. In a previous work, we applied this theory to galaxy rotation curves and demonstrated good fits to our sample of 101 low surface brightness, high surface brightness and dwarf galaxies including 58 galaxies that were fitted photometrically with the single-parameter mass-to-light ratio ( M / L )_stars. The results obtained there were qualitatively similar to those obtained using Milgrom's phenomenological Modified Newtonian Dynamics (MOND) model, although the determined galaxy masses were quantitatively different, and MOND does not show a return to Keplerian behaviour at extragalactic distances. The results obtained here are compared to those obtained using Milgrom's phenomenological MOND model which does not fit the X-ray galaxy cluster masses unless an auxiliary dark matter component is included.     

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_⊙.     

Multiwavelength observations of an evolved galaxy group: an end-point of galaxy merging?

The group of galaxies RXJ1340.6+4018 has approximately the same bolometric X-ray luminosity as other bright galaxy groups and poor clusters such as the Virgo cluster. However, 70 per cent of the optical luminosity of the group comes from a dominant giant elliptical galaxy, compared with 5 per cent from M87 in Virgo.The second brightest galaxy in RXJ1340.6+4018 is a factor of 10 fainter (Δ m ₁₂=2.5 mag) than the dominant elliptical, and the galaxy luminosity function has a gap at about L *.
We interpret the properties of the system as a result of galaxy merging within a galaxy group. We find that the central galaxy lies on the Fundamental Plane of ellipticals, has an undisturbed, non-cD morphology, and has no spectral features indicative of recent star formation, suggesting that the last major merger occurred ≳4 Gyr ago. The deviation of the system from the cluster L _X− T relation in the opposite sense to most groups may be caused by an early epoch of formation of the group or a strong cooling flow.
The unusual elongation of the X-ray isophotes and the similarity between the X-ray and optical ellipticities at large radii (∼230 kpc) suggest that both the X-ray gas and the outermost stars of the dominant galaxy are responding to an elongated dark matter distribution. RXJ1340.6+4018 may be part of a filamentary structure related to infall in the outskirts of the cluster A1774.     

The growth and assembly of a massive galaxy at z∼ 2

We study the stellar mass assembly of the Spiderweb galaxy  (MRC 1138−262)  , a massive   z = 2.2  radio galaxy in a protocluster and the probable progenitor of a brightest cluster galaxy. Nearby protocluster galaxies are identified and their properties are determined by fitting stellar population models to their rest-frame ultraviolet to optical spectral energy distributions. We find that within 150 kpc of the radio galaxy the stellar mass is centrally concentrated in the radio galaxy, yet most of the dust-uncorrected, instantaneous star formation occurs in the surrounding low-mass satellite galaxies. We predict that most of the galaxies within 150 kpc of the radio galaxy will merge with the central radio galaxy by   z = 0  , increasing its stellar mass by up to a factor of ≃2. However, it will take several hundred Myr for the first mergers to occur, by which time the large star formation rates are likely to have exhausted the gas reservoirs in the satellite galaxies. The tidal radii of the satellite galaxies are small, suggesting that stars and gas are being stripped and deposited at distances of tens of kpc from the central radio galaxy. These stripped stars may become intracluster stars or form an extended stellar halo around the radio galaxy, such as those observed around cD galaxies in cluster cores.     

Central mass-to-light ratios and dark matter fractions in early-type galaxies

Dynamical studies of local elliptical galaxies and the Fundamental Plane point to a strong dependence of the total mass-to-light ratio ( M / L ) on luminosity with a relation of the form   M / L ∝ L ^γ  . The 'tilt'γ may be caused by various factors, including stellar population properties (metallicity, age and star formation history), initial mass function, rotational support, luminosity profile non-homology and dark matter (DM) fraction. We evaluate the impact of all these factors using a large uniform data set of local early-type galaxies from Prugniel & Simien. We take particular care in estimating the stellar masses, using a general star formation history, and comparing different population synthesis models. We find that the stellar M / L contributes little to the tilt. We estimate the total M / L using simple Jeans dynamical models, and find that adopting accurate luminosity profiles is important but does not remove the need for an additional tilt component, which we ascribe to DM. We survey trends of the DM fraction within one effective radius, finding it to be roughly constant for galaxies fainter than   M _B∼−20.5  , and increasing with luminosity for the brighter galaxies; we detect no significant differences between S0s and fast- and slow-rotating ellipticals. We construct simplified cosmological mass models and find general consistency, where the DM transition point is caused by a change in the relation between luminosity and effective radius. A more refined model with varying galaxy star formation efficiency suggests a transition from total mass profiles (including DM) of faint galaxies distributed similarly to the light to near-isothermal profiles for the bright galaxies. These conclusions are sensitive to various systematic uncertainties which we investigate in detail, but are consistent with the results of dynamical studies at larger radii.     

Satellite kinematics – II. The halo mass–luminosity relation of central galaxies in SDSS

The kinematics of satellite galaxies reflect the masses of the extended dark matter haloes in which they orbit, and thus shed light on the mass–luminosity relation (MLR) of their corresponding central galaxies. In this paper, we select a large sample of centrals and satellites from the Sloan Digital Sky Survey and measure the kinematics (velocity dispersions) of the satellite galaxies as a function of the r -band luminosity of the central galaxies. Using the analytical framework presented in More, van den Bosch & Cacciato, we use these data to infer both the mean and the scatter of the MLR of central galaxies, carefully taking account of selection effects and biases introduced by the stacking procedure. As expected, brighter centrals on average reside in more massive haloes. In addition, we find that the scatter in halo masses for centrals of a given luminosity,  σ_{log  M}   , also increases with increasing luminosity. As we demonstrate, this is consistent with  σ_{log  L}   , which reflects the scatter in the conditional probability function   P ( L _c| M )  , being independent of halo mass. Our analysis of the satellite kinematics yields  σ_{log  L} = 0.16  ±  0.04  , in excellent agreement with constraints from clustering and group catalogues, and with predictions from a semi-analytical model of galaxy formation. We thus conclude that the amount of stochasticity in galaxy formation, which is characterized by  σ_{log  L}   , is well constrained, independent of halo mass and in a good agreement with current models of galaxy formation.     

The effects of photoionization on galaxy formation – II. Satellite galaxies in the Local Group

We use a self-consistent model of galaxy formation and the evolution of the intergalactic medium to study the effects of the reionization of the Universe at high redshift on the properties of satellite galaxies like those seen around the Milky Way. Photoionization suppresses the formation of small galaxies, so that surviving satellites are preferentially those that formed before the Universe reionized. As a result, the number of satellites expected today is about an order of magnitude smaller than the number inferred by identifying satellites with subhaloes of the same circular velocity in high-resolution simulations of the dark matter. The resulting satellite population has an abundance similar to that observed in the Local Group, although the distribution of circular velocities differs somewhat from the available data. We explore many other properties of satellite galaxies, including their gas content, metallicity and star formation rate, and find generally good agreement with available data. Our model predicts the existence of many as yet undetected satellites in the Local Group. We quantify their observability in terms of their apparent magnitude and surface brightness, and also in terms of their constituent stars. A near-complete census of the Milky Way's satellites would require imaging to   V ≈20  and to a surface brightness fainter than 26 V -band magnitudes per square arcsecond. Satellites with integrated luminosity   V =15  should contain of order 100 stars brighter than   B =26  , with central stellar densities of a few tens per square arcminute. Discovery of a large population of faint satellites would provide a strong test of current models of galaxy formation.     

Light curves for off-centre ignition models of Type Ia supernovae

We compare observations of the high-redshift galaxy population to the predictions of the galaxy formation model of Croton et al. and De Lucia & Blaizot. This model, implemented on the Millennium Simulation of the concordance Lambda cold dark matter cosmogony, introduces 'radio mode' feedback from the central galaxies of groups and clusters in order to obtain quantitative agreement with the luminosity, colour, morphology and clustering properties of the present-day galaxy population. Here we construct deep light cone surveys in order to compare model predictions to the observed counts and redshift distributions of distant galaxies, as well as to their inferred luminosity and mass functions out to redshift 5. With the exception of the mass functions, all these properties are sensitive to modelling of dust obscuration. A simple but plausible treatment agrees moderately well with most of the data. The predicted abundance of relatively massive  (∼ M _*)  galaxies appears systematically high at high redshift, suggesting that such galaxies assemble earlier in this model than in the real Universe. An independent galaxy formation model implemented on the same simulation matches the observed mass functions slightly better, so the discrepancy probably reflects incomplete or inaccurate galaxy formation physics rather than problems with the underlying cosmogony.     

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

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