Direct prediction of isotopic properties from molecular dynamics trajectories: Application to sulfide, sulfate and sulfur radical ions in hydrothermal fluids
Marc Blanchard , Elsa Desmaele , Gleb S. Pokrovski , Carlos Pinilla , Merlin Méheut , Rodolphe Vuilleumier
{"title":"Direct prediction of isotopic properties from molecular dynamics trajectories: Application to sulfide, sulfate and sulfur radical ions in hydrothermal fluids","authors":"Marc Blanchard , Elsa Desmaele , Gleb S. Pokrovski , Carlos Pinilla , Merlin Méheut , Rodolphe Vuilleumier","doi":"10.1016/j.chemgeo.2024.122202","DOIUrl":null,"url":null,"abstract":"<div><p>In hydrothermal fluids, disulfur (S<sub>2</sub><sup>•−</sup>) and trisulfur (S<sub>3</sub><sup>•−</sup>) radical anions have been observed to coexist with the major hydrogen sulfide and sulfate species. These radical ions have potentially important effects on the solubility, transport and fractionation of metals and sulfur, with consequences for ore deposit formation and, more generally, for geochemical cycles of metals and volatiles. It is therefore essential to know the intrinsic isotopic properties of these important sulfur species in order to use sulfur isotopes for tracing different geological processes. Here, the theoretical equilibrium isotopic properties of the disulfur and trisulfur radical ions are computed and compared to the hydrogen sulfide (H<sub>2</sub>S) and sulfate ion (SO<sub>4</sub><sup>2−</sup>), using, for the first time, a first-principles molecular dynamics (MD) approach. The isotopic properties are calculated directly from molecular dynamics trajectories using the vibrational density of states and the atomic kinetic energy, and then compared to the more established method based on sampling of several snapshots. This comparison allowed us to validate the new modelling method and to assess its advantages and limitations. The predicted equilibrium isotope fractionation in terms of <sup>34</sup>S/<sup>32</sup>S between S<sub>2</sub><sup>•−</sup> and S<sub>3</sub><sup>•−</sup> is small, i.e. <1‰, with a slight enrichment in the heavier isotope for S<sub>3</sub><sup>•−</sup>, over the temperature range 200–500 °C. Both radical ions are slightly depleted in the heavier isotope, by 1 to 2‰, relative to aqueous H<sub>2</sub>S. Our results help tuning sulfur isotope fractionation models used for tracing the origin and evolution of hydrothermal fluids. Our method opens large perspectives for using the rapidly growing body of MD simulation data in geosciences on structure and stability of aqueous complexes to assess in parallel element isotope fractionations from MD-generated trajectories.</p></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":null,"pages":null},"PeriodicalIF":3.6000,"publicationDate":"2024-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0009254124002821/pdfft?md5=080f8504403a7fea82aae9b7a1c5c368&pid=1-s2.0-S0009254124002821-main.pdf","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Chemical Geology","FirstCategoryId":"89","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0009254124002821","RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"GEOCHEMISTRY & GEOPHYSICS","Score":null,"Total":0}
引用次数: 0
Abstract
In hydrothermal fluids, disulfur (S2•−) and trisulfur (S3•−) radical anions have been observed to coexist with the major hydrogen sulfide and sulfate species. These radical ions have potentially important effects on the solubility, transport and fractionation of metals and sulfur, with consequences for ore deposit formation and, more generally, for geochemical cycles of metals and volatiles. It is therefore essential to know the intrinsic isotopic properties of these important sulfur species in order to use sulfur isotopes for tracing different geological processes. Here, the theoretical equilibrium isotopic properties of the disulfur and trisulfur radical ions are computed and compared to the hydrogen sulfide (H2S) and sulfate ion (SO42−), using, for the first time, a first-principles molecular dynamics (MD) approach. The isotopic properties are calculated directly from molecular dynamics trajectories using the vibrational density of states and the atomic kinetic energy, and then compared to the more established method based on sampling of several snapshots. This comparison allowed us to validate the new modelling method and to assess its advantages and limitations. The predicted equilibrium isotope fractionation in terms of 34S/32S between S2•− and S3•− is small, i.e. <1‰, with a slight enrichment in the heavier isotope for S3•−, over the temperature range 200–500 °C. Both radical ions are slightly depleted in the heavier isotope, by 1 to 2‰, relative to aqueous H2S. Our results help tuning sulfur isotope fractionation models used for tracing the origin and evolution of hydrothermal fluids. Our method opens large perspectives for using the rapidly growing body of MD simulation data in geosciences on structure and stability of aqueous complexes to assess in parallel element isotope fractionations from MD-generated trajectories.
期刊介绍:
Chemical Geology is an international journal that publishes original research papers on isotopic and elemental geochemistry, geochronology and cosmochemistry.
The Journal focuses on chemical processes in igneous, metamorphic, and sedimentary petrology, low- and high-temperature aqueous solutions, biogeochemistry, the environment and cosmochemistry.
Papers that are field, experimentally, or computationally based are appropriate if they are of broad international interest. The Journal generally does not publish papers that are primarily of regional or local interest, or which are primarily focused on remediation and applied geochemistry.
The Journal also welcomes innovative papers dealing with significant analytical advances that are of wide interest in the community and extend significantly beyond the scope of what would be included in the methods section of a standard research paper.