Uranophane dissolution and growth in CaCl2-SiO2(aq) test solutions

The dissolution and growth of uranophane [Ca(UO₂)₂(SiO₃OH)₂·5H₂O] have been examined in Ca- and Si-rich test solutions at low temperatures (20.5 ± 2.0 °C) and near-neutral pH (∼6.0). Uranium-bearing experimental solutions undersaturated and supersaturated with uranophane were prepared in matrices of ∼10⁻² M CaCl₂ and ∼10⁻³ M SiO₂(aq). The experimental solutions were reacted with synthetic uranophane and analyzed periodically over 10 weeks. Interpretation of the aqueous solution data permitted extraction of a solubility constant for the uranophane dissolution reaction and standard state Gibbs free energy of formation for uranophane ( kJ mol⁻¹).     

Dissolution of uranyl microprecipitates in subsurface sediments at Hanford Site, USA

The dissolution of uranium was investigated from contaminated sediments obtained at the US. Department of Energy (U.S. DOE) Hanford site. The uranium existed in the sediments as uranyl silicate microprecipitates in fractures, cleavages, and cavities within sediment grains. Uranium dissolution was studied in Na, Na-Ca, and NH₄ electrolytes with pH ranging from 7.0 to 9.5 under ambient CO₂ pressure. The rate and extent of uranium dissolution was influenced by uranyl mineral solubility, carbonate concentration, and mass transfer rate from intraparticle regions. Dissolved uranium concentration reached constant values within a month in electrolytes below pH 8.2, whereas concentrations continued to rise for over 200 d at pH 9.0 and above. The steady-state concentrations were consistent with the solubility of Na-boltwoodite and/or uranophane, which exhibit similar solubility under the experimental conditions. The uranium dissolution rate increased with increasing carbonate concentration, and was initially fast. It decreased with time as solubility equilibrium was attained, or dissolution kinetics or mass transfer rate from intraparticle regions became rate-limiting. Microscopic observations indicated that uranium precipitates were distributed in intragrain microfractures with variable sizes and connectivity to particle surfaces. Laser-induced fluorescence spectroscopic change of the uranyl microprecipitates was negligible during the long-term equilibration, indicating that uranyl speciation was not changed by dissolution. A kinetic model that incorporated mineral dissolution kinetics and grain-scale, fracture-matrix diffusion was developed to describe uranium release rate from the sediment. Model calculations indicated that 50-95% of the precipitated uranium was associated with fractures that were in close contact with the aqueous phase. The remainder of the uranium was deeply imbedded in particle interiors and exhibited effective diffusivities that were over three orders of magnitude lower than those in the fractures.     

Fluorescence spectroscopy of U(VI)-silicates and U(VI)-contaminated Hanford sediment

Time-resolved U(VI) laser fluorescence spectra (TRLFS) were recorded for a series of natural uranium-silicate minerals including boltwoodite, uranophane, soddyite, kasolite, sklodowskite, cuprosklodowskite, haiweeite, and weeksite, a synthetic boltwoodite, and four U(VI)-contaminated Hanford vadose zone sediments. Lowering the sample temperature from RT to ∼ 5.5 K significantly enhanced the fluorescence intensity and spectral resolution of both the minerals and sediments, offering improved possibilities for identifying uranyl species in environmental samples. At 5.5 K, all of the uranyl silicates showed unique, well-resolved fluorescence spectra. The symmetric O = U = O stretching frequency, as determined from the peak spacing of the vibronic bands in the emission spectra, were between 705 to 823 cm⁻¹ for the uranyl silicates. These were lower than those reported for uranyl phosphate, carbonate, or oxy-hydroxides. The fluorescence emission spectra of all four sediment samples were similar to each other. Their spectra shifted minimally at different time delays or upon contact with basic Na/Ca-carbonate electrolyte solutions that dissolved up to 60% of the precipitated U(VI) pool. The well-resolved vibronic peaks in the fluorescence spectra of the sediments indicated that the major fluorescence species was a crystalline uranyl mineral phase, while the peak spacing of the vibronic bands pointed to the likely presence of uranyl silicate. Although an exact match was not found between the U(VI) fluorescence spectra of the sediments with that of any individual uranyl silicates, the major spectral characteristics indicated that the sediment U(VI) was a uranophane-type solid (uranophane, boltwoodite) or soddyite, as was concluded from microprobe, EXAFS, and solubility analyses.     

An experimental study of the solubility and speciation of neodymium (III) fluoride in F-bearing aqueous solutions

The solubility of neodymium (III) fluoride was investigated at temperatures of 150, 200 and 250 °C, saturated water vapor pressure, and a total fluoride concentration (HF°_aq + F⁻) ranging from 2.0 × 10⁻³ to 0.23 mol/l. The results of the experiments show that Nd³⁺ and NdF²⁺ are the dominant species in solution at the temperatures investigated and were used to derive formation constants for NdF²⁺ and a solubility product for NdF₃. The solubility product of NdF₃(logK_sp=loga_Nd³⁺+3loga_F^-) is −24.4 ± 0.2, −22.8 ± 0.1, and −21.5 ± 0.2 at 250, 200 and 150 °C, respectively, and the formation constant of NdF²⁺(logβ=loga_NdF²⁺-loga_Nd³⁺-loga_F^-) is 6.8 ± 0.1, 6.2 ± 0.1, and 5.5 ± 0.2 at 250, 200 and 150 °C, respectively. The results of this study show that published theoretical predictions significantly overestimate the stability of NdF²⁺ and the solubility of NdF₃.The potential impact of the results on natural systems was evaluated for a hypothetical fluid with a composition similar to that responsible for REE mineralization in the Capitan pluton, New Mexico. In contrast to results obtained using the theoretical predictions of Haas [Haas J. R., Shock E. L., and Sassani D. C. (1995) Rare earth elements in hydrothermal systems: estimates of standard partial molal thermodynamic properties of aqueous complexes of the rare earth elements at high pressures and temperatures. Geochim. Cosmochim. Acta59, 4329-4350.], which indicate that NdF²⁺ is the dominant species in solution, calculations employing the data presented in this paper and previously published experimental data for chloride and sulfate species [Migdisov A. A., and Williams-Jones A. E. (2002) A spectrophotometric study of neodymium(III) complexation in chloride solutions. Geochim. Cosmochim. Acta66, 4311-4323; Migdisov A. A., Reukov V. V., and Williams-Jones A. E. (2006) A spectrophotometric study of neodymium(III) complexation in sulfate solutions at elevated temperatures. Geochim. Cosmochim. Acta70, 983-992.] show that neodymium chloride species predominate and that neodymium fluoride species are relatively unimportant. This suggests that accepted models for REE deposits that invoke fluoride complexation as the method of hydrothermal REE transport may need to be re-evaluated.     

An isopiestic study of aqueous NaBr and KBr at 50 °C: Chemical equilibrium model of solution behavior and solubility in the NaBr-H2O, KBr-H2O and Na-K-Br-H2O systems to high concentration and temperature

The isopiestic method has been used to determine the osmotic coefficients of the binary solutions NaBr-H₂O (from 0.745 to 5.953 mol kg⁻¹) and KBr-H₂O (from 0.741 to 5.683 mol kg⁻¹) at the temperature t = 50 °C. Sodium chloride solutions have been used as isopiestic reference standards. The isopiestic results obtained have been combined with all other experimental thermodynamic quantities available in literature (osmotic coefficients, water activities, bromide mineral’s solubilities) to construct a chemical model that calculates solute and solvent activities and solid-liquid equilibria in the NaBr-H₂O, KBr-H₂O and Na-K-Br-H₂O systems from dilute to high solution concentration within the 0-300 °C temperature range. The Harvie and Weare [Harvie C., and Weare J. (1980) The prediction of mineral solubilities in naturalwaters: the Na-K-Mg-Ca-Cl-SO₄-H₂O system from zero to high concentration at 25 °C. Geochim. Cosmochim. Acta44, 981-997] solubility modeling approach, incorporating their implementation of the concentration-dependent specific interaction equations of Pitzer [Pitzer K. (1973) Thermodynamics of electrolytes. I. Theoretical basis and general equations. J. Phys. Chem.77, 268-277] is employed. The model for binary systems is validated by comparing activity coefficient predictions with those given in literature, and not used in the parameterization process. Limitations of the mixed solutions model due to data insufficiencies are discussed. This model expands the variable temperature sodium-potassium model of Greenberg and Moller [Greenberg J., and Moller N. (1989) The prediction of mineral solubilities in natural waters: a chemical equilibrium model for the Na-K-Ca-Cl-SO₄-H₂O system to high concentration from 0 to 250 °C. Geochim. Cosmochim. Acta53, 2503-2518] by evaluating Br⁻ pure electrolyte and mixing solution parameters and the chemical potentials of three bromide solid phases: NaBr-2H₂O (cr), NaBr (cr) and KBr (cr).     

Hemimorphite as a natural sink for arsenic in zinc deposits and related mine tailings: Evidence from single-crystal EPR spectroscopy and hydrothermal synthesis

Hemimorphite is a refractory mineral in surface environments and occurs commonly in supergene non-sulfide Zn deposits and Zn mine tailings. Single-crystal electron paramagnetic resonance (EPR) spectra of gamma-ray-irradiated hemimorphite from Mapimi (Durango, Mexico) reveal two arsenic-associated oxyradicals: [AsO₄]⁴⁻ and [AsO₄]²⁻. Inductively coupled plasma mass spectrometry analyses confirm this sample to contain 270 ppm As and that hemimorphite from other Zn deposits has appreciable amounts of arsenic as well. Spin Hamiltonian parameters, including matrices g, A (⁷⁵As) and P(⁷⁵As), show that the [AsO₄]⁴⁻ radical formed from electron trapping by a locally uncompensated [AsO₄]³⁻ ion substituting for [SiO₄]⁴⁻. Matrices g, A(⁷⁵As) and P(⁷⁵As) of the [AsO₄]²⁻ radical show it to have the unpaired spin on the bridging oxygen of an [AsO₄]³⁻ ion at a Si site and linked to a monovalent impurity ion. This structural model for the [AsO₄]²⁻ radical is further supported by observed ²⁹Si and ¹H superhyperfine structures arising from interactions with a single Si atom (A/g_eβ_e = ∼1 mT at B//c) and two equivalent H atoms (A/g_eβ_e = ∼0.3 mT at B∧b = 10°), respectively. Hydrothermal experiments at 200 °C and ∼9.5 MPa show that hemimorphite contains up to ∼2.5 wt% As₂O₅ and suggest that both the arsenate concentration and the pH value in the solution affect the As content in hemimorphite. These results demonstrate that hemimorphite is capable of sequestering arsenate in its crystal lattice, hence is a natural sink for attenuating As in supergene non-sulfide Zn deposits and Zn mine tailings. Moreover, results from hemimorphite potentially have more far-reaching implications for major silicates such as zeolites in the immobilization and removal of arsenic in surface environments.     

Solubility of NaCl in CO2 at high pressure and temperature: First experimental measurements

NaCl solubility in gaseous carbon dioxide has been measured in the pressure range from 30 to 70 MPa at 623 and 673 K. Our originally-designed high pressure apparatus allows in situ sampling of a portion of the fluid phase for chemical analysis. The results indicate that the solubility of NaCl increases with both temperature and pressure, and is about 4-5 orders of magnitude higher than saturated NaCl pressure values at the same temperature conditions (6.02 × 10⁻¹² at 623 K and 1.51 × 10⁻¹⁰ at 673 K). It is also 1-2 orders of magnitude greater than predictions according to the Equation of State of the ternary H₂O-CO₂-NaCl system by Duan, Moeller and Weare [Duan, Z., Moller, N., and Weare, J. H. (1995) Equation of state for the NaCl-H₂O-CO₂ system: prediction of phase equilibria and volumetric properties. Geochim. Cosmochim. Acta59, 2869] and has the opposite pressure dependence. The activity values of NaCl in the vapor phase, calculated from the experiments (with pure molten NaCl as a standard state in the vapor), have been fitted to the Darken Quadratic Formalism: , where, x_NaCl,v is mole the fraction of NaCl in the vapor phase, , , where P is the pressure in MPa and T the absolute temperature. Caution should be exerted while extrapolating this empirical equation far beyond the experimental P-T-compositional range.     

Microscale controls on the fate of contaminant uranium in the vadose zone, Hanford Site, Washington

An alkaline brine containing uranyl (UO₂²⁺) leaked to the thick unsaturated zone at the Hanford Site. We examined samples from this zone at microscopic scale to determine the mode of uranium occurrence—microprecipitates of uranyl (UO₂²⁺) silicate within lithic-clast microfractures—and constructed a conceptual model for its emplacement, which we tested using a model of reactive diffusion at that scale. The study was driven by the need to understand the heterogeneous distribution of uranium and the chemical processes that controlled it. X-ray and electron microprobe imaging showed that the uranium was associated with a minority of clasts, specifically granitic clasts occupying less than four percent of the sediment volume. XANES analysis at micron resolution showed the uranium to be hexavalent. The uranium was precipitated in microfractures as radiating clusters of uranyl silicates, and sorbed uranium was not observed on other surfaces. Compositional determinations of the 1-3 μm precipitates were difficult, but indicated a uranyl silicate. These observations suggested that uranyl was removed from pore waters by diffusion and precipitation in microfractures, where dissolved silica within the granite-equilibrated solution would cause supersaturation with respect to sodium boltwoodite. This hypothesis was tested using a reactive diffusion model operating at microscale. Conditions favoring precipitation were simulated to be transient, and driven by the compositional contrast between pore and fracture space. Pore-space conditions, including alkaline pH, were eventually imposed on the microfracture environment. However, conditions favoring precipitation were prolonged within the microfracture by reaction at the silicate mineral surface to buffer pH in a solubility limiting acidic state, and to replenish dissolved silica. During this time, uranyl was additionally removed to the fracture space by diffusion from pore space. Uranyl is effectively immobilized within the microfracture environment within the presently unsaturated Vadose Zone.     

Dissolution of uranophane: An AFM, XPS, SEM and ICP study

Dissolution experiments on single crystals of uranophane and uranophane-β, Ca(H₂O)₅[(UO₂)(SiO₃(OH)]₂, from the Shinkolobwe mine of the Democratic Republic of Congo, were done in an aqueous HCl solution of pH 3.5 for 3 h, in HCl solutions of pH 2 for 5, 10 and 30 min, and in Pb²⁺-, Ba-, Sr-, Ca- and Mg-HCl solutions of pH 2 for 30 min. The basal surfaces of the treated uranophane crystals were examined using atomic-force microscopy (AFM), X-ray photoelectron spectroscopy (XPS) and scanning electron microscopy (SEM). Solutions after dissolution experiments on single crystals and synthetic powders were analysed with inductively coupled plasma-optical emission spectroscopy (ICP-OES) and mass spectroscopy (ICP-MS). The morphology of the observed etch pits (measured by AFM) were compared to the morphology, predicted on the basis of the bond-valence deficiency of polyhedron chains along the edges of the basal surface. Etch pits form in HCl solutions of pH 2. Their decrease in depth with the duration of the dissolution experiment is explained with the stepwave dissolution model, which describes the lowering of the surrounding area of an etch pit with continuous waves of steps emanated from the etch pit into the rest of the crystal surface. Hillocks form in an HCl solution of pH 3.5, and the chemical composition of the surface (as indicated by XPS) shows that these hillocks are the result of the precipitation of a uranyl-hydroxy-hydrate phase. Well-orientated hillocks form on the surface of uranophane in a SrCl₂-HCl solution of pH 2. They are part of an aged silica coating of composition Si₂O₂(OH)₄(H₂O)_n. An amorphous layer forms on the surface of uranophane in a MgCl₂-HCl solution of pH 2, which has a composition and structure similar to silicic acid. Small crystallites of uranyl-hydroxy-hydrate phases form on the surface of uranophane after treatment in Pb(NO₃)₂-HCl and BaCl₂-HCl solutions of pH 2. Dissolution experiments on synthetic uranophane powders show that in the early stage of the experiments, the dissolution rate of uranophane increase in the sequence Pb(NO₃)₂-HCl < BaCl₂-HCl < CaCl₂-HCl < HCl < SrCl₂-HCl < MgCl₂-HCl, indicating that the dissolution of uranophane is more enhanced in solutions containing divalent cations of small ionic radii and high Lewis acidity (Mg, MgCl⁺).     

The activity coefficients of Fe(III) hydroxide complexes in NaCl and NaClO4 solutions

The osmotic coefficients of FeCl₃ at 25 °C from 0.15 to 1.7 m [Rumyantsev et al., Z. Phys. Chem., 218, 1089-1127, 2004] have been used to determine the Pitzer parameters (β⁽⁰⁾, β⁽¹⁾ and C^?) for FeCl₃. Since the differences in the Pitzer coefficients of rare earths in NaCl and NaClO₄ are small, the values of Fe(ClO₄)₃ have been estimated using the differences between La(ClO₄)₃ and LaCl₃. The Pitzer coefficients for FeCl₃ combined with enthalpy and heat capacity data for the rare earths can be used to estimate the activity coefficients of Fe³⁺ in NaCl over a wide range of temperatures (0 to 50 °C) and ionic strength (0 to 6 m).The activity coefficients of Fe³⁺ in NaCl and NaClO₄ solutions have been used to determine the activity coefficients of Fe(OH)²⁺ in these solutions from the measured first hydrolysis constants of Fe³⁺ [Byrne et al., Mar. Chem., 97, 34-48, 2005]. The activity coefficients of , Fe(OH)₃ and from 0 to 50 °C have also been determined from the solubility measurements of Fe(III) in NaCl solutions [Liu and Millero, Geochim. Cosmochim Acta, 63, 3487-3497, 1999]. These activity coefficients have been fitted to the Pitzer equations. These results can be used to estimate the speciation of Fe(III) with OH⁻ in natural waters with high concentrations of NaCl from 0 to 50 °C.     

Organic matter and metamorphic history of CO chondrites

The metamorphic grades of a series of eight CO3 chondrites (ALHA77307, Colony, Kainsaz, Felix, Lancé, Ornans, Warrenton and Isna) have been quantified. The method used was based on the structural grade of the organic matter trapped in the matrix, which is irreversibly transformed by thermal metamorphism. The maturation of the organic matter is independent with respect to the mineralogical context and aqueous alteration. This metamorphic tracer is thus valid whatever the chemical class of chondrites. Moreover, it is sensitive to the peak metamorphic temperature.The structural grade of the organic matter was used along with other metamorphic tracers such as petrography of opaque minerals, Fa and Fs silicate composition in type I chondrules, presolar grains and noble gas (P3 component) abundance. The deduced metamorphic hierarchy and the attributed petrographic types are the following: ALHA77307 (3.03) < Colony (3.1) < Kainsaz (3.6) < Felix (3.6 (1)) < Ornans (3.6 (2)) < Lancé (3.6 (3)) < Warrenton (3.7 (1)) < Isna (3.7 (2)).For most metamorphosed objects, the peak metamorphic temperature can be estimated using a geothermometer calibrated with terrestrial metasediments [Beyssac O., Goffe B., Chopin C., and Rouzaud J. N. (2002) Raman spectrum of carbonaceous material in metasediments: a new geothermometer. J. Metamorph. Geol., 20, 859-871]. A value of 330 °C was obtained for Allende (CV chondrite), Warrenton and Isna, consistent with temperatures estimated from Fe diffusion [Weinbruch S., Armstrong J., and Palme H. (1994). Constraints on the thermal history of the Allende parent body as derive from olivine-spinel thermometry and Fe/Mg interdiffusion in olivine. Geochim. Cosmochim. Acta58(2), 1019-1030.], from the Ni content in sulfide-metal assemblages [Zanda B., Bourot-Denise M., and Hewins R. (1995) Condensate sulfide and its metamorphic transformations in primitive chondrites. Meteorit. Planet. Sci.30, A605.] and from the d₀₀₂ interlayer spacing in poorly graphitized carbon [Rietmeijer, F., and MacKinnon, I. (1985) Poorly graphitized carbon as a new cosmothermometer for primitive extraterrestrial materials. Nature, 315, 733-736]. The trapped noble gas and C content appear to be sensitive but not precise metamorphic tracers, indicating that the “Ornans paradox” does not exist. Major problems with the current petrologic types derived from Induced ThermoLuminescence are pointed out.     

Solubility of plutonium(VI) carbonate in saline solutions

Among the plutonium oxidation states found to form in the environment, mobile plutonium(VI) can exist under oxidizing conditions and in waters with high chloride content due to radiolysis effects. We are investigating the solubility and speciation of plutonium(VI) carbonate under conditions relevant to natural waters and brines such as those found near some geologic radioactive waste repositories. The solid Pu(VI) phase PuO₂CO₃(s) was prepared and its solubility was measured in NaCl and NaClO₄ solutions in a CO₂ atmosphere as a function of pH and ionic strength (0.1-5.6 m). The concentration of soluble plutonium in solution was calculated from spectroscopic data and liquid scintillation counting. Spectroscopic measurements also revealed the plutonium oxidation state. The apparent solubility product of PuO₂CO₃(s) was determined at selected electrolyte concentrations to be, log K_s,0 = −13.95 ± 0.07 (0.1 m NaCl), log K_s,0 = −14.07 ± 0.13 (5.6 m NaCl), and log K_s,0 = −15.26 ± 0.11 (5.6 m NaClO₄). Specific ion interaction theory was used to calculate the solubility product at zero ionic strength, .     

Hydrochemical characteristics of shallow groundwater in aquifer containing uranyl phosphate minerals, in the Köprübaşı (Manisa) area, Turkey

The concentrations of uranium, iron and the major constituents were determined in groundwater samples from aquifer containing uranyl phosphate minerals (meta-autunite, meta-torbernite and torbernite) in the Köprüba?? area. Groundwater samples from wells located at shallow depths (0.5–6 m) show usually near neutral pH values (6.2–7.1) and oxidizing conditions (Eh = 119–275 mV). Electrical conductivity (EC) values of samples are between 87 and 329 μS/cm^?1. They are mostly characterized by mixed cationic Ca dominating bicarbonate types. The main hydrogeochemical process is weathering of the silicates in the shallow groundwater system. All groundwater in the study area are considered undersaturated with respect to torbernite and autunite. PHREEQC predicted UO₂(HPO₄) ₂ ^2? as the unique species. The excellent positive correlation coefficient (r = 0.99) between U and PO₄ indicates the dissolved uranium in groundwater would be associated with the dissolution of uranyl phosphate minerals. The groundwater show U content in the range 1.71–70.45 μg/l but they are mostly lower than US EPA (2003) maximum contaminant level of 30 μg/l. This low U concentrations in oxic groundwater samples is attributed to the low solubility of U(VI) phosphate minerals under near neutral pH and low bicarbonate conditions. Iron closely associated with studied sediments, were also detected in groundwater. The maximum concentration of Fe in groundwater samples was 2837 μg/l, while the drinking water guidelines of Turkish (TSE 1997) and US EPA (2003) were suggested 200 and 300 μg/l, respectively. Furthermore, iron and uranium showed a significant correlation to each other with a correlation coefficient (r) of 0.94. This high correlation is probably related to the iron-rich sediments which contain also significant amounts of uranium mineralization. In addition to pH and bicarbonate controlling dissolution of uranyl phosphates, association of uranyl phosphates with iron (hydr) oxides seems to play important role in the amount of dissolved U in shallow groundwater.     

Thermodynamic characterization of boltwoodite and uranophane: Enthalpy of formation and aqueous solubility

Boltwoodite and uranophane are uranyl silicates common in oxidized zones of uranium ore deposits. An understanding of processes that impact uranium transport in the environment, especially pertaining to the distribution of uranium between solid phases and aqueous solutions, ultimately requires determination of thermodynamic parameters for such crystalline materials. We measured formation enthalpies of synthetic boltwoodites, K(UO₂)(HSiO₄)·H₂O and Na(UO₂)(HSiO₄)·H₂O, and uranophane, Ca(UO₂)₂(HSiO₄)₂·5H₂O, by high temperature oxide melt solution calorimetry. We also studied the aqueous solubility of these phases from both saturated and undersaturated conditions at a variety of pH. The combined data permit the determination of standard enthalpies, entropies and Gibbs free energies of formation for each phase and analysis of its potential geological impact from a thermodynamic point of view.     

Molecular H2O in armenite, BaCa2Al6Si9O30·2H2O, and epididymite, Na2Be2Si6O15·H2O: Heat capacity, entropy and local-bonding behavior of confined H2O in microporous silicates

Armenite, ideal formula BaCa₂Al₆Si₉O₃₀·2H₂O, and its dehydrated analog BaCa₂Al₆Si₉O₃₀ and epididymite, ideal formula Na₂Be₂Si₆O₁₅·H₂O, and its dehydrated analog Na₂Be₂Si₆O₁₅ were studied by low-temperature relaxation calorimetry between 5 and 300 K to determine the heat capacity, C_p, behavior of their confined H₂O. Differential thermal analysis and thermogravimetry measurements, FTIR spectroscopy, electron microprobe analysis and powder Rietveld refinements were undertaken to characterize the phases and the local environment around the H₂O molecule.The determined structural formula for armenite is Ba_0.88(0.01)Ca_1.99(0.02)Na_0.04(0.01)Al_5.89(0.03)Si_9.12(0.02)O₃₀·2H₂O and for epididymite Na_1.88(0.03)K_0.05(0.004)Na_0.01(0.004)Be_2.02(0.008)Si_6.00(0.01)O₁₅·H₂O. The infrared (IR) spectra give information on the nature of the H₂O molecules in the natural phases via their H₂O stretching and bending vibrations, which in the case of epididymite only could be assigned. The powder X-ray diffraction data show that armenite and its dehydrated analog have similar structures, whereas in the case of epididymite there are structural differences between the natural and dehydrated phases. This is also reflected in the lattice IR mode behavior, as observed for the natural phases and the H₂O-free phases. The standard entropy at 298 K for armenite is S° = 795.7 ± 6.2 J/mol K and its dehydrated analog is S° = 737.0 ± 6.2 J/mol K. For epididymite S° = 425.7 ± 4.1 J/mol K was obtained and its dehydrated analog has S° = 372.5 ± 5.0 J/mol K. The heat capacity and entropy of dehydration at 298 K are Δ = 3.4 J/mol K and ΔS^rxn = 319.1 J/mol K and Δ = −14.3 J/mol K and ΔS^rxn = 135.7 J/mol K for armenite and epididymite, respectively. The H₂O molecules in both phases appear to be ordered. They are held in place via an ion-dipole interaction between the H₂O molecule and a Ca cation in the case of armenite and a Na cation in epididymite and through hydrogen-bonding between the H₂O molecule and oxygen atoms of the respective silicate frameworks. Of the three different H₂O phases ice, liquid water and steam, the C_p behavior of confined H₂O in both armenite and epididymite is most similar to that of ice, but there are differences between the two silicates and from the C_p behavior of ice. Hydrogen-bonding behavior and its relation to the entropy of confined H₂O at 298 K is analyzed for various microporous silicates.The entropy of confined H₂O at 298 K in various silicates increases approximately linearly with increasing average wavenumber of the OH-stretching vibrations. The interpretation is that decreased hydrogen-bonding strength between a H₂O molecule and the silicate framework, as well as weak ion-dipole interactions, results in increased entropy of H₂O. This results in increased amplitudes of external H₂O vibrations, especially translations of the molecule, and they contribute strongly to the entropy of confined H₂O at T < 298 K.     

XPS spectra of uranyl minerals and synthetic uranyl compounds. II: The O 1s spectrum

The O 1s spectrum is examined for 19 uranyl minerals of different composition and structure. Spectra from single crystals were measured with a Kratos Axis Ultra X-ray Photoelectron Spectrometer with a magnetic-confinement charge-compensation system. Well-resolved spectra with distinct maxima, shoulders and inflections points, in combination with reported and measured binding energies for specific O²⁻ species and structural data of the uranyl minerals are used to resolve the fine structure of the O 1s envelope. The resolution of the O 1s spectra includes, for the first time, different O²⁻ bands, which are assigned to O atoms linking uranyl with uranyl polyhedra (UOU) and O atoms of uranyl groups (OUO). The resolved bands in the O 1s spectrum occur at distinct ranges in binding energy: bands for (UOU) occur at 529.6-530.4 eV, bands for (OUO) occur at 530.6-531.4 eV, bands for O²⁻ in the equatorial plane of the uranyl polyhedra linking uranyl polyhedra with (TO_n) groups (T = Si, S, C, P, Se) (TO) occur at 530.9-532.2 eV; bands for (OH) groups in the equatorial plane of the uranyl polyhedra (OH) occur at 532.0-532.5 eV, bands of (H₂O) groups in the interstitial complex of the uranyl minerals (H₂O_interst) occur at 533.0-533.8 eV and bands of physisorbed (H₂O) groups on the surface of uranyl minerals (H₂O_adsorb) occur at 534.8-535.2 eV. Treatment of uranyl minerals with acidic solutions results in a decrease in UOU and an increase in OH. Differences in the ratio of OH : OUO between the surface and bulk structure is larger for uranyl minerals with a high number of UOU and TO species in the bulk structure which is explained by protonation of underbonded UO, UOU and TO terminations on the surface. The difference in the ratio of H₂O_interst : OUO between the bulk and surface structures is larger for uranyl minerals with higher percentages of H₂O_interst as well as, with a higher number of interstitial H₂O groups that are not bonded to interstitial cations, resulting in easier dehydration of interstitial H₂O groups in uranyl minerals during exposure to a vacuum.     

Zinc complexation in aqueous sulfide solutions: Determination of the stoichiometry and stability of complexes via ZnS(cr) solubility measurements at 100 °C and 150 bars

The solubility of ZnS_(cr) was measured at 100 °C, 150 bars in sulfide solutions as a function of sulfur concentration (m(S_total) = 0.02-0.15) and acidity (pH_t = 2-11). The experiments were conducted using a Ti flow-through hydrothermal reactor enabling the sampling of large volumes of solutions at experimental conditions, with the subsequent concentration and determination of trace quantities of Zn. Prior to the experiments, a long-term in situ conditioning of the solid phase was performed in order to attain the reproducible Zn concentrations (i.e. solubilities). The ZnS_(cr) solubility product was monitored in the course of the experiment. The following species were found to account for Zn speciation in solution: Zn²⁺ (pH_t < 3), (pH_t 3-4.5), (pH_t 5-8), and ZnS(HS)⁻ (pH_t > 8) (pH_t predominance regions are given for m(S_total) = 0.1). Solubility data collected in this study at pH_t > 3 were combined with the ZnS_(cr) solubility product determined at lower pH to yield the following equilibrium constants (t = 100 °C, P = 150 bars):

     

A thermodynamic model for calculating methane solubility, density and gas phase composition of methane-bearing aqueous fluids from 273 to 523 K and from 1 to 2000 bar

A thermodynamic model is presented to calculate methane solubility, liquid phase density and gas phase composition of the H₂O-CH₄ and H₂O-CH₄-NaCl systems from 273 to 523 K (possibly up to 573 K), from 1 to 2000 bar and from 0 to 6 mol kg⁻¹ of NaCl with experimental accuracy. By a more strict theoretical approach and using updated experimental data, this model made substantial improvements over previous models: (1) the accuracy of methane solubility in pure water in the temperature range between 273 and 283 K is increased from about 10% to about 5%, but confirms the accuracy of the Duan model [Duan Z., Moller N., Weare J.H., 1992a. Prediction of methane solubilities in natural waters to high ionic strength from 0 to 250 °C and from 0 to 1600 bar. Geochim. Cosmochim. Acta56, 1451-1460] above 283 K up to 2000 bar; (2) the accuracy of methane solubility in the NaCl aqueous solutions is increased from >12% to about 6% on average from 273 K and 1 bar to 523 K and 2000 bar; (3) this model is able to calculate water content in the gas phase and liquid phase density, which cannot be calculated by previous models; and (4) it covers a wider range of temperature and pressure space. With a simple approach, this model is extended to predict CH₄ solubility in other aqueous salt solutions containing Na⁺, K⁺, Mg²⁺, Ca²⁺, Cl⁻ and , such as seawater and geothermal brines, with excellent accuracy. This model is also able to calculate homogenization pressure of fluid inclusions (CH₄-H₂O-NaCl) and CH₄ solubility in water at gas-liquid-hydrate phase equilibrium. A computer code is developed for this model and can be downloaded from the website: www.geochem-model.org/programs.htm.     

Solubilities of corundum, wollastonite and quartz in H2O-NaCl solutions at 800 °C and 10 kbar: Interaction of simple minerals with brines at high pressure and temperature

Solubilities of corundum (Al₂O₃) and wollastonite (CaSiO₃) were measured in H₂O-NaCl solutions at 800 °C and 10 kbar and NaCl concentrations up to halite saturation by weight-loss methods. Additional data on quartz solubility at a single NaCl concentration were obtained as a supplement to previous work. Single crystals of synthetic corundum, natural wollastonite or natural quartz were equilibrated with H₂O and NaCl at pressure (P) and temperature (T) in a piston-cylinder apparatus with NaCl pressure medium and graphite heater sleeves. The three minerals show fundamentally different dissolution behavior. Corundum solubility undergoes large enhancement with NaCl concentration, rising rapidly from Al₂O₃ molality (m_Al2O3) of 0.0013(1) (1σ error) in pure H₂O and then leveling off to a maximum of ∼0.015 at halite saturation (X_NaCl ≈ 0.58, where X is mole fraction). Solubility enhancement relative to that in pure H₂O, , passes through a maximum at X_NaCl ≈ 0.15 and then declines towards halite saturation. Quenched fluids have neutral pH at 25 °C. Wollastonite has low solubility in pure H₂O at this P and T(m_CaSiO3=0.0167(6)). It undergoes great enhancement, with a maximum solubility relative to that in H₂O at X_NaCl ≈ 0.33, and solubility >0.5 molal at halite saturation. Solute silica is 2.5 times higher than at quartz saturation in the system H₂O-NaCl-SiO₂, and quenched fluids are very basic (pH 11). Quartz shows monotonically decreasing solubility from m_SiO2=1.248 in pure H₂O to 0.202 at halite saturation. Quenched fluids are pH neutral. A simple ideal-mixing model for quartz-saturated solutions that requires as input only the solubility and speciation of silica in pure H₂O reproduces the data and indicates that hydrogen bonding of molecular H₂O to dissolved silica species is thermodynamically negligible. The maxima in for corundum and wollastonite indicate that the solute products include hydrates and Na⁺ and/or Cl⁻ species produced by molar ratios of reactant H₂O to NaCl of 6:1 and 2:1, respectively. Our results imply that quite simple mechanisms may exist in the dissolution of common rock-forming minerals in saline fluids at high P and T and allow assessment of the interaction of simple, congruently soluble rock-forming minerals with brines associated with deep-crustal metamorphism.     

Coupled phase and aqueous species equilibrium of the H2O-CO2-NaCl-CaCO3 system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite

A model is developed for the calculation of coupled phase and aqueous species equilibrium in the H₂O-CO₂-NaCl-CaCO₃ system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite. The vapor-liquid-solid (calcite, halite) equilibrium together with the chemical equilibrium of H⁺, Na⁺, Ca²⁺, , Ca(OH)⁺, OH⁻, Cl⁻, , , CO_2(aq) and CaCO_3(aq) in the aqueous liquid phase as a function of temperature, pressure, NaCl concentrations, CO_2(aq) concentrations can be calculated, with accuracy close to those of experiments in the stated T-P-m range, hence calcite solubility, CO₂ gas solubility, alkalinity and pH values can be accurately calculated. The merit and advantage of this model is its predictability, the model was generally not constructed by fitting experimental data.One of the focuses of this study is to predict calcite solubility, with accuracy consistent with the works in previous experimental studies. The resulted model reproduces the following: (1) as temperature increases, the calcite solubility decreases. For example, when temperature increases from 273 to 373 K, calcite solubility decreases by about 50%; (2) with the increase of pressure, calcite solubility increases. For example, at 373 K changing pressure from 10 to 500 bar may increase calcite solubility by as much as 30%; (3) dissolved CO₂ can increase calcite solubility substantially; (4) increasing concentration of NaCl up to 2 m will increase calcite solubility, but further increasing NaCl solubility beyond 2 m will decrease its solubility.The functionality of pH value, alkalinity, CO₂ gas solubility, and the concentrations of many aqueous species with temperature, pressure and NaCl_(aq) concentrations can be found from the application of this model. Online calculation is made available on www.geochem-model.org/models/h2o_co2_nacl_caco3/calc.php.