Hello,
I am working on a MC simulation of a solid using NVT ensemble with Metropolis algorithm. The botlzman factor is of the form exp(-βΔU(r1,r2...rN)), where U is just the potential energy.
From thermodyanmics, dU=dq+dw (dq is the amount of heat put in, and dw is the amount of work done...