BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Aseem Paranjape

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMhep-phhep-th

Catalog footprint

What is connected

34works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

34 published item(s)

preprint2026arXiv

Sahyadri: A simulation suite for the cosmology dependence of the Cosmic Web

We present Sahyadri, a suite of cosmological $N$-body simulations designed to enable precision studies of the low-redshift Universe with next-generation spectroscopic surveys. Sahyadri includes systematic variations of six cosmological parameters around Planck 2018 constraints, with seed-matched initial conditions enabling cosmological parameter derivatives. Each simulation evolves $2048^3$ particles in a periodic box of side length $200$ $h^{-1}$ Mpc, yielding a particle mass of $m_{\rm{p}} = 8.1 \times 10^{7}\,h^{-1}\,M_{\odot}$ in the fiducial Planck 2018 cosmology. This resolution enables robust identification of dark matter halos down to $M_{\rm min} = 3.2 \times 10^{9}$ $h^{-1}$ $M_\odot$, which represents a factor of $\sim$25 improvement over the AbacusSummit suite, and is over two orders of magnitude better than the Quijote and Aemulus suites. We estimate that approximately 40% of DESI BGS galaxies at redshift $z < 0.15$ - roughly 1.6 million objects - reside in halos accessible to Sahyadri but beyond the reach of existing parameter-varying simulation suites. We demonstrate Sahyadri's capabilities through measurements of the matter power spectrum, halo mass function and power spectrum, and beyond 2-point statistics such as the Voronoi volume function and $k^{\rm th}$ nearest neighbour statistics, showing excellent agreement with theoretical predictions and significant sensitivity to $Ω_{\rm m}$ variations. We implement a custom compression scheme reducing storage requirements by a factor of $\sim$3 while maintaining sub-percent clustering accuracy. Key data products will be made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A simulated annealing approach to parameter inference with expensive likelihoods

We present a new approach to parameter inference targeted on generic situations where the evaluation of the likelihood $\mathcal{L}$ (i.e., the probability to observe the data given a fixed model configuration) is numerically expensive. Inspired by ideas underlying simulated annealing, the method first evaluates $χ^2=-2\ln\mathcal{L}$ on a sparse sequence of Latin hypercubes of increasing density in parameter (eigen)space. The semi-stochastic choice of sampling points accounts for anisotropic gradients of $χ^2$ and rapidly zooms in on the minimum of $χ^2$. The sampled $χ^2$ values are then used to train an interpolator which is further used in a standard Markov Chain Monte Carlo (MCMC) algorithm to inexpensively explore the parameter space with high density, similarly to emulator-based approaches now popular in cosmological studies. Comparisons with example linear and non-linear problems show gains in the number of likelihood evaluations of factors of 10 to 100 or more, as compared to standard MCMC algorithms. As a specific implementation, we publicly release the code PICASA: Parameter Inference using Cobaya with Anisotropic Simulated Annealing, which combines the minimizer (of a user-defined $χ^2$) with Gaussian Process Regression for training the interpolator and a subsequent MCMC implementation using the COBAYA framework. Being agnostic to the nature of the observable data and the theoretical model, our implementation is potentially useful for a number of emerging problems in cosmology, astrophysics and beyond.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mass accretion rates and multi-scale halo environment in cold and warm dark matter cosmologies

We study the evolving environment dependence of mass accretion by dark haloes in simulations of cold and warm dark matter (CDM and WDM) cosmologies. The latter allows us to probe the nature of halo growth at scales below the WDM half-mode mass, which form an extreme regime of nonlinear collisionless dynamics and offer an excellent test-bed for ideas relating to hierarchical growth. As environmental proxies, we use the local halo-centric matter density $δ$ and tidal anisotropy $α$, as well as large-scale halo bias $b_1$. Our analysis, while reproducing known trends for environment-dependent accretion in CDM, as well as the comparison between accretion in CDM and WDM, reveals several interesting new features. As expected from excursion set models, WDM haloes have higher specific accretion rates, dominated by the accretion of diffuse mass, as compared to CDM haloes. For low-mass WDM haloes, we find that the environment-dependence of both diffuse mass accretion as well as accretion by mergers is almost fully explained by $α$. For the other cases, $δ$ plays at least a comparable role. We detect, for the first time, a significant and evolving assembly bias due to diffuse mass accretion for low-mass CDM and WDM haloes (after excluding splashback objects), with a $z=0$ strength higher than with almost all known secondary variables and largely explained by $α$. Our results place constraints on semi-analytical merger tree algorithms, which in turn could affect the predictions of galaxy evolution models based on them.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Analytical halo models of cosmic tidal fields

The non-linear cosmic web environment of dark matter haloes plays a major role in shaping their growth and evolution, and potentially also affects the galaxies that reside in them. We develop an analytical (halo model) formalism to describe the tidal field of anisotropic halo-centric density distributions, as characterised by the halo-centric tidal tensor $\langle T_{ij} \rangle(<R)$ spherically averaged on scale $R\sim4R_{\rm vir}$ for haloes of virial radius $R_{\rm vir}$. We focus on axisymmetric anisotropies, which allows us to explore simple and intuitive toy models of (sub)halo configurations that exemplify some of the most interesting anisotropies in the cosmic web. We build our models around the spherical Navarro-Frenk-White (NFW) profile after describing it as a Gaussian mixture, which leads to almost fully analytical expressions for the `tidal anisotropy' scalar $α(<4R_{\rm vir})$ extracted from the tidal tensor. Our axisymmetric examples include (i) a spherical halo at the axis of a cylindrical filament, (ii) an off-centred satellite in a spherical host halo and (iii) an axisymmetric halo. Using these, we demonstrate several interesting results. For example, the tidal tensor at the axis of a pure cylindrical filament gives $α^{\rm (fil)}(<R)=1/2$ exactly, for any $R$. Also, $α(<4R_{\rm vir,sat})$ for a satellite of radius $R_{\rm vir,sat}$ as a function of its host-centric distance is a sensitive probe of dynamical mass loss of the satellite in its host environment. Finally, we discuss a number of potentially interesting extensions and applications of our formalism that can deepen our understanding of the multi-scale phenomenology of the cosmic web.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A sea of tides: star formation and the central-satellite dichotomy in a continuum of tidal environments

The environment-dependent bimodality of the distribution of stellar mass ($M_\ast$) and specific star formation rate (sSFR) of galaxies, and its explanation in terms of the central-satellite dichotomy, form a cornerstone of our current understanding of galaxy evolution in the hierarchical structure formation paradigm. We revisit this framework in the IllustrisTNG simulation in the context of the most extreme local tidal anisotropy $α_{\rm peak}$ experienced by each galaxy over cosmic time, which is an excellent proxy for environmental influence. We show that, while sharing a common monotonic $M_\ast$-$v_{\rm peak}$ relation, central, satellite and `splashback' galaxies define a hierarchy of increasing $α_{\rm peak}$. We also find that the sSFR of objects in small haloes unaffected by feedback from an active nucleus typically decreases with increasing $α_{\rm peak}$. Our results support an alternate viewpoint in which a galaxy can be identified by the value of $α_{\rm peak}$; i.e., rather than being placed on the central-satellite dichotomy, a galaxy is better classified by its location in a continuum of tidal environments. This conceptual shift can potentially yield a more robust understanding of galaxy evolution and the galaxy-dark matter connection, e.g., in accurately modelling subtle effects such as sSFR-induced secondary clustering.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Studying the Lyman-$α$ optical depth fluctuations at $z \sim 5.5$ using fast semi-numerical methods

We present a computationally efficient and fast semi-numerical technique for simulating the Lyman-$α$ (Ly$α$) absorption optical depth in presence of neutral hydrogen "islands" left over from reionization at redshifts $5 \lesssim z \lesssim 6$. The main inputs to the analysis are (i) a semi-numerical photon-conserving model of ionized regions during reionization (named SCRIPT) along with a prescription for simulating the shadowing by neutral islands and (ii) the fluctuating Gunn-Peterson approximation to model the Ly$α$ absorption. Our model is then used for simulating the large-scale fluctuations in the effective optical depth as observed along sight lines towards high$-z$ quasars. Our model is fully described by five parameters. By setting two of them to default values and varying the other three, we obtain the constraints on reionization history at $5 \lesssim z \lesssim 6$ as allowed by the data. We confirm that reionization is \emph{not} complete before $z \sim 5.6$ at $\gtrsim 2σ$ confidence, with the exact confidence limits depending on how the non-detections of the flux in the data are treated. We also confirm that the completion of reionization can be as late as $z \sim 5.2$. With further improvements in the model and with more sight lines at $z \sim 6$, we can take advantage of the computational efficiency of our analysis to obtain more stringent constraints on the ionization fraction at the tail-end of reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Voronoi volume function: A new probe of cosmology and galaxy evolution

We study the Voronoi volume function (VVF) -- the distribution of cell volumes (or inverse local number density) in the Voronoi tessellation of any set of cosmological tracers (galaxies/haloes). We show that the shape of the VVF of biased tracers responds sensitively to physical properties such as halo mass, large-scale environment, substructure and redshift-space effects, making this a hitherto unexplored probe of both primordial cosmology and galaxy evolution. Using convenient summary statistics -- the width, median and a low percentile of the VVF as functions of average tracer number density -- we explore these effects for tracer populations in a suite of N-body simulations of a range of dark matter models. Our summary statistics sensitively probe primordial features such as small-scale oscillations in the initial matter power spectrum (as arise in models involving collisional effects in the dark sector), while being largely insensitive to a truncation of initial power (as in warm dark matter models). For vanilla cold dark matter (CDM) cosmologies, the summary statistics display strong evolution and redshift-space effects, and are also sensitive to cosmological parameter values for realistic tracer samples. Comparing the VVF of galaxies in the GAMA survey with that of abundance matched CDM (sub)haloes tentatively reveals environmental effects in GAMA beyond halo mass (modulo unmodelled satellite properties). Our exploratory analysis thus paves the way for using the VVF as a new probe of galaxy evolution physics as well as the nature of dark matter and dark energy.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

General relativistic 'screening' in cosmological simulations

We revisit the issue of interpreting the results of large volume cosmological simulations in the context of large scale general relativistic effects. We look for simple modifications to the nonlinear evolution of the gravitational potential $ψ$ that lead on large scales to the correct, fully relativistic description of density perturbations in the Newtonian gauge. We note that the relativistic constraint equation for $ψ$ can be cast as a diffusion equation, with a diffusion length scale determined by the expansion of the Universe. Exploiting the weak time evolution of $ψ$ in all regimes of interest, this equation can be further accurately approximated as a Helmholtz equation, with an effective relativistic 'screening' scale $\ell$ related to the Hubble radius. We demonstrate that it is thus possible to carry out N-body simulations in the Newtonian gauge by replacing Poisson's equation with this Helmholtz equation, involving a trivial change in the Green's function kernel. Our results also motivate a simple, approximate (but very accurate) gauge transformation - $δ_{\rm N}(\mathbf{k}) \approx δ_{\rm sim}(\mathbf{k})\times (k^2+\ell^{-2})/k^2$ - to convert the density field $δ_{\rm sim}$ of standard collisionless N-body simulations (initialised in the comoving synchronous gauge) into the Newtonian gauge density $δ_{\rm N}$ at arbitrary times. A similar conversion can also be written in terms of particle positions. Our results can be interpreted in terms of a Jeans stability criterion induced by the expansion of the Universe. The appearance of the screening scale $\ell$ in the evolution of $ψ$, in particular, leads to a natural resolution of the 'Jeans swindle' in the presence of super-horizon modes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modelling the 21 cm Signal From the Epoch of Reionization and Cosmic Dawn

Studying the cosmic dawn and the epoch of reionization through the redshifted 21 cm line are among the major science goals of the SKA1. Their significance lies in the fact that they are closely related to the very first stars in the universe. Interpreting the upcoming data would require detailed modelling of the relevant physical processes. In this article, we focus on the theoretical models of reionization that have been worked out by various groups working in India with the upcoming SKA in mind. These models include purely analytical and semi-numerical calculations as well as fully numerical radiative transfer simulations. The predictions of the 21 cm signal from these models would be useful in constraining the properties of the early galaxies using the SKA data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Photon Number Conserving Models of H II Bubbles during Reionization

Traditional excursion set based models of H II bubble growth during the epoch of reionization are known to violate photon number conservation, in the sense that the mass fraction in ionized bubbles in these models does not equal the ratio of the number of ionizing photons produced by sources and the number of hydrogen atoms in the intergalactic medium. E.g., for a Planck13 cosmology with electron scattering optical depth $τ\simeq0.066$, the discrepancy is $\sim15$ per cent for $x_{\rm HII}=0.1$ and $\sim5$ per cent for $x_{\rm HII}=0.5$. We demonstrate that this problem arises from a fundamental conceptual shortcoming of the excursion set approach (already recognised in the literature on this formalism) which only tracks average mass fractions instead of the exact, stochastic source counts. With this insight, we build an approximately photon number conserving Monte Carlo model of bubble growth based on partitioning regions of dark matter into halos. Our model, which is formally valid for white noise initial conditions (ICs), shows dramatic improvements in photon number conservation, as well as substantial differences in the bubble size distribution, as compared to traditional models. We explore the trends obtained on applying our algorithm to more realistic ICs, finding that these improvements are robust to changes in the ICs. Since currently popular semi-numerical schemes of bubble growth also violate photon number conservation, we argue that it will be worthwhile to pursue new, explicitly photon number conserving approaches. Along the way, we clarify some misconceptions regarding this problem that have appeared in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Simulating the Large-Scale Structure of HI Intensity Maps

Intensity mapping of neutral hydrogen (HI) is a promising observational probe of cosmology and large-scale structure. We present wide field simulations of HI intensity maps based on N-body simulations of a $2.6\, {\rm Gpc / h}$ box with $2048^3$ particles (particle mass $1.6 \times 10^{11}\, {\rm M_\odot / h}$). Using a conditional mass function to populate the simulated dark matter density field with halos below the mass resolution of the simulation ($10^{8}\, {\rm M_\odot / h} < M_{\rm halo} < 10^{13}\, {\rm M_\odot / h}$), we assign HI to those halos according to a phenomenological halo to HI mass relation. The simulations span a redshift range of 0.35 < z < 0.9 in redshift bins of width $Δz \approx 0.05$ and cover a quarter of the sky at an angular resolution of about 7'. We use the simulated intensity maps to study the impact of non-linear effects and redshift space distortions on the angular clustering of HI. Focusing on the autocorrelations of the maps, we apply and compare several estimators for the angular power spectrum and its covariance. We verify that these estimators agree with analytic predictions on large scales and study the validity of approximations based on Gaussian random fields, particularly in the context of the covariance. We discuss how our results and the simulated maps can be useful for planning and interpreting future HI intensity mapping surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Correlating galaxy colour and halo concentration: A tunable Halo Model of galactic conformity

We extend the Halo Occupation Distribution (HOD) framework to generate mock galaxy catalogs exhibiting varying levels of "galactic conformity", which has emerged as a potentially powerful probe of environmental effects in galaxy evolution. Our model correlates galaxy colours in a group with the concentration of the common parent dark halo through a "group quenching efficiency" $ρ$ which makes older, more concentrated halos $\textit{at fixed mass}$ preferentially host redder galaxies. We find that, for a specific value of $ρ$, this 1-halo conformity matches corresponding measurements in a group catalog based on the Sloan Digital Sky Survey. Our mocks also display conformity at large separations from isolated objects, potentially an imprint of halo assembly bias. A detailed study - using mocks with assembly bias erased while keeping 1-halo conformity intact - reveals a rather nuanced situation, however. At separations $\lesssim 4$Mpc, conformity is mainly a 1-halo effect dominated by the largest halos and is $\textit{not}$ a robust indicator of assembly bias. Only at very large separations ($\gtrsim 8$Mpc) does genuine 2-halo conformity, driven by the assembly bias of small halos, manifest distinctly. We explain all these trends in standard Halo Model terms. Our model opens the door to parametrized HOD analyses that self-consistently account for galactic conformity at all scales.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Information Gains from Cosmic Microwave Background Experiments

To shed light on the fundamental problems posed by Dark Energy and Dark Matter, a large number of experiments have been performed and combined to constrain cosmological models. We propose a novel way of quantifying the information gained by updates on the parameter constraints from a series of experiments which can either complement earlier measurements or replace them. For this purpose, we use the Kullback-Leibler divergence or relative entropy from information theory to measure differences in the posterior distributions in model parameter space from a pair of experiments. We apply this formalism to a historical series of Cosmic Microwave Background experiments ranging from Boomerang to WMAP, SPT, and Planck. Considering different combinations of these experiments, we thus estimate the information gain in units of bits and distinguish contributions from the reduction of statistical errors and the `surprise' corresponding to a significant shift of the parameters' central values. For this experiment series, we find individual relative entropy gains ranging from about 1 to 30 bits. In some cases, e.g. when comparing WMAP and Planck results, we find that the gains are dominated by the surprise rather than by improvements in statistical precision. We discuss how this technique provides a useful tool for both quantifying the constraining power of data from cosmological probes and detecting the tensions between experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Testing spherical evolution for modelling void abundances

We compare analytical predictions of void volume functions to those measured from N-body simulations, detecting voids with the zobov void finder. We push to very small, nonlinear voids, below few Mpc radius, by considering the unsampled DM density field. We also study the case where voids are identified using halos. We develop analytical formula for the void abundance of both the excursion set approach and the peaks formalism. These formula are valid for random walks smoothed with a top-hat filter in real space, with a large class of realistic barrier models. We test the extent to which the spherical evolution approximation, which forms the basis of the analytical predictions, models the highly aspherical voids that occur in the cosmic web, and are found by a watershed-based algorithm such as zobov. We show that the volume function returned by zobov is quite sensitive to the choice of treatment of sub-voids, a fact that has not been appreciated previously. For reasonable choices of sub-void exclusion, we find that the Lagrangian density delta_v of the zobov voids -- which is predicted to be a constant delta_v = -2.7 in the spherical evolution model -- is different from the predicted value, showing substantial scatter and scale dependence. This result applies to voids identified at z=0 with effective radius between 1 and 10 Mpc/h. Our analytical approximations are flexible enough to give a good description of the resulting volume function; however, this happens for choices of parameter values that are different from those suggested by the spherical evolution assumption. We conclude that analytical models for voids must move away from the spherical approximation in order to be applied successfully to observations, and we discuss some possible ways forward.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

3-dimensional spherical analyses of cosmological spectroscopic surveys

Spectroscopic redshift surveys offer great prospects for constraining the dark sector in cosmology. Future surveys will however be both deep and wide and will thus require an analysis in 3-dimensional spherical geometry. We review and compare several methods which have been proposed in the literature for this purpose, focusing in particular on implementations of the spherical harmonic tomography (SHT) power spectrum $C^{i j}_{l}$ and the spherical Fourier Bessel (SFB) power spectrum $C_{l} (k, k')$. Using a Fisher analysis, we compare the forecasted constraints on cosmological parameters using these statistics. These constraints typically rely on approximations such as the Limber approximation and make specific choices in the numerical implementation of each statistic. Using a series of toy models, we explore the applicability of these approximations and study the sensitivity of the SHT and SFB statistics to the details of their implementation. In particular, we show that overlapping redshift bins may improve cosmological constraints using the SHT statistic when the number of bins is small, and that the SFB constraints are quite robust to changes in the assumed distance-redshift relation. We also find that the SHT can be tailored to be more sensitive to modes at redshifts close to the survey boundary, while the SFB appears better suited to capture information beyond the smooth shape of the power spectrum. In this context, we discuss the pros and cons of the different techniques and their impact on the design and analysis of future wide field spectroscopic surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A simple model linking galaxy and dark matter evolution

We construct a simple phenomenological model for the evolving galaxy population by incorporating pre-defined baryonic prescriptions into a dark matter hierarchical merger tree. Specifically the model is based on the simple gas-regulator model introduced by Lilly et al. 2013 coupled with the empirical quenching rules of Peng et al. 2010/12. The simplest model already does quite well in reproducing, without re-adjusting the input parameters, many observables including the Main Sequence sSFR-mass relation, the faint end slope of the galaxy mass function and the shape of the star-forming and passive mass functions. Compared with observations and/or the recent phenomenological model of Behroozi et al. 2013 based on epoch-dependent abundance-matching, our model also qualitatively reproduces the evolution of the Main Sequence sSFR(z) and SFRD(z) star formation rate density relations, the $M_s - M_h$ stellar-to-halo mass relation and also the $SFR - M_h$ relation. Quantitatively the evolution of sSFR(z) and SFRD(z) is not steep enough, the $M_s - M_h$ relation is not quite peaked enough and, surprisingly, the ratio of quenched to star-forming galaxies around M* is not quite high enough. We show that these deficiencies can simultaneously be solved by ad hoc allowing galaxies to re-ingest some of the gas previously expelled in winds, provided that this is done in a mass-dependent and epoch-dependent way. These allow the model galaxies to reduce an inherent tendency to saturate their star-formation efficiency. This emphasizes how efficient galaxies around M* are in converting baryons into stars and highlights the fact that quenching occurs just at the point when galaxies are rapidly approaching the maximum possible efficiency of converting baryons into stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An improved model of HII bubbles during the epoch of reionization

The size distribution of ionized regions during the epoch of reionization -- a key ingredient in understanding the HI power spectrum observable by 21cm experiments -- can be modelled analytically using the excursion set formalism of random walks in the smoothed initial density field. To date, such calculations have been based on simplifying assumptions carried forward from the earliest excursion set models of two decades ago. In particular, these models assume that the random walks have uncorrelated steps and that haloes can form at arbitrary locations in the initial density field. We extend these calculations by incorporating recent technical developments that allow us to (a) include the effect of correlations in the steps of the walks induced by a realistic smoothing filter and (b) more importantly, account for the fact that dark matter haloes preferentially form near peaks in the initial density. A comparison with previous calculations shows that including these features, particularly the peaks constraint on halo locations, has large effects on the size distribution of the HII bubbles surrounding these haloes. For example, when comparing models at the same value of the globally averaged ionized volume fraction, the typical bubble sizes predicted by our model are more than a factor 2 larger than earlier calculations. Our results can potentially have a significant impact on estimates of the observable HI power spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmology with Galaxy Clusters: Systematic Effects in the Halo Mass Function

We investigate potential systematic effects in constraining the amplitude of primordial fluctuations σ_8 arising from the choice of halo mass function in the likelihood analysis of current and upcoming galaxy cluster surveys. We study the widely used N-body simulation fit of Tinker et al. (T08) and, as an alternative, the recently proposed analytical model of Excursion Set Peaks (ESP). We first assess the relative bias between these prescriptions when constraining σ_8 by sampling the ESP mass function to generate mock catalogs and using the T08 fit to analyse them, for various choices of survey selection threshold, mass definition and statistical priors. To assess the level of absolute bias in each prescription, we then repeat the analysis on dark matter halo catalogs in N-body simulations designed to mimic the mass distribution in the current data release of Planck SZ clusters. This N-body analysis shows that using the T08 fit without accounting for the scatter introduced when converting between mass definitions (alternatively, the scatter induced by errors on the parameters of the fit) can systematically over-estimate the value of σ_8 by as much as 2σ for current data, while analyses that account for this scatter should be close to unbiased in σ_8. With an increased number of objects as expected in upcoming data releases, regardless of accounting for scatter, the T08 fit could over-estimate the value of σ_8 by ~1.5σ. The ESP mass function leads to systematically more biased but comparable results. A strength of the ESP model is its natural prediction of a weak non-universality in the mass function which closely tracks the one measured in simulations and described by the T08 fit. We suggest that it might now be prudent to build new unbiased ESP-based fitting functions for use with the larger datasets of the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Measuring nonlocal Lagrangian peak bias

We investigate nonlocal Lagrangian bias contributions involving gradients of the linear density field, for which we have predictions from the excursion set peak formalism. We begin by writing down a bias expansion which includes all the bias terms, including the nonlocal ones. Having checked that the model furnishes a reasonable fit to the halo mass function, we develop a 1-point cross-correlation technique to measure bias factors associated with 2-distributed quantities. We validate the method with numerical realizations of peaks of Gaussian random fields before we apply it to N-body simulations. We focus on the lowest (quadratic) order nonlocal contributions. We can reproduce our measurement of χ_{10} if we allow for an offset between the Lagrangian halo center-of-mass and the peak position. The sign and magnitude of χ_{10} is consistent with Lagrangian haloes sitting near linear density maxima. The resulting contribution to the halo bias can safely be ignored for M = 10^13 Msun/h, but could become relevant at larger halo masses. For the second nonlocal bias χ_{01} however, we measure a much larger magnitude than predicted by our model. We speculate that some of this discrepancy might originate from nonlocal Lagrangian contributions induced by nonspherical collapse.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The locations of halo formation and the peaks formalism

We investigate the problem of predicting the halo mass function from the properties of the Lagrangian density field. We focus on a perturbation spectrum with a small-scale cut-off (as in warm dark matter cosmologies). This cut-off results in a strong suppression of low mass objects, providing additional leverage to rigorously test which perturbations collapse and to what mass. We find that all haloes are consistent with forming near peaks of the initial density field, with a strong correlation between proto-halo density and ellipticity. We demonstrate that, while standard excursion set theory with correlated steps completely fails to reproduce the mass function, the inclusion of the peaks constraint leads to the correct number of haloes but significantly underpredicts the masses of low-mass objects (with the predicted halo mass function at low masses behaving like dn/dln m ~ m^{2/3}). This prediction is very robust and cannot be easily altered within the framework of a single collapse barrier. The nature of collapse in the presence of a small-scale cut-off thus reveals that excursion set calculations require a more detailed understanding of the collapse-time of a general ellipsoidal perturbation to predict the ultimate collapsed mass of a peak -- a problem that has been hidden in the large abundance of small-scale structure in CDM. We demonstrate how this problem can be resolved within the excursion set framework.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Improved Calculation of the Non-Gaussian Halo Mass Function

The abundance of collapsed objects in the universe, or halo mass function, is an important theoretical tool in studying the effects of primordially generated non-Gaussianities on the large scale structure. The non-Gaussian mass function has been calculated by several authors in different ways, typically by exploiting the smallness of certain parameters which naturally appear in the calculation, to set up a perturbative expansion. We improve upon the existing results for the mass function by combining path integral methods and saddle point techniques (which have been separately applied in previous approaches). Additionally, we carefully account for the various scale dependent combinations of small parameters which appear. Some of these combinations in fact become of order unity for large mass scales and at high redshifts, and must therefore be treated non-perturbatively. Our approach allows us to do this, and to also account for multi-scale density correlations which appear in the calculation. We thus derive an accurate expression for the mass function which is based on approximations that are valid over a larger range of mass scales and redshifts than those of other authors. By tracking the terms ignored in the analysis, we estimate theoretical errors for our result and also for the results of others. We also discuss the complications introduced by the choice of smoothing filter function, which we take to be a top-hat in real space, and which leads to the dominant errors in our expression. Finally, we present a detailed comparison between the various expressions for the mass functions, exploring the accuracy and range of validity of each.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Excursion Sets and Non-Gaussian Void Statistics

Primordial non-Gaussianity (NG) affects the large scale structure (LSS) of the universe by leaving an imprint on the distribution of matter at late times. Much attention has been focused on using the distribution of collapsed objects (i.e. dark matter halos and the galaxies and galaxy clusters that reside in them) to probe primordial NG. An equally interesting and complementary probe however is the abundance of extended underdense regions or voids in the LSS. The calculation of the abundance of voids using the excursion set formalism in the presence of primordial NG is subject to the same technical issues as the one for halos, which were discussed e.g. in arXiv:1005.1203. However, unlike the excursion set problem for halos which involved random walks in the presence of one barrier $δ_c$, the void excursion set problem involves two barriers $δ_v$ and $δ_c$. This leads to a new complication introduced by what is called the "void-in-cloud" effect discussed in the literature, which is unique to the case of voids. We explore a path integral approach which allows us to carefully account for all these issues, leading to a rigorous derivation of the effects of primordial NG on void abundances. The void-in-cloud issue in particular makes the calculation conceptually rather different from the one for halos. However, we show that its final effect can be described by a simple yet accurate approximation. Our final void abundance function is valid on larger scales than the expressions of other authors, while being broadly in agreement with those expressions on smaller scales.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-Gaussian halo abundances in the excursion set approach with correlated steps

We study the effects of primordial non-Gaussianity on the large scale structure in the excursion set approach, accounting for correlations between steps of the random walks in the smoothed initial density field. These correlations are induced by realistic smoothing filters (as opposed to a filter that is sharp in k-space), but have been ignored by many analyses to date. We present analytical arguments -- building on existing results for Gaussian initial conditions -- which suggest that the effect of the filter at large smoothing scales is remarkably simple, and is in fact identical to what happens in the Gaussian case: the non-Gaussian walks behave as if they were smooth and deterministic, or "completely correlated". As a result, the first crossing distribution (which determines, e.g., halo abundances) follows from the single-scale statistics of the non-Gaussian density field -- the so-called "cloud-in-cloud" problem does not exist for completely correlated walks. Also, the answer from single-scale statistics is simply one half that for sharp-k walks. We explicitly test these arguments using Monte Carlo simulations of non-Gaussian walks, showing that the resulting first crossing distributions, and in particular the factor 1/2 argument, are remarkably insensitive to variations in the power spectrum and the defining non-Gaussian process. We also use our Monte Carlo walks to test some of the existing prescriptions for the non-Gaussian first crossing distribution. Since the factor 1/2 holds for both Gaussian and non-Gaussian initial conditions, it provides a theoretical motivation (the first, to our knowledge) for the common practice of analytically prescribing a ratio of non-Gaussian to Gaussian halo abundances.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Extreme Tail of the Non-Gaussian Mass Function

Number counts of massive high-redshift clusters provide a window to study primordial non-Gaussianity. The current quality of data, however, forces the statistical analysis to probe a region of parameter space -- the extreme tail of the mass function -- which is neither accessible in any of the currently available theoretical prescriptions for calculating the mass function, nor calibrated in N-body simulations. In this work we present a new analytical prescription for calculating a "resummed" non-Gaussian halo mass function, which is constructed to remain stable in the extreme tail. We show that the prescription works well in the parameter regime that has been currently explored in simulations. We then use Fisher matrix techniques to compare our prescription with an extrapolated fit to N-body simulations, which has recently been used to obtain constraints from data collected by the South Pole Telecope. We show that for the current data, both prescriptions would lead to statistically consistent constraints. As the data improve, however, there is a possibility of introducing a statistically significant bias in the constraints due to the choice of prescription, especially if non-Gaussianity is scale dependent and becomes relatively large on cluster scales. It would then be necessary to test the accuracy of the prescriptions in N-body simulations that can probe clusters with high masses and redshifts in the presence of large non-Gaussianity.

0 citations0 reviewsNo code linkedOpen access