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

Francisco-Shu Kitaura

Francisco-Shu Kitaura 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.GAApplicationsastro-ph.IMmath.STStatistics TheoryArtificial IntelligenceMethodology

Catalog footprint

What is connected

43works
8topics
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

43 published item(s)

preprint2022arXiv

Cosmic Void Baryon Acoustic Oscillation Measurement: Evaluation of Sensitivity to Selection Effects

Cosmic voids defined as a subset of Delaunay Triangulation (DT) circumspheres have been used to measure the Baryon Acoustic Oscillations (BAO) scale; providing tighter constraints on cosmological parameters when combined with matter tracers. These voids are defined as spheres larger than a given radius threshold, which is constant over the survey volume. However, the response of these void tracers to observational systematics has not yet been studied. In this work we analyse the response of void clustering to selection effects. We find for the case of moderate (<20 per cent) incompleteness, void selection based on a constant radius cut yields robust measurements. This is particularly true for BAO-reconstructed galaxy samples, where large-scale void exclusion effects are mitigated. Moreover, we observe for the case of severe (up to 90 per cent) incompleteness -- such as can be found at the edges of the radial selection function -- that an accurate estimation of the void distribution is necessary for unbiased clustering measurements. In addition, we find that without reconstruction, using a constant threshold under these conditions produces a stronger void exclusion effect that can affect the clustering on large scales. A new void selection criteria dependent on the (local) observed tracer density that maximises the BAO peak significance prevents the aforementioned exclusion features from contaminating the BAO signal. Finally, we verify, with large simulations including light cone evolution, that both void sample definitions (local and constant) yield unbiased and consistent BAO scale measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Predicted future fate of COSMOS galaxy protoclusters over 11 Gyr with constrained simulations

Cosmological simulations are crucial tools in studying the Universe, but they typically do not directly match real observed structures. Constrained cosmological simulations, on the other hand, are designed to match the observed distribution of galaxies. Here we present constrained simulations based on spectroscopic surveys at a redshift of z~2.3, corresponding to an epoch of nearly 11 Gyrs ago. This allows us to 'fast-forward' the simulation to our present-day and study the evolution of observed cosmic structures self-consistently. We confirm that several previously-reported protoclusters will evolve into massive galaxy clusters by our present epoch, including the 'Hyperion' structure that we predict will collapse into a giant filamentary supercluster spanning 100 Megaparsecs. We also discover previously unknown protoclusters, with lower final masses than typically detectable by other methods, that nearly double the number of known protoclusters within this volume. Constrained simulations, applied to future high-redshift datasets, represents a unique opportunity for studying early structure formation and matching galaxy properties between high and low redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from multi-tracer BAO analysis with galaxies and voids

We construct cosmic void catalogues with the DIVE void finder upon SDSS BOSS DR12 and eBOSS DR16 galaxy samples with BAO reconstruction applied, and perform a joint BAO analysis using different types of galaxies and the corresponding voids. The BAO peak is evident for the galaxy-galaxy, galaxy-void, and void-void correlation functions of all datasets, including the ones cross correlating luminous red galaxy and emission line galaxy samples. Two multi-tracer BAO fitting schemes are then tested, one combining the galaxy and void correlation functions with a weight applied to voids, and the other using a single BAO dilation parameter for all clustering measurements of different tracers. Both methods produce consistent results with mock catalogues, and on average ~10 per cent improvements of the BAO statistical uncertainties are observed for all samples, compared to the results from galaxies alone. By combining the clustering of galaxies and voids, the uncertainties of BAO measurements from the SDSS data are reduced by 5 to 15 per cent, yielding 0.9, 0.8, 1.1, 2.3, and 2.9 per cent constraints on the distance $D_{_{\rm V}}(z)$, at effective redshifts of 0.38, 0.51, 0.70, 0.77, and 0.85, respectively. When combined with BAO measurements from SDSS MGS, QSO, and Ly$α$ samples, as well as the BBN results, we obtain $H_0 = 67.58 \pm 0.91\,{\rm km}\,{\rm s}^{-1}\,{\rm Mpc}^{-1}$, $Ω_{\rm m} = 0.290 \pm 0.015$, and $Ω_Λh^2 = 0.3241 \pm 0.0079$ in the flat-$Λ$CDM framework, where the 1$\,σ$ uncertainties are around 6, 6, and 17 per cent smaller respectively, compared to constraints from the corresponding anisotropic BAO measurements without voids and LRG-ELG cross correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

UNITSIM-Galaxies: data release and clustering of emission-line galaxies

New surveys such as ESA's Euclid mission are planned to map with unprecedented precision the large-scale structure of the Universe by measuring the 3D positions of tens of millions of galaxies. It is necessary to develop theoretically modelled galaxy catalogues to estimate the expected performance and to optimise the analysis strategy of these surveys. We populate two pairs of (1 Gpc/h)^3 volume dark-matter-only simulations from the UNIT project with galaxies using the SAGE semi-analytic model of galaxy formation, coupled to the photoionisation model get_emlines to estimate their H_alpha emission. These catalogues represent a unique suite that includes galaxy formation physics and - thanks to the fixed-pair technique used - an effective volume of ~(5 Gpc/h)^3, which is several times larger than the Euclid survey. We present the performance of these data and create five additional emission-line galaxy (ELG) catalogues by applying a dust attenuation model as well as adjusting the flux threshold as a function of redshift in order to reproduce Euclid-forecast dN/dz values. As a first application, we study the abundance and clustering of those model H_alpha ELGs: for scales greater than ~5 Mpc/h, we find a scale-independent bias with a value of $b\sim 1$ at redshift z ~ 0.5, that can increase nearly linearly to $b\sim 4$ at $z\sim 2$, depending on the ELG catalogue. Model galaxy properties, including their emission-line fluxes (with and without dust extinction) are publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Constraining Primordial Non-Gaussianity with Post-reconstructed Galaxy Bispectrum in Redshift Space

Galaxy bispectrum is a promising probe of inflationary physics in the early universe as a measure of primordial non-Gaussianity (PNG), whereas its signal-to-noise ratio is significantly affected by the mode coupling due to non-linear gravitational growth. In this paper, we examine the standard reconstruction method of linear cosmic mass density fields from non-linear galaxy density fields to de-correlate the covariance in redshift-space galaxy bispectra. In particular, we evaluate the covariance of the bispectrum for massive-galaxy-sized dark matter halos with reconstruction by using 4000 independent $N$-body simulations. Our results show that the bispectrum covariance for the post-reconstructed field approaches the Gaussian prediction at scale of $k<0.2\, h\, {\rm Mpc}^{-1}$. We also verify the leading-order PNG-induced bispectrum is not affected by details of the reconstruction with perturbative theory. We then demonstrate the constraining power of the post-reconstructed bispectrum for PNG at redshift of $\sim0.5$. Further, we perform a Fisher analysis to make a forecast of PNG constraints by galaxy bispectra including anisotropic signals. Assuming a massive galaxy sample in the SDSS Baryon Oscillation Spectroscopic Survey, we find that the post-reconstructed bispectrum can constrain the local-, equilateral- and orthogonal-types of PNG with $Δf_{\rm NL} \sim$13, 90 and 42, respectively, improving the constraints with the pre-reconstructed bispectrum by a factor of $1.3-3.2$. In conclusion, the reconstruction plays an essential role in constraining various types of PNG signatures with a level of $Δf_{\rm NL}<1$ from the galaxy bispectrum based on upcoming galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Improving baryon acoustic oscillation measurement with the combination of cosmic voids and galaxies

We develop a methodology to optimise the measurement of Baryon Acoustic Oscillation (BAO) from a given galaxy sample. In our previous work, we demonstrated that one can measure BAO from tracers in under-dense regions (voids). In this study, we combine the over-dense and under-dense tracers (galaxies & voids) to obtain better constraints on the BAO scale. To this end, we modify the de-wiggled BAO model with an additional parameter to describe both the BAO peak and the underlying exclusion pattern of void 2PCFs. We show that after applying BAO reconstruction to galaxies, the BAO peak scale of both galaxies and voids are unbiased using the modified model. Furthermore, we use a new 2PCF estimator for a multi-tracer analysis with galaxies and voids. In simulations, the joint sample improves by about 10% the constraint for the post-reconstruction BAO peak position compared to the result from galaxies alone, which is equivalent to an enlargement of the survey volume by 20%. Applying this method to the BOSS DR12 data, we have an 18% improvement for the low-z sample (0.2<z<0.5), but a worse constraint for the high-z sample (0.5<z<0.75), which is consistent with statistical fluctuations for the current survey volume. Future larger samples will give more robust improvements due to less statistical fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The bias of dark matter tracers: assessing the accuracy of mapping techniques

We present a comparison between approximated methods for the construction of mock catalogs based on the halo-bias mapping technique. To this end, we use as reference a high resolution $N$-body simulation of 3840$^3$ dark matter particles on a 400$h^{-1}\rm{Mpc}$ cube box from the Multidark suite. In particular, we explore parametric versus non-parametric bias mapping approaches and compare them at reproducing the halo distribution in terms of the two and three point statistics down to $\sim 10^8\,{\rm M}_{\odot}\,h^{-1}$ halo masses. Our findings demonstrate that the parametric approach remains inaccurate even including complex deterministic and stochastic components. On the contrary, the non-parametric one is indistinguishable from the reference $N$-body calculation in the power-spectrum beyond $k=1\,h\,{\rm Mpc}^{-1}$, and in the bispectrum for typical configurations relevant to baryon acoustic oscillation analysis. We conclude, that approaches which extract the full bias information from $N$-body simulations in a non-parametric fashion are ready for the analysis of the new generation of large scale structure surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Signatures of the Primordial Universe from Its Emptiness: Measurement of Baryon Acoustic Oscillations from Minima of the Density Field

Sound waves from the primordial fluctuations of the Universe imprinted in the large-scale structure, called baryon acoustic oscillations (BAOs), can be used as standard rulers to measure the scale of the Universe. These oscillations have already been detected in the distribution of galaxies. Here we propose to measure BAOs from the troughs (minima) of the density field. Based on two sets of accurate mock halo catalogues with and without BAOs in the seed initial conditions, we demonstrate that the BAO signal cannot be obtained from the clustering of classical disjoint voids, but is clearly detected from overlapping voids. The latter represent an estimate of all troughs of the density field. We compute them from the empty circumsphere centers constrained by tetrahedra of galaxies using Delaunay triangulation. Our theoretical models based on an unprecedented large set of detailed simulated void catalogues are remarkably well confirmed by observational data. We use the largest recently publicly available sample of luminous red galaxies from SDSS-III BOSS DR11 to unveil for the first time a >3$σ$ BAO detection from voids in observations. Since voids are nearly isotropically expanding regions, their centers represent the most quiet places in the Universe, keeping in mind the cosmos origin and providing a new promising window in the analysis of the cosmological large-scale structure from galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in Fourier-space

We analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is performed in Fourier-space, using the power spectrum monopole and quadrupole. The dataset includes $1\,198\,006$ galaxies over the redshift range $0.2 < z < 0.75$. We divide this dataset into three (overlapping) redshift bins with the effective redshifts $\zeff = 0.38$, $0.51$ and $0.61$. We demonstrate the reliability of our analysis pipeline using N-body simulations as well as $\sim 1000$ MultiDark-Patchy mock catalogues, which mimic the BOSS-DR12 target selection. We apply density field reconstruction to enhance the BAO signal-to-noise ratio. By including the power spectrum quadrupole we can separate the line-of-sight and angular modes, which allows us to constrain the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ separately. We obtain two independent $1.6\%$ and $1.5\%$ constraints on $D_A(z)$ and $2.9\%$ and $2.3\%$ constraints on $H(z)$ for the low ($\zeff=0.38$) and high ($\zeff=0.61$) redshift bin, respectively. We obtain two independent $1\%$ and $0.9\%$ constraints on the angular averaged distance $D_V(z)$, when ignoring the Alcock-Paczynski effect. The detection significance of the BAO signal is of the order of $8σ$ (post-reconstruction) for each of the three redshift bins. Our results are in good agreement with the Planck prediction within $Λ$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in~\citet{Alam2016} to produce the final cosmological constraints from BOSS.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: combining correlated Gaussian posterior distributions

The cosmological information contained in anisotropic galaxy clustering measurements can often be compressed into a small number of parameters whose posterior distribution is well described by a Gaussian. We present a general methodology to combine these estimates into a single set of consensus constraints that encode the total information of the individual measurements, taking into account the full covariance between the different methods. We illustrate this technique by applying it to combine the results obtained from different clustering analyses, including measurements of the signature of baryon acoustic oscillations (BAO) and redshift-space distortions (RSD), based on a set of mock catalogues of the final SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). Our results show that the region of the parameter space allowed by the consensus constraints is smaller than that of the individual methods, highlighting the importance of performing multiple analyses on galaxy surveys even when the measurements are highly correlated. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The methodology presented here is used in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the configuration-space clustering wedges

We explore the cosmological implications of anisotropic clustering measurements in configuration space of the final galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We implement a new detailed modelling of the effects of non-linearities, galaxy bias and redshift-space distortions that can be used to extract unbiased cosmological information from our measurements for scales $s \gtrsim 20\,h^{-1}{\rm Mpc}$. We combined the galaxy clustering information from BOSS with the latest cosmic microwave background (CMB) observations and Type Ia supernovae samples and found no significant evidence for a deviation from the $Λ$CDM cosmological model. In particular, these data sets can constrain the dark energy equation of state parameter to $w_{\rm DE}=-0.996\pm0.042$ when assumed time-independent, the curvature of the Universe to $Ω_{k}=-0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum m_ν < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the growth rate of cosmic structures assuming $f(z)=Ω_{\rm m}(z)^γ$ and obtain $γ= 0.609\pm 0.079$, in good agreement with the predictions of general relativity of $γ=0.55$. We compress the information of our clustering measurements into constraints on the parameter combinations $D_{\rm V}(z)/r_{\rm d}$, $F_{\rm AP}(z)$ and $fσ_8(z)$ at the effective redshifts of $z=0.38$, $0.51$ and $0.61$ with their respective covariance matrices and find good agreement with the predictions for these parameters obtained from the best-fitting $Λ$CDM model to the CMB data from the Planck satellite. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: double-probe measurements from BOSS galaxy clustering \& Planck data -- towards an analysis without informative priors

We develop a new methodology called double-probe analysis with the aim of minimizing informative priors in the estimation of cosmological parameters. We extract the dark-energy-model-independent cosmological constraints from the joint data sets of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy sample and Planck cosmic microwave background (CMB) measurement. We measure the mean values and covariance matrix of $\{R$, $l_a$, $Ω_b h^2$, $n_s$, $log(A_s)$, $Ω_k$, $H(z)$, $D_A(z)$, $f(z)σ_8(z)\}$, which give an efficient summary of Planck data and 2-point statistics from BOSS galaxy sample, where $R=\sqrt{Ω_m H_0^2}\,r(z_*)$, and $l_a=πr(z_*)/r_s(z_*)$, $z_*$ is the redshift at the last scattering surface, and $r(z_*)$ and $r_s(z_*)$ denote our comoving distance to $z_*$ and sound horizon at $z_*$ respectively. The advantage of this method is that we do not need to put informative priors on the cosmological parameters that galaxy clustering is not able to constrain well, i.e. $Ω_b h^2$ and $n_s$. Using our double-probe results, we obtain $Ω_m=0.304\pm0.009$, $H_0=68.2\pm0.7$, and $σ_8=0.806\pm0.014$ assuming $Λ$CDM; and $Ω_k=0.002\pm0.003$ and $w=-1.00\pm0.07$ assuming o$w$CDM. The results show no tension with the flat $Λ$CDM cosmological paradigm. By comparing with the full-likelihood analyses with fixed dark energy models, we demonstrate that the double-probe method provides robust cosmological parameter constraints which can be conveniently used to study dark energy models. We extend our study to measure the sum of neutrino mass and obtain $Σm_ν<0.10/0.22$ (68\%/95\%) assuming $Λ$CDM and $Σm_ν<0.26/0.52$ (68\%/95\%) assuming $w$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Observational systematics and baryon acoustic oscillations in the correlation function

We present baryon acoustic oscillation (BAO) scale measurements determined from the clustering of 1.2 million massive galaxies with redshifts 0.2 < z < 0.75 distributed over 9300 square degrees, as quantified by their redshift-space correlation function. In order to facilitate these measurements, we define, describe, and motivate the selection function for galaxies in the final data release (DR12) of the SDSS III Baryon Oscillation Spectroscopic Survey (BOSS). This includes the observational footprint, masks for image quality and Galactic extinction, and weights to account for density relationships intrinsic to the imaging and spectroscopic portions of the survey. We simulate the observed systematic trends in mock galaxy samples and demonstrate that they impart no bias on baryon acoustic oscillation (BAO) scale measurements and have a minor impact on the recovered statistical uncertainty. We measure transverse and radial BAO distance measurements in 0.2 < z < 0.5, 0.5 < z < 0.75, and (overlapping) 0.4 < z < 0.6 redshift bins. In each redshift bin, we obtain a precision that is 2.7 per cent or better on the radial distance and 1.6 per cent or better on the transverse distance. The combination of the redshift bins represents 1.8 per cent precision on the radial distance and 1.1 per cent precision on the transverse distance. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: BAO measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies

[abridged] We present an anisotropic analysis of the baryonic acoustic oscillation (BAO) scale in the twelfth and final data release of the Baryonic Oscillation Spectroscopic Survey (BOSS). We independently analyse the LOWZ and CMASS galaxy samples: the LOWZ sample contains contains 361 762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$; the CMASS sample consists of 777 202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. We extract the BAO peak position from the monopole power spectrum moment, $α_0$, and from the $μ^2$ moment, $α_2$, where $μ$ is the cosine of the angle to the line-of-sight. The $μ^2$-moment provides equivalent information to that available in the quadrupole but is simpler to analyse. After applying a reconstruction algorithm to reduce the BAO suppression by bulk motions, we measure the BAO peak position in the monopole and $μ^2$-moment, which are related to radial and angular shifts in scale. We report $H(z_{\rm LOWZ})r_s(z_d)=(11.60\pm0.60)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.66\pm0.16$ with a cross-correlation coefficient of $r_{HD_A}=0.41$, for the LOWZ sample; and $H(z_{\rm CMASS})r_s(z_d)=(14.56\pm0.37)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.42\pm0.13$ with a cross-correlation coefficient of $r_{HD_A}=0.47$, for the CMASS sample. We combine these results with the measurements of the BAO peak position in the monopole and quadrupole correlation function of the same dataset \citep[][companion paper]{Cuestaetal2015} and report the consensus values: $H(z_{\rm LOWZ})r_s(z_d)=(11.63\pm0.69)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.67\pm0.15$ with $r_{HD_A}=0.35$ for the LOWZ sample; $H(z_{\rm CMASS})r_s(z_d)=(14.67\pm0.42)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.47\pm0.12$ with $r_{HD_A}=0.52$ for the CMASS sample.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: mock galaxy catalogues for the BOSS Final Data Release

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $σ$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Modeling the clustering and halo occupation distribution of BOSS-CMASS galaxies in the Final Data Release

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $Λ$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$σ$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$σ$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies

We measure and analyse the clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) relative to the line-of-sight (LOS), for LOWZ and CMASS galaxy samples drawn from the final Data Release 12 (DR12). The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm lowz}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm cmass}=0.57$. From the power spectrum monopole and quadrupole moments around the LOS, we measure the growth of structure parameter $f$ times the amplitude of dark matter density fluctuations $σ_8$ by modeling the Redshift-Space Distortion signal. When the geometrical Alcock-Paczynski effect is also constrained from the same data, we find joint constraints on $fσ_8$, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon $D_A(z)/r_s(z_d)$. We find $f(z_{\rm lowz})σ_8(z_{\rm lowz})=0.394\pm0.062$, $D_A(z_{\rm lowz})/r_s(z_d)=6.35\pm0.19$, $H(z_{\rm lowz})r_s(z_d)=(11.41\pm 0.56)\,{10^3\rm km}s^{-1}$ for the LOWZ sample, and $f(z_{\rm cmass})σ_8(z_{\rm cmass})=0.444\pm0.038$, $D_A(z_{\rm cmass})/r_s(z_d)=9.42\pm0.15$, $H(z_{\rm cmass})r_s(z_d)=(13.92 \pm 0.44)\, {10^3\rm km}s^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 measurements. Assuming the Hubble parameter and angular distance parameter are fixed at fiducial $Λ$CDM values, we find $f(z_{\rm lowz})σ_8(z_{\rm lowz})=0.485\pm0.044$ and $f(z_{\rm cmass})σ_8(z_{\rm cmass})=0.436\pm0.022$ for the LOWZ and CMASS samples, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the power spectrum and bispectrum of the DR12 BOSS galaxies

We measure and analyse the bispectrum of the final, Data Release 12, galaxy sample provided by the Baryon Oscillation Spectroscopic Survey, splitting by selection algorithm into LOWZ and CMASS galaxies. The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. Combining the power spectrum, measured relative to the line-of-sight, with the spherically averaged bispectrum, we are able to constrain the product of the growth of structure parameter, $f$, and the amplitude of dark matter density fluctuations, $σ_8$, along with the geometric Alcock-Paczynski parameters, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch, $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon, $D_A(z)/r_s(z_d)$. After combining pre-reconstruction RSD analyses of the power spectrum monopole, quadrupole and bispectrum monopole; with post-reconstruction analysis of the BAO power spectrum monopole and quadrupole, we find $f(z_{\rm LOWZ})σ_8(z_{\rm LOWZ})=0.427\pm 0.056$, $D_A(z_{\rm LOWZ})/r_s(z_d)=6.60 \pm 0.13$, $H(z_{\rm LOWZ})r_s(z_d)=(11.55\pm 0.38)10^3\,{\rm kms}^{-1}$ for the LOWZ sample, and $f(z_{\rm CMASS})σ_8(z_{\rm CMASS})=0.426\pm 0.029$, $D_A(z_{\rm CMASS})/r_s(z_d)=9.39 \pm 0.10$, $H(z_{\rm CMASS})r_s(z_d)=(14.02\pm 0.22)10^3\,{\rm kms}^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 and DR12 measurements. Combining our dataset with {\it Planck15} we perform a null test of General Relativity (GR) through the $γ$-parametrisation finding $γ=0.733^{+0.068}_{-0.069}$, which is $\sim2.7σ$ away from the GR predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: single-probe measurements from CMASS anisotropic galaxy clustering

With the largest spectroscopic galaxy survey volume drawn from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), we can extract cosmological constraints from the measurements of redshift and geometric distortions at quasi-linear scales (e.g. above 50 $h^{-1}$Mpc). We analyze the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS galaxy sample, at the effective redshift $z=0.59$, to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)σ_8(z)$, and the physical matter density $Ω_mh^2$. We obtain robust measurements by including a polynomial as the model for the systematic errors, and find it works very well against the systematic effects, e.g., ones induced by stars and seeing. We provide accurate measurements $\{D_A(0.59)r_{s,fid}/r_s$ $\rm Mpc$, $H(0.59)r_s/r_{s,fid}$ $km s^{-1} Mpc^{-1}$, $f(0.59)σ_8(0.59)$, $Ω_m h^2\}$ = $\{1427\pm26$, $97.3\pm3.3$, $0.488 \pm 0.060$, $0.135\pm0.016\}$, where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$ Mpc is the sound scale of the fiducial cosmology used in this study. The parameters which are not well constrained by our galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets, e.g., cosmic microwave background (CMB), are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e., CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. The uncertainty on the dark energy equation of state parameter, $w$, from CMB+CMASS is about 8 per cent. The uncertainty on the curvature fraction, $Ω_k$, is 0.3 per cent. We do not find deviation from flat $Λ$CDM.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bayesian redshift-space distortions correction from galaxy redshift surveys

We present a Bayesian reconstruction method which maps a galaxy distribution from redshift-space to real-space inferring the distances of the individual galaxies. The method is based on sampling density fields assuming a lognormal prior with a likelihood given by the negative binomial distribution function modelling stochastic bias. We assume a deterministic bias given by a power law relating the dark matter density field to the expected halo or galaxy field. Coherent redshift-space distortions are corrected in a Gibbs-sampling procedure by moving the galaxies from redshift-space to real-space according to the peculiar motions derived from the recovered density field using linear theory with the option to include tidal field corrections from second order Lagrangian perturbation theory. The virialised distortions are corrected by sampling candidate real-space positions (being in the neighbourhood of the observations along the line of sight), which are compatible with the bulk flow corrected redshift-space position adding a random dispersion term in high density collapsed regions. The latter are defined according to the eigenvalues of the Hessian. This approach presents an alternative method to estimate the distances to galaxies using the three dimensional spatial information, and assuming isotropy. Hence the number of applications is very broad. In this work we show the potential of this method to constrain the growth rate up to $k$ ~ 0.3 $h$ Mpc$^{-1}$. Furthermore it could be useful to correct for photo-metric redshift errors, and to obtain improved BAO reconstructions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraining the baryon fraction in the Warm Hot Intergalactic Medium at low redshifts with PLANCK data

We cross-correlate foreground cleaned Planck Nominal Cosmic Microwave Background (CMB) maps with two templates constructed from the Two-Micron All-Sky Redshift Survey of galaxies. The first template traces the large-scale filamentary distribution characteristic of the Warm-Hot Intergalactic Medium (WHIM) out to ~90 Mpc/h. The second traces preferentially the virialized gas in unresolved halos around galaxies. We find a marginal signal from the correlation of Planck data and the WHIM template with a signal-to-noise from 0.84 to 1.39 at the different Planck frequencies, and with a frequency dependence compatible with the thermal Sunyaev-Zel'dovich (tSZ) effect. When we restrict our analysis to the 60% of the sky outside the plane of the Galaxy and known point sources and galaxy clusters, the cross-correlation at zero lag is 0.064 +/- 0.051 muK. The correlation extends out to ~6 deg, which at the median depth of our template corresponds to a physical length of ~ 6-8 Mpc/h. On the same fraction of the sky, the cross-correlation of the CMB data with the second template is <0.17 muK (95% C.L.), providing no statistically significant evidence of a contribution from bound gas to the previous result. This limit translates into a physical constraint on the properties of the shock-heated WHIM of a log-normal model describing the weakly nonlinear density field. We find that our upper limit is compatible with a fraction of 45% of all baryons residing in filaments at overdensities ~1-100 and with temperatures in the range 10^4.5-10^7.5 K, in agreement with the detection at redshift z~0.5 of van Waerbeke et al. (2014).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraints on the missing baryons from the kinetic Sunyaev-Zeldovich effect in Planck data

We estimate the amount of the {\it missing baryons} detected by the \Planck\ measurements of the cosmic microwave background in the direction of Central Galaxies (CGs) identified in the Sloan galaxy survey. The peculiar motion of the gas inside and around the CGs unveils values of the Thomson optical depth $τ_{\rm T}$ in the range $0.2$--$2\times 10^{-4}$, indicating that the regions probed around CGs contain roughly half of the total amount of baryons in the Universe at the epoch where the CGs are found. If baryons follow dark matter, the measured $τ_{\rm T}$s are compatible with the detection all the baryons existing inside and around the CGs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmic flows and the expansion of the Local Universe from nonlinear phase-space reconstructions

This work investigates the impact of cosmic flows and density perturbations on Hubble constant $H_0$ measurements using nonlinear phase-space reconstructions of the Local Universe (LU). In particular, we rely on 25 constrained N-body simulations of the LU using the 2MRS galaxy sample within distances of 90 Mpc/h. These have been randomly extended up to volumes enclosing 360 Mpc/h with augmented Lagrangian perturbation theory (750 simulations in total), accounting in this way for gravitational mode coupling from larger scales, correcting for periodic boundary effects, and estimating systematics of missing attractors ($σ_{\rm large}=134$ km/s). We report on Local Group speed reconstructions, which are compatible with CMB-dipole measurements: $|v_{\rm LG}|=685\pm137$ km/s. The direction $(l,b)=(260.5\pm 13.3,39.1\pm 10.4)^\circ$ is compatible with observations after considering the variance of large scales. Accounting for large scales, our local bulk flow estimations assuming a $Λ$CDM model are compatible with estimates based on velocity data derived from the Tully-Fisher relation. We focus on low redshift supernova measurements ($0.01<z<0.025$), which have been found to disagree with probes at larger distances. Our analysis indicates that there are two effects contributing to this tension. First, the anisotropic distribution of supernovae aligns with the velocity dipole and induces a systematic boost in $H_0$. Second, a divergent region surrounding the Virgo Supercluster is responsible for an additional positive bias in $H_0$. Taking these effects into account yields a correction of $ΔH_0=-1.76 \pm 0.21 $ km/s/Mpc, thereby reducing the tension between local probes and more distant probes. Effectively $H_0$ is lower by about $2\%$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

DIVE in the cosmic web: voids with Delaunay Triangulation from discrete matter tracer distributions

We present a novel parameter-free cosmological void finder (\textsc{dive}, Delaunay TrIangulation Void findEr) based on Delaunay Triangulation (DT), which efficiently computes the empty spheres constrained by a discrete set of tracers. We define the spheres as DT voids, and describe their properties, including an universal density profile together with an intrinsic scatter. We apply this technique on 100 halo catalogues with volumes of 2.5\,$h^{-1}$Gpc side each, with a bias and number density similar to the BOSS CMASS Luminous Red Galaxies, performed with the \textsc{patchy} code. Our results show that there are two main species of DT voids, which can be characterised by the radius: they have different responses to halo redshift space distortions, to number density of tracers, and reside in different dark matter environments. Based on dynamical arguments using the tidal field tensor, we demonstrate that large DT voids are hosted in expanding regions, whereas the haloes used to construct them reside in collapsing ones. Our approach is therefore able to efficiently determine the troughs of the density field from galaxy surveys, and can be used to study their clustering. We further study the power spectra of DT voids, and find that the bias of the two populations are different, demonstrating that the small DT voids are essentially tracers of groups of haloes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Halo mass distribution reconstruction across the cosmic web

We study the relation between halo mass and its environment from a probabilistic perspective. We find that halo mass depends not only on local dark matter density, but also on non-local quantities such as the cosmic web environment and the halo-exclusion effect. Given these accurate relations, we have developed the HADRON-code (Halo mAss Distribution ReconstructiON), a technique which permits us to assign halo masses to a distribution of haloes in three-dimensional space. This can be applied to the fast production of mock galaxy catalogues, by assigning halo masses, and reproducing accurately the bias for different mass cuts. The resulting clustering of the halo populations agree well with that drawn from the BigMultiDark $N$-body simulation: the power spectra are within 1-$σ$ up to scales of $k=0.2\,h\,{\rm Mpc}^{-1}$, when using augmented Lagrangian perturbation theory based mock catalogues. Only the most massive haloes show a larger deviation. For these, we find evidence of the halo-exclusion effect. A clear improvement is achieved when assigning the highest masses to haloes with a minimum distance separation. We also compute the 2- and 3-point correlation functions, and find an excellent agreement with $N$-body results. Our work represents a quantitative application of the cosmic web classification. It can have further interesting applications in the multi-tracer analysis of the large-scale structure for future galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measuring Baryon Acoustic Oscillations from the clustering of voids

We investigate the necessary methodology to optimally measure the baryon acoustic oscillation (BAO) signal, from voids based on galaxy redshift catalogues. To this end, we study the dependency of the BAO signal on the population of voids classified by their sizes. We find for the first time the characteristic features of the correlation function of voids including the first robust detection of BAOs in mock galaxy catalogues. These show an anti-correlation around the scale corresponding to the smallest size of voids in the sample (the void exclusion effect), and dips at both sides of the BAO peak, which can be used to determine the significance of the BAO signal without any priori model. Furthermore, our analysis demonstrates that there is a scale dependent bias for different populations of voids depending on the radius, with the peculiar property that the void population with the largest BAO significance corresponds to tracers with approximately zero bias on the largest scales. We further investigate the methodology on an additional set of 1,000 realistic mock galaxy catalogues reproducing the SDSS-III/BOSS CMASS DR11 data, to control the impact of sky mask and radial selection function. Our solution is based on generating voids from randoms including the same survey geometry and completeness, and a post-processing cleaning procedure in the holes and at the boundaries of the survey. The methodology and optimal selection of void populations validated in this work have been used to perform the first BAO detection from voids in observations, presented in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

SDSS-III Baryon Oscillation Spectroscopic Survey Data Release 12: galaxy target selection and large scale structure catalogues

The Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey (SDSS) III project, has provided the largest survey of galaxy redshifts available to date, in terms of both the number of galaxy redshifts measured by a single survey, and the effective cosmological volume covered. Key to analysing the clustering of these data to provide cosmological measurements is understanding the detailed properties of this sample. Potential issues include variations in the target catalogue caused by changes either in the targeting algorithm or properties of the data used, the pattern of spectroscopic observations, the spatial distribution of targets for which redshifts were not obtained, and variations in the target sky density due to observational systematics. We document here the target selection algorithms used to create the galaxy samples that comprise BOSS. We also present the algorithms used to create large scale structure catalogues for the final Data Release (DR12) samples and the associated random catalogues that quantify the survey mask. The algorithms are an evolution of those used by the BOSS team to construct catalogues from earlier data, and have been designed to accurately quantify the galaxy sample. The code used, designated MKSAMPLE, is released with this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bayesian inference of cosmic density fields from non-linear, scale-dependent, and stochastic biased tracers

We present a Bayesian reconstruction algorithm to generate unbiased samples of the underlying dark matter field from halo catalogues. Our new contribution consists of implementing a non-Poisson likelihood including a deterministic non-linear and scale-dependent bias. In particular we present the Hamiltonian equations of motions for the negative binomial (NB) probability distribution function. This permits us to efficiently sample the posterior distribution function of density fields given a sample of galaxies using the Hamiltonian Monte Carlo technique implemented in the Argo code. We have tested our algorithm with the Bolshoi $N$-body simulation at redshift $z = 0$, inferring the underlying dark matter density field from sub-samples of the halo catalogue with biases smaller and larger than one. Our method shows that we can draw closely unbiased samples (compatible within 1-$σ$) from the posterior distribution up to scales of about $k$~1 h/Mpc in terms of power-spectra and cell-to-cell correlations. We find that a Poisson likelihood yields reconstructions with power spectra deviating more than 10% at $k$=0.2 h/Mpc. Our reconstruction algorithm is especially suited for emission line galaxy data for which a complex non-linear stochastic biasing treatment beyond Poissonity becomes indispensable.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining the halo bispectrum in real and redshift space from perturbation theory and non-linear stochastic bias

We present a method to produce mock galaxy catalogues with efficient perturbation theory schemes, which match the number density, power spectra and bispectra in real and in redshift space from N-body simulations. The essential contribution of this work is the way in which we constrain the bias parameters in the PATCHY-code. In addition of aiming at reproducing the two-point statistics, we seek the set of bias parameters, which constrain the univariate halo probability distribution function (PDF) encoding higher-order correlation functions. We demonstrate that halo catalogues based on the same underlying dark matter field with a fix halo number density, and accurately matching the power spectrum (within 2%), can lead to very different bispectra depending on the adopted halo bias model. A model ignoring the shape of the halo PDF can lead to deviations up to factors of 2. The catalogues obtained additionally constraining the shape of the halo PDF can significantly lower the discrepancy in the three-point statistics, yielding closely unbiased bispectra both in real and in redshift space; which are in general compatible with those corresponding to an N-body simulation within 10% (deviating at most up to 20%). Our calculations show that the constant linear bias of ~2 for Luminous Red Galaxy (LRG) like galaxies seen in the power spectrum, mainly comes from sampling halos in high density peaks, choosing a high density threshold rather than from a factor multiplying the dark matter density field. Our method contributes towards an efficient modelling of the halo/galaxy distribution required to estimate uncertainties in the clustering measurements from galaxy redshift surveys. We have also demonstrated that it represents a powerful tool to test various bias models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hunting down systematics in baryon acoustic oscillations after cosmic high noon

Future dark energy experiments will require better and more accurate theoretical predictions for the baryonic acoustic oscillations (BAO) signature in the spectrum of cosmological perturbations. Here, we use large N-body simulations of the \LambdaCDM Planck cosmology to study any possible systematic shifts and damping in BAO due to the impact of nonlinear gravitational growth of structure, scale dependent and non-local bias, and redshift-space distortions. The effect of cosmic variance is largely reduced by dividing the tracer power spectrum by that from a BAO-free simulation starting with the same phases. This permits us to study with unprecedented accuracy (better than 0.02% for dark matter and 0.07% for low-bias halos) small shifts of the pristine BAO wavenumbers towards larger k, and non-linear damping of BAO wiggles in the power spectrum of dark matter and halo populations in the redshift range z=0-1. For dark matter, we provide an accurate parametrization of the evolution of αas a function of the linear growth factor D(z). For halo samples, with bias ranging from 1.2 to 2.8, we measure a typical BAO shift of ~0.25%, observed in real-space, which does not show an appreciable evolution with redshift within the uncertainties. Moreover, we report a constant shift as a function of halo bias. We find a different evolution of the damping of the acoustic feature in all halo samples as compared to dark matter with haloes suffering less damping, and also find some weak dependence on bias. A larger BAO shift and damping is measured in redshift-space which can be well explained by linear theory due to redshift-space distortions. A clear modulation in phase with the acoustic scale is observed in the scale-dependent halo bias due to the presence of the baryonic acoustic oscillations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cosmological Structure Formation with Augmented Lagrangian Perturbation Theory

We present a new fast and efficient approach to model structure formation with Augmented Lagrangian Perturbation Theory (ALPT). Our method is based on splitting the displacement field into a long and a short-range component. The long-range component is computed by second order LPT (2LPT). This approximation contains a tidal nonlocal and nonlinear term. Unfortunately, 2LPT fails on small scales due to severe shell crossing and a crude quadratic behaviour in the low density regime. The spherical collapse (SC) approximation has been recently reported to correct for both effects by adding an ideal collapse truncation. However, this approach fails to reproduce the structures on large scales where it is significantly less correlated with the N-body result than 2LPT or linear LPT (the Zeldovich approximation). We propose to combine both approximations using for the short-range displacement field the SC solution. A Gaussian filter with a smoothing radius r_S is used to separate between both regimes. We use the result of 25 dark matter only N-body simulations to benchmark at z=0 the different approximations: 1st, 2nd, 3rd order LPT, SC and our novel combined ALPT model. This comparison demonstrates that our method improves previous approximations at all scales showing ~25% and ~75% higher correlation than 2LPT with the N-body solution at k = 1 and 2 h Mpc^-1, respectively. We conduct a parameter study to determine the optimal range of smoothing radii and find that the maximum correlation is achieved with r_S = 4 - 5 h^-1 Mpc. This structure formation approach could be used for various purposes, such as setting-up initial conditions for N-body simulations, generating mock galaxy catalogues, cosmic web analysis or for reconstructions of the primordial density fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modelling Baryon Acoustic Oscillations with Perturbation Theory and Stochastic Halo Biasing

In this work we investigate the generation of mock halo catalogues based on perturbation theory and nonlinear stochastic biasing with the novel PATCHY-code. In particular, we use Augmented Lagrangian Perturbation Theory (ALPT) to generate a dark matter density field on a mesh starting from Gaussian fluctuations and to compute the peculiar velocity field. ALPT is based on a combination of second order LPT (2LPT) on large scales and the spherical collapse model on smaller scales. We account for the systematic deviation of perturbative approaches from N-body simulations together with halo biasing adopting an exponential bias model. We then account for stochastic biasing by defining three regimes: a low, an intermediate and a high density regime, using a Poisson distribution in the intermediate regime and the negative binomial distribution to model over-dispersion in the high density regime. Since we focus in this study on massive halos, we suppress the generation of halos in the low density regime. The various nonlinear and stochastic biasing parameters, and density thresholds (five) are calibrated with the large BigMultiDark N-body simulation to match the power spectrum of the corresponding halo population. Our mock catalogues show power spectra, both in real- and redshift-space, which are compatible with N-body simulations within about 2% up to k ~ 1 h Mpc^-1 at z = 0.577 for a sample of halos with the typical BOSS CMASS galaxy number density. The corresponding correlation functions are compatible down to a few Mpc. We also find that neglecting over-dispersion in high density regions produces power spectra with deviations of 10% at k ~ 0.4 h Mpc^-1. These results indicate the need to account for an accurate statistical description of the galaxy clustering for precise studies of large-scale surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmic Structure and Dynamics of the Local Universe

We present a cosmography analysis of the Local Universe based on the recently released Two-Micron All-Sky Redshift Survey (2MRS). Our method is based on a Bayesian Networks Machine Learning algorithm (the Kigen-code) which self-consistently samples the initial density fluctuations compatible with the observed galaxy distribution and a structure formation model given by second order Lagrangian perturbation theory (2LPT). From the initial conditions we obtain an ensemble of reconstructed density and peculiar velocity fields which characterize the local cosmic structure with high accuracy unveiling nonlinear structures like filaments and voids in detail. Coherent redshift space distortions are consistently corrected within 2LPT. From the ensemble of cross-correlations between the reconstructions and the galaxy field and the variance of the recovered density fields we find that our method is extremely accurate up to k ~ 1 h Mpc^-1 and still yields reliable results down to scales of about 3-4 h^-1 Mpc. The motion of the local group we obtain within ~ 80 h^-1 Mpc (v_LG=522+-86 km s^-1, l_LG=291^o +- 16^o, b_LG=34^o+-8^o) is in good agreement with measurements derived from the CMB and from direct observations of peculiar motions and is consistent with the predictions of LambdaCDM.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Estimating cosmic velocity fields from density fields and tidal tensors

In this work we investigate the nonlinear and nonlocal relation between cosmological density and peculiar velocity fields. Our goal is to provide an algorithm for the recon- struction of the nonlinear velocity field from the fully nonlinear density. We find that including the gravitational tidal field tensor using second order Lagrangian perturba- tion theory (2LPT) based upon an estimate of the linear component of the nonlinear density field significantly improves the estimate of the cosmic flow in comparison to linear theory not only in the low density, but also and more dramatically in the high density regions. In particular we test two estimates of the linear component: the log- normal model and the iterative Lagrangian linearisation. The present approach relies on a rigorous higher order Lagrangian perturbation theory analysis which incorpo- rates a nonlocal relation. It does not require additional fitting from simulations being in this sense parameter free, it is independent of statistical-geometrical optimisation and it is straightforward and efficient to compute. The method is demonstrated to yield an unbiased estimator of the velocity field on scales ~> 5 Mpc/h with closely Gaussian distributed errors. Moreover, the statistics of the divergence of the peculiar velocity field is extremely well recovered showing a good agreement with the true one from N-body simulations. The typical errors of about 10 km/s (1 sigma confidence intervals) are reduced by more than 80% with respect to linear theory in the scale range between 5 and 10 Mpc/h in high density regions (δ > 2). We also find that iterative Lagrangian linearisation is significantly superior in the low density regime with respect to the lognormal model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multiscale Inference of Matter Fields and Baryon Acoustic Oscillations from the Ly-alpha Forest

We present a novel Bayesian method for the joint reconstruction of cosmological matter density fields, peculiar velocities and power-spectra in the quasi-nonlinear regime. We study its applicability to the Ly-alpha forest based on multiple quasar absorption spectra. Our approach to this problem includes a multiscale, nonlinear, two-step scheme since the statistics describing the matter distribution depends on scale, being strongly non-Gaussian on small scales (< 0.1 h^{-1} Mpc) and closely lognormal on scales >~10 h^{-1} Mpc. The first step consists on performing 1D highly resolved matter density reconstructions along the line-of-sight towards z~2-3 quasars based on an arbitrary non-Gaussian univariate model for matter statistics. The second step consists on Gibbs-sampling based on conditional PDFs. The matter density field is sampled in real space with Hamiltonian-sampling using the Poisson/Gamma-lognormal model, while redshift distortions are corrected with linear Lagrangian perturbation theory. The power-spectrum of the lognormal transformed variable which is Gaussian distributed (and thus close to the linear regime) can consistently be sampled with the inverse Gamma distribution function. We test our method through numerical N-body simulations with a computational volume large enough (> 1 h^{-3} Gpc^3) to show that the linear power-spectra are nicely recovered over scales larger than >~20 h^{-1} Mpc, i.e. the relevant range where features imprinted by the baryon-acoustics oscillations (BAOs) appear.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Initial Conditions of the Universe from Constrained Simulations

I present a new approach to recover the primordial density fluctuations and the cosmic web structure underlying a galaxy distribution. The method is based on sampling Gaussian fields which are compatible with a galaxy distribution and a structure formation model. This is achieved by splitting the inversion problem into two Gibbs-sampling steps: the first being a Gaussianisation step transforming a distribution of point sources at Lagrangian positions -which are not a priori given- into a linear alias-free Gaussian field. This step is based on Hamiltonian sampling with a Gaussian-Poisson model. The second step consists on a likelihood comparison in which the set of matter tracers at the initial conditions is constrained on the galaxy distribution and the assumed structure formation model. For computational reasons second order Lagrangian Perturbation Theory is used. However, the presented approach is flexible to adopt any structure formation model. A semi-analytic halo-model based galaxy mock catalog is taken to demonstrate that the recovered initial conditions are closely unbiased with respect to the actual ones from the corresponding N-body simulation down to scales of a ~ 5 Mpc/h. The cross-correlation between them shows a substantial gain of information, being at k ~ 0.3 h/Mpc more than doubled. In addition the initial conditions are extremely well Gaussian distributed and the power-spectra follow the shape of the linear power-spectrum being very close to the actual one from the simulation down to scales of k ~ 1 h/Mpc.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bayesian analysis of cosmic structures

We revise the Bayesian inference steps required to analyse the cosmological large-scale structure. Here we make special emphasis in the complications which arise due to the non-Gaussian character of the galaxy and matter distribution. In particular we investigate the advantages and limitations of the Poisson-lognormal model and discuss how to extend this work. With the lognormal prior using the Hamiltonian sampling technique and on scales of about 4 h^{-1} Mpc we find that the over-dense regions are excellent reconstructed, however, under-dense regions (void statistics) are quantitatively poorly recovered. Contrary to the maximum a posteriori (MAP) solution which was shown to over-estimate the density in the under-dense regions we obtain lower densities than in N-body simulations. This is due to the fact that the MAP solution is conservative whereas the full posterior yields samples which are consistent with the prior statistics. The lognormal prior is not able to capture the full non-linear regime at scales below ~ 10 h^{-1} Mpc for which higher order correlations would be required to describe the matter statistics. However, we confirm as it was recently shown in the context of Ly-alpha forest tomography that the Poisson-lognormal model provides the correct two-point statistics (or power-spectrum).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cosmic density field reconstruction from Ly-alpha forest data

We present a novel, fast method to recover the density field through the statistics of the transmitted flux in high redshift quasar absorption spectra. The proposed technique requires the computation of the probability distribution function of the transmitted flux (P_F) in the Ly-alpha forest region and, as a sole assumption, the knowledge of the probability distribution function of the matter density field (P_Delta). We show that the probability density conservation of the flux and matter density unveils a flux-density (F-Delta) relation which can be used to invert the Ly-alpha forest without any assumption on the physical properties of the intergalactic medium. We test our inversion method at z=3 through the following steps: [i] simulation of a sample of synthetic spectra for which P_Delta is known; [ii] computation of P_F; [iii] inversion of the Ly-alpha forest through the F-Delta relation. Our technique, when applied to only 10 observed spectra characterized by a signal-to noise ratio S/N >= 100 provides an exquisite (relative error epsilon_Delta <~ 12 % in >~ 50 % of the pixels) reconstruction of the density field in >~ 90 % of the line of sight. We finally discuss strengths and limitations of the method.

0 citations0 reviewsNo code linkedOpen access