Reference: Daniel V. Schroeder, *An Introduction to Thermal Physics*, (Addison-Wesley, 2000) – Problem 3.25.

The properties of an Einstein solid can be investigated analytically. The multiplicity of a solid with oscillators and energy quanta is

For large and , we can approximate this by

For very large systems, the square root factor can be neglected compared to the two power terms as it is merely “large” compared to the “very large” power terms. We thus get

From this we can get the entropy

The total energy of the system is

where is the energy of a single quantum, so we can get the temperature as defined from the entropy

We can invert this to obtain :

Finally, the heat capacity is obtained from

A plot of the dimensionless quantities versus looks like this:

To compare this with Schroeder’s Figure 1.14, we need to insert some actual numbers so we can plot versus . Figure 1.14 is for one mole of atoms, but an Einstein solid treats oscillators, where each oscillator is a one-dimensional simple harmonic oscillator. Each atom in a real solid can vibrate in three directions, so one mole of atoms in a real solid contains three moles of oscillators in an Einstein solid. Thus to compare the plots for one mole of “real” atoms, we should use .

Using and multiplying the result of by to convert the result from to (which are the units used in Figure 1.14), we then need to choose some values for to get the curves for lead, aluminum and diamond.

As a starting point, I used the values from the numerical solution done earlier. The values for lead and aluminum turned out (surprisingly, considering the numerical solution was worked out for very small systems) to work quite well, but the value for diamond needed to be adjusted a bit. The final values were

The resulting plots are:

The curves are lead (red), aluminum (green) and diamond (blue). They are a good match to the curves shown in Figure 1.14.

We can get an expression for the heat capacity in the limit of high temperatures by expanding the exponential. As a first approximation we’ll use

Then from 17

This is the expected result from the equipartition theorem, since the heat capacity for a system whose thermal energy is entirely in the form of quadratic degrees of freedom is

where is the number of degrees of freedom. For a harmonic oscillator, (one from the kinetic energy term and one from the potential energy), so we’d expect .

To get a more accurate approximation, we can expand the exponential out to a few more terms using the Taylor expansion for small :

To get the next term in the expansion beyond 22, we need to expand out to , so we get

The denominator can now be expanded again using a Taylor series. This gets quite tedious so I used Maple to do the expansion (you’re welcome to try it by hand by calculating the derivatives) with the result

Weckiai RannilaCan you explain a little further what is meant by expand the derivative and what steps would be necessary to recieve the expansion in question

gwrowePost authorI don’t know which equation(s) you’re referring to (please use the equation numbers!!) but all the expansions are just standard Taylor series. Since some of them get quite messy if you do them by hand, I’ve used Maple to do them instead.

DanteI think you missed a factor of k in the last two terms of eq. 5

gwrowePost authorFixed now. Thanks.