A photometric redshift sample of Luminous Red Galaxies (hereafter LRGs) obtained from The DECam Legacy Survey (DECaLS) is analysed to probe cosmic distances by exploiting the wedge approach of the two-point correlation function. Although the cosmological information is highly contaminated by the uncertainties existing in the photometric redshifts from the galaxy map, an angular diameter distance can be probed at the perpendicular configuration in which the measured correlation function is minimally contaminated. An ensemble of wedged correlation functions selected up to a given threshold based on having the least contamination was studied in the previous work (arXiv:1903.09651v2 [astro-ph.CO]) using simulations, and the extracted cosmological information was unbiased within this threshold. We apply the same methodology for analysing the LRG sample from DECaLS which will provide the optical imaging for targeting two-thirds of the DESI footprint and measure the angular diameter distances at $z=0.69$ and $z=0.87$ to be $D_{A}(0.697)=(1499 \pm 77\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ and $D_{A}(0.874)=(1680 \pm 109\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ with a fractional error of 5.14% and 6.48% respectively. We obtain a value of $H_{0}=67.59\pm5.52$ km/s/Mpc which supports the $H_0$ measured by all other BAO results and is consistent with $\Lambda$CDM model.

We present CosmoHub (https://cosmohub.pic.es), a web application based on Hadoop to perform interactive exploration and distribution of massive cosmological datasets. Recent Cosmology seeks to unveil the nature of both dark matter and dark energy mapping the large-scale structure of the Universe, through the analysis of massive amounts of astronomical data, progressively increasing during the last (and future) decades with the digitization and automation of the experimental techniques. CosmoHub, hosted and developed at the Port d'Informaci\'o Cient\'ifica (PIC), provides support to a worldwide community of scientists, without requiring the end user to know any Structured Query Language (SQL). It is serving data of several large international collaborations such as the Euclid space mission, the Dark Energy Survey (DES), the Physics of the Accelerating Universe Survey (PAUS) and the Marenostrum Institut de Ci\`encies de l'Espai (MICE) numerical simulations. While originally developed as a PostgreSQL relational database web frontend, this work describes the current version of CosmoHub, built on top of Apache Hive, which facilitates scalable reading, writing and managing huge datasets. As CosmoHub's datasets are seldomly modified, Hive it is a better fit. Over 60 TiB of catalogued information and $50 \times 10^9$ astronomical objects can be interactively explored using an integrated visualization tool which includes 1D histogram and 2D heatmap plots. In our current implementation, online exploration of datasets of $10^9$ objects can be done in a timescale of tens of seconds. Users can also download customized subsets of data in standard formats generated in few minutes.

Context. Studies of galaxy pairs can provide valuable information to jointly understand the formation and evolution of galaxies and galaxy groups. Consequently, taking into account the new high precision photo-z surveys, it is important to have reliable and tested methods that allow us to properly identify these systems and estimate their total masses and other properties. Aims. In view of the forthcoming Physics of the Accelerating Universe Survey (PAUS) we propose and evaluate the performance of an identification algorithm of projected close isolated galaxy pairs. We expect that the photometric selected systems can adequately reproduce the observational properties and the inferred lensing mass - luminosity relation of a pair of truly bound galaxies that are hosted by the same dark matter halo. Methods. We develop an identification algorithm that considers the projected distance between the galaxies, the projected velocity difference and an isolation criteria in order to restrict the sample to isolated systems. We apply our identification algorithm using a mock galaxy catalog that mimics the features of PAUS. To evaluate the feasibility of our pair finder, we compare the identified photometric samples with a test sample that considers that both members are included in the same halo. Also, taking advantage of the lensing properties provided by the mock catalog, we apply a weak lensing analysis to determine the mass of the selected systems. Results. Photometric selected samples tend to show high purity values, but tend to misidentify truly bounded pairs as the photometric redshift errors increase. Nevertheless, overall properties such as the luminosity and mass distributions are successfully reproduced. We also accurately reproduce the lensing mass - luminosity relation as expected for galaxy pairs located in the same halo.

What are the mass and galaxy profiles of cosmic voids? In this paper we use two methods to extract voids in the Dark Energy Survey (DES) Year 1 redMaGiC galaxy sample to address this question. We use either 2D slices in projection, or the 3D distribution of galaxies based on photometric redshifts to identify voids. For the mass profile, we measure the tangential shear profiles of background galaxies to infer the excess surface mass density. The signal-to-noise ratio for our lensing measurement ranges between 10.7 and 14.0 for the two void samples. We infer their 3D density profiles by fitting models based on N-body simulations and find good agreement for void radii in the range 15-85 Mpc. Comparison with their galaxy profiles then allows us to test the relation between mass and light at the 10%-level, the most stringent test to date. We find very similar shapes for the two profiles, consistent with a linear relationship between mass and light both within and outside the void radius. We validate our analysis with the help of simulated mock catalogues and estimate the impact of photometric redshift uncertainties on the measurement. Our methodology can be used for cosmological applications, including tests of gravity with voids. This is especially promising when the lensing profiles are combined with spectroscopic measurements of void dynamics via redshift-space distortions.

Photometric galaxy surveys constitute a powerful cosmological probe but rely on the accurate characterization of their redshift distributions using only broadband imaging, and can be very sensitive to incomplete or biased priors used for redshift calibration. S\'anchez & Bernstein (2019) presented a hierarchical Bayesian model which estimates those from the robust combination of prior information, photometry of single galaxies and the information contained in the galaxy clustering against a well-characterized tracer population. In this work, we extend the method so that it can be applied to real data, developing some necessary new extensions to it, especially in the treatment of galaxy clustering information, and we test it on realistic simulations. After marginalizing over the mapping between the clustering estimator and the actual density distribution of the sample galaxies, and using prior information from a small patch of the survey, we find the incorporation of clustering information with photo-$z$'s to tighten the redshift posteriors, and to overcome biases in the prior that mimic those happening in spectroscopic samples. The method presented here uses all the information at hand to reduce prior biases and incompleteness. Even in cases where we artificially bias the spectroscopic sample to induce a shift in mean redshift of $\Delta \bar z \approx 0.05,$ the final biases in the posterior are $\Delta \bar z \lesssim0.003.$ This robustness to flaws in the redshift prior or training samples would constitute a milestone for the control of redshift systematic uncertainties in future weak lensing analyses.

The largest structures in the cosmic web probe the dynamical nature of dark energy through their integrated Sachs-Wolfe imprints. In the strength of the signal, typical cosmic voids have shown good consistency with expectation $A_{\rm ISW}=\Delta T^{\rm data} / \Delta T^{\rm theory}=1$, given the substantial cosmic variance. Discordantly, large-scale hills in the gravitational potential, or supervoids, have shown excess signals. In this study, we mapped out 87 new supervoids in the total 5000 deg$^2$ footprint of the Dark Energy Survey at $0.2<z<0.9$ to probe these anomalous claims. We found an excess imprinted profile with $ A_{\rm ISW}\approx4.1\pm2.0$ amplitude. The combination with independent BOSS data reveals an ISW imprint of supervoids at the $3.3\sigma$ significance level with an enhanced $A_{\rm ISW}\approx5.2\pm1.6$ amplitude. The tension with $\Lambda$CDM predictions is equivalent to $2.6\sigma$ and remains unexplained.

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.

Small temperature anisotropies in the Cosmic Microwave Background can be sourced by density perturbations via the late-time integrated Sachs-Wolfe effect. Large voids and superclusters are excellent environments to make a localized measurement of this tiny imprint. In some cases excess signals have been reported. We probed these claims with an independent data set, using the first year data of the Dark Energy Survey in a different footprint, and using a different super-structure finding strategy. We identified 52 large voids and 102 superclusters at redshifts $0.2 < z < 0.65$. We used the Jubilee simulation to a priori evaluate the optimal ISW measurement configuration for our compensated top-hat filtering technique, and then performed a stacking measurement of the CMB temperature field based on the DES data. For optimal configurations, we detected a cumulative cold imprint of voids with $\Delta T_{f} \approx -5.0\pm3.7~\mu K$ and a hot imprint of superclusters $\Delta T_{f} \approx 5.1\pm3.2~\mu K$ ; this is $\sim1.2\sigma$ higher than the expected $|\Delta T_{f}| \approx 0.6~\mu K$ imprint of such super-structures in $\Lambda$CDM. If we instead use an a posteriori selected filter size ($R/R_{v}=0.6$), we can find a temperature decrement as large as $\Delta T_{f} \approx -9.8\pm4.7~\mu K$ for voids, which is $\sim2\sigma$ above $\Lambda$CDM expectations and is comparable to previous measurements made using SDSS super-structure data.

Using dark matter simulations we show how halo bias is determined by local density and not by halo mass. This is not totally surprising, as according to the peak-background split model, local density is the property that constraints bias at large scales. Massive haloes have a high clustering because they reside in high density regions. Small haloes can be found in a wide range of environments which determine their clustering amplitudes differently. This contradicts the assumption of standard Halo Occupation Distribution (HOD) models that the bias and occupation of haloes is determined solely by their mass. We show that the bias of central galaxies from semi-analytic models of galaxy formation as a function of luminosity and colour is not correctly predicted by the standard HOD model. Using local density instead of halo mass the HOD model correctly predicts galaxy bias. These results indicate the need to include information about local density and not only mass in order to correctly apply HOD analysis in these galaxy samples. This new model can be readily applied to observations and has the advantage that the galaxy density can be directly observed, in contrast with the dark matter halo mass.

We present a new method to measure the redshift-dependent galaxy bias by combining information from the galaxy density field and the weak lensing field. This method is based on Amara et al. (2012), where they use the galaxy density field to construct a bias-weighted convergence field kg. The main difference between Amara et al. (2012) and our new implementation is that here we present another way to measure galaxy bias using tomography instead of bias parameterizations. The correlation between kg and the true lensing field k allows us to measure galaxy bias using different zero-lag correlations, such as <kgk>/<kk> or <kgkg>/<kgk>. Our method measures the linear bias factor on linear scales under the assumption of no stochasticity between galaxies and matter. We use the MICE simulation to measure the linear galaxy bias for a flux-limited sample (i < 22.5) in tomographic redshift bins using this method. This paper is the first that studies the accuracy and systematic uncertainties associated with the implementation of the method, and the regime where it is consistent with the linear galaxy bias defined by projected 2-point correlation functions (2PCF). We find that our method is consistent with linear bias at the percent level for scales larger than 30 arcmin, while nonlinearities appear at smaller scales. This measurement is a good complement to other measurements of bias, since it does not depend strongly on sigma8 as the 2PCF measurements. We apply this method to the Dark Energy Survey Science Verification data in a follow-up paper.

It is well known that the probability distribution function (PDF) of galaxy density contrast is approximately lognormal; whether the PDF of mass fluctuations derived from weak lensing convergence (kappa_WL) is lognormal is less well established. We derive PDFs of the galaxy and projected matter density distributions via the Counts in Cells (CiC) method. We use maps of galaxies and weak lensing convergence produced from the Dark Energy Survey (DES) Science Verification data over 139 deg^2. We test whether the underlying density contrast is well described by a lognormal distribution for the galaxies, the convergence and their joint PDF. We confirm that the galaxy density contrast distribution is well modeled by a lognormal PDF convolved with Poisson noise at angular scales from 10-40 arcmin (corresponding to physical scales of 3-10 Mpc). We note that as kappa_WL is a weighted sum of the mass fluctuations along the line of sight, its PDF is expected to be only approximately lognormal. We find that the kappa_WL distribution is well modeled by a lognormal PDF convolved with Gaussian shape noise at scales between 10 and 20 arcmin, with a best-fit chi^2/DOF of 1.11 compared to 1.84 for a Gaussian model, corresponding to p-values 0.35 and 0.07 respectively, at a scale of 10 arcmin. Above 20 arcmin a simple Gaussian model is sufficient. The joint PDF is also reasonably fitted by a bivariate lognormal. As a consistency check we compare the variances derived from the lognormal modelling with those directly measured via CiC. Our methods are validated against maps from the MICE Grand Challenge N-body simulation.

We study halo clustering bias with second- and third-order statistics of halo and matter density fields in the MICE Grand Challenge simulation. We verify that two-point correlations deliver reliable estimates of the linear bias parameters at large scales, while estimations from the variance can be significantly affected by non-linear and possibly non-local contributions to the bias function. Combining three-point auto- and cross-correlations we find, for the first time in configuration space, evidence for the presence of such non-local contributions. These contributions are consistent with predicted second-order non-local effects on the bias functions originating from the dark matter tidal field. Samples of massive haloes show indications of bias (local or non-local) beyond second order. Ignoring non-local bias causes $20-30$\% and $5-10$\% overestimation of the linear bias from three-point auto- and cross-correlations respectively. We study two third-order bias estimators which are not affected by second-order non-local contributions. One is a combination of three-point auto- and cross- correlation. The other is a combination of third-order one- and two-point cumulants. Both methods deliver accurate bias estimations of the linear bias. Furthermore their estimations of second-order bias agree mutually. Ignoring non-local bias causes higher values of the second-order bias from three-point correlations. Our results demonstrate that third-order statistics can be employed for breaking the growth-bias degeneracy.

The Physics of the Accelerating Universe (PAU) survey at the William Herschel Telescope (WHT) will use a new optical camera (PAUCam) with a large set of narrow-band filters to perform a photometric galaxy survey with a quasi-spectroscopic redshift precision of \sigma(z)/(1 + z) ~ 0.0035 and map the large-scale structure of the universe in three dimensions up to i_AB < 22.5-23.0. In this paper we present a detailed photo-z performance study using photometric simulations for 40 equally-spaced 12.5-nm-wide (FWHM) filters with a ~25% overlap and spanning the wavelength range from 450 nm to 850 nm, together with a ugrizY broad-band filter system. We then present the migration matrix r_ij, containing the probability that a galaxy in a true redshift bin j is measured in a photo-z bin i, and study its effect on the determination of galaxy auto- and cross-correlations. Finally, we also study the impact on the photo-z performance of small variations of the filter set in terms of width, wavelength coverage, etc., and find a broad region where slightly modified filter sets provide similar results, with the original set being close to optimal.

The magnification effect of gravitational lensing is a powerful probe of the distribution of matter in the universe, yet it is frequently overlooked due to the fact that its signal to noise is smaller than that of lensing shear. Because its systematic errors are quite different from those of shear, magnification is nevertheless an important approach with which to study the distribution of large scale structure. We present lensing mass profiles of spectroscopic luminous red galaxies (LRGs) and galaxy clusters determined through measurements of the weak lensing magnification of photometric LRGs in their background. We measure the change in detected galaxy counts as well as the increased average galaxy flux behind the lenses. In addition, we examine the average change in source color due to extinction by dust in the lenses. By simultaneously fitting these three probes we constrain the mass profiles and dust-to-mass ratios of the lenses in six bins of lens richness. For each richness bin we fit an NFW halo mass, brightest cluster galaxy (BCG) mass, second halo term, and dust-to-mass ratio. The resulting mass-richness relation is consistent with previous analyses of the catalogs, and limits on the dust-to-mass ratio in the lenses are in agreement with expectations. We explore the effects of including the (low signal-to-noise) flux magnification and reddening measurements in the analysis compared to using only the counts magnification data; the additional probes significantly improve the agreement between our measured mass-richness relation and previous results.

Photometric galaxy surveys are an essential tool to further our understanding of the large-scale structure of the universe, its matter and energy content and its evolution. These surveys necessitate the determination of the galaxy redshifts using photometric techniques (photo-z). Oftentimes, it is advantageous to remove from the galaxy sample those for which one suspects that the photo-z estimation might be unreliable. In this paper, we show that applying these photo-z quality cuts blindly can grossly bias the measured galaxy correlations within and across photometric redshift bins. We then extend the work of Ho et al. (2012) and Ross et al. (2011) to develop a simple and effective method to correct for this using the data themselves. Finally, we apply the method to the Mega-Z catalog, containing about a million luminous red galaxies in the redshift range 0.45 < z < 0.65. After splitting the sample into four \Delta z = 0.05 photo-z bins using the BPZ algorithm, we see how our corrections bring the measured galaxy auto- and cross-correlations into agreement with expectations. We then look for the BAO feature in the four bins, with and without applying the photo-z quality cuts, and find a broad agreement between the BAO scales extracted in both cases. Intriguingly, we observe a correlation between galaxy density and photo-z quality even before any photo-z quality cuts are applied. This may be due to uncorrected observational effects that result in correlated gradients across the sky of the galaxy density and the galaxy photo-z precision. Our correction procedure could also help to mitigate some of these systematic effects.

Near future cosmology will see the advent of wide area photometric galaxy surveys, like the Dark Energy Survey (DES), that extent to high redshifts (z ~ 1 - 2) but with poor radial distance resolution. In such cases splitting the data into redshift bins and using the angular correlation function $w(\theta)$, or the $C_{\ell}$ power spectrum, will become the standard approach to extract cosmological information or to study the nature of dark energy through the Baryon Acoustic Oscillations (BAO) probe. In this work we present a detailed model for $w(\theta)$ at large scales as a function of redshift and bin width, including all relevant effects, namely nonlinear gravitational clustering, bias, redshift space distortions and photo-z uncertainties. We also present a model for the full covariance matrix characterizing the angular correlation measurements, that takes into account the same effects as for $w(\theta)$ and also the possibility of a shot-noise component and partial sky coverage. Provided with a large volume N-body simulation from the MICE collaboration we built several ensembles of mock redshift bins with a sky coverage and depth typical of forthcoming photometric surveys. The model for the angular correlation and the one for the covariance matrix agree remarkably well with the mock measurements in all configurations. The prospects for a full shape analysis of $w(\theta)$ at BAO scales in forthcoming photometric surveys such as DES are thus very encouraging.

Clustering statistics are compared in the Automatic Plate Machine (APM) and the Edinburgh/Durham Southern Galaxy Catalogue (EDSGC) angular galaxy surveys. Both surveys were independently constructed from scans of the same adjacent UK IIIa--J Schmidt photographic plates with the APM and COSMOS microdensitometers, respectively. The comparison of these catalogs is a rare practical opportunity to study systematic errors, which cannot be achieved via simulations or theoretical methods. On intermediate scales, $0.1^\circ < \theta < 0.5^\circ$, we find good agreement for the cumulants or reduced moments of counts in cells up to sixth order. On larger scales there is a small disagreement due to edge effects in the EDSGC, which covers a smaller area. On smaller scales, we find a significant disagreement that can only be attributed to differences in the construction of the surveys, most likely the dissimilar deblending of crowded fields. The overall agreement of the APM and EDSGC is encouraging, and shows that the results for intermediate scales should be fairly robust. On the other hand, the systematic deviations found at small scales are significant in a regime, where comparison with theory and simulations is possible. This is an important fact to bear in mind when planning the construction of future digitized galaxy catalogs.

We analyze the uncertainties in the amplitudes of the spatial correlation functions estimated from angular correlations in a sample from the APM Galaxy Survey, with $b_J=17-20$. We model the uncertainties in the selection function and in the evolution of clustering. In particular we estimate $\sigma_8^{APM}$, the rms galaxy number fluctuations in spheres of radius at $8 \Mpc$, from the measured angular variance in the APM. The uncertainty in $\sigma_8^{APM}$ has three main contributions: 8\% from sampling and selection function uncertainties, 7\% from the uncertainty in the evolution of clustering and 3\% from the uncertainty in the value of $\Omega_0$. Including all these contributions, we find $\sigma_8^{APM}$ is in the range $0.78-1.08$. If the galaxy clustering in the APM evolves as expected from gravitational clustering of matter fluctuations, then $\sigma_8^{APM}=0.95 \pm 0.07$ ($1.00 \pm 0.08$) for $\Omega_0 \simeq 1$ ($\Omega_0 \simeq 0$), close to the values for nearby optical samples. On the other hand, if we assume that clustering evolution is fixed in comoving coordinates $\sigma_8^{APM}=0.83 \pm 0.05$ ($0.87 \pm 0.06$), closer to the results for nearby IRAS samples. The final uncertainty in the range of values for the hierarchical amplitudes $S_J\equiv \xibar_J/\xibar_2^{J-1}$ is typically twice the estimated sampling errors, with the highest values for the case of less clustering evolution. We compare our estimates with other results and discuss the implications for models of structure formation.

We estimate the variance $\xibar_2$, the skewness $\xibar_3$ and the kurtosis $\xibar_4$ in the distribution of density fluctuations in a complete sample from the APM Cluster Redshift Survey with 339 clusters and a mean depth $ \sim 250\Mpc$. We are able to measure the statistics of fluctuations in spheres of radius $R \simeq 5-80 \Mpc$, with reasonable errorbars. The statistics in the cluster distribution follow the hierarchical pattern $\xibar_J=S_J~\xibar_2^{J-1}$ with $S_J$ roughly constant, $S_3 \simeq 2$ and $S_4 \sim 8$. We analyse the distribution of clusters taken from N-body simulations of different dark matter models. The results are compared with an alternative method of simulating clusters which uses the truncated Zel'dovich approximation. We argue that this alternative method is not reliable enough for making quantitative predictions of $\xibar$. The N-body simulation results follow similar hierarchical relations to the observations, with $S_J$ almost unaffected by redshift distortions from peculiar motions. The standard $\Omega=1$ Cold Dark Matter (CDM) model is inconsistent with either the second, third or fourth order statistics at all scales. However both a hybrid Mixed Dark Matter model and a low density CDM variant agree with the $\xibar_J$ observations.