In this work we utilise the most recent publicly available type Ia supernova (SN Ia) compilations and implement a well formulated cosmological model based on Lema\^{i}tre-Tolman-Bondi metric in presence of cosmological constant $\Lambda$ ($\Lambda$LTB) to test for signatures of large local inhomogeneities at $z\leq0.15$. Local underdensities in this redshift range have been previously found based on luminosity density (LD) data and galaxy number counts. Our main constraints on the possible local void using the Pantheon SN Ia dataset are: redshift size of $z_{\rm size}=0.068^{+0.021}_{-0.030}$; density contrast of $\delta\Omega_0/\Omega_0=-10.5_{-7.4}^{+9.3}\%$ between 16th and 84th percentiles. Investigating the possibility to alleviate the $\sim9\%$ disagreement between measurements of present expansion rate $H_0$ coming from calibrated local SN Ia and high-$z$ cosmic microwave background data, we find large local void to be a very unlikely explanation alone, consistently with previous studies. However, the level of matter inhomogeneity at a scale of $\sim$100Mpc that is allowed by SN Ia data, although not expected from cosmic variance calculations in standard model of cosmology, could be the origin of additonal systematic error in distance ladder measurements based on SN Ia. Fitting low-redshift Pantheon data with a cut $0.023<z<0.15$ to the $\Lambda$LTB model and to the Taylor expanded luminosity distance formula we estimate that this systematic error amounts to $1.1\%$ towards the lower $H_0$ value. A test for local anisotropy in Pantheon SN Ia data yields null evidence. Analysis of LD data provides a constraint on contrast of large isotropic void $\delta\Omega_0/\Omega_0=-51.9\%\pm6.3\%$, which is in $\sim4\sigma$ tension with SN Ia results. More data are necessary to better constrain the local matter density profile and understand the disagreement between SN and LD samples

The quest for understanding the late-time acceleration is haunted by an immense freedom in the analysis of dynamical models for dark energy in extended parameter spaces. Often-times having no prior knowledge at our disposal, arbitrary choices are implemented to reduce the degeneracies between parameters. We also encounter this issue in the case of quintessence fields, where a scalar degree of freedom drives the late-time acceleration. In this study, we implement a more physical prescription, the \textit{flow} condition, to fine-tune the quintessence evolution for several field potentials. We find that this prescription agrees well with the most recent catalogue of data, namely supernovae type Ia, baryon acoustic oscillations, cosmic clocks and distance to last scattering surface, and it enables us to infer the initial conditions for the field, both potential and cosmological parameters. At $2\sigma$ we find stricter bounds on the potential parameters $f/m_{pl}>0.26$ and $n<0.15$ for the PNGB and IPL potentials, respectively, while constraints on cosmological parameters remain extremely consistent across all assumed potentials. By implementing information criteria to assess their ability to fit the data, we do not find any evidence against thawing models, which in fact are statistically equivalent to $\Lambda$CDM, and the freezing ones are moderately disfavoured. Through our analysis we place upper bounds on the slope of quintessence potentials, consequently revealing a strong tension with the recently proposed swampland criterion, finding the 2$\sigma$ upper bound of $\lambda \sim 0.31$ for the exponential potential.

In the current work, we have implemented an extension of the standard Gaussian Process formalism, namely the Multi-Task Gaussian Process with the ability to perform a joint learning of several cosmological data simultaneously. We have utilised the "low-redshift" expansion rate data from Supernovae Type-Ia (SN), Baryon Acoustic Oscillations (BAO) and Cosmic Chronometers (CC) data in a joint analysis. We have tested several possible models of covariance functions and find very consistent estimates for cosmologically relevant parameters. In the current formalism, we also find provisions for heuristic arguments which allow us to select the best-suited kernel for the reconstruction of expansion rate data. We also utilised our method to account for systematics in CC data and find an estimate of $H_0 = 68.52^{+0.94 + 2.51 (sys)}_{-0.94} $ $\textrm{km/s Mpc}^{-1}$ and a corresponding $r_d = 145.61^{+2.82}_{ - 2.82 - 4.3 (sys)} $ Mpc as our primary result. Subsequently, we find constraints on the present deceleration parameter $q_0 = -0.52 \pm 0.06$ and the transition redshift $z_T = 0.64^{+0.12}_{-0.09}$. All the estimated cosmological parameters are found to be in good agreement with the standard $\Lambda$CDM scenario. Including the local model-independent $H_0$ estimate to the analysis we find $H_0 = 71.40^{ + 0.30 + 1.65 (sys)}_{- 0.30 } $ $\textrm{km/s Mpc}^{-1}$ and the corresponding $r_d = 141.29^{ + 1.31 }_{-1.31-2.63 (sys)}$ Mpc. Also, the constraints on $r_d H_0$ remain consistent throughout our analysis and also with the model-dependent CMB estimate. Using the $\mathcal{O}m(z)$ diagnostic, we find that the concordance model is very consistent within the redshift range $z \lesssim 2$ and mildly discrepant for $z \gtrsim 2$.

We conduct a selective analysis of the isotropic ($D_V$) and anisotropic ($AP$) components of the most recent Baryon Acoustic Oscillations (BAO) data. We find that these components provide significantly different constraints and could provide strong diagnostics for model selection, also in view of more precise data to arrive. For instance, in the $\Lambda$CDM model, we find a mild tension of $\sim 2 \sigma$ for the $\Omega_m$ estimates obtained using $D_V$ and $AP$ separately. Considering both $\Omega_k$ and $w$ as free parameters, we find that the concordance model is in tension with the best-fit values provided by the BAO data alone at 2.2$\sigma$. We complemented the BAO data with the Supernova Ia (SNIa) and Observational \textit{Hubble} datasets to perform a joint analysis on the $\Lambda$CDM model and its standard extensions. By assuming $\Lambda$CDM scenario, we find that these data provide $H_0 = 69.4 \pm 1.7$ \text{km/s Mpc$^{-1} $} as the best-fit value for the present expansion rate. In the $k\Lambda$CDM scenario we find that the evidence for acceleration using the BAO data alone is more than $\sim 5.8\sigma$, which increases to $8.4 \sigma$ in our joint analysis.

In this paper we summarise the constraints that low-redshift data --such as supernovae Ia (SN Ia), baryon acoustic oscillations (BAO) and cosmic chronometers (CC)-- are able to set on the concordance model and its extensions, as well as on inhomogeneous but isotropic models. We provide a broad overlook into these cosmological scenarios and several aspects of data analysis. In particular, we review a number of systematic issues of SN Ia analysis that include magnitude correction techniques, selection bias and their influence on the inferred cosmological constraints. Furthermore, we examine the isotropic and anisotropic components of the BAO data and their individual relevance for cosmological model-fitting. We extend the discussion presented in earlier works regarding the inferred dynamics of cosmic expansion and its present rate from the low-redshift data. Specifically, we discuss the cosmological constraints on the accelerated expansion and related model-selections. In addition, we extensively talk about the Hubble constant problem, then focus on the low-redshift data constraint on $H_0$ that is based on CC. Finally, we present the way in which this result compares the two of high-redshift $H_0$ estimate and local (redshift zero) measurements that are in tension.

A recent analysis of the Supernova Ia data claims a 'marginal' ($\sim3\sigma$) evidence for a cosmic acceleration. This result has been complemented with a non-accelerating $R_{h}=ct$ cosmology, which was presented as a valid alternative to the $\Lambda$CDM model. In this paper, we use the same analysis to show that a non-marginal evidence for acceleration is actually found. We compare the standard Friedmann models to the $R_{h}=ct$ cosmology by complementing SN Ia data with the Baryon Acoustic Oscillations, Gamma Ray Bursts and Observational Hubble datasets. We also study the power-law model which is a functional generalisation of $R_{h}=ct$. We find that the evidence for late-time acceleration is beyond refutable at a 4.56$\sigma$ confidence level from SN Ia data alone, and at an even stronger confidence level ($5.38\sigma$) from our joint analysis. Also, the non-accelerating $R_{h}=ct$ model fails to statistically compare with the $\Lambda$CDM having a $\Delta(\text{AIC})\sim30$.

We test the theoretical predictions of several cosmological models against different observables to compare the indirect estimates of the current expansion rate of the Universe determined from model fitting with the direct measurements based on Cepheids data published recently. We perform a statistical analysis of type Ia supernova (SN Ia), Hubble parameter, and baryon acoustic oscillation data. A joint analysis of these datasets allows us to better constrain cosmological parameters, but also to break the degeneracy that appears in the distance modulus definition between $H_0$ and the absolute B-band magnitude of SN Ia, $M_0$. From the theoretical side, we considered spatially flat and curvature-free $\Lambda$CDM, $w$CDM, and inhomogeneous Lema\^{i}tre-Tolman-Bondi (LTB) models. To analyse SN Ia we took into account the distributions of SN Ia intrinsic parameters. For the $\Lambda$CDM model we find that $\Omega_m=0.35\pm0.02$, $H_0=(67.8\pm1.0)\,$km$\,$s$^{-1}/$Mpc, while the corrected SN absolute magnitude has a normal distribution ${\cal N}(19.13,0.11)$. The $w$CDM model provides the same value for $\Omega_m$, while $H_0=(66.5\pm1.8)\,$km$\,$s$^{-1}/$Mpc and $w=-0.93\pm0.07$. When an inhomogeneous LTB model is considered, the combined fit provides $H_0=(64.2\pm1.9)\,$km$\,$s$^{-1}/$Mpc. Both the Akaike information criterion and the Bayes factor analysis cannot clearly distinguish between $\Lambda$CDM and $w$CDM cosmologies, while they clearly disfavour the LTB model. For the $\Lambda$CDM, our joint analysis of the SN Ia, the Hubble parameter, and the baryon acoustic oscillation datasets provides $H_0$ values that are consistent with CMB-only Planck measurements, but they differ by $2.5\sigma$ from the value based on Cepheids data.