Variational optimization of second order density matrices: the past, the present and the future