Monazite-(Ce) and xenotime-(Y) occur as secondary minerals in iron-oxide-apatite (IOA) deposits, and their stability and composition are important indicators of timing and conditions of metasomatism. Both of these minerals occur as replacement of apatite and display slight but important variations in light (e.g. La, Ce, Pr, Nd, etc.) and heavy (e.g. Y, Er, Dy, Yb, etc.) REE concentrations [1,2]. The causes for these chemical variations can be quantified by combining thermodynamic modeling with field observations. Major challenges for determining the stability of these minerals in hydrothermal solutions are the underlying models for calculating the thermodynamic properties of REE-bearing mineral solid solutions and aqueous species as a function of temperature and pressure. The thermodynamic properties of monazite and xenotime have been determined using several calorimetric methods [3], but only a few hydrothermal solubility studies have been undertaken, which test the reliability and compatibility of both the calorimetric data and thermodynamic properties of associated REE aqueous species [4,5]. Here, we evaluate the conditions of REE metasomatism in the Pea Ridge IOA-REE deposit in Missouri, and combine newly available experimental solubility data to simulate the speciation of LREE vs. HREE, and the partitioning of REE as a function of varying fluid compositions and temperatures. Our new experimental data will be implemented in the MINES thermodynamic database (http:// tdb.mines.edu) for modeling the chemistry of crustal fluid-rock equilibria [6]. [1] Harlov et al. (2016), Econ. Geol. 111, 1963-1984;[2] Hofstra et al. (2016), Econ. Geol. 111, 1985-2016; [3] Navrotsky et al. (2015), J. Chem. Thermodyn. 88, 126-141; [4] Gysi et al. (2015), Chem. Geol. 83-95; [5] Gysi et al. (2018), Geochim. Cosmochim. Acta 242, 143-164; [6] Gysi (2017), Pure and Appl. Chem. 89, 581-596.
more »
« less
Reaction calorimetry and structural crystal properties of non-ideal binary rhabdophane solid solutions (Ce1−xREExPO4·nH2O)
Rhabdophane is a hydrous phosphate that commonly replaces monazite as a weathering product in critical mineral deposits during the alteration of rare earth elements (REE) bearing carbonatites and alkaline igneous complexes. It is an important host to the light (L)REE (i.e., La to Gd) but the stability and structure of binary solid solutions between the Ce and the other LREE endmembers have not yet been determined experimentally. Here we present room temperature calorimetric experiments that were used to measure the enthalpy of precipitation of rhabdophane (Ce1−xREExPO4·nH2O; REE = La, Pr, Nd, Sm, Eu, and Gd). The solids were characterized using X-ray diffraction, scanning electron microscopy, Raman spectroscopy, and the role of water in the rhabdophane structure was further determined using thermogravimetric analysis coupled with differential scanning calorimetry. The calorimetric experiments indicate a non-ideal behavior for all of the binary solid solutions investigated with an excess enthalpy of mixing (ΔHex) described by a 2- to 3-term Guggenheim parameters equation. The solid solutions were categorized into three groups: (1) binary Ce-La and Ce-Pr which display positive ΔHex values with a slight asymmetry; (2) binary Ce-Nd and Ce-Sm which display negative ΔHex values with a nearly symmetric shape; (3) Ce-Eu and Ce-Gd which display both negative and positive ΔHex values with nearly symmetric shape. The excess Gibbs energy (ΔGex) of the solid solutions was further investigated using a thermodynamic analysis approach of aqueous-solid solution equilibria and the optimization programs GEMS and GEMSFITS. The resulting ΔGex values combined with the calorimetric ΔHex values indicate that there is likely an excess entropy contribution implying important short-range structural modifications in the solid solutions dependent on the deviation of the REE ionic radii from the size of Ce3+. These observations corroborate with the trends in the Raman v1 stretching bands of the PO4-site. The excess molar volumes determined from X-ray diffraction analysis further indicate an overall asymmetric behavior in all of the studied binary solid solutions, which becomes increasingly important from La to Gd. The pronounced short-range order–disorder occurring in groups 2 and 3 solid solutions mimics some of the behavior observed from previous studies in anhydrous monazite solid solutions. This study highlights the potential to use the chemistry and the structural modifications of rhabdophane as potential indicators of formation conditions in geologic systems and permits improving our modeling capabilities of REE partitioning in critical minerals systems.
more »
« less
- PAR ID:
- 10506366
- Publisher / Repository:
- Mendeley Data
- Date Published:
- Journal Name:
- Geochimica et Cosmochimica Acta
- ISSN:
- 0016-7037
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Charlier, B (Ed.)Quantifying the oxygen fugacity (fo2) of high temperature lithospheric fluids, including hydrothermal systems, presents a challenge because these fluids are difficult to capture and measure in the same manner as quenched glasses of silicate melts. The chemical properties of fluids can however be inferred through mineral proxies that interacted with the fluids through precipitation or recrystallization. Here, we present hydrothermal experiments to quantify the partition coefficients of rare earth elements (REEs) – including redox-sensitive Ce and Eu – between zircon and fluid. Experiments were conducted in a piston cylinder device at temperatures that range from 1200 to 800 ◦C under fo2-buffered conditions in a SiO2-ZrO2-NaCl-REE-oxide system, and similar experiments were performed in the absence of NaCl (31 total experiments). The fo2 was buffered to values that range from approximately 3 log units below to 7 log units above the fayalite magnetite quartz equilibrium. Zircon REE concentrations were quantified using laser ablation inductively coupled plasma mass spectrometry whereas the quenched fluids were extracted and measured by solution-based inductively coupled plasma mass spectrometry. Zircon Ce anomalies, quantified relative to La and Pr, exhibit sensitivity to oxygen fugacity and temperature and our preferred calibration is: log [ Ce Ce* ) D 1 ] = (0.237 ± 0.040)× log(fo2) + 9437±640 T(K) 5.02 ± 0.38 where the Ce anomalies are calculated from the partition coefficients for La, Ce, and Pr. Zircon Eu anomalies are also a function of oxygen fugacity though they exhibit no systematic dependence on T. Our preferred calibration is described by: Eu Eu* ) D = 1 1+100.30±0.04 [0.27±0.03]×ΔFMQ We performed additional calculations, in which lattice strain parabolas were fit to all non-redox sensitive rare earth elements that were added to the starting composition (i.e., La, Pr, Sm, Gd, Dy, Ho, Tm, Lu) as an alternate means to calculate anomalies. This method yields broadly similar results, though we prefer the La-Pr calibrations due to the non-systematic REE patterns frequently encountered with hydrothermal zircons; e.g., LREE zircon enrichment relative to other REEs. These experiments are applied to quantify the fo2 of fluids during mineralization of critical element-bearing systems, and separately to calculate the oxygen fugacity values of fluids formed during plate boundary processes.more » « less
-
Monazite is a light rare earth element (REE) phosphate found in REE mineral deposits, such as those formed in (per)alkaline and carbonatite magmatic-hydrothermal systems, where it occurs in association to the development of alteration zones and hydrothermal veins. Although it has been recognized that monazite may undergo replacement by coupled dissolution-precipitation processes, currently there is no model describing the compositional REE variations in monazite resulting from direct interaction with or precipitation from hydrothermal fluids. To develop such a model requires quantification of the thermodynamic properties of the aqueous REE species and the properties of the monazite endmembers and their solid solutions. The thermodynamic properties of monazite endmembers have been determined previously using calorimetric methods and low temperature solubility studies, but only a few solubility studies have been conducted at >100 °C. In this study, the solubility products (logKs0) of LaPO4, PrPO4, NdPO4, and EuPO4 monazite endmembers have been measured at temperatures between 100 and 250 °C and saturated water vapor pressure. The solubility products are reported with an uncertainty of ±0.2 (95% confidence) according to the reaction, REEPO4(s) = REE3+ + PO43−. (see table in manuscript) The REE phosphates display a retrograde solubility, with the measured Ks0 values varying several orders of magnitude over the experimental temperature range. Discrepancies were observed between the experimental solubility products and the calculated values resulting from combining calorimetric data of monazite with the properties of the aqueous REE3+ and PO43− species available in the literature. The differences between the calculated and measured standard Gibbs energy of reaction (ΔrG0) for PrPO4, NdPO4, and EuPO4 increased with higher temperatures (up to 15 kJ mol−1 at 250 °C), whereas for LaPO4 these differences increased at lower temperatures (up to 8 kJ mol−1 at 100 °C). To reconcile these discrepancies, the standard enthalpy of formation (ΔfH0) of monazite was optimized by fitting the experimental solubility data and extrapolating these fits to reference conditions of 25 °C and 1 bar. The optimized thermodynamic data provide the first internally consistent dataset for the solubility of all the monazite endmembers, and can be used to model REE partitioning between monazite and hydrothermal fluids at >100 °C.more » « less
-
Abstract Phase formation in multicomponent rare‐earth oxides is determined by a combination of composition, sintering atmosphere, and cooling rate. Polycrystalline ceramics comprising various combinations of Ce, Gd, La, Nd, Pr, Sm, and Y oxides in equiatomic proportions were synthesized using solid‐state sintering. The effects of composition, sintering atmosphere, and cooling rate on phase formation were investigated. Single cubic or monoclinic structures were obtained with a slow cooling of 3.3°C/min, confirming that rare‐earth oxides follow a different structure stabilization process than transition metal high‐entropy oxides. In an oxidizing atmosphere, both Ce and Pr induce a cubic structure, while only Ce plays that role in an inert or reducing atmosphere. Samples without Ce or Pr develop a single monoclinic structure. The structures formed at initial synthesis may be converted to a different one, when the ceramics are annealed in an additional atmosphere. Phase evolution of a five‐cation composition was also studied as a function of sintering temperature. The binary oxides used as raw materials completely dissolve into a single cubic structure at 1450°C in air.more » « less
-
The reduction potentials (reported vs. Fc + /Fc) for a series of Cp′ 3 Ln complexes (Cp′ = C 5 H 4 SiMe 3 , Ln = lanthanide) were determined via electrochemistry in THF with [ n Bu 4 N][BPh 4 ] as the supporting electrolyte. The Ln( iii )/Ln( ii ) reduction potentials for Ln = Eu, Yb, Sm, and Tm (−1.07 to −2.83 V) follow the expected trend for stability of 4f 7 , 4f 14 , 4f 6 , and 4f 13 Ln( ii ) ions, respectively. The reduction potentials for Ln = Pr, Nd, Gd, Tb, Dy, Ho, Er, and Lu, that form 4f n 5d 1 Ln( ii ) ions ( n = 2–14), fall in a narrow range of −2.95 V to −3.14 V. Only cathodic events were observed for La and Ce at −3.36 V and −3.43 V, respectively. The reduction potentials of the Ln( ii ) compounds [K(2.2.2-cryptand)][Cp′ 3 Ln] (Ln = Pr, Sm, Eu) match those of the Cp′ 3 Ln complexes. The reduction potentials of nine (C 5 Me 4 H) 3 Ln complexes were also studied and found to be 0.05–0.24 V more negative than those of the Cp′ 3 Ln compounds.more » « less
An official website of the United States government

