Search NASA⌕ Search

SEARCH · Search NASA

Results for “Aqueous speciation modeling”

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.

88 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Speciation and Organic Phase Structure in Nitric Acid Extraction with Trioctylamine

Understanding chemical speciation and intermolecular interactions in multicomponent liquids is essential to understanding their phase and chemical equilibria, which underpin chemical separation processes, including solvent extraction. Here we report on the extraction of nitric acid from its aqueous solutions into organic solutions of trioctylamine (TOA) in toluene, investigated with spectroscopic, X-ray scattering, and computational tools to understand molecular speciation in the organic phase and its relationship with the nanoscale structure of the organic phase. Trends in acid and water extraction clearly show two and three regimes, respectively, indicating different stoichiometric relationships, but speciation of HNO 3 , water, and amine in these regimes is not apparent. 1 H NMR of the organic phase shows that there are at least two distinct acidic protons in the organic phase while ATR-FTIR results show that the organic phase with excess acid extraction is a mixture of trioctylammonium-nitrate ion pairs (TOAHNO 3 ), and undissociated HNO 3 molecules. Comparison with DFT-computed IR spectra show that the chain-like configurations of TOAHNO 3 HNO 3 H 2 O are favored over TOAHNO 3 H 2 OHNO 3 , i.e., direct interaction between the nitrate and HNO 3 molecules is more favored compared to a water-mediated interaction. SAXS of the organic phases were modeled as sums of Ornstein-Zernike (O-Z) scattering and a prepeak feature in the higher Q region that corresponds to extractant packing. The extraction of undissociated HNO 3 by the ion pairs leads to an increased X-ray scattering contrast in the organic phase without any significant change in the correlation length. These results show that the organic phase nanostructure is more sensitive to the concentration of TOAHNO 3 and is relatively unaffected by excess acid extraction. Finally, these findings will enable a molecular understanding of the mechanisms behind metal extraction from acidic media with basic extractants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Quantitative Prediction of Uranium Speciation and Amidoxime Binding in Seawater from Advanced Simulation Techniques

The goal of this project was to quantitatively predict stability of uranium species and thermochemistry of uranium binding with amidoxime ligands in realistic seawater conditions. The specific aims of the project were: (i) determine the stable species of uranium and its solvation structure in realistic seawater conditions; also speciation and solvation of uranium's major competing species in seawater, vanadium; (ii) accurately predict the free energies (and thereby stability constants) of binding between uranium and amidoxime; (iii) accurately predict the reaction enthalpies and entropies of binding between uranium and amidoxime as well as between vanadium and amidoxime. Main results and conclusions are as follows. First-principles molecular dynamic simulations found that the structure of the Ca 2 UO 2 (CO 3 ) 3 complex is very stable and that one Ca ion binds to the center UO 2 (CO 3 ) 3 4- anion stronger than the other Ca ion. This finding suggests that using time-resolved EXAFS spectra may confirm the asymmetry in binding of the two Ca ions in the aqueous Ca 2 UO 2 (CO 3 ) 3 complex. To consider the common ions, solvation of the Ca 2 UO 2 (CO 3 ) 3 complex in seawater was simulated by classical molecular dynamics simulation. It was found that the structure of the Ca 2 UO 2 (CO 3 ) 3 complex is very stable in the model seawater. A Na+ ion was found to be closely associated with the Ca 2 UO 2 (CO 3 ) 3 complex by indirectly interacting with one axial oxygen atom of the UO 2 group bridged by a water molecule. In addition, the Na + ion interacts closely with one Ca 2+ ion than the other. The present simulations revealed the key role of common ions such as Na + in impacting the solvation, structure, and apparent charge of the Ca 2 UO 2 (CO 3 ) 3 complex in seawater. The sequential processes of UO 2 2+ binding with the three CO 3 2- groups and the two Ca 2+ ions in pure water and in 0.1 M NaCl were simulated by classical molecular dynamics with both the non-constrained brute-force approach and umbrella sampling. The simulated free energies show excellent agreement with the experiment. The coordination structures of UO 2 (CO 3 ) 3 4- and CaUO 2 (CO 3 ) 3 2- were found to be significantly affected by the presence of Na + ions, leading to a monodentate binding of a carbonate group to U and a water molecule entering the first coordination shell of U. To shed light on the binding between uranyl and glutardiamidoxime (H 2 B), an important model and state-of-the-art ligand for seawater uranium extraction, a suite of computational methods was employed. From molecular dynamics (MD) simulations, it was found that the binding configurations of B 2- with uranyl favor the twofold distorted η 2 binding between the oximate ends (C=N-O-) and U, while HB - prefers the chelating mode for the oximate end with the neutral end being solvated by water. The free energies of sequential ligand binding to form UO 2 B, [UO 2 B 2 ] 2- , and [UO 2 (HB)B] - were simulated with umbrella sampling and very good agreement with the experimental values was achieved, which corroborates the structural insights into the binding mode. State-of-the-art polymeric sorbents employ both amidoximate and carboxylate groups on the side chains to achieve optimal U uptake and selectivity, so we simulated the binding of a model amidoximate–carboxylate bifunctional ligand with uranyl. Classical MD and free-energy simulations in 0.5 M NaCl showed that the carboxylate group binds first to uranyl, leading to a loose intermediate state, and then, the amidoximate group binds, resulting in a more stable and tight chelate state. Binding of the second bifunctional ligand follows a similar process, and the two ligands prefer a trans configuration around the uranyl group. The simulated free energies indicate that the two bifunctional ligands bind with uranyl 14 kcal/mol stronger than the two ligands with only amidoximate groups, confirming an important synergy between amidoximate and carboxylate groups in binding uranyl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Resolving Local Configurational Contributions to X-ray and Neutron Radial Distribution Functions Within Solutions of Concentrated Electrolytes – A Case Study of Concentrated NaOH

Extreme conditions of complex materials often lead to a manifold of local environments that challenge characterization and require new advances at the intersection of modern experimental and theoretical techniques. In this contribution, highly caustic and viscous aqueous NaOD solutions were characterized with a combination of X-ray and neutron radial distribution function (RDF) analyses, molecular dynamics simulations and sub-ensemble analysis. The interpretation of the experimental RDFs extends beyond the traditional understanding of atom-pair interactions and to the individual contributions of specific molecular-scale configurations (below 3.5 Å). This integrated approach yields unique insight into the complex and significant breath of ion solvation and pairing structural motifs present, the experimental sensitivity of RDFs to changes in local geometries, and the challenge of experimentally differentiating the ensemble of all superimposed local environments — a feature of increasing importance within the extreme condition of high ionic strength. This study was facilitated by the development of a chemically resistant sample cell for in-situ neutron total scattering measurements of highly caustic materials at elevated temperature and pressure conditions. Advancing these integrated capabilities is critical for ongoing model validation efforts and for extending the predictive capabilities needed for the extremes of chemical composition, aiding a variety of industrial processes that include advanced nuclear waste remediation and metal refining strategies.

Semrouni, David↗
NASA NTRS1999

Trace Chemistry

The goals of the trace chemistry group were to identify the processes relevant to aerosol and aerosol precursor formation occurring within aircraft gas turbine engines; that is, within the combustor, turbine, and nozzle. The topics of discussion focused on whether the chemistry of aerosol formation is homogeneous or heterogeneous; what species are important for aerosol and aerosol precursor formation; what modeling/theoretical activities to pursue; what experiments to carry out that both support modeling activities and elucidate fundamental processes; and the role of particulates in aerosol and aerosol precursor formation. The consensus of the group was that attention should be focused on SO2, SO3, and aerosols. Of immediate concern is the measurement of the concentration of the species SO3, SO2, H2SO4 OH, HO2, H2O2, O, NO, NO2, HONO, HNO3, CO, and CO2 and particulates in various engines, both those currently in use and those in development. The recommendation was that concentration measurements should be made at both the combustor exit and the engine exit. At each location the above species were classified into one of four categories of decreasing importance, Priority I through IV, as follows: Combustor exit: Priority I species - SO3:SO2 ratio, SO3, SO2, and particulates; Priority II species: OH and O; Priority III species - NO and NO2; and Priority IV species - CO and CO2. For the Engine exit: Priority I species - SO3:SO2 ratio, SO3, SO2,H2SO4, and particulates; Priority II species: OH,HO2, H2O2, and O; Priority III species - NO, NO2, HONO, and HNO3; and Priority IV species - CO and CO2. Table I summarizes the anticipated concentration range of each of these species. For particulate matter, the quantities of interest are the number density, size distribution, and composition. In order to provide data for validating multidimensional reacting flow models, it would be desirable to make 2-D, time-resolved measurements of the concentrations of the above species and, in addition, of the pressure, temperature, and velocity. A near term goal of the experimental program should be to confirm the nonlinear effects of sulfur speciation, and if present, to provide an explanation for them. It is also desirable to examine if the particulate matter retains any sulfur. The recommendation is to examine the effects on SOx production of variations in fuel-bound sulfur and aromatic content (which may affect the amount of particulates formed). These experiments should help us to understand if there is a coupling between particulate formation and SO, concentration. Similarly, any coupling with NOx can be examined either by introducing NOx into the combustion air or by using fuel-bound nitrogen. Also of immediate urgency is the need to establish and validate a detailed mechanism for sulfur oxidation/aerosol formation, whose chemistry is concluded to be homogeneous, because there is not enough surface area for heterogeneous effects. It is envisaged that this work will involve both experimental and theoretical programs. The experimental work will require, in addition to the measurements described above, fundamental studies in devices such as flow reactors and shock tubes. Complementing this effort should be modeling and theoretical activities. One impediment to the successful modeling of sulfur oxidation is the lack of reliable data for thermodynamic and transport properties for several species, such as aqueous nitric acid, sulfur oxides, and sulfuric acid. Quantum mechanical calculations are recommended as a convenient means of deriving values for these properties. Such calculations would also help establish rate constants for several important reactions for which experimental measurements are inherently fraught with uncertainty. Efforts to implement sufficiently detailed chemistry into computational fluid dynamic codes should be continued. Zero- and one-dimensional flow models are also useful vehicles for elucidating the minimal set of species and reactions that must be included in two- and three-dimensional modeling studies.

Radhakrishnan, Krishnan↗
DOE OSTI2020

Current Needs in Radiochemistry

There are a number of current needs in the development of the Pu and Am separations that underpin larger interest areas including nuclear waste issues, nuclear weapons, nuclear forensics, nuclear fuels and similarly related fields. In the large areas of interest, there are a number of sub-fields that include corrosion chemistry, pyrochemistry, solution-phase chemistry, rapid separations, environmental speciation, and other similar efforts. However, in all of the these fields it is important to realize that the underlying chemistry that underpins these efforts is similar in that Pu and Am are bound to ligands that can be in the solid, aqueous, or gas phase and separations occurs due to the preferential binding of the ligand to Pu or Am based on a number of factors that include pH, concentration, oxidation state, etc. In the past, most of the Pu and Am chemistry has occurred in limited scope within universities where micro-scale and bench-top scale chemistry are the dominate investigative theme. At US national laboratory and their international counterparts, larger scale reactions have been attempted such that industrial scale chemistry has been completed. Some of this technology has been harnessed in radiopharmaceuticals and nuclear fuel reprocessing. However, such bench-top research is expensive to due to the safety and security that is concerned when dealing with Pu, and scale-up of reactions is even more resource intensive. In an effort to offset this expensive, it would be of interest to explore the ability to predict selectivity based on the computational modeling of actinide chemistry. This has historically been limited by the abilities of large supercomputers, but with the emphasis on developing exascale computers within the US and similar efforts by other countries this predictive modeling has potential to produce significant results that would yield improved understanding of actinide chemistry.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

An assessment of ocean alkalinity enhancement using aqueous hydroxides: kinetics, efficiency, and precipitation thresholds

Abstract. Ocean alkalinity enhancement (OAE) is a promising approach to marine carbon dioxide removal (mCDR) that leverages the large surface area and carbon storage capacity of the oceans to sequester atmospheric CO2 as dissolved bicarbonate (HCO3-). One OAE method involves the conversion of salt in seawater into aqueous alkalinity (NaOH), which is returned to the ocean. The resulting increase in seawater pH and alkalinity causes a shift in dissolved inorganic carbon (DIC) speciation toward carbonate and a decrease in the surface ocean pCO2. The shift in the pCO2 results in enhanced uptake of atmospheric CO2 by the seawater due to gas exchange. In this study, we systematically test the efficiency of CO2 uptake in seawater treated with NaOH at aquarium (15 L) and tank (6000 L) scales to establish operational boundaries for safety and efficiency in advance of scaling up to field experiments. CO2 equilibration occurred on the order of weeks to months, depending on circulation, air forcing, and air bubbling conditions within the test tanks. An increase of ∼0.7–0.9 mol DIC per mol added alkalinity (in the form of NaOH) was observed through analysis of seawater bottle samples and pH sensor data, consistent with the value expected given the values of the carbonate system equilibrium calculations for the range of salinities and temperatures tested. Mineral precipitation occurred when the bulk seawater pH exceeded 10.0 and Ωaragonite exceeded 30.0. This precipitation was dominated by Mg(OH)2 over hours to 1 d before shifting to CaCO3,aragonite precipitation. These data, combined with models of the dilution and advection of alkaline plumes, will allow the estimation of the amount of carbon dioxide removal expected from OAE pilot studies. Future experiments should better approximate field conditions including sediment interactions, biological activity, ocean circulation, air–sea gas exchange rates, and mixing zone dynamics.

Ringham, Mallory C. (ORCID:0000000348020185)↗
NASA NTRS1996

Numerical simulation of iodine speciation in relation to water disinfection aboard manned spacecraft I. Equilibria

Elemental iodine (I2) is currently used as the drinking water disinfectant aboard the Shuttle Orbiter and will also be incorporated into the water recovery and distribution system for the International Space Station Alpha. Controlled release of I2 is achieved using the Microbial Check Valve (MCV), a flow-through device containing an iodinated polymer which imparts a bacteriostatic residual concentration of approximately 2mg/L to the aqueous stream. During regeneration of MCV canisters, I2 concentrations of approximately 300 mg/L are used. Dissolved iodine undergoes a series of hydrolytic disproportionation and related reactions which result in the formation of an array of inorganic species including: I-, I3-, HOI, OI-, IO3-, HIO3, I2OH-, I2O(-2), and H2OI+. Numerical estimation of the steady-state distribution of inorganic iodine containing species in pure water at 25 degrees C has been achieved by simultaneous solution of the multiple equilibrium expressions as a function of pH. The results are reported herein.

NASA Discipline Environmental Health↗
DOE OSTI2020

Electrochemistry of PUREX: R is for Reduction and Ion Transfer

The PUREX (Plutonium, Uranium, Reduction, EXtraction) process is the foundation of all industrial activities that involve the recycling of plutonium and uranium from used nuclear fuel. With over 70 years of research, engineering, and operations, it can be argued that little opportunity is left for further discovery and new fundamental understandings regarding PUREX, in general, and tri-n-butyl phosphate (TBP), in particular. Through use of an advanced electroanalytical approach, molecular-level insights regarding the back-extraction of Pu(IV) by reductive separation is reported. Studies of model PUREX systems, with lanthanide ions extracted into 20 % TBP in n-dodecane, are used to demonstrate the approach. These studies reveal knowledge pertinent to mass transfer aspects of the biphasic reductive stripping reaction central to the separation of U and Pu in PUREX. Two lanthanide ions—Ce(IV) and Yb(III)—were extracted from their concentrated solutions in aqueous 3 M HNO 3 to expressly provide electrically conducting third phases to facilitate electrochemical data acquisition. To simplify the diabolically-complicated reductive separation aspects of PUREX, the controlled polarization of electrode surfaces was used to reduce Ce(IV) and Yb(III) in place of chemical reductants. Electroanalyses of the biphasic systems with Ce(IV) and Yb(III) using three-phase electrode differential pulse voltammetry demonstrate the effects of mass transfer across water-oil interfaces on the electrode potentials for the reduction to Ce(III) and Yb(II). These reductions trigger nitrate transfer out of the oil phase for Ce and proton transfer into the organic phase for Yb. This is in correlation with the transfer of the charge-neutral ion pairs Ce·3NO 3 and Yb·2NO 3 out of the oil phases. Here, the differences in ion transfer steps reflect significant variations in the equilibrium speciation of multinuclear Ce(IV) and mononuclear Yb(III) solvates of TBP in the oil phases and the loading of the organic phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Iodine Distribution During Evaporation of Hanford Waste Treatment Plant Direct Feed Low Activity Waste Effluent Management Facility Simulant

The Hanford Waste Treatment and Immobilization Plant (WTP) Low Activity Waste (LAW) vitrification facility will generate an aqueous condensate recycle stream from the off-gas system. The baseline plan for disposition of this stream during full WTP operations is to send it to the WTP Pretreatment Facility, where it will be blended with LAW, concentrated by evaporation, and recycled to the LAW vitrification facility. However, during the Direct Feed LAW (DFLAW) scenario, planned disposition of this stream involves concentrating the condensate in a new evaporator at the Effluent Management Facility (EMF) and returning it to the LAW melter. The LAW melter condensate stream will contain components, e.g. halides and sulfates, that are volatile at melter temperatures, have limited solubility in glass waste forms, and present a material corrosion concern. Further, some minor constituents in the melter condensate, such as volatile organic compounds and iodine, have been found to be at least partially volatile in the EMF. The condensate from the EMF will be sent to the Hanford Effluent Treatment Facility (ETF). The presence of some recalcitrant organics and iodine, as 129I, are expected to be problematic in the ETF. In order to better predict and prepare for the distribution of organics and 129I in the flowsheet, it is key to understand their chemistry and experimentally determine their behavior during evaporation. To do this, modeling and testing were performed to examine the speciation of iodine under relevant conditions, as well as determine partitioning of iodine and the organics in the EMF evaporator. This overall program examines the potential treatment and immobilization of the LAW melter condensate stream to enable alternative disposal. The objective of this task was to (1) perform modeling to determine the speciation and fate of iodine species, (2) prepare a simulant of the LAW Melter Off-gas Condensate expected during DFLAW operations, (3) spike in the key volatile organics and non-radioactive iodine, and (4) demonstrate evaporation in order to predict the final composition of the effluents from the EMF evaporator to aid in planning for their disposition. The results of the modeling of iodine speciation indicate that the iodine is expected to remain in the evaporator pot as iodide and possibly iodate ion if the pH is raised to 11.5-12, which is the planned condition for the EMF. Experimental testing confirmed that the iodide remained in the evaporator concentrate and was below detection limits in all condensate and knock-out pot samples, indicating it was not volatile. The organic nitrile compounds that were added to the simulant were absent from the evaporator feed solution, either lost through evaporation, or more likely, decomposition. The acetone and methylene chloride were below the target values in the feed and condensate samples, and were presumably lost to evaporation or decomposition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

NEWTS Integrated Dataset (version 2.0)

The National Energy Water Treatment and Speciation (NEWTS) Integrated Dataset v2.0 provides water researchers, community leaders, regulators, and industry stakeholders with a unified and standardized energy-process wastewater chemistry database. This resource is derived from 39 state and federal entities, and scientific publications, and contains more than 700,000 sample records, many of which also provide geospatial information. The dataset includes chemistry data for several different energy-process wastewater types including produced water, other oil and gas wastewaters, mine drainage, coal ash leachate, power plant wastewater, and geothermal fluids. The NEWTS Integrated Dataset was built to support prudent decision-making, characterization of potential critical mineral sources, and modeling of treatment and valorization options. A subset of this novel resource is also featured on the NEWTS State-Level Database Dashboard. Additional data can be found in the NEWTS EDX Group and the NEWTS Federal Database Dashboard.

AMD↗
DOE OSTI2025

Evaluation of New Additions to OLI Software in Predicting Mercuric and Mercurous Species in Liquid Waste Operations

Speciation of mercury during the pretreatment steps of tank waste processing is critical to successful mercury removal prior to vitrification during Liquid Waste Operations (LWO) at SRS. OLI software has been used to predict mercury speciation and activity throughout LWO. The OLI software operates based on a thermodynamic framework called the Mixed Solvent Electrolyte (MSE) framework. The MSE framework allows prediction in theoretically infinitely dilute to concentrated mixtures (e.g., purely solute solutions). Before modification to the MSE framework databanks, certain critical mercury species were missing in the MSE databank, and some thermodynamic data needed to be updated for the OLI software to accurately predict mercury chemical species in SRS waste tanks. To better reflect streams across LWO, new mercury species were integrated into the MSE database. To evaluate the changes to the OLI MSE framework per the Technical Task Request (TTR) and the Task Technical and Quality Assurance Plan (TTQAP), waste stream compositions from Tanks 38, 43, and Tank 50 decontaminated salt solution (DSS) were used as model inputs. Models were developed and executed using both the old and new databases. Compositional analyses from caustic Tank 50 DSS and caustic Tanks 38 and 43 were used as the input streams. These streams represent the most comprehensive chemical data sets where both mercury and tank constituents were measured together. Results for Tank 50 DSS predict HgO as the predominant species in both databases. Both methyl and dimethyl Hg species are present when the new database is ‘on’ and are not predicted with the new database turned ‘off’. The new database predicts a greater amount of HgO and a greater fraction of it in the solid phase. Pourbaix diagrams (potential vs. pH) generated for each Tank 50 DSS were identical regardless of which database was used. Elemental Hg and HgO were predicted in the water stable region under basic conditions. Tanks 38 and 43 follow similar trends as the Tank 50 DSS models. Unlike Tanks 38 and 50 DSS, the Tank 43 Pourbaix plot shows a region of stability for an aqueous HgOHCO3 - species between approximately pH 7-11. In all streams, when MeHg+ is included in the inputs, the new database predicts aqueous MeHgOH as the dominant species. If elemental or dimethyl mercury is in the waste stream, the new database model predicts they are unchanged and remain in those states and quantities. Additionally, the total mercury values are reported for both the measured input data and the OLI output data for all considered tanks. The summary indicates that the percentage error between the measured and calculated values is less than 1% in all cases The reconciliations and generation of the Pourbaix diagrams for Tank 50 DSS took approximately ten times longer with the new database ‘on’. In addition, over the course of that time, models with the new database ‘on’ were more likely to crash or display an error. Some modest performance improvements were noted when modeling with an i7 processor versus an i5. An example error is found in Appendix A. Furthermore, Appendix B provides V&V for two chemical systems analyzed with the OLI software, results were satisfactory. It is recommended to utilize the new databases (i.e., HCO.ddb and SR-Hg.ddb) in future Savannah River Mission Completion applications of OLI to represent pseudo steady-state. Furthermore, the integration and utilization of the new databases (i.e., HCO.ddb and SR-Hg.ddb) in modeling applications (e.g., Aspen) is also recommended.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Stability of aqueous neodymium complexes in carbonate-bearing solutions from 100–600 °C

Rare earth element exploration requires a quantitative understanding of factors governing their mobilization and economic concentration. However, the behavior of rare earth elements in carbonate- bearing hydrothermal fluids associated with carbonatite-hosted deposits is poorly understood, and conflicting mechanisms of rare earth transport by anionic ligands and alkali behavior have been described. Here, we report quantitative data to characterize the role of carbonate-bearing solutions in the hydrothermal mobilization of neodymium. Solubility studies of neodymium phosphate were performed at temperatures ranging from 100 to 600 °C in carbonate-bearing solutions. The thermodynamic data determined for the predominant complex were used to model the separation of neodymium from thorium in a simple flow-through system based on fluid and mineral compositions characteristic of carbonatite deposits. Our data suggest that neodymium transport is controlled by the stability of the carbonate species NdCO 3 OH o , and at temperatures of 500–600 °C, the concentrations of neodymium in solutions can reach ~1000 ppm.

58 GEOSCIENCES↗
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2025

Spectrophotometric determination of the stability of La hydroxyl complexes at near neutral to alkaline pH from 25 to 75 °C

The hydrolysis of rare earth elements (REE) potentially controls their mobility during fluid-rock interaction in a broad range of pH and temperature conditions. However, there is still a lack of thermodynamic data for modeling accurately the stability of REE hydroxyl complexes in hydrothermal aqueous fluids. Here, in this study, UV–Vis spectrophotometric experiments were conducted from 25 to 75 °C in near-neutral to alkaline NaOH-bearing aqueous solutions with varying lanthanum (La) concentrations (0 to ∼0.23 mmol/kg). The color indicator m-cresol purple was used to determine in situ pH and derive the average OH− ligand number ($\overrightarrow{n}$) and formation constants for the La hydroxyl complexes (LaOH 2+ , La(OH) 2 + , and La(OH) 3 0 ). From 25 to 50 °C, $\overrightarrow{n}$ ranges between ∼1 and 2 at pH from 7.0 to 9.3. At 75 °C, $\overrightarrow{n}$ ranges between ∼1.5 and 3 at pH from 6.3 to 8.8. These results suggest the predominance of LaOH 2+ and La(OH) 2 + complexes from 25 to 50 °C, and an increased predominance of La(OH) 3 0 at 75 °C. The cumulative formation constants (β n °, n = 1 to 3) are derived for the reaction La 3+ + nOH − = La(OH) n 3-n , and fitted between 25 and 250 °C by combining the UV–Vis and literature solubility data. The resulting logβ n ° are expressed as function of temperature (T in Kelvin): logβ 1 ° = −1.786 + 0.0133 T + 1.049·10 3 /T; logβ 2 ° = −5.797 + 0.0267 T + 2.713·10 3 /T; logβ 3 ° = 6.435 + 0.0223 T + 512.7/T. A comparison between these new fits and existing extrapolations using the Helgeson-Kirkham-Flowers equation of state indicates significant differences in the predicted hydrolysis of La. The latter extrapolations should therefore be updated for the hydrolysis of REE.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Mobility of heavy rare earth elements in magmatic-hydrothermal systems: Experimental determination of DyPO 4 solubility in supercritical fluids

Dysprosium (Dy) is part of the heavy rare earth elements (REE) and a critical component in the production of high-performance super magnets used in wind turbines, electric vehicles, and other green technologies. The HREE phosphate, xenotime, is commonly associated with hydrothermal systems, where the mobility, fractionation, and deposition of REE are controlled by the stability of aqueous complexes and solubility of REE minerals which depend on temperature, pressure, and fluid composition (salinity, pH, and ligand availability). Here, we conducted a series of synthetic DyPO 4 solubility experiments in NaCl-HCl-bearing aqueous solutions (pH of 2 at 25 °C) from 500 to 700 °C and 1.5 kbar using externally heated pressure vessels. DyPO 4 displays a strong prograde solubility with temperature across all studied fluid salinities (0–0.5 m NaCl). The DyPO₄ solubility isotherms at 500 °C, 600 °C, and 700 °C show a complex relationship with salinity where DyPO₄ solubility is highest (39.2–1552 ppm at 500–700 °C) in 0 m NaCl solutions, decreases strongly (21–331 ppm at 500–700 °C) in the 0.06 m NaCl solutions and remains relatively constant in the 0.25 m (42.8–248 ppm) and 0.5 m NaCl (40.2–124 ppm) solutions. Thermodynamic modeling suggests that pH increases as a function of temperature and salinity, and controls solubility due to the stability of chloride complexes (i.e., DyCl +2 and DyCl 2 + ) at 500–600 °C, transitioning to hydroxyl complexes (i.e., Dy(OH) 3 0 ) dominating at 700 °C and relatively alkaline pH. Comparison with NdPO 4 solubility data indicates that the light REE in monazite is significantly more soluble than the heavy REE in xenotime in saline fluids at high temperature (500–700 °C), implying efficient light/heavy REE fractionation in supercritical fluids. The lower solubility of DyPO 4 relative to NdPO 4 suggests that at high temperatures, acidic (pH 2–4), low-salinity (∼5 wt% NaCl) fluids will favor the precipitation of xenotime over monazite as it will cool from 700 °C to below 500 °C. These observations align with paragenetic mineral assemblages in NYF pegmatites (e.g., Baveno, Italy), where xenotime commonly immobilizes heavy REE during the high-temperature hydrothermal stages.

58 GEOSCIENCES↗
DOE OSTI2024

Hydrothermal solution calorimetry in acidic aqueous solutions and revisiting the standard partial molal thermodynamic properties of Nd 3+ from 25 to 300 °C

The mobility of rare earth elements (REE) can be predicted in aqueous fluids using geochemical modeling but the accuracy of these models strongly depends on the availability of robust thermodynamic properties for the REE aqueous species. The REE 3+ aqua ions are important in the derivation of the formation constants of all the major REE complexes including the chloride, sulfate, and fluoride species which predominate in many hydrothermal-magmatic systems. However, the thermodynamic properties of the REE 3+ aqua ions are still commonly derived from the Helgeson-Kirkham-Flowers (HKF) equation of state parameters tabulated several decades ago. The standard state thermodynamic properties at reference conditions (25 °C and 1bar) and their extrapolations to high temperature need to be verified, if not revised, based on hydrothermal experiments. In this study, the enthalpy of solution was measured for synthetic Nd hydroxide from 25 to 150 ºC to retrieve the standard partial molal thermodynamic properties of Nd 3+ as a function of temperature. The experiments were conducted in aqueous perchloric acid based solutions with starting pH of 2 and varying ionic strength (0.01 to 0.09 mol/kg NaClO 4 ). The standard partial molal enthalpy of formation (Δ f H°) of Nd 3+ derived from the experimental study displays differences of up to 10 kJ/mol compared to the enthalpy values derived from the HKF equation of state in the studied temperature range. These inaccuracies are resolved by adjusting the standard partial molal Gibbs energy of formation (Δ f G°) of Nd 3+ at 25 °C and 1 bar from -672.0 to -679.7 kJ/mol. The heat capacity function (C p °) derived between 25 and 150 ºC can be described by: C p ° = a 0 + a 1 ·T + a 2 ·T -2 , with a 0 = 1256, a 1 = -2.68, a 2 = -55.56·10 6 and T in Kelvin. A set of recommended thermodynamic properties is provided for the Nd 3+ aqua ions and corrections are provided for the chloride and fluoride species to remain internally consistent with the experimentally derived properties. These results allow predicting accurately the solubility of monazite between 25 and 300 ºC. Before these corrections, the properties for the Nd 3+ aqua ions derived from the HKF parameters resulted in up to ~1.5 orders of magnitude lower monazite solubility than determined experimentally. Here, a revision of the REE +3 aqua ions properties is necessary to accurately predict the mobility of REE in hydrothermal acidic solutions.

58 GEOSCIENCES↗