Use a computer to study the entropy, temperature, and heat capacity of an Einstein solid, as follows. Let the solid contain 50 oscillators (initially), and from 0 to 100 units of energy. Make a table, analogous to Table 3.2, in which each row represents a different value for the energy. Use separate columns for the energy, multiplicity, entropy, temperature, and heat capacity. To calculate the temperature, evaluate $\Delta U/\Delta S$for two nearby rows in the table. (Recall that $U=qϵ$for some constant $ϵ$.) The heat capacity $(\Delta U/\Delta T)$can be computed in a similar way. The first few rows of the table should look something like this:

(In this table I have computed derivatives using a "centered-difference" approximation. For example, the temperature $.28$is computed as $2/(7.15-0)$.) Make a graph of entropy vs. energy and a graph of heat capacity vs. temperature. Then change the number of oscillators to 5000 (to "dilute" the system and look at lower temperatures), and again make a graph of heat capacity vs. temperature. Discuss your prediction for the heat capacity, and compare it to the data for lead, aluminum, and diamond shown in Figure 1.14. Estimate the numerical value of $\epsilon$in electron-volts, for each of those real solids.

The given Einstein solid contains 50 oscillators from 0 to 100 units of energy.

That is,

$$\begin{gathered} N=50 \\ q=1\text{to}100 \end{gathered}$$

$U=qϵ$

The entropy of the system is given as:

$$\begin{gathered} \frac{S}{k}=\ln \Omega \\ \frac{S}{k}=\ln \left(\frac{q+N-1}{\lfloor q\lfloor N-1}\right) \end{gathered}$$

The temperature is given as:

$$\begin{gathered} \frac{1}{T}=\frac{\partial S}{\partial U} \\ T=\frac{\Delta U}{\Delta S} \\ T=\frac{\Delta \left(\epsilon q\right)}{\Delta \left(k\ln \Omega \right)} \\ \frac{kT}{ϵ}=\frac{\Delta \left(q\right)}{\Delta \left(\ln \Omega \right)} \end{gathered}$$

Heat capacity per oscillation is given as:

$$\begin{gathered} \frac{C}{N}=\frac{\Delta U}{N\Delta T} \\ \frac{{\displaystyle C}}{{\displaystyle N}}=\frac{\Delta \left(\epsilon q\right)}{N\Delta \left(\frac{\Delta \left(\epsilon q\right)}{\Delta \left(k\ln \Omega \right)}\right)} \\ \frac{C}{Nk}=\frac{\Delta \left(q\right)}{N\Delta \left(\frac{\Delta \left(q\right)}{\Delta \left(\ln \Omega \right)}\right)} \end{gathered}$$

The table consisting of entropy, temperature, and heat capacity of an Einstein solid for 50 oscillations with 0 to 100 units of energy can be prepared as follows:

The graphs of entropy vs energy and heat capacity vs temperature can be sketched as follows:

For 50 oscillations:

The graph of entropy vs energy can be made as:

The graph of heat capacity vs temperature can be made as:

For 5000 oscillations:

The graph of entropy vs energy can be made as:

The graph of heat capacity vs temperature can be made as:

The heat capacity is substantially smaller than the values in figure 1.14, since, for 5000 oscillations with the same units of energy, the energy is lowered, lowering the system's temperature.

For $q=1,\Omega =1$

Solving for $\epsilon$,we get,

$$\begin{gathered} U=ϵq \\ \epsilon =\frac{U}{q} \\ ϵ=\frac{Tk\left(\ln \Omega \right)}{q} \\ \epsilon =kT \end{gathered}$$

For lead, heat capacity is maximum at $T=80K$

$$\begin{gathered} {\epsilon }_{\text{lead}}=kT \\ {\epsilon }_{\text{lead}}=\left(1.38\times {10}^{-23}\right)\times 80 \\ {\epsilon }_{\text{lead}}=1.1\times {10}^{-21}\mathrm{J} \\ {\epsilon }_{\text{lead}}=6.875\times {10}^{-3}\mathrm{eV} \end{gathered}$$

For aluminum heat capacity is maximum at $T=300\mathrm{K}$

$$\begin{gathered} {\epsilon }_{alu\min um}=kT \\ {\epsilon }_{alu\min um}=\left(1.38\times {10}^{-23}\right)\times 300 \\ {\epsilon }_{alu\min um}=4.14\times {10}^{-21}\mathrm{J} \\ {\epsilon }_{alu\min um}=2.587\times {10}^{-2}\mathrm{eV} \end{gathered}$$

For diamond heat capacity is maximum at $T=2230\mathrm{K}$

$$\begin{gathered} {\epsilon }_{\text{diamond}}=kT \\ {\epsilon }_{\text{diamond}}=\left(1.38\times {10}^{-23}\right)\times 2230 \\ {\epsilon }_{\text{diamond}}=3.0774\times {10}^{-20}\mathrm{J} \\ {\epsilon }_{\text{diamond}}=0.19233\mathrm{eV} \end{gathered}$$

The table of the values can be prepared as:

The graphs can be prepared as:

For 50 oscillations:

For 5000 oscillations:

The value of $\epsilon$are:

$$\begin{gathered} {ϵ}_{\text{lead}}=6.875\times {10}^{-3}\mathrm{eV} \\ {ϵ}_{\text{aluminium}}=2.587\times {10}^{-2}\mathrm{eV} \\ {ϵ}_{\text{diamond}}=0.19233\mathrm{eV} \end{gathered}$$