11.2: Heat Capacity of Solids
- Page ID
- 372726
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)In the early 1800’s French scientists Pierre Dulong and Alexis Petit noticed that many solid elements had near-identical per molar heat capacities as shown in Figure 11.6. This led to the Dulong–Petit law: \(C_m=3R\). Previously in Ch. 2 we encountered gases with such high heat capacities due to their three translational and three rotational motions, which imparts six degrees of freedom that rendered \(C_m=3R\) due to the equipartition theorem (Table 2.1). While this is perfectly understandable for a gas, it’s hard to reconcile why heavy metallic elements such as lead and a gas such as methane have the same heat capacities. The answer is simple- solids have three (vibrational) potential energy and three kinetic energy degrees of freedom due to their bonding. Thus, the equipartition theorem still appears valid for solids, unless, those solids are cold, in which case something goes terribly wrong as shown in Figure 11.7.
From the data in Figure 11.7 we can see that, at low temperatures, solids of copper and tungsten have almost no heat capacity! And while it rises with temperature to the Dulong-Petit limit, it was quite the curiosity for many years why this behavior was observed in many materials. In fact, it was one of Albert Einstein’s early contributions in 1907 to approximately describe the increase in heat capacity as vibrations in the solid state, called phonons, that behave just like Planck’s photons. In this regard, instead of energy being related to the wavelength of light it was due to the frequency of an atomic vibration in the solid state. All one has to do is simply change the definition of average internal energy (U), as derived from the Boltzmann formula, from the Planck to the Einstein form: \(\frac{h\cdot \frac{c}{\lambda }}{e^{\frac{\hslash \cdot c}Callstack:
at (Under_Construction/Free_Energy_(Snee)/11:_Boltzmann_Statistics/11.02:_Heat_Capacity_of_Solids), /content/body/p[2]/span[1]/span, line 1, column 1
\[C=3N_A\cdot \frac{\partial }{\partial T}\left(\frac{\hslash \omega }{e^{\frac{\hslash \omega }{k_BT}}-1}\right)=3N_A\cdot \frac{-\hslash \omega }
Callstack:
at (Under_Construction/Free_Energy_(Snee)/11:_Boltzmann_Statistics/11.02:_Heat_Capacity_of_Solids), /content/body/p[3]/span[1], line 1, column 1
Callstack:
at (Under_Construction/Free_Energy_(Snee)/11:_Boltzmann_Statistics/11.02:_Heat_Capacity_of_Solids), /content/body/p[3]/span[2], line 1, column 1
This result correctly predicts that, as T\(\mathrm{\to}\)0 K, the heat capacity also \(\mathrm{\to}\)0 J/K/mol. At higher temperatures where: \(e^{\frac{\hslash \omega }{k_BT}}\sim 1+\frac{\hslash \omega }{k_BT}\) we find that:
\[C\sim 3N_Ak_B{\left(\frac{\hslash \omega }{k_BT}\right)}^2\frac{1+\frac{\hslash \omega }{k_BT}}{{\left(\frac{\hslash \omega }{k_BT}\right)}^2}\sim 3N_Ak_B{\left(\frac{\hslash \omega }{k_BT}\right)}^2{\left(\frac{k_BT}{\hslash \omega }\right)}^2=3N_Ak_B=3R \nonumber \]
All of these results make us believe that Einstein has the correct approach; however, this turns out not to be the case! The problem is that the heat capacity rises exponentially in the Einstein model whereas in reality most materials have a C\({}_{p,m}\)\(\mathrm{\sim}\)T\({}^{3}\) behavior. And as shown by Peter Debye in 1912 this was due to Einstein’s use of a single frequency for the solid-state vibrations, whereas there are a range of frequencies as shown in the inset of Figure 11.7.
The solution to address concerning the shortcoming of the Einstein model is to use a range of phonon frequencies to define the internal energy as so:
\[U=\int{}degeneracy\left(\omega \right)\cdot \langle energy\left(\omega \right)\rangle =\int{}degeneracy\left(\omega \right)\cdot \frac{\hslash \omega }{e^{\frac{\hslash \omega }{k_BT}}-1}\partial \omega \nonumber \]
So as not to labor the point one can show that the degeneracy of phonons as a function of frequency is \(\frac{V{\omega }^2}{2{\pi }^2v^3_s}\), where \(v_s\) is the speed of sound in the solid. As a result the Debye equation for the heat capacity of a solid is:\(C=\frac{\partial U}{\partial T}=\frac{\partial }{\partial T}\frac{V{\omega }^2}{2{\pi }^2v^3_s}\int^{{\omega }_D}_0{}\frac{V{\omega }^2}{2{\pi }^2v^3_s}\frac{\hslash \omega }{e^{\frac{\hslash \omega }{k_BT}}-1}\partial \omega \propto T^3\)
While we don’t need to worry further about the mathematics of the result, suffice to say the result is proportional to \(T^3\) as experimentally observed. It is also interesting to note that the upper limit of the Debye equation isn’t infinity as one might expect. This is because Debye realized that there was an upper limit to the range of phonon frequencies as shown in Figure 11.8, which we will call the Debye frequency \({\omega }_D\). Clearly, there cannot be a frequency higher than allowed by the atomic spacing. In practice, this upper limit is an empirical parameter that can be used to create a fit to experimental data.
Conclusion. The purpose of this chapter was to demonstrate that the Boltzmann formula has uses far beyond that encountered with the Maxwell-Boltzmann description of the velocities of gases. In terms of the development of physical chemistry as a science, it wasn’t long after the phenomena of Blackbody radiation and the Debye theory of heat capacity were described that scientists started to make the connection between matter and waves. This ultimately resulted in quantum theory, which is the subject of the next chapter.