Abstract
Creation and annihilation operators that replace one ion in a crystal by an impurity in one of a set of states are defined. To use statistical mechanics with these operators, an ensemble that is close to the grand canonical ensemble is introduced. From the model Hamiltonian (given at a fixed value of the volume per crystal site) and the empirical knowledge of the pressure, most of the thermodynamic quantities are derived. A simple example is worked out in detail.