Abstract
An efficient order of N molecular dynamics method for the simulation of two-body and three-body systems is presented. Due to its high speed it enables one to simulate large MD systems on a mainframe and investigate some complex processes in terms which are closer to experimental conditions (e.g. crystal growth or ion implantation on relatively large-by MD standards-substrates).