Einstein solid

Heat capacity of an Einstein solid as a function of temperature. Experimental value of 3Nk is recovered at high temperatures.

The heat capacity of an object at constant volume V is defined through the internal energy U as

${\displaystyle C_{V}=\left({\partial U \over \partial T}\right)_{V}.}$

${\displaystyle T}$, the temperature of the system, can be found from the entropy

${\displaystyle {1 \over T}={\partial S \over \partial U}.}$

To find the entropy consider a solid made of ${\displaystyle N}$ atoms, each of which has 3 degrees of freedom. So there are ${\displaystyle 3N}$ quantum harmonic oscillators (hereafter SHOs for "Simple Harmonic Oscillators").

${\displaystyle N^{\prime }=3N}$

Possible energies of an SHO are given by

${\displaystyle E_{n}=\hbar \omega \left(n+{1 \over 2}\right)}$

where the n of SHO is usually interpreted as the excitation state of the oscillating mass but here n is usually interpreted as the number of phonons (bosons) occupying that vibrational mode (frequency). The net effect is that the energy levels are evenly spaced, and one can define a quantum of energy due to a phonon as

${\displaystyle \varepsilon =\hbar \omega }$

which is the smallest and only amount by which the energy of an SHO is increased. Next, we must compute the multiplicity of the system. That is, compute the number of ways to distribute ${\displaystyle q}$ quanta of energy among ${\displaystyle N^{\prime }}$ SHOs. This task becomes simpler if one thinks of distributing ${\displaystyle q}$ pebbles over ${\displaystyle N^{\prime }}$ boxes

or separating stacks of pebbles with ${\displaystyle N^{\prime }-1}$ partitions

or arranging ${\displaystyle q}$ pebbles and ${\displaystyle N^{\prime }-1}$ partitions

The last picture is the most telling. The number of arrangements of ${\displaystyle n}$ objects is ${\displaystyle n!}$. So the number of possible arrangements of ${\displaystyle q}$ pebbles and ${\displaystyle N^{\prime }-1}$ partitions is ${\displaystyle \left(q+N^{\prime }-1\right)!}$. However, if partition #3 and partition #5 trade places, no one would notice. The same argument goes for quanta. To obtain the number of possible distinguishable arrangements one has to divide the total number of arrangements by the number of indistinguishable arrangements. There are ${\displaystyle q!}$ identical quanta arrangements, and ${\displaystyle (N^{\prime }-1)!}$ identical partition arrangements. Therefore, multiplicity of the system is given by

${\displaystyle \Omega ={\left(q+N^{\prime }-1\right)! \over q!(N^{\prime }-1)!}}$

which, as mentioned before, is the number of ways to deposit ${\displaystyle q}$ quanta of energy into ${\displaystyle N^{\prime }}$ oscillators. Entropy of the system has the form

${\displaystyle S/k=\ln \Omega =\ln {\left(q+N^{\prime }-1\right)! \over q!(N^{\prime }-1)!}.}$

${\displaystyle N^{\prime }}$ is a huge number—subtracting one from it has no overall effect whatsoever:

${\displaystyle S/k\approx \ln {\left(q+N^{\prime }\right)! \over q!N^{\prime }!}}$

With the help of Stirling's approximation, entropy can be simplified:

${\displaystyle S/k\approx \left(q+N^{\prime }\right)\ln \left(q+N^{\prime }\right)-N^{\prime }\ln N^{\prime }-q\ln q.}$

Total energy of the solid is given by

${\displaystyle U={N^{\prime }\varepsilon \over 2}+q\varepsilon ,}$

since there are q energy quanta in total in the system in addition to the ground state energy of each oscillator. Some authors, such as Schroeder, omit this ground state energy in their definition of the total energy of an Einstein solid.

We are now ready to compute the temperature

${\displaystyle {1 \over T}={\partial S \over \partial U}={\partial S \over \partial q}{dq \over dU}={1 \over \varepsilon }{\partial S \over \partial q}={k \over \varepsilon }\ln \left(1+N^{\prime }/q\right)}$

Elimination of q between the two preceding formulas gives for U:

${\displaystyle U={N^{\prime }\varepsilon \over 2}+{N^{\prime }\varepsilon \over e^{\varepsilon /kT}-1}.}$

The first term is associated with zero point energy and does not contribute to specific heat. It will therefore be lost in the next step.

Differentiating with respect to temperature to find ${\displaystyle C_{V}}$ we obtain:

${\displaystyle C_{V}={\partial U \over \partial T}={N^{\prime }\varepsilon ^{2} \over kT^{2}}{e^{\varepsilon /kT} \over \left(e^{\varepsilon /kT}-1\right)^{2}}}$

or

${\displaystyle C_{V}=3Nk\left({\varepsilon \over kT}\right)^{2}{e^{\varepsilon /kT} \over \left(e^{\varepsilon /kT}-1\right)^{2}}.}$

Although the Einstein model of the solid predicts the heat capacity accurately at high temperatures, and in this limit

${\displaystyle \lim _{T\rightarrow \infty }C_{V}=3Nk}$,

which is equivalent to Dulong–Petit law.

Nevertheless, the heat capacity noticeably deviates from experimental values at low temperatures. See Debye model for how to calculate accurate low-temperature heat capacities.

Heat capacity is obtained through the use of the canonical partition function of a simple quantum harmonic oscillator.

${\displaystyle Z=\sum _{n=0}^{\infty }e^{-E_{n}/kT}}$

where

${\displaystyle E_{n}=\varepsilon \left(n+{1 \over 2}\right)}$

substituting this into the partition function formula yields

${\displaystyle {\begin{aligned}Z=\sum _{n=0}^{\infty }e^{-\varepsilon \left(n+1/2\right)/kT}=e^{-\varepsilon /2kT}\sum _{n=0}^{\infty }e^{-n\varepsilon /kT}=e^{-\varepsilon /2kT}\sum _{n=0}^{\infty }\left(e^{-\varepsilon /kT}\right)^{n}\\={e^{-\varepsilon /2kT} \over 1-e^{-\varepsilon /kT}}={1 \over e^{\varepsilon /2kT}-e^{-\varepsilon /2kT}}={1 \over 2\sinh \left({\varepsilon \over 2kT}\right)}.\end{aligned}}}$

This is the partition function of one harmonic oscillator. Because, statistically, heat capacity, energy, and entropy of the solid are equally distributed among its atoms, we can work with this partition function to obtain those quantities and then simply multiply them by ${\displaystyle N^{\prime }}$ to get the total. Next, let's compute the average energy of each oscillator

${\displaystyle \langle E\rangle =U=-{1 \over Z}\partial _{\beta }Z}$

where

${\displaystyle \beta ={1 \over kT}.}$

Therefore,

${\displaystyle U=-2\sinh \left({\varepsilon \over 2kT}\right){-\cosh \left({\varepsilon \over 2kT}\right) \over 2\sinh ^{2}\left({\varepsilon \over 2kT}\right)}{\varepsilon \over 2}={\varepsilon \over 2}\coth \left({\varepsilon \over 2kT}\right).}$

Heat capacity of one oscillator is then

${\displaystyle c_{V}={\partial U \over \partial T}=-{\varepsilon \over 2}{1 \over \sinh ^{2}\left({\varepsilon \over 2kT}\right)}\left(-{\varepsilon \over 2kT^{2}}\right)=k\left({\varepsilon \over 2kT}\right)^{2}{1 \over \sinh ^{2}\left({\varepsilon \over 2kT}\right)}.}$

Up to now, we calculated the heat capacity of a unique degree of freedom, which has been modeled as a quantum harmonic. The heat capacity of the entire solid is then given by ${\displaystyle C_{V}=3Nc_{V}}$, where the total number of degree of freedom of the solid is three (for the three directional degree of freedom) times ${\displaystyle N}$, the number of atoms in the solid. One thus obtains

${\displaystyle C_{V}=3Nk\left({\varepsilon \over 2kT}\right)^{2}{1 \over \sinh ^{2}\left({\varepsilon \over 2kT}\right)}.}$

which is algebraically identical to the formula derived in the previous section.

The quantity ${\displaystyle T_{\rm {E}}=\varepsilon /k}$ has the dimensions of temperature and is a characteristic property of a crystal. It is known as the Einstein temperature.[2] Hence, the Einstein crystal model predicts that the energy and heat capacities of a crystal are universal functions of the dimensionless ratio ${\displaystyle T/T_{\rm {E}}}$. Similarly, the Debye model predicts a universal function of the ratio ${\displaystyle T/T_{\rm {D}}}$, where ${\displaystyle T_{\rm {D}}}$ is the Debye temperature.