Phonon densities of states (DOS) of bcc $\alpha$-$^{57}$Fe were measured from
room temperature through the 1044K Curie transition and the 1185K fcc
$\gamma$-Fe phase transition using nuclear resonant inelastic x-ray scattering.
At higher temperatures all phonons shift to lower energies (soften) with
thermal expansion, but the low transverse modes soften especially rapidly above
700K, showing strongly nonharmonic behavior that persists through the magnetic
transition. Interatomic force constants for the bcc phase were obtained by
iteratively fitting a Born-von K\'{a}rm\'{a}n model to the experimental phonon
spectra using a genetic algorithm optimization. The second-nearest-neighbor
fitted axial force constants weakened significantly at elevated temperatures.
An unusually large nonharmonic behavior is reported, which increases the
vibrational entropy and accounts for a contribution of 35 meV/atom in the free
energy at high temperatures. The nonharmonic contribution to the vibrational
entropy follows the thermal trend of the magnetic entropy, and may be coupled
to magnetic excitations. A small change in vibrational entropy across the
$\alpha$-$\gamma$ structural phase transformation is also reported.