Nonharmonic Phonons in α-Iron at High Temperatures

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.