The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

In this letter we investigate the connection between the cosmic web and the halo distribution through the gravitational potential. We combine three fields of research, cosmic web classification, perturbation theory expansions of the halo bias, and halo (galaxy) mock catalog making methods. In particular, we use the invariants of the tidal field tensor as generating functions (dubbed I-web), to reproduce the halo number counts of a reference catalog from full gravity calculations, populating the dark matter field on a mesh well into the non-linear regime ($\sim5$ Mpc scales). Our results show an unprecedented agreement with the reference power spectrum within 0.5% up to $k=0.72\,h$ Mpc$^{-1}$. By analysing the three point statistics on large scales (configurations of up to $k=0.2\,h$ Mpc$^{-1}$), we find evidence for non-local bias at the 4.8 $\sigma$ confidence level (an information gain of $\sim$ 3.4 $\sigma$ over the commonly used T-web), being fully compatible with the reference catalog. In particular, we find that a detailed description of tidal anisotropic clustering on large scales is crucial to achieve this accuracy. We conclude that the I-web can potentially be useful to study the cosmic web, to improve the generation of mock galaxy catalogs, to improve on halo mass reconstructions, to study primordial non-Gaussianities, to develop new effective Eulerian galaxy bias models at the field level, and to investigate galaxy evolution improving on environmental studies.

In this paper we demonstrate that the information encoded in \emph{one} single (sufficiently large) $N$-body simulation can be used to reproduce arbitrary numbers of halo catalogues, using approximated realisations of dark matter density fields with different initial conditions. To this end we use as a reference one realisation (from an ensemble of $300$) of the Minerva $N$-body simulations and the recently published Bias Assignment Method to extract the local and non-local bias linking the halo to the dark matter distribution. We use an approximate (and fast) gravity solver to generate $300$ dark matter density fields from the down-sampled initial conditions of the reference simulation and sample each of these fields using the halo-bias and a kernel, both calibrated from the arbitrarily chosen realisation of the reference simulation. We show that the power spectrum, its variance and the three-point statistics are reproduced within $\sim 2\%$ (up to $k\sim1.0\,h\,{\rm Mpc}^{-1}$), $\sim 5-10\%$ and $\sim 10\%$, respectively. Using a model for the real space power spectrum (with three free bias parameters), we show that the covariance matrices obtained from our procedure lead to parameter uncertainties that are compatible within $\sim 10\%$ with respect to those derived from the reference covariance matrix, and motivate approaches that can help to reduce these differences to $\sim 1\%$. Our method has the potential to learn from one simulation with moderate volumes and high-mass resolution and extrapolate the information of the bias and the kernel to larger volumes, making it ideal for the construction of mock catalogues for present and forthcoming observational campaigns such as Euclid or DESI.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

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.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Aiming to extract cosmological information from linear scales of the WISE$\times$SuperCOSMOS photometric redshift catalog, we perform a characterization of the systematic effects associated with stellar content, evidencing the presence of contamination and obscuration. We create an integrated model for these effects (which together we call `usurper contamination'), devise a method to remove both of them simultaneously and show its functionality by applying it to a set of mock catalogs. When administered to WISE$\times$SuperCOSMOS data, our method shows to improve the measurements of angular power spectra on scales $\ell\lesssim15$ and the extraction of cosmological parameters therefrom, even though a significant excess of power remains at these scales. When ignoring scales $\ell<15$, we still find strong indications of systematics, albeit these can be localized in the southern equatorial hemisphere. An independent analysis of the northern hemisphere at $\ell\geq 15$ agrees with a $\Lambda$CDM model with parameters from the Planck satellite and gives $\Omega_{\mathrm{c}}=0.254\pm0.020$ and $\Omega_{\mathrm{b}}<0.065$ at 95% confidence limit when combined with priors on $H_0$, $A_s$ and $n_s$.

This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.

We investigate the clustering properties of radio sources in the Alternative Data Release 1 of the TIFR GMRT Sky Survey (TGSS), focusing on large angular scales, where previous analyses have detected a large clustering signal. After appropriate data selection, the TGSS sample we use contains ~110,000 sources selected at 150 MHz over ~70% of the sky. The survey footprint is largely superimposed on that of the NRAO VLA Sky Survey (NVSS) with the majority of TGSS sources having a counterpart in the NVSS sample. These characteristics make TGSS suitable for large-scale clustering analyses and facilitate the comparison with the results of previous studies. In this analysis we focus on the angular power spectrum, although the angular correlation function is also computed to quantify the contribution of multiple-component radio sources. We find that on large angular scales, corresponding to multipoles $2 \leq \ell \leq 30$, the amplitude of the TGSS angular power spectrum is significantly larger than that of the NVSS. We do not identify any observational systematic effects that may explain this mismatch. We have produced a number of physically motivated models for the TGSS angular power spectrum and found that all of them fail to match observations, even when taking into account observational and theoretical uncertainties. The same models provide a good fit to the angular spectrum of the NVSS sources. These results confirm the anomalous nature of the TGSS large-scale power, which has no obvious physical origin and seems to indicate that unknown systematic errors are present in the TGSS dataset.

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.

We present BAM: a novel Bias Assignment Method envisaged to generate mock catalogs. Combining the statistics of dark matter tracers from a high resolution cosmological $N$-body simulation and the dark matter density field calculated from down-sampled initial conditions using efficient structure formation solvers, we extract the halo-bias relation on a mesh of a $3\,h^{-1}$ Mpc cell side resolution as a function of properties of the dark matter density field (e.g. local density, cosmic web type), automatically including stochastic, deterministic, local and non-local components. We use this information to sample the halo density field, accounting for ignored dependencies through an iterative process. By construction, our approach reaches $\sim 1\%$ accuracy in the majority of the $k$-range up to the Nyquist frequency without systematic deviations for power spectra (about $k \sim 1\, h$ Mpc$^{-1}$) using either particle mesh or Lagrangian perturbation theory based solvers. When using phase-space mapping to compensate the low resolution of the approximate gravity solvers, our method reproduces the bispectra of the reference within $10\%$ precision studying configurations tracing the quasi-nonlinear regime. BAM has the potential to become a standard technique to produce mock halo and galaxy catalogs for future galaxy surveys and cosmological studies being highly accurate, efficient and parameter free.

We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method.

Using the almost all-sky 2MASS Photometric Redshift catalogue (2MPZ) we perform for the first time a tomographic analysis of galaxy angular clustering in the local Universe ($z<0.24$). We estimate the angular auto- and cross-power spectra of 2MPZ galaxies in three photometric redshift bins, and use dedicated mock catalogues to assess their errors. We measure a subset of cosmological parameters, having fixed the others at their Planck values, namely the baryon fraction $f_{b}=0.14^{+0.09}_{-0.06}$, the total matter density parameter $\Omega_{\rm m}=0.30\pm 0.06$, and the effective linear bias of 2MPZ galaxies $b_{\rm eff}$, which grows from $1.1^{+0.3}_{-0.4}$ at $\langle z \rangle=0.05$ up to $2.1^{+0.3}_{-0.5}$ at $\langle z\rangle=0.2$, largely because of the flux-limited nature of the dataset. The results obtained here for the local Universe agree with those derived with the same methodology at higher redshifts, and confirm the importance of the tomographic technique for next-generation photometric surveys such as Euclid or LSST.

We aim to quantify the capability of the inhomogeneous distribution of galaxy clusters, represented by the two-point statistics in Fourier space, to retrieve information on the underlying scaling relations. We make a case study using the mass X-ray luminosity scaling relation for galaxy clusters and study its impact on the clustering pattern of these objects. We define the luminosity-weighted power spectrum and introduce the luminosity power spectrum as direct assessment of the clustering of the property of interest, in our case, the cluster X-ray luminosity. Using a suite of halo catalogs extracted from N-body simulations and realistic estimates of the mass X-ray luminosity relation, we measured these statistics with their corresponding covariance matrices. By carrying out a Fisher matrix analysis, we quantified the content of information (by means of a figure-of merit) encoded in the amplitude, shape, and full shape of our probes for two-point statistics. The full shape of the luminosity power spectrum, when analyzed up to scales of k~0.2 h/Mpc, yields a figure of merit which is two orders of magnitude above the figure obtained from the unweighted power spectrum, and only one order of magnitude below the value encoded in X-ray luminosity function estimated from the same sample. This is a significant improvement over the analysis developed with the standard (i.e., unweighted) clustering probes. The measurements of the clustering of galaxy clusters and its explicit dependence on the cluster intrinsic properties can contribute to improving the degree of knowledge regarding the underlying links between cluster observables and the cluster masses (Abridged).

We describe the construction of a suite of galaxy cluster mock catalogues from N-body simulations, based on the properties of the new ROSAT-ESO Flux-Limited X-Ray (REFLEX II) galaxy cluster catalogue. Our procedure is based on the measurements of the cluster abundance, and involves the calibration of the underlying scaling relation linking the mass of dark matter haloes to the cluster X-ray luminosity determined in the \emph{ROSAT} energy band $0.1-2.4$ keV. In order to reproduce the observed abundance in the luminosity range probed by the REFLEX II X-ray luminosity function ($0.01<L_{X}/(10^{44}{\rm erg}\,{\rm s}^{-1}h^{-2})<10$), a mass-X ray luminosity relation deviating from a simple power law is required. We discuss the dependence of the calibration of this scaling relation on the X-ray luminosity and the definition of halo masses and analyse the one- and two-point statistical properties of the mock catalogues. Our set of mock catalogues provides samples with self-calibrated scaling relations of galaxy clusters together with inherent properties of flux-limited surveys. This makes them a useful tool to explore different systematic effects and statistical methods involved in constraining both astrophysical and cosmological information from present and future galaxy cluster surveys.