Yin-Chuan Li , Hao-Long Li , Hai-Zhen Wei , Martin R. Palmer , Yong-Hui Li , Ya-Ru Zhou
{"title":"Potassium isotopic equilibrium fractionation induced by clay adsorption and incorporation: A density functional theory study","authors":"Yin-Chuan Li , Hao-Long Li , Hai-Zhen Wei , Martin R. Palmer , Yong-Hui Li , Ya-Ru Zhou","doi":"10.1016/j.gca.2025.08.026","DOIUrl":null,"url":null,"abstract":"<div><div>Chemical weathering of continental silicate rocks is a key process in controlling the concentration and isotope composition of many elements in seawater, including potassium. K isotopes show great potential for tracing the silicate weathering. However, their application remains constrained by a lack of understanding of K isotopic fractionation mechanisms during chemical weathering, wherein secondary clay minerals adsorb and incorporate K derived from silicate minerals. Using first-principles molecular dynamics simulations and density functional theory (DFT) calculation, the coordination states of K in aqueous fluid and at the kaolinite (001) and (010) surfaces, as well as the K isotopic equilibrium fractionation caused by the adsorption of kaolinite surface and the incorporation of illite and glauconite, have been investigated. The results show that five (<sup>V</sup>K) and six (<sup>VI</sup>K) coordinated K dominate in aqueous fluids (∼88 %) and on kaolinite (001) (∼84 %), while <sup>IV</sup>K and <sup>V</sup>K prevail on (010) surfaces (∼80 %). The 1000lnβ values of aqueous fluids and kaolinite (001) and (010) surface configurations do not correlate well with structural factors (e.g., K–O bond lengths, the distortions of K–O bond length and O–K–O angle of coordination polyhedral). Rather, the 1000lnβ values exhibit a strong linear correlation with the force constants acting on the K<sup>+</sup> ion in minerals and fluids. This process results in K isotopic equilibrium fractionation between the kaolinite (001) and (010) surfaces and aqueous fluids of 0.63 ± 0.37 ‰ and 0.58 ± 0.42 ‰, respectively, at 25 °C. The tetrahedral Al absence alters the K isotope equilibrium fractionation between clay minerals and aqueous fluids (0.25 ± 0.30 ‰ vs. −0.37 ± 0.24 ‰ for illite; −0.06 ± 0.22 ‰ vs. −0.35 ± 0.23 ‰; for glauconite; at 25 °C). These theoretical calculations may be used to model K isotope evolution in surface weathering environments to gain a more detailed understanding of K isotope fractionation under varying conditions of temperature, pH, ionic strength, and salinity.</div></div>","PeriodicalId":327,"journal":{"name":"Geochimica et Cosmochimica Acta","volume":"407 ","pages":"Pages 32-46"},"PeriodicalIF":5.0000,"publicationDate":"2025-08-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Geochimica et Cosmochimica Acta","FirstCategoryId":"89","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0016703725004302","RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"GEOCHEMISTRY & GEOPHYSICS","Score":null,"Total":0}
引用次数: 0
Abstract
Chemical weathering of continental silicate rocks is a key process in controlling the concentration and isotope composition of many elements in seawater, including potassium. K isotopes show great potential for tracing the silicate weathering. However, their application remains constrained by a lack of understanding of K isotopic fractionation mechanisms during chemical weathering, wherein secondary clay minerals adsorb and incorporate K derived from silicate minerals. Using first-principles molecular dynamics simulations and density functional theory (DFT) calculation, the coordination states of K in aqueous fluid and at the kaolinite (001) and (010) surfaces, as well as the K isotopic equilibrium fractionation caused by the adsorption of kaolinite surface and the incorporation of illite and glauconite, have been investigated. The results show that five (VK) and six (VIK) coordinated K dominate in aqueous fluids (∼88 %) and on kaolinite (001) (∼84 %), while IVK and VK prevail on (010) surfaces (∼80 %). The 1000lnβ values of aqueous fluids and kaolinite (001) and (010) surface configurations do not correlate well with structural factors (e.g., K–O bond lengths, the distortions of K–O bond length and O–K–O angle of coordination polyhedral). Rather, the 1000lnβ values exhibit a strong linear correlation with the force constants acting on the K+ ion in minerals and fluids. This process results in K isotopic equilibrium fractionation between the kaolinite (001) and (010) surfaces and aqueous fluids of 0.63 ± 0.37 ‰ and 0.58 ± 0.42 ‰, respectively, at 25 °C. The tetrahedral Al absence alters the K isotope equilibrium fractionation between clay minerals and aqueous fluids (0.25 ± 0.30 ‰ vs. −0.37 ± 0.24 ‰ for illite; −0.06 ± 0.22 ‰ vs. −0.35 ± 0.23 ‰; for glauconite; at 25 °C). These theoretical calculations may be used to model K isotope evolution in surface weathering environments to gain a more detailed understanding of K isotope fractionation under varying conditions of temperature, pH, ionic strength, and salinity.
期刊介绍:
Geochimica et Cosmochimica Acta publishes research papers in a wide range of subjects in terrestrial geochemistry, meteoritics, and planetary geochemistry. The scope of the journal includes:
1). Physical chemistry of gases, aqueous solutions, glasses, and crystalline solids
2). Igneous and metamorphic petrology
3). Chemical processes in the atmosphere, hydrosphere, biosphere, and lithosphere of the Earth
4). Organic geochemistry
5). Isotope geochemistry
6). Meteoritics and meteorite impacts
7). Lunar science; and
8). Planetary geochemistry.