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
Tuesday, May 5, 2009
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment