- Calculate the potential energy of this arrangement.

- Do the same calculation using and as units of energy and length, respectively. These parameters then vanish from the expression for the pair energy, and the calculation is done with quantities of order .

- With the above units for energy and length, together with the atomic mass unit, compute the metric value of the self-consistent unit of time? Let one of the particles have a metric speed , typical of the thermal velocities of atoms or small molecules. What is the value of in self-consistent units?

*Advice:* It is convenient to count the lower, left and front
face of the cube as belonging to the basic cell, while the three other
faces belong to the next periodic cells.

To test for correct arragement of the particles, compute the diagnostic

which is sometimes called ``melting factor''. For a

By scaling all lengths, adjust the
volume such that the reduced number density becomes .

__PROJECT MD:__
Augment the subroutine STARTCONF by a procedure that assigns random
velocities to the particles, making sure that the
sum total of each velocity component is zero.

__PROJECT MC/MD:__
The second subroutine will serve to compute the total potential
energy in the system, assuming a Lennard-Jones interaction and applying
the nearest image convention:

Write such a subroutine and call it ENERGY. Use it to compute the energy in the system created by STARTCONF.

F. J. Vesely / University of Vienna