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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Proton Chelating Ligands Drive Improved Chemical Separations for Rhodium

Current methods for the extraction of rhodium carry the highest carbon footprint and worst pollution metrics of all of the elements used in modern technological applications. Improving upon existing methods is made difficult by the limited understanding of the molecular-level chemistry occurring in extraction processes, particularly in the hydrometallurgical separation step. While many of the precious metals can be separated by solvent extraction, there currently exist no commercial extractants for Rh. This is due to its complicated mixed speciation upon leaching into hydrochloric acid, which gives rise to difficulties in designing effective reagents for solvent extraction. Herein we show that the diamidoamine reagent N-n-hexylbis(N-methyl-N-n-octylethylamide)amine transports Rh(III) from aqueous HCl into an organic phase as the monoaquated dianion [RhCl 5 (H 2 O)] 2– through the formation of an outer-sphere assembly; this assembly has been characterized by experimentation (slope analysis, FT-IR and NMR spectroscopy, EXAFS, SANS, and ESI-MS) and computational modeling. The paper demonstrates the importance of applying a broad range of techniques to obtain a convincing mode of action for the complex processes involved in anion recognition in the solution phase. A consistent and comprehensive understanding of how the ligand operates to achieve Rh(III) selectivity over the competitor anion Cl – has emerged. This knowledge will guide the design of extractants and thus offers promise for improving the sustainability of metal extraction from both traditional mining sources and the recycling of secondary source materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neptunium (V) Solubility in the Presence and Absence of WIPP Relevant Ligands

CRA (Compliance Recertification Application)-2024 will update the oxidation state specific model used to calculate actinide solubility in the WIPP (Waste Isolation Pilot Plant). The pentavalent Np(V) ion is the most stable ion in solution within the An(V) series. This valence stability of Np(V) in aqueous solution makes Np a useful representative element for studying the aqueous complex chemistry of all An(V), but in the predicted WIPP scenario, this is the only TRU component that increases in the inventory with time as it is the main daughter product of 241 Am decay (alpha, n) so its relative importance increases throughout the 10,000-year performance period of the WIPP, it is only used to represent itself. This work is focused on confirming organic, borate and carbonate effects on solubility of Np(V) in WIPP simulated brines and 5 M NaCl and reassessing the predicted model speciation to strengthen the overall WIPP An(V) model.

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

Distinct Mechanisms and Hydricities of Cp*Ir-Based CO 2 Hydrogenation Catalysts in Basic Water

Transition metal-catalyzed reversible hydrogenation of CO 2 to formate in aqueous solutions under ambient conditions is an attractive environmentally friendly strategy for the storage and transportation of H 2 in a liquid chemical form. Here, mechanistic details for the CO 2 hydrogenation by a series of pentamethylcyclopentadienyl iridium(III) (Cp*Ir) complexes with picolinamidate ligands are investigated through density functional theory and kinetic isotope effects (KIEs) calculations. Additionally, adopting a speciation approach at a pH of 8.2, CO 2 hydrogenation mechanisms of 37 distinct protonation states and conformers of Cp*Ir type complexes are investigated. This imposes a challenging test for the computational modeling in terms of providing a consistent correlation with experimental kinetics and KIE studies of multiple catalysts instead of a single complex. Overall, H 2 heterolysis to generate an iridium hydride intermediate was demonstrated to be the rate determining step, which proceeds via distinct pathways depending on the stereoelectronic properties of the ligand. A peculiar iridium dihydride route was also uncovered, which could be optimized to accelerate the H 2 heterolysis. We have further computed the thermodynamic and kinetic hydricities of 74 complexes with [Cp*Ir(L)(H)] q ( q = -1, 0, 1 depending on the ligand charge by deprotonation or protonation events) general formula. Hydricity values do not show any noticeable correlation with the thermodynamics (Δ G ) of [Cp*Ir(L)(HCO 2 – )] formation but exhibit linear correlation with the kinetics (Δ G ‡ ) of electrophilic CO 2 attack to iridium hydride species, especially when the charges and local structural effects of the metal hydrides are considered. Insights gained in this study on (i) the factors determining the preferred pathway for the rate-limiting H 2 heterolysis step, (ii) the correlation between thermodynamic hydricity and the kinetics of CO 2 insertion to iridium hydride species, and (iii) mechanistic analysis via combination of free energy profiles and KIE studies for a series of iridium catalysts will provide guidelines for the design of next-generation catalysts for the reversible H 2 storage through CO 2 utilization.

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

Mercury Complexation with Dissolved Organic Matter Released from Thirty-Six Types of Biochar

Previous studies show mercury (Hg) can be effectively removed from solution by biochar, but limited attention was paid on the complexation between Hg and components released from biochars, e.g. dissolved organic matter (DOM). Here, aqueous data from batch-style experiments were modeled using PHREEQC, incorporating thermodynamic constants between Hg and DOM, which was assumed to be composed of thiol, carboxylic, and phenolic functional groups. Modelling results suggest that > 99% Hg complexed with thiol groups in DOM. The modelled concentrations of Hg–DOM complexes from low-T (low-temperature, 300°C) biochars were greater than from high-T (600°C) biochars. The concentrations of Hg–DOM complexes were lower in wood-based than in agricultural residue- and manure-based biochars. Hg–DOM complexes may affect Hg speciation, bioavailability, transport, and methylation processes. This research describes a method to evaluate Hg–DOM interactions, and the results indicate extra caution regarding Hg–DOM complex formation is required in the selection of biochar for Hg remediation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Pb Sorption at the Barite (001)–Water Interface

Sorption of ions at the mineral-water interface is an important factor that determines the fate of toxic metals in the environment. Here, we use barite as a model substrate to understand the interaction of toxic-metal lead (Pb) with ionic crystals. The coverage and location of Pb sorbed to the (001) surface was measured as a function of aqueous Pb concentration ([Pb] aq ) using in situ specular resonant anomalous X-ray reflectivity (RAXR) to determine the sorption capacity and process. In this work, the results show that Pb sorption occurs via incorporation (primarily within the top barite layer similar to 3 angstrom in depth) and adsorption (mostly as an inner-sphere complex at similar to 2 angstrom in height) simultaneously. Both the incorporated and adsorbed Pb coverages increase with increasing [Pb] aq up to [Pb] aq similar to 200 mu M, above which the adsorbed fraction increases more rapidly than the incorporated fraction. This enhanced adsorption has a height distribution that is further extended (>= 15 angstrom from the surface) than that observed in lower [Pb] aq . This change in distribution is interpreted as arising from additional sorption of outer-sphere species or Pb-bearing phases precipitated on the surface. Desorption experiments in Pb-free solutions show that the incorporated fraction is more resistant to removal than the adsorbed fraction, consistent with the speciation-dependent stabilities premised in the classical sorption models.

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

Molecular insights into Yb(III) speciation in sulfate-bearing hydrothermal fluids from X-ray absorption spectra informed by ab initio molecular dynamics

Rare earth elements (REEs) are critical for advanced technologies, yet in hydrothermal aqueous solutions the molecular level details of their interaction with ligands that control their geochemical transport and deposition remain poorly understood. Here, this study elucidates the coordination behavior of Yb 3+ in sulfate-rich hydrothermal fluids using in situ extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulations. By integrating multi-angle EXAFS with AIMD-derived constraints, we precisely resolve Yb 3+ coordination structures and ligand interactions under hydrothermal conditions. At room temperature, Yb 3+ is coordinated by five water molecules and two sulfate ligands (coordination number, CN = 8), forming a distorted square antiprism geometry. Increasing temperature induces progressive dehydration, reducing the hydration shell and favoring stronger sulfate complexation. At 200°C, sulfate ligands reorganize around Yb 3+ , shifting its geometry to a capped octahedron (CN = 7). At 300 °C, sulfate binding dominates, leading to structural reorganization that parallels the onset of sulfate mineral precipitation, consistent with the retrograde solubility of REE sulfates. These findings provide direct molecular-scale evidence that sulfate acts as both a transport and deposition ligand, critically influencing REE mobility in geochemical environments. Our results can also help to refine thermodynamic models of REE speciation in high-temperature hydrothermal fluids and improve our understanding of REE ore formation processes in nature.

AIMD↗
DOE OSTI2020

Arsenic Fate in Peat Controlled by the pH-Dependent Role of Reduced Sulfur

Reduced sulfur (S) has a contrasting role in the fate of arsenic (As) in peatlands. Sulfur bridges provide efficient binding of As to organic carbon (C), but the formation of aqueous As–S species, so-called thioarsenates, leads to a low to no sorption tendency to organic C functional groups. Here, we studied how pH changes the role of reduced S in desorption and retention of presorbed As in model peat. Control desorption experiments without S addition revealed that As was mobilized, predominantly as arsenite, in all treatments with relative mobilization increasing with pH (4.5 < 7.0 < 8.5). Addition of sulfide or polysulfide caused substantial As retention at acidic conditions but significantly enhanced As desorption compared to controls at neutral to alkaline pH. Thioarsenates dominated As speciation at pH 7.0 and 8.5 (maximum, 79%) and remained in solution without (re)sorption to peat. Predominance of arsenite in control experiments and no evidence of surface-bound thioarsenates at pH 7.0 suggest mobilization to proceed via arsenite desorption, reaction with dissolved or surface-bound reduced S, and formation of thioarsenates. Our results suggest that natural or management-related increases in pH or increases in reduced S in near-neutral pH environments can turn organic matter from an As sink into a source.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of nanoconfinement and surface charge on iron adsorption on mesoporous silica

Here, we present a combined molecular dynamics (MD) simulation and X-ray absorption fine structure (XAFS) spectroscopic investigation of aqueous iron adsorption on nanoconfined amorphous silica surfaces. The simulation models examine the effects of pore size, pH (surface charge), iron valency, and counter-ion (chloride or hydroxide). The simulation methods were validated by comparing the coordination environment of adsorbed iron with coordination numbers and bond lengths derived from XAFS. In the MD models, nanoconfinement effects on local iron coordination were investigated by comparing results for unconfined silica surfaces and in confined domains within 2 nm, 4 nm, and 8 nm pores. Experimentally, coordination environments of iron adsorbed onto mesoporous silica with 4 nm and 8 nm pores at pH 7.5 were investigated. The effect of pH in the MD models was included by simulating Fe(II) adsorption onto negatively charged SiO 2 surfaces and Fe(III) adsorption on neutral surfaces. The simulation results show that iron adsorption depends significantly on silica surface charge, as expected based on electrostatic interactions. Adsorption on a negatively charged surface is an order of magnitude greater than on the neutral surface, and simulated surface coverages are consistent with experimental results. Pore size effects from the MD simulations were most notable in the adsorption of Fe(II) at deprotonated surface sites (SiO - ), but adsorption trends varied with concentration and aqueous Fe speciation. The coordination environment of adsorbed iron varied significantly with the type of anion. Considerable ion pairing with hydroxide anions led to the formation of oligomeric surface complexes and aqueous species, resulting in larger iron hydroxide clusters at higher surface loadings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding Ionizing Radiation-Induced Speciation, Chemistry, and Transport in Nuclear Materials

The research to be presented will cover an overview of three projects supported by the United States Department of Energy Office of Science Solar Photochemistry Program: (i) Radiation-Induced Late Actinide Redox Chemistry (2019–2022); (ii) Understanding Radiation-Induced Iodine Speciation, Chemistry, and Transport in High Temperature Molten Salts (2022–2025); and (iii) Radiation-Induced Chemistry of Nuclear Materials (2023–2026). The overarching aim of these distinct projects is to develop a deep fundamental mechanistic understanding of ionizing-radiation-induced processes occurring in solutions over multiple time and length scales. Investigated solutions range from aqueous to organic to molten salt liquids in the presence and absence of ligands and various redox-active ions and molecules, including the actinides. The ultimate goal of these projects is to utilize the new foundational knowledge gained to construct quantitative multiscale computer models with the capacity to predict radiation-induced behavior in complex systems of relevance to nuclear energy technologies and beyond.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2019

Analysis of equilibrium CO 2 solubility and thermodynamic models for aqueous 1‐(2‐hydoxyethyl)‐piperidine solution

Abstract 1‐(2‐Hydoxyethyl)‐piperidine (1‐(2HE)‐PP) is a new tertiary amine with desirable properties and can be potentially used to formulate superior absorbents for CO 2 capture. The equilibrium CO 2 solubility of 1‐(2HE)‐PP solution is measured over temperatures from 298 to 333 K, CO 2 partial pressures from 8.1 to 101.3 kPa and initial amine concentrations from 1 to 5 M. Two thermodynamic models, namely semiempirical model and activity coefficient model are developed for the system. The activity coefficient model shows better estimation solubility with an absolute average relative deviation (AARD) of 7.6%. In the comparison between the two models, a comprehensive analysis is presented. Some suggestions are provided for the similar model development. In addition, the speciation plot of CO 2 loaded 1‐(2HE)‐PP solution is predicted based on the activity coefficient model. The predictive pH values agree well with experimental data with AARD of 1.0%. Finally, the potential of 1‐(2‐HE)PP to be an alternative amine in CO 2 capture is evaluated.

Xiao, Min↗
DOE OSTI2025

The role of Ca-bridged organic matter in an alkaline soil, as revealed by multimodal chemical imaging

Mineral–organic matter (OM) studies have predominantly focused on acidic soils that are abundant in iron (Fe) oxides and aluminum (Al) oxides. We have probed mineral–OM interactions in an alkaline or calcareous soil of the Aridisols class. Unlike the role of Fe and Al, the role of Ca-minerals (particularly calcite), which are ubiquitous in alkaline soils, in OM sequestration is not well understood. Multiple recent model studies with aqueous Ca2+ or synthetic calcite and a suite of OM compounds have shown Ca-OM assemblages to be spatially correlated with calcite at the microscale. To study the chemical state of both Ca and Fe and their competing role in soil organic matter (SOM) stabilization, we performed laboratory characterization using x-ray diffraction, Mössbauer spectroscopy, x-ray photoelectron spectroscopy, scanning electron microscopy, and scanning transmission electron microscopy, alongside synchrotron-based microscale chemical imaging using scanning transmission x-ray microscopy combined with near-edge x-ray absorption fine structure. Ca mineral–organic associations were found to be ubiquitous in this system and are likely critical for understanding SOM stabilization/degradation in alkaline soils. From our findings on mineralogy, speciation, and the nature of Ca-OM bridging, we identified differences in C and Ca chemistry based on the relative location of OM to Ca minerals. The OM near the calcite crystal was enriched in lipid and protein moieties, Ca-OM next to Fe minerals displayed a strong contribution from aromatic compounds, while on the surface of microbes, the carbonate was believed to be of microbial in origin, as also suggested by preliminary works reporting on the formation of amorphous calcite or nano-calcite. In Ca-OM admixed with carbonate, it was difficult to distinguish Ca-associated OM from amorphous calcite or nano-calcite.

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

Impact of dissolved organic matter on mercury and methylmercury sorption to activated carbon in soils: implications for remediation

Activated carbon (AC) amendments have shown promise in reducing inorganic mercury (Hg(II) complexes, “Hg”) and methylmercury (MeHg) risk in contaminated soils. However, the effectiveness of AC in Hg and MeHg immobilization has varied among studies, suggesting that site biogeochemistry might dictate efficacy. In this study, we examined the effect of dissolved organic matter (DOM) on MeHg and Hg sorption to AC. We evaluated the impact of Suwannee River Humic Acid (SRHA) on sorption to AC directly using an isotherm approach and in a soil/AC mixture using slurry microcosms. Aqueous sorption coefficients to AC (log K AC ) for Hg-SRHA and MeHg-SRHA complexes were one to two orders of magnitude lower (Hg-SRHA = 4.53, MeHgSRHA = 4.35) than those for chloride complexes (HgCl 2 = 6.55, MeHgCl = 4.90) and more closely resembled the log K AC of SRHA (3.64). In anoxic, sulfidic soil slurries, the K AC for sulfide species appeared stronger than for chloride or SRHA species for both Hg and MeHg. AC significantly reduced porewater concentrations of both ambient MeHg and a fresh Me 199 Hg spike, and the addition of up to 60 mg L -1 SRHA did not reduce sorption to AC. The AC also reduced ambient Hg and 201 Hg porewater concentrations, but as SRHA concentration increased, the magnitude of solid phase sorption decreased. Speciation modeling revealed that SRHA may have impacted Hg distribution to the solid phase by reducing HgS precipitation. In conclusion, this study highlights the need for site-specific evaluation of AC efficacy and the value in developing biogeochemical models of AC performance for Hg control.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Acidities of MgO surface sites: implications for the formation mechanism of Mg(OH) 2

The hydroxylation of periclase (MgO) to brucite (Mg(OH) 2 ) is thought to be an important intermediate step when using MgO to capture CO 2 from the atmosphere. However, the mechanism of hydroxylation of MgO to form Mg(OH) 2 is poorly understood. In this work, we used atomic-scale density functional tight binding simulations coupled with the metadynamics rare event method to analyze the surface chemistry of MgO and the acid dissociation equilibrium constants (pK a ) of its surface sites. The method and parameters were validated by calculating the pK a for hydroxylation of the first shell water bound to aqueous Mg 2+ ion. The pK a value derived using a probabilistic method was 12.3, which is in fair agreement with the accepted value of 11.4, with the difference between them equal to a ∼5 kJ mol −1 error in the calculations. We then extended these pK a calculations to probe the hydroxylation reactions of the surface sites of the MgO(100)–water interface, arriving at pK a s of 5.4 to deprotonate terminal water molecules bound to the surface magnesium sites (η-OH 2 or 〉MgOH 2 ), and 13.9 to deprotonate hydroxylated bridging oxygen sites (μ 5 -oxo or 〉O). Hydroxide (OH − ) adsorption on the surface was also probed and found to be less thermodynamically favorable than deprotonation of the terminal water molecule. The plausibility of the computed pK a s was verified using an activity-based speciation model and compared to pH measurements of water equilibrated with MgO nanoparticles and single crystals. The model predicted a solution pH of 7.1 when surface sites buffered and the pH of 12.0 when MgO dissolution dominated. These are close to the experimental initial solution pHs of 7–7.5 and the long term pHs of ∼10.5. The similarity suggests that the calculated pK a values from the DFTB+/metadynamics simulations are plausible and that these methods can be a useful tool to probe reaction mechanisms involving covalent bonds.

Adapa, Sai Krishna Reddy [Oak Ridge National Labor↗
DOE OSTI2025

Unravelling the radiation-induced redox chemistry of plutonium ions in aqueous solution

Plutonium plays a critical role in nuclear fuel cycle technologies, but our understanding of its fundamental radiation-induced redox chemistry is limited. Changes in oxidation states affect the speciation and transport of plutonium ions in solution. For example, solvent extraction techniques used to separate and recover plutonium from used nuclear fuel rely on the selective formation, maintenance, and complexation of specific plutonium oxidation states. However, radiolytically generated radicals, ions, and molecules can drive the oxidation state distribution of plutonium ions far from equilibrium, ultimately changing the physical and chemical properties of the bulk system. These radiation-induced processes are inevitable due to the ionizing radiation fields generated by the radioactive decay of plutonium and its daughter nuclides. Therefore, mechanistically understanding how plutonium's various oxidation states respond to ionizing radiation is essential for predicting its behavior in solution. Here, we present significant advances in our understanding of radiation-induced plutonium redox chemistry by using time-resolved (electron pulse) and dose accumulation (alpha and gamma) irradiation techniques, along with quantitative multiscale modeling methods.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2019

Environmental Mercury Chemistry – In Silico

Mercury (Hg) is a global environmental contaminant. Major anthropogenic sources of Hg emission include gold mining and the burning of fossil fuels. Once deposited in aquatic environments, Hg can undergo redox reactions, form complexes with ligands, and adsorb onto particles. It can also be methylated by microorganisms. Mercury, especially its methylated form methylmercury, can be taken up by organisms, where it bioaccumulates and biomagnifies in the food chain, leading to detrimental effects on ecosystem and human health. In support of the recently enforced Minamata Convention on Mercury, a legally binding international convention aimed at reducing the anthropogenic emission of—and human exposure to—Hg, its global biogeochemical cycle must be understood. Thus, a detailed understanding of the molecular-level interactions of Hg is crucial.The ongoing rapid development of hardware and methods has brought computational chemistry to a point that it can usefully inform environmental science. This is particularly true for Hg, which is difficult to handle experimentally due to its ultratrace concentrations in the environment and its toxicity. The current account provides a synopsis of the application of computational chemistry to filling several major knowledge gaps in environmental Hg chemistry that have not been adequately addressed experimentally.Environmental Hg chemistry requires defining the factors that determine the relative affinities of different ligands for Hg species, as they are critical for understanding its speciation, transformation and bioaccumulation in the environment. Formation constants and the nature of bonding have been determined computationally for environmentally relevant Hg(II) complexes such as chlorides, hydroxides, sulfides and selenides, in various physical phases. Quantum chemistry has been used to determine the driving forces behind the speciation of Hg with hydrochalcogenide and halide ligands. Of particular importance is the detailed characterization of solvation effects. Indeed, the aqueous phase reverses trends in affinities found computationally in the gas phase. Computation has also been used to investigate complexes of methylmercury with (seleno)amino acids, providing a molecular-level understanding of the toxicological antagonism between Hg and selenium (Se). Furthermore, evidence is emerging that ice surfaces play an important role in Hg transport and transformation in polar and alpine regions. Therefore, the diffusion of Hg and its ions through an idealized ice surface has been characterized.Microorganisms are major players in environmental mercury cycling. Some methylate inorganic Hg species, whereas others demethylate methylmercury. Quantum chemistry has been used to investigate catalytic mechanisms of enzymatic Hg methylation and demethylation. The complex interplay between the myriad chemical reactions and transport properties both in and outside microbial cells determines net biogeochemical cycling. In conclusion, prospects for scaling up molecular work to obtain a mechanistic understanding of Hg cycling with comprehensive multiscale biogeochemical modeling are also discussed.

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

Genome-Resolved Metagenomics and Detailed Geochemical Speciation Analyses Yield New Insights into Microbial Mercury Cycling in Geothermal Springs

Geothermal systems emit substantial amounts of aqueous, gaseous, and methylated mercury, but little is known about microbial influences on mercury speciation. Here, we report results from genome-resolved metagenomics and mercury speciation analysis of acidic warm springs in the Ngawha Geothermal Field (< 55°C, pH < 4.5), Northland Region, Aotearoa New Zealand. Our aim was to identify the microorganisms genetically equipped for mercury methylation, demethylation, or Hg(II) reduction to volatile Hg(0) in these springs. Dissolved total and methylated mercury concentrations in two adjacent springs with different mercury speciation ranked among the highest reported from natural sources (250 to 16,000 ng liter -1 and 0.5 to 13.9 ng liter -1 , respectively). Total solid mercury concentrations in spring sediments ranged from 1,274 to 7,000 μg g -1 . In the context of such ultrahigh mercury levels, the geothermal microbiome was unexpectedly diverse and dominated by acidophilic and mesophilic sulfur- and iron-cycling bacteria, mercury- and arsenic-resistant bacteria, and thermophilic and acidophilic archaea. By integrating microbiome structure and metagenomic potential with geochemical constraints, we constructed a conceptual model for biogeochemical mercury cycling in geothermal springs. The model includes abiotic and biotic controls on mercury speciation and illustrates how geothermal mercury cycling may couple to microbial community dynamics and sulfur and iron biogeochemistry.

59 BASIC BIOLOGICAL SCIENCES↗
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↗