Theoretical Electron Density Distributions for Fe- and Cu-Sulfide Earth Materials: A Connection between Bond Length, Bond Critical Point Properties, Local Energy Densities, and Bonded Interactions

Bond critical point and local energy density properties together with net atomic charges were calculated for theoretical electron density distributions, rho(r), generated for a variety of Fe and Cu metal-sulfide materials with high- and low-spin Fe atoms in octahedral coordination and high-spin Fe atoms in tetrahedral coordination. The electron density, rho(rc), the Laplacian, triangle down2rho(rc), the local kinetic energy, G(rc), and the oxidation state of Fe increase as the local potential energy density, V(rc), the Fe-S bond lengths, and the coordination numbers of the Fe atoms decrease. The properties of the bonded interactions for the octahedrally coordinated low-spin Fe atoms for pyrite and marcasite are distinct from those for high-spin Fe atoms for troilite, smythite, and greigite. The Fe-S bond lengths are shorter and the values of rho(rc) and triangle down2rho(rc) are larger for pyrite and marcasite, indicating that the accumulation and local concentration of rho(r) in the internuclear region are greater than those involving the longer, high-spin Fe-S bonded interactions. The net atomic charges and the bonded radii calculated for the Fe and S atoms in pyrite and marcasite are also smaller than those for sulfides with high-spin octahedrally coordinated Fe atoms. Collectively, the Fe-S interactions are indicated to be intermediate in character with the low-spin Fe-S interactions having greater shared character than the high-spin interactions. The bond lengths observed for chalcopyrite together with the calculated bond critical point properties are consistent with the formula Cu+Fe3+S2. The bond length is shorter and the rho(rc) value is larger for the FeS4 tetrahedron displayed by metastable greigite than those displayed by chalcopyrite and cubanite, consistent with a proposal that the Fe atom in greigite is tetravalent. S-S bond paths exist between each of the surface S atoms of adjacent slabs of FeS6 octahedra comprising the layer sulfide smythite, suggesting that the neutral Fe3S4 slabs are linked together and stabilized by the pathways of electron density comprising S-S bonded interactions. Such interactions not only exist between the S atoms for adjacent S8 rings in native sulfur, but their bond critical point properties are similar to those displayed by the metal sulfides.