🚀 go-pugleaf

Does anyone know this problem about the molecular dynamic program for a
pure hard sphere fluid. This reads input positions and velovities  from
a file, and after simulation it writes out the final positions and
velocities to a data file(Config). When changing from one state to
another, it is necessary to equilibrate the configuration at the new
state before collecting data for final statistical analysis. First,
choose a reduced density between 0.0 and 0.943. Beginning with 100
collisions, see how the final value of PV/NKT - 1 approaches an
equilibrated value by running again and again and using the utput in
CONFN as input for the next run; once equilibration is attained there
will be statistical fluctuations but no trend. How many collisions were
needed to reach equilibrium? How many collision did you use to obtain an
correct averaged value of PV/NKT - 1?

If you can solve to this problem, I deeply appreciate you.

e-mail: ylin@utkvx.utk.edu


D
D
If you can sovle for this prob