Showing papers by "Ravi K. Sheth published in 2017"

The Halo Boundary of Galaxy Clusters in the SDSS

Abstract: Kavli Institute for Cosmological Physics at the University of Chicago [NSF PHY-1125897]; Kavli Foundation; US Department of Energy [DE-SC0007901]; DOE [DE-SC0015975]; Sloan Foundation [FG-2016-6443]; Alfred P. Sloan Foundation; U.S. Department of Energy Office of Science; Center for High-Performance Computing at the University of Utah

The high mass end of the stellar mass function: dependence on stellar population models and agreement between fits to the light profile

Abstract: We quantify the systematic effects on the stellar mass function that arise from assumptions about the stellar population, as well as how one fits the light profiles of the most luminous galaxies at z∼0.1. When comparing results from the literature, we are careful to separate out these effects. Our analysis shows that while systematics in the estimated comoving number density that arise from different treatments of the stellar population remain of the order of ≤0.5 dex, systematics in photometry are now about 0.1 dex, in contrast to some recent claims in the literature. Compared to these more recent analyses, previous work based on Sloan Digital Sky Survey pipeline photometry leads to underestimates of ρ∗(≥M∗)by factors of 3–10 in the mass range 1011–1011.6 M⊙, but up to a factor of 100 at higher stellar masses. This impacts studies that match massive galaxies to dark matter haloes. Although systematics that arise from different treatments of the stellar population remain of the order of ≤0.5 dex, our finding that systematics in photometry now amount to only about 0.1 dex in the stellar mass density is a significant improvement with respect to a decade ago. Our results highlight the importance of using the same stellar population and photometric models whenever low- and high-redshift samples are compared.

The Halo Boundary of Galaxy Clusters in the SDSS

Abstract: Mass around dark matter halos can be divided into "infalling" material and "collapsed" material that has passed through at least one pericenter. Analytical models and simulations predict a rapid drop in the halo density profile associated with the transition between these two regimes. Using data from SDSS, we explore the evidence for such a feature in the density profiles of galaxy clusters and investigate the connection between this feature and a possible phase space boundary. We first estimate the steepening of the outer galaxy density profile around clusters: the profiles show an abrupt steepening, providing evidence for truncation of the halo profile. Next, we measure the galaxy density profile around clusters using two sets of galaxies selected based on color. We find evidence of an abrupt change in the galaxy colors that coincides with the location of the steepening of the density profile. Since galaxies are likely to be quenched of star formation and turn red inside of clusters, this change in the galaxy color distribution can be interpreted as the transition from an infalling regime to a collapsed regime. We also measure this transition using a model comparison approach which has been used recently in studies of the "splashback" phenomenon, but find that this approach is not a robust way to quantify the significance of detecting a splashback-like feature. Finally, we perform measurements using an independent cluster catalog to test for potential systematic errors associated with cluster selection. We identify several avenues for future work: improved understanding of the small-scale galaxy profile, lensing measurements, identification of proxies for the halo accretion rate, and other tests. With upcoming data from the DES, KiDS and HSC surveys, we can expect significant improvements in the study of halo boundaries.

Comparing pymorph and SDSS photometry – II. The differences are more than semantics and are not dominated by intracluster light

Abstract: The Sloan Digital Sky Survey (SDSS) pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background, and (ii) single-component or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than the de Vaucouleurs model used by the SDSS pipeline. We use the pymorph photometric reductions to isolate effect (ii) and show that it is the same in the full sample as in small group environments, and for satellites in the most massive clusters as well. None of these are expected to be significantly affected by intracluster light (ICL). We only see an additional effect for centrals in the most massive haloes, but we argue that even this is not dominated by ICL. Hence, for the vast majority of galaxies, the differences between pymorph and SDSS pipeline photometry cannot be ascribed to the semantics of whether or not one includes the ICL when describing the stellar mass of massive galaxies. Rather, they likely reflect differences in star formation or assembly histories. Failure to account for the SDSS underestimate has significantly biased most previous estimates of the SDSS luminosity and stellar mass functions, and therefore halo model estimates of the z ∼ 0.1 relation between the mass of a halo and that of the galaxy at its centre. We also show that when one studies correlations, at fixed group mass, with a quantity that was not used to define the groups, then selection effects appear. We show why such effects arise and should not be mistaken for physical effects.

Selection bias in dynamically measured supermassive black hole samples: scaling relations and correlations between residuals in semi-analytic galaxy formation models

Abstract: Recent work has confirmed that the masses of supermassive black holes, estimated from scaling relations with global properties such as the stellar masses of their host galaxies, may be biased high. Much of this may be caused by the requirement that the gravitational sphere of influence of the black hole must be resolved for the black-hole mass to be reliably estimated. We revisit this issue by using a comprehensive galaxy evolution semi-analytic model, which self-consistently evolves supermassive black holes from high-redshift seeds via gas accretion and mergers, and also includes AGN feedback. Once tuned to reproduce the (mean) correlation of black-hole mass with velocity dispersion, the model is unable to also account for the correlation with stellar mass. This behaviour is independent of the model's parameters, thus suggesting an internal inconsistency in the data. The predicted distributions, especially at the low-mass end, are also much broader than observed. However, if selection effects are included, the model's predictions tend to align with the observations. We also demonstrate that the correlations between the residuals of the local scaling relations are more effective than the scaling relations themselves at constraining AGN feedback models. In fact, we find that our semi-analytic model, while in apparent broad agreement with the scaling relations when accounting for selection biases, yields very weak correlations between their residuals at fixed stellar mass, in stark contrast with observations. This problem persists when changing the AGN feedback strength, and is also present in the $z\sim 0$ outputs of the hydrodynamic cosmological simulation Horizon-AGN, which includes state-of-the-art treatments of AGN feedback. This suggests that current AGN feedback models may be too weak or are simply not capturing the effect of the black hole on the stellar velocity dispersion.

Selection bias in dynamically-measured super-massive black hole samples: dynamical masses and dependence on Sérsic index

Abstract: We extend the comparison between the set of local galaxies having dynamically measured black holes with galaxies in the Sloan Digital Sky Survey (SDSS). We first show that the most up-to-date local black hole samples of early-type galaxies with measurements of effective radii, luminosities, and Sersic indices of the bulges of their host galaxies, have dynamical mass and Sersic index distributions consistent with those of SDSS early-type galaxies of similar bulge stellar mass. The host galaxies of local black hole samples thus do not appear structurally different from SDSS galaxies, sharing similar dynamical masses, light profiles and light distributions. Analysis of the residuals reveals that velocity dispersion is more fundamental than Sersic index n in the scaling relations between black holes and galaxies. Indeed, residuals with Sersic index could be ascribed to the (weak) correlation with bulge mass or even velocity dispersion. Finally, targetted Monte Carlo simulations that include the effects of the sphere of influence of the black hole, and tuned to reproduce the observed residuals and scaling relations in terms of velocity dispersion and stellar mass, show that, at least for galaxies with Mbulge > 1e10 Msun and n>5, the observed mean black hole mass at fixed Sersic index is biased significantly higher than the intrinsic value.

Revisiting the bulge-halo conspiracy I: dependence on galaxy properties and halo mass

Abstract: We carry out a systematic investigation of the total mass density profile of massive ($\mathrm{log}\,{M}_{\mathrm{star}}/{M}_{\odot }\gtrsim 11.3$) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies with measured stellar masses, effective radii, and Sersic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical ΛCDM profile. Our main results are as follows. (1) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius γ′ is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (2) Models characterized by a Salpeter or variable initial mass function (IMF) and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier IMF and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (3) Currently available data on the mass density profiles of very massive galaxies ($\mathrm{log}\,{M}_{\mathrm{star}}/{M}_{\odot }\gtrsim 12$), with ${M}_{\mathrm{halo}}\gtrsim 3\times {10}^{14}\,{M}_{\odot }$, favor instead models with a stellar profile flatter than a Sersic one in the very inner regions (r ≲ 3–5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ∼2 or ≲1.3, respectively, higher than those typically predicted by N-body numerical simulations.

Effective window function for Lagrangian halos

Abstract: The window function for protohalos in Lagrangian space is often assumed to be a top hat in real space. We measure this profile directly and find that it is more extended than a top hat but less extended than a Gaussian; its shape is well described by rounding the edges of the top hat by convolution with a Gaussian that has a scale length about 5 times smaller. This effective window ${W}_{\mathrm{eff}}$ is particularly simple in Fourier space, and has an analytic form in real space. Together with the excursion set bias parameters, ${W}_{\mathrm{eff}}$ describes the scale-dependence of the Lagrangian halo-matter cross correlation up to $k{R}_{\text{Lag}}\ensuremath{\sim}10$, where ${R}_{\text{Lag}}$ is the Lagrangian size of the protohalo. Moreover, with this ${W}_{\mathrm{eff}}$, all the spectral moments of the power spectrum are finite, allowing a straightforward estimate of the excursion set peak mass function. This estimate requires a prescription of the critical overdensity enclosed within a protohalo if it is to collapse, which we calibrate from simulations. We find that the resulting estimate of halo abundances is only accurate to about 20%, and we discuss why: a top hat in ``infall time'' towards the protohalo center need not correspond to a top hat in the initial spatial distribution, so models in which infall rather than smoothed overdensity is the relevant variable may be more accurate.

Radial Acceleration Relation between Baryons and Dark or Phantom Matter in the Super-critical Acceleration Regime of Nearly Spherical Galaxies

Abstract: The central regions of nearby elliptical galaxies are dominated by baryons (stars) and provide interesting laboratories for studying the radial acceleration relation (RAR). We carry out exploratory analyses and discuss the possibility of constraining the RAR in the super-critical acceleration range $(10^{-9.5},\hspace{1ex}10^{-8})$~${\rm m}~{\rm s}^{-2}$ using a sample of nearly round pure-bulge (spheroidal, dispersion-dominated) galaxies including 24 ATLAS$^{\rm 3D}$ galaxies and 4201 SDSS galaxies covering a wide range of masses, sizes and luminosity density profiles. We consider a range of current possibilities for the stellar mass-to-light ratio ($M_\star/L$), its gradient and dark or phantom matter (DM/PM) halo profiles. We obtain the probability density functions (PDFs) of the parameters of the considered models via Bayesian inference based on spherical Jeans Monte Carlo modeling of the observed velocity dispersions. We then constrain the DM/PM-to-baryon acceleration ratio $a_{\rm X}/a_{\rm B}$ from the PDFs. Unless we ignore observed radial gradients in $M_\star/L$, or assume unreasonably strong gradients, marginalization over nuisance factors suggests $a_{\rm X}/a_{\rm B} = 10^{p} (a_{\rm B}/a_{+1})^q$ with $p = -1.00 \pm 0.03$ (stat) $^{+0.11}_{-0.06}$ (sys) and $q=-1.02 \pm 0.09$ (stat) $^{+0.16}_{-0.00}$ (sys) around a super-critical acceleration $a_{+1}\equiv 1.2\times 10^{-9}~{\rm m}~{\rm s}^{-2}$. In the context of the $\Lambda$CDM paradigm, this RAR suggests that the NFW DM halo profile is a reasonable description of galactic halos even after the processes of galaxy formation and evolution. In the context of the MOND paradigm, this RAR favors the Simple interpolating function but is inconsistent with the vast majority of other theoretical proposals and fitting functions motivated mainly from sub-critical acceleration data.

Constraints on halo formation from cross-correlations with correlated variables

Abstract: Cross-correlations between biased tracers and the dark matter field encode information about the physical variables which characterize these tracers. However, if the physical variables of interest are correlated with one another, then extracting this information is not as straightforward as one might naively have thought. We show how to exploit these correlations so as to estimate scale-independent bias factors of all orders in a model-independent way. We also show that failure to account for this will lead to incorrect conclusions about which variables matter and which do not. Morever, accounting for this allows one to use the scale dependence of bias to constrain the physics of halo formation; to date the argument has been phrased the other way around. We illustrate by showing that the scale dependence of linear and nonlinear bias, measured on nonlinear scales, can be used to provide consistent estimates of how the critical density for halo formation depends on halo mass. Our methods work even when the bias is nonlocal and stochastic, such as when, in addition to the spherically averaged density field and its derivatives, the quadrupolar shear field also matters for halo formation. In such models, the nonlocal bias factors are closely related to the more familiar local nonlinear bias factors, which are much easier to measure. Our analysis emphasizes the fact that biased tracers are biased because they do not sample fields (density, velocity, shear, etc.) at all positions in space in the same way that the dark matter does.

Consistency relations for the Lagrangian halo bias and their implications

Abstract: The protohalo patches from which halos form are defined by a number of constraints imposed on the Lagrangian dark matter density field. Each of these constraints contributes to biasing the spatial distribution of the protohalos relative to the matter. We show how measurements of this spatial distribution -- linear combinations of protohalo bias factors -- can be used to make inferences about the physics of halo formation. Our analysis exploits the fact that halo bias factors satisfy consistency relations which encode this physics, and that these relations are the same even for sub-populations in which assembly bias has played a role. We illustrate our methods using a model in which three parameters matter: a density threshold, the local slope and the curvature of the smoothed density field. The latter two are nearly degenerate; our approach naturally allows one to build an accurate effective two-parameter model for which the consistency relations still apply. This, with an accurate description of the smoothing window, allows one to describe the protohalo-matter cross-correlation very well, both in Fourier and configuration space. We then use our determination of the large scale bias parameters together with the consistency relations, to estimate the enclosed density and mean slope on the Lagrangian radius scale of the protohalos. Direct measurements of these quantities, made on smaller scales than those on which the bias parameters are typically measured, are in good agreement.

Bimodal Formation Time Distribution for Infall Dark Matter Halos

Abstract: We use a 200 $h^{-1}Mpc$ a side N-body simulation to study the mass accretion history (MAH) of dark matter halos to be accreted by larger halos, which we call infall halos. We define a quantity $a_{\rm nf}\equiv (1+z_{\rm f})/(1+z_{\rm peak})$ to characterize the MAH of infall halos, where $z_{\rm peak}$ and $z_{\rm f}$ are the accretion and formation redshifts, respectively. We find that, at given $z_{\rm peak}$, their MAH is bimodal. Infall halos are dominated by a young population at high redshift and by an old population at low redshift. For the young population, the $a_{\rm nf}$ distribution is narrow and peaks at about $1.2$, independent of $z_{\rm peak}$, while for the old population, the peak position and width of the $a_{\rm nf}$ distribution both increases with decreasing $z_{\rm peak}$ and are both larger than those of the young population. This bimodal distribution is found to be closely connected to the two phases in the MAHs of halos. While members of the young population are still in the fast accretion phase at $z_{\rm peak}$, those of the old population have already entered the slow accretion phase at $z_{\rm peak}$. This bimodal distribution is not found for the whole halo population, nor is it seen in halo merger trees generated with the extended Press-Schechter formalism. The infall halo population at $z_{\rm peak}$ are, on average, younger than the whole halo population of similar masses identified at the same redshift. We discuss the implications of our findings in connection to the bimodal color distribution of observed galaxies and to the link between central and satellite galaxies.

Radial Acceleration Relation by Dark Matter and Baryons in the Intermediate Acceleration Regime

Abstract: We investigate the relation between radial accelerations due to the baryonic mass distribution ($a_{\rm B}$) and the dark (or phantom) matter distribution ($a_{\rm X}$) in the intermediate acceleration range $10^{-10} 10^{-9.5}~{\rm m}~{\rm s}^{-2}$. Using a generalized MOND model of the form $a_{\rm B}/a= (a/a_0)/\left[1+(a/a_0)^ u\right]^{1/ u}$ between the empirical acceleration $a$ and the Newtonian prediction $a_{\rm B}$ by the baryons, we find $a_0=(1.22\pm 0.08)\times 10^{-10}~{\rm m}~{\rm s}^{-2}$ and $ u=0.88 \pm 0.15$ for $10^{-10} < a_{\rm B} /[{\rm m}~{\rm s}^{-2}] < 10^{-8.5}$. Our results constrain theories of dark matter or modified gravity.

Revisiting the bulge-halo conspiracy I: Dependence on galaxy properties and halo mass

Abstract: We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.