There are only 16 particles, so the fluctuations in the kinetic energy are quite large, which make the equilbration hard to distinguish. To help distinguish, we did 16000 time steps. The data appears to be equilibrated after 200 data points. That's 2000 timesteps (the data is taken every 10 steps) and 2000 x.006 = 12 in simulation time units.
So throwing away the first 200 data points, we get an average kinetic energy of 1.493 +/- .025 and an average potential energy of -0.705 +/- .025.
The kinetic energy is higher than the desired temperature because the initial state (a lattice) had high potential energy, which converted to kinetic energy and raised the temperature during equilbration.
h = 0.01 | h = 0.1 | h = 1.0 |
Return to Exercise 1
Return to Main Page