We propose a new method, called MonteCarlo Posterior Fit, to boost the MonteCarlo sampling of likelihood (posterior) functions. The idea is to approximate the posterior function by an analytical multidimensional non-Gaussian fit. The many free parameters of this fit can be obtained by a smaller sampling than is needed to derive the full numerical posterior. In the examples that we consider, based on supernovae and cosmic microwave background data, we find that one needs an order of magnitude smaller sampling than in the standard algorithms to achieve comparable precision. This method can be applied to a variety of situations and is expected to significantly improve the performance of the MonteCarlo routines in all the cases in which sampling is very time-consuming. Finally, it can also be applied to Fisher matrix forecasts, and can help solve various limitations of the standard approach.

In this work we provide updated constraints on coupled dark energy (CDE) cosmology with Peebles-Ratra (PR) potential and constant coupling strength $\beta$. This modified gravity scenario introduces a fifth force between dark matter particles, mediated by a scalar field that plays the role of dark energy. The mass of the dark matter particles does not remain constant, but changes with time as a function of the scalar field. Here we focus on the phenomenological behavior of the model, and assess its ability to describe updated cosmological data sets that include the Planck 2018 cosmic microwave background (CMB) temperature, polarization and lensing, baryon acoustic oscillations, the Pantheon compilation of supernovae of Type Ia, data on $H(z)$ from cosmic chronometers, and redshift-space distortions. We also study which is the impact of the local measurement of $H_0$ from SH0ES and the strong-lensing time delay data from the H0LICOW collaboration on the parameter that controls the strength of the interaction in the dark sector. We find a peak corresponding to a coupling $\beta > 0$ and to a potential parameter $\alpha > 0$, more or less evident depending on the data set combination. We show separately the impact of each data set and remark that it is especially CMB lensing the one data set that shifts the peak the most towards $\Lambda$CDM. When a model selection criterion based on the full Bayesian evidence is applied, however, $\Lambda$CDM is still preferred in all cases, due to the additional parameters introduced in the CDE model.

We reconstruct in this paper the deceleration and jerk parameters as functions of the cosmological redshift from data on cosmic chronometers (CCH), baryon acoustic oscillations (BAOs), and the Pantheon+MCT compilation of supernovae of Type Ia (SnIa). The reconstruction is carried out with the Weighted Function Regression method, previously introduced by G\'omez-Valent & Amendola (2018). It improves the usual cosmographic approach by automatically implementing Occam's razor criterion. This makes our procedure to be more free of model and parametrization dependencies than many other analyses in the literature. The reconstructed functions are fully compatible with the predictions for the concordance model. In addition, we also discuss the confidence level at which we can claim that the Universe (assumed to be flat, homogeneous and isotropic) is currently accelerating. According to Jeffreys' scale and jargon, we find moderate evidence in favor of such speed-up using the data on SnIa+CCH, and very strong one when we also use data on BAOs. The measured current value of the deceleration parameter in the latter case reads $q_0\sim -0.60\pm 0.10$, and for the deceleration-acceleration transition redshift we find $z_t\sim 0.8\pm 0.10$. The former is $\sim 6\sigma$ away from $0$. This is in stark contrast, for instance, with the $\sim 17\sigma$ that are found in the context of the flat $\Lambda$CDM even without including the BAOs data. This indicates that cosmography and Occam's razor criterion play a crucial role in this discussion, and that estimating the evidence for positive acceleration only in the framework of a particular cosmological model or parametrization is clearly insufficient.

We derive new constraints on the Hubble parameter $H_0$ using the available data on $H(z)$ from cosmic chronometers (CCH), and the Hubble rate data points from the supernovae of Type Ia (SnIa) of the Pantheon compilation and the Hubble Space Telescope (HST) CANDELS and CLASH Multy-Cycle Treasury (MCT) programs. We employ two alternative techniques, Gaussian Processes (GPs) and the Weighted Polynomial Regression (WPR) method, to reconstruct the Hubble function, determine the derived values of $H_0$, and compare them with the local HST measurement provided by Riess {\it et al.} (2018), $H_0^{{\rm HST}}=(73.48\pm1.66)$ km/s/Mpc, and with the Planck+$\Lambda$CDM value, $H_0^{{\rm P18}}=(66.88\pm 0.92)$ km/s/Mpc. With GPs we obtain $H_0=(67.99\pm 1.94)$ km/s/Mpc, and with the WPR method $H_0=(68.90\pm 1.96)$ km/s/Mpc. Both are fully compatible at $<1\sigma$ c.l., and also with $H_0^{{\rm P18}}$. In contrast, they are in $\sim 2\sigma$ tension with $H_0^{{\rm HST}}$.

Recent studies suggest that dynamical dark energy (DDE) provides a better fit to the rising affluence of modern cosmological observations than the concordance model ($\Lambda$CDM) with a rigid cosmological constant, $\Lambda$. Such is the case with the running vacuum models (RVMs) and to some extent also with a simple XCDM parametrization. Apart from the cosmic microwave background (CMB) anisotropies, the most crucial datasets potentially carrying the DDE signature are: i) baryonic acoustic oscillations (BAO), and ii) direct large scale structure (LSS) formation data (i.e. the observations on $f(z)\sigma_8(z)$ at different redshifts). As it turns out, analyses mainly focusing on CMB and with insufficient BAO+LSS input, or those just making use of gravitational weak-lensing data for the description of structure formation, generally fail to capture the DDE signature, whereas the few existing studies using a rich set of CMB+BAO+LSS data (see in particular Sol\`a, G\'omez-Valent & de Cruz P\'erez 2015,2017; and Zhao et al. 2017) do converge to the remarkable conclusion that DDE might well be encoded in the current cosmological observations. Being the issue so pressing, here we explain both analytically and numerically the origin of the possible hints of DDE in the context of RVMs, which arise at a significance level of $3-4\sigma$. By performing a detailed study on the matter and vacuum perturbations within the RVMs, and comparing with the XCDM, we show why the running vacuum fully relaxes the existing $\sigma_8$-tension and accounts for the LSS formation data much better than the concordance model.

In this paper we present new constraints on the Hubble parameter $H_0$ using: (i) the available data on $H(z)$ obtained from cosmic chronometers (CCH); (ii) the Hubble rate data points extracted from the supernovae of Type Ia (SnIa) of the Pantheon compilation and the Hubble Space Telescope (HST) CANDELS and CLASH Multy-Cycle Treasury (MCT) programs; and (iii) the local HST measurement of $H_0$ provided by Riess et al. (2018), $H_0^{{\rm HST}}=(73.45\pm1.66)$ km/s/Mpc. Various determinations of $H_0$ using the Gaussian processes (GPs) method and the most updated list of CCH data have been recently provided by Yu, Ratra and Wang (2018). Using the Gaussian kernel they find $H_0=(67.42\pm 4.75)$ km/s/Mpc. Here we extend their analysis to also include the most released and complete set of SnIa data, which allows us to reduce the uncertainty by a factor $\sim 3$ with respect to the result found by only considering the CCH information. We obtain $H_0=(67.06\pm 1.68)$ km/s/Mpc, which favors again the lower range of values for $H_0$ and is in tension with $H_0^{{\rm HST}}$. The tension reaches the $2.71\sigma$ level. We round off the GPs determination too by taking also into account the error propagation of the kernel hyperparameters when the CCH with and without $H_0^{{\rm HST}}$ are used in the analysis. In addition, we present a novel method to reconstruct functions from data, which consists in a weighted sum of polynomial regressions (WPR). We apply it from a cosmographic perspective to reconstruct $H(z)$ and estimate $H_0$ from CCH and SnIa measurements. The result obtained with this method, $H_0=(68.90\pm 1.96)$ km/s/Mpc, is fully compatible with the GPs ones. Finally, a more conservative GPs+WPR value is also provided, $H_0=(68.45\pm 2.00)$ km/s/Mpc, which is still almost $2\sigma$ away from $H_0^{{\rm HST}}$.