8.8 Exercises 281
Estimateg(r)forr∈( 0 ,L 2 ]in your argon system. The easiest way is to divide the distance
interval into bins, loop over all pairs of particles and count how many distances belong in
each bin. Time-averaging the function gives a better description of the system’s general
behaviour. Plotg(r)for temperatures where the system is in solid and liquid phases. Does
it appear as expected? How would the exactg(r)look for a perfect crystal?