Density-functional calculation of the electronic structure and equilibrium geometry of iron pyrite (FeS2)

Development of molecular cluster models to probe pyrite surface reactivity

The recent discovery that anaerobic methanogens can reductively dissolve pyrite and utilize dissolution products as a source of iron and sulfur to meet their biosynthetic demands for these elements prompted the development of atomic-scale nanoparticle models, as maquettes of reactive surface sites, for describing the fundamental redox steps that take place at the mineral surface during reduction. The given report describes our computational approach for modeling n(FeS2 ) nanoparticles originated from mineral bulk structure. These maquettes contain a comprehensive set of coordinatively unsaturated Fe(II) sites that are connected via a range of persulfide (S2 2- ) ligation. In addition to the specific maquettes with n = 8, 18, and 32 FeS2 units, we established guidelines for obtaining low-energy structures by considering the pattern of ionic, covalent, and magnetic interactions among the metal and ligand sites. The developed models serve as computational nano-reactors that can be used to describe the reductive dissolution mechanism of pyrite to better understand the reactive sites on the mineral, where microbial extracellular electron-transfer reactions can occur.

Effects of sulfur-deficient defect and water on rearrangements of formamide on pyrite (100) surface

The efficient formation of HCN/HNC from formamide (FM) combining the advantages of water-assistance, self-catalyzed reactions, and the mineral surfaces was investigated. Periodic density functional theory calculations with plane-wave pseudopotential basis sets were performed to study the interaction of FM with pyrite (100) ideal and defect surfaces. Effects of sulfur vacancy defect and water on tautomerization and rearrangement barriers of FM on the (100) surface were evaluated. Calculated results show that FM adsorbs more strongly on the defect surface than on the ideal surface, with the lowest adsorption energy on the defect surface being -22 kcal/mol. The energy barriers for rearrangements of FM on these two surfaces being close to each other suggests that the adsorptions on the surfaces have small effects on the energy barriers. The energy barriers for formimic acid isomer formations are 44.5 and 46.0 kcal/mol, and those of aminohydroxymethylene formations are 72.6 and 71.9 kcal/mol on the ideal and defect surfaces, respectively. A reduction of ∼30 kcal/mol in tautomerization energy barriers is observed in water-assisted process on the defect surface. Because this reduction is close to that of the gas-phase reactions, the catalytic effect is clearly due to the presence of water molecule instead of the interaction with the surface. In this case, the pyrite surfaces with the ability to accumulate reactive species only play the role of connecting bridges between the two steps of the proposed reaction mechanism: the water-assisted rearrangement and the self-catalyzed dehydration.

Atomic properties and chemical bonding in the pyrite and marcasite polymorphs of FeS2: a combined experimental and theoretical electron density study

The chemical bonding in the pyrite (left) and marcasite (right) polymorphs of FeS₂is investigated by charge density analysis.

Electronic Structure, Pressure Dependence and Optical Properties of FeS2

A revisited electronic structure study of iron pyrite, FeS2, has been performed using a new Tight-Binding Linear Muffin-Tin Orbital (TB-LMTO) technique in which the radii of overlapping MT spheres are determined from a full potential construction. The interstitial spheres were chosen to provide an efficient packing of space while ensuring that the overlap between the spheres remain small. We have found that this treatment of interstitial spheres results in a dramatic improvement in the description of the electronic structure and the binding energy curves for FeS2 in comparison with a previous LMTO calculation. In particular, the energy band gap, the equilibrium lattice constant and the bulk modulus are all in much better agreement with experimental observations. Moreover, the calculated equation of state is in excellent accord with recent measured P- V data up to pressures of 15GPa with overall deviations of less than 10%. The predicted reflectivity spectrum of FeS2 as a function of pressure gives the observed behaviour of the optical edge. The bonding behaviour the orthorhombic marcasite phase of FeS2 is also discussed within this new TB-LMTO formalism.

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.

Glycine at the Pyrite−Water Interface: The Role of Surface Defects

Ab initio molecular dynamics simulations were performed in order to study chemisorption, electronic properties, and desorption of glycine at wet pyrite surfaces focusing on the role of surface point defects. The change in the electronic structure and its influence on the chemical reactivity of the free FeS(2)(100) surface due to sulfur vacancies was studied in detail yielding several adsorption modes of glycine and water molecules. Energetically preferred adsorption modes were furthermore investigated in the presence of hot pressurized water mimicking "Iron Sulfur World" prebiotic conditions. The metadynamics Car-Parrinello technique was employed to map the free energy landscape including paths and barriers for desorption of glycine from such wet defective surfaces. The ubiquitous sulfur vacancies are found to increase the retention time of the adsorbed amino acid by many orders of magnitudes in comparison to the ideal pyrite-water interface. The importance of these findings in terms of a possible two-dimensional primordial chemistry on mineral surfaces is discussed.

New Solid−Gas Metathetical Synthesis of Binary Metal Polysulfides and Sulfides at Intermediate Temperatures: Utilization of Boron Sulfides

A new simple synthetic method for binary metal polysulfides and sulfides was developed by utilizing an in situ formation of boron sulfides and their subsequent reactions with metal-source oxides in a closed container at intermediate temperatures above 350 degrees C at which the boron sulfides react in a gaseous form. The versatility of the new method is demonstrated with oxides of various transition metals (Ti, V, Mn, Fe, Ni, Nb, Mo, Ru, and W) and rare-earth metals (Y, Ce, Nd, Sm, Eu, Tb, and Er) as starting materials that exhibit different chemical characteristics. Regardless of the oxidation states of metals in the starting materials, the sulfidation reactions occurred quantitatively with stoichiometric mixtures of boron and sulfur, and within 24 h the reactions yielded pure products of TiS(2), TiS(3), VS(4), FeS(2), NiS(2), NbS(3), MoS(2), RuS(2), WS(2), Y(2)S(3), and RS(2) (R = Ce, Nd, Sm, Eu, Tb, and Er) which were the thermodynamically stable phases under the reaction conditions. The scope and implications of the new sulfidation method are also discussed.