Several molecular dynamics codes are routinely executed on large scale Cray vector and scalar MPP systems, with excellent scalability and performance results. Here we present some applications performance profiles and how the core algorithms involved can be optimally implemented in order to exploit the different architectural features of the Cray X1E and XT3.