We consider stability properties of spherically symmetric spacetimes of stars in metric f(R) gravity. We stress that these not only depend on the particular model, but also on the specific physical configuration. Typically configurations giving the desired $\gamma_{\rm PPN} \approx 1$ are strongly constrained, while those corresponding to $\gamma_{\rm PPN} \approx 1/2$ are less affected. Furthermore, even when the former are found strictly stable in time, the domain of acceptable static spherical solutions typically shrinks to a point in the phase space. Unless a physical reason to prefer such a particular configuration can be found, this poses a naturalness problem for the currently known metric f(R) models for accelerating expansion of the Universe.

We study both analytically and numerically the gravitational fields of stars in f(R) gravity theories. We derive the generalized Tolman-Oppenheimer-Volkov equations for these theories and show that in metric f(R) models the Parameterized Post-Newtonian parameter $\gamma_{\rm PPN} = 1/2$ is a robust outcome for a large class of boundary conditions set at the center of the star. This result is also unchanged by introduction of dark matter in the Solar System. We find also a class of solutions with $\gamma_{\rm PPN} \approx 1$ in the metric $f(R)=R-\mu^4/R$ model, but these solutions turn out to be unstable and decay in time. On the other hand, the Palatini version of the theory is found to satisfy the Solar System constraints. We also consider compact stars in the Palatini formalism, and show that these models are not inconsistent with polytropic equations of state. Finally, we comment on the equivalence between f(R) gravity and scalar-tensor theories and show that many interesting Palatini f(R) gravity models can not be understood as a limiting case of a Jordan-Brans-Dicke theory with $\omega \to -3/2$.

We study the interior spacetimes of stars in the Palatini formalism of f(R) gravity and derive a generalized Tolman-Oppenheimer-Volkoff and mass equation for a static, spherically symmetric star. We show that matching the interior solution with the exterior Schwarzschild-De Sitter solution in general gives a relation between the gravitational mass and the density profile of a star, which is different from the one in General Relativity. These modifications become neglible in models for which $\delta F(R) \equiv \partial f/\partial R - 1$ is a decreasing function of R however. As a result, both Solar System constraints and stellar dynamics are perfectly consistent with $f(R) = R - \mu^4/R$.

We explore in detail a dilatonic scalar-tensor theory of gravity inspired by large extra dimensions, where a radion field from compact extra dimensions gives rise to quintessence in our 4-dimensional world. We show that the model can give rise to other types of cosmologies as well, some more akin to $k$-essence and possibly variants of phantom dark energy. In our model the field (or radius) stabilization arises from quantum corrections to the effective 4D Ricci scalar. We then show that various constraints nearly determine the model parameters, and give an example of a quintessence-type cosmology consistent with observations. We show that the upcoming SNAP-experiment would easily distinguish the present model from a constant $\Lambda$ model with an \emph{equal} amount of dark energy, but that the SNAP-data alone will not be able distinguish it from a $\Lambda$ model with about 5% less dark energy.