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  .     

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.     

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.     

Clustering Property of Wolf-Rayet Galaxies in the SDSS

We have analysed, for the first time, the clustering properties of Wolf-Rayet (W-R) galaxies, using a large sample of 846 W-R galaxies selected from the Data Release 4 (DR4) of the Sloan Digital Sky Survey (SDSS). We compute the cross-correlation function between W- R galaxies and a reference sample of galaxies drawn from the DR4. We compare the function to the results for control samples of non-W-R star-forming galaxies that are matched closely in redshift, luminosity, concentration, 4000-A break strength and specific star formation rate (SSFR). On scales larger than a few Mpc, W-R galaxies have almost the same clustering amplitude as the control samples, indicating that W-R galaxies and non-W-R control galax- ies populate dark matter haloes of similar masses. On scales between 0.1-1 h-1 Mpc, W-R galaxies are less clustered than the control samples, and the size of the difference depends on the SSFR. Based on both observational and theoretical considerations, we speculate that this negative bias can be interpreted by W-R galaxies residing preferentially at the centers of their dark matter haloes. We examine the distribution of W-R galaxies more closely using the SDSS galaxy group catalogue of Yang et al., and find that ～82% of our W-R galaxies are the central galaxies of groups, compared to ～74% for the corresponding control galaxies. We find that W-R galaxies are hosted, on average, by dark matter haloes of masses of 1012,3 M☉, compared to 1012,1 M? For centrally-located W-R galaxies and 1012,7 M☉ For satellite ones. We would like to point out that this finding, which provides a direct observational support to our conjecture, is really very crude due to the small number of W-R galaxies and the incom- pleteness of the group catalogue, and needs more work in future with larger samples.     

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.     

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.     

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.     

Spatial clustering of Ultra Steep Spectrum sources and galaxies

We present measurements of the clustering properties of galaxies in the field of redshift range 0.5 ≲ z ≲ 1.5 Ultra Steep Spectrum radio sources selected from the Sydney University Molonglo Sky Survey and the National Radio Astronomy Observatories Very Large Array Sky Survey. Galaxies in these USS fields were identified in deep near-infrared observations, complete down to   K _s= 20  , using the IRIS2 instrument at the Anglo-Australian Telescope. We used the redshift distribution of   K _s < 20  galaxies taken from Cimatti et al. (2002) to constrain the correlation length r ₀. We find a strong correlation signal of galaxies with   K _s < 20  around our USS sample. A comoving correlation length   r ₀= 14.0 ± 2.8  h ⁻¹ Mpc  and γ= 1.98 ± 0.15 are derived in a flat cosmological model universe.
We compare our findings with those obtained in a cosmological N -body simulation populated with galform semi-analytic galaxies. We find that clusters of galaxies with masses in the range   M = 10^13.4–14.2  h ⁻¹ M_⊙  have a cluster–galaxy cross-correlation amplitude comparable to those found between the USS hosts and galaxies. These results suggest that distant radio galaxies are excellent tracers of galaxy overdensities and pinpoint the progenitors of present day rich clusters of galaxies.     

The three-dimensional power spectrum of dark and luminous matter from the VIRMOS-DESCART cosmic shear survey

We present the first optimal power spectrum estimation and three-dimensional deprojections for the dark and luminous matter and their cross-correlations. The results are obtained using a new optimal fast estimator, deprojected using minimum variance and Singular Value Decomposition (SVD) techniques. We show the resulting 3D power spectra for dark matter and galaxies, and their covariance for the VIRMOS-DESCART weak lensing shear and galaxy data. The survey is most sensitive to non-linear scales   k _NL∼ 1 h Mpc⁻¹  . On these scales, our 3D power spectrum of dark matter is in good agreement with the RCS 3D power spectrum found by Tegmark & Zaldarriaga. Our galaxy power is similar to that found by the 2MASS survey, and larger than that of SDSS, APM and RCS, consistent with the expected difference in galaxy population.
We find an average bias   b = 1.24 ± 0.18  for the I -selected galaxies, and a cross-correlation coefficient   r = 0.75 ± 0.23  . Together with the power spectra, these results optimally encode the entire two point information about dark matter and galaxies, including galaxy–galaxy lensing. We address some of the implications regarding galaxy haloes and mass-to-light ratios. The best-fitting 'halo' parameter   h ≡ r / b = 0.57 ± 0.16  , suggesting that dynamical masses estimated using galaxies systematically underestimate total mass.
Ongoing surveys, such as the Canada–France–Hawaii Telescope Legacy Survey, will significantly improve on the dynamic range, and future photometric redshift catalogues will allow tomography along the same principles.     

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

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

The assembly bias of dark matter haloes to higher orders

We use an extremely large volume  (2.4  h ⁻³ Gpc³)  , high-resolution N -body simulation to measure the higher order clustering of dark matter haloes as a function of mass and internal structure. As a result of the large simulation volume and the use of a novel 'cross-moment' counts-in-cells technique which suppresses discreteness noise, we are able to measure the clustering of haloes corresponding to rarer peaks than was possible in previous studies; the rarest haloes for which we measure the variance are 100 times more clustered than the dark matter. We are able to extract, for the first time, halo bias parameters from linear up to fourth order. For all orders measured, we find that the bias parameters are a strong function of mass for haloes more massive than the characteristic mass   M _*  . Currently, no theoretical model is able to reproduce this mass dependence closely. We find that the bias parameters also depend on the internal structure of the halo up to fourth order. For haloes more massive than   M _*  , we find that the more concentrated haloes are more weakly clustered than the less concentrated ones. We see no dependence of clustering on concentration for haloes with masses   M < M _*  ; this is contrary to the trend reported in the literature when segregating haloes by their formation time. Our results are insensitive to whether haloes are labelled by the total mass returned by the friends-of-friends group finder or by the mass of the most massive substructure. This implies that our conclusions are not an artefact of the particular choice of group finding algorithm. Our results will provide important input to theoretical models of galaxy clustering.     

Empirical distributions of galactic λ spin parameters from the Sloan Digital Sky Survey

Using simple dimensional arguments for both spiral and elliptical galaxies, we present formulae to derive an estimate of the halo spin parameter λ for any real galaxy, in terms of common observational parameters. This allows a rough estimate of λ, which we apply to a large volume-limited sample of galaxies taken from the Sloan Digital Sky Survey data base. The large numbers involved (11 597) allow the derivation of reliable λ distributions, as signal adds up significantly in spite of the errors in the inferences for particular galaxies. We find that if the observed distribution of λ is modelled with a lognormal function, as often done for this distribution in dark matter haloes that appear in cosmological simulations, we obtain parameters  λ₀= 0.04 ± 0.005  and  σ_λ= 0.51 ± 0.05  , interestingly consistent with values derived from simulations. For spirals, we find a good correlation between empirical values of λ and visually assigned Hubble types, highlighting the potential of this physical parameter as an objective classification tool.     

Forming a large disc galaxy from a z < 1 major merger

Using high-resolution SPH simulations in a fully cosmological Λ cold dark matter context, we study the formation of a bright disc-dominated galaxy that originates from a 'wet' major merger at   z = 0.8  . The progenitors of the disc galaxy are themselves disc galaxies that formed from early major mergers between galaxies with blue colours. A substantial thin stellar disc grows rapidly following the last major merger and the present-day properties of the final remnant are typical of early-type spiral galaxies, with an i -band bulge-to-disc ratio ∼0.65, a disc scalelength of 7.2 kpc,   g − r = 0.5 mag  , an H  i linewidth ( W ₂₀/2) of 238 km s⁻¹ and total magnitude   i =−22.4  . The key ingredients for the formation of a dominant stellar disc component after a major merger are (i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, (ii) supernova feedback that is able to partially suppress star formation during mergers and (iii) relative fading of the spheroidal component. The gas fraction of the progenitors' discs does not exceed 25 per cent at   z < 3  , emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of discs and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disc galaxies by the high likelihood of interactions and mergers for galaxy-sized haloes at relatively low z .     

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.     

Completely dark galaxies: their existence, properties and strategies for finding them

There are a number of theoretical and observational hints that large numbers of low-mass galaxies composed entirely of dark matter exist in the field. The theoretical considerations follow from the prediction of cold dark matter theory that there exist many low-mass galaxies for every massive one. The observational considerations follow from the observed paucity of these low-mass galaxies in the field but not in dense clusters of galaxies; this suggests that the lack of small galaxies in the field is due to the inhibition of star formation in the galaxies as opposed to the fact that their small dark matter haloes do not exist. In this work we outline the likely properties of low-mass dark galaxies, and describe observational strategies for finding them, and where in the sky to search. The results are presented as a function of the global properties of dark matter, in particular the presence or absence of a substantial baryonic dark matter component. If the dark matter is purely cold and has a Navarro, Frenk & White density profile, directly detecting dark galaxies will only be feasible with present technology if the galaxy has a maximum velocity dispersion in excess of 70 km s⁻¹, in which case the dark galaxies could strongly lens background objects. This is much higher than the maximum velocity dispersions in most dwarf galaxies. If the dark matter in galaxy haloes has a baryonic component close to the cosmic ratio, the possibility of directly detecting dark galaxies is much more realistic; the optimal method of detection will depend on the nature of the dark matter. A number of more indirect methods are also discussed.