Motivated by the fact that Fe2O3 nanoparticles are used in the treatment of cancer, we have examined the role of ligands on the magnetic properties of these particles by focusing on (Fe2O3)(4) as a prototype system with H as ligands. Using the Broken-Symmetry Density Functional Theory, we observed a strong collective superexchange in the hydrogenated Fe8O12H8 cluster. The average antiferromagnetic exchange coupling constant between the four iron-iron oxo-bridged pairs was found to be-178 cm(-1), whereas coupling constants between hydroxo-bridged pairs were much smaller. We found that despite the apparent symmetry of the iron atom framework, it is not reasonable to assume this symmetry when fitting the exchange coupling constants. We also analyzed the geometrical and magnetic properties of Fe8O12Hn for n = 0-12 and found that hydrogenating oxo-bridges would generally inhibit the Fe-O-Fe antiferromagnetic superexchange interactions. Antiferromagnetic lowest total energy states become favorable only when specific distributions of hydrogen atoms are realized. The (HO)(4)-Fe-4(all spin-up)-O-4-Fe-4(all spindown)-(OH)(4) configuration in Fe8O12H8 presents such an example. This symmetric configuration can be diatomic system. considered a super