Monte‐Carlo simulation is an efficient means of investigating electron transport in semiconductors. Nevertheless, the prevailing physical phenomena at low electric fields are usually neglected because a Monte‐Carlo algorithm is needed only at rather high fields. It is shown that reaching thermodynamic equilibrium between electrons and crystal lattice in a Monte‐Carlo simulation requires only including inelastic interactions with acoustic phonons at the center of the Brillouin zone. Electron energy and velocity distributions at equilibrium resulting from this umproved model are compared with the theoretical Fermi–Dirac distribution.