  Non-abelian symmetries are thought to be incompatible with many-body
localization, but have been argued to produce in certain disordered systems a
broad non-ergodic regime distinct from many-body localization. In this context,
we present a numerical study of properties of highly-excited eigenstates of
disordered chains with SU(3) symmetry. We find that while weakly disordered
systems rapidly thermalize, strongly-disordered systems indeed exhibit
non-thermal signatures over a large range of system sizes, similar to the one
found in previously studied SU(2) systems. Our analysis is based on the
spectral, entanglement, and thermalization properties of eigenstates obtained
through large-scale exact diagonalization exploiting the full SU(3) symmetry.
