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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A Uranium Bioremediation Reactive Transport Benchmark

A reactive transport benchmark problem set has been developed based on in situ uranium bio-immobilization experiments that have been performed at a former uranium mill tailings site in Rifle, Colorado, USA. Acetate-amended groundwater stimulates indigenous microorganisms to catalyze the reduction of U(VI) to a sparingly soluble U(IV) mineral. The interplay between the flow, acetate loading periods and rates, microbially-mediated and geochemical reactions leads to dynamic behavior in metal- and sulfate-reducing bacteria, pH, alkalinity, and reactive mineral surfaces. The benchmark is based on an 8.5 m long one-dimensional model domain with constant saturated flow and uniform porosity. The 159-day simulation introduces acetate and bromide through the upgradient boundary in 14-day and 85-day pulses separated by a 10 day interruption. Acetate loading is tripled during the second pulse, which is followed by a 50 day recovery period. Terminal electron accepting processes for goethite, phyllosilicate Fe(III), U(VI), and sulfate are modeled using Monod-type rate laws. Major ion geochemistry modeled includes mineral reactions, as well as aqueous and surface complexation reactions for UO2++, Fe++, and H+. In addition to the dynamics imparted by the transport of the acetate pulses, U(VI) behavior involves the interplay between bioreduction, which is dependent on acetate availability, and speciation-controlled surface complexation, which is dependent on pH, alkalinity and available surface complexation sites. The general difficulty of this benchmark is the large number of reactions (74), multiple rate law formulations, a multisite uranium surface complexation model, and the strong interdependency and sensitivity of the reaction processes. Results are presented for three simulators: HYDROGEOCHEM, PHT3D, and PHREEQC.

bioremediation↗
DOE OSTI2024

UV-Vis spectrophotometric determination of rare earth elements (REE) speciation at near-neutral to alkaline pH. Part I: m-cresol purple properties from 25-75 °C and Er hydrolysis

The speciation and mobility of rare earth elements (REE) strongly depends on pH which controls the formation of charged aqueous hydroxyl species. The latter potentially play an important role in controlling heavy REE adsorption on clay minerals in near-neutral to alkaline waters such as in regolith-hosted REE mineral deposits. However, accurate REE hydrolysis constants are needed for developing geochemical models that can predict the role of these charged species in natural systems. Here, we develop a robust experimental UV-Vis spectrophotometric method using m-cresol purple to determine in situ pH from 25 to 75 °C. This method is used to derive the average ligand number and hydrolysis constants of erbium (Er) at 25 °C in aqueous solutions with low ionic strength (≤ 0.001 mol/L) at pH from ~7 to 9.5 and in the presence of Er concentrations from 0 to 0.057 mM. The average ligand number ranges between 1 and 3 indicating that Er(OH) 2+ , Er(OH) 2 + and Er(OH) 3 0 control speciation in the experiments. The logarithm of the Er hydrolysis constants (log*β n ° , n= 1 to 3) derived at infinite dilution for the reaction Er 3+ + nH 2 O = Er(OH) n 3-n + nH + are: *β 1 ° = –7.22 ± 0.10, *β 2 ° = –14.52 ± 0.08, *β 3 ° = –23.24 ± 0.04. Implementation of these experimental data into a geochemical model indicates that the Er(OH) 2 + and Er(OH) 3 0 species are both stable in a much wider pH range than previously predicted. Consequently, the positively charged REE hydroxyl complexes can potentially control the fractionation of light vs. heavy REE via adsorption as observed in the formation of certain regolith-hosted REE deposits.

58 GEOSCIENCES↗
DOE OSTI2018

Geochemical Simulations of Plutonium Redox Speciation at Increasing Ionic Strengths - 18567

In salt repositories like the Waste Isolation Pilot Plant (WIPP), reducing conditions are expected, supporting two aqueous species of plutonium: Pu(III) and Pu(IV). Of these states, Pu(III) is more soluble than Pu(IV). Because Pu behavior is uncertain in WIPP brines under reducing conditions, performance assessment calculations assume a 50/50 partitioning between the two, justified through calculations using dilute solutions. Recent experimental evidence suggests an increasing importance for Pu(III). To reexamine this assumption, we modeled Pu speciation in NaCl and dilute WIPP brine solutions equilibrated with a Pu solid, including PuO{sub 2} and Pu(OH){sub 3}. To maintain reducing conditions, sulfate and sulfide were added. Models were performed using PHREEQC and the Specific Ion Theory (SIT) database on dilutions ranging from 0.1 M to 3 M. Models were compared to calculations using the Pitzer database with PuO{sub 2} and with analogs Am(III) and Th(IV), all up to 11 M. Results indicate that under reducing conditions, Pu(III) concentrations increase with ionic strength and can be up to 80 fold higher than Pu(IV). While concentrations varied due to differences in solubility for each solid, Pu(IV) concentrations remained nearly constant in all solutions. Analogs showed similar trends, though Am(III) is 3 times more soluble than Pu(III) while Th(IV) is 6 times more soluble than Pu(IV). Pu(IV) concentrations using Pitzer calculations were similar to those with SIT. Our study indicates that assuming equivalent partitioning of Pu species may under-represent concentrations during a release. They also show that in PA calculations there is a need to consider redox state at high ionic strength rather than adopting an assumption of equal partitioning between the two charges. (authors)

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

Polynuclear Speciation of Trivalent Cations near the Surface of an Electrolyte Solution

Despite long-standing efforts, there is no agreed upon structural model for electrolyte solutions at air-liquid interfaces. We report the simultaneous detection of the near-surface and bulk coordination environments of a trivalent metal cation (europium) in an aqueous solution by use of X-ray absorption spectroscopy. Within the first few nanometers of the liquid surface, the cations exhibit oxygen coordination typical of inner-sphere hydration of an aquated Eu3+ cation. Beyond that, outer-sphere ion-ion correlations are observed that are otherwise not present in the bulk electrolyte. The combination of near-surface and bulk sensitivities to probe metal ion speciation in electrolyte solutions is achieved by detecting electron-yield and X-ray fluorescence signals from an inverted pendant drop. The results provide new knowledge about the near-surface chemistry of aqueous solutions of relevance to aerosols and ion transport processes in chemical separations and biological systems.

Electrolytes↗
DOE OSTI2026

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗
DOE OSTI2025

National Energy Water Treatment & Speciation (NEWTS): A Water & Critical Mineral Database and Dashboard

The scarcity of water resources, the need for beneficial water reuse, and the challenges of wastewater treatment are becoming increasingly pressing in economic, social, and environmental domains. Addressing these concerns requires effective treatment strategies to manage wastewater streams and tackle environmental and economic issues. Furthermore, the recovery of critical minerals from the waste streams associated with energy production holds the promise of offsetting treatment costs and securing local sources of valuable minerals. However, relevant data on these waste streams are dispersed and challenging to locate. The process of ingesting such data into modeling software often involves multiple steps, requiring data restructuring to meet software-input requirements. The non-standardized reporting of water data makes data aggregation and reformatting a time-consuming process. Additionally, essential attributes necessary for modeling water treatment and mineral scale formation are frequently missing. Moreover, data gaps vary depending on the region of interest. Consequently, there is a pressing need for high-quality energy-water composition data that can be easily imported into water chemistry modeling software. To address this need, the National Energy Technology Laboratory has created the National Energy Water Treatment and Speciation (NEWTS) Database and Dashboard—a free online tool catering to community leaders and water researchers. NEWTS facilitates a comprehensive understanding of the composition of energy-related wastewater streams in the United States. The datasets provide detailed concentrations and speciation of major and minor aqueous compounds in energy-related wastewater streams, including power plant leachate, acid mine drainage, brackish water, and oil and gas produced water across the United States. Many of the aqueous species are critical minerals (Li, REEs) in high demand to modernize the world’s energy infrastructure. Many of the datasets also contain volumetric flow-rates needed to model the treatment and reuse scenarios in advanced aqueous chemistry software programs. The NEWTS Database and Dashboard offer public access to hitherto challenging-to-access datasets, presented in a standardized format that is tailored for easy input into aqueous chemistry modeling software. By performing the work needed to transform dispersed, disparate data sources into unified, model-ready datasets, NEWTS serves as an essential resource in advancing water treatment research and sustainable water resource management.

produced water management↗
DOE OSTI2023

National Energy Water Treatment & Speciation (NEWTS): A Water & Critical Minerals Database and Dashboard

The poster, presented at the American Geophysical Union (AGU) 2023 Conference, highlights NETL's National Energy Water Treatment & Speciation (NEWTS) Database & Dashboard. This poster reviews the work done to aggregate disparate water quality data for energy-related wastewater streams, to process and integrate the datasets, and enable easy integration with aqueous chemistry modeling software to support water management and beneficial reuse decision-making.

Wenzlick, Madison↗
DOE OSTI2015

Effects of aqueous uranyl speciation on the kinetics of microbial uranium reduction

The ability to predict the success of the microbial reduction of soluble U(VI) to highly insoluble U(IV) as an in situ bioremediation strategy is complicated by the wide range of geochemical conditions at contaminated sites and the strong influence of aqueous uranyl speciation on the bioavailability and toxicity of U(VI) to metal-reducing bacteria. In order to determine the effects of aqueous uranyl speciation on uranium bioreduction kinetics, incubations and viability assays with Shewanella putrefaciens strain 200 were conducted over a range of pH and dissolved inorganic carbon (DIC), Ca 2+ , and Mg 2+ concentrations. A speciation-dependent kinetic model was developed to reproduce the observed time series of total dissolved uranium concentration over the range of geochemical conditions tested. The kinetic model yielded the highest rate constant for the reduction of uranyl non-carbonate species (i.e., the ‘free’ hydrated uranyl ion, uranyl hydroxides, and other minor uranyl complexes), indicating that they represent the most readily reducible fraction of U(VI) despite being the least abundant uranyl species in solution. In the presence of DIC, Ca 2+ , and Mg 2+ is suppressed during the formation of more bioavailable uranyl non-carbonate species and resulted in slower bioreduction rates. At high concentrations of bioavailable U(VI), however, uranium toxicity to S. putrefaciens inhibited bioreduction, and viability assays confirmed that the concentration of non-carbonate uranyl species best predicts the degree of toxicity. The effect of uranium toxicity was accounted for by incorporating the free ion activity model of metal toxicity into the bioreduction rate law. These results demonstrate that, in the absence of competing terminal electron acceptors, uranium bioreduction kinetics can be predicted over a wide range of geochemical conditions based on the bioavailability and toxicity imparted on U(VI) by solution composition. Finally, these findings also imply that the concentration of uranyl non-carbonate species, despite being extremely low, is a determining factor controlling uranium bioreduction at contaminated sites.

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

Optimizing CO 2 -Loaded Aqueous Amine Solutions for Higher Electrocatalytic CO 2 Reduction Activity

The activity of aqueous-based carbon dioxide reduction (CO 2 R) reactions is often limited by the solubility of CO 2 . The addition of amines can increase the total dissolved carbon in water through the formation of bicarbonate and carbamate species, which has been used to a great effect to capture CO 2 from dilute streams. Here, in this study, we explore the effect of 12 primary and secondary amines of varying Brønsted basicity, steric profile, and hydrogen-bonding capabilities on the aqueous CO 2 R to CO activity of a molecular Ni(cyclam)Cl 2 catalyst with a Hg electrode. Addition of some of the amines results in greater activity and selectivity for CO production compared to equivalent aqueous solutions without added amines. Under optimal conditions (0.4 M 3-amino-propionitrile), there is an over sevenfold increase in partial current density and greater selectivity for CO compared to equivalent conditions with no amine. Interestingly, the increase in activity did not correlate to any single property across the 12 amines. To elucidate the effect of the amine additives on catalysis, we used vapor–liquid equilibrium modeling (VLE), 13 C NMR spectroscopy, and computational analysis to determine the carbon speciation of the solutions. These results indicate that for amines without ethylalcohol functionalities, CO 2 R activity correlates with carbamate concentration, which is in turn governed by amine basicity and steric effects. However, this correlation does not persist for amines with ethylalcohol functionalities, which can form more stable carbamates through intramolecular-hydrogen bonding. These studies demonstrate that amine additives can enhance aqueous CO 2 R activity and selectivity and describe amine properties that lead to these higher performance metrics.

Amines↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗
DOE OSTI2018

Solubility of U(VI) in chloride solutions. III. The stable oxides/hydroxides in MgCl 2 systems: Pitzer activity model for the system UO 2 2+ –Na + –K + –Mg 2+ –H + –OH – –Cl – –H 2 O(l)

Here we have developed new chemical, thermodynamic and activity models for the system UO 2 2+ –Na + –K + –Mg 2+ –H + –OH – –Cl – –H 2 O(l) within the Pitzer approach. The new thermodynamic model is based on previously reported data treated within the SIT approach for NaCl and KCl systems, as well as on new experimental data determined in this work for the MgCl 2 system. The solubility of uranium(VI) was studied in 0.01–5.15 mol·kg w –1 MgCl 2 solutions at pH m = 4.1–9.7 (with pH m = –log [H + ]). Experiments were performed under Ar atmosphere at T = (22 ± 2) °C. In all investigated systems, the solubility of U(VI) is controlled by metaschoepite, UO 3 ·2H 2 O(cr). In contrast to previously investigated NaCl and KCl systems, no ternary Mg–U(VI)–OH(s) solid phases formed in alkaline MgCl 2 solutions within the timeframe of this study (≤ 200 days). A very significant increase in the solubility (up to 3 log 10 -units) is observed in acidic to near-neutral pH m conditions when increasing MgCl 2 concentration from 0.01 to 5.15 mol·kg w –1 , which reflects the strong ion interaction processes taking place in concentrated MgCl 2 brines. The solubility of UO 3 ·2H 2 O(cr) in the investigated NaCl, KCl and MgCl 2 solutions is well described with the solubility and hydrolysis constants recommended by Altmaier et al., (2017) and NEA–TDB, and a Pitzer activity model derived in the present work. The latter model considers experimental data reported in the present study and available in the literature for NaCl, KCl and MgCl 2 systems (solubility, potentiometric and spectroscopic data), in combination with well-stablished estimation methods and correlations with SIT coefficients. Finally, chemical, thermodynamic and Pitzer activity models provided in this work for the system UO 2 2+ –Na + –K + –Mg 2+ –H + –OH – –Cl – –H 2 O(l) accurately describe all evaluated datasets, and represent an adequate tool for the calculation of U(VI) solubility and aqueous speciation in a variety of geochemical conditions including concentrated brine systems of relevance in salt-based repositories for nuclear waste disposal.

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

Probing the Mechanism of Selective Phosphate Adsorption from Wastewater Using Aqueous and Synchrotron X-ray Characterization

Ion exchange shows promise for recovering phosphate from wastewater as value-added products, but requires high phosphate selectivity to compete with conventional treatment. Hybrid anion exchange (HAIX) resins, which contain non-selective basic functional groups and selective iron oxide nanoparticles (FeOnp), can effectively remove phosphate from wastewater. However, knowledge gaps remain regarding the mechanisms of phosphate selectivity and influence of competing ions, hindering needed efforts to model adsorption dynamics and design scalable adsorption processes for varying wastewaters. To address these gaps, we integrated aqueous-phase adsorption analysis with solid-phase, synchrotron-based X-ray characterization; this integration facilitated elucidation of the distribution and speciation of iron, phosphate, and competing anions on HAIX resins. We compared a quaternary ammonium-functionalized HAIX resin (SBA) to a tertiary amine version (WBA) to determine the role of functional groups. X-ray radiography revealed differences in FeOnp speciation (goethite vs. ferrihydrite) and distribution (peripheral vs. homogeneous) between the resins, resulting in varied phosphate affinity and intraparticle diffusion resistance. Using micro-X-ray fluorescence (μ-XRF) and micro-X-ray absorption near-edge structure (μ-XANES) spectroscopy, we identified differences in where and how phosphate binds across resin types and wastewaters. Across wastewater compositions, FeOnp sites in WBA contribute more to phosphate adsorption than in SBA, possibly due to variations in Fe distribution and speciation. Phosphate adsorption densities calculated from quantitative μ-XRF maps matched those from aqueous analysis, demonstrating the effectiveness of this integrated approach. Altogether, results demonstrate the use of synchrotron-based X-ray characterization for investigating adsorption mechanisms and advance HAIX as a phosphate recovery technology from complex wastewaters.

Nutrient recovery↗
DOE OSTI2021

Raman Spectroscopy Coupled with Chemometric Analysis for Speciation and Quantitative Analysis of Aqueous Phosphoric Acid Systems

Complex chemical systems that exhibit varied and matrix-dependent speciation are notoriously difficult to monitor and characterize on-line and in real-time. Optical spectroscopy is an ideal tool for in situ characterization of chemical species that can enable quantification as well as species identification. Chemometric modeling, a multivariate method, has been successfully paired with optical spectroscopy to enable measurement of analyte concentrations even in complex solutions where univariate methods such as Beer’s law analysis fail. Here, Raman spectroscopy is used to quantify the concentration of phosphoric acid and its three deprotonated forms during a titration. In this system, univariate approaches would be difficult to apply due to multiple species being present simultaneously within the solution as pH is varied. Locally Weighted Regression (LWR) modeling was used to determine phosphate concentration from spectral signature. LWR results, in tandem with Multivariate Curve Resolution modeling, provide direct measurement of the concentration of each phosphate species using only the Raman signal. Furthermore, results are presented within the context of fundamental solution chemistry, including Pitzer equations to compensate for activity coefficients and non-idealities associated with high ionic strength systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Photoreduction of gaseous oxidized mercury changes global atmospheric mercury speciation, transport and deposition

Anthropogenic mercury (Hg(0)) emissions oxidize to gaseous Hg(II) compounds, before deposition to Earth surface ecosystems. Atmospheric reduction of Hg(II) competes with deposition, thereby modifying the magnitude and pattern of Hg deposition. Global Hg models have postulated that Hg(II) reduction in the atmosphere occurs through aqueous-phase photoreduction that may take place in clouds. Here we report that experimental rainfall Hg(II) photoreduction rates are much slower than modelled rates. We compute absorption cross sections of Hg(II) compounds and show that fast gas-phase Hg(II) photolysis can dominate atmospheric mercury reduction and lead to a substantial increase in the modelled, global atmospheric Hg lifetime by a factor two. Models with Hg(II) photolysis show enhanced Hg (0) deposition to land, which may prolong recovery of aquatic ecosystems long after Hg emissions are lowered, due to the longer residence time of Hg in soils compared with the ocean. Fast Hg(II) photolysis substantially changes atmospheric Hg dynamics and requires further assessment at regional and local scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Development of a Trace Detection Method for Mercury Speciation in Hanford Tank Waste - 19014

The U.S. Department of Energy's (DOE) Hanford Site in southeast Washington State stores 56 million gallons of chemical and radioactive waste stored in underground tanks. Mercury species in liquid and solid fractions of Hanford tank waste pose a risk to the implementation and execution of DoE's River Protection Project mission due to mercury concentration limits imposed under waste acceptance criteria, toxicity characteristic leaching procedure restrictions for secondary waste disposal, and air emission permits. Although there is currently no processing of tank waste at Hanford, waste processing is imminent, and representative mercury inventory tracking will be important to anticipate and mitigate stack emissions and secondary liquid waste disposal issues. A methodology to accurately determine total mercury and mercury speciation in Hanford tank waste, while maintaining low mercury detection limits, currently does not exist at Hanford. Moreover, the mercury species in the waste solids are unknown and may be waste type-specific. The development of an analytical method to separate and speciate mercury in tank waste was undertaken by the Pacific Northwest National Laboratory (PNNL) while leveraging an existing non-radiological, ultra-trace mercury detection technique frequently performed at PNNL's Marine Sciences Laboratory in Sequim, Washington. The ultimate goal of this work is to develop an analytical method to obtain detailed mercury speciation information and update flowsheet modeling assumptions related to mercury partitioning at Hanford. The first phase of testing was completed in the development of an analytical methodology for the low-level determination of total mercury, elemental mercury, methylmercury, and dimethylmercury in Hanford tank waste. This paper presents Phase I testing which focused on determining whether mercury species can be preconcentrated from aqueous Hanford tank waste solution onto a solid substrate, while simultaneously eliminating or very significantly reducing the parent radionuclide activity. This was demonstrated by developing protocols for the separation and preconcentration of total mercury and methylmercury onto solid substrates, then measuring the substrates for radiological content at PNNL's Radiochemical Processing Laboratory on the Hanford Site. The mercury content was successfully separated onto stable substrates and the radionuclide activity reduction was significant enough that the preconcentrated samples could be subjected to quantification using far less restrictive handling and radiological controls than would be required for the parent sample (i.e., contamination area fume hood vs. hot cell). This mercury separation method required no sample dilution, which enables maximum mercury detectability, while still allowing mercury speciation. (authors)

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

Uranium (VI) Adsorbate Structures on Portlandite [Ca(OH)2] Type Surfaces Determined by Computational Modelling and X-ray Absorption Spectroscopy

Portlandite [Ca(OH) 2 ] is a potentially dominant solid phase in the high pH fluids expected within the cementitious engineered barriers of Geological Disposal Facilities (GDF). This study combined X-ray Absorption Spectroscopy with computational modelling in order to provide atomic-scale data which improves our understanding of how a critically important radionuclide (U) will be adsorbed onto this phase under conditions relevant to a GDF environment. Such data are fundamental for predicting radionuclide mass transfer. Surface coordination chemistry and speciation of uranium with portlandite [Ca(OH) 2 ] under alkaline groundwater conditions (ca. pH 12) were determined by both in situ and ex situ grazing incidence extended X-ray absorption fine structure analysis (EXAFS) and by computational modelling at the atomic level. Free energies of sorption of aqueous uranyl hydroxides, [UO 2 (OH) n ] 2–n (n = 0–5) with the (001), (100) and (203) or (101) surfaces of portlandite are predicted from the potential of mean force using classical molecular umbrella sampling simulation methods and the structural interactions are further explored using fully periodic density functional theory computations. Although uranyl is predicted to only weakly adsorb to the (001) and (100) clean surfaces, there should be significantly stronger interactions with the (203/101) surface or at hydroxyl vacancies, both prevalent under groundwater conditions. The uranyl surface complex is typically found to include four equatorially coordinated hydroxyl ligands, forming an inner-sphere sorbate by direct interaction of a uranyl oxygen with surface calcium ions in both the (001) and (203/101) cases. In contrast, on the (100) surface, uranyl is sorbed with its axis more parallel to the surface plane. The EXAFS data are largely consistent with a surface structural layer or film similar to calcium uranate, but also show distinct uranyl characteristics, with the uranyl ion exhibiting the classic dioxygenyl oxygens at 1.8 Å and between four and five equatorial oxygen atoms at distances between 2.28 and 2.35 Å from the central U absorber. These experimental data are wholly consistent with the adsorbate configuration predicted by the computational models. These findings suggest that, under the strongly alkaline conditions of a cementitious backfill engineered barrier, there would be significant uptake of uranyl by portlandite to inhibit the mobility of U(VI) from the near field of a geological disposal facility.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tonian Carbonates Record Phosphate-Rich Shallow Seas

The early-middle Neoproterozoic is thought to have witnessed significant perturbations to marine P cycling, in turn facilitating the rise of eukaryote-dominated primary production. However, with few robust constraints on aqueous P concentrations, current understanding of Neoproterozoic P cycling is generally model-dependent. To provide new geochemical constraints, we combined microanalytical data sets with solid-state Nuclear Magnetic Resonance, synchrotron-based X-ray Absorption Near Edge Structure spectroscopy, and micro-X-ray Fluorescence imaging to characterize the speciation and distribution of P in Tonian shallow-water carbonate rocks. These data reflect shallow water phosphate concentrations 10–100× higher than modern systems, supporting the hypothesis that tectonically-driven influxes in P periodically initiated kinetically-controlled CaCO 3 deposition, in turn destabilizing marine carbonate chemistry, climate, and nutrient inventories. Alongside these observations, a new compilation and statistical analysis of mudstone geochemistry data indicates that, in parallel, C org and P burial increased across later Tonian continental margins until becoming decoupled at the close of the Tonian, implicating widespread N-limitation triggered by increasing atmospheric O 2 .

58 GEOSCIENCES↗
DOE OSTI2018

Quantum Chemical Approach for Calculating Stability Constants of Mercury Complexes

Stability constants are central to the multi-scale modeling of the thermodynamic speciation, cycling, and transport of mercury (Hg) and other contaminants in aquatic environments. However, for Hg, experimental values for many relevant complexes are not available, and for others can span ranges in excess of 10 log units. The missing data and large uncertainties mean there are significant knowledge gaps in predictions of thermodynamic speciation. As an alternative to experimental measurements, thermodynamic quantities can be calculated with quantum chemical methods. Among these, density functional theory (DFT) with a polarizable continuum solvent combines accuracy with practicability. In this paper, we present an accurate and quick approach in which we use DFT with continuum solvation to calculate stability constants of Hg complexes with inorganic and low molecular-weight organic ligands in aqueous solution. Specifically, we use the M06/[SDD]6-31+G(d,p) level of theory in combination with a modified version of the SMD solvent model in which the solute radii are re-optimized with a scaled solvent-accessible surface (sSAS) approach. For the set of 37 Hg complexes used for optimization, which contain environmentally relevant functional groups and have reliable experimental stability constants, we obtain a mean unsigned error of 1.4 log units. Testing the method on an independent set of 12 Hg complexes reproduces the experimental stability constants to a mean unsigned error of 1.6 log units. Finally, this approach is a substantial step toward generally applicable rapid stability constant derivation for a wide range of Hg complexes, including those present in dissolved organic matter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗