A description of many-particle systems, which is more fundamental than the fluid approach, is to consider them as a kinetic gas. In this approach the dynamical variable in which the properties of the system are encoded, is the distribution of the gas particles in position and velocity space, called 1-particle distribution function (1PDF). However, when the gravitational field of a kinetic gas is derived via the Einstein-Vlasov equations, the information about the velocity distribution of the gas particles is averaged out and therefore lost. We propose to derive the gravitational field of a kinetic gas directly from its 1PDF, taking the velocity distribution fully into account. We conjecture that this refined approach could possibly account for the observed dark energy phenomenology.

We present a systematic analysis of the dynamics of flat Friedmann-Lema\^{i}tre-Robertson-Walker cosmological models with radiation and dust matter in generalized teleparallel $f(T)$ gravity. We show that the cosmological dynamics of this model is fully described by a function $W(H)$ of the Hubble parameter, which is constructed from the function $f(T)$. After reducing the phase space to two dimensions we derive the conditions on $W(H)$ for the occurrence of de Sitter fixed points, accelerated expansion, crossing the phantom divide, and finite time singularities. Depending on the model parameters it is possible to have a bounce (from contraction to expansion) or a turnaround (from expansion to contraction), but cyclic or oscillating scenarios are prohibited. As an illustration of the formalism we consider power law $f(T) = T + \alpha(-T)^n$ models, and show that these allow only one period of acceleration and no phantom divide crossing.

We compute the parametrized post-Newtonian parameter $\gamma$ in the case of a static point source for multiscalar-tensor gravity with completely general nonderivative couplings and potential in the Jordan frame. Similarly to the single massive field case $\gamma$ depends exponentially on the distance from the source and is determined by the length of a vector of non-minimal coupling in the space of scalar fields and its orientation relative to the mass eigenvectors. Using data from the Cassini tracking experiment, we estimate bounds on a general theory with two scalar fields. Our formalism can be utilized for a wide range of models, which we illustrate by applying it to nonminimally coupled Higgs SU(2) doublet, general hybrid metric-Palatini gravity, linear ($\Box^{-1}$) and quadratic ($\Box^{-2}$) nonlocal gravity.

We calculate the PPN parameters gamma and beta for scalar-tensor gravity with a generic coupling function omega and scalar potential V in the Jordan conformal frame in the case of a static spherically symmetric source. Since the potential generally introduces a radial dependence to the effective gravitational constant as well as to gamma and beta, we discuss the issue of defining these PPN parameters and compare our expressions with previous calculations in simpler cases. We confront our results with current observational constraints on the values of gamma and beta and thus draw restrictions on the form of the functions omega and V around their asymptotic background values.