####################################################################### Ethasham Ahmed, December 2016 This is a simple molecular dynamics program that simulates the collison of two clusters. Argon has been chosen arbitrarily as the particle within the clusters. The program uses Lennard Jones potential. Velocity verlet algorithm is used to solve newtons equation of motion Equipartition theorem is used to calculate temperature Spatial dimentions (DIM) is set to 3 by default Time step (dt) is set to 1.7e-4 by default Within the code, calculations are done in reduced units which means mass, sigma and epsiolon are 1 However, by the end of the code, all values are converted back to SI units. Initial positions of the particles are determined from two input files The program writes to a .xyz file every 10 steps by defult. #######################################################################