Source author record

Ravi K. Sheth

Ravi K. Sheth appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph hep-th gr-qc astro-ph.GA hep-ph astro-ph.HE physics.data-an

Catalog footprint

What is connected

56works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Protohalos and their connection to halo assembly, shape and structure

Protohalos, primordial regions in the initial cosmic density field that evolve into dark matter halos, are crucial for understanding cosmic structure formation. Motivated by the potential to reconstruct protohalo positions and shapes from observed galaxies using a novel approach grounded in optimal transport theory, we revisit the relationship between the structural properties of protohalos and the assembly histories, concentrations, and final morphologies of their associated dark matter halos. To better understand halo assembly, we introduce a new estimator defined by an integral over redshifts and compare its performance to $z_{50}$, the commonly used redshift at which half of the final halo mass is formed. We quantify protohalo structure using the three invariants of the inertia, deformation, and energy shear tensors. Although past research has correlated the first two invariants of the deformation and energy tensors with halo formation, our findings reveal that the third invariant also significantly correlates with halo assembly and final shape.

preprint2023arXiv

Optimal Transport Reconstruction of Baryon Acoustic Oscillations

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$.

preprint2021arXiv

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

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.

preprint2021arXiv

The smearing scale in Laguerre reconstructions of the correlation function

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.

preprint2020arXiv

A universal threshold for primordial black hole formation

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$.

preprint2020arXiv

Analytical thresholds for black hole formation in general cosmological backgrounds

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]$.

preprint2020arXiv

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

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.

preprint2019arXiv

Linear Point and Sound Horizon as Purely Geometric standard rulers

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.

preprint2019arXiv

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

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.

preprint2019arXiv

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

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.

preprint2016arXiv

Excursion set peaks: the role of shear

Recent analytical work on the modelling of dark halo abundances and clustering has demonstrated the advantages of combining the excursion set approach with peaks theory. We extend these ideas and introduce a model of excursion set peaks that incorporates the role of initial tidal effects or shear in determining the gravitational collapse of dark haloes. The model -- in which the critical density threshold for collapse depends on the tidal influences acting on protohaloes -- is well motivated from ellipsoidal collapse arguments and is also simple enough to be analytically tractable. We show that the predictions of this model are in very good agreement with measurements of the halo mass function and traditional scale dependent halo bias in N-body simulations across a wide range of masses and redshift. The presence of shear in the collapse threshold means that halo bias is naturally predicted to be nonlocal, and that protohalo densities at fixed mass are naturally predicted to have Lognormal-like distributions. We present the first direct estimate of Lagrangian nonlocal bias in N-body simulations, finding broad agreement with the model prediction. Finally, the simplicity of the model (which has essentially a single free parameter) opens the door to building efficient and accurate non-universal fitting functions of halo abundances and bias for use in precision cosmology.

preprint2015arXiv

The universality of the virial halo mass function and models for non-universality of other halo definitions

The abundance of galaxy clusters can constrain both the geometry and growth of structure in our Universe. However, this probe could be significantly complicated by recent claims of nonuniversality -- non-trivial dependences with respect to the cosmological model and redshift. In this work we analyse the dependance of the mass function on the way haloes are identified and establish if this can cause departures from universality. In order to explore this dependance, we use a set of different N-body cosmological simulations (Le SBARBINE simulations), with the latest cosmological parameters from the Planck collaboration; this first suite of simulations is followed by a lower resolution set, carried out with different cosmological parameters. We identify dark matter haloes using a Spherical Overdensity algorithm with varying overdensity thresholds (virial, 2000rho_c, 1000rho_c, 500rho_c, 200rho_c and 200rho_b) at all redshifts. We notice that, when expressed in term of the rescaled variable nu, the mass functionfor virial haloes is a nearly universal as a function of redshift and cosmology, while this is clearly not the case for the other overdensities we considered. We provide fitting functions for the halo mass function parameters as a function of overdensity, that allow to predict, to within a few percent accuracy, the halo mass function for a wide range of halo definitions, redshifts and cosmological models. We then show how the departures from universality associated with other halo definitions can be derived by combining the universality of the virial definition with the expected shape of the density profile of halos.

preprint2014arXiv

On the Markovian assumption in the excursion set approach: the approximation of Markov Velocities

The excursion set approach uses the statistics of the density field smoothed on a wide range of scales, to gain insight into a number of interesting processes in nonlinear structure formation, such as cluster assembly, merging and clustering. The approach treats the curve defined by the height of the overdensity fluctuation field when changing the smoothing scale as a random walk. The steps of the walks are often assumed to be uncorrelated, so that the walk heights are a Markov process, even though this assumption is known to be inaccurate for physically relevant filters. We develop a model in which the walk steps, rather than heights, are a Markov process, and correlations between steps arise because of nearest neighbour interactions. This model is a particular case of a general class, which we call Markov Velocity models. We show how these can approximate the walks generated by arbitrary power spectra and filters, and, unlike walks with Markov heights, provide a very good approximation to physically relevant models. We define a Markov Velocity Monte Carlo algorithm to generate walks whose first crossing distribution is very similar to that of TopHat-smoothed LCDM walks. Finally, we demonstrate that Markov Velocity walks generically exhibit a simple but realistic form of assembly bias, so we expect them to be useful in the construction of more realistic merger history trees.

preprint2014arXiv

Stochasticity in halo formation and the excursion set approach

The simplest stochastic halo formation models assume that the traceless part of the shear field acts to increase the initial overdensity (or decrease the underdensity) that a protohalo (or protovoid) must have if it is to form by the present time. Equivalently, it is the difference between the overdensity and (the square root of the) shear that must be larger than a threshold value. To estimate the effect this has on halo abundances using the excursion set approach, we must solve for the first crossing distribution of a barrier of constant height by the random walks associated with the difference, which is now (even for Gaussian initial conditions) a non-Gaussian variate. The correlation properties of such non-Gaussian walks are inherited from those of the density and the shear, and, since they are independent processes, the solution is in fact remarkably simple. We show that this provides an easy way to understand why earlier heuristic arguments about the nature of the solution worked so well. In addition to modelling halos and voids, this potentially simplifies models of the abundance and spatial distribution of filaments and sheets in the cosmic web.

preprint2014arXiv

The excursion set approach in non-Gaussian random fields

Insight into a number of interesting questions in cosmology can be obtained from the first crossing distributions of physically motivated barriers by random walks with correlated steps. We write the first crossing distribution as a formal series, ordered by the number of times a walk upcrosses the barrier. Since the fraction of walks with many upcrossings is negligible if the walk has not taken many steps, the leading order term in this series is the most relevant for understanding the massive objects of most interest in cosmology. This first term only requires knowledge of the bivariate distribution of the walk height and slope, and provides an excellent approximation to the first crossing distribution for all barriers and smoothing filters of current interest. We show that this simplicity survives when extending the approach to the case of non-Gaussian random fields. For non-Gaussian fields which are obtained by deterministic transformations of a Gaussian, the first crossing distribution is simply related to that for Gaussian walks crossing a suitably rescaled barrier. Our analysis shows that this is a useful way to think of the generic case as well. Although our study is motivated by the possibility that the primordial fluctuation field was non-Gaussian, our results are general. In particular, they do not assume the non-Gaussianity is small, so they may be viewed as the solution to an excursion set analysis of the late-time, nonlinear fluctuation field rather than the initial one. They are also useful for models in which the barrier height is determined by quantities other than the initial density, since most other physically motivated variables (such as the shear) are usually stochastic and non-Gaussian. We use the Lognormal transformation to illustrate some of our arguments.

preprint2013arXiv

Bias deconstructed: Unravelling the scale dependence of halo bias using real space measurements

We explore the scale dependence of halo bias using real space cross-correlation measurements in N-body simulations and in Pinocchio, an algorithm based on Lagrangian Perturbation Theory. Recent work has shown how to interpret such real space measurements in terms of k-dependent bias in Fourier space, and how to remove the k-dependence to reconstruct the k-independent peak-background split halo bias parameters. We compare our reconstruction of the linear bias, which requires no free parameters, with previous estimates from N-body simulations which were obtained directly in Fourier space at large scales, and find very good agreement. Our reconstruction of the quadratic bias is similarly parameter-free, although in this case there are no previous Fourier space measurements to compare with. Our analysis of N-body simulations explicitly tests the predictions of the excursion set peaks (ESP) formalism of Paranjape et al. (2013) for the scale dependence of bias; we find that the ESP predictions accurately describe our measurements. In addition, our measurements in Pinocchio serve as a useful, successful consistency check between Pinocchio and N-body simulations that is not accessible to traditional measurements.

preprint2013arXiv

Ellipsoidal halo finders and implications for models of triaxial halo formation

We describe an algorithm for identifying ellipsoidal haloes in numerical simulations, and quantify how the resulting estimates of halo mass and shape differ with respect to spherical halo finders. Haloes become more prolate when fit with ellipsoids, the difference being most pronounced for the more aspherical objects. Although the ellipsoidal mass is systematically larger, this is less than 10% for most of the haloes. However, even this small difference in mass corresponds to a significant difference in shape. We quantify these effects also on the initial mass and deformation tensors, on which most models of triaxial collapse are based. By studying the properties of protohaloes in the initial conditions, we find that models in which protohaloes are identified in Lagrangian space by three positive eigenvalues of the deformation tensor are tenable only at the masses well-above $M_*$. The overdensity $δ$ within almost any protohalo is larger than the critical value associated with spherical collapse (increasing as mass decreases); this is in good qualitative agreement with models which identify haloes requiring that collapse have occured along all three principal axes, each axis having turned around from the universal expansion at a different time. The distributions of initial values are in agreement with the simplest predictions associated with ellipsoidal collapse, assuming initially spherical protohaloes, collapsed around random positions which were sufficiently overdense. However, most protohaloes are not spherical and departures from sphericity increase as protohalo mass decreases. [Abridged]

preprint2013arXiv

Excursion set peaks: a self-consistent model of dark halo abundances and clustering

We describe how to extend the excursion set peaks framework so that its predictions of dark halo abundances and clustering can be compared directly with simulations. These extensions include: a halo mass definition which uses the TopHat filter in real space; the mean dependence of the critical density for collapse delta_c on halo mass m; and the scatter around this mean value. All three of these are motivated by the physics of triaxial rather than spherical collapse. A comparison of the resulting mass function with N-body results shows that, if one uses delta_c(m) and its scatter as determined from simulations, then all three are necessary ingredients for obtaining ~10% accuracy. E.g., assuming a constant value of delta_c with no scatter, as motivated by the physics of spherical collapse, leads to many more massive halos than seen in simulations. The same model is also in excellent agreement with N-body results for the linear halo bias, especially at the high mass end where the traditional peak-background split argument applied to the mass function fit is known to underpredict the measured bias by ~10%. In the excursion set language, our model is about walks centered on special positions (peaks) in the initial conditions -- we discuss what it implies for the usual calculation in which all walks contribute to the statistics.

preprint2013arXiv

Lagrangian bias in the local bias model

It is often assumed that the halo-patch fluctuation field can be written as a Taylor series in the initial Lagrangian dark matter density fluctuation field. We show that if this Lagrangian bias is local, and the initial conditions are Gaussian, then the two-point cross-correlation between halos and mass should be linearly proportional to the mass-mass auto-correlation function. This statement is exact and valid on all scales; there are no higher order contributions, e.g., from terms proportional to products or convolutions of two-point functions, which one might have thought would appear upon truncating the Taylor series of the halo bias function. In addition, the auto-correlation function of locally biased tracers can be written as a Taylor series in the auto-correlation function of the mass; there are no terms involving, e.g., derivatives or convolutions. Moreover, although the leading order coefficient, the linear bias factor of the auto-correlation function is just the square of that for the cross-correlation, it is the same as that obtained from expanding the mean number of halos as a function of the local density only in the large-scale limit. In principle, these relations allow simple tests of whether or not halo bias is indeed local in Lagrangian space. We discuss why things are more complicated in practice. We also discuss our results in light of recent work on the renormalizability of halo bias, demonstrating that it is better to renormalize than not. We use the Lognormal model to illustrate many of our findings.

preprint2013arXiv

Measures of galaxy environment -- II. Rank-ordered mark correlations

We analyze environmental correlations using mark clustering statistics with the mock galaxy catalogue constructed by Muldrew et al. (Paper I). We find that mark correlation functions are able to detect even a small dependence of galaxy properties on the environment, quantified by the overdensity $1+δ$, while such a small dependence would be difficult to detect by traditional methods. We then show that rank ordering the marks and using the rank as a weight is a simple way of comparing the correlation signals for different marks. With this we quantify to what extent fixed-aperture overdensities are sensitive to large-scale halo environments, nearest-neighbor overdensities are sensitive to small-scale environments within haloes, and colour is a better tracer of overdensity than is luminosity.

preprint2013arXiv

Non-local Lagrangian bias

Halos are biased tracers of the dark matter distribution. It is often assumed that the patches from which halos formed are locally biased with respect to the initial fluctuation field, meaning that the halo-patch fluctuation field can be written as a Taylor series in that of the dark matter. If quantities other than the local density influence halo formation, then this Lagrangian bias will generically be nonlocal; the Taylor series must be performed with respect to these other variables as well. We illustrate the effect with Monte-Carlo simulations of a model in which halo formation depends on the local shear (the quadrupole of perturbation theory), and provide an analytic model which provides a good description of our results. Our model, which extends the excursion set approach to walks in more than one dimension, works both when steps in the walk are uncorrelated, as well as when there are correlations between steps. For walks with correlated steps, our model includes two distinct types of nonlocality: one is due to the fact that the initial density profile around a patch which is destined to form a halo must fall sufficiently steeply around it -- this introduces k-dependence to even the linear bias factor, but otherwise only affects the monopole of the clustering signal. The other is due to the surrounding shear field; this affects the quadratic and higher order bias factors, and introduces an angular dependence to the clustering signal. In both cases, our analysis shows that these nonlocal Lagrangian bias terms can be significant, particularly for massive halos; they must be accounted for in analyses of higher order clustering such as the halo bispectrum in Lagrangian or Eulerian space. Although we illustrate these effects using halos, our analysis and conclusions also apply to the other constituents of the cosmic web -- filaments, sheets and voids.

preprint2013arXiv

Peak-peak correlations in the cosmic background radiation from cosmic strings

We examine the two-point correlation function of local maxima in temperature fluctuations at the last scattering surface when this stochastic field is modified by the additional fluctuations produced by straight cosmic strings via the Kaiser-Stebbins effect. We demonstrate that one can detect the imprint of cosmic strings with tension $Gμ\gtrsim 1.2 \times 10^{-8}$ on noiseless $1^\prime$ resolution cosmic microwave background (CMB) maps at 95% confidence interval. Including the effects of foregrounds and anticipated systematic errors increases the lower bound to $Gμ\gtrsim 9.0\times 10^{-8}$ at $2σ$ confidence level. Smearing by beams of order 4' degrades the bound further to $Gμ\gtrsim 1.6 \times 10^{-7}$. Our results indicate that two-point statistics are more powerful than 1-point statistics (e.g. number counts) for identifying the non-Gaussianity in the CMB due to straight cosmic strings.

preprint2013arXiv

Stochastic bias in multi-dimensional excursion set approaches

We describe a simple fully analytic model of the excursion set approach associated with two Gaussian random walks: the first walk represents the initial overdensity around a protohalo, and the second is a crude way of allowing for other factors which might influence halo formation. This model is richer than that based on a single walk, because it yields a distribution of heights at first crossing. We provide explicit expressions for the unconditional first crossing distribution which is usually used to model the halo mass function, the progenitor distributions, and the conditional distributions from which correlations with environment are usually estimated. These latter exhibit perhaps the simplest form of what is often called nonlocal bias, and which we prefer to call stochastic bias, since the new bias effects arise from `hidden-variables' other than density, but these may still be defined locally. We provide explicit expressions for these new bias factors. We also provide formulae for the distribution of heights at first crossing in the unconditional and conditional cases. In contrast to the first crossing distribution, these are exact, even for moving barriers, and for walks with correlated steps. The conditional distributions yield predictions for the distribution of halo concentrations at fixed mass and formation redshift. They also exhibit assembly bias like effects, even when the steps in the walks themselves are uncorrelated. Finally, we show how the predictions are modified if we add the requirement that halos form around peaks: these depend on whether the peaks constraint is applied to a combination of the overdensity and the other variable, or to the overdensity alone. Our results demonstrate the power of requiring models to reproduce not just halo counts but the distribution of overdensities at fixed protohalo mass as well.

preprint2013arXiv

The importance of stepping up in the excursion set approach

Recently, we provided a simple but accurate formula which closely approximates the first crossing distribution associated with random walks having correlated steps. The approximation is accurate for the wide range of barrier shapes of current interest and is based on the requirement that, in addition to having the right height, the walk must cross the barrier going upwards. Therefore, it only requires knowledge of the bivariate distribution of the walk height and slope, and is particularly useful for excursion set models of the massive end of the halo mass function. However, it diverges at lower masses. We show how to cure this divergence by using a formulation which requires knowledge of just one other variable. While our analysis is general, we use examples based on Gaussian initial conditions to illustrate our results. Our formulation, which is simple and fast, yields excellent agreement with the considerably more computationally expensive Monte-Carlo solution of the first crossing distribution, for a wide variety of moving barriers, even at very low masses.

preprint2012arXiv

Formation times, mass growth histories and concentrations of dark matter haloes

We develop a simple model for estimating the mass growth histories of dark matter halos. The model is based on a fit to the formation time distribution, where formation is defined as the earliest time that the main branch of the merger tree contains a fraction f of the final mass M . Our analysis exploits the fact that the median formation time as a function of f is the same as the median of the main progenitor mass distribution as a function of time. When coupled with previous work showing that the concentration c of the final halo is related to the formation time tf associated with f \sim 0.04, our approach provides a simple algorithm for estimating how the distribution of halo concentrations may be expected to depend on mass, redshift and the expansion history of the background cosmology. We also show that one can predict log10 c with a precision of about 0.13 and 0.09 dex if only its mass, or both mass and tf are known. And, conversely, one can predict log10 tf from mass or c with a precision of 0.12 and 0.09 dex, approximately independent of f. Adding the mass to the c-based estimate does not result in further improvement. These latter results may be useful for studies which seek to compare the age of the stars in the central galaxy in a halo with the time the core was first assembled.

preprint2012arXiv

Gravity and Large-Scale Non-local Bias

The relationship between galaxy and matter overdensities, bias, is most often assumed to be local. This is however unstable under time evolution, we provide proofs under several sets of assumptions. In the simplest model galaxies are created locally and linearly biased at a single time, and subsequently move with the matter (no velocity bias) conserving their comoving number density (no merging). We show that, after this formation time, the bias becomes unavoidably non-local and non-linear at large scales. We identify the non-local gravitationally induced fields in which the galaxy overdensity can be expanded, showing that they can be constructed out of the invariants of the deformation tensor (Galileons). In addition, we show that this result persists if we include an arbitrary evolution of the comoving number density of tracers. We then include velocity bias, and show that new contributions appear, a dipole field being the signature at second order. We test these predictions by studying the dependence of halo overdensities in cells of fixed matter density: measurements in simulations show that departures from the mean bias relation are strongly correlated with the non-local gravitationally induced fields identified by our formalism. The effects on non-local bias seen in the simulations are most important for the most biased halos, as expected from our predictions. The non-locality seen in the simulations is not fully captured by assuming local bias in Lagrangian space. Accounting for these effects when modeling galaxy bias is essential for correctly describing the dependence on triangle shape of the galaxy bispectrum, and hence constraining cosmological parameters and primordial non-Gaussianity. We show that using our formalism we remove an important systematic in the determination of bias parameters from the galaxy bispectrum, particularly for luminous galaxies. (abridged)

preprint2012arXiv

On the anisotropic density distribution on large scales

Motivated by the recent detection of an enhanced clustering signal along the major axis of haloes in N-body simulations, we derive a formula for the anisotropic density distribution around haloes and voids on large scales. Our model, which assumes linear theory and that the formation and orientation of nonlinear structures are strongly correlated with the Lagrangian shear, is in good agreement with measurements. We also show that the measured amplitude is inconsistent with a model in which the alignment is produced by the initial inertia rather than shear tensor.

preprint2012arXiv

One step beyond: The excursion set approach with correlated steps

We provide a simple formula that accurately approximates the first crossing distribution of barriers having a wide variety of shapes, by random walks with a wide range of correlations between steps. Special cases of it are useful for estimating halo abundances, evolution, and bias, as well as the nonlinear counts in cells distribution. We discuss how it can be extended to allow for the dependence of the barrier on quantities other than overdensity, to construct an excursion set model for peaks, and to show why assembly and scale dependent bias are generic even at the linear level.

preprint2012arXiv

Peaks theory and the excursion set approach

We describe a model of dark matter halo abundances and clustering which combines the two most widely used approaches to this problem: that based on peaks and the other based on excursion sets. Our approach can be thought of as addressing the cloud-in-cloud problem for peaks and/or modifying the excursion set approach so that it averages over a special subset, rather than all possible walks. In this respect, it seeks to account for correlations between steps in the walk as well as correlations between walks. We first show how the excursion set and peaks models can be written in the same formalism, and then use this correspondence to write our combined excursion set peaks model. We then give simple expressions for the mass function and bias, showing that even the linear halo bias factor is predicted to be k-dependent as a consequence of the nonlocality associated with the peak constraint. At large masses, our model has little or no need to rescale the variable delta_c from the value associated with spherical collapse, and suggests a simple explanation for why the linear halo bias factor appears to lie above that based on the peak-background split at high masses when such a rescaling is assumed. Although we have concentrated on peaks, our analysis is more generally applicable to other traditionally single-scale analyses of large-scale structure.

preprint2012arXiv

Plain fundamentals of Fundamental Planes: Analytics and algorithms

The coefficients a and b of the Fundamental Plane relation R ~ Sigma^a I^b depend on whether one minimizes the scatter in the R direction or orthogonal to the Plane. We provide explicit expressions for a and b (and confidence limits) in terms of the covariances between logR, logSigma and logI. Our analysis is more generally applicable to any other correlations between three variables: e.g., the color-magnitude-Sigma relation, the L-Sigma-Mbh relation, or the relation between the X-ray luminosity, Sunyaev-Zeldovich decrement and optical richness of a cluster, so we provide IDL code which implements these ideas, and we show how our analysis generalizes further to correlations between more than three variables. We show how to account for correlated errors and selection effects, and quantify the difference between the direct, inverse and orthogonal fit coefficients. We show that the three vectors associated with the Fundamental Plane can all be written as simple combinations of a and b because the distribution of I is much broader than that of Sigma, and Sigma and I are only weakly correlated. Why this should be so for galaxies is a fundamental open question about the physics of early-type galaxy formation. If luminosity evolution is differential, and Rs and Sigmas do not evolve, then this is just an accident: Sigma and I must have been correlated in the past. On the other hand, if the (lack of) correlation is similar to that at the present time, then differential luminosity evolution must have been accompanied by structural evolution. A model in which the luminosities of low-L galaxies evolve more rapidly than do those of higher-L galaxies is able to produce the observed decrease in a (by a factor of 2 at z~1) while having b decrease by only about 20 percent. In such a model, the Mdyn/L ratio is a steeper function of Mdyn at higher z.

preprint2012arXiv

Scale dependent halo bias in the excursion set approach

If one accounts for correlations between scales, then nonlocal, k-dependent halo bias is part and parcel of the excursion set approach, and hence of halo model predictions for galaxy bias. We present an analysis that distinguishes between a number of different effects, each one of which contributes to scale-dependent bias in real space. We show how to isolate these effects and remove the scale dependence, order by order, by cross-correlating the halo field with suitably transformed versions of the mass field. These transformations may be thought of as simple one-point, two-scale measurements that allow one to estimate quantities which are usually constrained using n-point statistics. As part of our analysis, we present a simple analytic approximation for the first crossing distribution of walks with correlated steps which are constrained to pass through a specified point, and demonstrate its accuracy. Although we concentrate on nonlinear, nonlocal bias with respect to a Gaussian random field, we show how to generalize our analysis to more general fields.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: a large sample of mock galaxy catalogues

We present a fast method of producing mock galaxy catalogues that can be used to compute covariance matrices of large-scale clustering measurements and test the methods of analysis. Our method populates a 2nd-order Lagrangian Perturbation Theory (2LPT) matter field, where we calibrate masses of dark matter halos by detailed comparisons with N-body simulations. We demonstrate the clustering of halos is recovered at ~10 per cent accuracy. We populate halos with mock galaxies using a Halo Occupation Distribution (HOD) prescription, which has been calibrated to reproduce the clustering measurements on scales between 30 and 80 Mpc/h. We compare the sample covariance matrix from our mocks with analytic estimates, and discuss differences. We have used this method to make catalogues corresponding to Data Release 9 of the Baryon Oscillation Spectroscopic Survey (BOSS),producing 600 mock catalogues of the "CMASS" galaxy sample. These mocks enabled detailed tests of methods and errors that formed an integral part of companion analyses of these galaxy data.

preprint2012arXiv

The luminosities of the brightest cluster galaxies and brightest satellites in SDSS groups

We show that the distribution of luminosities of Brightest Cluster Galaxies in an SDSS-based group catalog suggests that BCG luminosities are just the statistical extremes of the group galaxy luminosity function. This latter happens to be very well approximated by the all-galaxy luminosity function (restricted to Mr<-19.9), provided one uses a parametrization of this function that is accurate at the bright end. A similar analysis of the luminosity distribution of the Brightest Satellite Galaxies suggests that they are best thought of as being the second brightest pick from the same luminosity distribution of which BCGs are the brightest. I.e., BSGs are not the brightest of some universal satellite luminosity function, in contrast to what Halo Model analyses of the luminosity dependence of clustering suggest. However, we then use mark correlations to provide a novel test of these order statistics, showing that the hypothesis of a universal luminosity function (i.e. no halo mass dependence) from which the BCGs and BSGs are drawn is incompatible with the data, despite the fact that there was no hint of this in the BCG and BSG luminosity distributions themselves. We also discuss why, since extreme value statistics are explicitly a function of the number of draws, the consistency of BCG luminosities with extreme value statistics is most clearly seen if one is careful to perform the test at fixed group richness N. Tests at, e.g., fixed total group luminosity Ltot, will generally be biased and may lead to erroneous conclusions.

preprint2011arXiv

A hierarchy of voids: More ado about nothing

We extend earlier work on the problem of estimating the void-volume function -- the abundance and evolution of large voids which grow gravitationally in an expanding universe -- in two ways. The first removes an ambiguity about how the void-in-cloud process, which erases small voids, should be incorporated into the excursion set approach. The main technical change here is to think of voids within a fully Eulerian, rather than purely Lagrangian, framework. The second accounts for correlations between different spatial scales in the initial conditions. We provide numerical and analytical arguments showing how and why both changes modify the predicted abundances substantially. In particular, we show that the predicted importance of the void-in-cloud process depends strongly on whether or not one accounts for correlations between scales. With our new formulation, the void-in-cloud process dramatically reduces the predicted abundances of voids if such correlations are ignored, but only matters for the smallest voids in the more realistic case in which the spatial correlations are included.

preprint2011arXiv

Halo abundances and counts-in-cells: The excursion set approach with correlated steps

The Excursion Set approach has been used to make predictions for a number of interesting quantities in studies of nonlinear hierarchical clustering. These include the halo mass function, halo merger rates, halo formation times and masses, halo clustering, analogous quantities for voids, and the distribution of dark matter counts in randomly placed cells. The approach assumes that all these quantities can be mapped to problems involving the first crossing distribution of a suitably chosen barrier by random walks. Most analytic expressions for these distributions ignore the fact that, although different k-modes in the initial Gaussian field are uncorrelated, this is not true in real space: the values of the density field at a given spatial position, when smoothed on different real-space scales, are correlated in a nontrivial way. As a result, the problem is to estimate first crossing distribution by random walks having correlated rather than uncorrelated steps. In 1990, Peacock & Heavens presented a simple approximation for the first crossing distribution of a single barrier of constant height by walks with correlated steps. We show that their approximation can be thought of as a correction to the distribution associated with what we call smooth completely correlated walks. We then use this insight to extend their approach to treat moving barriers, as well as walks that are constrained to pass through a certain point before crossing the barrier. For the latter, we show that a simple rescaling, inspired by bivariate Gaussian statistics, of the unconditional first crossing distribution, accurately describes the conditional distribution, independently of the choice of analytical prescription for the former. In all cases, comparison with Monte-Carlo solutions of the problem shows reasonably good agreement. (Abridged)

preprint2011arXiv

Halo bias in the excursion set approach with correlated steps

In the Excursion Set approach, halo abundances and clustering are closely related. This relation is exploited in many modern methods which seek to constrain cosmological parameters on the basis of the observed spatial distribution of clusters. However, to obtain analytic expressions for these quantities, most Excursion Set based predictions ignore the fact that, although different k-modes in the initial Gaussian field are uncorrelated, this is not true in real space: the values of the density field at a given spatial position, when smoothed on different real-space scales, are correlated in a nontrivial way. We show that when the excursion set approach is extended to include such correlations, then one must be careful to account for the fact that the associated prediction for halo bias is explicitly a real-space quantity. Therefore, care must be taken when comparing the predictions of this approach with measurements in simulations, which are typically made in Fourier-space. We show how to correct for this effect, and demonstrate that ignorance of this effect in recent analyses of halo bias has led to incorrect conclusions and biased constraints.

preprint2011arXiv

Modelling the shapes of the largest gravitationally bound objects

We combine the physics of the ellipsoidal collapse model with the excursion set theory to study the shapes of dark matter halos. In particular, we develop an analytic approximation to the nonlinear evolution that is more accurate than the Zeldovich approximation; we introduce a planar representation of halo axis ratios, which allows a concise and intuitive description of the dynamics of collapsing regions and allows one to relate the final shape of a halo to its initial shape; we provide simple physical explanations for some empirical fitting formulae obtained from numerical studies. Comparison with simulations is challenging, as there is no agreement about how to define a non-spherical gravitationally bound object. Nevertheless, we find that our model matches the conditional minor-to-intermediate axis ratio distribution rather well, although it disagrees with the numerical results in reproducing the minor-to-major axis ratio distribution. In particular, the mass dependence of the minor-to-major axis distribution appears to be the opposite to what is found in many previous numerical studies, where low-mass halos are preferentially more spherical than high-mass halos. In our model, the high-mass halos are predicted to be more spherical, consistent with results based on a more recent and elaborate halo finding algorithm, and with observations of the mass dependence of the shapes of early-type galaxies. We suggest that some of the disagreement with some previous numerical studies may be alleviated if we consider only isolated halos.

preprint2011arXiv

Scale-dependent halo bias from scale-dependent growth

We derive a general expression for the large-scale halo bias, in theories with a scale-dependent linear growth, using the excursion set formalism. Such theories include modified gravity models, and models in which the dark energy clustering is non-negligible. A scale dependence is imprinted in both the formation and evolved biases by the scale-dependent growth. Mergers are accounted for in our derivation, which thus extends earlier work which focused on passive evolution. There is a simple analytic form for the bias for those theories in which the nonlinear collapse of perturbations is approximately the same as in general relativity. As an illustration, we apply our results to a simple Yukawa modification of gravity, and use SDSS measurements of the clustering of luminous red galaxies to constrain the theory's parameters.

preprint2010arXiv

Convolution and deconvolution based estimates of galaxy scaling relations from photometric redshift surveys

In addition to the maximum likelihood approach, there are two other methods which are commonly used to reconstruct the true redshift distribution from photometric redshift datasets: one uses a deconvolution method, and the other a convolution. We show how these two techniques are related, and how this relationship can be extended to include the study of galaxy scaling relations in photometric datasets. We then show what additional information photometric redshift algorithms must output so that they too can be used to study galaxy scaling relations, rather than just redshift distributions. We also argue that the convolution based approach may permit a more efficient selection of the objects for which calibration spectra are required.

preprint2010arXiv

Halo model description of the non-linear dark matter power spectrum at $k \gg 1$ Mpc$^{-1}$

Accurate knowledge of the non-linear dark-matter power spectrum is important for understanding the large-scale structure of the Universe, the statistics of dark-matter haloes and their evolution, and cosmological gravitational lensing. We analytically model the dark-matter power spectrum and its cross-power spectrum with dark-matter haloes. Our model extends the halo-model formalism, including realistic substructure population within individual dark-matter haloes and the scatter of the concentration parameter at fixed halo mass. We consider three prescriptions for the mass-concentration relation and two for the substructure distribution in dark-matter haloes. We show that this extension of the halo model mainly increases the predicted power on the small scales, and is crucial for proper modeling the cosmological weak-lensing signal due to low-mass haloes. Our extended formalism shows how the halo model approach can be improved in accuracy as one increases the number of ingredients that are calibrated from n-body simulations.

preprint2010arXiv

Modeling scale-dependent bias on the baryonic acoustic scale with the statistics of peaks of Gaussian random fields

Models of galaxy and halo clustering commonly assume that the tracers can be treated as a continuous field locally biased with respect to the underlying mass distribution. In the peak model pioneered by BBKS, one considers instead density maxima of the initial, Gaussian mass density field as an approximation to the formation site of virialized objects. In this paper, the peak model is extended in two ways to improve its predictive accuracy. Firstly, we derive the two-point correlation function of initial density peaks up to second order and demonstrate that a peak-background split approach can be applied to obtain the k-independent and k-dependent peak bias factors at all orders. Secondly, we explore the gravitational evolution of the peak correlation function within the Zel'dovich approximation. We show that the local (Lagrangian) bias approach emerges as a special case of the peak model, in which all bias parameters are scale-independent and there is no statistical velocity bias. We apply our formulae to study how the Lagrangian peak biasing, the diffusion due to large scale flows and the mode-coupling due to nonlocal interactions affect the scale dependence of bias from small separations up to the baryon acoustic oscillation (BAO) scale. For 2-sigma density peaks collapsing at z=0.3, our model predicts a ~ 5% residual scale-dependent bias around the acoustic scale that arises mostly from first-order Lagrangian peak biasing (as opposed to second-order gravity mode-coupling). We also search for a scale dependence of bias in the large scale auto-correlation of massive halos extracted from a very large N-body simulation provided by the MICE collaboration. For halos with mass M>10^{14}Msun/h, our measurements demonstrate a scale-dependent bias across the BAO feature which is very well reproduced by a prediction based on the peak model.

preprint2010arXiv

Redshift space correlations and scale-dependent stochastic biasing of density peaks

We calculate the redshift space correlation function and the power spectrum of density peaks of a Gaussian random field. In the linear regime k < 0.1 h/Mpc, the redshift space power spectrum is P^s_{pk}(k,u) = exp(-f^2 s_{vel}^2 k^2 u^2) * [b_{pk}(k) + b_{vel}(k) f u^2]^2 * P_m(k), where u is the angle with respect to the line of sight, s_{vel} is the one-dimensional velocity dispersion, f is the growth rate, and b_{pk}(k) and b_{vel}(k) are k-dependent linear spatial and velocity bias factors. For peaks, the value of s_{vel} depends upon the functional form of b_{vel}. The peaks model is remarkable because it has unbiased velocities -- peak motions are driven by dark matter flows -- but, in order to achieve this, b_{vel} is k-dependent. We speculate that this is true in general: k-dependence of the spatial bias will lead to k-dependence of b_{vel} even if the biased tracers flow with the dark matter. Because of the k-dependence of the linear bias parameters, standard manipulations applied to the peak model will lead to k-dependent estimates of the growth factor that could erroneously be interpreted as a signature of modified dark energy or gravity. We use the Fisher formalism to show that the constraint on the growth rate f is degraded by a factor of two if one allows for a k-dependent velocity bias of the peak type. We discuss a simple estimate of nonlinear evolution and illustrate the effect of the peak bias on the redshift space multipoles. For k < 0.1 h/Mpc, the peak bias is deterministic but k-dependent, so the configuration space bias is stochastic and scale dependent, both in real and redshift space. We provide expressions for this stochasticity and its evolution (abridged).

preprint2010arXiv

The nonlinear redshift space probability distribution function in models with local primordial non-Gaussianity

We use the ellipsoidal collapse approximation to investigate the nonlinear redshift space evolution of the density field with primordial non-Gaussianity of the local f_{nl}-type. We utilize the joint distribution of eigenvalues of the initial non-Gaussian shear field and evaluate the evolved redshift space probability distribution function (PDF). It is shown that, similar to the real space analysis, the underdense tail of the nonlinear redshift space PDF differs significantly from that for Gaussian initial conditions. We also derive the lowest order correction of the Kaiser's formulain the presence of a non-zero f_{nl}.

preprint2010arXiv

The Substructure Hierarchy in Dark Matter Haloes

We present a new algorithm for identifying the substructure within simulated dark matter haloes. The method is an extension of that proposed by Tormen et al. (2004) and Giocoli et al. (2008a), which identifies a subhalo as a group of self-bound particles that prior to being accreted by the main progenitor of the host halo belonged to one and the same progenitor halo (hereafter satellite). However, this definition does not account for the fact that these satellite haloes themselves may also have substructure, which thus gives rise to sub-subhaloes, etc. Our new algorithm identifies substructures at all levels of this hierarchy, and we use it to determine the mass function of all substructure (counting sub-haloes, sub-subhaloes, etc.). On average, haloes which formed more recently tend to have a larger mass fraction in substructure and to be less concentrated than average haloes of the same mass. We provide quantitative fits to these correlations. Even though our algorithm is very different from that of Gao et al. (2004), we too find that the subhalo mass function per unit mass at redshift z = 0 is universal. This universality extends to any redshift only if one accounts for the fact that host haloes of a given mass are less concentrated at higher redshifts, and concentration and substructure abundance are anti-correlated. This universality allows a simple parametrization of the subhalo mass function integrated over all host halo masses, at any given time. We provide analytic fits to this function which should be useful in halo model analyses which equate galaxies with halo substructure when interpreting clustering in large sky surveys. Finally, we discuss systematic differences in the subhalo mass function that arise from different definitions of (host) halo mass.

preprint2009arXiv

Halo abundances in the f_{nl} model

We show how the excursion set moving barrier model for halo abundances may be generalized to the local non-Gaussian f_{nl} model. Our estimate assumes that the distribution of step sizes depends on f_{nl}, but that they are otherwise uncorrelated. Our analysis is consistent with previous results for the case of a constant barrier, and highlights some implicit assumptions. It also clarifies the basis of an approximate analytic solution to the moving barrier problem in the Gaussian case, and shows how it might be improved.

preprint2009arXiv

Reconstructing galaxy fundamental distributions and scaling relations from photometric redshift surveys. Applications to the SDSS early-type sample

Noisy distance estimates associated with photometric rather than spectroscopic redshifts lead to a mis-estimate of the luminosities, and produce a correlated mis-estimate of the sizes. We consider a sample of early-type galaxies from the SDSS DR6 for which both spectroscopic and photometric information is available, and apply the generalization of the V_max method to correct for these biases. We show that our technique recovers the true redshift, magnitude and size distributions, as well as the true size-luminosity relation. We find that using only 10% of the spectroscopic information randomly spaced in our catalog is sufficient for the reconstructions to be accurate within about 3%, when the photometric redshift error is dz = 0.038. We then address the problem of extending our method to deep redshift catalogs, where only photometric information is available. In addition to the specific applications outlined here, our technique impacts a broader range of studies, when at least one distance-dependent quantity is involved. It is particularly relevant for the next generation of surveys, some of which will only have photometric information.

preprint2009arXiv

Sizes and ages of SDSS ellipticals: Comparison with hierarchical galaxy formation models

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that the shape, normalization, and dispersion around the mean size-stellar mass relation is the same for young and old systems, provided the stellar mass is greater than 3*10^10 Msun. This is difficult to reproduce in pure passive evolution models, which generically predict older galaxies to be much more compact than younger ones of the same stellar mass. However, this aspect of our measurements is well reproduced by hierarchical models of galaxy formation. Whereas the models predict more compact galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the more massive objects, resulting in a flat size-age relation at the present time. At lower masses, the models predict that mergers are less frequent, so that the expected anti-correlation between age and size is not completely erased. This is in good agreement with our data: below 3*10^10 Msun, the effective radius R_e is a factor of ~2 lower for older galaxies. These successes of the models are offset by the fact that the predicted sizes have other serious problems, which we discuss.

preprint2009arXiv

The initial shear field in models with primordial local non-Gaussianity and implications for halo and void abundances

We generalize Doroshkevich's celebrated formulae for the eigenvalues of the initial shear field associated with Gaussian statistics to the local non-Gaussian f_{nl} model. This is possible because, to at least second order in f_{nl}, distributions at fixed overdensity are unchanged from the case f_{nl}=0. We use this generalization to estimate the effect of f_{nl}\ne 0 on the abundance of virialized halos. Halo abundances are expected to be related to the probability that a certain quantity in the initial fluctuation field exceeds a threshold value, and we study two choices for this variable: it can either be the sum of the eigenvalues of the initial deformation tensor (the initial overdensity), or its smallest eigenvalue. The approach based on a critical overdensity yields results which are in excellent agreement with numerical measurements. We then use these same methods to develop approximations describing the sensitivity of void abundances on f_{nl}. While a positive f_{nl} produces more extremely massive halos, it makes fewer extremely large voids. Its effect thus is qualitatively different from a simple rescaling of the normalisation of the density fluctuation field σ_8. Therefore, void abundances furnish complementary information to cluster abundances, and a joint comparison of both might provide interesting constraints on primordial non-Gaussianity.

preprint2006arXiv

Bispectrum and Nonlinear Biasing of Galaxies: Perturbation Analysis, Numerical Simulation and SDSS Galaxy Clustering

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

preprint2003arXiv

Composite Luminosity Functions Based on the Sloan Digital Sky Survey "Cut & Enhance" Galaxy Cluster Catalog

We present here results on the composite luminosity functions of galaxies in the clusters of galaxies selected from the SDSS Cut and Enhance cluster catalog (CE; Goto et al. 2001). We construct the composite luminosity function in the five SDSS bands, u,g,r,i and z, using 204 CE clusters ranging from z=0.02 to z=0.25. We use photometric redshifts to construct composite luminosity functions. The robustness of the weighting scheme was tested using Monte Carlo simulation. We find the slope of composite LFs become flatter toward redder color band. Comparing with field LFs of SDSS (Blanton et al. 2001), cluster LFs have brighter characteristic magnitude and flatter slopes in g, r, i and z bands. These results are consistent with the hypothesis that the cluster LF has two distinct underlying populations i.e. the bright end of the LF is dominated by bright early types that follow a gaussian--like luminosity distribution, while the faint--end of the cluster LF is a steep power--law like function dominated by star--forming (bluer) galaxies. We also study the composite luminosity functions for early--type and late--type galaxies using profile fits, a concentration parameter and u-r color to classify galaxy morphology. The strong dependence of LF on galaxy morphology is found.The faint end slope of the LF is always flatter for early--type galaxies than late--type regardless of passband and methodology. These results are consistent with the hypothesis that the cluster region are dominated by bright elliptical galaxies.

preprint2003arXiv

The Environment of Passive Spiral Galaxies in the SDSS

In previous work on galaxy clusters, several authors reported a discovery of an unusual population of galaxies, which have spiral morphologies, but do not show any star formation activity. These galaxies are called ``passive spirals'', and have been interesting since it has been difficult to understand the existence of such galaxies. Using a volume limited sample (0.05<z<0.1 and Mr<-20.5; 25813 galaxies) of the Sloan Digital Sky Survey data, we have found 73 (0.28$\pm$0.03%) passive spiral galaxies and studied their environments. It is found that passive spiral galaxies live in local galaxy density 1-2 Mpc$^{-2}$ and 1-10 cluster-centric virial radius. Thus the origins of passive spiral galaxies are likely to be cluster related. These characteristic environments coincide with the previously reported environment where galaxy star formation rate suddenly declines and the so-called morphology-density relation turns. It is likely that the same physical mechanism is responsible for all of these observational results. The existence of passive spiral galaxies suggests that a physical mechanism that works calmly is preferred to dynamical origins such as major merger/interaction since such a mechanism can destroy spiral arm structures. Compared with observed cluster galaxy evolution such as the Butcher-Oemler effect and the morphological Butcher-Oemler effect, passive spiral galaxies are likely to be a galaxy population in transition between red, elliptical/S0 galaxies in low redshift clusters and blue, spiral galaxies more numerous in higher redshift clusters.

preprint2003arXiv

The Morphological Butcher-Oemler effect in the SDSS Cut&Enhance Galaxy Cluster Catalog

We investigate the evolution of the fractions of late type cluster galaxies as a function of redshift, using one of the largest, most uniform cluster samples available. The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created using a single automated cluster finding algorithm on uniform data from a single telescope, with accurate CCD photometry, thus, minimizing selection biases. We use four independent methods to analyze the evolution of the late type galaxy fraction. Specifically, we select late type galaxies based on: restframe g-r color, u-r color, galaxy profile fitting and concentration index. The first criterion corresponds to the one used in the classical Butcher-Oemler analyses. The last three criteria are more sensitive to the morphological type of the galaxies. In all four cases, we find an increase in the fraction of late type galaxies with increasing redshift, significant at the 99.9% level. The results confirm that cluster galaxies do change colors with redshift (the Butcher-Oemler effect) and, in addition, they change their morphology to later-type toward higher redshift -- indicating a morphological equivalent of the Butcher-Oemler effect. We also find a tendency of richer clusters to have lower fractions of late type galaxies. The trend is consistent with a ram pressure stripping model, where richer clusters have more effective ram pressure due to their higher temperature.

preprint1996arXiv

Galton--Watson branching processes and the growth of gravitational clustering

The Press--Schechter description of gravitational clustering from an initially Poisson distribution is shown to be equivalent to the well studied Galton--Watson branching process. This correspondence is used to provide a detailed description of the evolution of hierarchical clustering, including a complete description of the merger history tree. The relation to branching process epidemic models means that the Press--Schechter description can also be understood using the formalism developed in the study of queues. The queueing theory formalism, also, is used to provide a complete description of the merger history of any given Press--Schechter clump. In particular, an analytic expression for the merger history of any given Poisson Press--Schechter clump is obtained. This expression allows one to calculate the partition function of merger history trees. It obeys an interesting scaling relation; the partition function for a given pair of initial and final epochs is the same as that for certain other pairs of initial and final epochs. The distribution function of counts in randomly placed cells, as a function of time, is also obtained using the branching process and queueing theory descriptions. Thus, the Press--Schechter description of the gravitational evolution of clustering from an initially Poisson distribution is now complete. All these interrelations show why the Press--Schechter approach works well in a statistical sense, but cannot provide a detailed description of the dynamics of the clustering particles themselves. One way to extend these results to more general Gaussian initial conditions is discussed.

preprint1996arXiv

The non-linear correlation function and the shapes of virialized halos

The correlation function xi(r) of matter in the non-linear regime is assumed to be determined by the density profiles rho(r) and the mass distribution n(M) of virialized halos. The Press--Schechter approach is used to compute n(M), and the stable clustering hypothesis is used to determine the density profiles of these Press--Schechter halos. Thus, the shape and amplitude of xi(r) on small scales is related to the initial power spectrum of density fluctuations. The case of clustering from scale-free initial conditions is treated in detail. If n is the slope of the initial power spectrum of density fluctuations, then stable clustering requires that xi(r)\propto r^{-gamma}, where gamma is a known function of n. If halo--halo correlations can be neglected, then rho(r)\propto r^{-epsilon}, where epsilon = (gamma+3)/2 = 3(4+n)/(5+n). For all values of n of current interest, this slope is steeper than the value 3(3+n)/(4+n) that was obtained by Hoffman & Shaham in their treatment of the shapes of the outer regions of collapsed halos. Our main result is a prediction for the amplitude of the non-linear correlation function. The predicted amplitude and its dependence on n are in good quantitative agreement with N-body simulations of self-similar clustering. If stable clustering is a good approximation only inside the half-mass radii of Press--Schechter halos, then the density contrast required for the onset of stable clustering can be estimated. This density contrast is in the range ~300-600 and increases with the initial slope n, in agreement with estimates from N-body simulations.

preprint1995arXiv

Constrained realizations and minimum variance reconstruction of non-Gaussian random fields

With appropriate modifications, the Hoffman--Ribak algorithm that constructs constrained realizations of Gaussian random fields having the correct ensemble properties can also be used to construct constrained realizations of those non-Gaussian random fields that are obtained by transformations of an underlying Gaussian field. For example, constrained realizations of lognormal, generalized Rayleigh, and chi-squared fields having $n$ degrees of freedom constructed this way will have the correct ensemble properties. The lognormal field is considered in detail. For reconstructing Gaussian random fields, constrained realization techniques are similar to reconstructions obtained using minimum variance techniques. A comparison of this constrained realization approach with minimum variance, Wiener filter reconstruction techniques, in the context of lognormal random fields, is also included. The resulting prescriptions for constructing constrained realizations as well as minimum variance reconstructions of lognormal random fields are useful for reconstructing masked regions in galaxy catalogues on smaller scales than previously possible, for assessing the statistical significance of small-scale features in the microwave background radiation, and for generating certain non-Gaussian initial conditions for $N$-body simulations.

preprint1995arXiv

The distribution of pairwise peculiar velocities in the nonlinear regime

The distribution of pairwise, relative peculiar velocities, $f(u;r)$, on small nonlinear scales, $r$, is derived from the Press--Schechter approach. This derivation assumes that Press--Schechter clumps are virialized and isothermal. The virialized assumption requires that the circular velocity, $V_c \propto M^{1/3}$, where $M$ denotes the mass of the clump. The isothermal assumption means that the circular velocity is independent of radius. Further, it is assumed that the velocity distribution within a clump is Maxwellian, that the pairwise relative velocity distribution is isotropic, and that on nonlinear scales clump-clump motions are unimportant when calculating the distribution of velocity differences. Comparison with $N$-body simulations shows that, on small nonlinear scales, all these assumptions are accurate. For most power spectra of interest, the resulting line of sight, pairwise, relative velocity distribution, $f(u_{\rm r})$, is well approximated by an exponential, rather than a Gaussian distribution. This simple Press--Schechter model is also able to provide a natural explanation for the observed, non-Maxwellian shape of $f(v)$, the distribution of peculiar velocities.

Ravi K. Sheth

What is connected

Connect this record

See the researcher in context

Building this map preview

56 published item(s)

Protohalos and their connection to halo assembly, shape and structure

Optimal Transport Reconstruction of Baryon Acoustic Oscillations

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

The smearing scale in Laguerre reconstructions of the correlation function

A universal threshold for primordial black hole formation

Analytical thresholds for black hole formation in general cosmological backgrounds

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

Linear Point and Sound Horizon as Purely Geometric standard rulers

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

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

Excursion set peaks: the role of shear

The universality of the virial halo mass function and models for non-universality of other halo definitions

On the Markovian assumption in the excursion set approach: the approximation of Markov Velocities

Stochasticity in halo formation and the excursion set approach

The excursion set approach in non-Gaussian random fields

Bias deconstructed: Unravelling the scale dependence of halo bias using real space measurements

Ellipsoidal halo finders and implications for models of triaxial halo formation

Excursion set peaks: a self-consistent model of dark halo abundances and clustering

Lagrangian bias in the local bias model

Measures of galaxy environment -- II. Rank-ordered mark correlations

Non-local Lagrangian bias

Peak-peak correlations in the cosmic background radiation from cosmic strings

Stochastic bias in multi-dimensional excursion set approaches

The importance of stepping up in the excursion set approach

Formation times, mass growth histories and concentrations of dark matter haloes

Gravity and Large-Scale Non-local Bias

On the anisotropic density distribution on large scales

One step beyond: The excursion set approach with correlated steps

Peaks theory and the excursion set approach

Plain fundamentals of Fundamental Planes: Analytics and algorithms

Scale dependent halo bias in the excursion set approach

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: a large sample of mock galaxy catalogues

The luminosities of the brightest cluster galaxies and brightest satellites in SDSS groups

A hierarchy of voids: More ado about nothing

Halo abundances and counts-in-cells: The excursion set approach with correlated steps

Halo bias in the excursion set approach with correlated steps

Modelling the shapes of the largest gravitationally bound objects

Scale-dependent halo bias from scale-dependent growth

Convolution and deconvolution based estimates of galaxy scaling relations from photometric redshift surveys

Halo model description of the non-linear dark matter power spectrum at $k \gg 1$ Mpc$^{-1}$

Modeling scale-dependent bias on the baryonic acoustic scale with the statistics of peaks of Gaussian random fields

Redshift space correlations and scale-dependent stochastic biasing of density peaks

The nonlinear redshift space probability distribution function in models with local primordial non-Gaussianity

The Substructure Hierarchy in Dark Matter Haloes

Halo abundances in the f_{nl} model

Reconstructing galaxy fundamental distributions and scaling relations from photometric redshift surveys. Applications to the SDSS early-type sample

Sizes and ages of SDSS ellipticals: Comparison with hierarchical galaxy formation models

The initial shear field in models with primordial local non-Gaussianity and implications for halo and void abundances

Bispectrum and Nonlinear Biasing of Galaxies: Perturbation Analysis, Numerical Simulation and SDSS Galaxy Clustering

Composite Luminosity Functions Based on the Sloan Digital Sky Survey "Cut & Enhance" Galaxy Cluster Catalog

The Environment of Passive Spiral Galaxies in the SDSS

The Morphological Butcher-Oemler effect in the SDSS Cut&Enhance Galaxy Cluster Catalog

Galton--Watson branching processes and the growth of gravitational clustering

The non-linear correlation function and the shapes of virialized halos

Constrained realizations and minimum variance reconstruction of non-Gaussian random fields

The distribution of pairwise peculiar velocities in the nonlinear regime