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 generation of simulated convergence maps is of key importance in fully exploiting weak lensing by Large Scale Structure (LSS) from which cosmological parameters can be derived. In this paper we present an extension of the PINOCCHIO code which produces catalogues of dark matter haloes so that it is capable of simulating weak lensing by LSS. Like WL-MOKA, the method starts with a random realisation of cosmological initial conditions, creates a halo catalogue and projects it onto the past-light-cone, and paints in haloes assuming parametric models for the mass density distribution within them. Large scale modes that are not accounted for by the haloes are constructed using linear theory. We discuss the systematic errors affecting the convergence power spectra when Lagrangian Perturbation Theory at increasing order is used to displace the haloes within PINOCCHIO, and how they depend on the grid resolution. Our approximate method is shown to be very fast when compared to full ray-tracing simulations from an N-Body run and able to recover the weak lensing signal, at different redshifts, with a few percent accuracy. It also allows for quickly constructing weak lensing covariance matrices, complementing PINOCCHIO's ability of generating the cluster mass function and galaxy clustering covariances and thus paving the way for calculating cross covariances between the different probes. This work advances these approximate methods as tools for simulating and analysing surveys data for cosmological purposes.

Observational cosmology is passing through a unique moment of grandeur with the amount of quality data growing fast. However, in order to better take advantage of this moment, data analysis tools have to keep up the pace. Understanding the effect of baryonic matter on the large-scale structure is one of the challenges to be faced in cosmology. In this work, we have thoroughly studied the effect of baryonic physics on different lensing statistics. Making use of the Magneticum Pathfinder suite of simulations we show that the influence of luminous matter on the 1-point lensing statistics of point sources is significant, enhancing the probability of magnified objects with $\mu>3$ by a factor of $2$ and the occurrence of multiple-images by a factor $5-500$ depending on the source redshift and size. We also discuss the dependence of the lensing statistics on the angular resolution of sources. Our results and methodology were carefully tested in order to guarantee that our uncertainties are much smaller than the effects here presented.

The abundances of dark matter halos in the universe are described by the halo mass function (HMF). It enters most cosmological analyses and parametrizes how the linear growth of primordial perturbations is connected to these abundances. Interestingly, this connection can be made approximately cosmology independent. This made it possible to map in detail its near-universal behavior through large-scale simulations. However, such simulations may suffer from systematic effects, especially if baryonic physics is included. In this paper we ask how well observations can constrain directly the HMF. The observables we consider are galaxy cluster number counts, galaxy cluster power spectrum and lensing of type Ia supernovae. Our results show that DES is capable of putting the first meaningful constraints on the HMF, while both Euclid and J-PAS can give stronger constraints, comparable to the ones from state-of-the-art simulations. We also find that an independent measurement of cluster masses is even more important for measuring the HMF than for constraining the cosmological parameters, and can vastly improve the determination of the halo mass function. Measuring the HMF could thus be used to cross-check simulations and their implementation of baryon physics. It could even, if deviations cannot be accounted for, hint at new physics.

The supernova (SN) Hubble diagram residual contains valuable information on both the present matter power spectrum and its growth history. In this paper we show that this information can be retrieved with precision by combining both peculiar velocity and weak-lensing analysis on the data. To wit, peculiar velocity induces correlations on the nearby SN while lensing induces a non-Gaussian dispersion in faraway objects. We show that both effects have almost orthogonal degeneracies and discuss how they can be extracted simultaneously from the data. We analyze the JLA supernova catalog in a 14-dimensional parameter space, assuming a flexible growth-rate index $\gamma$. We arrive at the following marginalized constraints: $\sigma_8 = 0.65^{+0.23}_{-0.37}$ and $\gamma = 1.38^{+1.7}_{-0.65}$. Assuming instead GR as the correct gravitation theory (and thus $\gamma \equiv 0.55$), the constraints in $\sigma_8$ tighten further: $\sigma_8 = 0.40^{+0.21}_{-0.23}$. We show that these constraints complement well the ones obtained from other datasets and that they could improve substantially with more SNe.

A recently proposed technique allows one to constrain both the background and perturbation cosmological parameters through the distribution function of supernova Ia apparent magnitudes. Here we extend this technique to alternative cosmological scenarios, in which the growth of structure does not follow the $\Lambda$CDM prescription. We apply the method first to the supernova data provided by the JLA catalog combined with all the current independent redshift distortion data and with low-redshift cluster data from Chandra and show that although the supernovae alone are not very constraining, they help in reducing the confidence regions. Then we apply our method to future data from LSST and from a survey that approximates the Euclid satellite mission. In this case we show that the combined data are nicely complementary and can constrain the normalization $\sigma_8$ and the growth rate index $\gamma$ to within $0.6\%$ and $7\%$, respectively. In particular, the LSST supernova catalog is forecast to give the constraint $\gamma (\sigma_8/0.83)^{6.7} = 0.55 \pm 0.1$. We also report on constraints relative to a step-wise parametrization of the growth rate of structures. These results show that supernova lensing serves as a good cross-check on the measurement of perturbation parameters from more standard techniques.

A method was recently proposed which allows the conversion of the weak-lensing effects in the supernova Hubble diagram from noise into signal. Such signal is sensitive to the growth of structure in the universe, and in particular can be used as a measurement of $\sigma_8$ which is independent from more traditional methods such as those based on the CMB, cosmic shear or cluster abundance. We extend here that analysis to allow for intrinsic non-Gaussianities in the supernova PDF, and discuss how this can be best modelled using the Bayes Factor. Although it was shown that a precise measurement of $\sigma_8$ requires ~$10^5$ supernovae, current data already allows an important proof of principle. In particular we make use of the 732 supernovae with z < 1 of the recent JLA catalog and show that a simple treatment of intrinsic non-Gaussianities with a couple of nuisance parameters is enough for our method to yield the values $\sigma_8 = 0.84^{+0.28}_{-0.65}$ or $\sigma_8 < 1.45$ at a $2\sigma$ confidence level. This result is consistent with mock simulations and it is also in agreement with independent measurements and presents the first ever measurement of $\sigma_8$ using supernova magnitudes alone.