Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrate”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Optimized thermodynamic properties of REE aqueous species (REE 3+ and REEOH 2+ ) and experimental database for modeling the solubility of REE phosphate minerals (monazite, xenotime, and rhabdophane) from 25 to 300 °C

Rare earth elements (REE) are critical elements found in monazite, xenotime, and hydrated REE phosphates which typically form in hydrothermal mineral deposits. Accurate predictions of the solubility of these REE phosphates and the speciation of REE in aqueous fluids are both key to understanding the controls on the transport, fractionation, and deposition of REE in natural systems. Previous monazite and xenotime solubility experiments indicate the presence of large discrepancies between experimentally derived solubility constants versus calculated solubilities by combining different data sources for the thermodynamic properties of minerals and aqueous species at hydrothermal conditions. In this study, these discrepancies were resolved by using the program GEMSFITS to optimize the standard partial molal Gibbs energy of formation (Δ f G° 298 ) of REE aqueous species (REE 3+ and REE hydroxyl complexes) at 298.15 K and 1 bar while keeping the thermodynamic properties fixed for the REE phosphates. A comprehensive experimental database was compiled using solubility data available between 25 and 300 °C. The latter permits conducting thermodynamic parameter optimization of Δ f G° 298 for REE aqueous species. Optimal matching of the rhabdophane solubility data between 25 and 100 °C requires modifying the Δ f G° 298 values of REE 3+ by 1–6 kJ/mol, whereas matching of the monazite solubility data between 100 and 300 °C requires modifying the Δ f G° 298 values of both REE 3+ and REEOH 2+ by ~15–31 kJ/mol and ~2–10 kJ/mol, respectively. For xenotime, adjustments of Δ f G° 298 values by 1–26 kJ/mol are only necessary for the REE 3+ species. The optimizations indicate that the solubility of monazite in acidic solutions is controlled by the light (L)REE 3+ species at <150 °C and the LREEOH 2+ species at >150 °C, whereas the solubility of xenotime is controlled by the heavy (H)REE 3+ species between 25 and 300 °C. Based on the optimization results, we conclude that the revised Helgeson-Kirkham-Flowers equation of state does not reliably predict the thermodynamic properties of REE 3+ , REEOH 2+ , and likely other REE hydroxyl species at hydrothermal conditions. We therefore provide an experimental database (ThermoExp_REE) as a basic framework for future updates, extensions with other ligands, and optimizations as new experimental REE data become available. As a result, the optimized thermodynamic properties of aqueous species and minerals are available open access to accurately predict the solubility of REE phosphates in fluid-rock systems.

58 GEOSCIENCES↗

Quantification of the REE 3+ aqua ions and chloride species in aqueous fluids by in situ Raman spectroscopy using perturbations of the water band

Acidic NaCl-rich aqueous fluids play a crucial role in forming hydrothermal rare earth elements (REE) mineral deposits. Aqueous REE mobility is mostly controlled by the stabilities of REE 3+ and REE chloride species. Our current knowledge of REE speciation is based on solubility data, thermodynamic models and in situ spectroscopic measurements, sometimes coupled with molecular simulations. Here, in this study, we investigate Nd and Yb speciation in pH2 Cl-bearing solutions at 25 °C and 0.1 MPa with variable Cl/REE ratios using Raman Spectroscopy in solutions with 0.1 to 0.6 mol/kg NdCl 3 or YbCl 3 and 0.2 to 3.2 mol/kg NaCl. Due to the challenges in resolving the REE-Cl band, we developed a new method using the water vibrational mode and multivariate curve resolution (MCR) analysis. The Raman spectra for the vibrational band of water (2700 to 3900 cm –1 ) were collected at 25 °C and fitted by three Gaussian sub-peaks, then quantified using MCR analysis to de-convolute the water band into bulk H 2 O and the perturbations caused by of Cl – , REE 3+ , and REE chloride species. REE speciation based on the perturbations of the water band indicates that REE 3+ aqua ions dominate acidic solutions at 25 °C, but up to ~20 mol% YbCl 2+ forms at high YbCl 3 concentrations. The new method is promising for quantifying in situ speciation of the REE 3+ aqua ions and REE chloride species in aqueous fluids while providing information on the hydration of ions. This method improves our molecular level understanding of REE aqueous species stability and their role in REE mobilization during fluid-rock interaction.

58 GEOSCIENCES↗

Effects of sodium and calcium chloride ionic stresses on model yeast membranes revealed by molecular dynamics simulation

As efforts to move a renewable economy grow, it will be necessary to make use of microbial conversion strategies for the production of novel materials or the upgrading of waste to high-value products. One critical technical challenge currently limiting waste upgrading remains the difficulty in obtaining single-pot conversion techniques where physical, chemical, and biological conversion are performed in a single step. To overcome this challenge, a detailed understanding of how different stresses impact microbial membrane stability will be necessary. In this work, using all-atom molecular dynamics simulations, we examine the impacts of moderate concentrations of NaCl and CaCl 2 on a model yeast plasma membrane. Weak, though statistically significant, changes in membrane morphology and dynamics functions are observed that are consistent with swelling and stiffening. Additionally, an examination of the ion-lipid contacts and the behavior of water at the water-membrane interface suggests that the impacts of these common salts may, in part, be mediated through changes to water-membrane hydrogen-bonding and hydration water dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Colloidal phenomena reflect the interplay between interfacial solution structure, interparticle forces, and dynamical response

The physical and chemical properties of colloidal systems are central to both natural and synthetic processes that enable a wide range of energy and environmental applications. These properties are an emergent outcome of the structure that colloidal systems adopt in response to the interplay of electrical, van der Waals, hydration, and steric forces coupled to the inherent fluctuations of the system. Here, in this study, we review the nature of those forces and show how variations in their relative importance lead to distinct free energy landscapes and consequent dynamic behaviors, illustrated for colloidal dispersions, nanocrystal superlattices, and mesocrystlline solids formed through oriented attachment. We end by discussing the challenges to future progress arising from the need to seamless couple from molecular to particle to ensemble scales.

42 ENGINEERING↗

Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale surface features of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. Additionally, a simple analytical model for considering the average effect of such surface features producing nanoscale roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to unambiguously determine the zeta potential in the presence of nanoscale surface features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plasma decay in hydrocarbons and hydrocarbon- and H 2 O-containing mixtures excited by high-voltage nanosecond discharge at elevated gas temperatures

Plasma decay after a high-voltage nanosecond discharge was experimentally and numerically studied in pure hydrocarbons (C 2 H 6 and C 3 H 8 ), and H 2 O:N 2 and C 3 H 8 :O 2 mixtures for pressures in the range 2–4 Torr and gas temperatures from 300 to 600 K. In a stoichiometric C 3 H 8 :O 2 mixture, plasma decay was also studied in a repetitively pulsed discharge for varing numbers of discharge pulses (varying degrees of fuel oxidation). The rate of plasma decay was determined from the temporal evolution of electron density measured using the microwave interferometer. It was observed that gas heating to 600 K leads to a decrease in the rate of plasma decay in all cases. The effect of heating on plasma decay was most profound in the H 2 O:N 2 mixture (after a single discharge pulse) and in the C 3 H 8 :O 2 mixture for high degrees of fuel oxidation. A kinetic scheme was developed to numerically simulate the plasma decay in hydrocarbons and combustible mixtures. Numerical analysis showed that, under the conditions studied, plasma decay was controlled by dissociative electron recombination with simple molecular and cluster ions. Gas heating led to a decrease in the rate of the electron-ion recombination and the rate of conversion of molecular ions to cluster ions. As a result, the gas temperature increase caused a decrease in the fraction of cluster ions for which the recombination coefficients are an order of magnitude higher than the recombination coefficients for molecular ions. The influence of gas heating on the decrease of the amount of cluster ions was more important when the ion composition was dominated by hydrated H 3 O + (H 2 O) k ions. The rates of the formation of these ions are extremely sensitive to any variations in gas temperature. A result, in agreement with our observations, gas heating led to an anomalous decrease in the rates of plasma decay in the H 2 O:N 2 mixture, as well as in the C 3 H 8 :O 2 mixture when H 2 O molecules were produced due to fuel oxidation.

20 FOSSIL-FUELED POWER PLANTS↗

Microstructural characterization and assessment of mechanical properties of concrete based on combined elemental analysis techniques and Fast-Fourier transform-based simulations

This paper presents a pixel-based modeling approach of concrete which combines an experimental characterization of concrete and the Fast-Fourier transform simulations. High-resolution phase maps created from experimental characterization by micro X-ray fluorescence, energy dispersive X-ray analysis, and X-ray diffraction contain 9 different phases, including 22.17 vol.% of hydrated cement paste, 54.21 vol.% of minerals, and 23.62 vol.% of interfaces. These phases provide the input for determining the effective elastic properties, coupled with Fast-Fourier transform-based simulation. The simulation results show that the effective range of Young’s modulus of concrete is comparable with the range of experimental values ~ 37 ± 4 GPa with the assumption of realistic properties of interfaces.

36 MATERIALS SCIENCE↗

Experimental study of mechanistic factors influencing solvent-driven fractional crystallization of calcium sulfate

To advance dimethyl ether-driven fractional crystallization (DME-FC), a more sustainable method of water treatment and mineral recovery, a range of chemical equilibria were measured. These include varying concentrations of miscible organic solvents (MOS) used to experimentally measure the solvent-induced solid-liquid equilibrium (SLE) of calcium sulfate (CaSO 4 ) in water. Seven MOS, including dimethyl ether (DME), acetonitrile (MeCN), 1,4-dioxane, tetrahydrofuran (THF), acetone, ethanol, and diethylamine, were screened to establish trends associated with molecular volume, functional groups, and physical properties. The effect of MOS on CaSO 4 removal differed at concentrations <0.15 mol fraction MOS; MOS with greater molecular volume (THF, 1,4-dioxane, and diethylamine) induced greater CaSO 4 precipitation on a per mole basis. The solvent-induced SLE for all MOS converged between 0.15 and 0.2 mol fraction MOS, reaching a CaSO 4 concentration consistent with a water to MOS hydration ratio of 5:1 to 6:1, which may correspond to the solvent generating a solution-based pseudo-clathrate structure with continuity within the solution. Finally, solution pseudo-clathrate structures provide a mechanistic basis for DME-FC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stochastic generation of electrolyzer anode catalyst layers

Here, we introduce a stochastic methodology to reproduce the complex pore structure observed in commercial iridium catalyst layers. This method preserves the α pore (pores smaller than 250 nm) and β pore (pores greater than or equal to 250 nm) regions of the catalyst layer. The morphology of the generated materials was validated by comparing the pore size distributions of generated materials against those obtained from commercial materials imaged using x-ray nano computed tomography. We further demonstrate that the pore size distributions of the generated materials are statistically indistinguishable from the imaged catalyst layers, indicating that the stochastic methodology is capable of accurately reproducing catalyst layer morphology. Pore network modelling was conducted on the generated catalyst materials to simulate single-phase permeability, electrical conductivity, and ionic conductivity, and these properties were found to be within experimentally measured ranges for electrolyzer catalyst layers. Additionally, simulations were performed on the generated materials with varying ionomer and iridium catalyst loadings. As the ionomer loading is added, proton conductivity increases exponentially, which demonstrates the importance of optimizing ionomer loading, considering that these effects will be exacerbated in the hydration and temperature conditions of operating electrolyzers. The stochastic material generation method presented in this work is a powerful tool for the development of novel low loading catalyst layers, where the effect of various structural parameters on electrolyzer performance characteristics can be explored.

36 MATERIALS SCIENCE↗

Tuning and high throughput experimental screening of eutectic electrolytes with co-solvents for redox flow batteries

Eutectic solvents, with high salt concentrations and suppressed volatility, are promising alternatives to aqueous and volatile organic electrolytes for grid-scale energy storage devices such as redox flow batteries (RFBs). However, it is not known a priori which hydrogen bond acceptors and donors would form a eutectic solvent and at which compositions of these mixtures would yield desirable properties for RFBs. Specifically, the known deep eutectic solvents generally have high viscosities and low ionic conductivities, especially in the landscape of RFB electrolytes. We report a high throughput experimental (HTE) study was carried out to assess the properties of eutectic solvents based on choline chloride as the hydrogen bond acceptor and ethylene glycol and aniline as the hydrogen bond donors. In addition, water, acetonitrile, and dimethyl sulfoxide as co-solvents were examined in terms of tuning the viscosity, conductivity, and solubility of two redox-active species: methyl viologen dichloride hydrate (MVCl 2 ·H 2 O, or MV for simplicity) and 2,1,3-benzothiadiazole (BTZ) that are relevant to RFBs. The HTE framework enabled the development of a rapid and comprehensive understanding of the design of new eutectic electrolytes incorporating co-solvents to improve the physicochemical and electrochemical properties for RFBs.

25 ENERGY STORAGE↗

Effect of cell compression on the performance and the structure of proton exchange membrane water electrolyzer (PEMWE) assembly

Here, in the field of water electrolysis, the proton exchange membrane water electrolyzer (PEMWE) is currently the most advanced technology for producing hydrogen without emitting CO 2 . Although PEMWE plants are already in operation, further research is needed to improve cell efficiency and reduce the use of rare materials, such as iridium oxide catalysts for the oxygen evolution reaction (OER). One of the main causes of performance loss in PEMWE is the relatively low electric conductivity of the porous transport layer (PTL) and of the anode catalyst layer, which results in ohmic losses and low catalyst utilization during high current density operation. The objective of this study is to investigate how optimization of the PTL and electrode interface can increase the cell performance. To this end, we tested different cell assemblies using fibrous and sintered PTLs, decreasing membrane thickness, reducing iridium loading, and inserting a microporous layer to increase contact surface area. Electrochemical characterization of each cell configuration was systematically performed at various compression levels as the pressure is a crucial parameter influencing the electrode/PTL contact area. In parallel, X-ray microcomputed tomography (micro-CT) was employed to investigate the effects of cell hydration and compression on the structure of PEMWE components. This study combining electrochemistry and micro-CT imaging presents how optimizing the electrode/PTL contact surface area, minimizes ohmic losses, and enables PEMWE operation with low iridium loading at high current densities.

Catalyst - PTL interface↗

Influence of alkali-silica reaction on the shear capacity of reinforced concrete beams with minimum transverse reinforcement

A growing number of critical concrete infrastructure are affected by alkali-silica reaction (ASR) damage such as the Seabrook Nuclear Power Plant in New Hampshire, Parker arch-gravity dam in Arizona, and several highway bridges in California and Texas. ASR causes expansion and cracking and degrades the concrete mechanical properties. Despite a wealth of material level studies there is still limited large-scale experimental data regarding the effects of ASR on reinforced concrete (RC) members. Due to the brittle nature of the shear failure in RC structures, this study focuses on the shear response of full-scale ASR damaged RC beams with minimum shear reinforcement. Six RC beams were built with different levels of ASR susceptibility and conditioned in different environments during which continuous expansion monitoring was performed. The beams all contained reactive fine aggregate (sand) and two of them had additional alkali, 1.25% by weight of cement, to accelerate ASR. The highest expansion rate of the beams happened during the first 150 days and the expansions stayed constant after 240 days and 330 days for beams conditioned in outdoor and laboratory conditions, respectively, until the last measurement at 575 days. The highest expansion, 0.4%, was seen in the beams with additional alkali and conditioned outside with regular water spray. Out of the six beams, three were selected at different levels of ASR damage and two shear tests were performed on the minimally reinforced spans close to the ends of each beam. Results indicated that beams gain shear strength from ongoing cement hydration in the presence of moisture for ASR expansions less than 0.2%. Compared to one of the beams with 0.2% ASR expansion, one of the other samples with 0.4% expansion lost 6% of its shear strength, 25% of its shear stiffness, and showed about two times larger shear cracks and shear deformations at peak load. Finally, the shear reinforcement yielded at 20% less load in the beam with 0.4% expansion compared to the beam with 0.2% expansion.

42 ENGINEERING↗

Synchrotron-based X-ray fluorescence microscopy mapping the ionome of a toxic freshwater cyanobacterium

Harmful algal blooms (HABs) pose a major environmental concern across the globe. In abundance, cyanobacteria, or so-called green-blue algae can produce extremely dangerous cyanotoxins that harm humans and animals. This study focused on the mapping and distribution of intracellular macro-and micronutrients of the widespread freshwater cyanobacteria Microcystis aeruginosa (M. aeruginosa). Towards a better understanding of trace metal uptake and homeostasis throughout the cell cycle, we quantitatively mapped the spatial distribution of the elements P, K, Fe, Ca, Zn, Mn, and Cu across the ultrastructure of frozen-hydrated single cells using state-of-theart X-ray nanofluorescence imaging at the Advanced Photon Source (APS) at Argonne National Laboratory. Further, bulk cellular nutrient and trace metal content correlated well with the total intracellular elemental content in individual cells obtained by quantitative synchrotron X-ray fluorescence measurements. Multi-dimensional mappings showed P and K atoms colocalized as discrete semicircular hotspots that were analyzed with respect to their stoichiometry. Elevated Cu and Ca concentrations were detected along division plane of cells. P and K were found to have similar spatial elemental distribution with about 65% and 69% of the total cellular P and K, respectively, located at the hotspots. The P and K colocalization were refined further using nanotomography, showing a K envelope surrounding the P core. Inorganic P and organic P compounds were specified using solution-state 31 P nuclear magnetic resonance (NMR) spectroscopy from M. aeruginosa. Of the total extracted P determined by 31 P NMR spectroscopy, 47% were found to be nucleotides while only 11% were polyphosphates. Multimodal X-ray imaging provides a better understanding of intracellular biochemical processes in cyanobacteria, helping us monitor and combat an emerging environmental threat.

31P NMR spectroscopy↗

A comprehensive review of diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) techniques in protonic ceramic cells (PCCs): Current status and future perspective

Protonic ceramic cells (PCCs) have emerged as a promising technology for power generation, energy storage, and value-added chemical synthesis, offering benefits such as fuel flexibility, low emissions, and efficient operation at intermediate temperatures (300–600 ​°C). Recently, significant breakthroughs in materials and manufacturing methods have markedly enhanced the performance of PCCs. However, establishing a fundamental understanding of their electrocatalytic reactions has gained less attention. As a fast and cost-effective method for physicochemical fingerprinting, diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) has proven to be a surface-sensitive analytical tool for structural and functional studies. This review critically examines the most up-to-date applications of DRIFTS for characterizing key components of PCCs, including oxygen electrodes, protonic electrolytes, and hydrogen electrodes for different applications, with a focus on revealing hydration properties and catalytic reactions, and guiding rational material design. The challenges for advancing DRIFTS, including quantitative capabilities and operando applications for PCC investigations, are highlighted and strategies to tackle these challenges are discussed. Ultimately, this review underscores the critical role of DRIFTS in accelerating the development of high-performance and durable PCCs for next-generation energy solutions, offering methodologies and insights broadly applicable to a wide range of electrochemical energy conversion and storage technologies.

Diffuse Reflectance Infrared Fourier Transform Spe↗

Experimental apparatus and methodology to test and quantify thermal performance of micro and macro-encapsulated phase change materials in building envelope applications

Thermal energy storage (TES) is used as a viable technology to shift peak electricity demand caused by the space cooling requirements in buildings. Passive TES is implemented in building envelope via micro and macroencapsulation methods. This study describes the use of a state-of-the-art laboratory to test different PCM inclusions in simplified building walls. A microencapsulated PCM and two macroencapsulated PCMs are tested in a controlled environment to gather data for validation purposes of PCM modelling algorithms in building energy modelling programs. Data indicates that the chamber environment and the heating and cooling system can conduct full-cycle tests of wall panels with PCM inclusions. This study also generates data from parallel tests on 4 wall panels which can be used in building energy modelling programs to validate the PCM modelling algorithms. The cyclic tests also capture the thermal effects of PCMs and complex PCM behaviors like sub-cooling in PCM hydrate-salts.

25 ENERGY STORAGE↗

Thermodynamic modeling of aqueous lithium salt solutions with association electrolyte nonrandom two-liquid activity coefficient model

The high charge density of lithium ion and the resulting strong association phenomena make thermodynamic modeling of aqueous lithium electrolyte solutions extremely challenging. In this study, the association electrolyte nonrandom two-liquid activity coefficient model of Lin et al. (AIChE J. 2022, 68(2), e17422) is utilized to correlate and predict thermodynamic properties and solubility behavior of aqueous single electrolyte solutions of LiCl, LiBr, LiI, and LiNO 3 , and their mixed electrolyte solutions. Capturing self-association of water, cross-association of ion and water for hydration, and cross-association of cation and anion for ion-pairing, the association model accurately represents the literature experimental data up to saturation concentrations and at the temperature ranging from 263 K to 523 K. Here, this study further investigated the effect of anions of the lithium salts, and re-confirmed that the order of solution non-ideality as LiI > LiBr > LiCl > LiNO 3 because the anions with stronger association strengths are more likely to form ion pairs and thus lower the mean ionic activity coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrogen adsorption and transport in clay-rich geomaterials: Implications for large-volume underground hydrogen storage

Depleted oil and gas reservoirs, characterized by impermeable clay-rich caprocks, are promising sites for large-scale underground hydrogen storage (UHS), which is a key strategy to support hydrogen-based energy systems. However, experimental data on hydrogen storage in clay-rich geomaterials remain scarce. In this work, we experimentally investigated hydrogen adsorption and migration in clay-rich geomaterials in the presence of nitrogen and water under controlled temperatures. Experimental observations showed that hydrogen was adsorbed in dry illite. A dual-porosity transport model was developed to interpret hydrogen transport between large-pore and small-pore domains in illite. The large-pore domain is the space between clay particles (i.e., inter-particle space), whereas the small-pore domain is the nanoscale pore space between clay mineral layers (i.e., inter-layer or intra-particle space). In contrast, nitrogen showed no evidence of adsorption in dry illite because it cannot move into the small-pore domain due to the relatively large kinetic diameter, referred to as the molecular sieving effect. Here, we found that 0.7–1.3 nm is the length scale regulating this molecular sieving effect, matching the interlayer spacing in illite, suggesting that nitrogen is a promising cushion gas in UHS, which aims to maintain adequate pressure in the reservoir for economic operations. In wetted illite, hydrogen was not adsorbed into the interlayer space due to the occupation of adsorption sites by interlayer water, which highlights the critical role of the clay hydration state in controlling hydrogen-clay interactions. Additionally, hydrogen adsorption experiments on crushed shale indicated that the shale surface possessed adsorption sites more favorable for hydrogen than for nitrogen. Through these experiments, we provide new insights into hydrogen storage mechanisms in clay-rich geomaterials and offer valuable laboratory data for evaluating the performance of large-scale UHS systems.

Adsorption↗

An Experimental Study of the Solubility of Rare Earth Chloride Salts (La, Nd, Er) in HCl Bearing Water Vapor from 350 – 425 °C

In this work, the solubilities of the rare earth chlorides REECl 3 , where REE = (La, Nd, Er), were measured in HCl bearing water vapor from 350 – 425°C with water partial pressures ranging from 8 – 170 bar. Solubility data were fit to the Pitzer-Pabalan quasi-chemical model in order to extract thermodynamic parameters for the formation of the gaseous REE-chloride-water clusters REECl 3 (H 2 O) n . The data show that the solubility of the REE chlorides are orders of magnitude higher than salts such as NaCl or CuCl at low water fugacities, despite their sublimation energies being substantially higher. This enhanced solubility is likely due to the high enthalpy associated with binding a single water molecule to form the species REECl 3 (H 2 O), with derived enthalpies ranging from –378 to –465 kJ/mol. Addition of further water molecules to form higher order clusters (n > 1) involves enthalpy changes of ~ -20 kJ/mol, and are in effect thermodynamically suppressed over the temperature range 350 – 425°C. Despite the enhanced solubility of small REECl 3 water clusters, simulations of boiling processes demonstrate that the REE show highly conservative behavior, partitioning strongly into the dense aqueous phase. Not surprisingly, the presence of phosphates in the system makes this effect even more pronounced, completely immobilizing the REE. This would reduce transport in both the vapor and aqueous phase to negligible levels. However, we suggest that vapor phase transport of the REE may play a significant role in systems having a relatively low partial pressure of water (below the saturation point), where the relatively high stability of the first hydrated REE chloride clusters (REECl 3 (H 2 O) and REECl 3 (H 2 O) 2 ) will give a preference for gas transport of the REE relative to other elements. This can likely happen in systems involving a gas/melt exchange in fumarolic exhalations, where water vapor discharges at relatively low (close to atmospheric) pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗