A constant dark matter halo surface density in galaxies

We confirm and extend the recent finding that the central surface density  μ_0D≡ r ₀ρ₀  of galaxy dark matter haloes, where r ₀ and  ρ₀  are the halo core radius and central density, is nearly constant and independent of galaxy luminosity. Based on the co-added rotation curves (RCs) of ∼1000 spiral galaxies, the mass models of individual dwarf irregular and spiral galaxies of late and early types with high-quality RCs, and the galaxy–galaxy weak-lensing signals from a sample of spiral and elliptical galaxies, we find that  log μ_0D= 2.15 ± 0.2  in units of  log(M_⊙ pc⁻²)  . We also show that the observed kinematics of Local Group dwarf spheroidal galaxies are consistent with this value. Our results are obtained for galactic systems spanning over 14 mag, belonging to different Hubble types and whose mass profiles have been determined by several independent methods. In the same objects, the approximate constancy of  μ_0D  is in sharp contrast to the systematical variations, by several orders of magnitude, of galaxy properties, including  ρ₀  and central stellar surface density.     

Satellite kinematics – I. A new method to constrain the halo mass–luminosity relation of central galaxies

Satellite kinematics can be used to probe the masses of dark matter haloes of central galaxies. In order to measure the kinematics with sufficient signal-to-noise ratio, one uses the satellite galaxies of a large number of central galaxies stacked according to similar properties (e.g. luminosity). However, in general, the relation between the luminosity of a central galaxy and the mass of its host halo will have non-zero scatter. Consequently, this stacking results in combining the kinematics of satellite galaxies in haloes of different masses, which complicates the interpretation of the data. In this paper, we present an analytical framework to model satellite kinematics, properly accounting for this scatter and for various selection effects. We show that in the presence of scatter in the halo mass–luminosity relation, the commonly used velocity dispersion of satellite galaxies can not be used to infer a unique halo mass–luminosity relation. In particular, we demonstrate that there is a degeneracy between the mean and the scatter of the halo mass–luminosity relation. We present a new technique that can break this degeneracy, and which involves measuring the velocity dispersions using two different weighting schemes: host weighting (each central galaxy gets the same weight) and satellite weighting (each central galaxy gets a weight proportional to its number of satellites). The ratio between the velocity dispersions obtained using these two weighting schemes is sensitive to the scatter in the halo mass–luminosity relation, and can thus be used to infer a unique relation between light and mass from the kinematics of satellite galaxies.     

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.     

Mass-to-light ratio gradients in early-type galaxy haloes

Owing to the fact that the near future should see a rapidly expanding set of probes of the halo masses of individual early-type galaxies, we introduce a convenient parameter for characterizing the halo masses from both observational and theoretical results:  ∇_ℓϒ  , the logarithmic radial gradient of the mass-to-light ratio. Using halo density profiles from Λ-cold dark matter (CDM) simulations, we derive predictions for this gradient for various galaxy luminosities and star formation efficiencies  ε_SF  . As a pilot study, we assemble the available  ∇_ℓϒ  data from kinematics in early-type galaxies – representing the first unbiased study of halo masses in a wide range of early-type galaxy luminosities – and find a correlation between luminosity and  ∇_ℓϒ  , such that the brightest galaxies appear the most dark-matter dominated. We find that the gradients in most of the brightest galaxies may fit in well with the ΛCDM predictions, but that there is also a population of fainter galaxies whose gradients are so low as to imply an unreasonably high star formation efficiency  ε_SF > 1  . This difficulty is eased if dark haloes are not assumed to have the standard ΛCDM profiles, but lower central concentrations.     

Halo masses for optically selected and for radio-loud AGN from clustering and galaxy–galaxy lensing

We compute two-point correlation functions and measure the shear signal due to galaxy–galaxy lensing for 80 000 optically identified and 5700 radio-loud active galactic nuclei (AGN) from Data Release 4 of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed   M _*  . Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive haloes than optical AGN: in our samples, their mean halo masses are  1.6 × 10¹³  and  8 × 10¹¹  h ⁻¹ M_⊙  , respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter haloes of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in haloes where gas cools quasi-statically.     

H i and dark matter in the windy starburst dwarf galaxy NGC 1705

We present 21-cm H  i line observations of the blue compact dwarf galaxy NGC 1705. Previous optical observations show a strong outflow powered by an ongoing starburst dominating the H  ii morphology and kinematics. In contrast, most of the H  i lies in a rotating disc. An extraplanar H  i spur accounts for ∼8 per cent of the total H  i mass, and is possibly associated with the H  ii outflow. The inferred mass loss rate out of the core of the galaxy is significant, ∼0.2 − 2 M_⊙ yr⁻¹, but does not dominate the H  i dynamics. Mass model fits to the rotation curve show that the dark matter (DM) halo is dominant at nearly all radii and has a central density ρ₀ ≈ 0.1 M_⊙ pc⁻³: ten times higher than typically found in dwarf irregular galaxies, but similar to the only other mass-modelled blue compact dwarf, NGC 2915. This large difference strongly indicates that there is little evolution between dwarf irregular and blue compact dwarf types. Instead, dominant DM haloes may regulate the morphology of dwarf galaxies by setting the critical surface density for disc star formation. Neither our data nor catalogue searches reveal any likely external trigger to the starburst in NGC 1705.     

Constraints on dark matter physics from dwarf galaxies through galaxy cluster haloes

One of the predictions of the standard cold dark matter model is that dark haloes have centrally divergent density profiles. An extensive body of rotation curve observations of dwarf and low surface brightness galaxies shows the dark haloes of those systems to be characterized by soft constant-density central cores. Several physical processes have been proposed to produce soft cores in dark haloes, each one with different scaling properties. With the aim of discriminating among them we have examined the rotation curves of dark-matter-dominated dwarf and low surface brightness galaxies and the inner mass profiles of two clusters of galaxies lacking a central cD galaxy and with evidence of soft cores in the centre. The core radii and central densities of these haloes scale in a well-defined manner with the depth of their potential wells, as measured through the maximum circular velocity. As a result of our analysis we identify self-interacting cold dark matter as a viable solution to the core problem, where a non-singular isothermal core is formed in the halo centre surrounded by a Navarro, Frenk & White profile in the outer parts. We show that this particular physical situation predicts core radii in agreement with observations. Furthermore, using the observed scalings, we derive an expression for the minimum cross-section ( σ ) which has an explicit dependence with the halo dispersion velocity ( v ). If m _x is the mass of the dark matter particle: σ m _x ≈4×10⁻²⁵ (100 km s⁻¹  v ⁻¹) cm² GeV⁻¹.     

The growth of supermassive black holes in pseudo-bulges, classical bulges and elliptical galaxies

Using results from structural analysis of a sample of nearly 1000 local galaxies from the Sloan Digital Sky Survey, we estimate how the mass in central black holes is distributed amongst elliptical galaxies, classical bulges and pseudo-bulges, and investigate the relation between their stellar masses and central stellar velocity dispersion σ. Assuming a single relation between elliptical galaxy/bulge mass, M _Bulge, and central black hole mass, M _BH, we find that  55⁺⁸₋₄  per cent of the mass in black holes in the local universe is in the centres of elliptical galaxies,  41⁺⁴₋₂  per cent in classical bulges and  4^+0.9_−0.4  per cent in pseudo-bulges. We find that ellipticals, classical bulges and pseudo-bulges follow different relations between their stellar masses and σ, and the most significant offset occurs for pseudo-bulges in barred galaxies. This structural dissimilarity leads to discrepant black hole masses if single   M _BH– M _Bulge  and   M _BH–σ  relations are used. Adopting relations from the literature, we find that the   M _BH–σ  relation yields an estimate of the total mass density in black holes that is roughly 55 per cent larger than if the   M _BH– M _Bulge  relation is used.     

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.     

Where do 'red and dead' early-type void galaxies come from?

Void regions of the Universe offer a special environment for studying cosmology and galaxy formation, which may expose weaknesses in our understanding of these phenomena. Although galaxies in voids are observed to be predominately gas rich, star forming and blue, a subpopulation of bright-red void galaxies can also be found, whose star formation was shutdown long ago. Are the same processes that quench star formation in denser regions of the Universe also at work in voids?
We compare the luminosity function of void galaxies in the 2dF Galaxy Redshift Survey, to those from a galaxy formation model built on the Millennium simulation. We show that a global star formation suppression mechanism in the form of low-luminosity 'radio-mode' active galactic nuclei (AGN) heating is sufficient to reproduce the observed population of void early types. Radio-mode heating is environment independent other than its dependence on dark matter halo mass, where, above a critical mass threshold of approximately   M _vir∼ 10^12.5 M_⊙  , gas cooling on to the galaxy is suppressed and star formation subsequently fades. In the Millennium simulation, the void halo mass function is shifted with respect to denser environments, but still maintains a high-mass tail above this critical threshold. In such void haloes, radio-mode heating remains efficient and red galaxies are found; collectively these galaxies match the observed space density without any modification to the model. Consequently, galaxies living in vastly different large-scale environments but hosted by haloes of similar mass are predicted to have similar properties, consistent with observations.     

Simulating subhaloes at high redshift: merger rates, counts and types

Galaxies are believed to be in one-to-one correspondence with simulated dark matter subhaloes. We use high-resolution N -body simulations of cosmological volumes to calculate the statistical properties of subhalo (galaxy) major mergers at high redshift ( z = 0.6–5). We measure the evolution of the galaxy merger rate, finding that it is much shallower than the merger rate of dark matter host haloes at   z > 2.5  , but roughly parallels that of haloes at   z < 1.6  . We also track the detailed merger histories of individual galaxies and measure the likelihood of multiple mergers per halo or subhalo. We examine satellite merger statistics in detail: 15–35 per cent of all recently merged galaxies are satellites, and satellites are twice as likely as centrals to have had a recent major merger. Finally, we show how the differing evolution of the merger rates of haloes and galaxies leads to the evolution of the average satellite occupation per halo, noting that for a fixed halo mass, the satellite halo occupation peaks at   z ∼ 2.5  .     

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.     

The black hole mass–stellar velocity dispersion correlation: bulges versus pseudo-bulges

We investigate the correlation between the supermassive black holes (SMBHs) mass ( M _bh) and the stellar velocity dispersion  (σ_*)  in two types of host galaxies: the early-type bulges (disc galaxies with classical bulges or elliptical galaxies) and pseudo-bulges. In the form  log ( M _bh/M_⊙) =α+β log (σ_*/200 km s⁻¹)  , the best-fitting results for the 39 early-type bulges are the slope  β= 4.06 ± 0.28  and the normalization  α= 8.28 ± 0.05  ; the best-fitting results for the nine pseudo-bulges are  β= 4.5 ± 1.3  and  α= 7.50 ± 0.18  . Both relations have intrinsic scatter in  log  M _bh  of ≲0.27 dex. The   M _bh–σ_*  relation for pseudo-bulges is different from the relation in the early-type bulges over the 3σ significance level. The contrasting relations indicate the formation and growth histories of SMBHs depend on their host type. The discrepancy between the slope of the   M _bh–σ_*  relations using different definition of velocity dispersion vanishes in our sample, a uniform slope will constrain the coevolution theories of the SMBHs and their host galaxies more effectively. We also find the slope for the 'core' elliptical galaxies at the high-mass range of the relation appears steeper  (β≃ 5–6)  , which may be the imprint of their origin of dissipationless mergers.     

Satellite systems around galaxies in hydrodynamic simulations

We investigate the properties of satellite galaxies formed in N -body/SPH simulations of galaxy formation in the ΛCDM cosmology. The simulations include the main physical effects thought to be important in galaxy formation and, in several cases, produce realistic spiral discs. In total, a sample of nine galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes brighter than the resolution limit,   M_V =−12  , the luminosity function of the satellite galaxies in the simulations is in excellent agreement with data for the Local Group. The radial number density profile of the model satellites, as well as their gas fractions also match observations very well. In agreement with previous N -body studies, we find that the satellites tend to be distributed in highly flattened configurations whose major axis is aligned with the major axis of the (generally triaxial) dark halo. In two out of three systems with sufficiently large satellite populations, the satellite system is nearly perpendicular to the plane of the galactic disc, a configuration analogous to that observed in the Milk Way. The discs themselves are perpendicular to the minor axis of their host haloes in the inner parts, and the correlation between the orientation of the galaxy and the shape of the halo persists even out to the virial radius. However, in one case the disc's minor axis ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular momenta of the galaxies and their host halo tend to be well aligned.     

Modelling galaxy–galaxy weak lensing with Sloan Digital Sky Survey groups

We use galaxy groups selected from the Sloan Digital Sky Survey (SDSS) together with mass models for individual groups to study the galaxy–galaxy lensing signals expected from galaxies of different luminosities and morphological types. We compare our model predictions with the observational results obtained from the SDSS by Mandelbaum et al. for the same samples of galaxies. The observational results are well reproduced in a Λ cold dark matter (ΛCDM) model based on the Wilkinson Microwave Anisotropy Probe ( WMAP ) 3-yr data, but a ΛCDM model with higher σ₈, such as the one based on the WMAP 1-yr data, significantly overpredicts the galaxy–galaxy lensing signal. We model, separately, the contributions to the galaxy–galaxy lensing signals from different galaxies: central versus satellite, early type versus late type and galaxies in haloes of different masses. We also examine how the predicted galaxy–galaxy lensing signal depends on the shape, density profile and the location of the central galaxy with respect to its host halo.     

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.     

Kinematic properties and stellar populations of faint early-type galaxies – I. Velocity dispersion measurements of central Coma galaxies

We present velocity dispersion measurements for 69 faint early-type galaxies in the core of the Coma cluster, spanning  −22.0 ≲ M_R ≲−17.5 mag  . We examine the   L –σ  relation for our sample and compare it to that of bright elliptical galaxies (Es) from the literature. The distribution of the the faint early-type galaxies in the   L –σ  plane follows the relation   L ∝σ^2.01±0.36  , which is significantly shallower from   L ∝σ⁴  as defined for the bright Es. While increased rotational support for fainter early-type galaxies could account for some of the difference in slope, we show that it cannot explain it. We also investigate the colour–σ relation for our Coma galaxies. Using the scatter in this relation, we constrain the range of galaxy ages as a function of their formation epoch for different formation scenarios. Assuming a strong coordination in the formation epoch of faint early-type systems in Coma, we find that most had to be formed at least 6 Gyr ago and over a short 1-Gyr period.     

An 80-kpc Lyα halo around a high-redshift type-2 quasi-stellar object

We announce the discovery of an extended emission-line region associated with a high-redshift type-2 quasi-stellar object (QSO). The halo, which was discovered in our new wide-field narrow-band survey, resides at   z = 2.85  in the Spitzer First Look Survey region and is extended over ∼80 kpc. Deep very long baseline interferometry (VLBI) observations imply that approximately 50 per cent of the radio emission is extended on scales >200 pc. The inferred active galactic nuclei (AGN) luminosity is sufficient to ionize the extended halo, and the optical emission is consistent with being triggered coevally with the radio source. The Lyα halo is as luminous as those found around high-redshift radio galaxies; however, the active nucleus is several orders of magnitude less luminous at radio wavelengths than those Fanarof–Riley type II (FRIIs) more commonly associated with extended emission-line regions. AMS05 appears to be a high-redshift analogue to the radio-quiet quasar E1821+643 which is core dominated, but which also exhibits extended Fanarof–Riley type I (FRI)-like structure and contains an optically powerful AGN. We also find evidence for more quiescent kinematics in the Lyα emission line in the outer regions of the halo, reminiscent of the haloes around the more powerful FRIIs. The optical to mid-infrared spectral energy distribution is well described by a combination of an obscured QSO  ( L _bol∼ 3.4 ± 0.2 × 10¹³ L_⊙)  and a 1.4 Gyr old simple stellar population with mass  ∼3.9 ± 0.3 × 10¹¹ M_⊙  .     

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.