# 1.92: Debye's Heat Capacity Equation Fit to Experimental Data - Another Algorithm

$$\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}}$$

$\mathrm{n} :=30 \qquad \mathrm{i} :=1 \ldots \mathrm{n} \nonumber$

 $$\mathrm{T}_{\mathrm{i}} :=$$ $$\mathrm{C}_{\mathrm{i}} :=$$ 1 0.000818 3 0.0065 5 0.0243 8 0.0927 10 0.183 15 0.670 20 1.647 25 3.066 30 4.774 35 6.612 40 8.419 45 10.11 50 11.66 55 13.04 60 14.27 65 15.35 70 16.30 80 17.87 90 19.11 100 20.10 120 21.54 140 22.52 160 23.22 180 23.75 200 24.16 220 24.49 240 24.76 260 24.99 280 25.19 300 25.37

The heat capacity data were taken from the Handbook of Physics and Chemistry ‐ 72nd Edition, page 5‐71. The data are presented in units of Joules/mole/K.

Gas law constant:

$\mathrm{R} :=8.31451 \nonumber$

Define Einstein function for heat capacity and first derivative with respect $$\Theta$$:

$F(T, \Theta)=\left[\begin{array}{c}{9 \cdot R \cdot\left(\frac{T}{\Theta}\right)^{3} \cdot \int_{0}^{\frac{\Theta}{T}} \frac{x^{4} \cdot \exp (x)}{(\exp (x)-1)^{2}} d x} \\ {\frac{d}{d \Theta}\left[9 . R \cdot\left(\frac{T}{\Theta}\right)^{3} \cdot \int_{0}^{\frac{\Theta}{T}} \frac{x^{4} \cdot \exp (x)}{(\exp (x)-1)^{2}} d x\right]}\end{array}\right] \nonumber$

Call genfit to do nonlinear regression analysis.

$P :=\text { genfit }(T, C, 200, F) \qquad P=210.985 \nonumber$

Plot data and fit:

$\mathrm{t} :=1 \ldots 300 \nonumber$

This page titled 1.92: Debye's Heat Capacity Equation Fit to Experimental Data - Another Algorithm is shared under a CC BY 4.0 license and was authored, remixed, and/or curated by Frank Rioux via source content that was edited to the style and standards of the LibreTexts platform; a detailed edit history is available upon request.