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.     

Cluster strong lensing in the Millennium simulation: the effect of galaxies and structures along the line-of-sight

We use ray-tracing through the Millennium simulation to study how secondary matter structures along the line-of-sight and the stellar mass in galaxies affect strong cluster lensing, in particular the cross-section for giant arcs. Furthermore, we investigate the distribution of the cluster Einstein radii and the radial distribution of giant arcs. We find that additional structures along the line-of-sight increase the strong-lensing optical depth by  ∼10–25 per cent  , while strong-lensing cross-sections of individual clusters are frequently boosted by as much as  ∼50 per cent  . The enhancement is mainly due to structures that are not correlated with the lens. Cluster galaxies increase the strong-lensing optical depth by up to a factor of 2, while interloping galaxies are not significant. We conclude that these effects need to be taken into account for predictions of the giant arc abundance, but they are not large enough to fully account for the reported discrepancy between predicted and observed abundances.
Furthermore, we find that Einstein radii defined via the area enclosed by the critical curve are 10–30 per cent larger than those defined via radial surface mass density profiles. The contributions of radial and tangential arcs to the radial distribution of arcs can be clearly distinguished. The radial distribution of tangential arcs is very broad and extends out to several Einstein radii. Thus, individual arcs are not well suited for constraining Einstein radii.     

Galaxy density profiles and shapes – I. Simulation pipeline for lensing by realistic galaxy models

Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, which we aim to investigate in a consistent way. In this first paper in a series, we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii.     

Galaxy density profiles and shapes – II. Selection biases in strong lensing surveys

Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio and image separation distribution, with a flexible treatment of magnification bias to mimic different survey strategies. We present our results for two families of density profiles: cusped and deprojected Sérsic models. While we use fixed lens and source redshifts for most of the analysis, we show that the results are applicable to other redshift combinations, and we also explore the physics of how our results change for very different redshifts. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sérsic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.     

Strong lensing optical depths in a ΛCDM universe

We investigate strong gravitational lensing in the concordance ΛCDM cosmology by carrying out ray tracing along past light cones through the Millennium Simulation, the largest simulation of cosmic structure formation ever carried out. We extend previous ray-tracing methods in order to take full advantage of the large volume and the excellent spatial and mass resolution of the simulation. As a function of source redshift we evaluate the probability that an image will be highly magnified, will be highly elongated or will be one of a set of multiple images. We show that such strong lensing events can almost always be traced to a single dominant lensing object and we study the mass and redshift distribution of these primary lenses. We fit analytic models to the simulated dark haloes in order to study how our optical depth measurements are affected by the limited resolution of the simulation and of the lensing planes that we construct from it. We conclude that such effects lead us to underestimate total strong lensing cross-sections by about 15 per cent. This is smaller than the effects expected from our neglect of the baryonic components of galaxies. Finally we investigate whether strong lensing is enhanced by material in front of or behind the primary lens. Although strong lensing lines of sight are indeed biased towards higher than average mean densities, this additional matter typically contributes only a few per cent of the total surface density.     

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.     

The formation and survival of discs in a ΛCDM universe

We study the formation of galaxies in a Λ cold dark matter (ΛCDM) universe using high-resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight isolated haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at   z = 0  , none of our galaxies contains a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The   z = 0  spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at   z = 0  nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at   z ≳ 2  , regardless of their   z = 0  morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic ΛCDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.     

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.     

X-ray and strong lensing mass estimate of MS2137.3–2353

We present new mass estimates of the cluster of galaxies MS2137.3–2353, inferred from X-ray and strong lensing analyses. This cluster exhibits an outstanding strong lensing configuration and indicates a well-relaxed dynamical state, being most suitable for a mass reconstruction which combines both techniques. Despite this, several previous studies have claimed a significant discrepancy between the X-ray and the strong lensing mass estimates. The primary aim of this paper is to address and explain this mismatch. For this purpose, we have analysed Chandra observations to recover the profiles of the intracluster medium properties and, assuming a functional form for the matter density, the total mass distribution. The notable strong-lensing features of MS2137.3 allow us to reconstruct its projected mass in the central regions with good accuracy, by taking advantage of the lensing inversion code lenstool . We compare the results obtained for both methods. Our mass estimates for MS2137.3 are in agreement within errors, leading to a mean, extrapolated value of   M ₂₀₀≃ 4.4 ± 0.3 × 10¹⁴ M_⊙  , under the assumption of the Navarro–Frenk–White (NFW) mass profile. However, the strong lensing mass estimate is affected by the details of the brightest cluster galaxy mass modelling, since the radial arc is a very sensitive probe of the total mass derivative in the central region. In particular, we do not find evidence for a high concentration for the NFW density profile, as reported in some earlier works.     

Constraining the Mass Distribution of Cluster Galaxies by Weak Lensing

Analysing the weak lensing distortions of the images of faint background galaxies provides a means to constrain the average mass distribution of cluster galaxies and potentially to test the extent of their dark matter haloes as a function of the density of their environment. The observable image distortions are a consequence of the interplay between the effects of a global cluster mass distribution and the perturbations resulting from individual cluster galaxies. Starting from a reconstruction of the cluster mass distribution with conventional techniques, we apply a maximum likelihood method to infer the average properties of an ensemble of cluster galaxies. From simulations this approach is found to be reliable as long as the galaxies including their dark matter haloes only contribute a small fraction to the total mass of the system. If their haloes are extended, the galaxies contain a substantial mass fraction. In this case our method is still applicable in the outer regions of clusters, where the surface mass density is low, but yields biased estimates of the parameters describing the mass profiles of the cluster galaxies in the central part of the cluster. In that case it will be necessary to resort to more sophisticated strategies by modelling cluster galaxies and an underlying global mass distribution simultaneously. We conclude that galaxy–galaxy lensing in clusters provides a unique means to probe the presence and extent of dark haloes of cluster galaxies.     

The dark matter distribution in disc galaxies

We use high-quality optical rotation curves of nine low-luminosity disc galaxies to obtain the velocity profiles of the surrounding dark matter haloes. We find that they increase linearly with radius at least out to the edge of the stellar disc, implying that, over the entire stellar region, the density of the dark halo is about constant.
The properties of the mass structure of these haloes are similar to those found for a number of dwarf and low surface brightness galaxies, but provide a more substantial evidence of the discrepancy between the halo mass distribution predicted in the cold dark matter scenario and those actually detected around galaxies. We find that the density law proposed by Burkert reproduces the halo rotation curves, with halo central densities ( ρ ₀∼1–4×10⁻²⁴ g cm⁻³) and core radii ( r ₀∼5–15 kpc) scaling as ρ ₀∝ r ₀^−2/3.     

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.     

Dwarf galaxy rotation curves and the core problem of dark matter haloes

The standard cold dark matter (CDM) model has recently been challenged by the claim that dwarf galaxies have dark matter haloes with constant-density cores, whereas CDM predicts haloes with steeply cusped density distributions. Consequently, numerous alternative dark matter candidates have recently been proposed. In this paper we scrutinize the observational evidence for the incongruity between dwarf galaxies and the CDM model. To this end, we analyse the rotation curves of 20 late-type dwarf galaxies studied by Swaters. Taking the effects of beam smearing and adiabatic contraction into account, we fit mass models to these rotation curves with dark matter haloes with different cusp slopes, ranging from constant-density cores to r ⁻² cusps. Even though the effects of beam smearing are small for these data, the uncertainties in the stellar mass-to-light ratio and the limited spatial sampling of the halo's density distribution hamper a unique mass decomposition. Consequently, the rotation curves in our sample cannot be used to discriminate between dark haloes with constant-density cores and r ⁻¹ cusps. We show that the dwarf galaxies analysed here are consistent with CDM haloes in a ΛCDM cosmology, and that there is thus no need to abandon the idea that dark matter is cold and collisionless. However, the data are also consistent with any alternative dark matter model that produces dark matter haloes with central cusps less steep than r ^−1.5. In fact, we argue that based on existing H  i rotation curves alone, at best weak limits can be obtained on cosmological parameters and/or the nature of the dark matter. In order to make progress, rotation curves with higher spatial resolution and independent measurements of the mass-to-light ratio of the disc are required.     

Galaxy star formation in different environments

We use a semi-analytic model of galaxy formation to study signatures of large-scale modulations in the star formation (SF) activity in galaxies. In order to do this, we carefully define local and global estimators of the density around galaxies. The former are computed using a voronoi tessellation technique and the latter are parametrized by the normalized distance to haloes and voids, in terms of the virial and void radii, respectively. As a function of local density, galaxies show a strong modulation in their SF, a result that is in agreement with those from several authors. When taking subsamples of equal local density at different large-scale environments, we find relevant global effects whereby the fraction of red galaxies diminishes for galaxies in equal local density environments farther away from clusters and closer to voids. In general, the semi-analytic simulation is in good agreement with the available observational results, and offers the possibility to disentangle many of the processes responsible for the variation of galaxy properties with the environment; we find that the changes found in samples of galaxies with equal local environment but different distances to haloes or voids come from the variations in the underlying mass function of dark matter (DM) haloes. There is an additional possible effect coming from the host DM halo ages, indicating that halo assembly also plays a small but significant role (1.14σ) in shaping the properties of galaxies, and in particular, hints at a possible spatial correlation in halo/stellar mass ages. An interesting result comes from the analysis of the coherence of flows in different large-scale environments of fixed local densities; the neighbourhoods of massive haloes are characterized by lower coherences than control samples, except for galaxies in filament-like regions, which show highly coherent motions.     

Probing galactic dark matter in dense environments: on the strong lensing efficiency of galaxies in rich clusters

The recent detection by Limousin et al. of five new strong lensing events dominated by galaxy cluster members in Abell 1689, and outside the critical regime of the cluster itself, offers a way to obtain constraints on the cluster mass distribution in a region inaccessible to standard lensing analysis. In addition, modelling such systems will provide another window on the dark matter haloes of galaxies in very dense environments. Here, it is shown that the boost in image separation due to the external shear and convergence from a smooth cluster component means that more numerous, less massive galaxies have the potential to create multiple images with detectable separations, relative to isolated field galaxies. This comes in addition to a potential increase in their lensing (source plane) cross-section. To gain insight into the factors involved and as a precursor to a numerical study using N -body simulations, a simple analytic model of a cluster at   z = 0.3  lensing background galaxies at   z = 2  is considered here. The fiducial model has cluster members with isothermal density profiles and luminosities L , distributed in a Schechter function (faint-end slope  ν=−1.25  ), related to their velocity dispersions σ via the Faber–Jackson scaling L ∝σ⁴. Just outside the critical regime of the cluster, the scale of galaxy-dominated image separations is significantly increased. Folding in the fact that less massive galaxies present a lower lensing cross-section, and that the cross-section can itself be enhanced in an external field leads to a factor of a few times more detected events relative to field galaxies. These values will be higher closer to the critical curve. Given that the events in Abell 1689 were detected over a very small region of the cluster where ACS data were available, this motivates the search for such events in other clusters.     

Dissipationless mergers of elliptical galaxies and the evolution of the fundamental plane

We carry out numerical simulations of dissipationless major mergers of elliptical galaxies using initial galaxy models that consist of a dark matter haloes and a stellar bulge with properties consistent with the observed fundamental plane. By varying the density profile of the dark matter haloes [standard Navarro, Frenk & White (NFW) profile versus adiabatically contracted NFW profile], the global stellar to dark matter mass ratio and the orbit of the merging galaxies, we are able to assess the impact of each of these factors on the structure of the merger remnant. Our results indicate that the properties of the remnant bulge depend primarily on the angular momentum and energy of the orbit; for a cosmologically motivated orbit, the effective radius and velocity dispersion of the remnant bulge remain approximately on the fundamental plane. This indicates that the observed properties of elliptical galaxies are consistent with significant growth via late dissipationless mergers. We also find that the dark matter fraction within the effective radius of our remnants increases after the merger, consistent with the hypothesis that the tilt of the fundamental plane from the virial theorem is due to a varying dark matter fraction as a function of galaxy mass.     

Understanding the halo-mass and galaxy-mass cross-correlation functions

We use the Millennium Simulation (MS) to measure the cross-correlation between halo centres and mass (or equivalently the average density profiles of dark haloes) in a Lambda cold dark matter (ΛCDM) cosmology. We present results for radii in the range  10  h ⁻¹ kpc < r < 30  h ⁻¹ Mpc  and for halo masses in the range  4 × 10¹⁰ < M ₂₀₀ < 4 × 10¹⁴  h ⁻¹ M_⊙  . Both at   z = 0  and at   z = 0.76  these cross-correlations are surprisingly well fitted if the inner region is approximated by a density profile of NFW or Einasto form, the outer region by a biased version of the linear mass autocorrelation function, and the maximum of the two is adopted where they are comparable. We use a simulation of galaxy formation within the MS to explore how these results are reflected in cross-correlations between galaxies and mass. These are directly observable through galaxy–galaxy lensing. Here also we find that simple models can represent the simulation results remarkably well, typically to ≲10 per cent. Such models can be used to extend our results to other redshifts, to cosmologies with other parameters, and to other assumptions about how galaxies populate dark haloes. Our galaxy formation simulation already reproduces current galaxy–galaxy lensing data quite well. The characteristic features predicted in the galaxy–galaxy lensing signal should provide a strong test of the ΛCDM cosmology as well as a route to understanding how galaxies form within it.     

Numerical modelling of the vertical structure and dark halo parameters in disc galaxies

The non‐linear dynamics of bending instability and vertical structure of a galactic stellar disc embedded into a spherical halo are studied with N‐body numerical modelling. Development of the bending instability in stellar galactic disc is considered as the main factor that increases the disc thickness. Correlation between the disc vertical scale height and the halo‐to‐disc mass ratio is predicted from the simulations. The method of assessment of the spherical‐to‐disc mass ratio for edge‐on spiral galaxies with a small bulge is considered. Modelling of eight edge‐on galaxies: NGC 891, NGC 4738, NGC 5170, UGC 6080, UGC 7321, UGC 8286, UGC 9422 and UGC 9556 is performed. Parameters of stellar discs, dark haloes and bulges are estimated. The lower limit of the dark‐to‐luminous mass ratio in our galaxies is of the order of one within the limits of their stellar discs. The dark haloes dominate by mass in the galaxies with very thin stellar discs (NGC 5170, UGC 7321 and UGC 8286) (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Measuring the distribution of galaxies between haloes

We develop a method to measure the probability, P ( N;   M ), of finding N galaxies in a dark matter halo of mass M from the theoretically determined clustering properties of dark matter haloes and the observationally measured clustering properties of galaxies. Knowledge of this function and the distribution of the dark matter completely specifies all clustering properties of galaxies on scales larger than the size of dark matter haloes. Furthermore, P ( N;   M ) provides strong constraints on models of galaxy formation, since it depends upon the merger history of dark matter haloes and the galaxy–galaxy merger rate within haloes. We show that measurements from a combination of the Two Micron All Sky Survey and Sloan Digital Sky Survey or Two-degree Field Galaxy Redshift Survey data sets will allow P ( N;   M ) averaged over haloes occupied by bright galaxies to be accurately measured for N =0–2 .