In this work, we will be testing four different general \textit{f(R)}-gravity models, two of which are the more realistic models (namely the Starobinsky and the Hu-Sawicki models), to determine if they are viable alternative models to pursue a more vigorous constraining test upon them. For the testing of these models, we use 359 low- and intermediate-redshift Supernovae Type 1A data obtained from thRede SDSS-II/SNLS2 Joint Light-curve Analysis (JLA). We develop a Markov Chain Monte Carlo (MCMC) simulation to find a best-fitting function within reasonable ranges for each \textit{f(R)}-gravity model, as well as for the Lambda Cold Dark Matter ($\Lambda$CDM) model. For simplicity, we assume a flat universe with a negligible radiation density distribution. Therefore, the only difference between the accepted $\Lambda$CDM model and the \textit{f(R)}-gravity models will be the dark energy term and the arbitrary free parameters. By doing a statistical analysis and using the $\Lambda$CDM model as our "true model", we can obtain an indication whether or not a certain \textit{f(R)}-gravity model shows promise and requires a more in-depth view in future studies. In our results, we found that the Starobinsky model obtained a larger likelihood function value than the $\Lambda$CDM model, while still obtaining the cosmological parameters to be $\Omega_{m} = 0.268^{+0.027}_{-0.024}$ for the matter density distribution and $\bar{h} = 0.690^{+0.005}_{-0.005}$ for the Hubble uncertainty parameter. We also found a reduced Starobinsky model that are able to explain the data, as well as being statistically significant.

We explore exact f(R) gravity solutions that mimic Chaplygin-gas inspired LCDM cosmology. Starting with the original and modified Chaplygin gas equations of state, we reconstruct the forms of f(R) Lagrangians. The resulting solutions are generally quadratic in the Ricci scalar, but have appropriate LCDM solutions in limiting cases. These solutions, given appropriate initial conditions, can be potential candidates for scalar field-driven early universe expansion (inflation) and dark energy-driven late-time cosmic acceleration.

Over the last few years much attention has been given to the study of modified gravity theories in order to find a more natural explanation for the late time acceleration of the Universe. Nevertheless, a comparison of the matter power spectrum predictions made by these theories with available data has not yet been subjected to a detailed analysis. In the context of f(R) theories of gravity we study the predicted power spectra using both a dynamical systems approach for the background and solving for the matter perturbations without using the quasi-static approximation, comparing the theoretical results with several SDSS data. The importance of studying the first order perturbed equations by assuming the correct background evolution and the relevance of the initial conditions are also stressed. We determine the statistical significance in relation to the observational data and demonstrate their conflict with existing observations.