Methylmercury is the environmental form of neurotoxic mercury that is biomagnified in the food chain. Methylation rates are reduced when the metal is sequestered in crystalline mercury sulfides or bound to thiol groups in macromolecular natural organic matter. Mercury sulfide minerals are known to nucleate in anoxic zones, by reaction of the thiol-bound mercury with biogenic sulfide, but not in oxic environments. We present experimental evidence that mercury sulfide forms from thiol-bound mercury alone in aqueous dark systems in contact with air. The maximum amount of nanoparticulate mercury sulfide relative to thiol-bound mercury obtained by reacting dissolved mercury and soil organic matter matches that detected in the organic horizon of a contaminated soil situated downstream from Oak Ridge, TN, in the United States. The nearly identical ratios of the two forms of mercury in field and experimental systems suggest a common reaction mechanism for nucleating the mineral. We identified a chemical reaction mechanism that is thermodynamically favorable in which thiol-bound mercury polymerizes to mercury-sulfur clusters. The clusters form by elimination of sulfur from the thiol complexes via breaking of mercury-sulfur bonds as in an alkylation reaction. Addition of sulfide is not required. This nucleation mechanism provides one explanation for how mercury may be immobilized, and eventually sequestered, in oxygenated surface environments.

We present results obtained from high energy-resolution L3-edge XANES spectroscopy and first-principles calculations for the structure, bonding, and stability of mercury(II) complexes with thiolate and thioether ligands in crystalline compounds, aqueous solution, and macromolecular natural organic matter (NOM). Core-to-valence XANES features that vary in intensity differentiate with unprecedented sensitivity the number and identity of Hg ligands and the geometry of the ligand environment. Post-Hartree-Fock XANES calculations, coupled with natural population analysis, performed on MP2-optimized Hg[(SR)2···(RSR)n] complexes show that the shape, position, and number of electronic transitions observed at high energy-resolution are directly correlated to the Hg and S (l,m)-projected empty densities of states and occupations of the hybridized Hg 6s and 5d valence orbitals. Linear two-coordination, the most common coordination geometry in mercury chemistry, yields a sharp 2p to 6s + 5d electronic transition. This transition varies in intensity for Hg bonded to thiol groups in macromolecular NOM. The intensity variation is explained by contributions from next-nearest, low-charge, thioether-type RSR ligands at 3.0-3.3 Å from Hg. Thus, Hg in NOM has two strong bonds to thiol S and k additional weak Hg···S contacts, or 2 + k coordination. The calculated stabilization energy is -5 kcal/mol per RSR ligand. Detection of distant ligands beyond the first coordination shell requires precise measurement of, and comparison to, spectra of reference compounds as well as accurate calculation of spectra for representative molecular models. The combined experimental and theoretical approaches described here for Hg can be applied to other closed-shell atoms, such as Ag(I) and Au(I). To facilitate further calculation of XANES spectra, experimental data, a new crystallographic structure of a key mercury thioether complex, Cartesian coordinates of the computed models, and examples of input files are provided as Supporting Information .