I am trying to equilibrate a spherical cluster of 141 Silver atoms (FCC) inside
a cubic box. I have used eam potential.
The code has run without any errors. But I observe from the 3D positions of the
atoms that the atoms move all over the spherical cluster approximately in
circles. I do not expect this. Other than the atoms at the surface, I expect
only oscillatory motion about mean positions for the atoms.
Is there a problem in my code (given below)? Any help would be much appreciated.