Variational density study of the one-dimensional Hubbard model: the need for three-particle constraints