Reference: Daniel V. Schroeder, An Introduction to Thermal Physics, (Addison-Wesley, 2000) – Problem 2.11.
As a more-or-less real world example of an Einstein solid, we can look at a two-state paramagnet which is a material composed of magnetic dipoles which can align parallel or antiparallel to an external magnetic field. The restriction to two states comes from quantum mechanics, since the magnetic moment can take on only certain discrete values. The two-state paramagnet is the simplest case; in other magnets, the dipole can take on more than two states relative to the applied field.
We can model the two-state paramagnet as an Einstein solid where the role of the oscillators is played by the dipoles, and an energy quantum is the difference in energy between a parallel and antiparallel dipole. The antiparallel dipole has higher energy since a torque must be applied to twist the dipole against the field.
Using this model, we can picture two interacting paramagnets as two interacting Einstein solids. As before, we have two solids, and , containing and dipoles and and quanta of energy, with . For any particular partition of the quanta, that is, for particular values of and , the total number of microstates available to the compound system is
The overall number of microstates is
Consider a system with and . Using Maple to calculate the binomial coefficients (Maple has a ‘binomial’ function that does this automatically) and produce the plot, we have
A plot of the probabilities looks like this:
As we’d expect, the most probable state is , with a probability of 0.075. The least probable state is either or , both of which have probability .