We study the scalar field potential $V(\phi)$ in the scalar-tensor gravity with self-consistent polytropic stellar configurations. Without choosing a particular potential, we numerically derive the potential inside various stellar objects. We restrict the potential to conform to general relativity or to $f(R)$ gravity inside and require the solution to arrive at SdS vacuum at the surface. The studied objects are required to obtain observationally valid masses and radii corresponding to solar type stars, white dwarfs and neutron stars. We find that the resulting scalar-tensor potential $V(\phi)$ for the numerically derived polytrope that conforms to general relativity, in each object class, is highly dependent on the matter configuration as well as on the vacuum requirement at the boundary. As a result, every stellar configuration arrives at a potential $V(\phi)$ that is not consistent with the other stellar class potentials. Therefore, a general potential that conforms to all these polytropic stellar classes could not be found.

It is stated that a class of $f(R)$ gravity models seem to obtain $\Lambda$CDM behaviour for high redshifts and general relativistic behaviour locally at high curvatures. In the present paper, we numerically study polytropic configurations that resemble stars like young sun with Hu and Sawicki $f(R)$ gravity field equations and compare the spacetime at the boundary to the general relativistic counterpart. These polytropes are stationary spherically symmetric configurations and have regular metrics at the origin. Since Birkhoff's theorem does not apply for modified gravity, the solution outside may deviate from Schwarzschild-de Sitter spacetime. At the boundary, Post-Newtonian parametrization was used to determine how much the studied model deviates from the general relativistic $\Lambda$CDM model.

We study stellar configurations and the space-time around them in metric $f(R)$ theories of gravity. In particular, we focus on the polytropic model of the Sun in the $f(R)=R-\mu^4/R$ model. We show how the stellar configuration in the $f(R)$ theory can, by appropriate initial conditions, be selected to be equal to that described by the Lane-Emden -equation and how a simple scaling relation exists between the solutions. We also derive the correct solution analytically near the center of the star in $f(R)$ theory. Previous analytical and numerical results are confirmed, indicating that the space-time around the Sun is incompatible with Solar System constraints on the properties of gravity. Numerical work shows that stellar configurations, with a regular metric at the center, lead to $\gamma_{PPN}\simeq1/2$ outside the star ie. the Schwarzschild-de Sitter -space-time is not the correct vacuum solution for such configurations. Conversely, by selecting the Schwarzschild-de Sitter -metric as the outside solution, we find that the stellar configuration is unchanged but the metric is irregular at the center. The possibility of constructing a $f(R)$ theory compatible with the Solar System experiments and possible new constraints arising from the radius-mass -relation of stellar objects is discussed.

A class of supergravity inspired quintessence models is studied by comparing to cosmological data. The set of considered models includes several previously studied quintessential potentials, as well as the $\Lambda$CDM model. We find that even though the commonly studied supergravity inspired quintessence models fit the data better than the $\Lambda$CDM model, they are a relatively poor fit when compared to the best fit model in the studied class. Our results suggest a low energy scale, less than ${\cal{M}}\sim 1$ TeV, for the effective supergravity potential.