Nonequilibrium molecular dynamics simulations are reported for the shear viscosity of xenon using accurate two- and three-body potentials. The dependence of the hydrostatic pressure, and energy are observed to vary linearly with the square of the strain rate. This is contrast to the non-analytic three-halves power dependency on strain rate predicted by mode coupling theory. This result is attributed solely to the two-body potential. The main effect of the three-body potential is to alter the magnitude of the pressure, energy and viscosity profiles.