  The deuteron quadrupole moment is calculated using two CST 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 that 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 $\chi$EFT predictions to order N$^3$LO.
