First, the system must be initialized: let β β = 1 / k b T be the system's Boltzmann temperature and initialize the system with an initial state (which can be anything since the final result should not depend on it).