Update figure comparing Debye and Einstein models. Added experimental data for silver.

src/2_debye_model.md

@@ -102,8 +102,8 @@ where $x=\frac{\hbar\omega}{k_{\rm B}T}$ and $\Theta_{\rm D}\equiv\frac{\hbar\om
 
 ```python
 
-T = [1.35,2.,3.,4.,5.,6.,7.,8.,10.,12.,14.,16.,20.,28.56,36.16,47.09,55.88,65.19,74.56,83.91,103.14,124.2,144.38,166.78,190.17,205.3]
-c = [0.,0.,0.,0.,0.,0.,0.0719648,0.1075288,0.2100368,0.364008,0.573208,0.866088,1.648496,4.242576,7.07096,10.8784,13.47248,15.60632,17.27992,18.6188,20.33424,21.63128,22.46808,23.05384,23.47224,23.68144]
+T = np.array([1.35,2.,3.,4.,5.,6.,7.,8.,10.,12.,14.,16.,20.,28.56,36.16,47.09,55.88,65.19,74.56,83.91,103.14,124.2,144.38,166.78,190.17,205.3])
+c = np.array([0.,0.,0.,0.,0.,0.,0.0719648,0.1075288,0.2100368,0.364008,0.573208,0.866088,1.648496,4.242576,7.07096,10.8784,13.47248,15.60632,17.27992,18.6188,20.33424,21.63128,22.46808,23.05384,23.47224,23.68144])
 
 def c_einstein(T, T_E):
     x = T_E / T
@@ -120,7 +120,7 @@ def c_debye(T, T_D):
 fig, ax = pyplot.subplots()
 
 ax.scatter(T, c)
-ax.set_title('Heat capacity of silver compared to the Debye and Einstein models')
+#ax.set_title('Heat capacity of silver compared to the Debye and Einstein models')
 ax.plot(T, c_einstein(T, 151), label='Einstein model')
 ax.plot(T, c_debye(T, 215), label='Debye model')
 ax.set_ylim(bottom=0, top=26)
@@ -128,7+128,7 @@
 ax.set_xlabel('T(K)')
 ax.set_ylabel(r'C(J/mol-K)')
 ax.set_xticks([0, 100, 200])
 ax.set_yticks([24.945])
 ax.set_yticklabels(['$3R$'])
 ax.legend(loc='lower right')
 pyplot.hlines([24.945], 0, T[-1], linestyles='dashed')