From a Parametrized Post-Newtonian (PPN) [46-51] perspective, we address the question of whether or not the new degrees of freedom, represented by the PPN potentials, can lead to significant modifications in the dynamics of galaxies in the direction of rendering dark matter obsolete. Here, we focus on the study of rotation curves associated with spherically symmetric configurations. The values for the post-Newtonian parameters, which help us to classify the different metric theories of gravity, are tightly constrained mainly by solar system experiments [47,48]. However, in recent years, metric theories characterized by screening mechanisms [4,10,22,25] have become popular, due to the fact that they lead to the possibility of modifications in larger scales than the solar system while retaining the success of GR on it, allowing for violations of the constraints of the Post-Newtonian parameters. In such a context, we consider here two kinds of solutions for field equations: (i) Vacuum solutions and (ii) in the presence of a politropic distribution of matter. For the case (i) we find that the post-Newtonian corrections do not lead to modifications significant enough to be considered an alternative to the dark matter hypothesis. In the case (ii) we find that for a wide range of values for the PPN parameters {\gamma}, {\beta}, {\xi}, {\alpha}_3 , {\zeta}_1 and {\zeta}_2 , the need for dark matter is unavoidable, in order to find flat rotation curves. It is only for theories in which {\zeta}_3>0 that some resemblance of flat rotation curves is found. The latter suggests, at least for the models considered, that these are the only theories capable of replacing dark matter as a possible explanation for the dynamics of galaxies.

By considering the Einstein-Vlasov system for static spherically symmetric distributions of matter, we show that configurations with constant anisotropy parameter $\beta$ have, necessarily, a distribution function (DF) of the form $F=l^{-2\beta}\xi(\varepsilon)$, where $\varepsilon=E/m$ and $l=L/m$ are the relativistic energy and angular momentum per unit rest mass, respectively. We exploit this result to obtain DFs for the general relativistic extension of the Hypervirial family introduced by Nguyen and Lingam (2013), which Newtonian potential is given by $\phi(r)=-\phi_o /[1+(r/a)^{n}]^{1/n}$ ($a$ and $\phi_o$ are positive free parameters, $n=1,2,...$). Such DFs can be written in the form $F_{n}=l^{n-2}\xi_{n}(\varepsilon)$. For odd $n$, we find that $\xi_n$ is a polynomial of order $2n+1$ in $\varepsilon$, as in the case of the Hernquist model ($n=1$), for which $F_1\propto l^{-1}\left(2\varepsilon-1\right)\left(\varepsilon-1\right)^2$. For even $n$, we can write $\xi_n$ in terms of incomplete beta functions (Plummer model, $n=2$, is an example). Since we demand that $F\geq 0$ throughout the phase space, the particular form of each $\xi_n$ leads to restrictions for the values of $\phi_o$. For example, for the Hernquist model we find that $0\leq \phi_o \leq2/3$, i.e. an upper bounding value less than the one obtained for Nguyen and Lingam ($0\leq \phi_o \leq1$), based on energy conditions.