The Fe vibrational density of states (VDOS) has been determined for the heme proteins deoxymyoglobin, metmyoglobin, and cytochrome f in the oxidized and reduced states, using nuclear resonance vibrational spectroscopy (NRVS). For cytochrome f in particular, the NRVS spectrum is compared with multiwavelength resonance Raman spectra to identify those Raman modes with significant Fe displacement. Modes not seen by Raman due to optical selection rules appear in the NRVS spectrum. The mean Fe force constant extracted from the VDOS illustrates how Fe dynamics varies among these four monoheme proteins, and is correlated with oxidation and spin state trends seen in model heme compounds. The protein's contribution to Fe motion is dominant at low frequencies, where coupling to the backbone tightly constrains Fe displacements in cytochrome f, in contrast to enhanced heme flexibility in myoglobin.