In order to see all of the important effects, consider a single block of 35 atoms (105 oscillators) with up to 300 quanta of energy. Note that in this analysis you are calculating quantities for a single block, not two blocks in contact. To make specific comparison with experimental data, consider the cases of aluminum (Al) and lead (Pb). For each metal, plot the theoretical heat capacity C per atom (not per oscillator) vs. T, with dots showing the actual experimental data. By trial and error, adjust the interatomic spring stiffness ks until your calculations approximately fit the experimental data given in the attached vpython shell. In trying to fit the data, ignore the lowest-temperature data, where electrons in the metal make a significant contribution that is not included in the Einstein model. Required:
What values of ks, obtained by trial and error, give good fits for aluminum and for lead?