As highlighted in the Introduction, the Lambda Cold Dark Matter ($\mathsf{\Lambda }$CDM) model is the standard or concordance model of a spatially flat universe. In this model, dark energy is represented by the cosmological constant $\mathsf{\Lambda }$, which is assumed to be associated with the vacuum energy density of quantum fields, and its energy density is a constant here $\mathsf{\Lambda }=4.33·{10}^{-66}{\mathrm{eV}}^2$. The pressure and the energy density in the $\mathsf{\Lambda }$CDM model are related as leading to the constant EoS parameter

The action describing the cosmological constant $\mathsf{\Lambda }$ is [43] where $g\equiv \det \left(g_{\mu \nu }\right)$ is the determinant of the metric tensor $g_{\mu \nu }$, R is the Ricci scalar, and $S_{\mathrm{M}}$ is the action describing matter. The spatially flat $\mathsf{\Lambda }$CDM model is typically characterized by six independent parameters [159]: the physical baryon density parameter, ${\mathsf{\Omega }}_{\mathrm{b}}{h}^2$; the dark matter physical density parameter, ${\mathsf{\Omega }}_{\mathrm{c}}{h}^2$; the age of the universe, $t_0$; the scalar spectral index, $n_{\mathrm{s}}$; the amplitude of the curvature fluctuations, ${\Delta }_{\mathrm{R}}^2$; and the optical depth during the reionization period for $z\in (6,20)$, $\tau$. In addition to these parameters, the $\mathsf{\Lambda }$CDM model is described by six extended fixed parameters: the total density parameter, ${\mathsf{\Omega }}_{\mathrm{tot}}$; the EoS parameter, $w_{\mathsf{\Lambda }}$; the total mass of three types of neutrinos, $\sum m_{\nu }$; the effective number of the relativistic degrees of freedom, $N_{\mathrm{eff}}$; the tensor/scalar ratio, r; the scalar spectral index running, $a_s$.

The extension of the spatially flat $\mathsf{\Lambda }$CDM model to spatially non-flat hypersurfaces is the oCDM model. The first Friedmann equations describing the evolution of the universe for the spatially non-flat oCDM model and for the spatially flat $\mathsf{\Lambda }$CDM model are defined, respectively, as

, here ${\mathsf{\Omega }}_{\mathrm{r}0}={\rho }_{\mathrm{r}0}/{\rho }_{\mathrm{cr}}$, ${\mathsf{\Omega }}_{\mathrm{m}0}={\rho }_{\mathrm{m}0}/{\rho }_{\mathrm{cr}}$, ${\mathsf{\Omega }}_{\mathsf{\Lambda }}={\rho }_{\mathsf{\Lambda }}/{\rho }_{\mathrm{cr}}$ are density parameters at present epoch for radiation, matter and vacuum, respectively, with ${\rho }_{\mathrm{r}0}$ and ${\rho }_{\mathrm{m}0}$ are energy densities for radiation and matter at present epoch, respectively. The value of the critical energy density at present epoch is equal to ${\rho }_{\mathrm{cr}}=3M_{\mathrm{pl}}^2{H}_0^2/8\pi =1.8791h^2·{10}^{-29}{\mathrm{g}\mathrm{cm}}^{-3}$; ${\mathsf{\Omega }}_{\mathrm{k}0}=-\mathrm{k}/H_0^2$ is a spatial curvature density parameter at present epoch, ${\rho }_{\mathrm{k}0}=-3\mathrm{k}{M}_{\mathrm{pl}}^2/8\pi$ is a spatial curvature density, $\mathrm{k}$ is a curvature parameter, and $E\left(a\right)\equiv H\left(a\right)/H_0$ is a normalized Hubble parameter.

Observational constraints on cosmological parameters ${\mathsf{\Omega }}_{m0}$ and ${\mathsf{\Omega }}_{\mathsf{\Lambda }}$, obtained from different cosmological datasets for the $\mathsf{\Lambda }$CDM model and for the oCDM model, are represented in Figure 1.

As we mentioned above, the $\mathsf{\Lambda }$CDM model is the fiducial model against which all alternative models are compared regarding their fit to observational data. Its predictions agree with the observational data pertaining to the accelerated expansion of the universe, the statistical distribution of LSS, the CMB temperature and polarization anisotropies, and the abundance of light elements in the universe [76,160].

There are numerous physically motivated alternative models for the $\mathsf{\Lambda }$CDM model [43,56,164,165,166,167,168,169,170,171,172,173,174]. One of the prominent alternatives to the $\mathsf{\Lambda }$CDM model are dynamical scalar field $\varphi$CDM models [36,91,92]. In these models, the role of dark energy is played by a slowly varying uniform self-interacting scalar field $\varphi$, minimally interacting with gravity. These $\varphi$CDM models involving a dynamical scalar field do not suffer from the fine tuning problem of the $\mathsf{\Lambda }$CDM model, and have a more natural explanation for the observed low-energy scale of dark energy. In addition, the EoS parameter time dependent in the $\varphi$CDM model, as opposed to it being a constant $w=-1$ in the $\mathsf{\Lambda }$CDM model. When the energy density of the scalar field begins to dominate over the energy density of both radiation and matter, the universe begins the stage of the accelerated expansion. At early times during the evolution of the universe, the behavior of the dynamical scalar field is different from that of the cosmological constant $\mathsf{\Lambda }$, but is almost indistinguishable from that of the cosmological constant during later times.

Dynamical scalar field $\varphi$CDM models are divided into two classes: the quintessence models [112] and phantom models [102,175]. These two classes of models differ from each other by the following attributes:

The action describing a scalar field in the presence of matter is given by where ${\mathcal{L}}_{\varphi }$ is the Lagrangian density of the scalar field, the form of which depends on the type of the chosen model, and as mentioned above, R is the Ricci scalar, and $S_{\mathrm{M}}$ is the action describing matter. We describe the form of ${\mathcal{L}}_{\varphi }$ for the quintessence and the phantom fields below.

The action for quintessential inflation models has the general form where $S_{\mathrm{R}}$ is the action describing radiation, $S_{\mathrm{I}}$ is the action describing the interactions of the inflaton field with the fermion ($\psi$), scalar ($\chi$), vector ($B_{\mu }$) degrees of freedom in the Standard Model and beyond.

To maintain inflation over a long period of time, it is necessary that the acceleration caused by the inflaton field be sufficiently small compared to its velocity over the Hubble time. Under these conditions, the first Friedmann’s and Klein-Gordon’s equations for the inflaton field in the spatially flat universe take the form [122,187,188]

The slow-roll regime of the inflaton field is provided by potential of the inflaton field $V\left(\varphi \right)$ with certain shapes: exponential [131], power-law [128,129], and plateau-like [132,133]. The slow-roll parameters, which determine the curvature and slope of the potential of the inflaton field, should remain small for some period of time to sustain the inflationary behavior The condition $ϵ\ll 1$ is necessary for the accelerated expansion of the universe to take place, while the condition $ϵ\simeq 1$ defines the end of inflation. In the case of $ϵ\approx 0$, the expansion rate of the universe is approximately constant, i.e., the scale factor grows exponentially with time, $a\propto e^{Ht}$, resulting in an almost de Sitter expansion of the universe.

The scalar spectral index ($n_s$), tensor spectral index ($n_t$), scalar spectral index running ($a_s$), and tensor-to-scalar ratio (r) are defined, respectively, as [131]

During the inflationary epoch of the universe, scalar and tensor perturbations are created from quantum vacuum fluctuations and are spatially stretched to superhorizon scales, where they become classical, and the almost scale-invariant tilted primordial power spectrum is formed [189]. The tilted primordial scalar ${\mathcal{P}}_s\left(k\right)$ and tensor ${\mathcal{P}}_t\left(k\right)$ power spectra for spatially flat tilted quintessential inflation models are defined in terms of the wave number k as [122,187,188] where $A_s$ and $A_t$ are the curvature perturbations amplitude and tensor amplitude at the pivot scale $k_0=0.05{\mathrm{Mpc}}^{-1}$ [51].

The untilted primordial power spectrum for untilted spatially non-flat quintessential inflation models is defined as [130,190] here $q=\sqrt{k^2+K}$ is the wavenumber for scalar perturbations. In the spatially flat limit $K=0$, $P\left(q\right)$ reduces to the $n_s=1$ spectrum.

As mentioned above, one of the major unresolved problems of modern cosmology is the so called coincidence problem in the $\mathsf{\Lambda }$CDM model, in which the energy densities of dark energy and dark matter are of the same order of magnitude at the present epoch. One way to resolve this problem is to assume that these components somehow interact with each other. IDE models consider the transformation of dark energy and dark matter into each other, with their interaction between described by the following modified continuity equations for dark energy and matter, respectively where ${\rho }_{\mathrm{m}}$ is the matter energy density, ${\rho }_{Œ}$ and $p_{Œ}$ are respectively the energy density and pressure of the scalar field, and ${\delta }_{\mathrm{couple}}$ is the coupling coefficient between dark energy and matter. In IDE models, the following forms of the coupling coefficient ${\delta }_{\mathrm{couple}}$ are typically used [136,137] where $n=\sqrt{8\pi /M_{\mathrm{pl}}^2}$, and $\beta$ and Q are dimensionless constants. The IDE models are subdivided into two types, as described below [88,136,137].

The observed magnitudes of type Ia supernovas are one of the best data for constraining the distance redshift relationship through the determination of the luminosity distance. In the $\varphi$CDM models, the distances tend to be smaller compared to the $\mathsf{\Lambda }$CDM predictions at the same redshift. This provides an opportunity for differentiating these models from each other.

One of the first studies in this direction was performed in Podariu and Ratra [193]. They used three datasets of SNe Ia apparent magnitude versus redshift: (i) R98 data [162], both including and excluding the unclassified SN Ia 1997ck at $z=0.97$ (with 50 and 49 SNe Ia apparent magnitude data, respectively), (ii) P99 data [2], and (iii) a third set with the corrected/effective stretch factor magnitudes for the 54 Fit C SNe Ia of P99 apparent magnitude data, and obtained constraints on the $\varphi$CDM model with RP potential ($\varphi$CDM-RP model) (see Figure 4).

Caresia et al. [194] obtained constraints on parameters of the $\varphi$CDM with the RP and Sugra potentials [177,195], and also the extended quintessence models with the inverse power-law RP potential [196,197], from datasets of apparent magnitude versus redshift measurements of 176 SNe Ia [2,162,198], and data of the SuperNovae Acceleration Probe (SNAP) satellite [199]1. Obtained constraints on model parameters are shown in Figure 5 and Figure 6. No useful constraints on the model parameters were found for the Sugra potential, while $1\sigma$ constraints of $\alpha <0.8$ and $\alpha <0.6$, for both the extended and ordinary quintessence models using the RP potential, were obtained using the SNe Ia apparent magnitude and SNAP satellite data respectively.

Doran et al. [200] considered a DE model parameterized as [200] where $R\left(a\right)=ln({\mathsf{\Omega }}_{\varphi }\left(a\right)/(1-{\mathsf{\Omega }}_{\varphi }\left(a\right))$ corresponds to the constant b is defined by the EoS parameter at present epoch $w_0$, the dark energy density parameter at present epoch ${\mathsf{\Omega }}_{\varphi 0}$, and the parameter ${\mathsf{\Omega }}_e$ characterizing the amount of dark energy at early times to which it asymptotes for very large redshifts, as Using a combination of datasets from SNe Ia [162], Wilkinson Microwave Anisotropy Probe (WMAP) [201], Cosmic Background Imager (CBI) [202], Very Small Array (VSA) [203], SDSS [204], and HST [205], the authors find $w_0<-0.8$ and ${\mathsf{\Omega }}_e<0.03$ at the $1\sigma$ confidence level; the contours are shown in Figure 7. It should be noted that the SNe Ia apparent magnitude data are most sensitive to $w_0$, while CMB temperature anisotropies and LSS growth rate are the best constraints of ${\mathsf{\Omega }}_e$ (see Figure 8).

Fuzfa & Alim [206] studied the $\varphi$CDM model with the RP and Sugra potentials in a spatially closed universe. The estimated values of ${\mathsf{\Omega }}_{\mathrm{m}0}$ and ${\mathsf{\Omega }}_{\varphi 0}$, using SNe Ia apparent magnitude data from the SNLS collaboration [207], are quite different than those for the standard spatially flat $\mathsf{\Lambda }$CDM model (Figure 9). Such a result is expected due to the different cosmic acceleration and dark matter clustering predicted between quintessence models and the standard $\mathsf{\Lambda }$CDM model, arising from differences in cosmological parameters, even at $z=0$. The quintessence scalar field creates more structures outside the filaments, lighter halos with higher internal velocity dispersion, as seen from N-body simulations performed by the authors to study the influence of quintessence on the distribution of matter on large scales.

Pavlov et al. [208] found that for the $\varphi$CDM-RP model in a spacetime with nonzero spatial curvature, the dynamical scalar field has an attractor solution in the curvature dominated epoch, while the energy density of the scalar field increases relative to that of the spatial curvature. In the left panel of Figure 10, we see that the values ${\mathsf{\Omega }}_{\mathrm{m}0}=0.27$ and $\alpha =3$ are consistent with these constraints for a range of values ${\mathsf{\Omega }}_{\mathrm{k}0}$, for a set of values of $H_0{t}_0$2. The right panel of Figure 10 shows a similar analysis for several values of the cosmological test parameter $▵({\mathsf{\Omega }}_{\mathrm{m}0},{\mathsf{\Omega }}_{\mathrm{k}0},\alpha )=\delta \left(t_0\right)/(1+z_i)\delta \left(t_i\right)$, where $\delta \left(t_0\right)$ and $\delta \left(t_i\right)$ are the values of the matter density contrast at respectively at the present time $t_0$, and an arbitrary time $t_i$ such that $a\left(t_i\right)\ll a\left(t_0\right)$, i.e., a time when the universe is well approximated by the Einstein-de Sitter model in the matter-dominated epoch.

Farooq et al. [209] constrained the $\varphi$CDM-RP model in a spacetime with non-zero spatial curvature, as well as the XCDM model, using the Union2.1 compilation of 580 SNe Ia apparent magnitude measurements of Suzuki et al. [210], Hubble parameter observations [28,30,211,212], and the $0.1\le z\le 0.75$ BAO peak length scale measurements [21,22,213] (see Figure 11). They constrain the spatial curvature density parameter today to be $|{\mathsf{\Omega }}_{\mathrm{k}0}|\le 0.15$ at $1\sigma$ confidence level, and more precise data are required to tighten the bounds on the parameters.

Assuming that the Hubble constant $H_0$ tension of the $\mathsf{\Lambda }$CDM model is actually a tension on the SNe Ia absolute magnitude $M_B$, Nunes & Di Valentino [139] assessed the $M_B$ tension by comparing of spatially flat $\mathsf{\Lambda }$CDM model, wCDM and IDE models using two datasets compilation SNe Ia Pantheon sample [163]+ BAO [22,24,25,215,216] + big bang nucleosynthesis (BBN) [217] and SNe Ia Pantheon sample + BAO [22,24,25,215,216] + BBN [217] + $M_B$ [218] (see Figure 12). They found that the IDE model can alleviate both the $M_B$ and $H_0$ tensions with a coupling different from zero at 2$\sigma$ confidence level with a preference for a compilation of Pantheon + BAO + BBN + $M_B$ datasets.

The CMB provides a very accurate determination of the angular diameter distance at a redshift of $z\sim 1000$. This measurement is sensitive to the entire expansion of the universe over this wide range of redshifts. As pointed out before, $\varphi$CDM models tend to predict smaller distances and can therefore be constrained with the CMB geometric measurements.

In one of the first such studies, Doran et al. [219] used CMB temperature anisotropy data from the BOOMERANG and MAXIMA experiments [220,221] to distinguish quintessential inflation models with different EoS parameters, described by a kinetic term $k\left(\varphi \right)$ of the cosmon field3: (i) the RP potential with $k\left(\varphi \right)=1$, (ii) the leaping kinetic term model which has $V\left(\varphi \right)=M_{\overline{\mathrm{pl}}}^{4}exp(-\varphi /M_{\overline{\mathrm{pl}}})$, $M_{\overline{\mathrm{pl}}}=\sqrt{8\pi }{M}_{\mathrm{pl}}$ is the reduced Plank mass, $k\left(\varphi \right)=k_{\min }+tanh[(\varphi -{\varphi }_1)/M_{\overline{\mathrm{pl}}}]+1$, ${\varphi }_1\approx 277$ eV and $k_{\min }=[0.05,0.1,0.2,0.26]$ [222], and (iii) the exponential potential with $V\left(\varphi \right)=M_{\overline{pl}}^{4}exp(-\sqrt{2}\alpha \varphi /M_{\overline{\mathrm{pl}}})$, $\alpha =\sqrt{3/2{\mathsf{\Omega }}_{\varphi 0}}$, and $k\left(\varphi \right)=1$ [134]. The dark energy density parameters today and at last scattering epoch, ${\mathsf{\Omega }}_{\varphi 0}$ and ${\mathsf{\Omega }}_{Œ\mathrm{ls}}$, and the averaged EoS parameter of the field $\varphi$ are used to parameterize the separation of peaks in CMB temperature anisotropies, which can be used to measure the value of ${\mathsf{\Omega }}_{\varphi }$ before the last scattering (see Figure 13).

Caldwell et al. [223] investigated how early quintessence dark energy, i.e., a non-negligible quintessence energy density during the recombination and structure formation epochs, affects the baryon-photon fluid and the clustering of dark matter, and thus the CMB temperature anisotropy and the matter power spectra. They showed that early quintessence is characterized by a suppressed ability to cluster at small length scales, as suggested by the compilation of data from WMAP [224,225], CBI [226,227], Arcminute Cosmology Bolometer Array Receiver (ACBAR) [228], the LSS growth rate dataset of Two degree Field (2dF) Galaxy Redshift Survey [229,230,231], and $L_{y-\alpha }$ forest [232,233]; these are shown in Figure 14. Furthermore, quintessential inflation models are compatible with these data for a constant spectral index of primordial density perturbations, as seen in the left panel of Figure 14.

Pettorino et al. [234] studied a class of extended $\varphi$CDM models, where the scalar field is exponentially coupled to the Ricci scalar and is described by the RP potential. The projection of the ISW effect on the CMB temperature anisotropy4 is found to be considerably larger in the exponential case with respect to a quadratic non-minimal coupling as seen in Figure 15. This reflects the fact that the effective gravitational constant depends exponentially on the dynamics of the scalar field.

Mukherjee et al. [235] conducted a likelihood analysis of the Cosmic Background Explorer - Differential Microwave Radiometers (COBE-DMR) sky maps to normalize the $\varphi$CDM-RP model in flat space5. As seen from Figure 16, this model remains an observationally viable alternative to the standard spatially flat $\mathsf{\Lambda }$CDM model.

Samushia and Ratra [238] constrainted model parameters of the $\mathsf{\Lambda }$CDM model, the XCDM model, and the $\varphi$CDM-RP model using galaxy cluster gas mass fraction data [239], for this, they introduced an auxiliary random variable as opposed to integrating over nuisance parameters of the Markov Chain Monte Carlo (MCMC) method. Two different sets of priors were chosen to study the influence of the type of priors on the obtained results – one set has [7] $h=0.73\pm 0.03$, ${\mathsf{\Omega }}_{\mathrm{b}}{h}^2=0.0223\pm 0.0008$ (1$\sigma$ errors), and the other set has $h=0.68\pm 0.04$ [240,241], and ${\mathsf{\Omega }}_{\mathrm{b}}{h}^2=0.0205\pm 0.0018$ [242]. The obtained constraints on the $\varphi$CDM model with the RP potential are shown in Figure 17. We see that ${\mathsf{\Omega }}_{\mathrm{m}}$ is better constrained than $\alpha$, whose best-fit value is $\alpha =0$, corresponding to the standard spatially flat $\mathsf{\Lambda }$CDM model, however, the scalar field $\varphi$CDM model is not excluded.

Chen et al. [243] constrained the scalar field $\varphi$CDM-RP model and the $\mathsf{\Lambda }$CDM model with massive neutrinos assuming two different neutrino mass hierarchies in both the spatially flat and non flat universes, using a joint dataset comprising of CMB temperature anisotropy data [12,244], BAO peak length scale data from 6dF Galaxy Survey (6dFGS), from SDSS - Main Galaxy Sample (MGS), Baryon Oscillation Spectroscopic Survey (BOSS)-LOWZ (galaxies within the redshift range $0.2<z<0.43$), BOSS CMASS-DR11 (galaxies within the redshift range $0.43<z<0.7$) [23], the joint light-curve analysis (JLA) compilation from SNe Ia apparent magnitude measurements [245], and the Hubble Space Telescope $H_0$ prior observations [29]. Assuming three species of degenerate massive neutrinos, they found the $2\sigma$ upper bounds of $\sum m_{\nu }<0.165\left(0.299\right)$ eV and $\sum m_{\nu }<0.164\left(0.301\right)$ eV respectively for the spatially flat (spatially non-flat) $\mathsf{\Lambda }$CDM model and the spatially flat (spatially non-flat) $\varphi$CDM model (see Figure 18). The inclusion of spatial curvature as a free parameter leads to a significant expansion of the confidence regions for $\sum m_{\nu }$ and other parameters in spatially flat $\varphi$CDM models, but the corresponding differences are larger for both the spatially non-flat $\mathsf{\Lambda }$CDM and spatially non-flat $\varphi$CDM models.

Park & Ratra [246] constrained the spatially flat tilted and spatially non-flat untilted $\varphi$CDM-RP inflation model by analyzing CMB temperature anisotropy angular power spectrum data from the Planck 2015 mission [247], BAO peak length scale measurements [26], a Pantheon collection of 1048 SNe Ia apparent magnitude measurements over the broader redshift range of $0.01<z<2.3$ [163], Hubble parameter observation [21,25,28,30,31,32,33,34,211,248], and LSS growth rate measurements [25] (see Figure 19 and Figure 20). Constraints on parameters of the spatially non-flat model was improved from $1.8\sigma$ to more than $3.1\sigma$ confidence level by combining CMB temperature anisotropy data with other datasets. Present observations favor a spatially closed universe with the spatial curvature contributing about two-thirds of a percent of the current total cosmological energy budget. The spatially flat tilted $\varphi$CDM model is a $0.4\sigma$ better fit to the observational data than is the standard spatially flat tilted $\mathsf{\Lambda }$CDM model, i.e., current observational data allow for the possibility of dynamical dark energy in the universe. The spatially non flat tilted $\varphi$CDM model better fits the DES bounds on the rms amplitude of mass fluctuations ${\sigma }_8$ as a function of the matter density parameter at present epoch ${\mathsf{\Omega }}_{\mathrm{m}0}$ but it does not provide such a good agreement with the larger multipoles of Planck 2015 CMB temperature anisotropy data as the spatially flat tilted $\mathsf{\Lambda }$CDM model.

Constraints on model parameters in the XCDM and $\varphi$CDM-RP (spatially flat tilted) inflation models using the compilation of CMB [247] and BAO data [22,23,24,249,250,251] were derived by Ooba et al. [249]. The authors calculated the angular power spectra of the CMB temperature anisotropy using the CLASS code of Blas et al. (2011) [250], and executed the MCMC analysis with Monte Python (Audren et al. [251]). Results of this analysis are shown in Figure 21.

Having one additional parameter compared to the standard spatially-flat $\mathsf{\Lambda }$CDM model, both $\varphi$CDM and XCMB models better fit the TT + lowP + lensing + BAO peak length scale data than does the standard spatially-flat $\mathsf{\Lambda }$CDM model. For the $\varphi$CDM model, $▵{\chi }^2=-1.60$, and for the XCDM model, $▵{\chi }^2=-1.26$ relative to the $\mathsf{\Lambda }$CDM model. The improvement over the standard spatially-flat $\mathsf{\Lambda }$CDM model in $1.3\sigma$ and in $1.1\sigma$ for the XCDM model are not significant, but these dynamical dark energy models cannot be ruled out. Both the $\varphi$CDM and XCMB dynamical dark energy models reduce the tension between the Planck 2015 (Aghanim et al. [247]) CMB temperature anisotropy and the weak lensing constraints of the rms amplitude of mass fluctuations ${\sigma }_8$.

Park & Ratra [252] also constrained the Hubble constant $H_0$ value in the spatially flat and spatially non-flat $\mathsf{\Lambda }$CDM, XCDM, $\varphi$CDM-RP models using various combinations of datasets: BAO peak length scale measurements [26], a Pantheon collection of 1048 SNe Ia apparent magnitude measurements over the broader redshift range of $0.01<z<2.3$ [163], and Hubble parameter observations [21,25,28,30,31,32,33,34,211,248]. The resulting constraints are shown in Figure 22. According to this analysis, the dataset slightly favors to the untilted spatially non-flat dynamical XCDM and $\varphi$CDM quintessential inflation models, as well as smaller Hubble constant $H_0$ values.

The compilation of the South Pole Telescope polarization (SPTpol) CMB temperature anisotropy data [253], alone and in combination with Planck 2015 CMB temperature anisotropy data [247] and non-CMB temperature anisotropy data, consisting of the Pantheon Type SNe Ia apparent magnitude measurements [163], BAO peak length scale measurements [22,24,25,26,248], Hubble parameter $H\left(z\right)$ data [21,28,30,31,32,33,34,211], and LSS growth rate data [25] was used by Park & Ratra [254] to obtain constraints on parameters of the spatially flat and untilted spatially non-flat $\mathsf{\Lambda }$CDM, XCDM, scalar field $\varphi$CDM-RP quintessential quintessential inflation models The results obtained from constraints on parameters of the untilted spatially non-flat $\varphi$CDM model with the inverse power-law RP potential using only SPTpol CMB temperature anisotropy data and with combination of other datasets are presented in Figure 23. In each dark energy model, constraints on cosmological parameters from SPTpol measurements, Planck CMB temperature anisotropy and non-CMB temperature anisotropy measurements are largely consistent with one another. Smaller angular scale SPTpol measurements (used jointly with only Planck CMB temperature anisotropy data or with the combination of Planck CMB temperature anisotropy data and non-CMB temperature anisotropy data) favor the untilted spatially closed models.

Di Valentino et al. [88] explored the IDE models to find out if these models can resolve both the Hubble constant $H_0$ tension problem of the standard spatially flat $\mathsf{\Lambda }$CDM model and resolve the contradictions between observations of Hubble constant in high and low redshifts in the spatially non-flat $\mathsf{\Lambda }$CDM scenario.

The authors constrained on parameters of the spatially flat IDE and $\mathsf{\Lambda }$CDM models as well as the spatially non-flat IDE and $\mathsf{\Lambda }$CDM models applying CMB Planck 2018 data [13], BAO [22,24,25] measurements, 1048 data points in the redshift range $z\in (0.01,2.3)$ of the Pantheon SNe Ia luminosity distance data [163], a Gaussian prior of the Hubble constant ($H_0=74.03\pm 1.42\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$ at 1$\sigma$ CL), obtained from a reanalysis of HST data by the SH0ES collaboration [81].

Based on the results of this observational analysis, it was found that the Planck 2018 CMB data favor spatially closed hypersurfaces at more than 99% CL (with a significance of 5$\sigma$). While a larger value of the Hubble constant, i.e., alleviation of the Hubble constant tension (with a significance of 3.6$\sigma$) has been obtained for the spatially non-flat IDE models. The authors concluded that searches for other forms of the interaction function $\xi$ and the EoS for the dark energy component in IDE models are needed, which may further ease the tension of the Hubble constant. 1$\sigma$ and 2$\sigma$ confidence level contours on parameters of the spatially non-flat IDE model are shown in Figure 24.

Another potentially powerful probe of $\varphi$CDM signatures is the growth rate in low redshift LSS. The growth rate is expected to be stronger in $\varphi$CDM models compared to their $\mathsf{\Lambda }$CDM counterparts.

Pavlov et al. [255] constrained the spatially flat $\varphi$CDM-RP, the XCDM, the wCDM, and the $\mathsf{\Lambda }$CDM models from future LSS growth rate data, by considering that the full sky space-based survey will observe $H_{\alpha }$-emitter galaxies over 15000 ${\mathrm{deg}}^2$ of the sky. For the bias and density of observed galaxies, they applied the predictions of Orsi et al. [256] and Geach et al. [257], respectively. They also assumed that half of the galaxies would be detected within the reliable redshift range, which roughly reflects the expected outcomes of proposed space missions, such as the ESA Euclidean Space Telescope (EUCLID) mission and the NASA Wide-Field Infrared Telescope (WFIRST) mission. The obtained results are shown in Figure 25, where we see that measurements of the LSS growth rate in the near future will be able to constrain scalar field $\varphi$CDM models with an accuracy of about 10%, considering the fiducial spatially flat $\mathsf{\Lambda }$CDM model, an improvement of almost an order of the magnitude compared to those from currently available datasets [213,258,259,260,261,262,263]. Constraints on the growth index parameter $\gamma$ are the most restrictive in the $\mathsf{\Lambda }$CDM model than in other models. In the $\varphi$CDM model, constraints on the growth index parameter $\gamma$ are about a third more tighter than in the wCDM and XCDM models.

Pavlov et al. [264] also obtained constraints on the above DE models from Hubble parameter $H\left(z\right)$ observations [28,30,211,212], from Union2.1 compilation of 580 SNe Ia apparent magnitude measurements [210], and a compilation of 14 independent LSS growth rate measurements within the redshift range $0.067\le z\le 0.8$ [21,22,213,265]. The authors performed two joint analyses, first for the combination of $H\left(z\right)$ and SNe Ia apparent magnitude data, and the other for measurements of LSS growth rate, Hubble parameter $H\left(z\right)$ and SNe Ia apparent magnitude; the results of these analyses are presented in Figure 26. Constraints on cosmological parameters of the spatially flat $\varphi$CDM model from LSS growth rate data are quite restrictive. In combination with SNe Ia apparent magnitude versus redshift data and Hubble parameter measurements, LSS growth rate data are consistent with the standard spatially flat $\mathsf{\Lambda }$CDM model, as well as with the spatially flat $\varphi$CDM model.

Avsajanishvili et al. [266] constrained the parameters ${\mathsf{\Omega }}_{\mathrm{m}}$ and $\alpha$ in the spatially flat $\varphi$CDM-RP model. Applying only measurements of the LSS growth rate [267], the authors obtained a strong degeneracy between the model parameters ${\mathsf{\Omega }}_{\mathrm{m}}$ and $\alpha$, Figure 27 (Left panel). This was followed by obtaining constraints from a compilation of data from the LSS growth rate measurements [267], and the distance-redshift ratio of the BAO peak length scale observations and prior distance from CMB temperature anisotropy [268], which eliminated the degeneracy between ${\mathsf{\Omega }}_{\mathrm{m}}$ and $\alpha$, giving ${\mathsf{\Omega }}_{\mathrm{m}}=0.30\pm 0.04$ and $0\le \alpha \le 1.30$ at $1\sigma$ confidence level (the best-fit value for the model parameter $\alpha$ is $\alpha =0$). Constraints on ${\mathsf{\Omega }}_{\mathrm{m}}$ and $\alpha$ from data compilation of Gupta et al. (2012) and Giostri et al. (2012) are presented in Figure 27 (Right panel).

Avsajanishvili et al. [269] also constrained various quintessence and phantom scalar field $\varphi$CDM models presented in Table 1 and Table 2 using observational data predicted for the Dark Energy Spectroscopic Instrument (DESI) [248]. The parameters of these models were constrained using MCMC methods by comparing measurements of the expansion rate of the universe $H\left(z\right)$, the angular diameter distance, and the LSS growth rate predicted for the standard spatially flat $\mathsf{\Lambda }$CDM model with corresponding values calculated for the $\varphi$CDM models. Results of constraints for the Zlatev-Wang-Steinhardt potential, the phantom pNGb potential, and the inverse power-law RP potential are shown in Figure 28 and Figure 29. To compare quintessence and phantom models, Bayesian statistical tests were conducted, namely the Bayes factor, as well as the $AIC$ and $BIC$ information criteria, were calculated. The $\varphi$CDM scalar field models could not be unambiguously preferred, from the DESI predictive data, over the standard $\mathsf{\Lambda }$CDM spatially flat model, the latter still being the most preferred dark energy model. The authors also investigated how the $\varphi$CDM models can be approximated by the CPL parametrization, by plotting the CPL-$\mathsf{\Lambda }$CDM 3 sigma confidence level contours, using MCMC techniques, and displayed on them the largest ranges of the current EoS parameters for each $\varphi$CDM model. These ranges were obtained for different values of model parameters or initial conditions from the prior ranges. The authors classified the scalar field models based on whether they can or cannot be distinguished from the standard spatially flat $\mathsf{\Lambda }$CDM model at the present epoch, as seen in Figure 31. They found that all studied models can be divided into two classes: models that have attractor solutions and models whose evolution depends on initial conditions.

Peracaula et al. [270] constrained the spatially flat $\mathsf{\Lambda }$CDM, XCDM, and $\varphi$CDM-RP models by constructing three datasets: DS1/SP consisting of SNe Ia apparent magnitude+$H\left(z\right)$+BAO peak length scale+LSS growth rate+CMB temperature anisotropy data with matter power spectrum SP; DS1/BSP consisting of SNe Ia apparent magnitude+$H\left(z\right)$+BAO peak length scale+LSS growth rate+CMB temperature anisotropy data with both matter power spectrum and bispectrum; and DS2/BSP, which involves BAO peak length scale+LSS growth rate+CMB temperature anisotropy data with both matter power spectrum and bispectrum. These datasets include 1063 SNe Ia apparent magnitude data [79,163], 31 measurements of $H\left(z\right)$ from cosmic chronometers [35,211], 16 BAO peak length scale data [271,272], LSS growth rate data, specifically 18 points from data [21,272,273], one point from the weak lensing observable $S_8$ [274], full CMB likelihood from Planck 2015 TT+lowP+lensing [12]. The obtained constraints are shown in Figure 32 and Figure 33. The authors tested the effect of separating the expansion history data (SNe Ia apparent magnitude+$H\left(z\right)$) from CMB temperature anisotropy characteristics and LSS formation data (BAO peak length scale+LSS), where LSS includes redshift-space distortions (RSD) and weak lensing measurements, and found that the expansion history data are not particularly sensitive to the dynamic effects of dark energy, while compilation of data BAO peak length scale+LSS+CMB temperature anisotropy is more sensitive. Also the influence of the bispectral component of the matter correlation function on the dynamics of dark energy is studied. For this, BAO peak length scale+LSS data were considered, including both the conventional power spectrum and the bispectrum. As a result, when the bispectral component is excluded, the results obtained are consistent with previous studies by other authors, which means that no clear signs of dynamical dark energy have been found in this case. On the contrary, when the bispectrum component was included in the BAO peak length scale+LSS growth rate dataset for the $\varphi$CDM model, a significant dynamical dark energy signal was achieved at $2.5-3\sigma$ confidence level. The bispectrum can therefore be a very useful tool for tracking and examining the possible dynamical features of dark energy and their influence on the LSS formation in the linear regime.

Park & Ratra [275] constrained the tilted spatially flat and untilted spatially non-flat XCDM model by applying the Planck 2015 CMB temperature anisotropy data [247], BAO peak length scale measurements [26], a Pantheon collection of 1048 SNe Ia apparent magnitude measurements over the broader redshift range $0.01<z<2.3$ [163], Hubble parameter observations [21,25,28,30,31,32,33,34,211,248], and LSS growth rate measurements [25], and obtained results as shown in Figure 34 and Figure 35. These data slightly favor the spatially closed XCDM model over the spatially flat $\mathsf{\Lambda }$CDM model at $1.2\sigma$ confidence level, while also being in better agreement with the untilted spatially flat XCDM model than with the spatially flat $\mathsf{\Lambda }$CDM model at the $0.3\sigma$ confidence level. Current observational data is unable to rule out dynamical dark energy models. The dynamical untilted spatially nonflat XCDM model is compatible with the Dark Energy Survey (DES) limits on the current value of the rms mass fluctuations amplitude ${\sigma }_8$ as a function of the matter density parameter at present epoch ${\mathsf{\Omega }}_{\mathrm{m}0}$, but it does not give such a good agreement with higher multipoles of CMB temperature anisotropy data, as the standard spatially flat $\mathsf{\Lambda }$CDM model.

Samushia & Ratra [276] constrained the standard spatially flat $\mathsf{\Lambda }$CDM, the XCDM, and the$\varphi$ CDM-RP models from BAO peak length scale measurements [17,20], in conjunction with WMAP measurements of the apparent acoustic horizon angle, and galaxy cluster gas mass fraction measurements [239]. These constraints are presented in Figure 36. It is seen that the measurements of Percival et al. (2007) constrain the $\varphi$CDM model less effectively (left panel of Figure 36), while measurements of joint BAO peak length scale and galaxy cluster gas mass of the $\varphi$CDM model give consistent and more accurate constraints than those derived from other data, i.e., $\alpha <3.5$ (right panel of Figure 36).

The above models were also constrained by Samushia et al. [277] using lookback time versus redshift data [278], passively evolving galaxies data [211], current BAO peak length scale data, and SNe Ia apparent magnitude measurements. Applying a bayesian prior on the total age of the universe based on WMAP data, the authors obtained constraints on the $\varphi$CDM model as shown in Figure 37. Constraints on the $\varphi$CDM model by joint datasets consisting of measurements of the age of the universe, SNe Ia Union apparent magnitude and BAO peak length scale are more tighter than those obtained from datasets consisting of data of the lookback time and the age of the universe.

The quintessential inflation model with the generalized exponential potential $V\left(\varphi \right)\propto exp(-\lambda {\varphi }^n/M_{\mathrm{pl}}^n),n>1$ was studied by Geng et al. [131]. The authors extended this model including massive neutrinos that are non-minimally coupled with a scalar field, obtaining observational constraints on parameters using combinations of data: CMB temperature anisotropy [244,244], BAO peak length scale from BOSS [23,263], and 11 SNe Ia apparent magnitudes from Supernova Legacy Survey (SNLS) [207]. It was found that the upper bound on possible values of the sum of neutrino masses $\sum m_{\nu }<2.5$ eV is significantly larger than in the spatially flat $\mathsf{\Lambda }$CDM model (Figure 38). The authors concluded that the model under consideration is in good agreement with observations and represents a successful scheme for the unification of primordial inflaton field causing inflation in the very early universe and dark energy causing the accelerated expansion of the universe at the present epoch.

The compilation of CMB angular power spectrum data from the Planck 2015 mission [247], and BAO peak length scale measurements from the matter power spectra obtained by missions: 6dFGS [22], BOSS, LOWZ and CMASS [23], and SDSS-MGS [24] was applied by Ooba et al. [279] to obtain constraints on the spatially non-flat quintessential inflation $\varphi$CDM-RP model. The theoretical angular power spectra of the CMB temperature anisotropy were calculated using the Cosmic Linear Anisotropy Solving System (CLASS) code of Blas et al. [250], and the MCMC analysis was performed with Monte Python of Audren et al. [251]. The results of this analysis are presented in Figure 39. The authors also used a physically consistent power spectrum for energy density inhomogeneities in the spatially non-flat (spatially closed) quintessential inflation $\varphi$CDM model, and found that the spatially closed $\varphi$CDM model provides a better fit to the lower multipole region of CMB temperature anisotropy data compared to that provided by the tilted spatially flat $\mathsf{\Lambda }$CDM model. The former reduces the tension between the Planck and the weak lensing ${\sigma }_8$ constraints, while the higher multipole region of the CMB temperature anisotropy data is in better agreement with the tilted spatially flat $\mathsf{\Lambda }$CDM model than with the spatially closed $\varphi$CDM model (see Figure 40).

Ryan et al. [280] constrained the parameters of the $\varphi$CDM-RP, the XCDM, and the $\mathsf{\Lambda }$CDM models from BAO peak length scale measurements [22,24,25,26,248], and the Hubble parameter $H\left(z\right)$ data [21,28,30,31,32,33,34,211]. The results obtained for the $\varphi$CDM model are presented in Figure 41, which shows that this dataset is consistent with the standard spatially flat $\mathsf{\Lambda }$CDM model. Depending on the value of the Hubble constant $H_0$ as a prior and the cosmological model under consideration, data provides evidence in favor of the spatially non-flat scalar field $\varphi$CDM model.

Chudaykin et al. [281] obtained constraints on the parameters of the oCDM, XCDM (here $w_0$CDM), and wCDM models by using the joint analysis from data of BAO peak length scale, BBN and SNe Ia apparent magnitude. The resulting constraints are completely independent of the CMB temperature anisotropy data but compete with the CMB temperature anisotropy constraints in terms of parameter error bars. The authors consequently obtained the value of the spatial curvature density parameter at present epoch ${\mathsf{\Omega }}_{\mathrm{k}0}=-0.{043}_{-0.036}^{+0.036}$ at 1$\sigma$ confidence level, which is consistent with the spatially flat universe; in the spatially flat XCDM model, the value of the dark energy EoS parameter at present epoch $w_0=-1.{031}_{-0.048}^{+0.052}$ at 1$\sigma$ confidence level, which approximately equals to the value of the EoS parameter for the $\mathsf{\Lambda }$CDM model; values of the $w_0$ and $w_a$ in the CPL parameterization of the EoS parameter of the wCDM model $w_0=-0.{98}_{-0.11}^{+0.099}$ and $w_a=-0.{33}_{-0.48}^{+0.63}$ at 1$\sigma$ confidence level. The authors also found that the exclusion of the SNe Ia apparent magnitude data from the joint data analysis does not significantly weaken the resulting constraints. It means that when using a single external BBN prior, full-shape and BAO peak length scale data can provide reliable constraints independent of CMB temperature anisotropy constraints. The resulting constraints on model parameters of the spatially flat XCDM and wCDM models are shown in Figure 42. The authors also tightened the observational constraints on cosmological parameters with the inclusion of the hexadecapole ($l=4$) moment of the redshift-space power spectrum (see Figure 43).

Bernui et al. [67] investigated the effect of BAO measurements on IDE models that have significantly different dynamic behavior compared to the prediction of the standard $\mathsf{\Lambda }$CDM model. The authors used the compilation of 15 transversal 2D BAO measurements [282,283] and CMB data [86] to constrain IDE models. It was found that transversal 2D BAO and traditional 3D BAO measurements can generate completely different observational constraints on the coupling parameter in IDE models. Moreover, in contrast to the joint Planck + BAO analysis, where it is not possible to solve the Hubble constant $H_0$ tension, the joint Planck + BAO (transversal) analysis agrees well with the measurements made by the SH0ES team and applied to the IDE models, solves the Hubble constant $H_0$ tension. 1$\sigma$ and 2$\sigma$ confidence level contours constraints on the coupling parameter $\xi$ in IDE model using the 2D transversal 2D BAO are shown in Figure 44.

Samushia & Ratra [284] used the Simon, Verde & Jimenez (SVJ) [211] definition of the redshift dependence of the Hubble parameter $H\left(z\right)$ (so-called SVJ $H\left(z\right)$ data) to constrain cosmological parameters in the scalar field $\varphi$CDM-RP model. According to the results obtained (see Figure 45), using the $H\left(z\right)$ data, the constraints on the matter density parameter ${\mathsf{\Omega }}_{\mathrm{m}}$ are more stringent than those on the model parameter $\alpha$. Constraints on the matter density ${\mathsf{\Omega }}_{\mathrm{m}}$ are approximately as tight as the ones derived from the galaxy cluster gas mass fraction data [285] and from the SNe Ia apparent magnitude data [286].

Chen & Ratra [287] analyzed constrains on the model parameters of the $\varphi$CDM-RP, the XCDM, and the $\mathsf{\Lambda }$CDM models, using 13 Hubble parameter $H\left(z\right)$ data versus redshift [28,212]. The authors showed (see Figure 46) that the Hubble parameter $H\left(z\right)$ data yield quite strong constraints on the parameters of the $\varphi$CDM model. The constraints derived from the $H\left(z\right)$ measurements are almost as restrictive as those implied by the currently available lookback time observations, and the GRB luminosity data, but more stringent than those based on the currently available galaxy cluster angular size data. However, they are less restrictive than those following from the joint analysis of SNe Ia apparent magnitude and BAO peak length scale data. The joint analysis of the Hubble parameter $H\left(z\right)$ data with SNe Ia apparent magnitude and BAO peak length scale data favor the standard spatially flat $\mathsf{\Lambda }$CDM model but do not exclude the dynamical scalar field $\varphi$CDM model.

In [288], Farooq et al. obtained constrains on the parameters of the $\varphi$CDM-RP, the XCDM, the wCDM, and the $\mathsf{\Lambda }$CDM models from analysis of measurements of the BAO peak length scale, SNe Ia apparent magnitude [210], 21 Hubble parameter $H\left(z\right)$ [28,30,211,212]. The results of this analysis are shown in Figure 47. Constraints are more restrictive with the inclusion of 8 new $H\left(z\right)$ measurements [30] than those derived by Chen & Ratra [287]. This analysis favors the standard spatially flat $\mathsf{\Lambda }$CDM model but does not exclude the scalar field $\varphi$CDM model.

Farooq & Ratra [289] worked out constrains on the parameters of the $\varphi$CDM-RP, the XCDM, and the $\mathsf{\Lambda }$CDM models from measurements of the Hubble parameter $H\left(z\right)$ at redshift $z=2.3$ [290] and 21 lower redshift measurements [28,30,211,212]. Constraints with the inclusion of the new $H\left(z\right)$ measurement of Busca et al. are more restrictive than those derived by Farooq et al., Figure 48. As seen in this figure, the $H\left(z\right)$ constraints depend on the Hubble constant prior to $H_0$ used in the analysis. The resulting constraints are more stringent than those which follow from measurements of the SNe Ia apparent magnitude of Suzuki et al. (2012). This joint analysis consisting of measurements of $H\left(z\right)$, SNe Ia apparent magnitude, and BAO peak length scale favor the standard spatially flat $\mathsf{\Lambda }$CDM model, but the dynamical scalar field $\varphi$CDM model is not excluded as well.

Farooq & Ratra [291] found constraints on the parameters of the $\varphi$CDM-RP model from the compilation of 28 independent measurements of the Hubble parameter $H\left(z\right)$ within the range of redshift $0.07\le z\le 2.3$. Measurements of $H\left(z\right)$ require a currently accelerating cosmological expansion at $3\sigma$ confidence level. The authors determined the deceleration-acceleration transition redshift $z_{\mathrm{da}}=0.74\pm 0.05$. This result is in good agreement with the result obtained by Busca et al. [290], which is $z_{\mathrm{da}}=0.82\pm 0.08$ based on 11 measurements of $H\left(z\right)$ from BAO peak length scale data within the range of redshift $0.2\le z\le 2.3$. The resulting constraints with different priors of $H_0$ are demonstrated in Figure 49.

Farooq et al. [214] analyzed constraints on the parameters of the spatially flat $\varphi$CDM-RP, the XCDM, and the $\mathsf{\Lambda }$CDM models from a compilation of measurements of the Hubble parameter $H\left(z\right)$. To get this compilation, the authors used weighted mean and median statistics techniques to combine 23 independent lower redshifts $z<1.04$, Hubble parameter $H\left(z\right)$ measurements, and define binned forms of them. Then this compilation was combined with 5 $H\left(z\right)$ measurements at the higher redshifts $1.3\le z\le 2.3$. The resulting constraints are shown in Figure 50. As seen from the figure, the weighted mean binned data are almost identical to those derived from analysis using 28 independent measurements of $H\left(z\right)$. Binned weighted-mean values of $H\left(z\right)/(1+z)$ versus redshift data are presented in Figure 51. These results are consistent with a moment of the deceleration-acceleration transition at redshift $z_{da}=0.74\pm 0.05$ derived by Farooq & Ratra [291], which corresponds to the standard spatially flat $\mathsf{\Lambda }$CDM model.

Chen et al. [292] used 28 measurements of Hubble parameter $H\left(z\right)$ within the redshift range $0.07\le z\le 2.3$ [21,28,30,31,211,290,293] to define the value of the Hubble constant $H_0$ in the $\varphi$CDM-RP, the wCDM, and the spatially flat and spatially non-flat $\mathsf{\Lambda }$CDM models. The result obtained for the $\varphi$CDM-RP model is shown in Figure 52. The value of the Hubble constant $H_0$ is found as: for the spatially flat and spatially non-flat $\mathsf{\Lambda }$CDM model, $H_0=68.3_{-3.3}^{+2.7}\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$ and $H_0=68.4_{-3.3}^{+2.9}\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$; for the wCDM model, $H_0=65.0_{-6.6}^{+6.5}\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$; for the $\varphi$CDM model, $H_0=67.9_{-2.4}^{+2.4}\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$ (at 1$\sigma$ confidence level). The obtained $H_0$ values are more consistent with the smaller values determined from the recent CMB temperature anisotropy and BAO peak length scale data and with the values derived from the median statistics analysis of Huchra’s compilation of $H_0$ data.

Farooq et al. [142] determined constraints on the parameters of the $\varphi$CDM-RP, the XCDM, the wCDM, and the $\mathsf{\Lambda }$CDM models in the spatially flat and spatially non-flat universe. The authors used the updated compilation of 38 measurements of the Hubble parameter $H\left(z\right)$ within the redshift range $0.07\le z\le 2.36$ [21,25,28,30,31,32,33,34,211,248]. The result for these constraints is shown in Figure 53. The authors determined the redshift of the cosmological deceleration-acceleration transition, $z_{\mathrm{da}}$, and the value of the Hubble constant $H_0$ from the $H\left(z\right)$ measurements. The determined values of $z_{\mathrm{da}}$ are insensitive to the chosen model, and depend only on the assumed value of the Hubble constant $H_0$. The weighted mean of these measurements is $z_{\mathrm{da}}=0.72\pm 0.05(0.84\pm 0.03)$ for $H_0=68\pm 2.8(73.24\pm 1.74)\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$. The authors proposed a model-independent method to determine the value of the Hubble constant $H_0$. The $H\left(z\right)$ data are consistent with the standard spatially flat $\mathsf{\Lambda }$CDM model while do not rule out the spatially non-flat XCDM and spatially non-flat $\varphi$CDM models.

Ryan et al. [143] determined constraints on the parameters of the spatially flat and spatially non-flat $\mathsf{\Lambda }$CDM, XCDM, $\varphi$CDM-RP models using BAO peak length scale measurements [22,24,25,26,248], the Hubble parameter $H\left(z\right)$ data [21,30,31,32,33,34,211], and quasar (QSO) angular size data [294,295]. 1$\sigma$, 2$\sigma$ and 3$\sigma$ confidence level contours constraints on the parameters of the spatially non-flat $\varphi$CDM model with the RP potential from $H\left(z\right)$, QSO and BAO peak length scale datasets are presented in Figure 54. Depending on the chosen model and dataset, the observational data slightly favor both the spatially closed hypersurfaces with ${\mathsf{\Omega }}_{\mathrm{k}0}<0$ at $1.7\sigma$ confidence level, and the dynamical dark energy models over the standard spatially flat $\mathsf{\Lambda }$CDM model at a slightly higher than $2\sigma$ confidence level. Furthermore, depending on the dataset and the model, the observational data favor a lower Hubble constant value than the one measured by the local data at $1.8\sigma$ confidence level to $3.4\sigma$ confidence level.

Cao et al. [296] found constraints on the parameters of the spatially flat and non-flat $\mathsf{\Lambda }$CDM, XCDM, and $\varphi$CDM-RP models using $H_{II}$ starburst galaxy apparent magnitude measurements [297,298], the compilation of 1598 X-ray and UV flux measurements of QSO 2015 data within the redshift range $0.036\le z\le 5.1003$ and 2019 QSO data [299,300] only and in conjunction with BAO peak length scale measurements [22,24,25,26,248], Hubble parameter $H\left(z\right)$ data [21,28,30,31,32,33,34,211]. The constraints on the parameters of the spatially flat and spatially non-flat $\varphi$CDM model with the RP potential obtained from datasets mentioned above are shown in Figure 55. A combined analysis of all datasets leads to the relatively model-independent and restrictive estimates for the values of matter density parameter at present epoch ${\mathsf{\Omega }}_{\mathrm{m}0}$ and the Hubble constant $H_0$. Depending on the cosmological model, these estimates are consistent with a lower value of $H_0$ in the range of $2.0\sigma$ to $3.4\sigma$ confidence level. Combined datasets favor the spatially flat $\mathsf{\Lambda }$CDM, while at the same time do not rule out dynamical dark energy models.

The compilation of 1598 X-ray and UV flux measurements of QSO 2015 data within the redshift range $0.036\le z\le 5.1003$ and 2019 QSO data [299,300] only and in conjunction with BAO peak length scale measurements [22,24,25,26,248], Hubble parameter $H\left(z\right)$ data[21,28,30,31,32,33,34,211] was applied by Khadka & Ratra [145] to impose constraints on the parameters of the tilted spatially flat and untilted spatially non-flat $\mathsf{\Lambda }$CDM, XCDM, and $\varphi$CDM-RP quintessential inflation models. Obtained constraints for the untilted spatially non-flat $\varphi$CDM-RP model from the combination of various datasets and extended QSO data only are presented in Figure 56. In most of the models, the QSO data favor the values of the matter density parameter ${\mathsf{\Omega }}_{\mathrm{m}0}\sim 0.5-0.6$, while in a combined analysis of QSO data with $H\left(z\right)$ + BAO peak length scale dataset, the values of the matter density parameter at present epoch ${\mathsf{\Omega }}_{\mathrm{m}0}$ are shifted slightly towards larger values. A combined set of data QSO + BAO peak length scale + $H\left(z\right)$ is consistent with the standard spatially flat $\mathsf{\Lambda }$CDM model, but favors slightly both the spatially closed hypersurfaces and the dynamical dark energy models.

Khadka & Ratra [146] obtained constraints on the parameters of the tilted spatially flat and untilted spatially non-flat $\mathsf{\Lambda }$CDM, XCDM, $\varphi$CDM-RP quintessential inflation models from a compilation of 808 X-ray and UV flux measurements of QSOs (quasi-stellar objects) within the redshift range $0.061\le z\le 6.28$ alone [299] and in conjunction with BAO peak length scale measurements [22,24,25,26,248], Hubble parameter $H\left(z\right)$ data [21,28,30,31,32,33,34,211]. 1$\sigma$, 2$\sigma$ and 3$\sigma$ confidence level contours constraints on the parameters of the untilted spatially non-flat $\varphi$CDM model with the RP potential from the combination of various datasets are presented in Figure 57. The constraints using only the QSO data are significantly weaker but consistent with those from the combination of the $H\left(z\right)$ + BAO peak length scale data. Combined analysis from QSO + $H\left(z\right)$ + BAO peak length scale data is consistent with the standard spatially flat $\mathsf{\Lambda }$CDM model but slightly favors both closed spatial hypersurfaces and the untilted spatially non-flat $\varphi$CDM model.

Cao et al. [301] found constraints on the parameters of the spatially flat and non-flat $\mathsf{\Lambda }$CDM, XCDM, and $\varphi$CDM-RP models using the higher-redshift GRB data [302,303], starburst galaxy ($H_{II}$G) measurements [297,298,304], and QSO angular size (QSO-AS) data [294,295]. Constraints from the combined analysis of cosmological parameters of the spatially flat and non-flat $\varphi$CDM-RP model are presented in Figure 58. The constraints from the combined analysis of these datasets are consistent with the currently accelerating cosmological expansion, as well as with the constraints obtained from the analysis of the Hubble parameter $H\left(z\right)$ data and the measurements of the BAO peak length scale. From the analysis of the $H\left(z\right)$ + BAO peak length scale + QSO-AS + $H_{II}$G + GRB dataset, the model-independent values of the matter density parameter at present epoch ${\mathsf{\Omega }}_{\mathrm{m}0}=0.313\pm 0.013$ and the Hubble constant $H_0=69.3\pm 1.2\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$ are obtained. This analysis favors the spatially flat $\mathsf{\Lambda }$CDM model but also does not rule out dynamical dark energy models.

Khadka & Ratra [147] determined constraints on the parameters of the spatially flat and non-flat $\mathsf{\Lambda }$CDM, XCDM, and $\varphi$CDM-RP models from the compilation of X-ray and UV flux measurements of 2038 QSOs which span the redshift range $0.009\le z\le 7.5413$ [300,305]. The obtained results are shown in Figure 59. The authors found that for the full QSO dataset, parameters of the X-ray and UV luminosities $L_X-L_{UV}$ relation used to standardize these QSO data depends on the cosmological model, and therefore cannot be used to constrain the cosmological parameters in these models. Subsets of these QSOs, limited by redshift $z\le 1.5-1.7$, obey the $L_X-L_{UV}$ relation in a way that is independent of the cosmological model and can therefore be used to constrain the cosmological parameters. Constraints from these smaller subsets of lower redshift QSO data are generally consistent, but much weaker than those inferred from the Hubble parameter $H\left(z\right)$ and the BAO peak length scale measurements.