Cosmological behavior of Bekenstein's modified theory of gravity

We study the background cosmology governed by the tensor-vector-scalar theory of gravity proposed by Bekenstein. We consider a broad family of functions that lead to modified gravity and calculate the evolution of the field variables both numerically and analytically. We find a range of possible behaviors, from scaling to the late time domination of either the additional gravitational degrees of freedom or the background fluid.