Galaxy Density Field Research Articles

Recovering dark-matter clustering from galaxies with Gaussianization

The Gaussianization transform has been proposed as a method to remove the issues of scale-dependent galaxy bias and nonlinearity from galaxy clustering statistics, but these benefits have yet to be thoroughly tested for realistic galaxy samples. In this paper, we test the effectiveness of the Gaussianization transform for different galaxy types by applying it to realistic simulated blue and red galaxy samples. We show that in real space, the shapes of the Gaussianized power spectra of both red and blue galaxies agree with that of the underlying dark matter, with the initial power spectrum, and with each other to smaller scales than do the statistics of the usual (untransformed) density field. However, we find that the agreement in the Gaussianized statistics breaks down in redshift space. We attribute this to the fact that red and blue galaxies exhibit very different fingers of god in redshift space. After applying a finger-of-god compression, the agreement on small scales between the Gaussianized power spectra is restored. We also compare the Gaussianization transform to the clipped galaxy density field and find that while both methods are effective in real space, they have more complicated behaviour in redshift space. Overall, we find that Gaussianization can be useful in recovering the shape of the underlying dark matter power spectrum to k ~ 0.5 h/Mpc and of the initial power spectrum to k ~ 0.4 h/Mpc in certain cases at z = 0.

Reconstructing the galaxy density field with photometric redshifts

Context. Measuring environment for large numbers of galaxies in the distant Universe is an open problem in astrophysics, as environment is important in determining many properties of galaxies during their formation and evolution. In order to measure galaxy environments, we need galaxy positions and redshifts. Photometric redshifts are more easily available for large numbers of galaxies, but at the price of larger uncertainties than spectroscopic redshifts.

Voids in massive neutrino cosmologies

Cosmic voids are a promising environment to characterize neutrino-induced effects on the large-scale distribution of matter in the universe. We perform a comprehensive numerical study of the statistical properties of voids, identified both in the matter and galaxy distributions, in massive and massless neutrino cosmologies. The matter density field is obtained by running several independent N-body simulations with cold dark matter and neutrino particles, while the galaxy catalogs are modeled by populating the dark matter halos in simulations via a halo occupation distribution (HOD) model to reproduce the clustering properties observed by the Sloan Digital Sky Survey (SDSS) II Data Release 7. We focus on the impact of massive neutrinos on the following void statistical properties: number density, ellipticities, two-point statistics, density and velocity profiles. Considering the matter density field, we find that voids in massive neutrino cosmologies are less evolved than those in the corresponding massless neutrinos case: there is a larger number of small voids and a smaller number of large ones, their profiles are less evacuated, and they present a lower wall at the edge. Moreover, the degeneracy between σ8 and Ων is broken when looking at void properties. In terms of the galaxy density field, we find that differences among cosmologies are difficult to detect because of the small number of galaxy voids in the simulations. Differences are instead present when looking at the matter density and velocity profiles around these voids.

The VIMOS Public Extragalactic Redshift Survey

Aims. Using the VIMOS Public Extragalactic Redshift Survey (VIPERS) we aim to jointly estimate the key parameters that describe the galaxy density field and its spatial correlations in redshift space. Methods. We use the Bayesian formalism to jointly reconstruct the redshift-space galaxy density field, power spectrum, galaxy bias and galaxy luminosity function given the observations and survey selection function. The high-dimensional posterior distribution is explored using the Wiener filter within a Gibbs sampler. We validate the analysis using simulated catalogues and apply it to VIPERS data taking into consideration the inhomogeneous selection function. Results. We present joint constraints on the anisotropic power spectrum as well as the bias and number density of red and blue galaxy classes in luminosity and redshift bins as well as the measurement covariances of these quantities. We find that the inferred galaxy bias and number density parameters are strongly correlated although these are only weakly correlated with the galaxy power spectrum. The power spectrum and redshift-space distortion parameters are in agreement with previous VIPERS results with the value of the growth rate $f\sigma_8 = 0.38$ with 18% uncertainty at redshift 0.7.

ECO AND RESOLVE: GALAXY DISK GROWTH IN ENVIRONMENTAL CONTEXT

We study the relationships between galaxy environments and galaxy properties related to disk (re)growth, considering two highly complete samples that are approximately baryonic mass limited into the high-mass dwarf galaxy regime, the Environmental COntext (ECO) catalog (data release herein) and the B-semester region of the REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. We quantify galaxy environments using both group identification and smoothed galaxy density field methods. We use by-eye and quantitative morphological classifications plus atomic gas content measurements and estimates. We find that blue early-type (E/S0) galaxies, gas-dominated galaxies, and UV-bright disk host galaxies all become distinctly more common below group halo mass ~10^11.5 Msun, implying that this low group halo mass regime may be a preferred regime for significant disk growth activity. We also find that blue early-type and blue late-type galaxies inhabit environments of similar group halo mass at fixed baryonic mass, consistent with a scenario in which blue early types can regrow late-type disks. In fact, we find that the only significant difference in the typical group halo mass inhabited by different galaxy classes is for satellite galaxies with different colors, where at fixed baryonic mass red early and late types have higher typical group halo masses than blue early and late types. More generally, we argue that the traditional morphology-environment relation (i.e., that denser environments tend to have more early types) can be largely attributed to the morphology-galaxy mass relation for centrals and the color-environment relation for satellites.

On fast generation of cosmological random fields

The statistical translation invariance of cosmological random fields is broken by a finite survey boundary, correlating the observable Fourier modes. Standard methods for generating Gaussian fields either neglect these correlations, or are costly, or both. Here we report on a fast and exact simulation method applicable to a wide class of two-point statistics that requires the simulation of a periodic grid of only twice the survey side with fast Fourier transforms. Super-survey modes, dominating the covariance of power spectra beyond linear scales in galaxy surveys and causing the correlation of large and small scales, "beat coupling", or "super-sample" covariance, are precisely accounted for in non-linear transformations of the Gaussian field. As an application, we simulate the CFHTLS $\sim 7^\circ \times 7^\circ$ W1 galaxy density field, modeled as a Poisson sampling of a lognormal density field. We show that our simulations produce power spectra, $A^*$-power spectra, counts-in-cells probability distributions as well as covariances perfectly consistent with the data. In addition, our technique reproduces the information plateau beyond linear scales as observed previously in Sloan Digital Sky Survey galaxy catalogs and in $N$-body simulations. Our method is thus an efficient yet powerful simulation and prediction tool for galaxy survey data and covariances.

Measuring the cosmic bulk flow with 6dFGSv

AbstractWhile recent years have seen rapid growth in the number of galaxy peculiar velocity measurements, disagreements remain about the extent to which the peculiar velocity field - a tracer of the large-scale distribution of mass - agrees with both ΛCDM expectations and with velocity field models derived from redshift surveys. The 6dF Galaxy Survey includes peculiar velocities for nearly 9 000 early-type galaxies (6dFGSv), making it the largest and most homogeneous galaxy peculiar velocity sample to date. We have used the 6dFGS velocity field to determine the amplitude and scale of large-scale cosmic flows in the local universe and test standard cosmological models. We also compare the galaxy density and peculiar velocity fields to establish the distribution of dark and luminous matter and better constrain key cosmological parameters such as the redshift-space distortion parameter.

Non-Gaussianity of the cosmic infrared background anisotropies – I. Diagrammatic formalism and application to the angular bispectrum

We present the first halo model based description of the Cosmic Infrared Background (CIB) non-Gaussianity (NG) that is fully parametric. To this end, we introduce, for the first time, a diagrammatic method to compute high order polyspectra of the 3D galaxy density field. It allows an easy derivation and visualisation of the different terms of the polyspectrum. We apply this framework to the power spectrum and bispectrum, and we show how to project them on the celestial sphere in the purpose of the application to the CIB angular anisotropies. Furthermore, we show how to take into account the particular case of the shot noise terms in that framework. Eventually, we compute the CIB angular bispectrum at 857 GHz and study its scale and configuration dependencies, as well as its variations with the halo occupation distribution parameters. Compared to a previously proposed empirical prescription, such physically motivated model is required to describe fully the CIB anisotropies bispectrum. Finally, we compare the CIB bispectrum with the bispectra of other signals potentially present at microwave frequencies, which hints that detection of CIB NG should be possible above 220 GHz.

More than just halo mass: modelling how the red galaxy fraction depends on multiscale density in an HOD framework

The fraction of galaxies with red colours depends sensitively on environment, and on the way in which environment is measured. To distinguish competing theories for the quenching of star formation, a robust and complete description of environment is required, to be applied to a large sample of galaxies. The environment of galaxies can be described using the density field of neighbours on multiple scales - the multiscale density field. We are using the Millennium simulation and a simple HOD prescription which describes the multiscale density field of Sloan Digital Sky Survey DR7 galaxies to investigate the dependence of the fraction of red galaxies on the environment. Using a volume limited sample where we have sufficient galaxies in narrow density bins, we have more dynamic range in halo mass and density for satellite galaxies than for central galaxies. Therefore we model the red fraction of central galaxies as a constant while we use a functional form to describe the red fraction of satellites as a function of halo mass which allows us to distinguish a sharp from a gradual transition. While it is clear that the data can only be explained by a gradual transition, an analysis of the multiscale density field on different scales suggests that colour segregation within the haloes is needed to explain the results. We also rule out a sharp transition for central galaxies, within the halo mass range sampled.

Extragalactic filament detection with a layer smoothing method

Filaments are clearly visible in galaxy distributions, but they are difficult to detect by computer algorithms. Most methods of filament detection can be used only with numerical simulations of a large-scale structure. New simple and effective methods for the real filament detection should be developed. The method of a smoothed galaxy density field was applied in this work to SDSS data of galaxy positions. Five concentric radial layers of 100 Mpc are appropriate for filaments detection. Two methods were tested for the first layer and one more method is proposed.

Weak lensing galaxy cluster field reconstruction

In this paper, we compare three methods to reconstruct galaxy cluster density fields with weak lensing data. The first method called FLens integrates an inpainting concept to invert the shear field with possible gaps, and a multi-scale entropy denoising procedure to remove the noise contained in the final reconstruction, that arises mostly from the random intrinsic shape of the galaxies. The second and third methods are based on a model of the density field made of a multi-scale grid of radial basis functions. In one case, the model parameters are computed with a linear inversion involving a singular value decomposition. In the other case, the model parameters are estimated using a Bayesian MCMC optimization implemented in the lensing software Lenstool. Methods are compared on simulated data with varying galaxy density fields. We pay particular attention to the errors estimated with resampling. We find the multi-scale grid model optimized with MCMC to provide the best results, but at high computational cost, especially when considering resampling. The SVD method is much faster but yields noisy maps, although this can be mitigated with resampling. The FLens method is a good compromise with fast computation, high signal to noise reconstruction, but lower resolution maps. All three methods are applied to the MACS J0717+3745 galaxy cluster field, and reveal the filamentary structure discovered in Jauzac et al. 2012. We conclude that sensitive priors can help to get high signal to noise, and unbiased reconstructions.

Cosmological information in the intrinsic alignments of luminous red galaxies

The intrinsic alignments of galaxies are usually regarded as a contaminant to weak gravitational lensing observables. The alignment of Luminous Red Galaxies, detected unambiguously in observations from the Sloan Digital Sky Survey, can be reproduced by the linear tidal alignment model of Catelan, Kamionkowski & Blandford (2001) on large scales. In this work, we explore the cosmological information encoded in the intrinsic alignments of red galaxies. We make forecasts for the ability of current and future spectroscopic surveys to constrain local primordial non-Gaussianity and Baryon Acoustic Oscillations (BAO) in the cross-correlation function of intrinsic alignments and the galaxy density field. For the Baryon Oscillation Spectroscopic Survey, we find that the BAO signal in the intrinsic alignments is marginally significant with a signal-to-noise ratio of 1.8 and 2.2 with the current LOWZ and CMASS samples of galaxies, respectively, and increasing to 2.3 and 2.7 once the survey is completed. For the Dark Energy Spectroscopic Instrument and for a spectroscopic survey following the EUCLID redshift selection function, we find signal-to-noise ratios of 12 and 15, respectively. Local type primordial non-Gaussianity, parametrized by fNL = 10, is only marginally significant in the intrinsic alignments signal with signal-to-noise ratios < 2 for the three surveys considered.

Independent constraints on local non-Gaussianity from the peculiar velocity and density fields

Primordial, non-Gaussian perturbations can generate scale-dependent bias in the galaxy distribution. This in turn will modify correlations between galaxy positions and peculiar velocities at late times, since peculiar velocities reflect the underlying matter distribution, whereas galaxies are a biased tracer of the same. We study this effect, and show that non-Gaussianity can be constrained by comparing the observed peculiar velocity field to a model velocity field reconstructed from the galaxy density field assuming linear bias. The amplitude of the spatial correlations in the residual map obtained after subtracting one velocity field from the other is directly proportional to the strength of the primordial non-Gaussianity. We construct the corresponding likelihood function use it to constrain the amplitude of the linear flow $\beta$ and the amplitude of local non-Gaussianity $f^{\rm NL}_{\rm local}$. Applying our method to two observational data sets, the Type-Ia supernovae (A1SN) and Spiral Field \textit{I}-band (SFI++) catalogues, we obtain constraints on the linear flow parameter consistent with the values derived previously assuming Gaussianity. The marginalised 1-D distribution of $|f^{\rm NL}_{\rm local}|$ does not show strong evidence for non-zero $f^{\rm NL}_{\rm local}$, and we set 95% upper limits $|f^{\rm NL}_{\rm local}|<51.4$ from A1SN and $|f^{\rm NL}_{\rm local}|<92.6$ from SFI++. These limits on $f^{\rm NL}_{\rm local}$ are as tight as any set by previous large-scale structure measurements. Our method can be applied to any survey with radial velocities and density field data, and provides an independent check of recent CMB constraints on $f^{\rm NL}_{\rm local}$, extending these to smaller spatial scales.

Simulating structure formation of the Local Universe

In this work we present cosmological N-body simulations of the Local Universe with initial conditions constrained by the Two-Micron Redshift Survey (2MRS) within a cubic volume of 180 Mpc/h side-length centred at the Local Group. We use a self-consistent Bayesian based approach to explore the joint parameter space of primordial density fluctuations and peculiar velocity fields, which are compatible with the 2MRS galaxy distribution after cosmic evolution. This method (the KIGEN-code) includes the novel ALPT (Augmented Lagrangian Perturbation Theory) structure formation model which combines second order LPT (2LPT) on large scales with the spherical collapse model on small scales. Furthermore we describe coherent flows with 2LPT and include a dispersion term to model fingers-of-god (fogs). These implementations are crucial to avoid artificial filamentary structures, which appear when using a structure formation model with 2LPT and data with compressed fogs. We assume LCDM cosmology throughout our method. The recovered initial Gaussian fields are used to perform a set of 25 constrained simulations. Statistically this ensemble of simulations is in agreement with a reference set of 25 simulations based on randomly seeded Gaussian fluctuations in terms of matter statistics, power-spectra and mass functions. Considering the entire volume, we obtain correlation coefficients of about 98.3% for the comparison between the simulated density fields and the galaxy density field in log-space with Gaussian smoothing scales of r_S=3.5 Mpc/h (74% for r_S=1.4 Mpc/h). The cross power-spectra show correlations with the galaxy distribution which weaken towards smaller length scales until they vanish at scales of 2.2-3.0 Mpc/h. The simulations we present provide a fully nonlinear density and velocity field with a high level of correlation with the observed galaxy distribution at scales of a few Mpc.

Relativistic effects in galaxy clustering in a parametrized post-Friedmann universe

We explore the signatures of quintessence and modified gravity theories in the relativistic description of galaxy clustering within a parametrized post-Friedmann framework. For this purpose, we develop a calibration method to consistently account for horizon-scale effects in the linear parametrized Post-Friedmann perturbations of minimally and nonminimally coupled scalar-tensor theories and test it against the full model-specific fluctuations. We further study the relativistic effects in galaxy clustering for the normal and self-accelerating branches of the Dvali-Gabadadze-Porrati braneworld model as well as for phenomenological modifications of gravity. We quantify the impact of modified gravity and dark energy models on galaxy clustering by computing the velocity-to-matter density ratio F, the velocity contribution R, and the potential contribution P and give an estimate of their detectability in future galaxy surveys. Our results show that, in general, the relativistic correction contains additional information on gravity and dark energy, which needs to be taken into account in consistent horizon-scale tests of departures from LCDM using the galaxy-density field.

SUZAKUOBSERVATIONS OF THE OUTSKIRTS OF A1835: DEVIATION FROM HYDROSTATIC EQUILIBRIUM

We present results of four-pointing Suzaku X-ray observations (total ~200 ks) of the intracluster medium (ICM) in the Abell 1835 galaxy cluster (kT ~ 8 keV, z = 0.253) out to the virial radius (r_vir ~ 2.9 Mpc) and beyond. Faint X-ray emission from the ICM out to r_vir is detected. The temperature gradually decreases with radius from ~8 keV in the inner region to ~2 keV at r_vir. The entropy profile is shown to flatten beyond r_500, in disagreement with the r_1.1 dependence predicted from the accretion shock heating model. The thermal pressure profile in the range 0.3r_500 < r < r_vir agrees well with that obtained from the stacked Sunyaev-Zel'dovich effect observations with the Planck satellite. The hydrostatic mass profile in the cluster outskirts (r_500 < r < r_vir) falls well short of the weak lensing one derived from Subaru/Suprime-Cam observations, showing an unphysical decrease with radius. The gas mass fraction at r_vir defined with the lensing total mass agrees with the cosmic baryon fraction from the WMAP 7-year data. All these results indicate, rather than the gas-clumping effect, that the bulk of the ICM in the cluster outskirts is far from hydrostatic equilibrium and infalling matter retained some of its kinetic energy. Finally, combining with our recent Suzaku and lensing analysis of Abell 1689, a cluster of similar mass, temperature, and redshift, we show that the cluster temperature distribution in the outskirts is significantly correlated with the galaxy density field in the surrounding large-scale environment at (1-2)r_vir.

Cross-bispectra and trispectra of the non-linear integrated Sachs–Wolfe effect and the tracer galaxy density field

In order to investigate possibilities to measure non-Gaussian signatures of the non-linear iSW effect, we study in this work the family of mixed bispectra <tau^q gamma^(3-q)> and trispectra <tau^q gamma^(4-q)> between the integrated Sachs-Wolfe (iSW) temperature perturbation tau and the galaxy over-density gamma. We use standard Eulerian perturbation theory restricted to tree level expansion for predicting the cosmic matter field. As expected, the spectra are found to decrease in amplitude with increasing q. The transition scale between linear domination and the scales, on which non-linearities take over, moves to larger scales with increasing number of included iSW source fields q. We derive the cumulative signal-to-noise ratios for a combination of Planck CMB data and the galaxy sample of a Euclid-like survey. Including scales down to l_max = 1000 we find sobering values of sigma = 0.83 for the mixed bispectrum and sigma = 0.19 in case of the trispectrum for q=1. For higher values of q the polyspectra <tau^2 gamma> and <tau^3 gamma> are found to be far below the detection limit.

A method to search for correlations of UHECR masses with the large scale structures in the local galaxy density field

One of the main goals of investigations using present and future giant extensive air shower (EAS) arrays is the mass composition of ultra-high energy cosmic rays (UHECRs). A new approach to the problem is presented, combining analysis of arrival directions with the statistical test of the pairs of EAS samples. An idea of the method is to search for possible correlations of UHECR masses with their presumably separate sources. The method is based on a non-parametric statistical test, specifically Wilcoxon rank sum routine, which does not depend on the populations fitting any parameterized distributions.

The clustering of galaxies in the SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: constraints on primordial non-Gaussianity

We analyze the density field of 264,283 galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) and included in the SDSS data release nine (DR9). In total, the SDSS DR9 BOSS data includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of more than 3,000 deg^2. We measure the power spectrum of these galaxies with redshifts 0.43 < z < 0.7 in order to constrain the amount of local non-Gaussianity, f_NL,local, in the primordial density field, paying particular attention to the impact of systematic uncertainties. The BOSS galaxy density field is systematically affected by the local stellar density and this influences the ability to accurately measure f_NL,local. In the absence of any correction, we find (erroneously) that the probability that f_NL,local is greater than zero, P(f_NL,local >0), is 99.5%. After quantifying and correcting for the systematic bias and including the added uncertainty, we find -45 < f_NL,local < 195 at 95% confidence, and P(f_NL,local >0) = 91.0%. A more conservative approach assumes that we have only learned the k-dependence of the systematic bias and allows any amplitude for the systematic correction; we find that the systematic effect is not fully degenerate with that of f_NL,local, and we determine that -82 < f_NL,local < 178 (at 95% confidence) and P(f_NL,local >0) = 68%. This analysis demonstrates the importance of accounting for the impact of Galactic foregrounds on f_NL,local measurements. We outline the methods that account for these systematic biases and uncertainties. We expect our methods to yield robust constraints on f_NL,local for both our own and future large-scale-structure investigations.

Second-order matter fluctuations via higher order galaxy correlators

We provide a formula for extracting the value of the rms of the linear matter fluctuations on a scale R directly from redshift surveys data. It allows to constrain the real-space amplitude of sigma_R without requiring any modeling of the nature and power spectrum of the matter distribution. Furthermore, the formalism is completely insensitive to the character of the bias function, namely its eventual scale or non-linear dependence. By contrasting measurements of sigma_R with predictions from linear perturbation theory, one can test for eventual departures from the standard description of gravity on large cosmological scales. The proposed estimator exploits the information contained in the 1-point moments and 2-point correlators of the matter and galaxy density fields, and it can be applied on cosmic scales where linear and semi-linear perturbative approximations of the evolution of matter overdensities offer a satisfactory description of the full underlying theory. We implement the test with N-body simulations to quantify potential systematics and successfully show that we are able to recover the present day value of sigma_8 `hidden' in the simulation. We also design a consistency check to gauge the soundness of the results inferred when the formalism is applied to real (as opposed to simulated) data. We expect that this approach will provide a sensitive probe of the clustering of matter when applied to future large redshift survey such as BigBOSS and EUCLID.