In [4] we proved that all Gauss methods areNτ(0)-compatible for neutral delay differential equations (NDDEs) of the form:y'(t) = ay(t) + by(t - τ) + cy'(t - τ), t > 0,y(t) = g(t), -τ≤ t ≤ 0, (0.1)where a, b, c are real, τ> 0, g(t) is a continuous real valued function. In this paper we are going to use the theory of order stars to characterize the asymptotic stability properties of Gauss methods for NDDEs. And then proved that all Gauss methods are Nτ(0)-stable.