This paper starts from a toy model for inflation in a class of modified theories of gravity in the metric formalism. Instead of the standard procedure -- assuming a non-linear Lagrangian $f(R)$ in the Jordan frame -- we start from a simple $\phi^2$ potential in the Einstein frame and investigate the corresponding $f(R)$ in the former picture. The addition of an ad-hoc Cosmological Constant in the Einstein frame leads to a Thermodynamical interpretation of this physical system, which allows further insight on its (meta)stability and evolution.

We address the issue of the existence of inequivalent definitions of gravitational mass in $R^{2}$-gravity. We present several definitions of gravitational mass, and discuss the formal relations between them. We then consider the concrete case of a static and spherically symmetric neutron star, and solve numerically the equations of motion for several values of the free parameter of the model. We compare the features of the mass-radius relations obtained for each definition of gravitational mass, and we comment on their dependence on the free parameter. We then argue that $R^{2}$-gravity is a valuable proxy to discuss the existence of inequivalent definitions of gravitational mass in a generic modified gravity theory, and present some comments on the general case.

Understanding the influence of dark energy on the formation of structures is currently a major challenge in Cosmology, since it can distinguish otherwise degenerated viable models. In this work we consider the Top-Hat Spherical-Collapse (SC) model with dark energy, which can partially (or totally) cluster, according to a free parameter $\gamma$. The {\it lack of} energy conservation has to be taken into account accordingly, as we will show. We determine characteristic quantities for the SC model, such as the critical contrast density and radius evolution, with particular emphasis on their dependence on the clustering parameter $\gamma$.

We investigate the linear regime of f(R) = R + alpha R^2 gravity for static, spherically symmetric and asymptotically flat configurations of matter. We show that, in vacuum and deep inside the range of the extra scalar degree of freedom, the post-Newtonian parameter gamma is not equal to 1/2, as established in the literature, but it assumes larger values depending on the pressure of the star. We provide an explicit expression for gamma in terms of the mass, of the integrated pressure of the star and of the ratio between the star's radius and the range of the extra degree of freedom. We corroborate our results by providing numerical solutions for the case of a neutron star.

Critical overdensity $\delta_c$ is a key concept in estimating the number count of halos for different redshift and halo-mass bins, and therefore, it is a powerful tool to compare cosmological models to observations. There are currently two different prescriptions in the literature for its calculation, namely, the differential-radius and the constant-infinity methods. In this work we show that the latter yields precise results {\it only} if we are careful in the definition of the so-called numerical infinities. Although the subtleties we point out are crucial ingredients for an accurate determination of $\delta_c$ both in general relativity and in any other gravity theory, we focus on $f(R)$ modified-gravity models in the metric approach; in particular, we use the so-called large ($F=1/3$) and small-field ($F=0$) limits. For both of them, we calculate the relative errors (between our method and the others) in the critical density $\delta_c$, in the comoving number density of halos per logarithmic mass interval $n_{\ln M}$ and in the number of clusters at a given redshift in a given mass bin $N_{\rm bin}$, as functions of the redshift. We have also derived an analytical expression for the density contrast in the linear regime as a function of the collapse redshift $z_c$ and $\Omega_{m0}$ for any $F$.

By using the Sloan Digital Sky Survey (SDSS) first year type Ia supernova (SN Ia) compilation, we compare two different approaches (traditional \chi^2 and complete likelihood) to determine parameter constraints when the magnitude dispersion is to be estimated as well. We consider cosmological constant + Cold Dark Matter (\Lambda CDM) and spatially flat, constant w Dark Energy + Cold Dark Matter (FwCDM) cosmological models and show that, for current data, there is a small difference in the best fit values and $\sim$ 30% difference in confidence contour areas in case the MLCS2k2 light-curve fitter is adopted. For the SALT2 light-curve fitter the differences are less significant ($\lesssim$ 13% difference in areas). In both cases the likelihood approach gives more restrictive constraints. We argue for the importance of using the complete likelihood instead of the \chi^2 approach when dealing with parameters in the expression for the variance.