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, F(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, F(rc), the Laplacian, 32F(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 F(rc) and 32F(rc) are larger for pyrite and marcasite, indicating that the accumulation and local concentration of F(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 F(rc) value is large