Phonon density of states of a fcc crystal

Student Project by Christian Neubauer and Gernot Kraberger (2010)

The form below generates a table of where the first column is the angular frequency ω in rad/s and the second column is the density of states D(ω) in units of s/(rad m³).

D(ω) [1015 s/(rad m³)] 

ω [1012 rad/s]

Mass: m =

[10-24 kg]

Spring constant: C =


Lattice constant: a =


The Matlab files used to calculate the data are base.m, kvals.m, disp_solver.m and analy.m.

The density of states can be used to calculate the temperature dependence of thermodynamic quantities.