This can happen only if the gradient
vanishes everywhere. Thus, the scalar field must be
constant. Adding a time variable, the same argument shows that the scalar field must depend only on time, i.e. it must be homogeneous in surfaces of "constant time". In a cosmological model, these surfaces will be the hyperslices orthogonal to the world lines of the matter.
They don't. There are plenty of cosmological models (dust solutions) which are
homogeneous but not
isotropic, and there are models which are not even homogeneous. In another thread I recently gave a simple explicit example of a homogeneous but anisotropic model, the
Bianchi II dust, an exact dust solution which exhibits Kasner epochs similar to those exhibited by the
Mixmaster model, aka the
Bianchi IX dust. I've also given the line element defining the
Stephani dust, an simple exact dust solution which is an inhomogeneous model (a perturbation of an FRW model). The
Kantowski-Sachs dusts provide further popular examples of anisotropic but homogeneous cosmological models.
See
http://xxx.lanl.gov/abs/gr-qc/9812046 for an excellent review of cosmological models, including the question of symmetry. See the website in my sig (happy b'day to John B, happy b'day to he!) for many more citations to good review papers and textbooks.