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

Robert E. Smith

Robert E. Smith 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.GAgr-qcastro-phhep-phphysics.soc-phSocial and Information Networkscs.CYhep-th

Catalog footprint

What is connected

35works
9topics
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

35 published item(s)

preprint2020arXiv

A network perspective on intermedia agenda-setting

In Communication Theory, intermedia agenda-setting refers to the influence that different news sources may have on each other, and how this subsequently affects the breadth of information that is presented to the public. Several studies have attempted to quantify the impact of intermedia agenda-setting in specific countries or contexts, but a large-scale, data-driven investigation is still lacking. Here, we operationalise intermedia agenda-setting by putting forward a methodology to infer networks of influence between different news sources on a given topic, and apply it on a large dataset of news articles published by globally and locally prominent news organisations in 2016. We find influence to be significantly topic-dependent, with the same news sources acting as agenda-setters (i.e., central nodes) with respect to certain topics and as followers (i.e., peripheral nodes) with respect to others. At the same time, we find that the influence networks associated with most topics exhibit small world properties, which we find to play a significant role towards the overall diversity of sentiment expressed about the topic by the news sources in the network. In particular, we find clustering and density of influence networks to act as competing forces in this respect, with the former increasing and the latter reducing diversity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

When facts fail: Bias, polarisation and truth in social networks

Online social networks provide users with unprecedented opportunities to engage with diverse opinions. At the same time, they enable confirmation bias on large scales by empowering individuals to self-select narratives they want to be exposed to. A precise understanding of such tradeoffs is still largely missing. We introduce a social learning model where most participants in a network update their beliefs unbiasedly based on new information, while a minority of participants reject information that is incongruent with their preexisting beliefs. This simple mechanism generates permanent opinion polarization and cascade dynamics, and accounts for the aforementioned tradeoff between confirmation bias and social connectivity through analytic results. We investigate the model's predictions empirically using US county-level data on the impact of Internet access on the formation of beliefs about global warming. We conclude by discussing policy implications of our model, highlighting the downsides of debunking and suggesting alternative strategies to contrast misinformation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Matter power spectrum and the challenge of percent accuracy

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for $k\leq 2$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than $L=0.5$ $h^{-1}\rm Gpc$ to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An exploration of galaxy-galaxy lensing and galaxy clustering in the Millennium-XXL simulation

The combination of galaxy-galaxy lensing and galaxy clustering data has the potential to simultaneously constrain both the cosmological and galaxy formation models. In this paper we perform a comprehensive exploration of these signals and their covariances through a combination of analytic and numerical approaches. First, we derive analytic expressions for the projected galaxy correlation function and stacked tangential shear profile and their respective covariances, which include Gaussian and discreteness noise terms. Secondly, we measure these quantities from mock galaxy catalogues obtained from the Millennium-XXL simulation and semi-analytic models of galaxy formation. We find that on large scales (R>10 Mpc/h), the galaxy bias is roughly linear and deterministic. On smaller scales (R<5 Mpc/h) the bias is a complicated function of scale and luminosity, determined by the different spatial distribution and abundance of satellite galaxies present when different magnitude cuts are applied, as well as by the mass dependence of the host haloes on magnitude. Our theoretical model for the covariances provides a reasonably good description of the measured ones on small and large scales. However, on intermediate scales (1<R<10 Mpc/h), the predicted errors are ~2-3 times smaller, suggesting that the inclusion of higher-order, non-Gaussian terms in the covariance will be required for further improvements. Importantly, both our theoretical and numerical methods show that the galaxy-galaxy lensing and clustering signals have a non-zero cross-covariance matrix with significant bin-to-bin correlations. Future surveys aiming to combine these probes must take this into account in order to obtain unbiased and realistic constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modified Gravity N-body Code Comparison Project

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $Λ$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Anisotropic Line Correlation Function as a Probe of Anisotropies in Galaxy Surveys

We propose an anisotropic generalisation of the line correlation function (ALCF) to separate and quantify phase information in the large-scale structure of galaxies. The line correlation function probes the strictly non-linear regime of structure formation and since phase information drops out of the power spectrum, the line correlation function provides a complementary tool to commonly used techniques based on two-point statistics. Furthermore, it is independent of linear bias as well as the Gaussian variance on the modulus of the density field and thus may also prove to be advantageous compared to the bispectrum or similar higher-order statistics for certain cases. For future applications it is vital, though, to be able to account for observational effects that cause anisotropies in the distribution of galaxies. Based on a number of numerical studies, we find that our ALCF is well suited to accomplish this task and we demonstrate how the Alcock-Paczynski effect and kinematical redshift-space distortions can in principle be measured via the ALCF.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The same with less: The cosmic web of warm versus cold dark matter dwarf galaxies

We explore fundamental properties of the distribution of low mass dark matter halos within the cosmic web using warm dark matter (WDM) and cold dark matter (CDM) cosmological simulations. Using self abundance-matched mock galaxy catalogs, we show that the distribution of dwarf galaxies in a WDM universe, wherein low mass halo formation is heavily suppressed, is nearly indistinguishable to that of a CDM universe whose low mass halos are not seen because galaxy formation is suppressed below some threshold halo mass. However, if the scatter between dwarf galaxy luminosity and halo properties is large enough, low mass CDM halos would sometimes host relatively bright galaxies thereby populating CDM voids with the occasional isolated galaxy and reducing the numbers of completely empty voids. Otherwise, without high mass to light scatter, all mock galaxy clustering statistics that we consider--the auto-correlation function, the numbers and radial profiles of satellites, the numbers of isolated galaxies, and the PDF of small voids--are nearly identical in CDM and WDM. WDM voids are neither larger nor emptier than CDM voids, when constructed from abundance-matched halo catalogs. It is thus a challenge to determine whether the CDM problem of the over-abundance of small halos with respect to the number density of observed dwarf galaxies has a cosmological solution or an astrophysical solution. However, some clues about the dark matter particle and the scatter between the properties of dwarf galaxies and their dark matter halo hosts might be found in the cosmic web of galaxies in future surveys of the local volume.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Toward optimal cluster power spectrum analysis

The power spectrum of galaxy clusters is an important probe of the cosmological model. In this paper we determine the optimal weighting scheme for maximizing the signal-to-noise ratio for such measurements. We find a closed form analytic expression for the optimal weights. Our expression takes into account: cluster mass, finite survey volume effects, survey masking, and a flux limit. The implementation of this weighting scheme requires knowledge of the measured cluster masses, and analytic models for the bias and space-density of clusters as a function of mass and redshift. Recent studies have suggested that the optimal method for reconstruction of the matter density field from a set of clusters is mass-weighting (Seljak et al 2009, Hamaus et al 2010, Cai et al 2011). We compare our optimal weighting scheme with this approach and also with the original power spectrum scheme of Feldman et al (1994). We show that our optimal weighting scheme outperforms these approaches for both volume- and flux-limited cluster surveys. Finally, we present a new expression for the Fisher information matrix for cluster power spectrum analysis. Our expression shows that for an optimally weighted cluster survey the cosmological information content is boosted, relative to the standard approach of Tegmark (1997).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards optimal estimation of the galaxy power spectrum

The galaxy power spectrum encodes a wealth of information about cosmology and the matter fluctuations. Its unbiased and optimal estimation is therefore of great importance. In this paper we generalise the framework of Feldman et al. (1994) to take into account the fact that galaxies are not simply a Poisson sampling of the underlying dark matter distribution. Besides finite survey-volume effects and flux-limits, our optimal estimation scheme incorporates several of the key tenets of galaxy formation: galaxies form and reside exclusively in dark matter haloes; a given dark matter halo may host several galaxies of various luminosities; galaxies inherit part of their large-scale bias from their host halo. Under these broad assumptions, we prove that the optimal weights "do not" explicitly depend on galaxy luminosity, other than through defining the maximum survey volume and effective galaxy density at a given position. Instead, they depend on the bias associated with the host halo; the first and second factorial moments of the halo occupation distribution; a selection function, which gives the fraction of galaxies that can be observed in a halo of mass M at position r in the survey; and an effective number density of galaxies. If one wishes to reconstruct the matter power spectrum, then, provided the model is correct, this scheme provides the only unbiased estimator. The practical challenges with implementing this approach are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

What is the optimal way to measure the galaxy power spectrum?

Measurements of the galaxy power spectrum contain a wealth of information about the Universe. Its optimal extraction is vital if we are to truly understand the micro-physical nature of dark matter and dark energy. In Smith & Marian (2015) we generalized the power spectrum methodology of Feldman et al. (1994) to take into account the key tenets of galaxy formation: galaxies form and reside exclusively in dark matter haloes; a given dark matter halo may host galaxies of various luminosities; galaxies inherit the large-scale bias associated with their host halo. In this paradigm we derived the optimal weighting and reconstruction scheme for maximizing the signal-to-noise on a given band power estimate. For a future all-sky flux-limited galaxy redshift survey of depth b_J ~22, we now demonstrate that the optimal weighting scheme does indeed provide improved S/N at the level of ~20% when compared to Feldman et al. (1994) and ~60% relative to Percival et al. (2003), for scales of order k~0.5 Mpc/h. Using a Fisher matrix approach, we show that the cosmological information yield is also increased relative to these alternate methods -- especially the primordial power spectrum amplitude and dark energy equation of state.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A new method to measure galaxy bias

We present a new approach for modelling galaxy/halo bias that utilizes the full non-linear information contained in the moments of the matter density field, which we derive using a set of numerical simulations. Although our method is general, we perform a case study based on the local Eulerian bias scheme truncated to second-order. Using 200 N-body simulations covering a total comoving volume of 675 (Gpc/h)^3, we measure several 2- and 3-point statistics of the halo distribution to unprecedented accuracy. We use the bias model to fit the halo-halo power spectrum, the halo-matter cross spectrum and the corresponding three bispectra for wavenumbers in the range 0.04 < k < 0.12 h Mpc^-1. We find the constraints on the bias parameters obtained using the full non-linear information differ significantly from those derived using standard perturbation theory at leading order. Hence, neglecting the full non-linear information leads to biased results for this particular scale range. We also test the validity of the second-order Eulerian local biasing scheme by comparing the parameter constraints derived from different statistics. Analysis of the halo-matter cross-correlation coefficients defined for the 2- and 3-point statistics reveals further inconsistencies contained in the second-order Eulerian bias scheme, suggesting it is too simple a model to describe halo bias with high accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Precision cosmology in muddy waters: Cosmological constraints and N-body codes

Future large-scale structure surveys of the Universe will aim to constrain the cosmological model and the true nature of dark energy with unprecedented accuracy. In order for these surveys to achieve their designed goals, they will require predictions for the nonlinear matter power spectrum to sub-percent accuracy. Through the use of a large ensemble of cosmological N-body simulations, we demonstrate that if we do not understand the uncertainties associated with simulating structure formation, i.e. knowledge of the `true' simulation parameters, and simply seek to marginalize over them, then the constraining power of such future surveys can be significantly reduced. However, for the parameters {n_s, h, Om_b, Om_m}, this effect can be largely mitigated by adding the information from a CMB experiment, like Planck. In contrast, for the amplitude of fluctuations sigma8 and the time-evolving equation of state of dark energy {w_0, w_a}, the mitigation is mild. On marginalizing over the simulation parameters, we find that the dark-energy figure of merit can be degraded by ~2. This is likely an optimistic assessment, since we do not take into account other important simulation parameters. A caveat is our assumption that the Hessian of the likelihood function does not vary significantly when moving from our adopted to the 'true' simulation parameter set. This paper therefore provides strong motivation for rigorous convergence testing of N-body codes to meet the future challenges of precision cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cosmological parameter constraints from galaxy-galaxy lensing and galaxy clustering with the SDSS DR7

Recent studies have shown that the cross-correlation coefficient between galaxies and dark matter is very close to unity on scales outside a few virial radii of galaxy halos, independent of the details of how galaxies populate dark matter halos. This finding makes it possible to determine the dark matter clustering from measurements of galaxy-galaxy weak lensing and galaxy clustering. We present new cosmological parameter constraints based on large-scale measurements of spectroscopic galaxy samples from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We generalise the approach of Baldauf et al. (2010) to remove small scale information (below 2 and 4 Mpc/h for lensing and clustering measurements, respectively), where the cross-correlation coefficient differs from unity. We derive constraints for three galaxy samples covering 7131 sq. deg., containing 69150, 62150, and 35088 galaxies with mean redshifts of 0.11, 0.28, and 0.40. We clearly detect scale-dependent galaxy bias for the more luminous galaxy samples, at a level consistent with theoretical expectations. When we vary both σ_8 and Ω_m (and marginalise over non-linear galaxy bias) in a flat LCDM model, the best-constrained quantity is σ_8 (Ω_m/0.25)^{0.57}=0.80 +/- 0.05 (1-sigma, stat. + sys.), where statistical and systematic errors have comparable contributions, and we fixed n_s=0.96 and h=0.7. These strong constraints on the matter clustering suggest that this method is competitive with cosmic shear in current data, while having very complementary and in some ways less serious systematics. We therefore expect that this method will play a prominent role in future weak lensing surveys. When we combine these data with WMAP7 CMB data, constraints on σ_8, Ω_m, H_0, w_{de} and \sum m_ν become 30--80 per cent tighter than with CMB data alone, since our data break several parameter degeneracies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Halo Mass Function and the Free Streaming Scale

The nature of structure formation around the particle free streaming scale is still far from understood. Many attempts to simulate hot, warm, and cold dark matter cosmologies with a free streaming cutoff have been performed with cosmological particle-based simulations, but they all suffer from spurious structure formation at scales below their respective free streaming scales -- i.e. where the physics of halo formation is most affected by free streaming. We perform a series of high resolution numerical simulations of different WDM models, and develop an approximate method to subtract artificial structures in the measured halo mass function. The corrected measurements are then used to construct and calibrate an extended Press-Schechter (EPS) model with sharp-$k$ window function and adequate mass assignment. The EPS model gives accurate predictions for the low redshift halo mass function of CDM and WDM models, but it significantly under-predicts the halo abundance at high redshifts. By taking into account the ellipticity of the initial patches and connecting the characteristic filter scale to the smallest ellipsoidal axis, we are able to eliminate this inconsistency and obtain an accurate mass function over all redshifts and all dark matter particle masses covered by the simulations. As an additional application we use our model to predict the microhalo abundance of the standard neutralino-CDM scenario and we give the first quantitative prediction of the mass function over the full range of scales of CDM structure formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Halo Stochasticity from Exclusion and non-linear Clustering

The clustering of galaxies in ongoing and upcoming galaxy surveys contains a wealth of cosmological information, but extracting this information is a non-trivial task since galaxies and their host haloes are stochastic tracers of the matter density field. This stochasticity is usually modeled as Poisson shot noise, which is constant as a function of wavenumber with amplitude given by 1/n, where n is the number density of galaxies. Here we use dark matter haloes in N-body simulations to show evidence for deviations from this simple behaviour and develop models that explain the behaviour of the large scale stochasticity. First, haloes are extended, non-overlapping objects, i.e., their correlation function needs to go to -1 on small scales. This leads to a negative correction to the stochasticity relative to the Poisson value at low wavenumber k, decreasing to zero for wavenumbers large compared to the inverse exclusion scale. Second, haloes show a non-linear enhancement of clustering outside the exclusion scale, leading to a positive stochasticity correction. Both of these effects go to zero for high-k, making the stochasticity scale dependent even for k<0.1 h/Mpc. We show that the corrections in the low-k regime are the same in Eulerian and Lagrangian space, but that the transition scale is pushed to smaller scales for haloes observed at present time, relative to the initial conditions. These corrections vary with halo mass and redshift. We also discuss simple applications of these effects to the galaxy samples with non-vanishing satellite fraction, where the stochasticity can again deviate strongly from the fiducial Poisson expectation. Overall these effects affect the clustering of galaxies at a level of a few percent even on very large scales and need to be modelled properly if we want to extract high precision cosmological information from the upcoming galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The cosmological information of shear peaks: beyond the abundance

We study the cosmological information of weak lensing (WL) peaks, focusing on two other statistics besides their abundance: the stacked tangential-shear profiles and the peak-peak correlation function. We use a large ensemble of simulated WL maps with survey specifications relevant to future missions like Euclid and LSST, to explore the three peak probes. We find that the correlation function of peaks with high signal-to-noise (S/N) measured from fields of size 144 sq. deg. has a maximum of ~0.3 at an angular scale ~10 arcmin. For peaks with smaller S/N, the amplitude of the correlation function decreases, and its maximum occurs on smaller angular scales. We compare the peak observables measured with and without shape noise and find that for S/N~3 only ~5% of the peaks are due to large-scale structures, the rest being generated by shape noise. The covariance matrix of the probes is examined: the correlation function is only weakly covariant on scales < 30 arcmin, and slightly more on larger scales; the shear profiles are very correlated for theta > 2 arcmin, with a correlation coefficient as high as 0.7. Using the Fisher-matrix formalism, we compute the cosmological constraints for {Om_m, sig_8, w, n_s} considering each probe separately, as well as in combination. We find that the correlation function of peaks and shear profiles yield marginalized errors which are larger by a factor of 2-4 for {Om_m, sig_8} than the errors yielded by the peak abundance alone, while the errors for {w, n_s} are similar. By combining the three probes, the marginalized constraints are tightened by a factor of ~2 compared to the peak abundance alone, the least contributor to the error reduction being the correlation function. This work therefore recommends that future WL surveys use shear peaks beyond their abundance in order to constrain the cosmological model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Toward an accurate mass function for precision cosmology

Cosmological surveys aim to use the evolution of the abundance of galaxy clusters to accurately constrain the cosmological model. In the context of LCDM, we show that it is possible to achieve the required percent level accuracy in the halo mass function with gravity-only cosmological simulations, and we provide simulation start and run parameter guidelines for doing so. Some previous works have had sufficient statistical precision, but lacked robust verification of absolute accuracy. Convergence tests of the mass function with, for example, simulation start redshift can exhibit false convergence of the mass function due to counteracting errors, potentially misleading one to infer overly optimistic estimations of simulation accuracy. Percent level accuracy is possible if initial condition particle mapping uses second order Lagrangian Perturbation Theory, and if the start epoch is between 10 and 50 expansion factors before the epoch of halo formation of interest. The mass function for halos with fewer than ~1000 particles is highly sensitive to simulation parameters and start redshift, implying a practical minimum mass resolution limit due to mass discreteness. The narrow range in converged start redshift suggests that it is not presently possible for a single simulation to capture accurately the cluster mass function while also starting early enough to model accurately the numbers of reionisation era galaxies, whose baryon feedback processes may affect later cluster properties. Ultimately, to fully exploit current and future cosmological surveys will require accurate modeling of baryon physics and observable properties, a formidable challenge for which accurate gravity-only simulations are just an initial step.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cluster Density Profiles as a Test of Modified Gravity

We present a new test of gravitational interactions at the r\sim(0.2-20)Mpc scale, around the virial radius of dark matter halos measured through cluster-galaxy lensing of maxBCG clusters from the Sloan Digital Sky Survey (SDSS). We employ predictions from self-consistent simulations of f(R) gravity to find an upper bound on the background field amplitude of f_R0<3.5x10^-3 at the 1D-marginalized 95% confidence level. As a model-independent assessment of the constraining power of cluster profiles measured through weak gravitational lensing, we also constrain the amplitude F_0 of a phenomenological modification based on the profile enhancement induced by f(R) gravity when not including effects from the increased cluster abundance in f(R). In both scenarios, dark-matter-only simulations of the concordance model corresponding to f_R0=0 and F_0=0 are consistent with the lensing measurements, i.e., at the 68% and 95% confidence level, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How covariant is the galaxy luminosity function?

We investigate the error properties of certain galaxy luminosity function (GLF) estimators. Using a cluster expansion of the density field, we show how, for both volume and flux limited samples, the GLF estimates are covariant. The covariance matrix can be decomposed into three pieces: a diagonal term arising from Poisson noise; a sample variance term arising from large-scale structure in the survey volume; an occupancy covariance term arising due to galaxies of different luminosities inhabiting the same cluster. To evaluate the theory one needs: the mass function and bias of clusters, and the conditional luminosity function (CLF). We use a semi-analytic model (SAM) galaxy catalogue from the Millennium run N-body simulation and the CLF of Yang et al. (2003) to explore these effects. The GLF estimates from the SAM and the CLF qualitatively reproduce results from the 2dFGRS. We also measure the luminosity dependence of clustering in the SAM and find reasonable agreement with 2dFGRS results for bright galaxies. However, for fainter galaxies, L<L*, the SAM overpredicts the relative bias by ~10-20%. We use the SAM data to estimate the errors in the GLF estimates for a volume limited survey of volume V~0.13 [Gpc/h]^3. We find that different luminosity bins are highly correlated: for L<L* the correlation coefficient is r>0.5. Our theory is in good agreement with these measurements. These strong correlations can be attributed to sample variance. For a flux-limited survey of similar volume, the estimates are only slightly less correlated. We explore the importance of these effects for GLF model parameter estimation. We show that neglecting to take into account the bin-to-bin covariances can lead to significant systematic errors in best-fit parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measuring primordial non-Gaussianity through weak lensing peak counts

We explore the possibility of detecting primordial non-Gaussianity of the local type using weak lensing peak counts. We measure the peak abundance in sets of simulated weak lensing maps corresponding to three models f_NL={0, +100, -100}. Using survey specifications similar to those of Euclid and without assuming any knowledge of the lens and source redshifts, we find the peak functions of the non-Gaussian models with f_NL=+-100 to differ by up to 15% from the Gaussian peak function at the high-mass end. For the assumed survey parameters, the probability of fitting an f_NL=0 peak function to the f_NL=+-100 peak functions is less than 0.1%. Assuming the other cosmological parameters known, f_NL can be measured with an error Δf_NL ~ 13. It is therefore possible that future weak lensing surveys like Euclid and LSST may detect primordial non-Gaussianity from the abundance of peak counts, and provide complementary information to that obtained from the cosmic microwave background.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measuring primordial non-Gaussianity with weak-lensing surveys

We study the ability of future weak lensing (WL) surveys to constrain primordial non-Gaussianity of the local type. We use a large ensemble of simulated WL maps with survey specifications relevant to Euclid and LSST. The simulations assume Cold Dark Matter cosmologies that vary certain parameters around fiducial values: the non-Gaussianity parameter f_NL, the matter density parameter Omega_m, the amplitude of the matter power spectrum sigma_8, the spectral index of the primordial power spectrum n_s, and the dark-energy equation-of-state parameter w_0. We assess the sensitivity of the cosmic shear correlation functions, the third-order aperture mass statistics, and the abundance of shear peaks to these parameters. We find that each of the considered probes provides unmarginalized constraints of Delta f_NL ~ 20 on f_NL. Marginalized constraints from any individual WL probe are much weaker due to strong correlations between parameters. However, the parameter errors can be substantially reduced by combining information from different WL probes. Combining all WL probes yields the following marginal (68% confidence level) uncertainties: Delta f_NL ~ 50, Delta Omega_m ~ 0.002, Delta sigma_8 ~ 0.004, Delta n_s ~ 0.007, and Delta w_0 ~ 0.03. We examine the bias induced by neglecting f_NL on the constraints on the other parameters. We find sigma_8 and w_0 to be the most affected. Moreover, neglecting non-Gaussianity leads to a severe underestimation of the uncertainties in the other cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear Evolution of Cosmological Structures in Warm Dark Matter Models

The dark energy dominated warm dark matter (WDM) model is a promising alternative cosmological scenario. We explore large-scale structure formation in this paradigm. We do this in two different ways: with the halo model approach and with the help of an ensemble of high resolution N-body simulations. Combining these quasi-independent approaches, leads to a physical understanding of the important processes which shape the formation of structures. We take a detailed look at the halo mass function, the concentrations and the linear halo bias of WDM. In all cases we find interesting deviations with respect to CDM. In particular, the concentration-mass relation displays a turnover for group scale dark matter haloes, for the case of WDM particles with masses of the order ~0.25 keV. This may be interpreted as a hint for top-down structure formation on small scales. We implement our results into the halo model and find much better agreement with simulations. On small scales the WDM halo model now performs as well as its CDM counterpart.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimized detection of shear peaks in weak lensing maps

We present a new method to extract cosmological constraints from weak lensing (WL) peak counts, which we denote as `the hierarchical algorithm'. The idea of this method is to combine information from WL maps sequentially smoothed with a series of filters of different size, from the largest down to the smallest, thus increasing the cosmological sensitivity of the resulting peak function. We compare the cosmological constraints resulting from the peak abundance measured in this way and the abundance obtained by using a filter of fixed size, which is the standard practice in WL peak studies. For this purpose, we employ a large set of WL maps generated by ray-tracing through N-body simulations, and the Fisher matrix formalism. We find that if low-S/N peaks are included in the analysis (S/N ~ 3), the hierarchical method yields constraints significantly better than the single-sized filtering. For a large future survey such as Euclid or LSST, combined with information from a CMB experiment like Planck, the results for the hierarchical (single-sized) method are: Δn=0.0039 (0.004); ΔΩm=0.002 (0.0045); Δσ8=0.003 (0.006); Δw=0.019 (0.0525). This forecast is conservative, as we assume no knowledge of the redshifts of the lenses, and consider a single broad bin for the redshifts of the sources. If only peaks with S/N >= 6 are considered, then there is little difference between the results of the two methods. We also examine the statistical properties of the hierarchical peak function: Its covariance matrix has off-diagonal terms for bins with S/N <= 6 and aperture mass of M < 3 x 1e+14 Ms/h, the higher bins being largely uncorrelated and therefore well described by a Poisson distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modelling large-scale halo bias using the bispectrum

We study the relation between the halo and matter density fields -- commonly termed bias -- in the LCDM framework. In particular, we examine the local model of biasing at quadratic order in the matter density. This model is characterized by parameters b_1 and b_2. Using an ensemble of N-body simulations, we apply several statistical methods to estimate the parameters. We measure halo and matter fluctuations smoothed on various scales and find that the parameters vary with smoothing scale. We argue that, for real-space measurements, owing to the mixing of wavemodes, no scale can be found for which the parameters are independent of smoothing. However, this is not the case in Fourier space. We measure halo power spectra and construct estimates for an effective large-scale bias. We measure the configuration dependence of the halo bispectra B_hhh and reduced bispectra Q_hhh for very large-scale k-space triangles. From this we constrain b_1 and b_2. Using the lowest-order perturbation theory, we find that for B_hhh the best-fit parameters are in reasonable agreement with one another as the triangle scale is varied, but that the fits become poor as smaller scales are included. The same is true for Q_hhh. The best-fit parameters depend on the discreteness correction. This led us to consider halo-mass cross-bispectra. The results from these statistics support our earlier findings. We develop a test to explore the importance of missing higher-order terms in the models. We prove that low-order expansions are not able to correctly model the data, even on scales k_1~0.04 h/Mpc. If robust inferences are to be drawn from galaxy surveys, then accurate models for the full nonlinear matter bispectrum and trispectrum will be essential.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonlinear clustering in models with primordial non-Gaussianity: the halo model approach

We develop the halo model of large-scale structure as an accurate tool for probing primordial non-Gaussianity. In this study we focus on understanding the matter clustering at several redshifts. The primordial non-Gaussianity is modeled as a quadratic correction to the local Gaussian potential, and is characterized by the parameter f_NL. In our formulation of the halo model we pay special attention to the effect of halo exclusion, and show that this can potentially solve the long standing problem of excess power on large scales in this model. The model depends on the mass function, clustering and density profiles of halos. We test these ingredients using a large ensemble of high-resolution Gaussian and non-Gaussian numerical simulations. In particular, we provide a first exploration of how density profiles change in the presence of primordial non-Gaussianities. We find that for f_NL positive/negative high mass halos have an increased/decreased core density, so being more/less concentrated than in the Gaussian case. We also examine the halo bias and show that, if the halo model is correct, then there is a small asymmetry in the scale-dependence of the bias on very large scales, which arises because the Gaussian bias must be renormalized. We show that the matter power spectrum is modified by ~2.5% and ~3.5% on scales k~1.0 h/Mpc at z=0 and z=1, respectively. Our halo model calculation reproduces the absolute amplitude to within 10% and the ratio of non-Gaussian to Gaussian spectra to within 1%. We also measure the matter correlation functions and find similarly good agreement between the model and the data. We anticipate that this modeling will be useful for constraining f_NL from measurements of the shear correlation function in future weak lensing surveys such as Euclid.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing the Warm Dark Matter paradigm with large-scale structures

We explore the impact of a LWDM cosmological scenario on the clustering properties of large-scale structure in the Universe. We do this by extending the halo model. The new development is that we consider two components to the mass density: one arising from mass in collapsed haloes, and the second from a smooth component of uncollapsed mass. Assuming that the nonlinear clustering of dark matter haloes can be understood, then from conservation arguments one can precisely calculate the clustering properties of the smooth component and its cross-correlation with haloes. We then explore how the three main ingredients of the halo calculations, the mass function, bias and density profiles are affected by WDM. We show that, relative to CDM: the mass function is suppressed by ~50%, for masses ~100 times the free-streaming mass-scale; the bias of low mass haloes can be boosted by up to 20%; core densities of haloes can be suppressed. We also examine the impact of relic thermal velocities on the density profiles, and find that these effects are constrained to scales r<1 kpc/h, and hence of little importance for dark matter tests, owing to uncertainties in the baryonic physics. We use our modified halo model to calculate the non-linear matter power spectrum, and find significant small-scale power in the model. However, relative to the CDM case, the power is suppressed. We then calculate the expected signal and noise that our set of LWDM models would give for a future weak lensing mission. We show that the models should in principle be separable at high significance. Finally, using the Fisher matrix formalism we forecast the limit on the WDM particle mass for a future full-sky weak lensing mission like Euclid or LSST. With Planck priors and using multipoles l<5000, we find that a lower limit of 2.6 keV should be easily achievable.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

What do cluster counts really tell us about the Universe?

We study the covariance matrix of the cluster mass function in cosmology. We adopt a two-line attack: firstly, we employ the counts-in-cells framework to derive an analytic expression for the covariance of the mass function. Secondly, we use a large ensemble of N-body simulations in the LCDM framework to test this. Our theoretical results show that the covariance can be written as the sum of two terms: a Poisson term, which dominates in the limit of rare clusters; and a sample variance term, which dominates for more abundant clusters. Our expressions are analogous to those of Hu & Kravtsov (2003) for multiple cells and a single mass tracer. Calculating the covariance depends on: the mass function and bias of clusters, and the variance of mass fluctuations within the survey volume. The predictions show that there is a strong bin-to-bin covariance between measurements. In terms of the cross-correlation coefficient, we find r~0.5 for haloes with M<3e14 Msol at z=0. Comparison of these predictions with estimates from simulations shows excellent agreement. We use the Fisher matrix formalism to explore the cosmological information content of the counts. We compare the Poisson likelihood model, with the more realistic likelihood model of Lima & Hu (2004), and all terms entering the Fisher matrices are evaluated using the simulations. We find that the Poisson approximation should only be used for the rarest objects, M>3e14 Msol, otherwise the information content of a survey of size V~13.5 [Gpc/h]^3 would be overestimated, resulting in errors that are ~2 times smaller. As an auxiliary result, we show that the bias of clusters, obtained from the cluster-mass cross-variance, is linear on scales >50 Mpc/h, whereas that obtained from the auto-variance is nonlinear.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An algorithm for the direct reconstruction of the dark matter correlation function from weak lensing and galaxy clustering

The clustering of matter on cosmological scales is an essential probe for studying the physical origin and composition of our Universe. To date, most of the direct studies have focused on shear-shear weak lensing correlations, but it is also possible to extract the dark matter clustering by combining galaxy-clustering and galaxy-galaxy-lensing measurements. In this study we develop a method that can constrain the dark matter correlation function from galaxy clustering and galaxy-galaxy-lensing measurements, by focusing on the correlation coefficient between the galaxy and matter overdensity fields. To generate a mock galaxy catalogue for testing purposes, we use the Halo Occupation Distribution approach applied to a large ensemble of N-body simulations to model pre-existing SDSS Luminous Red Galaxy sample observations. Using this mock catalogue, we show that a direct comparison between the excess surface mass density measured by lensing and its corresponding galaxy clustering quantity is not optimal. We develop a new statistic that suppresses the small-scale contributions to these observations and show that this new statistic leads to a cross-correlation coefficient that is within a few percent of unity down to 5 Mpc/h. Furthermore, the residual incoherence between the galaxy and matter fields can be explained using a theoretical model for scale-dependent bias, giving us a final estimator that is unbiased to within 1%. We also perform a comprehensive study of other physical effects that can affect the analysis, such as redshift space distortions and differences in radial windows between galaxy clustering and weak lensing observations. We apply the method to a range of cosmological models and show the viability of our new statistic to distinguish between cosmological models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Confirmation of general relativity on large scales from weak lensing and galaxy velocities

Although general relativity underlies modern cosmology, its applicability on cosmological length scales has yet to be stringently tested. Such a test has recently been proposed, using a quantity, EG, that combines measures of large-scale gravitational lensing, galaxy clustering and structure growth rate. The combination is insensitive to 'galaxy bias' (the difference between the clustering of visible galaxies and invisible dark matter) and is thus robust to the uncertainty in this parameter. Modified theories of gravity generally predict values of EG different from the general relativistic prediction because, in these theories, the 'gravitational slip' (the difference between the two potentials that describe perturbations in the gravitational metric) is non-zero, which leads to changes in the growth of structure and the strength of the gravitational lensing effect3. Here we report that EG = 0.39 +/- 0.06 on length scales of tens of megaparsecs, in agreement with the general relativistic prediction of EG $\approx$ 0.4. The measured value excludes a model within the tensor-vector-scalar gravity theory, which modifies both Newtonian and Einstein gravity. However, the relatively large uncertainty still permits models within f(R) theory, which is an extension of general relativity. A fivefold decrease in uncertainty is needed to rule out these models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Minimizing the stochasticity of halos in large-scale structure surveys

In recent work (Seljak, Hamaus and Desjacques 2009) it was found that weighting central halo galaxies by halo mass can significantly suppress their stochasticity relative to the dark matter, well below the Poisson model expectation. In this paper we extend this study with the goal of finding the optimal mass-dependent halo weighting and use $N$-body simulations to perform a general analysis of halo stochasticity and its dependence on halo mass. We investigate the stochasticity matrix, defined as $C_{ij}\equiv<(δ_i -b_iδ_m)(δ_j-b_jδ_m)>$, where $δ_m$ is the dark matter overdensity in Fourier space, $δ_i$ the halo overdensity of the $i$-th halo mass bin and $b_i$ the halo bias. In contrast to the Poisson model predictions we detect nonvanishing correlations between different mass bins. We also find the diagonal terms to be sub-Poissonian for the highest-mass halos. The diagonalization of this matrix results in one large and one low eigenvalue, with the remaining eigenvalues close to the Poisson prediction $1/\bar{n}$, where $\bar{n}$ is the mean halo number density. The eigenmode with the lowest eigenvalue contains most of the information and the corresponding eigenvector provides an optimal weighting function to minimize the stochasticity between halos and dark matter. We find this optimal weighting function to match linear mass weighting at high masses, while at the low-mass end the weights approach a constant whose value depends on the low-mass cut in the halo mass function. Finally, we employ the halo model to derive the stochasticity matrix and the scale-dependent bias from an analytical perspective. It is remarkably successful in reproducing our numerical results and predicts that the stochasticity between halos and the dark matter can be reduced further when going to halo masses lower than we can resolve in current simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precision cluster mass determination from weak lensing

Weak gravitational lensing has been used extensively in the past decade to constrain the masses of galaxy clusters, and is the most promising observational technique for providing the mass calibration necessary for precision cosmology with clusters. There are several challenges in estimating cluster masses, particularly (a) the sensitivity to astrophysical effects and observational systematics that modify the signal relative to the theoretical expectations, and (b) biases that can arise due to assumptions in the mass estimation method, such as the assumed radial profile of the cluster. All of these challenges are more problematic in the inner regions of the cluster, suggesting that their influence would ideally be suppressed for the purpose of mass estimation. However, at any given radius the differential surface density measured by lensing is sensitive to all mass within that radius, and the corrupted signal from the inner parts is spread out to all scales. We develop a new statistic that is ideal for estimation of cluster masses because it completely eliminates mass contributions below a chosen scale (which we suggest should be about 20 per cent of the virial radius), and thus reduces sensitivity to systematic and astrophysical effects. We use simulated and analytical profiles to quantify systematic biases on the estimated masses for several standard methods of mass estimation, finding that these can lead to significant mass biases that range from ten to over fifty per cent. The mass uncertainties when using our new statistic are reduced by up to a factor of ten relative to the standard methods, while only moderately increasing the statistical errors. This new method of mass estimation will enable a higher level of precision in future science work with weak lensing mass estimates for galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cosmological perturbation theory for baryons and dark matter I: one-loop corrections in the RPT framework

We generalize the `renormalized' perturbation theory (RPT) formalism of Crocce & Scoccimarro (2006a) to deal with multiple fluids in the Universe and here we present calculations up to the one-loop level. We apply the approach to the non-linear evolution of baryon and cold dark matter (CDM) perturbations, evolving from distinct sets of initial conditions. In current models of structure formation, it is standard practice to treat baryons and CDM as an effective single component fluid. In this approximation, one uses a weighed sum of late-time baryon and CDM transfer functions to set initial conditions. Here, we explore whether this approach can be used for high precision work. We show that, even for a pure linear treatment, there is a large-scale scale-dependent bias between baryons and CDM for WMAP5 cosmology. This bias is >1% until the present day, when it is driven towards unity through gravitational relaxation. Using the RPT formalism we test this approximation in the non-linear regime, and show that the CDM power spectrum in the 2-component fluid differs from that obtained from a 1-component fluid by ~3% on scales of order k~0.05 h/Mpc at z=10, and by ~0.5% at z=0. However, for the case of baryons the situation is worse and we find that the power spectrum is suppressed by ~15% on scales k~0.05 h/Mpc at z=10, and by ~3-5% at z=0, relative to the total matter. Importantly, besides the suppression of the spectrum, baryonic acoustic oscillations (BAO) are amplified for baryon and damped for CDM spectra. Thus, high precision modeling of baryons can not be probed through CDM only simulations; detection significance of BAO will be amplified in probes that study baryonic matter. Total mass can be modeled accurately using a 1-component fluid approach at all times. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Covariance of cross-correlations: towards efficient measures for large-scale structure

We study the covariance of the cross-power spectrum of different tracers for the large-scale structure. We develop the counts-in-cells framework for the multi-tracer approach, and use this to derive expressions for the full non-Gaussian covariance matrix. We show, that for the usual auto-power statistic, besides the off-diagonal covariance generated through gravitational mode-coupling, the discreteness of the tracers and their associated sampling distribution can generate strong off-diagonal covariance, and that this becomes the dominant source of covariance as k>>k_f=2 pi/L. On comparison with the derived expressions for the cross-power covariance, we show that the off-diagonal terms can be suppressed, if one cross-correlates a high tracer-density sample with a low one. Taking the effective estimator efficiency to be proportional to the signal-to-noise ratio (SN), we show that, to probe clustering as a function of physical properties of the sample, i.e. cluster mass or galaxy luminosity, then the cross-power approach can out perform the auto-power one by factors of a few. We confront the theory with measurements of the mass-mass, halo-mass, and halo-halo power spectra from a large ensemble of N-body simulations. We show that there is a significant SN advantage to be gained from using the cross-power approach when studying the bias of rare haloes. The analysis is repeated in configuration space and again SN improvement is found. We estimate the covariance matrix for these samples, and find strong off-diagonal contributions. The covariance depends on halo mass, with higher mass samples having stronger covariance. In agreement with theory, we show that the covariance is suppressed for the cross-power. This work points the way towards improved estimators for clustering studies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The cosmology dependence of weak lensing cluster counts

We present the main results of a numerical study of weak lensing cluster counting. We examine the scaling with cosmology of the projected-density-peak mass function. Our main conclusion is that the projected-peak and the three-dimensional mass functions scale with cosmology in an astonishingly close way. This means that, despite being derived from a two-dimensional field, the weak lensing cluster abundance can be used to constrain cosmology in the same way as the three-dimensional mass function probed by other types of surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The impact of correlated projections on weak lensing cluster counts

Large-scale structure projections are an obstacle in converting the shear signal of clusters detected in weak-lensing maps into virial masses. However, this step is not necessary for constraining cosmology with the shear-peak abundance, if we are able to predict its amplitude. We generate a large ensemble of N-body simulations spanning four cosmological models, with total volume V~1 (Gpc/h)^3 per model. Variations to the matter density parameter and amplitude of fluctuations are considered. We measure the abundance of peaks in the mass density projected in ~100 Mpc/h slabs to determine the impact of structures spatially correlated with the simulation clusters, identified by the 3D friends-of-friends algorithm. The halo model shows that the choice of the smoothing filter for the density field is important in reducing the contribution of correlated projections to individual halo masses. Such contributions are less than 2% in the case of the optimal, compensated filter used throughout this analysis. We measure the change in the mass of peaks when projected in slabs of various thicknesses. Peaks in slabs of 26 Mpc/h and 102 Mpc/h suffer an average mass change of less than 2% compared to their mass in slabs of 51 Mpc/h. We then explore the cosmology dependence of the projected-peak mass function, and find that, for a wide range of slab thicknesses (<500 Mpc/h), it scales with cosmology in exactly the same way as the 3D friends-of-friends mass function and the Sheth-Tormen formula. This extends the earlier result of Marian et al. (2009). Finally, we show that for all cosmological models considered, the low and intermediate mass bins of the peak abundance can be described using a modified Sheth-Tormen functional form to within 10%-20% accuracy.

0 citations0 reviewsNo code linkedOpen access