  We present a new reduced-order computational method for the molecular
dynamics simulation of entangled polymer systems. The polymer chains are
modeled as continuous Gaussian chains. Our algorithm is based on the
application of the molecular dynamics simulation method to the pseudospectral
representation of the Fourier modes of the chains.
