Details

At low temperatures, the Dulong–Petit law fails and the specific heat of solid is no longer constant. The Einstein quantum model for a solid correctly explains this phenomenon. We can make a simplified model of a solid with a 2D cell array (Black and Ogborn, 1976) with one oscillator in each cell with an energy of a multiple of ℏω, and thermal interaction reproduced through a transfer of a quantum of energy between a pair of these cells chosen randomly. The number of points assigned to each cell is the number of quanta possessed by the oscillator (occupancy). The initial state assigns one quantum to every cell. After about 1,000 steps, the distribution of a thermal equilibrium state and Boltzmann distribution is obtained; it is stable after more steps. Blue bars represent instantaneous values , purple bars represent the mean values given by the canonical distribution (equal to ), and the red curve represents the fit of the occupancies.