CO adsorption structures and energetics on the iron (100), (110), (111), (210), (211), and (310) surfaces from the lowest coverage up to saturation have been computed using spin-polarized density functional theory and ab initio thermodynamics. It is found that different adsorption configurations on each of these surfaces at high coverage can coexist. The stepwise adsorption energies and dissociation barriers at different coverage reveal equilibriums between desorption and dissociation of adsorbed CO molecules. Only molecular CO adsorption is possible at very high coverage and only dissociative CO adsorption at very low coverage, whereas mixed molecular and dissociative CO adsorption becomes possible at medium coverage. The computed stable adsorption configurations and the respective C–O and Fe–C stretching frequencies as well as desorption temperatures on the (100), (110), and (111) surfaces agree very well with the available experimental data. Such agreements between theory and experiment validate our co...