  The detection of the GW170817/GRB170817A event improved the constraints on
the propagation speed of gravitational waves, thus placing possible variations
caused by dark energy under restraint. For models based on scalar fields
belonging to the family of Horndeski Lagrangians, non-minimal derivative
couplings are now severely constrained, entailing a substantially limited
phenomenology. In this work we want to stress that there is still a plethora of
dark energy models that get around this obstacle while still providing
interesting phenomenologies able to distinguish them from the standard
cosmology. We focus on a class involving vector fields as a proxy, but our
discussion is extensible to a broader class of models. In particular, we show
the possibility of having a non-minimal derivative coupling giving a
non-trivial effect on scalar modes without affecting gravitational waves and
the possibility of having a second tensor mode that can oscillate into
gravitational waves. We also present a novel class of configurations breaking
rotational invariance but with an energy-momentum tensor that is isotropic
on-shell. This peculiar feature makes the scalar and vector sectors of the
perturbations mix so that, even in a perfectly isotropic background cosmology,
preferred direction effects can appear in the perturbations. We also comment on
models that give rise to isotropic solutions when averaging over rapid
oscillations of the vector fields. The explored models are classified according
to distinctive field configurations that provide inequivalent realisations of
the Cosmological Principle.
