Arsenic in hydrothermal apatite : 2 oxidation state , mechanism of uptake , and comparison 3 between experiments and nature 4 5 6

23 Element substitution that occurs during fluid-rock interaction permits assessment of fluid 24 composition and interaction conditions in ancient geological systems, and provides a way to fix 25 contaminants from aqueous solutions. We conducted a series of hydrothermal mineral replacement 26 experiments to determine whether a relationship can be established between arsenic (As) 27 distribution in apatite and fluid chemistry. Calcite crystals were reacted with phosphate solutions 28 spiked with As(V), As(III), and mixed As(III)/As(V) species at 250 ̊C and water-saturated 29 pressure. Arsenic-bearing apatite rims formed in several hours, and within 48 hours the calcite 30 grains were fully replaced. X-ray Absorption Near-edge Spectroscopy (XANES) data show that 31 As retained the trivalent oxidation state in the fully-reacted apatite grown from solutions 32 containing only As(III). Extended X-ray Fine Spectroscopy (EXAFS) data reveal that these As(III) 33 ions are surrounded by about three oxygen atoms at an As-O bond length close to that of an 34 arsenate group (AsO4 ), indicating that they occupy tetrahedral phosphate sites. The three35 coordinated As(III)-O3 structure, with three oxygen atoms and one lone electron pair around 36 As(III), was confirmed by geometry optimization using ab initio molecular simulations. 37 The micro-XANES imaging data show that apatite formed from solutions spiked with 38 mixed As(III) and As(V) retained only As(V) after completion of the replacement reaction; in 39 contrast, partially reacted samples revealed a complex distribution of As(V)/As(III) ratios, with 40 As(V) concentrated in the center of the grain and As(III) towards the rim. Most natural apatites 41 from the Ernest Henry Iron Oxide Copper Gold deposit, Australia, show predominantly As(V), 42 but two grains retained some As(III) in their core. The As-anomalous amphibolite-facies gneiss 43 from Binntal, Switzerland, only revealed As(V), despite the fact that these apatites in both cases 44 formed under conditions where As(III) is expected to be the dominant As form in hydrothermal 45 fluids. 46 These results show that incorporation of As in apatite is a complicated process, and 47 sensitive to the local fluid composition during crystallization, and that some of the complexity in 48 As zoning in partially reacted apatite may be due to local fluctuations of As(V)/As(III) ratios in 49 the fluid and to kinetic effects during the mineral replacement reaction. Our study shows for the 50 first time that As(III) can be incorporated into the apatite structure, although not as efficiently as 51 As(V). Uptake of As(III) is probably highly dependent on the reaction mechanism. As(III)O3 352 moieties replace phosphate groups, but cause a high strain on the lattice; as a result, As(III) is 53 easily exchanged (or oxidized) for As(V) during hydrothermal recrystallization, and the fully 54 Arsenic in apatite Page 3 reacted grains only record the preferred oxidation state (i.e., As(V)) from mixed-oxidation state 55 solutions. Overall this study shows that the observed oxidation state of As in apatite may not reflect 56 the original As(III)/As(V) ratio of the parent fluid, due to the complex nature of As(III) uptake and 57 possible in-situ oxidation during recrystallization. 58 59

[1]  L. Coogan,et al.  Apatite Trace Element Compositions: A Robust New Tool for Mineral Exploration , 2016 .

[2]  J. Brugger,et al.  Textural and compositional complexities resulting from coupled dissolution–reprecipitation reactions in geomaterials , 2015 .

[3]  J. Brugger,et al.  Palladium complexation in chloride- and bisulfide-rich fluids: Insights from ab initio molecular dynamics simulations and X-ray absorption spectroscopy , 2015 .

[4]  John M. Hughes,et al.  Structurally Robust, Chemically Diverse: Apatite and Apatite Supergroup Minerals , 2015 .

[5]  T. Cheng,et al.  Silicate Minerals as a Source of Arsenic Contamination in Groundwater , 2014, Water, Air, & Soil Pollution.

[6]  C. Ryan,et al.  Speciation mapping of environmental samples using XANES imaging , 2014 .

[7]  C. MacRae,et al.  Complex mineral zoning patterns caused by ultra-local equilibrium at reaction interfaces , 2014 .

[8]  J. Hazemann,et al.  Speciation and thermodynamic properties of manganese(II) chloride complexes in hydrothermal fluids: In situ XAS study , 2014 .

[9]  F. di Benedetto,et al.  Arsenic-bearing calcite in natural travertines: evidence from sequential extraction, μXAS, and μXRF. , 2013, Environmental science & technology.

[10]  Kazuya Tanaka,et al.  Differences in the immobilization of arsenite and arsenate by calcite , 2012 .

[11]  A. Berger,et al.  Allanite behaviour during incipient melting in the southern Central Alps , 2012 .

[12]  S. Tofail,et al.  Unravelling the specific site preference in doping of calcium hydroxyapatite with strontium from ab initio investigations and Rietveld analyses. , 2012, Physical chemistry chemical physics : PCCP.

[13]  L. Valenzano,et al.  Arsenic uptake by natural calcite: An XAS study , 2011 .

[14]  A. Williams-Jones,et al.  Reading the mineral record of fluid composition from element partitioning , 2010 .

[15]  D. P. Siddons,et al.  The Maia Spectroscopy Detector System: Engineering for Integrated Pulse Capture, Low-Latency Scanning and Real-Time Processing , 2010 .

[16]  D. P. Siddons,et al.  The new Maia detector system : Methods for High Definition Trace Element Imaging of natural material , 2010 .

[17]  A. Slepko First principles study of biomineral hydroxyapatite , 2010 .

[18]  Yinian Zhu,et al.  Synthesis and characterization of arsenate/phosphate hydroxyapatite solid solution , 2009 .

[19]  J. Hazemann,et al.  Metal speciation in fluid inclusions using microbeam X-ray absorption spectroscopy , 2009 .

[20]  Yu. V. Shvarov,et al.  HCh: New potentialities for the thermodynamic simulation of geochemical systems offered by windows , 2008 .

[21]  Michael Bühl,et al.  Geometries of Transition-Metal Complexes from Density-Functional Theory. , 2006, Journal of chemical theory and computation.

[22]  N. Oliver,et al.  Comparing closed system, flow‐through and fluid infiltration geochemical modelling: examples from K‐alteration in the Ernest Henry Fe‐oxide–Cu–Au system , 2005 .

[23]  T. White,et al.  Structural derivation and crystal chemistry of apatites. , 2003, Acta crystallographica. Section B, Structural science.

[24]  R. Bateman,et al.  Origins of Nd–Sr–Pb isotopic variations in single scheelite grains from Archaean gold deposits, Western Australia , 2002 .

[25]  S. Sutton,et al.  Evidence of heterovalent europium in zoned Llallagua apatite using wavelength dispersive XANES , 2001 .

[26]  P. Brown,et al.  Arsenic immobilization by calcium arsenate formation , 1999 .

[27]  J.C. Elliott,et al.  Rietveld refinement of the crystallographic structure of human dental enamel apatites , 1999 .

[28]  J. Brugger,et al.  As, Sb, Be and Ce enrichment in minerals from a metamorphosed Fe-Mn deposit, Val Ferrera, eastern Swiss Alps , 1999 .

[29]  David N. Jamieson,et al.  Dynamic analysis: on-line quantitative PIXE microanalysis and its use in overlap-resolved elemental mapping , 1993 .

[30]  Martins,et al.  Efficient pseudopotentials for plane-wave calculations. , 1991, Physical review. B, Condensed matter.

[31]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[32]  Everett L. Shock,et al.  Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Correlation algorithms for ionic species and equation of state predictions to 5 kb and 1000°C , 1988 .

[33]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[34]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[35]  I. D. Brown,et al.  Bond‐valence parameters obtained from a systematic analysis of the Inorganic Crystal Structure Database , 1985 .

[36]  M. Catti,et al.  On the topotactic dehydration Ca3(AsO4)2.11H2O (phaunouxite) → Ca3(AsO4)2.10H2O (rauenthalite), and the structures of both minerals , 1983 .

[37]  R. D. Shannon Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides , 1976 .

[38]  G. Chiari,et al.  The crystal structure of Na2HAsO4.7H2O , 1970 .

[39]  W. Skinner,et al.  Formation of As(II)-pyrite during experimental replacement of magnetite under hydrothermal conditions , 2013 .

[40]  A. Putnis Mineral Replacement Reactions , 2009 .

[41]  P. Ballirano,et al.  Refinement of the crystal structure of arsenolite, AS2O3 , 2002 .

[42]  M. Fleet,et al.  Compositions of the Apatite-Group Minerals: Substitution Mechanisms and Controlling Factors , 2002 .

[43]  C. G. Ryan,et al.  Quantitative trace element imaging using PIXE and the nuclear microprobe , 2000, Int. J. Imaging Syst. Technol..

[44]  J. Brugger Inhomogeneous Distribution of REE and Initial Nd Isotopic Compositions in Scheelite Crystals from an Archaean Gold Deposit (Mt Charlotte, Western Australia) , 1998 .

[45]  R. Morton,et al.  Elemental X-Ray Imaging of Fossils , 1992 .

[46]  W. L. Jolly A Modern Inorganic Chemistry , 2022, Nature.