Trimeron-phonon coupling in magnetite

Using density functional theory, we study the lattice dynamical properties of magnetite (${\mathrm{Fe}}_{3}{\mathrm{O}}_{4}$) in the high-temperature cubic and low-temperature monoclinic phases. The calculated phonon dispersion curves and density of states are compared with the available experimental data obtained by inelastic neutron, inelastic x-ray, and nuclear inelastic scattering. We find a very good agreement between the theoretical and experimental results for the monoclinic $Cc$ structure revealing the strong coupling between the charge-orbital (trimeron) order and specific phonon modes. For the cubic phase, clear discrepancies arise due to fluctuation effects, which are not included in the calculation method. Despite this shortcoming, we argue that the main spectral features can be understood assuming that the strong trimeron-phonon coupling is extended above the Verwey transition, with lattice dynamics influenced by the short-range order instead of the average cubic structure. Our results indicate the validity of trimerons (and trimeron-phonon coupling) to explain the physics of magnetite much beyond their original formulation.