Tuesday, May 5, 2009

Even more on Ex4

There is a slight typo in the lecture notes that I believe may cause confusion regarding your homework. When we perform a simulation in reduced units, such as the monatomic Lennard-Jones system, Boltzmann's constant kB is absorbed into the dimensionless temperature:

T* = kB T / epsilon

Thus, the acceptance criterion does not explicitly contain the quantity kB. Another way of thinking about this is that kB = 1 for reduced units:

Pacc = min{ 1, exp[ -(U2* - U1*) / T*] }

The point is, you should not need to specify kB in your simulation runs.

MSS

No comments: