Fe(CO)(n) (n=1-5) complexes have been studied using density functional theory (DFT) methods. Several functionals have been used in the geometry optimizations, harmonic frequencies computation and calculation of the iron-carbonyl bond dissociation energies. Coupled-cluster single double (triple) bond dissociation energies have also been computed for the smaller systems. The obtained results show that DFT methods yield reasonable geometries and vibrational frequencies. Regarding the bond dissociation energies, it is shown that the validity of the results depends on whether there is a change in the atomic state of the metal during the dissociation; When the atomic state is the same for both complexes, the bond dissociation energy computed using gradient corrected functionals is within the range of the experimental values, while when the atomic state changes, DFT methods overestimate the bond dissociation energy due to a poor description of the atomic multiplets.