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 large amount of cosmological data already available (and in the near future) makes necessary the development of efficient numerical codes. Many software products have been implemented to perform cosmological analyses considering one or few probes. The need of multi-task software is rapidly increasing, in order to combine numerous cosmological probes along with their specificity (e.g., astrophysical descriptions and systematic errors). In this work we mention some of these libraries, bringing out some challenges they will face in the few-percent error era (on the cosmological parameters). We review some concepts of the standard cosmological model, and examine some specific topics on their implementation, bringing, for example, the discussion on how some quantities are numerically defined in different codes. We also consider implementation differences between public codes, mentioning their advantages/disadvantages.

The attempt to describe the recent accelerated expansion of the Universe includes different propositions for dark energy models and modified gravity theories. Establish their features in order to discriminate and even rule out part of these models using observational data is a fundamental issue of cosmology. In the present work we consider a class of extended theories of gravity (ETG) that are minimally coupled to the ordinary matter fields. In this context, and assuming a homogeneous and isotropic spacetime, we derive the energy conditions for this ETG class. We then put constraints on these conditions using a model-independent approach to reconstruct the deceleration function along with the Joint Light-curve Analysis (JLA) supernova sample, 11 baryon acoustic oscillation and 22 cosmic-chronometer measurements. We also consider an additional condition imposing the strong energy condition only on the ordinary matter. This is to guarantee the presence of only attractive matter in the energy-momentum tensor, at least in the redshift range of the observations, i.e., the recent accelerated expansion of the Universe is due solely to the modifications in the gravity theory. The main result of this work is a general reconstruction of the energy conditions valid for every considered ETG.

The standard model of cosmology, {\Lambda}CDM, is the simplest model that matches the current observations, but it relies on two hypothetical components, to wit, dark matter and dark energy. Future galaxy surveys and cosmic microwave background (CMB) experiments will independently shed light on these components, but a joint analysis that includes cross-correlations will be necessary to extract as much information as possible from the observations. In this paper, we carry out a multi-probe analysis based on pseudo-spectra and test it on publicly available data sets. We use CMB temperature anisotropies and CMB lensing observations from Planck as well as the spectroscopic galaxy and quasar samples of SDSS-III/BOSS, taking advantage of the large areas covered by these surveys. We build a likelihood to simultaneously analyse the auto and cross spectra of CMB lensing and tracer overdensity maps before running Monte-Carlo Markov Chains (MCMC) to assess the constraining power of the combined analysis. We then add the CMB temperature anisotropies likelihood and obtain constraints on cosmological parameters ($H_0$, $\omega_b$, $\omega_c$, ${\ln10^{10}A_s}$, $n_s$ and $z_{re}$) and galaxy biases. We demonstrate that the joint analysis can additionally constrain the total mass of neutrinos ${\Sigma m_{\nu}}$ as well as the dark energy equation of state $w$ at once (for a total of eight cosmological parameters), which is impossible with either of the data sets considered separately. Finally, we discuss limitations of the analysis related to, e.g., the theoretical precision of the models, particularly in the non-linear regime.

Among the various possibilities to probe the theory behind the recent accelerated expansion of the universe, the energy conditions (ECs) are of particular interest, since it is possible to confront and constrain the many models, including different theories of gravity, with observational data. In this context, we use the ECs to probe any alternative theory whose extra term acts as a cosmological constant. For this purpose, we apply a model-independent approach to reconstruct the recent expansion of the universe. Using Type Ia supernova, baryon acoustic oscillations and cosmic-chronometer data, we perform a Markov Chain Monte Carlo analysis to put constraints on the effective cosmological constant $\Omega^0_{\rm eff}$. By imposing that the cosmological constant is the only component that possibly violates the ECs, we derive lower and upper bounds for its value. For instance, we obtain that $0.59 < \Omega^0_{\rm eff} < 0.91$ and $0.40 < \Omega^0_{\rm eff} < 0.93$ within, respectively, $1\sigma$ and $3\sigma$ confidence levels. In addition, about 30\% of the posterior distribution is incompatible with a cosmological constant, showing that this method can potentially rule it out as a mechanism for the accelerated expansion. We also study the consequence of these constraints for two particular formulations of the bimetric massive gravity. Namely, we consider the Visser's theory and the Hassan and Roses's massive gravity by choosing a background metric such that both theories mimic General Relativity with a cosmological constant. Using the $\Omega^0_{\rm eff}$ observational bounds along with the upper bounds on the graviton mass we obtain constraints on the parameter spaces of both theories.

Approximate Bayesian Computation (ABC) enables parameter inference for complex physical systems in cases where the true likelihood function is unknown, unavailable, or computationally too expensive. It relies on the forward simulation of mock data and comparison between observed and synthetic catalogues. Here we present cosmoabc, a Python ABC sampler featuring a Population Monte Carlo (PMC) variation of the original ABC algorithm, which uses an adaptive importance sampling scheme. The code is very flexible and can be easily coupled to an external simulator, while allowing to incorporate arbitrary distance and prior functions. As an example of practical application, we coupled cosmoabc with the numcosmo library and demonstrate how it can be used to estimate posterior probability distributions over cosmological parameters based on measurements of galaxy clusters number counts without computing the likelihood function. cosmoabc is published under the GPLv3 license on PyPI and GitHub and documentation is available at http://goo.gl/SmB8EX

We determine the mass scale of Planck galaxy clusters using gravitational lensing mass measurements from the Cluster Lensing And Supernova survey with Hubble (CLASH). We have compared the lensing masses to the Planck Sunyaev-Zeldovich (SZ) mass proxy for 21 clusters in common, employing a Bayesian analysis to simultaneously fit an idealized CLASH selection function and the distribution between the measured observables and true cluster mass. We used a tiered analysis strategy to explicitly demonstrate the importance of priors on weak lensing mass accuracy. In the case of an assumed constant bias, $b_{SZ}$, between true cluster mass, $M_{500}$, and the Planck mass proxy, $M_{PL}$, our analysis constrains $1- b_{SZ} = 0.73\pm 0.10$ when moderate priors on weak lensing accuracy are used, including a zero-mean Gaussian with standard deviation of 8% to account for possible bias in lensing mass estimations. Our analysis explicitly accounts for possible selection bias effects in this calibration sourced by the CLASH selection function. Our constraint on the cluster mass scale is consistent with recent results from the Weighing the Giants program and the Canadian Cluster Comparison Project. It is also consistent, at 1.34$\sigma$, with the value needed to reconcile the Planck SZ cluster counts with Planck's base $\Lambda$CDM model fit to the primary cosmic microwave background anisotropies.

Evidences for late-time acceleration of the Universe are provided by multiple probes, such as Type Ia supernovae, the cosmic microwave background (CMB) and large-scale structure (LSS). In this work, we focus on the integrated Sachs--Wolfe (ISW) effect, i.e., secondary CMB fluctuations generated by evolving gravitational potentials due to the transition between, e.g., the matter and dark energy (DE) dominated phases. Therefore, assuming a flat universe, DE properties can be inferred from ISW detections. We present a Bayesian approach to compute the CMB--LSS cross-correlation signal. The method is based on the estimate of the likelihood for measuring a combined set consisting of a CMB temperature and a galaxy contrast maps, provided that we have some information on the statistical properties of the fluctuations affecting these maps. The likelihood is estimated by a sampling algorithm, therefore avoiding the computationally demanding techniques of direct evaluation in either pixel or harmonic space. As local tracers of the matter distribution at large scales, we used the Two Micron All Sky Survey (2MASS) galaxy catalog and, for the CMB temperature fluctuations, the ninth-year data release of the Wilkinson Microwave Anisotropy Probe (WMAP9). The results show a dominance of cosmic variance over the weak recovered signal, due mainly to the shallowness of the catalog used, with systematics associated with the sampling algorithm playing a secondary role as sources of uncertainty. When combined with other complementary probes, the method presented in this paper is expected to be a useful tool to late-time acceleration studies in cosmology.

Distance measurements are currently the most powerful tool to study the expansion history of the universe without specifying its matter content nor any theory of gravitation. Assuming only an isotropic, homogeneous and flat universe, in this work we introduce a model-independent method to reconstruct directly the deceleration function via a piecewise function. Including a penalty factor, we are able to vary continuously the complexity of the deceleration function from a linear case to an arbitrary $(n+1)$-knots spline interpolation. We carry out a Monte Carlo (MC) analysis to determine the best penalty factor, evaluating the bias-variance trade-off, given the uncertainties of the SDSS-II and SNLS supernova combined sample (JLA), compilations of baryon acoustic oscillation (BAO) and $H(z)$ data. The bias-variance analysis is done for three fiducial models with different features in the deceleration curve. For each fiducial model, we test different reconstructions using, in each case, more than $10^4$ catalogs in a total of about $5\times 10^5$. This investigation proved to be essential in determining the best reconstruction to study these data. We show that, evaluating a single fiducial model, the conclusions about the bias-variance ratio are misleading. We determine the reconstruction method in which the bias represents at most $10\%$ of the total uncertainty. Finally, we apply the Ensemble Sampler Markov Chain Monte Carlo (ESMCMC) method to explore the posterior of the deceleration function up to redshift $1.3$ (using only JLA) and $2.3$ (JLA+BAO+$H(z)$). We obtain that the standard cosmological model agrees within $3\sigma$ level with the reconstructed results in the whole studied redshift intervals. Since our method is calibrated to minimize the bias, the error bars of the reconstructed functions are a good approximation for the total uncertainty.

Sunyaev-Zel'dovich (SZ) surveys are promising probes of cosmology - in particular for Dark Energy (DE) -, given their ability to find distant clusters and provide estimates for their mass. However, current SZ catalogs contain tens to hundreds of objects and maximum likelihood estimators may present biases for such sample sizes. In this work we use the Monte Carlo approach to determine the presence of bias on cosmological parameter estimators from cluster abundance as a function of the area and depth of the survey, and the number of cosmological parameters fitted. Assuming perfect knowledge of mass and redshift some estimators have non-negligible biases. For example, the bias of $\sigma_8$ corresponds to about $40%$ of its statistical error bar when fitted together with $\Omega_c$ and $w_0$. Including a SZ mass-observable relation decreases the relevance of the bias, for the typical sizes of current surveys. The biases become negligible when combining the SZ data with other cosmological probes. However, we show that the biases from SZ estimators do not go away with increasing sample sizes and they may become the dominant source of error for an all sky survey at the South Pole Telescope (SPT) sensitivity. The results of this work validate the use of the current maximum likelihood methods for present SZ surveys, but highlight the need for further studies for upcoming experiments. [abridged]