The mass-varying neutrino scenario is analyzed for three trial quintessence potentials (Ferreira-Joyce, inverse exponential, and thawing oscillating). The neutrino mass is generated via Yukawa coupling to the scalar field which represents dark energy. The inverse exponential and oscillating potentials are shown to successfully generate the neutrino masses in the range $m \sim 10^{-2}-10^{-3}~$eV and to yield the current dark energy density in the regime of the late-time acceleration of the Universe. Depending on the choice of potentials, the acceleration could occur in two different regimes: (1) the regime of instability, and (2) the stable regime. The first regime of instability is after the Universe underwent a first-order transition and is rolling toward the new stable vacuum. The imaginary sound velocity $c^2_s < 0$ in this regime implies growing fluctuations of the neutrino density (clustering). In the second regime, the Universe smoothly changes its stable states via a continuous transition. Since $c^2_s > 0$, the neutrino density is stable. For all cases the predicted late-time acceleration of the Universe is asymptotically very close to that of the $\Lambda$CDM model. Further extensions of the theory to modify the neutrino sector of the Standard Model and to incorporate inflation are also discussed.

We investigated the evolution of the background expansion and the growth rate of the matter density fluctuations in the scalar field $\phi$CDM Ratra-Peebles model. We constrained the model parameter $\alpha$ and the matter density parameter $\Omega_{\rm m}$ using the recent measurements of the growth rate of the matter density fluctuations and baryon acoustic oscillation peak position. In addition, we studied a number of the $\phi$CDM scalar field models in order to determine whether or not it would be possible to discriminate them from the $\Lambda$CDM model using predicted data for the future Dark Energy Spectroscopic Instrument (DESI) observations. For this purpose, we carried out the statistical Bayesian analysis, and computed Bayes coefficients, as well as Akaike and Bayesian information criteria. We found that projected DESI results could provide compelling when comparing most $\phi$CDM models with the $\Lambda$CDM model. We also conducted the MCMC analysis and obtained the constraints on the parameters of the scalar field models, comparing the observational data for: the universe expansion rate, the angular diameter distance and the growth rate function, with the corresponding data generated for the $\Lambda$CDM model. We investigated how well the Chevallier-Polarsky-Linder (CPL) parametrization approximates the various scalar field models. We determined the location of scalar field model in the phase space of the CPL parameters.

Most dark energy models have the $\Lambda$CDM as their limit, and if future observations constrain our universe to be close to $\Lambda$CDM Bayesian arguments about the evidence and the fine-tuning will have to be employed to discriminate between the models. Assuming a baseline $\Lambda$CDM model we investigate a number of quintessence and phantom dark energy models, and we study how they would perform when compared to observational data, such as the expansion rate, the angular distance, and the growth rate measurements, from the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. We sample posterior likelihood surfaces of these dark energy models with Monte Carlo Markov Chains while using central values consistent with the Planck $\Lambda$CDM universe and covariance matrices estimated with Fisher information matrix techniques. We find that for this setup the Bayes factor provides a substantial evidence in favor of the $\Lambda$CDM model over most of the alternatives. We also investigated how well the CPL parametrization approximates various scalar field dark energy models, and identified the location for each dark energy model in the CPL parameter space.

We explore the scalar field quintessence freezing model of dark energy with the inverse Ratra-Peebles potential. We study the cosmic expansion and the large scale structure growth rate. We use recent measurements of the growth rate and the baryon acoustic oscillation peak positions to constrain the matter density $\Omega_\mathrm{m}$ parameter and the model parameter $\alpha$ that describes the steepness of the scalar field potential. We solve jointly the equations for the background expansion and for the growth rate of matter perturbations. The obtained theoretical results are compared with the observational data. We perform the Baysian data analysis to derive constraints on the model parameters.

Dark Energy models with slowly-rolling cosmological scalar field provide a popular alternative to the standard, time-independent cosmological constant model. We study simultaneous evolution of background expansion and growth in the scalar field model with the Ratra-Peebles self-interaction potential. We use recent measurements of the linear growth rate and the baryon acoustic oscillation peak positions to constrain the model parameter $\alpha$ that describes the steepness of the scalar field potential.