Near-infrared luminosity function and colours of dwarf galaxies in the Coma cluster

We present K -band observations of the low-luminosity galaxies in the Coma cluster, which are responsible for the steep upturn in the optical luminosity function at M _R∼−16, discovered recently. The main results of this study are as follows.
(i) The optical–near-infrared colours of these galaxies imply that they are dwarf spheroidal galaxies. The median B − K colour for galaxies with −19.3< M_K <−16.3 is 3.6 mag.
(ii) The K -band luminosity function in the Coma cluster is not well constrained, because of the uncertainties due to the field-to-field variance of the background. However, within the estimated large errors, this is consistent with the R -band luminosity function, shifted by ∼3 mag.
(iii) Many of the cluster dwarfs lie in a region of the B − K versus B − R colour–colour diagram where background galaxies are rare ( B − K <5; 1.2< B − R <1.6). Local dwarf spheroidal galaxies lie in this region too. This suggests that a better measurement of the K -band cluster luminosity can be made if the field-to-field variance of the background can be measured as a function of colour, even if it is large.
(iv) If we assume that none of the galaxies in the region of the B − K versus B − R plane given in (iii) in our cluster fields are background, and that all the cluster galaxies with 15.5< K <18.5 lie in this region of the plane, then we measure α=−1.41^+0.34_−0.37 for −19.3< M_K −16.3, where α is the logarithmic slope of the luminosity function. The uncertainties in this number come from counting statistics.     

Evolution of the u-band luminosity function from redshift 1.2 to 0

We produce and analyse u -band (  λ≈ 355  nm) luminosity functions (LFs) for the red and blue populations of galaxies using data from the Sloan Digital Sky Survey (SDSS) u -band Galaxy Survey ( u GS) and Deep Evolutionary Exploratory Probe 2 (DEEP2) survey. From a spectroscopic sample of 41 575 SDSS u GS galaxies and 24 561 DEEP2 galaxies, we produce colour magnitude diagrams and make use of the colour bimodality of galaxies to separate red and blue populations. LFs for eight redshift slices in the range  0.01 < z < 1.2  are determined using the  1/ V _max  method and fitted with Schechter functions showing that there is significant evolution in   M *  , with a brightening of 1.4 mag for the combined population. The integration of the Schechter functions yields the evolution in the u -band luminosity density (LD) out to   z ∼ 1  . By parametrizing the evolution as  ρ∝ (1 + z )^β  , we find that  β= 1.36 ± 0.2  for the combined populations and  β= 2.09 ± 0.2  for the blue population. By removing the contribution of the old stellar population to the u -band LD and correcting for dust attenuation, we estimate the evolution in the star formation rate (SFR) of the Universe to be  β_SFR= 2.5 ± 0.3  . Discrepancies between our result and higher evolution rates measured using the infrared and far-UV can be reconciled by considering possibilities such as an underestimated dust correction at high redshifts or evolution in the stellar initial mass function.     

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).     

The quasar epoch and the stellar ages of early-type galaxies

We investigate the hypothesis that quasars formed together with the stellar populations of early-type galaxies. This hypothesis – in conjunction with the stellar ages of early-type galaxies from population synthesis models, the relation of black hole mass to bulge velocity dispersion, and the velocity dispersion distribution of spheroids from the Sloan Digital Sky Survey – completely determines the cosmic accretion history of supermassive black holes and the redshift evolution of the characteristic luminosity. On the other hand, the precise shape of the luminosity function of quasars depends on the light curve of quasars and – in the optical, but not so much in X-rays – on the covering factor of the dust surrounding the active nucleus. We find a plausible set of assumptions for which the coeval formation of supermassive black holes and elliptical galaxies is in good agreement with the observed B -band and X-ray luminosity functions of quasars.     

Lyman break galaxies at z= 4–6 in cosmological smoothed particle hydrodynamics simulations

We perform a spectrophotometric analysis of galaxies at redshifts z = 4–6 in cosmological smoothed particle hydrodynamics simulations of a Λ cold dark matter universe. Our models include radiative cooling and heating by a uniform ultraviolet (UV) background, star formation, supernova feedback, and a phenomenological model for galactic winds. Analysing a series of simulations of varying box size and particle number allows us to isolate the impact of numerical resolution on our results. Specifically, we determine the luminosity functions in B , V , R , i ' and z ' filters, and compare the results with observational surveys of Lyman break galaxies (LBGs) performed with the Subaru telescope and the Hubble Space Telescope . We find that the simulated galaxies have UV colours consistent with observations and fall in the expected region of the colour–colour diagrams used by the Subaru group. The stellar masses of the most massive galaxies in our largest simulation increase their stellar mass from   M _★∼ 10¹¹ M_⊙  at z = 6 to   M _★∼ 10^11.7 M_⊙  at z = 3. Assuming a uniform extinction of E ( B − V ) = 0.15, we also find reasonable agreement between simulations and observations in the space density of UV bright galaxies at z = 3–6, down to the magnitude limit of each survey. For the same moderate extinction level of E ( B − V ) ∼ 0.15, the simulated luminosity functions match observational data, but have a steep faint-end slope with α∼−2.0. We discuss the implications of the steep faint-end slope found in the simulations. Our results confirm the generic conclusion from earlier numerical studies that UV bright LBGs at z ≥ 3 are the most massive galaxies with E ( B − V ) ∼ 0.15 at each epoch.     

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.     

Luminosity and surface brightness distribution of K-band galaxies from the UKIDSS Large Area Survey

We present luminosity and surface-brightness distributions of 40 111 galaxies with K -band photometry from the United Kingdom Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS), Data Release 3 and optical photometry from Data Release 5 of the Sloan Digital Sky Survey (SDSS). Various features and limitations of the new UKIDSS data are examined, such as a problem affecting Petrosian magnitudes of extended sources. Selection limits in K - and r -band magnitude, K -band surface brightness and K -band radius are included explicitly in the  1/ V _max  estimate of the space density and luminosity function. The bivariate brightness distribution in K -band absolute magnitude and surface brightness is presented and found to display a clear luminosity–surface brightness correlation that flattens at high luminosity and broadens at low luminosity, consistent with similar analyses at optical wavelengths. Best-fitting Schechter function parameters for the K -band luminosity function are found to be   M *− 5 log  h =−23.19 ± 0.04, α=−0.81 ± 0.04  and  φ*= (0.0166 ± 0.0008)  h ³ Mpc⁻³  , although the Schechter function provides a poor fit to the data at high and low luminosity, while the luminosity density in the K band is found to be   j = (6.305 ± 0.067) × 10⁸ L_⊙  h  Mpc⁻³  . However, we caution that there are various known sources of incompleteness and uncertainty in our results. Using mass-to-light ratios determined from the optical colours, we estimate the stellar mass function, finding good agreement with previous results. Possible improvements are discussed that could be implemented when extending this analysis to the full LAS.     

Mass profiles and anisotropies of early-type galaxies

We discuss the problem of using stellar kinematics of early-type galaxies to constrain the orbital anisotropies and radial mass profiles of galaxies. We demonstrate that compressing the light distribution of a galaxy along the line of sight produces approximately the same signature in the line-of-sight velocity profiles as radial anisotropy. In particular, fitting spherically symmetric dynamical models to apparently round, isotropic face-on flattened galaxies leads to a spurious bias towards radial orbits in the models, especially if the galaxy has a weak face-on stellar disc. Such face-on stellar discs could plausibly be the cause of the radial anisotropy found in spherical models of intermediate luminosity ellipticals such as NGC 2434, 3379 and 6703.
In the light of this result, we use simple dynamical models to constrain the outer mass profiles of a sample of 18 round, early-type galaxies. The galaxies follow a Tully–Fisher relation parallel to that for spiral galaxies, but fainter by at least 0.8 mag ( I -band) for a given mass. The most luminous galaxies show clear evidence for the presence of a massive dark halo, but the case for dark haloes in fainter galaxies is more ambiguous. We discuss the observations that would be required to resolve this ambiguity.     

The high-mass end of the Tully–Fisher relation

We study the location of massive disc galaxies on the Tully–Fisher (TF) relation. Using a combination of K -band photometry and high-quality rotation curves, we show that in traditional formulations of the TF relation (using the width of the global H  i profile or the maximum rotation velocity), galaxies with rotation velocities larger than 200 km s⁻¹ lie systematically to the right of the relation defined by less massive systems, causing a characteristic 'kink' in the relations. Massive, early-type disc galaxies in particular have a large offset, up to 1.5 mag, from the main relation defined by less massive and later-type spirals.
The presence of a change in slope at the high-mass end of the TF relation has important consequences for the use of the TF relation as a tool for estimating distances to galaxies or for probing galaxy evolution. In particular, the luminosity evolution of massive galaxies since z ≈ 1 may have been significantly larger than estimated in several recent studies.
We also show that many of the galaxies with the largest offsets have declining rotation curves and that the change in slope largely disappears when we use the asymptotic rotation velocity as kinematic parameter. The remaining deviations from linearity can be removed when we simultaneously use the total baryonic mass (stars + gas) instead of the optical or near-infrared luminosity. Our results strengthen the view that the TF relation fundamentally links the mass of dark matter haloes with the total baryonic mass embedded in them.     

Dark matter in elliptical galaxies – I. Is the total mass density profile of the NFW form or even steeper?

Elliptical galaxies are modelled as Sérsic luminosity distributions with density profiles (DPs) for the total mass adopted from the DPs of haloes within dissipationless ΛCDM (cold dark matter) N -body simulations. Ellipticals turn out to be inconsistent with cuspy low-concentration NFW models representing the total mass distribution, neither are they consistent with a steeper −1.5 inner slope, nor with the shallower models proposed by Navarro et al., nor with NFW models 10 times more concentrated than predicted, as deduced from several X-ray observations – the mass models, extrapolated inwards, lead to local mass-to-light ratios that are smaller than the stellar value inside an effective radius ( R _e), and to central aperture velocity dispersions that are much smaller than observed. This conclusion remains true as long as there is no sharp steepening (slope < −2) of the dark matter DPs just inside 0.01 virial radii.
The very low total mass and velocity dispersion produced within R _e by an NFW-like total mass profile suggests that the stellar component should dominate the dark matter component out to at least R _e. It should then be difficult to kinematically constrain the inner slope of the DP of ellipticals. The high-concentration parameters deduced from X-ray observations appear to be a consequence of fitting an NFW model to the total mass DP made up of a stellar component that dominates inside and a dark matter component that dominates outwards.
An appendix gives the virial mass dependence of the concentration parameter, central density and total mass of the Navarro et al. model. In a second appendix are given single integral expressions for the velocity dispersions averaged along the line of sight, in circular apertures and in thin slits, for general luminosity density and mass distributions, with isotropic orbits.     

Galaxy evolution in the high-redshift, colour-selected cluster RzCS 052 at z= 1.02

We present deep I and z ' imaging of the colour-selected cluster RzCS 052 and study the colour–magnitude relation of this cluster, its scatter, the morphological distribution on the red sequence, the luminosity and stellar mass functions of red galaxies and the cluster blue fraction. We find that the stellar populations of early-type galaxies in this cluster are uniformly old and that their luminosity function does not show any sign of evolution other than the passive evolution of their stellar populations. We rule out a significant contribution from mergers in the build-up of the red sequence of RzCS 052. The cluster has a large (∼30 per cent) blue fraction and we infer that the evolution of the blue galaxies is faster than an exponentially declining star formation model and that these objects have probably experienced starburst episodes. Mergers are unlikely to be the driver of the observed colour evolution, because of the measured constancy of the mass function, as derived from near-infrared photometry of 32 clusters, including RzCS 052, presented in a related paper. Mechanisms with clustercentric radial dependent efficiencies are disfavoured as well, because of the observed constant blue fraction with clustercentric distance.     

A sample of 6C radio sources designed to find objects at redshift z>4– III. Imaging and the radio galaxy K–z relation

In this paper, the third and final of a series, we present complete K -band imaging and some complementary I -band imaging of the filtered 6C* sample. We find no systematic differences between the K – z relation of 6C* radio galaxies and those from complete samples, so the near-infrared properties of luminous radio galaxies are not obviously biased by the additional 6C* radio selection criteria (steep spectral index and small angular size). The 6C* K – z data significantly improve delineation of the K – z relation for radio galaxies at high redshift ( z >2) . Accounting for non-stellar contamination, and for correlations between radio luminosity and stellar mass, we find little support for previous claims that the underlying scatter in the stellar luminosity of radio galaxies increases significantly at z >2 . In a particular spatially flat universe with a cosmological constant (Ω_M=0.3 and Ω_Λ=0.7) , the most luminous radio sources appear to be associated with galaxies with a luminosity distribution with a high mean (≈5  L *), and a low dispersion ( σ ∼0.5 mag) which formed their stars at epochs corresponding to z ≳2.5 . This result is in line with recent submillimetre studies of high-redshift radio galaxies and the inferred ages of extremely red objects from faint radio samples.     

The relation between stellar mass and weak lensing signal around galaxies: implications for modified Newtonian dynamics

We study the amplitude of the weak gravitational lensing signal as a function of stellar mass around a sample of relatively isolated galaxies. This selection of lenses simplifies the interpretation of the observations, which consist of data from the Red-Sequence Cluster Survey and the Sloan Digital Sky Survey. We find that the amplitude of the lensing signal as a function of stellar mass is well described by a power law with a best-fitting slope  α= 0.74 ± 0.08  . This result is inconsistent with modified Newtonian dynamics (MOND), which predicts  α= 0.5  (we find  α > 0.5  with 99.7 per cent confidence). As a related test, we determine the MOND mass-to-light ratio as a function of luminosity. Our results require dark matter for the most luminous galaxies ( L ≳ 10¹¹ L_⊙). We rule out an extended halo of gas or active neutrinos as a way of reconciling our findings with MOND. Although we focus on a single alternative gravity model, we note that our results provide an important test for any alternative theory of gravity.     

Further constraints on the evolution of Ks-selected galaxies in the GOODS/CDFS field

We have selected and analysed the properties of a sample of  2905 K_s < 21.5  galaxies in  ∼131 arcmin²  of the Great Observatories Origins Deep Survey (GOODS) Chandra Deep Field South (CDFS), to obtain further constraints on the evolution of K_s -selected galaxies with respect to the results already obtained in previous studies. We made use of the public deep multiwavelength imaging from the optical B through the infrared (IR) 4.5-μm bands, in conjunction with available spectroscopic and COMBO17 data in the CDFS, to construct an optimized redshift catalogue for our galaxy sample. We computed the K_s -band luminosity function and determined that its characteristic magnitude has a substantial brightening and a decreasing total density from   z = 0  to  〈 z 〉= 2.5  . We also analysed the colours and number density evolution of galaxies with different stellar masses. Within our sample, and in contrast to what is observed for less massive systems, the vast majority (∼85–90 per cent) of the most massive  ( M > 2.5 × 10¹¹ M_⊙)  local galaxies appear to be in place before redshift   z ∼ 1  . Around 65–70 per cent of the total assemble between redshifts   z = 1  and 3 and most of them display extremely red colours, suggesting that plausible star formation in these very massive systems should mainly proceed in obscured, short-time-scale bursts. The remaining fraction (up to ∼20 per cent) could be in place at even higher redshifts   z = 3–4  , pushing the first epoch of formation of massive galaxies beyond the limits of current near-IR surveys.     

On the formation of massive galaxies: a simultaneous study of number density, size and intrinsic colour evolution in GOODS

The evolution of number density, size and intrinsic colour is determined for a volume-limited sample of visually classified early-type galaxies selected from the Hubble Space Telescope /Advanced Camera for Surveys images of the Great Observatories Origins Deep Survey (GOODS) North and South fields (version 2). The sample comprises 457 galaxies over 320 arcmin² with stellar masses above  3 × 10¹⁰ M_⊙  in the redshift range  0.4 < z < 1.2  . Our data allow a simultaneous study of number density, intrinsic colour distribution and size. We find that the most massive systems  (≳3 × 10¹¹ M_⊙)  do not show any appreciable change in comoving number density or size in our data. Furthermore, when including the results from 2dF galaxy redshift survey, we find that the number density of massive early-type galaxies is consistent with no evolution between   z = 1.2  and 0, i.e. over an epoch spanning more than half of the current age of the Universe. We find large discrepancies between the predictions of semi-analytic models. Massive galaxies show very homogeneous intrinsic colour distributions, with nearly flat radial colour gradients, but with a significant negative correlation between stellar mass and colour gradient, such that red cores appear predominantly in massive galaxies. The distribution of half-light radii – when compared to   z ∼ 0  and   z > 1  samples – is compatible with the predictions of semi-analytic models relating size evolution to the amount of dissipation during major mergers.     

The globular clusters–stellar haloes connection in early-type galaxies

This paper explores if, and to what an extent, the stellar populations of early-type galaxies can be traced through the colour distribution of their globular cluster (GC) systems. The analysis, based on a galaxy sample from the Virgo Advanced Camera for Surveys data, is an extension of a previous approach that has been successful in the cases of the giant ellipticals NGC 1399 and NGC 4486, and assumes that the two dominant GC populations form along diffuse stellar populations sharing the cluster chemical abundances and spatial distributions. The results show that (a) integrated galaxy colours can be matched to within the photometric uncertainties and are consistent with a narrow range of ages; (b) the inferred mass to luminosity ratios and stellar masses are within the range of values available in the literature; (c) most GC systems occupy a thick plane in the volume space defined by the cluster formation efficiency, total stellar mass and projected surface mass density. The formation efficiency parameter of the red clusters shows a dependency with projected stellar mass density that is absent for the blue globulars. In turn, the brightest galaxies appear clearly detached from that plane as a possible consequence of major past mergers; (d) the stellar mass–metallicity relation is relatively shallow but shows a slope change at   M _*≈ 10¹⁰ M_⊙  . Galaxies with smaller stellar masses show predominantly unimodal GC colour distributions. This result may indicate that less massive galaxies are not able to retain chemically enriched interstellar matter.     

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 effect of dwarf galaxy disruption in semi-analytic models

We present results for a galaxy formation model that includes a simple treatment for the disruption of dwarf galaxies by gravitational forces and galaxy encounters within galaxy clusters. This is implemented a posteriori in a semi-analytic model by considering the stability of cluster dark matter subhaloes at   z = 0  . We assume that a galaxy whose dark matter substructure has been disrupted will itself disperse, while its stars become part of the population of intracluster stars responsible for the observed intracluster light. Despite the simplicity of this assumption, our results show a substantial improvement over previous models and indicate that the inclusion of galaxy disruption is indeed a necessary ingredient of galaxy formation models. We find that galaxy disruption suppresses the number density of dwarf galaxies by about a factor of 2. This makes the slope of the faint end of the galaxy luminosity function shallower, in agreement with observations. In particular, the abundance of faint, red galaxies is strongly suppressed. As a result, the luminosity function of red galaxies and the distinction between the red and the blue galaxy populations in colour–magnitude relationships are correctly predicted. Finally, we estimate a fraction of intracluster light comparable to that found in clusters of galaxies.