Lopsided and Bulging Distribution of Satellites around Paired Halos. II. 3D Analysis and Dependence on Projection and Selection Effects

ABSTRACTWe present a method to flexibly and self-consistently determine individual galaxies’ star formation rates (SFRs) from their host haloes’ potential well depths, assembly histories, and redshifts. The method is constrained by galaxies’ observed stellar mass functions, SFRs (specific and cosmic), quenched fractions, ultraviolet (UV) luminosity functions, UV–stellar mass relations, IRX–UV relations, auto- and cross-correlation functions (including quenched and star-forming subsamples), and quenching dependence on environment; each observable is reproduced over the full redshift range available, up to 0 < z < 10. Key findings include the following: galaxy assembly correlates strongly with halo assembly; quenching correlates strongly with halo mass; quenched fractions at fixed halo mass decrease with increasing redshift; massive quenched galaxies reside in higher-mass haloes than star-forming galaxies at fixed galaxy mass; star-forming and quenched galaxies’ star formation histories at fixed mass differ most at z < 0.5; satellites have large scatter in quenching time-scales after infall, and have modestly higher quenched fractions than central galaxies; Planck cosmologies result in up to 0.3 dex lower stellar – halo mass ratios at early times; and, none the less, stellar mass–halo mass ratios rise at z > 5. Also presented are revised stellar mass – halo mass relations for all, quenched, star-forming, central, and satellite galaxies; the dependence of star formation histories on halo mass, stellar mass, and galaxy SSFR; quenched fractions and quenching time-scale distributions for satellites; and predictions for higher-redshift galaxy correlation functions and weak lensing surface densities. The public data release (DR1) includes the massively parallel (>105 cores) implementation (the UniverseMachine), the newly compiled and remeasured observational data, derived galaxy formation constraints, and mock catalogues including lightcones.

We describe the luminosity function, based on Sersic fits to the light profiles, of CMASS galaxies at z ~ 0.55. Compared to previous estimates, our Sersic-based reductions imply more luminous, massive galaxies, consistent with the effects of Sersic- rather than Petrosian or de Vaucouleur-based photometry on the Sloan Digital Sky Survey (SDSS) main galaxy sample at z ~ 0.1. This implies a significant revision of the high mass end of the correlation between stellar and halo mass. Inferences about the evolution of the luminosity and stellar mass functions depend strongly on the assumed, and uncertain, k+e corrections. In turn, these depend on the assumed age of the population. Applying k+e corrections taken from fitting the models of Maraston et al. (2009) to the colors of both SDSS and CMASS galaxies, the evolution of the luminosity and stellar mass functions appears impressively passive, provided that the fits are required to return old ages. However, when matched in comoving number- or luminosity-density, the SDSS galaxies are less strongly clustered compared to their counterparts in CMASS. This rules out the passive evolution scenario, and, indeed, any minor merger scenarios which preserve the rank ordering in stellar mass of the population. Potential incompletenesses in the CMASS sample would further enhance this mismatch. Our analysis highlights the virtue of combining clustering measurements with number counts.

We investigate the quenching properties of central and satellite galaxies, utilizing the halo masses and central–satellite identifications from the Sloan Digital Sky Survey galaxy group catalog of Yang et al. We find that the quenched fractions of centrals and satellites of similar stellar masses have similar dependence on host halo mass. The similarity of the two populations is also found in terms of specific star formation rate and 4000 Å break. The quenched fractions of centrals and satellites of similar masses show similar dependencies on bulge-to-total light ratio, central velocity dispersion, and halo-centric distance in halos of given halo masses. The prevalence of optical/radio-loud active galactic nuclei is found to be similar for centrals and satellites at given stellar masses. All these findings strongly suggest that centrals and satellites of similar masses experience similar quenching processes in their host halos. We discuss implications of our results for the understanding of galaxy quenching.

We analyse galaxies in groups in the Sloan Digital Sky Survey (SDSS) and find a weak but significant assembly-type bias, where old central galaxies have a higher clustering amplitude (61 $\pm$ 9 per cent) at scales > 1 Mpc than young central galaxies of equal host halo mass ($M_{h} \sim 10^{11.8} h^{-1}$ $M_{\odot}$). The observational sample is volume-limited out to z=0.1 with $M_r -$ 5 log$(h) \le -19.6$. We construct a mock catalogue of galaxies that shows a similar signal of assembly bias (46 $\pm$ 9 per cent) at the same halo mass. We then adapt the model presented by Lacerna & Padilla (Paper I) to redefine the overdensity peak height, which traces the assembly bias such that galaxies in equal density peaks show the same clustering regardless of their stellar age, but this time using observational features such as a flux limit. The proxy for peak height, which is proposed as a better alternative than the virial mass, consists in the total mass given by the mass of neighbour host haloes in cylinders centred at each central galaxy. The radius of the cylinder is parametrized as a function of stellar age and virial mass. The best-fitting set of parameters that make the assembly bias signal lower than 5$-$15 per cent for both SDSS and mock central galaxies are similar. The idea behind the parametrization is not to minimize the bias, but it is to use this method to understand the physical features that produce the assembly bias effect. Even though the tracers of the density field used here differ significantly from those used in paper I, our analysis of the simulated catalogue indicates that the different tracers produce correlated proxies, and therefore the reason behind this assembly bias is the crowding of peaks in both simulations and the SDSS.

We explore the feasibility of learning the connection between Sloan Digital Sky Survey (SDSS) galaxies and ELUCID subhaloes with random forest (RF). ELUCID is a constrained N-body simulation constructed using the matter density field of SDSS. Based on a SDSS-ELUCID matched catalogue, we build RF models that predict Mr magnitude, colour, stellar mass M*, and specific star formation rate (sSFR) with several subhalo properties. While the RF can predict Mr and M* with reasonable accuracy, the prediction accuracy of colour and sSFR is low, which could be due to the mismatch between galaxies and subhaloes. To test this, we shuffle the galaxies in subhaloes of narrow mass bins in the local neighbourhood using galaxies of a semi-analytic model (SAM) and the TNG hydrodynamic simulation. We find that the shuffling only slightly reduces the colour prediction accuracy in SAM and TNG, which is still considerably higher than that of the SDSS. This suggests that the true connection between SDSS colour and subhalo properties could be weaker than that in the SAM and TNG without the mismatch effect. We also measure the Pearson correlation coefficient between the galaxy and subhalo properties in SDSS, SAM, and TNG. Similar to the RF results, we find that the colour–subhalo correlation in SDSS is lower than in both SAM and TNG. We also show that the galaxy–subhalo correlations depend on subhalo mass in the galaxy formation models. Advanced surveys with fainter galaxies will provide new insights into the galaxy–subhalo relation in the real Universe.

We have applied computer analysis to classify the broad morphological types of ∼3 · 106 Sloan Digital Sky Survey (SDSS) galaxies. For each galaxy, the catalog provides the DR8 object ID, the R.A., the decl., and the certainty for the automatic classification as either spiral or elliptical. The certainty of the classification allows us to control the accuracy of a subset of galaxies by sacrificing some of the least certain classifications. The accuracy of the catalog was tested using galaxies that were classified by the manually annotated Galaxy Zoo catalog. The results show that the catalog contains ∼900,000 spiral galaxies and ∼600,000 elliptical galaxies with classification certainty that has a statistical agreement rate of ∼98% with the Galaxy Zoo debiased “superclean” data set. The catalog also shows that objects assigned by the SDSS pipeline with a relatively high redshift (z > 0.4) can have clear visual spiral morphology. The catalog can be downloaded at http://vfacstaff.ltu.edu/lshamir/data/morph_catalog. The image analysis software that was used to create the catalog is also publicly available.

We present a new suite of mock galaxy catalogues mimicking the low-redshift Universe, based on an updated halo occupation distribution (HOD) model and a scaling relation between optical properties and the neutral hydrogen (H i) content of galaxies. Our algorithm is constrained by observations of the luminosity function and luminosity- and colour-dependent clustering of Sloan Digital Sky Survey (SDSS) galaxies, as well as the H i mass function and H i-dependent clustering of massive H i-selected galaxies in the Arecibo Legacy Fast ALFA (ALFALFA) survey. Mock central and satellite galaxies with realistic values of r-band luminosity, g − r and u − r colour, stellar mass and H i mass are populated in an N-body simulation, inheriting a number of properties of the density and tidal environment of their host haloes. The host halo of each central galaxy is also ‘baryonified’ with realistic spatial distributions of stars as well as hot and cold gas, along with the corresponding rotation curve. Our default HOD assumes that galaxy properties are a function of group halo mass alone, and can optionally include effects such as galactic conformity and colour-dependent galaxy assembly bias. The mocks predict the relation between the stellar mass and H i mass of massive H i galaxies, as well as the 2-point cross-correlation function of spatially co-located optical and H i-selected samples. They enable novel null tests for galaxy assembly bias, provide predictions for the H i velocity width function, and clarify the origin and universality of the radial acceleration relation in the Lambda cold dark matter framework.

We consider nonlinear biasing models of galaxies with particular attention to a correlation between the linear and quadratic biasing coefficients, $b_1$ and $b_2$. We first derive perturbative expressions for $b_1$ and $b_2$ in halo and peak biasing models. We then discuss our computations of the power spectra and bispectra of dark matter particles and halos using $N$-body simulation data and of volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their $b_1$ and $b_2$. We find that the values of those coefficients at linear regimes ($k<0.2$$h$ Mpc$^{-1}$) are fairly insensitive to the redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of the bispectra, $Q$, for equilateral triangles are insensitive to the values of $b_1$, implying that $b_2$ indeed correlates with $b_1$. The present results explain the previous finding of Kayo et al. (2004, PASJ, 56, 415) for the hierarchical relation of three-point correlation functions of SDSS galaxies. While the relations between $b_1$ and $b_2$ are quantitatively different for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the biasing due to the gravity in primordial Gaussian density fields.

We present distributions of orbital parameters of infalling satellites of $\Lambda$CDM haloes in the mass range $10^{12}-10^{14}$M$_\odot$, which represent the initial conditions for the subsequent evolution of substructures within the host halo. We use merger trees constructed in a high resolution cosmological N-body simulation to trace satellite haloes, and identify the time of infall. We find signficant trends in the distribution of orbital parameters with both the host halo mass and the ratio of satellite-to-host halo masses. For all host halo masses, satellites whose infall mass is a larger fraction of the host halo mass have more eccentric, radially biased orbits. At fixed satellite-to-host halo mass ratio, high mass haloes are biased towards accreting satellites on slightly more radial orbits. To charactise the orbital distributions fully requires fitting the correlated bivariate distribution of two chosen orbital parameters (e.g. radial and tangential velocity or energy and angular momentum). We provide simple fits to one choice of the bivariate distributions, which when transformed faithfully, captures the behaviour of any of the projected one-dimensional distributions.

We use a sample of ~6000 galaxies detected by the Arecibo Legacy Fast ALFA (ALFALFA) 21cm survey, to measure the clustering properties of HI-selected galaxies. We find no convincing evidence for a dependence of clustering on the galactic atomic hydrogen (HI) mass, over the range M_HI ~ 10^{8.5} - 10^{10.5} M_sun. We show that previously reported results of weaker clustering for low-HI mass galaxies are probably due to finite-volume effects. In addition, we compare the clustering of ALFALFA galaxies with optically selected samples drawn from the Sloan Digital Sky Survey (SDSS). We find that HI-selected galaxies cluster more weakly than even relatively optically faint galaxies, when no color selection is applied. Conversely, when SDSS galaxies are split based on their color, we find that the correlation function of blue optical galaxies is practically indistinguishable from that of HI-selected galaxies. At the same time, SDSS galaxies with red colors are found to cluster significantly more than HI-selected galaxies, a fact that is evident in both the projected as well as the full two-dimensional correlation function. A cross-correlation analysis further reveals that gas-rich galaxies "avoid" being located within ~3 Mpc of optical galaxies with red colors. Next, we consider the clustering properties of halo samples selected from the Bolshoi LambdaCDM simulation. A comparison with the clustering of ALFALFA galaxies suggests that galactic HI mass is not tightly related to host halo mass, and that a sizable fraction of subhalos do not host HI galaxies. Lastly, we find that we can recover fairly well the correlation function of HI galaxies by just excluding halos with low spin parameter. This finding lends support to the hypothesis that halo spin plays a key role in determining the gas content of galaxies.

We explore observational constraints on possible deviations from Newtonian gravity by means of large-scale clustering of galaxies. We measure the power spectrum and the bispectrum of Sloan Digital Sky Survey galaxies and compare the result with predictions in an empirical model of modified gravity. Our model assumes an additional Yukawa-like term with two parameters that characterize the amplitude and the length scale of the modified gravity. The model predictions are calculated using two methods; second-order perturbation theory and direct $N$-body simulations. These methods allow us to study nonlinear evolution of large-scale structure. Using the simulation results, we find that perturbation theory provides reliable estimates for the power spectrum and the bispectrum in the modified Newtonian model. We also construct mock galaxy catalogs from the simulations, and derive constraints on the amplitude and the length scale of deviations from Newtonian gravity. The resulting constraints from the power spectrum are consistent with those obtained in our earlier, indicating the validity of the previous empirical modeling of gravitational nonlinearity in the modified Newtonian model. If linear biasing is adopted, the bispectrum of the Sloan Digital Sky Survey galaxies yields constraints very similar to those from the power spectrum. If we allow for the nonlinear biasing instead, we find that the ratio of the quadratic to linear biasing coefficients, ${b}_{2}/{b}_{1}$, should satisfy $\ensuremath{-}0.25<{b}_{2}/{b}_{1}<0.19$ for $\ensuremath{\lambda}=5{h}^{\ensuremath{-}1}\text{ }\text{ }\mathrm{Mpc}$ and $\ensuremath{-}0.19<{b}_{2}/{b}_{1}<0.13$ for $\ensuremath{\lambda}=10{h}^{\ensuremath{-}1}\text{ }\text{ }\mathrm{Mpc}$ in the modified Newtonian model.

We develop empirical methods for modeling the galaxy population and populating cosmological N-body simulations with mock galaxies according to the observed properties of galaxies in survey data. We use these techniques to produce a new set of mock catalogs for the DEEP2 Galaxy Redshift Survey based on the output of the high-resolution Bolshoi simulation, as well as two other simulations with different cosmological parameters, all of which we release for public use. The mock-catalog creation technique uses subhalo abundance matching to assign galaxy luminosities to simulated dark-matter halos. It then adds color information to the resulting mock galaxies in a manner that depends on the local galaxy density, in order to reproduce the measured color–environment relation in the data. In the course of constructing the catalogs, we test various models for including scatter in the relation between halo mass and galaxy luminosity, within the abundance-matching framework. We find that there is no constant-scatter model that can simultaneously reproduce both the luminosity function and the autocorrelation function of DEEP2. This result has implications for galaxy-formation theory, and it restricts the range of contexts in which the mock catalogs can be usefully applied. Nevertheless, careful comparisons show that our new mock catalogs accurately reproduce a wide range of the other properties of the DEEP2 catalog, suggesting that they can be used to gain a detailed understanding of various selection effects in DEEP2.

In a wide-area spectroscopic survey of galaxies, it is nearly impossible to obtain a homogeneous sample of galaxies with respect to galaxy properties such as stellar mass and host halo mass across a range of redshifts. Despite the selection effect, theoretical templates in most analyses assume single tracers when compared with the measured clustering quantities. We demonstrate analytically that the selection effect inevitably introduces a bias in the redshift-space power spectrum on scales from linear to nonlinear scales. To quantitatively assess the impact of the selection effect, we construct mock galaxy catalogs from halos in N-body simulations by selecting halos above redshift-dependent mass thresholds such that the resulting redshift distribution of the halos, n(z), matches that of SDSS-like galaxies. We find that the selection effect causes fractional changes of up to only 1% and 2% in the monopole and quadrupole moments of the redshift-space power spectrum at k ≲ 0.3 hMpc-1, respectively, compared to the moments for the single mass-threshold (therefore single tracer) sample, for n g(z) of the SDSS-like galaxy samples. We also argue that the selection effect is unlikely to cause a significant bias in the estimation of cosmological parameters using the Fisher matrix method, provided that the redshift-dependent selection effect is modest.

We explore the dependence of the radial alignment of subhaloes on the mass of the host halo they orbit in. As the effect is seen on a broad range of scales including massive clusters as well as galactic systems it only appears natural to explore this phenomenon by means of cosmological simulations covering the same range in masses. We have 25 well resolved host dark matter haloes at our disposal ranging from 1015h−1 M⊙ down to 1012h−1 M⊙ each consisting of order of a couple of million particles within the virial radius. We observe that subhaloes tend to be more spherical than isolated objects. Both the distributions of sphericity and triaxiality of subhaloes are Gaussian-distributed with peak values of 〈s〉≈ 0.80 and 〈T〉≈ 0.56, irrespective of host mass. Interestingly, we note that the radial alignment is independent of host halo mass and the distribution of cos θ (i.e. the angle between the major-axis Ea of each subhalo and the radius vector of the subhalo in the reference frame of the host) is well fitted by a simple power law P(cos θ) ∝ cos4θ with the same fitting parameters for all host haloes.

This is the third paper in a series that reports on our investigation of the clustering properties of AGNs identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this paper, we extend the redshift range to 0.07<z<0.50 and measure the clustering amplitudes of both X-ray and optically-selected SDSS broad-line AGNs with and without radio detections as well as for X-ray selected narrow-line RASS/SDSS AGNs. We measure the clustering amplitude through cross-correlation functions (CCFs) with SDSS galaxies and derive the bias by applying a halo occupation distribution (HOD) model directly to the CCFs. We find no statistically convincing difference in the clustering of X-ray and optically-selected broad-line AGNs, as well as with samples in which radio-detected AGNs are excluded. This is in contrast to low redshift optically-selected narrow-line AGNs, where radio-loud AGNs are found in more massive halos than optical AGNs without a radio-detection. The typical dark matter halo masses of our broad-line AGNs are log M_DMH/[h^(-1) M_SUN] ~ 12.4-13.4, consistent with the halo mass range of typical non-AGN galaxies at low redshifts. We find no significant difference between the clustering of X-ray selected narrow-line AGNs and broad-line AGNs. We confirm the weak dependence of the clustering strength on AGN X-ray luminosity at a ~2 sigma level. Finally, we summarize the current picture of AGN clustering to z~1.5 based on three dimensional clustering measurements.