A molecular dynamics study of the organic ionic plastic crystal tetramethylammonium dicyanamide is reported here. An all atom force field, based on CHARMM parameters, has been used, and calculations were carried out under NPT conditions at several temperatures ranging from the rigid lattice (200 K) up to well above the melt (600 K). The volume expansion, radial distribution functions and the onset of different rotational and translational motions has been simulated and are discussed both in terms of plastic crystals in general as well as compared to experimental results previously obtained for this particular system. (c) 2006 Published by Elsevier B.V.