Showing 1–50 of 204 results for author: Sheth, R K

Eigen-decomposition of Covariance matrices: An application to the BAO Linear Point

Authors: Jaemyoung Jason Lee, Farnik Nikakhtar, Aseem Paranjape, Ravi K. Sheth

Abstract: The Baryon Acoustic Oscillation (BAO) feature in the two-point correlation function (TPCF) of discrete tracers such as galaxies is an accurate standard ruler. The covariance matrix of the TPCF plays an important role in determining how the precision of this ruler depends on the number density and clustering strength of the tracers, as well as the survey volume. An eigen-decomposition of this matri… ▽ More The Baryon Acoustic Oscillation (BAO) feature in the two-point correlation function (TPCF) of discrete tracers such as galaxies is an accurate standard ruler. The covariance matrix of the TPCF plays an important role in determining how the precision of this ruler depends on the number density and clustering strength of the tracers, as well as the survey volume. An eigen-decomposition of this matrix provides an objective way to separate the contributions of cosmic variance from those of shot-noise to the statistical uncertainties. For the signal-to-noise levels that are expected in ongoing and next-generation surveys, the cosmic variance eigen-modes dominate. These modes are smooth functions of scale, meaning that: they are insensitive to the modest changes in binning that are allowed if one wishes to resolve the BAO feature in the TPCF; they provide a good description of the correlated residuals which result from fitting smooth functional forms to the measured TPCF; they motivate a simple but accurate approximation for the uncertainty on the Linear Point (LP) estimate of the BAO distance scale. This approximation allows one to quantify the precision of the BAO distance scale estimate without having to generate a large ensemble of mock catalogs and explains why: the uncertainty on the LP does not depend on the functional form fitted to the TPCF or the binning used; the LP is more constraining than the peak or dip scales in the TPCF; the evolved TPCF is less constraining than the initial one, so that reconstruction schemes can yield significant gains in precision. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 10 pages, 8 figures, submitted to Physical Review D

The Stability of the BAO Linear Point under Modified Gravity

Authors: Jaemyoung Jason Lee, Bartolomeo Fiorini, Farnik Nikakhtar, Ravi K. Sheth

Abstract: Baryon Acoustic Oscillations (BAOs) are crucial in cosmological analysis, providing a standard ruler, as well as constraints on dark energy. In General Relativity models, the BAO Linear Point - the midpoint between the dip and the peak in the correlation function - has been shown to be rather robust to evolution and redshift space distortions. We show that this remains true even when the gravity m… ▽ More Baryon Acoustic Oscillations (BAOs) are crucial in cosmological analysis, providing a standard ruler, as well as constraints on dark energy. In General Relativity models, the BAO Linear Point - the midpoint between the dip and the peak in the correlation function - has been shown to be rather robust to evolution and redshift space distortions. We show that this remains true even when the gravity model is not General Relativity, at least for $f(R)$ and DGP gravity models which have the same expansion history as the standard $Λ$CDM. For the Linear Point to be able to distinguish between modified gravity (MG) and $Λ$CDM, survey volumes of order tens of cubic Gpc are required. △ Less

Submitted 19 June, 2024; v1 submitted 13 June, 2024; originally announced June 2024.

Comments: 9 pages, 5 figures, submitted to Physical Review D, v2

The energy shear of protohaloes

Authors: Marcello Musso, Giulia Despali, Ravi K. Sheth

Abstract: As it collapses to form a halo, the shape of a protohalo patch is deformed by the initial shear field. This deformation is often modeled using the "deformation" tensor, constructed from second derivatives of the gravitational potential, whose trace gives the initial overdensity. However, especially for lower mass protohalos, this matrix is not always positive definite: one of its eigenvalues has a… ▽ More As it collapses to form a halo, the shape of a protohalo patch is deformed by the initial shear field. This deformation is often modeled using the "deformation" tensor, constructed from second derivatives of the gravitational potential, whose trace gives the initial overdensity. However, especially for lower mass protohalos, this matrix is not always positive definite: one of its eigenvalues has a different sign from the others. We show that the evolution of a patch is better described by the "energy shear" tensor, which is positive definite and plays a direct role in the evolution. This positive-definiteness simplifies models of halo abundances, assembly and of the cosmic web. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 9 pages, 10 figures

Displacement Field Analysis via Optimal Transport: Multi-Tracer Approach to Cosmological Reconstruction

Authors: Farnik Nikakhtar, Ravi K. Sheth, Nikhil Padmanabhan, Bruno Lévy, Roya Mohayaee

Abstract: We demonstrate the effectiveness of one of the many multi-tracer analyses enabled by Optimal Transport (OT) reconstruction. Leveraging a semi-discrete OT algorithm, we determine the displacements between initial and observed positions of biased tracers and the remaining matter field. With only redshift-space distorted final positions of biased tracers and a simple premise for the remaining mass di… ▽ More We demonstrate the effectiveness of one of the many multi-tracer analyses enabled by Optimal Transport (OT) reconstruction. Leveraging a semi-discrete OT algorithm, we determine the displacements between initial and observed positions of biased tracers and the remaining matter field. With only redshift-space distorted final positions of biased tracers and a simple premise for the remaining mass distribution as input, OT solves the displacement field. This extracted field, assuming asymptotically uniform density and a gradient flow displacement, enables reconstruction of the initial overdensity fluctuation field. We show that the divergence of the OT displacement field is a good proxy of the linear density field, even though the method never assumes the linear theory growth. Additionally, this divergence field can be combined with the reconstructed protohalos to provide a higher signal-to-noise measurement of the BAO standard ruler than was possible with either measurement individually. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 12 pages, 10 figures

The statistics of primordial black holes in a radiation dominated Universe -- recent and new results

Authors: Cristiano Germani, Ravi K. Sheth

Abstract: We review the non-linear statistics of Primordial Black Holes that form from the collapse of over-densities in a radiation dominated Universe. We focus on the scenario in which large over-densities are generated by rare and Gaussian curvature perturbations during inflation. As new results, we show that the mass spectrum follows a power law determined by the critical exponent of the self-similar co… ▽ More We review the non-linear statistics of Primordial Black Holes that form from the collapse of over-densities in a radiation dominated Universe. We focus on the scenario in which large over-densities are generated by rare and Gaussian curvature perturbations during inflation. As new results, we show that the mass spectrum follows a power law determined by the critical exponent of the self-similar collapse up to a power spectrum dependent cut-off, and that the abundance related to very narrow power spectra is exponentially suppressed. Related to this, we discuss and explicitly show that the Press-Schechter approximation, as well as the statistics of mean profiles, lead to wrong conclusions for the abundance and mass spectrum. Finally, we clarify that the transfer function in the statistics of initial conditions for primordial black holes formation (the abundance) does not play a significant role. △ Less

Submitted 14 September, 2023; v1 submitted 5 August, 2023; originally announced August 2023.

Comments: v2 typos corrected and clarifications added. Results unchanged. Version accepted for publication as a review in the journal Universe, special issue on PBHs. 18 pages; 7 figures

Optimal Transport Reconstruction of Biased Tracers in Redshift Space

Authors: Farnik Nikakhtar, Nikhil Padmanabhan, Bruno Lévy, Ravi K. Sheth, Roya Mohayaee

Abstract: Recent research has emphasized the benefits of accurately reconstructing the initial Lagrangian positions of biased tracers from their positions at a later time, to gain cosmological information. A weighted semi-discrete optimal transport algorithm can achieve the required accuracy, provided the late-time positions are known, with minimal information about the background cosmology. The algorithm's… ▽ More Recent research has emphasized the benefits of accurately reconstructing the initial Lagrangian positions of biased tracers from their positions at a later time, to gain cosmological information. A weighted semi-discrete optimal transport algorithm can achieve the required accuracy, provided the late-time positions are known, with minimal information about the background cosmology. The algorithm's performance relies on knowing the masses of the biased tracers, and depends on how one models the distribution of the remaining mass that is not associated with these tracers. We demonstrate that simple models of the remaining mass result in accurate retrieval of the initial Lagrangian positions, which we quantify using pair statistics and the void probability function. This is true even if the input positions are affected by redshift-space distortions. The most sophisticated models assume that the masses of the tracers, and the amount and clustering of the missing mass are known; we show that the method is robust to realistic errors in the masses of the tracers and remains so as the model for the missing mass becomes increasingly crude. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Comments: 15 pages, 18 figures

Model-agnostic cosmological constraints from the baryon acoustic oscillation feature in redshift space

Authors: Aseem Paranjape, Ravi K. Sheth

Abstract: We develop a framework for self-consistently extracting cosmological information from the clustering of tracers in redshift space, $\textit{without}$ relying on model-dependent templates to describe the baryon acoustic oscillation (BAO) feature. Our approach uses the recently proposed Laguerre reconstruction technique for the BAO feature and its linear point $r_{\rm LP}$, and substantially extends… ▽ More We develop a framework for self-consistently extracting cosmological information from the clustering of tracers in redshift space, $\textit{without}$ relying on model-dependent templates to describe the baryon acoustic oscillation (BAO) feature. Our approach uses the recently proposed Laguerre reconstruction technique for the BAO feature and its linear point $r_{\rm LP}$, and substantially extends it to simultaneously model the multipoles $\ell=0,2,4$ of the anisotropic galaxy 2-point correlation function (2pcf). The approach is `model-agnostic': it assumes that the non-linear growth of structure smears the BAO feature by an approximately Gaussian kernel with a smearing scale $σ_{\rm v}$, but does not assume any fiducial cosmology for describing the shape of the feature itself. Using mock observations for two realistic survey configurations assuming $Λ$ cold dark matter ($Λ$CDM), combined with Bayesian parameter inference, we show that the linear point $r_{\rm LP}$ and smearing scale $σ_{\rm v}$ can be accurately recovered by our method in both existing and upcoming surveys. The precision of the recovery of $r_{\rm LP}$ is always better than $1\%$, while $σ_{\rm v}$ can be recovered with $\lesssim10\%$ uncertainty provided the linear galaxy bias $b$ is separately constrained, e.g., using weak lensing observations. Our method is also sensitive to the linear growth rate $f$, albeit with larger uncertainties and systematic errors, especially for upcoming surveys such as DESI. We discuss how our model can be modified to improve the recovery of $f$, such that the resulting constraints on $\{f,σ_{\rm v},r_{\rm LP}\}$ can potentially be used as a test of cosmological models including and beyond $Λ$CDM. △ Less

Submitted 8 September, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 18 pages, 6 figures; v2: minor changes, accepted in MNRAS

Getting in shape with minimal energy. A variational principle for protohaloes

Authors: Marcello Musso, Ravi K. Sheth

Abstract: In analytical models of structure formation, protohalos are routinely assumed to be peaks of the smoothed initial density field, with the smoothing filter being spherically symmetric. This works reasonably well for identifying a protohalo's center of mass, but not its shape. To provide a more realistic description of protohalo boundaries, one must go beyond the spherical picture. We suggest that t… ▽ More In analytical models of structure formation, protohalos are routinely assumed to be peaks of the smoothed initial density field, with the smoothing filter being spherically symmetric. This works reasonably well for identifying a protohalo's center of mass, but not its shape. To provide a more realistic description of protohalo boundaries, one must go beyond the spherical picture. We suggest that this can be done by looking for regions of fixed volume, but arbitrary shape, that minimize the enclosed energy. Such regions are surrounded by surfaces over which (a slightly modified version of) the gravitational potential is constant. We show that these equipotential surfaces provide an excellent description of protohalo shapes, orientations and associated torques. △ Less

Submitted 22 February, 2023; originally announced March 2023.

Comments: 5 pages, 6 figures

Stellar population analysis of MaNGA early-type galaxies: IMF dependence and systematic effects

Authors: M. Bernardi, H. Dominguez Sanchez, R. K. Sheth, J. R. Brownstein, R. R. Lane

Abstract: We study systematics associated with estimating simple stellar population (SSP) parameters -- age, metallicity [M/H], $α$-enhancement [$α$/Fe] and IMF shape -- and associated $M_*/L$ gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs) and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) h… ▽ More We study systematics associated with estimating simple stellar population (SSP) parameters -- age, metallicity [M/H], $α$-enhancement [$α$/Fe] and IMF shape -- and associated $M_*/L$ gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs) and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H$_β$, $\langle$Fe$\rangle$, TiO$_{\rm 2SDSS}$ and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius $R_e$: (a) ages are younger and [$α$/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about $R_e/2$; (c) this makes $M_*/L$ about $2\times$ larger in the central regions than beyond $R_e/2$. While the models of Conroy et al. (2018) return similar [M/H] and [$α$/Fe] profiles, the age and (hence) $M_*/L$ profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs and S0s: younger ages and less enhanced [$α$/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred $M_*/L$ gradients significantly. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 22 pages, 23 figures, accepted for publication in MNRAS

Bayesian evidence comparison for distance scale estimates

Authors: Aseem Paranjape, Ravi K. Sheth

Abstract: Constraints on cosmological parameters are often distilled from sky surveys by fitting templates to summary statistics of the data that are motivated by a fiducial cosmological model. However, recent work has shown how to estimate the distance scale using templates that are more generic: the basis functions used are not explicitly tied to any one cosmological model. We describe a Bayesian framewor… ▽ More Constraints on cosmological parameters are often distilled from sky surveys by fitting templates to summary statistics of the data that are motivated by a fiducial cosmological model. However, recent work has shown how to estimate the distance scale using templates that are more generic: the basis functions used are not explicitly tied to any one cosmological model. We describe a Bayesian framework for (i) determining how many basis functions to use and (ii) comparing one basis set with another. Our formulation provides intuition into how (a) one's degree of belief in different basis sets, (b) the fact that the choice of priors depends on basis set, and (c) the data set itself, together determine the derived constraints. We illustrate our framework using measurements in simulated datasets before applying it to real data. △ Less

Submitted 14 October, 2022; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: 10 pages, 7 figures; v2: minor clarifications, fixed some typos, accepted in MNRAS

Optimal Transport Reconstruction of Baryon Acoustic Oscillations

Authors: Farnik Nikakhtar, Ravi K. Sheth, Bruno Lévy, Roya Mohayaee

Abstract: A weighted, semi-discrete, fast optimal transport (OT) algorithm for reconstructing the Lagrangian positions of proto-halos from their evolved Eulerian positions is presented. The algorithm makes use of a mass estimate of the biased tracers and of the distribution of the remaining mass (the `dust'), but is robust to errors in the mass estimates. Tests with state-of-art cosmological simulations sho… ▽ More A weighted, semi-discrete, fast optimal transport (OT) algorithm for reconstructing the Lagrangian positions of proto-halos from their evolved Eulerian positions is presented. The algorithm makes use of a mass estimate of the biased tracers and of the distribution of the remaining mass (the `dust'), but is robust to errors in the mass estimates. Tests with state-of-art cosmological simulations show that if the dust is assumed to have a uniform spatial distribution, then the shape of the OT-reconstructed pair correlation function of the tracers is very close to linear theory, enabling sub-percent precision in the BAO distance scale that depends weakly, if at all, on a cosmological model. With a more sophisticated model for the dust, OT returns an estimate of the displacement field which yields superb reconstruction of the proto-halo positions, and hence of the shape and amplitude of the initial pair correlation function of the tracers. This enables direct and independent determinations of the bias factor $b$ and the smearing scale $Σ$, potentially providing new methods for breaking the degeneracy between $b$ and $σ_8$. △ Less

Submitted 9 January, 2023; v1 submitted 3 March, 2022; originally announced March 2022.

Comments: matches published version (Physical Review Letters)

The half mass radius of MaNGA galaxies: Effect of IMF gradients

Authors: M. Bernardi, R. K. Sheth, H. Dominguez Sanchez, B. Margalef-Bentabol, D. Bizyaev, R. R. Lane

Abstract: Gradients in the stellar populations (SP) of galaxies -- e.g., in age, metallicity, stellar Initial Mass Function (IMF) -- can result in gradients in the stellar mass to light ratio, $M_*/L$. Such gradients imply that the distribution of the stellar mass and light are different. For old SPs, e.g., in early-type galaxies at $z\sim 0$, the $M_*/L$ gradients are weak if driven by variations in age an… ▽ More Gradients in the stellar populations (SP) of galaxies -- e.g., in age, metallicity, stellar Initial Mass Function (IMF) -- can result in gradients in the stellar mass to light ratio, $M_*/L$. Such gradients imply that the distribution of the stellar mass and light are different. For old SPs, e.g., in early-type galaxies at $z\sim 0$, the $M_*/L$ gradients are weak if driven by variations in age and metallicity, but significantly larger if driven by the IMF. A gradient which has larger $M_*/L$ in the center increases the estimated total stellar mass ($M_*$) and reduces the scale which contains half this mass ($R_{e,*}$), compared to when the gradient is ignored. For the IMF gradients inferred from fitting MILES simple SP models to the H$_β$, $\langle$Fe$\rangle$, [MgFe] and TiO$_{\rm 2SDSS}$ absorption lines measured in spatially resolved spectra of early-type galaxies in the MaNGA survey, the fractional change in $R_{e,*}$ can be significantly larger than that in $M_*$, especially when the light is more centrally concentrated. The $R_{e,*}-M_*$ correlation which results is offset by 0.3 dex to smaller sizes compared to when these gradients are ignored. Comparisons with `quiescent' galaxies at higher-$z$ must account for evolution in SP gradients (especially age and IMF) and the light profile before drawing conclusions about how $R_{e,*}$ and $M_*$ evolve. The implied merging between higher-$z$ and the present is less contrived if $R_{e,*}/R_e$ at $z\sim 0$ is closer to our IMF-driven gradient calibration than to unity. △ Less

Submitted 15 November, 2022; v1 submitted 19 January, 2022; originally announced January 2022.

Comments: 16 pages, 15 figures, accepted for publication in MNRAS

The phenomenology of the external field effect in cold dark matter models

Authors: Aseem Paranjape, Ravi K. Sheth

Abstract: In general relativity (GR), the internal dynamics of a self-gravitating system under free-fall in an external gravitational field should not depend on the external field strength. Recent work has claimed a statistical detection of an `external field effect' (EFE) using galaxy rotation curve data. We show that large uncertainties in rotation curve analyses and inaccuracies in published simulation-b… ▽ More In general relativity (GR), the internal dynamics of a self-gravitating system under free-fall in an external gravitational field should not depend on the external field strength. Recent work has claimed a statistical detection of an `external field effect' (EFE) using galaxy rotation curve data. We show that large uncertainties in rotation curve analyses and inaccuracies in published simulation-based external field estimates compromise the significance of the claimed EFE detection. We further show analytically that a qualitatively similar statistical signal is, in fact, expected in a $Λ$-cold dark matter ($Λ$CDM) universe without any violation of the strong equivalence principle. Rather, such a signal arises simply because of the inherent correlations between galaxy clustering strength and intrinsic galaxy properties. We explicitly demonstrate the effect in a baryonified mock catalog of a $Λ$CDM universe. Although the detection of an EFE-like signal is not, by itself, evidence for physics beyond GR, our work shows that the $\textit{sign}$ of the EFE-like correlation between the external field strength and the shape of the radial acceleration relation can be used to probe new physics: e.g., in MOND, the predicted sign is opposite to that in our $Λ$CDM mocks. △ Less

Submitted 20 September, 2022; v1 submitted 30 November, 2021; originally announced December 2021.

Comments: 11 pages, 6 figures; v2: minor changes and added discussion, accepted in MNRAS

The smearing scale in Laguerre reconstructions of the correlation function

Authors: Farnik Nikakhtar, Ravi K. Sheth, Idit Zehavi

Abstract: To a good approximation, on large cosmological scales the evolved two-point correlation function of biased tracers is related to the initial one by a convolution. For Gaussian initial conditions, the smearing kernel is Gaussian, so if the initial correlation function is parametrized using simple polynomials then the evolved correlation function is a sum of generalized Laguerre functions of half-in… ▽ More To a good approximation, on large cosmological scales the evolved two-point correlation function of biased tracers is related to the initial one by a convolution. For Gaussian initial conditions, the smearing kernel is Gaussian, so if the initial correlation function is parametrized using simple polynomials then the evolved correlation function is a sum of generalized Laguerre functions of half-integer order. This motivates an analytic Laguerre reconstruction algorithm which previous work has shown is fast and accurate. This reconstruction requires as input the width of the smearing kernel. We show that the method can be extended to estimate the width of the smearing kernel from the same dataset. This estimate, and associated uncertainties, can then be used to marginalize over the distribution of reconstructed shapes, and hence provide error estimates on the value of the distance scale which are not tied to a particular cosmological model. We also show that if, instead, we parametrize the evolved correlation function using simple polynomials, then the initial one is a sum of Hermite polynomials, again enabling fast and accurate deconvolution. If one is willing to use constraints on the smearing scale from other datasets, then marginalizing over its value is simpler for Hermite reconstruction, potentially providing further speed-up in cosmological analyses. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 11 pages, 6 figures

Laguerre reconstruction of the BAO feature in halo-based mock galaxy catalogues

Authors: Farnik Nikakhtar, Ravi K. Sheth, Idit Zehavi

Abstract: Fitting half-integer generalized Laguerre functions to the evolved, real-space dark matter and halo correlation functions provides a simple way to reconstruct their initial shapes. We show that this methodology also works well in a wide variety of realistic, assembly biased, velocity biased and redshift-space distorted mock galaxy catalogs. We use the linear point feature in the monopole of the re… ▽ More Fitting half-integer generalized Laguerre functions to the evolved, real-space dark matter and halo correlation functions provides a simple way to reconstruct their initial shapes. We show that this methodology also works well in a wide variety of realistic, assembly biased, velocity biased and redshift-space distorted mock galaxy catalogs. We use the linear point feature in the monopole of the redshift-space distorted correlation function to quantify the accuracy of our approach. We find that the linear point estimated from the mock galaxy catalogs is insensitive to the details of the biasing scheme at the sub-percent level. However, the linear point scale in the nonlinear, biased, and redshift-space distorted field is systematically offset from its scale in the unbiased linear density fluctuation field by more than 1%. In the Laguerre reconstructed correlation function, this is reduced to sub-percent values, so it provides comparable accuracy and precision to methods that reconstruct the full density field before estimating the distance scale. The linear point in the reconstructed density fields provided by these other methods is likewise precise, accurate, and insensitive to galaxy bias. All reconstructions depend on some input parameters, and marginalizing over uncertainties in the input parameters required for reconstruction can degrade both accuracy and precision. The linear point simplifies the marginalization process, enabling more realistic estimates of the precision of the distance scale estimate for negligible additional computational cost. We show this explicitly for Laguerre reconstruction. △ Less

Submitted 26 July, 2021; originally announced July 2021.

Comments: Accepted for publication in PRD, 11 pages, 11 figures

Journal ref: Phys. Rev. D 104, 063504 (2021)

The distribution of HI velocity profiles in a $Λ$CDM universe

Authors: Aseem Paranjape, R. Srianand, Tirthankar Roy Choudhury, Ravi K. Sheth

Abstract: We model the distribution of the observed profiles of 21 cm line emission from neutral hydrogen (HI) in central galaxies selected from a statistically representative mock catalog of the local Universe in the Lambda-cold dark matter framework. The distribution of these HI velocity profiles (specifically, their widths $W_{50}$) has been observationally constrained, but has not been systematically st… ▽ More We model the distribution of the observed profiles of 21 cm line emission from neutral hydrogen (HI) in central galaxies selected from a statistically representative mock catalog of the local Universe in the Lambda-cold dark matter framework. The distribution of these HI velocity profiles (specifically, their widths $W_{50}$) has been observationally constrained, but has not been systematically studied theoretically. Our model profiles derive from rotation curves of realistically baryonified haloes in an N-body simulation, including the quasi-adiabatic relaxation of the dark matter profile of each halo in response to its baryons. We study the predicted $W_{50}$ distribution using a realistic pipeline applied to noisy profiles extracted from our luminosity-complete mock catalog with an ALFALFA-like survey geometry and redshift selection. Our default mock is in good agreement with observed ALFALFA results for $W_{50}\gtrsim700$ km/s, being incomplete at lower widths due to the intrinsic threshold of $M_r\leq-19$. Variations around the default model show that the velocity width function at $W_{50}\gtrsim300$ km/s is most sensitive to a possible correlation between galaxy inclination and host concentration, followed by the physics of quasi-adiabatic relaxation. We also study the excess kurtosis of noiseless velocity profiles, obtaining a distribution which tightly correlates with $W_{50}$, with a shape and scatter that depend on the properties of the turbulent HI disk. Our results open the door towards using the shapes of HI velocity profiles as a novel statistical probe of the baryon-dark matter connection. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 22 pages, 14 figures, submitted to MNRAS

The radial acceleration relation in a $Λ$CDM universe

Authors: Aseem Paranjape, Ravi K. Sheth

Abstract: We study the radial acceleration relation (RAR) between the total ($a_{\rm tot}$) and baryonic ($a_{\rm bary}$) centripetal acceleration profiles of central galaxies in the cold dark matter (CDM) paradigm. We analytically show that the RAR is intimately connected with the physics of the quasi-adiabatic relaxation of dark matter in the presence of baryons in deep potential wells. This cleanly demon… ▽ More We study the radial acceleration relation (RAR) between the total ($a_{\rm tot}$) and baryonic ($a_{\rm bary}$) centripetal acceleration profiles of central galaxies in the cold dark matter (CDM) paradigm. We analytically show that the RAR is intimately connected with the physics of the quasi-adiabatic relaxation of dark matter in the presence of baryons in deep potential wells. This cleanly demonstrates how the mean RAR and its scatter emerge in the low-acceleration regime ($10^{-12}\,{\rm m\,s}^{-2}\lesssim a_{\rm bary}\lesssim10^{-10}\,{\rm m\,s}^{-2}$) from an interplay between baryonic feedback processes and the distribution of CDM in dark halos. Our framework allows us to go further and study both higher and lower accelerations in detail, using analytical approximations and a realistic mock catalog of $\sim342,000$ low-redshift central galaxies with $M_r\leq-19$. We show that, while the RAR in the baryon-dominated, high-acceleration regime ($a_{\rm bary}\gtrsim10^{-10}\,{\rm m\,s}^{-2}$) is very sensitive to details of the relaxation physics, a simple `baryonification' prescription matching the relaxation results of hydrodynamical CDM simulations is remarkably successful in reproducing the observed RAR without any tuning. And in the (currently unobserved) ultra-low-acceleration regime ($a_{\rm bary}\lesssim 10^{-12}\,{\rm m\,s}^{-2}$), the RAR is sensitive to the abundance of diffuse gas in the halo outskirts, with our default model predicting a distinctive break from a simple power-law-like relation for HI-deficient, diffuse gas-rich centrals. Our mocks also show that the RAR provides more robust, testable predictions of the $Λ$CDM paradigm at galactic scales, with implications for alternative gravity theories, than the baryonic Tully-Fisher relation. △ Less

Submitted 23 July, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: 20 pages, 16 figures; v2 - added discussion, clarifications and references, conclusions unchanged, accepted in MNRAS

Laguerre reconstruction of the correlation function on Baryon Acoustic Oscillation scales

Authors: Farnik Nikakhtar, Ravi K. Sheth, Idit Zehavi

Abstract: The baryon acoustic oscillation feature can be used as a standard cosmological ruler. In practice, for sub-percent level accuracy on the distance scale, it must be standardized. The physical reason why is understood, so we use this to develop an algorithm which improves the estimated scale. The algorithm exploits the fact that, over the range of scales where the initial correlation function is wel… ▽ More The baryon acoustic oscillation feature can be used as a standard cosmological ruler. In practice, for sub-percent level accuracy on the distance scale, it must be standardized. The physical reason why is understood, so we use this to develop an algorithm which improves the estimated scale. The algorithm exploits the fact that, over the range of scales where the initial correlation function is well-fit by a polynomial, the leading order effects which distort the length of the ruler can be accounted for analytically. Tests of the method in numerical simulations show that it provides simple and fast reconstruction of the full shape of the BAO feature, as well as sub-percent determination of the linear point in the correlation function of biased tracers with minimal assumptions about the underlying cosmological model or the nature of the observed tracers. Our results also suggest that, for least squares estimators of the correlation function, half-integer generalized Laguerre functions are a particularly useful choice. △ Less

Submitted 26 July, 2021; v1 submitted 20 January, 2021; originally announced January 2021.

Comments: Accepted for publication in PRD, 16 pages, 17 figures

Journal ref: Phys. Rev. D 104, 043530 (2021)

Multi-wavelength mock galaxy catalogs of the low-redshift Universe

Authors: Aseem Paranjape, Tirthankar Roy Choudhury, Ravi K. Sheth

Abstract: We present a new suite of mock galaxy catalogs 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 (HI) content of galaxies. Our algorithm is constrained by observations of the luminosity function and luminosity- and colour-dependent clustering of SDSS galaxies, as well as the H… ▽ More We present a new suite of mock galaxy catalogs 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 (HI) content of galaxies. Our algorithm is constrained by observations of the luminosity function and luminosity- and colour-dependent clustering of SDSS galaxies, as well as the HI mass function and HI-dependent clustering of massive HI-selected galaxies in the ALFALFA survey. Mock central and satellite galaxies with realistic values of $r$-band luminosity, $g-r$ and $u-r$ colour, stellar mass and HI mass are populated in an $N$-body simulation, inheriting a number of properties of the density and tidal environment of their host halos. 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 HI mass of massive HI galaxies, as well as the 2-point cross-correlation function of spatially co-located optical and HI-selected samples. They enable novel null tests for galaxy assembly bias, provide predictions for the HI velocity width function, and clarify the origin and universality of the radial acceleration relation in the $Λ$CDM framework. △ Less

Submitted 8 March, 2021; v1 submitted 7 January, 2021; originally announced January 2021.

Comments: 17 pages, 13 figures; v2: minor changes and clarifications, accepted in MNRAS

Mock halo catalogs: assigning unresolved halo properties using correlations with local halo environment

Authors: Sujatha Ramakrishnan, Aseem Paranjape, Ravi K. Sheth

Abstract: Large-scale sky surveys require companion large volume simulated mock catalogs. To ensure precision cosmology studies are unbiased, the correlations in these mocks between galaxy properties and their large-scale environments must be realistic. Since galaxies are embedded in dark matter halos, an important first step is to include such correlations -- sometimes called assembly bias -- for dark matt… ▽ More Large-scale sky surveys require companion large volume simulated mock catalogs. To ensure precision cosmology studies are unbiased, the correlations in these mocks between galaxy properties and their large-scale environments must be realistic. Since galaxies are embedded in dark matter halos, an important first step is to include such correlations -- sometimes called assembly bias -- for dark matter halos. However, galaxy properties correlate with smaller scale physics in halos which large simulations struggle to resolve. We describe an algorithm which addresses and largely mitigates this problem. Our algorithm exploits the fact that halo assembly bias is unchanged as long as correlations between halo property $c$ and the intermediate-scale tidal environment $α$ are preserved. Therefore, knowledge of $α$ is sufficient to assign small-scale, otherwise unresolved properties to a halo in a way which preserves its large-scale assembly bias accurately. We demonstrate this explicitly for halo internal properties like formation history (concentration $c_{\rm 200b}$), shape $c/a$, dynamics $c_{v}/a_{v}$, velocity anisotropy $β$ and angular momentum (spin $λ$). Our algorithm increases a simulation's reach in halo mass and number density by an order of magnitude, with improvements in the bias signal as large as 45% for 30-particle halos, thus significantly reducing the cost of mocks for future weak lensing and redshift space distortion studies. △ Less

Submitted 22 February, 2021; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: 13 pages, 10 figures, minor revisions, accepted in MNRAS

On the Presence of a Universal Acceleration Scale in Elliptical Galaxies

Authors: Kyu-Hyun Chae, Mariangela Bernardi, Helena Dominguez Sanchez, Ravi K. Sheth

Abstract: Dark matter phenomena in rotationally supported galaxies exhibit a characteristic acceleration scale of $g_\dagger \approx 1.2\times 10^{-10}$ m s$^{-2}$. Whether this acceleration is a manifestation of a universal scale, or merely an emergent property with an intrinsic scatter, has been debated in the literature. Here we investigate whether a universal acceleration scale exists in dispersion-supp… ▽ More Dark matter phenomena in rotationally supported galaxies exhibit a characteristic acceleration scale of $g_\dagger \approx 1.2\times 10^{-10}$ m s$^{-2}$. Whether this acceleration is a manifestation of a universal scale, or merely an emergent property with an intrinsic scatter, has been debated in the literature. Here we investigate whether a universal acceleration scale exists in dispersion-supported galaxies using two uniform sets of integral field spectroscopy (IFS) data from SDSS-IV MaNGA and ATLAS$^{\rm 3D}$. We apply the spherical Jeans equation to 15 MaNGA and 4 ATLAS$^{\rm 3D}$ slow-rotator E0 (i.e., nearly spherical) galaxies. Velocity dispersion profiles for these galaxies are well determined with observational errors under control. Bayesian inference indicates that all 19 galaxies are consistent with a universal acceleration of $g_\dagger=1.5_{-0.6}^{+0.9}\times 10^{-10}$ m s$^{-2}$. Moreover, all 387 data points from the radial bins of the velocity dispersion profiles are consistent with a universal relation between the radial acceleration traced by dynamics and that predicted by the observed distribution of baryons. This universality remains if we include 12 additional non-E0 slow-rotator elliptical galaxies from ATLAS$^{\rm 3D}$. Finally, the universal acceleration from MaNGA and ATLAS$^{\rm 3D}$ is consistent with that for rotationally supported galaxies, so our results support the view that dark matter phenomenology in galaxies involves a universal acceleration scale. △ Less

Submitted 17 November, 2020; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: Published in ApJ Letters

Journal ref: 2020 November 10, 903, L31

Analytical thresholds for black hole formation in general cosmological backgrounds

Authors: Albert Escrivà, Cristiano Germani, Ravi K. Sheth

Abstract: We consider black holes which form from an initially spherically symmetric super-Hubble perturbation of a cosmological background filled by a perfect fluid $p=w ρ$ with $w\in (0,1]$. Previous work has shown that when $w = 1/3$ (radiation), there is a critical threshold for black hole formation ($δ_c$), which, to a very good approximation, only depends upon the curvature of the compaction function… ▽ More We consider black holes which form from an initially spherically symmetric super-Hubble perturbation of a cosmological background filled by a perfect fluid $p=w ρ$ with $w\in (0,1]$. Previous work has shown that when $w = 1/3$ (radiation), there is a critical threshold for black hole formation ($δ_c$), which, to a very good approximation, only depends upon the curvature of the compaction function around its peak value. We find that this generalizes to all $w\gtrsim 1/3$; for smaller $w$s the knowledge of the full shape of the compaction function is necessary. We provide analytic approximations for $δ_c$ which are accurate for $w\in [1/3,1]$. △ Less

Submitted 16 November, 2020; v1 submitted 10 July, 2020; originally announced July 2020.

Comments: 19 pages and 15 figures. v2: small changes. Accepted for publication in JCAP

Report number: ICCUB-20-016

Galaxy properties as revealed by MaNGA. III. Kinematic profiles and stellar population gradients in S0s

Authors: H. Domínguez Sánchez, M. Bernardi, F. Nikakhtar, B. Margalef-Bentabol, R. K. Sheth

Abstract: This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, $α$-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at $z\le 0.08$ from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by m… ▽ More This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, $α$-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at $z\le 0.08$ from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by measuring Lick indices and comparing them to stellar population synthesis models. This requires spectra with high signal-to-noise which we achieve by stacking in bins of luminosity (L$_r$) and central velocity dispersion ($σ_0$). We find that: 1) There is a bimodality in the S0 population: S0s more massive than $3\times 10^{10}M_\odot$ show stronger velocity dispersion and age gradients (age and $σ_r$ decrease outwards) but little or no metallicity gradient, while the less massive ones present relatively flat age and velocity dispersion profiles, but a significant metallicity gradient (i.e. [M/H] decreases outwards). Above $2\times10^{11}M_\odot$ the number of S0s drops sharply. These two mass scales are also where global scaling relations of ETGs change slope. 2) S0s have steeper velocity dispersion profiles than fast rotating elliptical galaxies (E-FRs) of the same luminosity and velocity dispersion. The kinematic profiles and stellar population gradients of E-FRs are both more similar to those of slow rotating ellipticals (E-SRs) than to S0s, suggesting that E-FRs are not simply S0s viewed face-on. 3) At fixed $σ_0$, more luminous S0s and E-FRs are younger, more metal rich and less $α$-enhanced. Evidently for these galaxies, the usual statement that 'massive galaxies are older' is not true if $σ_0$ is held fixed. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS. 15 pages, 20 figures

The Stellar Mass Fundamental Plane: The virial relation and a very thin plane for slow-rotators

Authors: M. Bernardi, H. Domínguez Sanchez, B. Margalef-Bentabol, F. Nikakhtar, R. K. Sheth

Abstract: Early-type galaxies -- slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars -- define a Fundamental Plane (FP) in the space of half-light radius $R_e$, enclosed surface brightness $I_e$ and velocity dispersion $σ_e$. Since $I_e$ and $σ_e$ are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which $I_e$ and $σ_e$ can be… ▽ More Early-type galaxies -- slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars -- define a Fundamental Plane (FP) in the space of half-light radius $R_e$, enclosed surface brightness $I_e$ and velocity dispersion $σ_e$. Since $I_e$ and $σ_e$ are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which $I_e$ and $σ_e$ can be used to estimate sizes $R_e$. We show that: 1) The thickness of the FP depends strongly on morphology. If the sample only includes E-SRs, then the observed scatter in $R_e$ is $\sim 16\%$, of which only $\sim 9\%$ is intrinsic. Removing galaxies with $M_*<10^{11}M_\odot$ further reduces the observed scatter to $\sim 13\%$ ($\sim 4\%$ intrinsic). The observed scatter increases to the $\sim 25\%$ usually quoted in the literature if E-FRs and S0s are added. If the FP is defined using the eigenvectors of the covariance matrix of the observables, then the E-SRs again define an exceptionally thin FP, with intrinsic scatter of only $5\%$ orthogonal to the plane. 2) The structure within the FP is most easily understood as arising from the fact that $I_e$ and $σ_e$ are nearly independent, whereas the $R_e-I_e$ and $R_e-σ_e$ correlations are nearly equal and opposite. 3) If the coefficients of the FP differ from those associated with the virial theorem the plane is said to be `tilted'. If we multiply $I_e$ by the global stellar mass-to-light ratio $M_*/L$ and we account for non-homology across the population by using Sérsic photometry, then the resulting stellar mass FP is less tilted. Accounting self-consistently for $M_*/L$ gradients will change the tilt. The tilt we currently see suggests that the efficiency of turning baryons into stars increases and/or the dark matter fraction decreases as stellar surface brightness increases. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: 13 pages, 9 figures, 3 tables, accepted for publication in MNRAS

Non-linear statistics of primordial black holes from gaussian curvature perturbations

Authors: Cristiano Germani, Ravi K. Sheth

Abstract: We develop the non-linear statistics of primordial black holes generated by a gaussian spectrum of primordial curvature perturbations. This is done by employing the compaction function as the main statistical variable under the constraints that: a) the over-density has a high peak at a point $\vec{x}_0$, b) the compaction function has a maximum at a smoothing scale $R$, and finally, c) the compact… ▽ More We develop the non-linear statistics of primordial black holes generated by a gaussian spectrum of primordial curvature perturbations. This is done by employing the compaction function as the main statistical variable under the constraints that: a) the over-density has a high peak at a point $\vec{x}_0$, b) the compaction function has a maximum at a smoothing scale $R$, and finally, c) the compaction function amplitude at its maximum is higher than the threshold necessary to trigger a gravitational collapse into a black hole of the initial over-density. Our calculation allows for the fact that the patches which are destined to form PBHs may have a variety of profile shapes and sizes. The predicted PBH abundances depend on the power spectrum of primordial fluctuations. For a very peaked power spectrum, our non-linear statistics, the one based on the linear over-density and the one based on the use of curvature perturbations, all predict a narrow distribution of PBH masses and comparable abundance. For broader power spectra the linear over-density statistics over-estimate the abundance of primordial black holes while the curvature-based approach under-estimates it. Additionally, for very large smoothing scales, the abundance is no longer dominated by the contribution of a mean over-density but rather by the whole statistical realisations of it. △ Less

Submitted 10 February, 2020; v1 submitted 15 December, 2019; originally announced December 2019.

Comments: v2: improved discussions and citations. Constraint on the abundance related to minimums added. Discussion on the broad power spectrum and related appendix added. 19 pages, 5 figures

Report number: ICCUB-19-021

Journal ref: Phys. Rev. D 101, 063520 (2020)

Probing black hole accretion tracks, scaling relations and radiative efficiencies from stacked X-ray active galactic nuclei

Authors: Francesco Shankar, David H. Weinberg, Christopher Marsden, Philip J. Grylls, Mariangela Bernardi, Guang Yang, Benjamin Moster, Rosamaria Carraro, David M. Alexander, Viola Allevato, Tonima T. Ananna, Angela Bongiorno, Giorgio Calderone, Francesca Civano, Emanuele Daddi, Ivan Delvecchio, Federica Duras, Fabio La Franca, Andrea Lapi, Youjun Lu, Nicola Menci, Mar Mezcua, Federica Ricci, Giulia Rodighiero, Ravi K. Sheth , et al. (3 additional authors not shown)

Abstract: The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises f… ▽ More The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 12 pages, 8 Figures. MNRAS, accepted. A discussion around the points raised by arXiv:1909.10821 is included in the Appendix

Linear Point and Sound Horizon as Purely Geometric standard rulers

Authors: Márcio O'Dwyer, Stefano Anselmi, Glenn D. Starkman, Pier-Stefano Corasaniti, Ravi K. Sheth, Idit Zehavi

Abstract: The Baryon Acoustic Oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the Purely Geometric BAO methods. To date two different Purely Geometric BAO approache… ▽ More The Baryon Acoustic Oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the Purely Geometric BAO methods. To date two different Purely Geometric BAO approaches have been proposed. The first exploits the linear-point standard ruler. The second, called correlation-function model-fitting, exploits the sound-horizon standard ruler. A key difference between them is that, when estimated from clustering data, the linear point makes use of a cosmological-model-independent procedure to extract the ratio of the ruler to the cosmic distance, while the correlation-function model-fitting relies on a phenomenological cosmological model for the correlation function. Nevertheless the two rulers need to be precisely defined independently of any specific observable. We define the linear point and sound horizon and we characterize and compare the two rulers' cosmological-parameter dependence. We find that they are both geometrical within the required accuracy, and they have the same parameter dependence for a wide range of parameter values. We estimate the rulers' best-fit values and errors given the cosmological constraints obtained by the Planck Satellite team from the CMB measurements. We do this for three different cosmological models encompassed by the Purely Geometric BAO methods. In each case we find that the relative errors of the two rulers coincide and they are insensitive to the assumed cosmological model. Interestingly both the linear point and the sound horizon shift by $0.5σ$ when we do not fix the spatial geometry to be flat in LCDM. This points toward a sensitivity of the rulers to different cosmological models when they are estimated from the CMB. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. D 101, 083517 (2020)

Constraining black hole-galaxy scaling relations from the large-scale clustering of Active Galactic Nuclei and implied mean radiative efficiency

Authors: Francesco Shankar, Viola Allevato, Mariangela Bernardi, Christopher Marsden, Andrea Lapi, Nicola Menci, Philip J. Grylls, Mirko Krumpe, Lorenzo Zanisi, Federica Ricci, Fabio La Franca, Ranieri D. Baldi, Jorge Moreno, Ravi K. Sheth

Abstract: A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demo… ▽ More A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 31 pages, 7 Figures. Accepted

A universal threshold for primordial black hole formation

Authors: Albert Escrivà, Cristiano Germani, Ravi K. Sheth

Abstract: In this letter, we argue and show numerically that the threshold to form primordial black holes from an initial spherically symmetric perturbation is, to an excellent approximation, universal, whenever given in terms of the compaction function averaged over a sphere of radius $r_m$, where $r_m$ is the scale on which the compaction function is maximum. This can be understood as the requirement that… ▽ More In this letter, we argue and show numerically that the threshold to form primordial black holes from an initial spherically symmetric perturbation is, to an excellent approximation, universal, whenever given in terms of the compaction function averaged over a sphere of radius $r_m$, where $r_m$ is the scale on which the compaction function is maximum. This can be understood as the requirement that, for a black hole to form, each shell of the averaged compaction function should have an amplitude exceeding the so-called Harada-Yoo-Kohri limit. For a radiation dominated universe we argued, supported by the numerical simulations, that this limit is $δ_c = 0.40$, which is slightly below the one quoted in the literature. Additionally, we show that the profile dependence of the threshold for the compaction function is only sensitive to its curvature at the maximum. We use these results to provide an analytic formula for the threshold amplitude of the compaction function at its maximum in terms of the normalised compaction function curvature at $r_m$. △ Less

Submitted 27 January, 2020; v1 submitted 31 July, 2019; originally announced July 2019.

Comments: 5 pages and 5 figures, v3: clarifications and more numerical evidences for the lower bound added. Version accepted in PRD

Report number: ICC-19-013

Journal ref: Phys. Rev. D 101, 044022 (2020)

Excursion set peaks in energy as a model for haloes

Authors: Marcello Musso, Ravi K. Sheth

Abstract: The simplest models of dark matter halo formation are based on the heuristic assumption, motivated by spherical collapse, that virialized haloes originate from initial regions that are maxima of the smoothed density field. Here, we replace this notion with the dynamical requirement that protohalo patches be regions where the local gravitational flow converges to a point. For this purpose, we look… ▽ More The simplest models of dark matter halo formation are based on the heuristic assumption, motivated by spherical collapse, that virialized haloes originate from initial regions that are maxima of the smoothed density field. Here, we replace this notion with the dynamical requirement that protohalo patches be regions where the local gravitational flow converges to a point. For this purpose, we look for spheres whose gravitational acceleration at the boundary -- relative to their center of mass -- points towards their geometric center: that is, spheres with null dipole moment. We show that these configurations are minima of the total energy, i.e. the most energetically bound spheres. For this reason, we study peaks of the energy overdensity field, and argue that the approach shows considerable promise. This change simply requires that one modify the standard top-hat filter, with the added important benefit that, for power spectra of cosmological interest, the resulting model is no longer plagued by divergences. Although the formalism is no more complicated than the overdensity based approach, the model is richer in the sense that it naturally predicts scatter in the overdensities of protohalo patches that are destined to form haloes of the same mass, in qualitative agreement with simulations of halo formation. △ Less

Submitted 27 October, 2021; v1 submitted 22 July, 2019; originally announced July 2019.

Comments: 9 pages, 5 figures

Journal ref: MNRAS (2021) 508-3, 3634-3648

Exact enumeration approach to first-passage time distribution of non-Markov random walks

Authors: Shant Baghram, Farnik Nikakhtar, M. Reza Rahimi Tabar, Sohrab Rahvar, Ravi K. Sheth, Klaus Lehnertz, Muhammad Sahimi

Abstract: We propose an analytical approach to study non-Markov random walks by employing an exact enumeration method. Using the method, we derive an exact expansion for the first-passage time (FPT) distribution for any continuous, differentiable non-Markov random walk with Gaussian or non-Gaussian multivariate distribution. As an example, we study the FPT distribution of a fractional Brownian motion with a… ▽ More We propose an analytical approach to study non-Markov random walks by employing an exact enumeration method. Using the method, we derive an exact expansion for the first-passage time (FPT) distribution for any continuous, differentiable non-Markov random walk with Gaussian or non-Gaussian multivariate distribution. As an example, we study the FPT distribution of a fractional Brownian motion with a Hurst exponent $H\in(1/2,1)$ that describes numerous non-Markov stochastic phenomena in physics, biology and geology, and for which the limit $H=1/2$ represents a Markov process. △ Less

Submitted 5 June, 2019; originally announced June 2019.

Comments: 23 pages, 4 figures, 1 table and 5 appendices. Published version

Journal ref: Physical Review E 99, 062101 (2019)

Constraints on the redshift evolution of astrophysical feedback with Sunyaev-Zeldovich effect cross-correlations

Authors: S. Pandey, E. J. Baxter, Z. Xu, J. Orlowski-Scherer, N. Zhu, A. Lidz, J. Aguirre, J. DeRose, M. Devlin, J. C. Hill, B. Jain, R. K. Sheth, S. Avila, E. Bertin, D. Brooks, E. Buckley-Geer, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, F. J. Castander, R. Cawthon, L. N. da Costa, J. De Vicente, S. Desai, H. T. Diehl , et al. (34 additional authors not shown)

Abstract: An understanding of astrophysical feedback is important for constraining models of galaxy formation and for extracting cosmological information from current and future weak lensing surveys. The thermal Sunyaev-Zel'dovich effect, quantified via the Compton-$y$ parameter, is a powerful tool for studying feedback, because it directly probes the pressure of the hot, ionized gas residing in dark matter… ▽ More An understanding of astrophysical feedback is important for constraining models of galaxy formation and for extracting cosmological information from current and future weak lensing surveys. The thermal Sunyaev-Zel'dovich effect, quantified via the Compton-$y$ parameter, is a powerful tool for studying feedback, because it directly probes the pressure of the hot, ionized gas residing in dark matter halos. Cross-correlations between galaxies and maps of Compton-$y$ obtained from cosmic microwave background surveys are sensitive to the redshift evolution of the gas pressure, and its dependence on halo mass. In this work, we use galaxies identified in year one data from the Dark Energy Survey and Compton-$y$ maps constructed from Planck observations. We find highly significant (roughly $12σ$) detections of galaxy-$y$ cross-correlation in multiple redshift bins. By jointly fitting these measurements as well as measurements of galaxy clustering, we constrain the halo bias-weighted, gas pressure of the Universe as a function of redshift between $0.15 \lesssim z \lesssim 0.75$. We compare these measurements to predictions from hydrodynamical simulations, allowing us to constrain the amount of thermal energy in the halo gas relative to that resulting from gravitational collapse. △ Less

Submitted 30 April, 2019; originally announced April 2019.

Comments: 21 pages, 12 figures, comments welcome

Journal ref: Phys. Rev. D 100, 063519 (2019)

Galaxy properties as revealed by MaNGA II. Differences in stellar populations of slow and fast rotator ellipticals and dependence on environment

Authors: M. Bernardi, H. Domínguez Sánchez, J. R. Brownstein, N. Drory, R. K. Sheth

Abstract: We present estimates of stellar population (SP) gradients from stacked spectra of slow (SR) and fast (FR) rotator elliptical galaxies from the MaNGA-DR15 survey. We find that: 1) FRs are $\sim 5$ Gyrs younger, more metal rich, less $α$-enhanced and smaller than SRs of the same luminosity $L_r$ and central velocity dispersion $σ_0$. This explains why when one combines SRs and FRs, objects which are… ▽ More We present estimates of stellar population (SP) gradients from stacked spectra of slow (SR) and fast (FR) rotator elliptical galaxies from the MaNGA-DR15 survey. We find that: 1) FRs are $\sim 5$ Gyrs younger, more metal rich, less $α$-enhanced and smaller than SRs of the same luminosity $L_r$ and central velocity dispersion $σ_0$. This explains why when one combines SRs and FRs, objects which are small for their $L_r$ and $σ_0$ tend to be younger. Their SP gradients are also different. 2) Ignoring the FR/SR dichotomy leads one to conclude that compact galaxies are older than their larger counterparts of the same mass, even though almost the opposite is true for FRs and SRs individually. 3) SRs with $σ_0\le 250$ km s$^{-1}$ are remarkably homogeneous within $\sim R_e$: they are old, $α$-enhanced and only slightly super-solar in metallicity. These SRs show no gradients in age and $M_*/L_r$, negative gradients in metallicity, and slightly positive gradients in [$α$/Fe] (the latter are model dependent). SRs with $σ_0\ge 250$ km $s^{-1}$ are slightly younger and more metal rich, contradicting previous work suggesting that age increases with $σ_0$. They also show larger $M_*/L_r$ gradients. 4) Self-consistently accounting for $M_*/L$ gradients yields $M_{\rm dyn}\approx M_*$ because gradients reduce $M_{\rm dyn}$ by $\sim 0.2$ dex while only slightly increasing the $M_*$ inferred using a Kroupa (not Salpeter) IMF. 5) The FR population all but disappears above $M_*\ge 3\times 10^{11}M_\odot$; this is the same scale at which the size-mass correlation and other scaling relations change. Our results support the finding that this is an important mass scale which correlates with the environment and above which mergers matter. △ Less

Submitted 27 August, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

Comments: 22 pages, 28 figures. Accepted for publication in MNRAS

Galaxy properties as revealed by MaNGA. I. Constraints on Initial Mass Function and M$_{*}$/L gradients in ellipticals

Authors: H. Domínguez Sánchez, M. Bernardi, J. R. Brownstein, N. Drory, R. K. Sheth

Abstract: We estimate ages, metallicities, $α$-element abundance ratios and stellar initial mass functions of elliptical (E) and S0 galaxies from the MaNGA-DR15 survey. We stack spectra and use a variety of single stellar population synthesis models to interpret the absorption line strengths in these spectra. We quantify how these properties vary across the population, as well as with galactocentric distanc… ▽ More We estimate ages, metallicities, $α$-element abundance ratios and stellar initial mass functions of elliptical (E) and S0 galaxies from the MaNGA-DR15 survey. We stack spectra and use a variety of single stellar population synthesis models to interpret the absorption line strengths in these spectra. We quantify how these properties vary across the population, as well as with galactocentric distance. This paper is the first of a series and is based on a sample of pure elliptical galaxies at z $\le$ 0.08. We show that the properties of the inner regions of Es with the largest luminosity (L$_r$) and central velocity dispersion ($σ_0$) are consistent with those associated with the commonly used Salpeter IMF, whereas a Kroupa-like IMF is a better description at $\sim$ 0.8R/Re (assuming [Ti/Fe] variations are limited). For these galaxies the stellar mass-to-light ratio decreases at most by a factor of 2 from the central regions to Re. In contrast, for lower L$_r$ and $σ_0$ galaxies, the IMF is shallower and M$_{*}$/L$_r$ in the central regions is similar to the outskirts. Although a factor of 2 is smaller than previous reports based on a handful of galaxies, it is still large enough to matter for dynamical mass estimates. Accounting self-consistently for these gradients when estimating both M$_{*}$ and M$_{dyn}$ brings the two into good agreement: gradients reduce M$_{dyn}$ by $\sim$ 0.2 dex while only slightly increasing the M$_{*}$ inferred using a Kroupa IMF. This is a different resolution of the M$_{*}$-M$_{dyn}$ discrepancy than has been followed in the recent literature where M$_{*}$ of massive galaxies is increased by adopting a Salpeter IMF while leaving Mdyn unchanged. A companion paper discusses how stellar population differences are even more pronounced if one separates slow from fast rotators. △ Less

Submitted 27 August, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

Comments: 22 pages, 29 figures, accepted for publication in MNRAS

Cosmic voids: a novel probe to shed light on our Universe

Authors: Alice Pisani, Elena Massara, David N. Spergel, David Alonso, Tessa Baker, Yan-Chuan Cai, Marius Cautun, Christopher Davies, Vasiliy Demchenko, Olivier Doré, Andy Goulding, Mélanie Habouzit, Nico Hamaus, Adam Hawken, Christopher M. Hirata, Shirley Ho, Bhuvnesh Jain, Christina D. Kreisch, Federico Marulli, Nelson Padilla, Giorgia Pollina, Martin Sahlén, Ravi K. Sheth, Rachel Somerville, Istvan Szapudi , et al. (4 additional authors not shown)

Abstract: Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust… ▽ More Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White Paper for the Astro2020 Decadal Survey, 5 pages text, 14 pages with references and affiliations

Cosmic web anisotropy is the primary indicator of halo assembly bias

Authors: Sujatha Ramakrishnan, Aseem Paranjape, Oliver Hahn, Ravi K. Sheth

Abstract: The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $α$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represen… ▽ More The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $α$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represent two distinct statistical links: one between the internal property and $α$, and the other between $α$ and large-scale ( $>30h^{-1}$Mpc) halo bias $b_1$. We focus on scalar internal properties of haloes related to formation time (concentration $c_{\rm vir}$), shape (mass ellipsoid asphericity $c/a$), velocity dispersion structure (velocity ellipsoid asphericity $c_v/a_v$ and velocity anisotropy $β$) and angular momentum (dimensionless spin $λ$) in the mass range $8\times10^{11}< M_{\rm vir}/(h^{-1}M_\odot)<5\times10^{14}$. Using conditional correlation coefficients and other detailed tests, we show that the joint distribution of $α$, $b_1$ and any of the internal properties $c\in\{β,c_v/a_v,c/a,c_{\rm vir},λ\}$ is consistent with $p(α,b_1,c)\simeq p(α)p(b_1|α)p(c|α)$, at all but the largest masses. $\textit{Thus, the assembly bias trends $c-b_1$ reflect the two fundamental correlations $c-α$ and $b_1-α$.}$ Our results are unaffected by the exclusion of haloes with recent major merger events or splashback objects, although the latter are distinguished by the fact that $α$ does not explain their assembly bias trends. The overarching importance of $α$ provides a new perspective on the nature of assembly bias of distinct haloes, with potential ramifications for incorporating realistic assembly bias effects into mock catalogs of future large-scale structure surveys and for detecting galaxy assembly bias. △ Less

Submitted 20 August, 2019; v1 submitted 5 March, 2019; originally announced March 2019.

Comments: 21 pages, 14 figures; v2: Added discussion, results strengthened. Conclusions unchanged. Accepted by MNRAS

Modeling Nearly Spherical Pure-Bulge Galaxies with a Stellar Mass-to-Light Ratio Gradient under the $Λ$CDM and MOND Paradigms: II. The Orbital Anisotropy of Slow Rotators within the Effective Radius

Authors: Kyu-Hyun Chae, Mariangela Bernardi, Ravi K. Sheth

Abstract: We investigate the anisotropy of the stellar velocity dispersions within the effective radius, $R_{\rm e}$, in 24 ATLAS$^{\rm 3D}$ pure-bulge galaxies, 16 of which are kinematic slow rotators (SRs). We allow the spherical anisotropy parameter $β$ to be radially varying and allow a radial gradient in the stellar mass-to-light ratio ($M_\star/L$) through the parameter $K$ introduced earlier. The med… ▽ More We investigate the anisotropy of the stellar velocity dispersions within the effective radius, $R_{\rm e}$, in 24 ATLAS$^{\rm 3D}$ pure-bulge galaxies, 16 of which are kinematic slow rotators (SRs). We allow the spherical anisotropy parameter $β$ to be radially varying and allow a radial gradient in the stellar mass-to-light ratio ($M_\star/L$) through the parameter $K$ introduced earlier. The median anisotropy for SRs depends on $K$ as follows: $\langleβ_{\rm m}\rangle = a + b K$ with $a=0.19\pm 0.05$, $b=-0.13\pm 0.07$ ($Λ$CDM) or $a=0.21\pm 0.05$, $b=-0.26\pm 0.08$ (MOND), where $β_{\rm m}$ refers to the radially averaged quantity. Under the $Λ$CDM paradigm this scaling is tied to a scaling of $\langle f_{\rm DM}\rangle = (0.16\pm 0.03) +(0.31\pm 0.06) K$, where $f_{\rm DM}$ refers to the DM fraction within a sphere of $r=R_{\rm e}$. For $K=0$ (constant $M_\star/L$), we obtain radially biased results with $\langleβ_{\rm m}\rangle \approx 0.2$ consistent with previous results. However, marginalizing over $0 < K < 1.5$ yields $\langleβ_{\rm m}\rangle = 0.06 ^{+0.11}_{-0.14}$ with $\langle f_{\rm DM}\rangle = 0.35 \pm 0.08$: isotropy is preferred. This isotropy hides the fact that $β_{\rm m}$ is correlated with kinematic features such as counter rotating cores (CRCs), kinematically distinct cores (KDCs), and low-level velocities (LVs): SRs with LVs are likely to be radially biased while SRs with CRCs are likely to be tangentially biased, and SRs with KDCs are intermediate. Existing cosmological simulations allow us to understand these results qualitatively in terms of their dynamical structures and formation histories although there exist quantitative tensions. More realistic cosmological simulations, particularly allowing for $M_\star/L$ gradients, may be required to better understand SRs. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 11 figures, 3 tables, ApJ, accepted for publication

Black hole scaling relations of active and quiescent galaxies: Addressing selection effects and constraining virial factors

Authors: Francesco Shankar, Mariangela Bernardi, Kayleigh Richardson, Christopher Marsden, Ravi K. Sheth, Viola Allevato, Luca Graziani, Mar Mezcua, Federica Ricci, Samantha J. Penny, Fabio La Franca, Fabio Pacucci

Abstract: Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-… ▽ More Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-sigma (velocity dispersion) relation. Local samples of active galactic nuclei (AGN) do not suffer from this selection effect, as in these samples Mbh is estimated from megamasers and/or reverberation mapping-based techniques. With the exception of megamasers, Mbh-estimates in these AGN samples are proportional to a virial coefficient fvir. Direct modelling of the broad line region suggests that fvir~3.5. However, this results in a Mbh-Mstar relation for AGN which lies below and is steeper than the one observed for quiescent black hole samples. A similar though milder trend is seen for the Mbh-sigma relation. Matching the high-mass end of the Mbh-Mstar and Mbh-sigma relations observed in quiescent samples requires fvir~15 and fvir~7, respectively. On the other hand, fvir~3.5 yields Mbh-sigma and Mbh-Mstar relations for AGN which are remarkably consistent with the expected `intrinsic' correlations for quiescent samples (i.e., once account has been made of the angular resolution-related selection effect), providing additional evidence that the sample of local quiescent black holes is biased. We also show that, as is the case for quiescent black holes, the Mbh-Mstar scaling relation of AGN is driven by velocity dispersion, thus providing additional key constraints to black hole-galaxy co-evolution models. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: 15 pages, 5 Figures. MNRAS, accepted

Cosmic distance inference from purely geometric BAO methods: Linear Point standard ruler and Correlation Function Model Fitting

Authors: Stefano Anselmi, Pier-Stefano Corasaniti, Ariel G. Sanchez, Glenn D. Starkman, Ravi K. Sheth, Idit Zehavi

Abstract: Leveraging the Baryon Acoustic Oscillations (BAO) feature present in clustering 2-point statistics, we aim to measure cosmological distances independently of the underlying background cosmological model. However this inference is complicated by late-time non-linearities that introduce model and tracer dependencies in the clustering correlation function and power spectrum, which must be properly ac… ▽ More Leveraging the Baryon Acoustic Oscillations (BAO) feature present in clustering 2-point statistics, we aim to measure cosmological distances independently of the underlying background cosmological model. However this inference is complicated by late-time non-linearities that introduce model and tracer dependencies in the clustering correlation function and power spectrum, which must be properly accounted for. With this in mind, we introduce the "Purely-Geometric-BAO," which provides a rigorous tool to measure cosmological distances without assuming a specific background cosmology. We focus on the 2-point clustering correlation function monopole, and show how to implement such an inference scheme employing two different methodologies: the Linear Point standard ruler (LP) and correlation-function model-fitting (CF-MF). For the first time we demonstrate how, by means of the CF-MF, we can measure very precisely the sound-horizon/isotropic-volume-distance ratio, $r_{d}/D_{V}(\bar{z})$, while correctly propagating all the uncertainties. Using synthetic data, we compare the outcomes of the two methodologies, and find that the LP provides up to $50\%$ more precise measurements than the CF-MF. Finally, we test a procedure widely employed in BAO analyses: fitting the 2-point function while fixing the cosmological and the non-linear-damping parameters at fiducial values. We find that this underestimates the distance errors by nearly a factor of $2$. We thus recommend that this practice be reconsidered, whether for parameter determination or model selection. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: 12 pages, 2 tables

Journal ref: Phys. Rev. D 99, 123515 (2019)

Density and velocity profiles around cosmic voids

Authors: Elena Massara, Ravi K. Sheth

Abstract: We study the evolution of the cross-correlation between voids and the mass density field - i.e. of void profiles. We show that approaches based on the spherical model alone miss an important contribution to the evolution on large scales of most interest to cosmology: they fail to capture the well-known fact that the large-scale bias factor of conserved tracers evolves. We also show that the operat… ▽ More We study the evolution of the cross-correlation between voids and the mass density field - i.e. of void profiles. We show that approaches based on the spherical model alone miss an important contribution to the evolution on large scales of most interest to cosmology: they fail to capture the well-known fact that the large-scale bias factor of conserved tracers evolves. We also show that the operations of evolution and averaging do not commute, but this difference is only significant within about two effective radii. We show how to include a term which accounts for the evolution of bias, which is directly related to the fact that voids move. The void motions are approximately independent of void size, so they are more significant for smaller voids that are typically more numerous. This term also contributes to void-matter pairwise velocities: including it is necessary for modeling the typical outflow speeds around voids. It is, therefore, important for void redshift space distortions. Finally, we show that the excursion set peaks/troughs approach provides a useful, but not perfect framework for describing void profiles and their evolution. △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: 16 pages, 18 figures

Modeling Nearly Spherical Pure-bulge Galaxies with a Stellar Mass-to-Light Ratio Gradient under the $Λ$CDM and MOND Paradigms: I. Methodology, Dynamical Stellar Mass, and Fundamental Mass Plane

Authors: Kyu-Hyun Chae, Mariangela Bernardi, Ravi K. Sheth

Abstract: We carry out spherical Jeans modeling of nearly round pure-bulge galaxies selected from the ATLAS$^{\rm 3D}$ sample. Our modeling allows for gradients in the stellar mass-to-light ratio ($M_\star/L$) through analytic prescriptions parameterized with a `gradient strength' $K$ introduced to accommodate any viable gradient. We use a generalized Osipkov-Merritt model for the velocity dispersion (VD) a… ▽ More We carry out spherical Jeans modeling of nearly round pure-bulge galaxies selected from the ATLAS$^{\rm 3D}$ sample. Our modeling allows for gradients in the stellar mass-to-light ratio ($M_\star/L$) through analytic prescriptions parameterized with a `gradient strength' $K$ introduced to accommodate any viable gradient. We use a generalized Osipkov-Merritt model for the velocity dispersion (VD) anisotropy. We produce Monte Carlo sets of models based on the stellar VD profiles under both the $Λ$CDM and MOND paradigms. Here, we describe the galaxy data, the empirical inputs, and the modeling procedures of obtaining the Monte Carlo sets. We then present the projected dynamical stellar mass, $M_{\rm \star e}$, within the effective radius $R_{\rm e}$, and the fundamental mass plane (FMP) as a function of $K$. We find the scaling of the $K$-dependent mass with respect to the ATLAS$^{\rm 3D}$ reported mass as: $\log_{10} \left[M_{\star{\rm e}}(K)/M_{\star{\rm e}}^{\rm A3D} \right] = a' + b' K$ with $a'=-0.019\pm 0.012$ and $b'=-0.18\pm 0.02$ ($Λ$CDM), or $a'=-0.023\pm 0.014$ and $b'=-0.23\pm 0.03$ (MOND), for $0\le K < 1.5$. The FMP has coefficients consistent with the virial expectation and only the zero-point scales with $K$. The median value of $K$ for the ATLAS$^{\rm 3D}$ galaxies is $\langle K\rangle =0.53^{+0.05}_{-0.04}$. We perform a similar analysis of the much larger SDSS DR7 spectroscopic sample. In this case, only the VD within a single aperture is available, so we impose the additional requirement that the VD slope be similar to that in the ATLAS$^{\rm 3D}$ galaxies. Our analysis of the SDSS galaxies suggests a positive correlation of $K$ with stellar mass. △ Less

Submitted 31 May, 2018; v1 submitted 30 March, 2018; originally announced April 2018.

Comments: ApJ, in press (32 pages in preprint style including 22 figures; Fig 21 corrected)

The Excursion set approach: Stratonovich approximation and Cholesky decomposition

Authors: Farnik Nikakhtar, Mohammadreza Ayromlou, Shant Baghram, Sohrab Rahvar, M. Reza Rahimi Tabar, Ravi K. Sheth

Abstract: The excursion set approach is a framework for estimating how the number density of nonlinear structures in the cosmic web depends on the expansion history of the universe and the nature of gravity. A key part of the approach is the estimation of the first crossing distribution of a suitably chosen barrier by random walks having correlated steps: The shape of the barrier is determined by the physic… ▽ More The excursion set approach is a framework for estimating how the number density of nonlinear structures in the cosmic web depends on the expansion history of the universe and the nature of gravity. A key part of the approach is the estimation of the first crossing distribution of a suitably chosen barrier by random walks having correlated steps: The shape of the barrier is determined by the physics of nonlinear collapse, and the correlations between steps by the nature of the initial density fluctuation field. We describe analytic and numerical methods for calculating such first up-crossing distributions. While the exact solution can be written formally as an infinite series, we show how to approximate it efficiently using the Stratonovich approximation. We demonstrate its accuracy using Monte-Carlo realizations of the walks, which we generate using a novel Cholesky-decomposition based algorithm, which is significantly faster than the algorithm that is currently in the literature. △ Less

Submitted 29 May, 2018; v1 submitted 12 February, 2018; originally announced February 2018.

Comments: 6 pages, 2 figures, Accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 478 (2018) no.4, 5296-5300

The dependence of galaxy clustering on tidal environment in the Sloan Digital Sky Survey

Authors: Aseem Paranjape, Oliver Hahn, Ravi K. Sheth

Abstract: The influence of the Cosmic Web on galaxy formation and evolution is of great observational and theoretical interest. We investigate whether the Cosmic Web leaves an imprint in the spatial clustering of galaxies in the Sloan Digital Sky Survey (SDSS), using the group catalog of Yang et al. and tidal field estimates at $\sim2h^{-1}$Mpc scales from the Mass-Tides-Velocity data set of Wang et al. We… ▽ More The influence of the Cosmic Web on galaxy formation and evolution is of great observational and theoretical interest. We investigate whether the Cosmic Web leaves an imprint in the spatial clustering of galaxies in the Sloan Digital Sky Survey (SDSS), using the group catalog of Yang et al. and tidal field estimates at $\sim2h^{-1}$Mpc scales from the Mass-Tides-Velocity data set of Wang et al. We use the $\textit{tidal anisotropy}$ $α$ (Paranjape et al.) to characterise the tidal environment of groups, and measure the redshift-space 2-point correlation function (2pcf) of group positions and the luminosity- and colour-dependent clustering of group galaxies using samples segregated by $α$. We find that all the 2pcf measurements depend strongly on $α$, with factors of $\sim20$ between the large-scale 2pcf of objects in the most and least isotropic environments. To test whether these strong trends imply `beyond halo mass' effects for galaxy evolution, we compare our results with corresponding 2pcf measurements in mock catalogs constructed using a halo occupation distribution that only uses halo mass as an input. We find that this prescription qualitatively reproduces $\textit{all}$ observed trends, and also quantitatively matches many of the observed results. Although there are some statistically significant differences between our `halo mass only' mocks and the data -- in the most and least isotropic environments -- which deserve further investigation, our results suggest that if the tidal environment induces additional effects on galaxy properties other than those inherited from their host halos, then these must be weak. △ Less

Submitted 16 January, 2018; v1 submitted 14 January, 2018; originally announced January 2018.

Comments: 11 pages, 6 figures, submitted to MNRAS

$M_*/L$ gradients driven by IMF variation: Large impact on dynamical stellar mass estimates

Authors: M. Bernardi, R. K. Sheth, H. Dominguez-Sanchez, J. -L. Fischer, K. -H. Chae, M. Huertas-Company, F. Shankar

Abstract: Within a galaxy the stellar mass-to-light ratio $Υ_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $Υ_*$ gradients than if the IMF is held fixed. If $Υ_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as… ▽ More Within a galaxy the stellar mass-to-light ratio $Υ_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $Υ_*$ gradients than if the IMF is held fixed. If $Υ_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as a factor of two for the most massive galaxies, though stellar population estimates $M_*^{\rm SP}$ are also affected. Large $Υ_*$-gradients have four main consequences: First, $M_*^{\rm dyn}$ cannot be estimated independently of stellar population synthesis models. Second, if there is a lower limit to $Υ_*$ and gradients are unknown, then requiring $M_*^{\rm dyn}=M_*^{\rm SP}$ constrains them. Third, if gradients are stronger in more massive galaxies, then $M_*^{\rm dyn}$ and $M_*^{\rm SP}$ can be brought into agreement, not by shifting $M_*^{\rm SP}$ upwards by invoking constant bottom-heavy IMFs, as advocated by a number of recent studies, but by revising $M_*^{\rm dyn}$ estimates in the literature downwards. Fourth, accounting for $Υ_*$ gradients changes the high-mass slope of the stellar mass function $φ(M_*^{\rm dyn})$, and reduces the associated stellar mass density. These conclusions potentially impact estimates of the need for feedback and adiabatic contraction, so our results highlight the importance of measuring $Υ_*$ gradients in larger samples. △ Less

Submitted 21 March, 2018; v1 submitted 14 December, 2017; originally announced December 2017.

Comments: 13 pages, 7 figures, MNRAS in press

Bimodal Formation Time Distribution for Infall Dark Matter Halos

Authors: Jingjing Shi, Huiyuan Wang, H. J. Mo, Lizhi Xie, Xiaoyu Wang, Andrea Lapi, Ravi K. Sheth

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 give… ▽ More 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. △ Less

Submitted 14 March, 2018; v1 submitted 1 December, 2017; originally announced December 2017.

Comments: 21 pages, 7 figures, accepted for publication by ApJ

Journal ref: Astrophys. J. 857 (2018), 127

The Linear Point Standard Ruler for galaxy survey data: validation with mock catalogues

Authors: Stefano Anselmi, Pier-Stefano Corasaniti, Glenn D. Starkman, Ravi K. Sheth, Idit Zehavi

Abstract: Due to late time non-linearities, the location of the acoustic peak in the two-point galaxy correlation function is a redshift-dependent quantity, thus it cannot be simply employed as a cosmological standard ruler. This has motivated the recent proposal of a novel ruler, also located in the Baryon Acoustic Oscillation range of scales of the correlation function, dubbed the "linear point". Unlike t… ▽ More Due to late time non-linearities, the location of the acoustic peak in the two-point galaxy correlation function is a redshift-dependent quantity, thus it cannot be simply employed as a cosmological standard ruler. This has motivated the recent proposal of a novel ruler, also located in the Baryon Acoustic Oscillation range of scales of the correlation function, dubbed the "linear point". Unlike the peak, it is insensitive at the $0.5\%$ level to many of the non-linear effects that distort the clustering correlation function and shift the peak. However, this is not enough to make the linear point a useful standard ruler. In addition, we require a model-independent method to estimate its value from real data, avoiding the need to deploy a poorly known non-linear model of the correlation function. In this manuscript, we precisely validate a procedure for model-independent estimation of the linear point. We also identify the optimal set-up to estimate the linear point from the correlation function using galaxy catalogs. The methodology developed here is of general validity, and can be applied to any galaxy correlation-function data. As a working example, we apply this procedure to the LOWZ and CMASS galaxy samples of the Twelfth Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS), for which the estimates of cosmic distances using the linear point have been presented in Anselmi et al. (2017) [1]. △ Less

Submitted 24 November, 2017; originally announced November 2017.

Comments: The real data distance constraints presented in the companion paper are improved: arXiv:1703.01275

Journal ref: Phys. Rev. D 98, 023527 (2018)

Stellar mass functions and implications for a variable IMF

Authors: M. Bernardi, R. K. Sheth, J. -L. Fischer, A. Meert, K. -H. Chae, H. Dominguez-Sanchez, M. Huertas-Company, F. Shankar, V. Vikram

Abstract: Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $σ_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with… ▽ More Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $σ_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with large $σ$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, $M_*^{\rm α_{JAM}}$, from $M_{*}$ and $σ_e$ for a sample of $6 \times 10^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R)\,R_e\,σ_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS$^{\rm 3D}$ samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS$^{\rm 3D}$ sample varies weakly with $n$: $(σ_R/σ_e) = (R/R_e)^{-γ(n)}$. The resulting stellar mass functions, based on $M_*^{\rm α_{JAM}}$ and the recalibrated virial mass, are in good agreement. If the $M_*^{\rm α_{JAM}}/M_* - σ_e$ correlation is indeed due to the IMF, and stellar M/L gradients can be ignored, then our $φ(M_*^{\rm α_{JAM}})$ is an estimate of the stellar mass function in which $σ_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $σ_e$ produces comparable results. By demonstrating that cheaper proxies are sufficiently accurate, our analysis should enable a more reliable census of the mass in stars for large galaxy samples, at a fraction of the cost. Our results are provided in tabular form. △ Less

Submitted 15 December, 2017; v1 submitted 19 October, 2017; originally announced October 2017.

Comments: 17 pages, 19 figures, 4 tables. Accepted for publication by MNRAS. Tables 1, C1 and C2 are provided as ancillary files

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

Authors: Kyu-Hyun Chae, Mariangela Bernardi, Ravi K. Sheth, In-Taek Gong

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 (… ▽ More 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 $Λ$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. △ Less

Submitted 12 April, 2019; v1 submitted 25 July, 2017; originally announced July 2017.

Comments: ApJ, accepted (22 pages including 19 figures)

Dependence of halo bias on mass and environment

Authors: Jingjing Shi, Ravi K. Sheth

Abstract: The simplest analyses of halo bias assume that halo mass alone determines halo clustering. However, if the large scale environment is fixed, then halo clustering is almost entirely determined by environment, and is almost completely independent of halo mass. We show why. Our analysis is useful for studies which use the environmental dependence of clustering to constrain cosmological and galaxy for… ▽ More The simplest analyses of halo bias assume that halo mass alone determines halo clustering. However, if the large scale environment is fixed, then halo clustering is almost entirely determined by environment, and is almost completely independent of halo mass. We show why. Our analysis is useful for studies which use the environmental dependence of clustering to constrain cosmological and galaxy formation models. It also shows why many correlations between galaxy properties and environment are merely consequences of the underlying correlations between halos and their environments, and provides a framework for quantifying such inherited correlations. △ Less

Submitted 31 August, 2017; v1 submitted 13 July, 2017; originally announced July 2017.

Comments: 9 pages, 6 figures, accepted by MNRAS

Halo assembly bias and the tidal anisotropy of the local halo environment

Authors: Aseem Paranjape, Oliver Hahn, Ravi K. Sheth

Abstract: We study the role of the local tidal environment in determining the assembly bias of dark matter haloes. Previous results suggest that the anisotropy of a halo's environment (i.e, whether it lies in a filament or in a more isotropic region) can play a significant role in determining the eventual mass and age of the halo. We statistically isolate this effect using correlations between the large-sca… ▽ More We study the role of the local tidal environment in determining the assembly bias of dark matter haloes. Previous results suggest that the anisotropy of a halo's environment (i.e, whether it lies in a filament or in a more isotropic region) can play a significant role in determining the eventual mass and age of the halo. We statistically isolate this effect using correlations between the large-scale and small-scale environments of simulated haloes at $z=0$ with masses between $10^{11.6}\lesssim (m/h^{-1}M_{\odot})\lesssim10^{14.9}$. We probe the large-scale environment using a novel halo-by-halo estimator of linear bias. For the small-scale environment, we identify a variable $α_R$ that captures the $\textit{tidal anisotropy}$ in a region of radius $R=4R_{\textrm{200b}}$ around the halo and correlates strongly with halo bias at fixed mass. Segregating haloes by $α_R$ reveals two distinct populations. Haloes in highly isotropic local environments ($α_R\lesssim0.2$) behave as expected from the simplest, spherically averaged analytical models of structure formation, showing a $\textit{negative}$ correlation between their concentration and large-scale bias at $\textit{all}$ masses. In contrast, haloes in anisotropic, filament-like environments ($α_R\gtrsim0.5$) tend to show a $\textit{positive}$ correlation between bias and concentration at any mass. Our multi-scale analysis cleanly demonstrates how the overall assembly bias trend across halo mass emerges as an average over these different halo populations, and provides valuable insights towards building analytical models that correctly incorporate assembly bias. We also discuss potential implications for the nature and detectability of galaxy assembly bias. △ Less

Submitted 25 February, 2018; v1 submitted 29 June, 2017; originally announced June 2017.

Comments: 19 pages, 15 figures; v2: revised in response to referee comments, added references and discussion, conclusions unchanged. Accepted in MNRAS