The deuteron quadrupole moment is calculated using two covariant spectator theory model wave functions obtained from the 2007 high-precision fits to np scattering data. Included in the calculation are a new class of isoscalar np interaction currents automatically generated by the nuclear force model used in these fits. The prediction for model WJC-1, with larger relativistic P-state components, is 2.5% smaller than the experimental result, in common with the inability of models prior to 2014 to predict this important quantity. However, model WJC-2, with very small P-state components, gives agreement to better than 1%, similar to the results obtained recently from chiral effective field theory predictions to next-to-next-to-next-to-leading order
