Correlation between geometry, electronic structure, and magnetism of solids is both intriguing and elusive. This is particularly strongly manifested in small clusters, where a vast number of unusual structures appear. Here, we employ density functional theory in combination with a genetic search algorithm and a hybrid functional to determine the structure of gas phase clusters. For cation clusters we also calculate the corresponding vibration spectra and compare them with experiments. We successfully identify and propose structures for . Within the triangular geometric structure of , a noncollinear, ferrimagnetic, and ferromagnetic state are comparable in energy. and are ferrimagnetic with a residual magnetic moment of due to ionization. is ferrimagnetic due to the odd number of Fe atoms. We compare the electronic structure with bulk magnetite and find to be mixed valence clusters. In contrast, in and , all Fe are found to be trivalent.