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.     

Collisional baryonic dark matter haloes

If dark haloes are composed of dense gas clouds, as has recently been inferred, then collisions between clouds lead to galaxy evolution. Collisions introduce a core in an initially singular dark matter distribution, and can thus help to reconcile scale-free initial conditions – such as are found in simulations – with observed haloes, which have cores. A pseudo-Tully–Fisher relation, between halo circular speed and visible mass (not luminosity), emerges naturally from the model: M _vis∝ V ^7/2.
Published data conform astonishingly well to this theoretical prediction. For our sample of galaxies, the mass–velocity relationship has much less scatter than the Tully–Fisher relation, and holds as well for dwarf galaxies (where diffuse gas makes a sizeable contribution to the total visible mass) as it does for giants. It seems very likely that this visible-mass/velocity relationship is the underlying physical basis for the Tully–Fisher relation, and this discovery in turn suggests that the dark matter is both baryonic and collisional.     

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.     

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 image separation distribution of strong lenses: halo versus subhalo populations

We present a halo model prediction of the image separation distribution of strong lenses. Our model takes into account the subhalo population, which has been ignored in previous studies, as well as the conventional halo population. Haloes and subhaloes are linked to central and satellite galaxies by adopting a universal scaling relation between masses of (sub)haloes and luminosities of galaxies. Our model predicts that 10–20 per cent of lenses should be caused by the subhalo population. The fraction of lensing by satellite galaxies (subhaloes) peaks at ∼1 arcsec and decreases rapidly with increasing image separations. We compute fractions of lenses which lie in groups and clusters and find them to be ∼14 and ∼4 per cent, respectively; nearly half of such lenses are expected to be produced by satellite galaxies, rather than central parts of haloes. We also study mass distributions of lensing haloes and find that, even at image separations of ∼3 arcsec, the deviation of lens mass distributions from isothermal profiles is large; at or beyond ∼3 arcsec, image separations are enhanced significantly by surrounding haloes. Our model prediction agrees reasonably well with observed image separation distributions from galaxy to cluster scales.     

Massive dark matter haloes around bright isolated galaxies in the 2dFGRS

We identify a large sample of isolated bright galaxies and their fainter satellites in the 2dF Galaxy Redshift Survey (2dFGRS). We analyse the dynamics of ensembles of these galaxies selected according to luminosity and morphological type by stacking the positions of their satellites and estimating the velocity dispersion of the combined set. We test our methodology using realistic mock catalogues constructed from cosmological simulations. The method returns an unbiased estimate of the velocity dispersion provided that the isolation criterion is strict enough to avoid contamination and that the scatter in halo mass at fixed primary luminosity is small. Using a maximum likelihood estimator that accounts for interlopers, we determine the satellite velocity dispersion within a projected radius of 175  h ⁻¹ kpc. The dispersion increases with the luminosity of the primary and is larger for elliptical galaxies than for spiral galaxies of similar b _J luminosity. Calibrating the mass–velocity dispersion relation using our mock catalogues, we find a dynamical mass within 175  h ⁻¹ kpc of     for elliptical galaxies and     for spiral galaxies. Finally, we compare our results with recent studies and investigate their limitations using our mock catalogues.     

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.     

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.     

GHASP: an Hα kinematic survey of spiral and irregular galaxies – V. Dark matter distribution in 36 nearby spiral galaxies

The results obtained from a study of the mass distribution of 36 spiral galaxies are presented. The galaxies were observed using Fabry–Perot interferometry as part of the GHASP survey. The main aim of obtaining high-resolution Hα 2D velocity fields is to define more accurately the rising part of the rotation curves which should allow to better constrain the parameters of the mass distribution. The Hα velocities were combined with low resolution H  i data from the literature, when available. Combining the kinematical data with photometric data, mass models were derived from these rotation curves using two different functional forms for the halo: an isothermal sphere (ISO) and a Navarro–Frenk–White (NFW) profile. For the galaxies already modelled by other authors, the results tend to agree. Our results point at the existence of a constant density core in the centre of the dark matter haloes rather than a cuspy core, whatever the type of the galaxy from Sab to Im. This extends to all types the result already obtained by other authors studying dwarf and low surface brightness galaxies but would necessitate a larger sample of galaxies to conclude more strongly. Whatever model is used (ISO or NFW), small core radius haloes have higher central densities, again for all morphological types. We confirm different halo scaling laws, such as the correlations between the core radius and the central density of the halo with the absolute magnitude of a galaxy: low-luminosity galaxies have small core radius and high central density. We find that the product of the central density with the core radius of the dark matter halo is nearly constant, whatever the model and whatever the absolute magnitude of the galaxy. This suggests that the halo surface density is independent from the galaxy type.     

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.     

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.     

Testing the modified Press–Schechter model against N-body simulations

A modified version of the extended Press–Schechter model for the growth of dark-matter haloes was introduced in two previous papers, with the aim of explaining the mass–density relation shown by haloes in high-resolution cosmological simulations. In this model, major mergers are well separated from accretion, thereby allowing a natural definition of halo formation and destruction. This makes it possible to derive analytic expressions for halo formation and destruction rates, the mass accretion rate and the probability distribution functions of halo formation times and progenitor masses. The stochastic merger histories of haloes can be readily derived and easily incorporated into semi-analytical models of galaxy formation, thus avoiding the usual problems encountered in the construction of Monte Carlo merger trees from the original extended Press–Schechter formalism. Here we show that the predictions of the modified Press–Schechter model are in good agreement with the results of N -body simulations for several scale-free cosmologies.     

The haloes of merger remnants

We perform collisionless N -body simulations of 1:1 galaxy mergers, using models which include a galaxy halo, disc and bulge, focusing on the behaviour of the halo component. The galaxy models are constructed without recourse to a Maxwellian approximation. We investigate the effect of varying the galaxies' orientation, their mutual orbit and the initial velocity anisotropy or cusp strength of the haloes upon the remnant halo density profiles and shape, as well as on the kinematics. We observe that the halo density profile (determined as a spherical average, an approximation we find appropriate) is exceptionally robust in mergers, and that the velocity anisotropy of our remnant haloes is nearly independent of the orbits or initial anisotropy of the haloes. The remnants follow the halo anisotropy – local density slope (β–γ) relation suggested by Hansen & Moore in the inner parts of the halo, but β is systematically lower than this relation predicts in the outer parts. Remnant halo axis ratios are strongly dependent on the initial parameters of the haloes and on their orbits. We also find that the remnant haloes are significantly less spherical than those described in studies of simulations which include gas cooling.     

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  .     

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.     

Dynamical friction and galaxy merging time-scales

The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N -body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M _sat/ M _host≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M _sat/ M _host≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halo's centre for a wide range of M _sat/ M _host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≲10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.     

The correlation of star formation quenching with internal galaxy properties and environment

We investigate the correlation of star formation quenching with internal galaxy properties and large-scale environment (halo mass) in empirical data and theoretical models. We make use of the halo-based group catalogue of Yang and collaborators, which is based on the Sloan Digital Sky Survey. Data from the Galaxy evolution explorer are also used to extract the recent star formation rate. In order to investigate the environmental effects, we examine the properties of 'central' and 'satellite' galaxies separately. For central galaxies, we are unable to conclude whether star formation quenching is primarily connected with halo mass or stellar mass, because these two quantities are themselves strongly correlated. For satellite galaxies, a nearly equally strong dependence on halo mass and stellar mass is seen. We make the same comparison for five different semi-analytic models based on three independently developed codes. We find that the models with active galactic nuclei feedback reproduce reasonably well the dependence of the fraction of central red and passive galaxies on halo mass and stellar mass. However, for satellite galaxies, the same models badly overproduce the fraction of red/passive galaxies and do not reproduce the empirical trends with stellar mass or halo mass. This satellite overquenching problem is caused by the too-rapid stripping of the satellites' hot gas haloes, which leads to rapid strangulation of star formation.     

The X-ray emission in post-merger ellipticals

The evolution in X-ray properties of early-type galaxies is largely unconstrained. In particular, little is known about how, and if, remnants of mergers generate hot gas haloes. Here we examine the relationship between X-ray luminosity and galaxy age for a sample of early-type galaxies. Comparing normalized X-ray luminosity to three different age indicators, we find that L _X L _B increases with age, suggesting an increase in X-ray halo mass with time after the last major star formation episode of a galaxy. The long-term nature of this trend, which appears to continue across the full age range of our sample, poses a challenge for many models of hot halo formation. We conclude that models involving a declining rate of type Ia supernovae, and a transition from outflow to inflow of the gas originally lost by galactic stars, offer the most promising explanation for the observed evolution in X-ray luminosity.