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³).
The Matlab file used to calculate the data are bccdos.m.
In the long wavelength limit (near ω = 0) we expect that the density of states should increase like ω². For the case of atoms on a bcc lattice connected to only their nearest neighbors by linear springs, the density of states increases linearly in ω for small ω. If linear springs to next nearest neighbor atoms are included, the density of states increases like ω² as it does for the phonon dispersion curve shown below. The matlab script also handles next nearest neighbors.
The density of states can be used to calculate the temperature dependence of thermodynamic quantities.