57 Fe Mössbauer spectroscopy can probe several local structural, electronic and magnetic properties of Fe-containing systems. However, to establish a direct relationship between these properties and a system's geometric structure, the experimental Mössbauer parameters need to be analyzed via electronic structure calculations. Herein, structural, electronic and magnetic effects of iron substituents in the topological insulator Bi2Se3, as uniquely probed by 57 Fe Mössbauer spectroscopy, have been determined via spin-polarized electronic structure calculations. The iron ion substituents, of nominal Fe 3+ (S = 5/2) oxidation and spin state, are unequivocally shown to substitute Bi 3+ sites in epitaxial Bi2Se3 thin-films used for Mössbauer measurements. Concomitant with iron substitution, localized structural rearrangements take place whereby the longer Bi-Se bonds of the native system are replaced by significantly shorter Fe-Se counterparts in the Fe-containing system. The resulting distorted-octahedral environment about substituent iron ions gives rise to characteristic Mössbauer parameters (δFe ≈ 0.51 mm/s, ∆EQ ≈ 0.20 mm/s) which have been calculated in excellent agreement with measured values for Fe-doped Bi2Se3 thin films. Consistent with a substituent Fe 3+ ion's nominal high-spin electronic configuration (t ↑↑↑ 2g e ↑↑ g ), an Fe-centered spin density has been established which, nevertheless, extends towards neighboring Se atoms via direct Fe-Se bonding and concomitant Fe(d)-Se(p) hybridization.