Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Application of community data to surface complexation modeling framework development: Iron oxide protolysis

This study presents a comprehensive community data-driven surface complexation modeling framework for simulating potentiometric titration of mineral surfaces. Compiled community data for ferrihydrite, goethite, hematite, and magnetite are fit to produce representative protolysis constants that can reproduce potentiometric titration data collected from multiple literature sources. Using this framework, the impact of surface complexation model type and surface site density (SSD) on the fit quality and protolysis constants can be readily evaluated. For example, the non-electrostatic model yielded a poor data fit compared to diffuse double layer model and constant capacitance models due to the absence of known surface charge effects. Regardless of the choice of iron oxide mineral, pK a1 decreased with increasing SSD while the opposite tendency was observed for pK a2 . This newly developed framework demonstrates a method to reconcile community data-wide potentiometric titration data using Findable, Accessible, Interoperable, Reusable data principles to produce mineral protolysis constants that improve robustness of surface complexation models for applications in metal sorption and reactive transport modeling. The framework is readily expandable (as community data increase) and extensible (as the number of minerals increase). The framework provides a path forward for developing self-consistent, comprehensive, and updateable surface complexation databases for surface complexation and reactive transport modeling.

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

Mechanisms of Zn removal from water by amorphous geopolymer: Molecular-level insights from X-ray absorption spectroscopy, isotope fractionation, and surface complexation modeling

Porous geopolymers have attracted widespread attention as promising heavy metal adsorbents that can be synthesized from aluminosilicate solid wastes. However, the precise microstructural evidence for adsorbed heavy metals on geopolymers remains unclear due to the insensitivity of conventional characterization techniques on minerals with amorphous structure and surface disorder. Batch adsorption and column experiments coupled with X-ray absorption spectroscopy (XAS), Zn stable isotope, and surface complexation model (SCM) were employed to reveal the Zn removal mechanisms with coal fly ash porous geopolymer (CFAPG) at a molecular scale. The macroscopic kinetic and isothermal adsorption of Zn on CFAPG were well described by the pseudo-second-order model and Bi_Langmuir equation, respectively, indicating the presence of abundant heterogeneous active sites on the CFAPG surface. Further, two types of active sites on the CFAPG surface were identified by XAS coupled with Zn isotopes in batch experiments at pH <= 6.0: one is pH-dependent and associated with tetrahedral zinc coordi-nation, and the other is pH-insensitive and associated with octahedral zinc coordination; these sites were confirmed by the bidentate SCM as the variable charge site (surface complexation, >S-OH) and the permanent negative charge site (cation exchange, >X - ), respectively. Furthermore, the important contribution of surface co -precipitation besides surface complexation and cation exchange to the Zn adsorption on CFAPG was identified by XAS coupled with SCM in a flow-through column experiment at pH >6.0. These investigations provide a systemic understanding of the Zn adsorption mechanisms on CFAPG and an SCM reference for the application and prediction of geopolymers in heavy metal-contaminated water remediations.

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

Machine learning surrogates for surface complexation model of uranium sorption to oxides

Abstract The safety assessments of the geological storage of spent nuclear fuel require understanding the underground radionuclide mobility in case of a leakage from multi-barrier canisters. Uranium, the most common radionuclide in non-reprocessed spent nuclear fuels, is immobile in reduced form (U(IV) and highly mobile in an oxidized state (U(VI)). The latter form is considered one of the most dangerous environmental threats in the safety assessments of spent nuclear fuel repositories. The sorption of uranium to mineral surfaces surrounding the repository limits their mobility. We quantify uranium sorption using surface complexation models (SCMs). Unfortunately, numerical SCM solvers often encounter convergence problems due to the complex nature of convoluted equations and correlations between model parameters. This study explored two machine learning surrogates for the 2-pK Triple Layer Model of uranium retention by oxide surfaces if released as U(IV) in the oxidizing conditions: random forest regressor and deep neural networks. Our surrogate models, particularly DNN, accurately reproduce SCM model predictions at a fraction of the computational cost without any convergence issues. The safety assessment of spent fuel repositories, specifically the migration of leaked radioactive waste, will benefit from having ultrafast AI/ML surrogates for the computationally expensive sorption models that can be easily incorporated into larger-scale contaminant migration models. One such model is presented here.

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

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗
DOE OSTI2024

Adsorption of Neodymium, Dysprosium, and Ytterbium to Goethite under Varying Aqueous Chemistry Conditions

The adsorption of rare earth elements (REEs) to iron oxides can regulate the mobility of REEs in the environment and is heavily influenced by water chemistry. This study utilized batch experiments to examine the adsorption of Nd, Dy, and Yb to goethite under varying pH, electrolyte (type and concentration), and concentrations of dissolved inorganic carbon and citrate. REE adsorption was strongly influenced by pH, with an increase from essentially no adsorption at pH 3.0 to nearly complete adsorption at pH 6.5 and higher. Citrate enhanced the adsorption of REEs at low pH (<5.0), likely by forming goethite-REE-citrate ternary surface complexes. However, citrate inhibited the adsorption of REEs at higher pH (>5.0) by forming aqueous REE-citrate complexes. Ionic strength had a small influence on REE adsorption, and the presence of dissolved inorganic carbon had no discernible effect. Equilibrium adsorption was interpreted with a triple-layer surface complexation model (SCM). The selection of surface complexation reactions was guided by extended X-ray absorption fine structure spectra. An SCM with a single bidentate inner-sphere surface complexation reaction for Nd and two inner-sphere surface complexation reactions (one monodentate and one bidentate reaction) for Dy and Yb effectively simulated adsorption across a broad range of conditions in the absence of citrate. Accounting for the effects of citrate on REE adsorption required the addition of up to two ternary REE-citrate-goethite surface complexes. The SCM can enable predictions of REE transport in subsurface environments that have goethite as an important adsorbent mineral. Furthermore, this predictive capability could contribute to identifying potential REE sources and facilitating efficient extraction of REEs.

58 GEOSCIENCES↗
DOE OSTI2023

LLNL Automized Surface Titration Model

The LLNL Automized Surface Titration Model (L-ASTM) is a community data-driven surface complexation modeling workflow for simulating potentiometric titration of mineral surfaces. The model accepts raw experimental potentiometric titration data formatted in a findable, accessible, interoperable, and reusable (FAIR) structure. The workflow was coded in Python and coupled to PHREEQC for surface complexation modeling and PEST for data fitting and parameter estimation.

Solchan, Han↗
DOE OSTI2024

M3SF-24LL010302052-Comprehensive Analysis of Radionuclide Interaction with Hydrothermally Altered Repository Materials

This progress report (Level 3 Milestone Number M3SF-24LL010302052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline Work Package Number SF-24LL01030205. The research is focused on actinide and radionuclide sequestration in hydrothermally altered repository materials. In FY24, we completed a rigorous analysis of Se sorption to iron oxide phases using our L-SCIE sorption database. This effort explicitly accounts for surface titration behavior of oxide surfaces that was the subject of a recently published manuscript (Han et al., 2023). With this in mind, our Se sorption analysis now yields a more robust workflow for developing self-consistent surface complexation modeling approaches that can be adapted to specific SCM conceptual and numerical approaches (i.e. non-electrostatic, diffuse layer, triple layer models). In late FY24, we plan to publish the results of our comprehensive surface complexation modeling of Se(IV) and Se(VI) sorption to iron oxide mineral phases and provide a path forward to developing robust radionuclide sorption models for use in performance assessment. In FY24, we also submitted a manuscript summarizing our approach to integrating radionuclide sorption and coprecipitation phenomena and evaluation of radionuclide partitioning values across a range of radionuclides relevant to performance assessment. We demonstrated our approach in detail using Se sorption and coprecipitation with iron oxide minerals as a test case. This manuscript was recently accepted for publication in Applied Geochemistry (Balboni et al., Accepted). We also continued experiments to identify radionuclide interaction with hydrothermally altered crystalline repository and backfill materials. Recent research performed at Los Alamos National Laboratory (LANL) and Sandia National Laboratory (SNL) has provided key insights regarding the hydrothermal alteration behavior of bentonite backfill in the presence of repository materials (steel, concrete, etc.). We are now examining how mineral alteration affects retardation behavior of a suite of radionuclides of interest to repository performance assessment. These experiments also allow us to test the predictive ability of our component additivity approach to surface complexation and ion exchange. Our guiding hypothesis is that a robust surface complexation/ion exchange model and associated database, developed using our L-SCIE approach, can effectively predict changes in radionuclide sorption behavior resulting from the hydrothermal alteration of mineralogy in a repository near field. A short update of results to date is presented below.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Management and International Sorption Model Collaboration (M4SF-23LL010302062-NEA-TDB)

This progress report (Level 4 Milestone Number M4SF-23LL010302062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline International Collaborations Activity Number SF-23LL01030206. The activity is focused on our long-term commitment of engaging our partners in international nuclear waste repository research. This includes participation in the Nuclear Energy Agency Thermochemical Database (NEA-TDB) Project and development of methodologies for integrating US and international thermodynamic databases for use in SFWST Generic Disposal System Assessment (GDSA) efforts. A continuing focus for FY23 efforts has been to support the US participation in the NEA-TDB effort (Mavrik Zavarin replaced Cindy Atkins-Duffin on the NEA-TDB Management Board (MB) and Executive Group (EG)) and developing mechanisms for integration of NEA-TDB thermochemical data with LLNL’s SUPCRTNE thermodynamic database that supports the SFWST GDSA activities. This effort is coordinated with the Argillite work package SUPCRTNE database development efforts. The goal is to provide a downloadable database that will be hosted on a LLNL website which integrates NEA-TDB data into the LLNL SUPCRTNE database where appropriate. As part of our international activities, we continue our effort to integrate international sorption databases into L-SCIE (Zavarin et al., 2022b). We presented opportunities to include sorption in the next phase of NEA-TDB efforts at the April 2023 EG meeting in Paris. FY23 efforts focused on ensuring interoperable database development across multiple international database development activities. The overall goal is to produce an open source database that can be shared and integrated with multiple nuclear waste programs internationally and harness modern data science workflows and algorithms to incorporate these new approaches into reactive transport and performance assessment models. In collaboration with our Helmholtz Zentrum Dresden Rossendorf partners, we recently demonstrated the power of FAIR open source databases by fitting iron oxide (hydrous ferric oxide, goethite, hematite, and magnetite) protolysis constants to all available L-SCIE data. The results were submitted as a manuscript to J. Colloid Interface Science. This work will inform future metal sorption studies on a variety of iron oxides in order to discern the most appropriate acidity constants and surface complexation modeling constructs to account for pH-dependent mineral surface charge behavior. This work also explored automated surface complexation model development workflows in order to generate higher throughput model input files for a more facile incorporation into GDSA activities.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Adsorption of neodymium, dysprosium, and ytterbium to goethite under varying aqueous chemistry conditions [dataset]

The adsorption of rare earth elements (REEs) to iron oxides can regulate the mobility of REEs in the environment and is heavily influenced by water chemistry. This study utilized batch experiments to examine the adsorption of Nd, Dy, and Yb to goethite under varying pH, electrolyte (type and concentration), and concentrations of dissolved inorganic carbon and citrate. REE adsorption was strongly influenced by pH, with an increase from essentially no adsorption at pH 3.0 to nearly complete adsorption at pH 6.5 and higher. Citrate enhanced the adsorption of REEs at low pH (< 5.0), likely by forming goethite-REE-citrate ternary surface complexes. However, citrate inhibited the adsorption of REEs at higher pH (> 5.0) by forming aqueous REE-citrate complexes. Ionic strength had a small influence on REE adsorption, and the presence of dissolved inorganic carbon had no discernible effect. Equilibrium adsorption was interpreted with a triple layer surface complexation model (SCM). The selection of surface complexation reactions was guided by extended X-ray absorption fine structure spectra. A SCM with a single bidentate inner-sphere surface complexation reaction for Nd and two inner-sphere surface complexation reactions (one monodentate and one bidentate reaction) for Dy and Yb effectively simulated adsorption across a broad range of conditions in the absence of citrate. Accounting for the effects of citrate on REE adsorption required the addition of up to two ternary REE-citrate-goethite surface complexes. The SCM can enable predictions of REE transport in subsurface environments that have goethite as an important adsorbent mineral. This predictive capability could contribute to identifying potential REE sources and facilitating efficient extraction of REEs.

58 GEOSCIENCES↗
DOE OSTI2023

A chemistry-informed hybrid machine learning approach to predict metal adsorption onto mineral surfaces

Historically, surface complexation model (SCM) constants and distribution coefficients (K d ) have been employed to quantify mineral-based retardation effects controlling the fate of metals in subsurface geologic systems. Our recent SCM development workflow, based on the Lawrence Livermore National Laboratory Surface Complexation/Ion Exchange (L-SCIE) database, illustrated a community FAIR data approach to SCM development by predicting uranium(VI)-quartz adsorption for a large number of literature-mined data. Here, we present an alternative hybrid machine learning (ML) approach that shows promise in achieving equivalent high-quality predictions compared to traditional surface complexation models. At its core, the hybrid random forest (RF) ML approach is motivated by the proliferation of incongruent SCMs in the literature that limit their applicability in reactive transport models. Our hybrid ML approach implements PHREEQC-based aqueous speciation calculations; values from these simulations are automatically used as input features for a random forest (RF) algorithm to quantify adsorption and avoid SCM modeling constraints entirely. Named the LLNL Speciation Updated Random Forest (L-SURF) model, this hybrid approach is shown to have applicability to U(VI) sorption cases driven by both ion-exchange and surface complexation, as is shown for quartz and montmorillonite cases. The approach can be applied to reactive transport modeling and may provide an alternative to the costly development of self-consistent SCM reaction databases.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Prediction of alkaline earth metal ion adsorption on goethite for various background electrolytes with the CD-MUSIC model

As water scarcity drives the use of more saline water sources, contaminant fate and transport models must capture the impact of high concentrations of alkaline earth metal ions (AEMs) and background electrolytes in these more complex waters. By utilizing macroscopic adsorption data from various electrolyte systems, a Charge Distribution – Multisite Complexation (CD-MUSIC) model, capable of incorporating electrolyte adsorption, was able to accurately simulate the adsorption behavior of alkaline earth metal ions onto goethite. The modeling effort was guided by previous spectroscopic and surface complexation modeling of alkaline earth metal adsorption and built on previous CD-MUSIC modeling that accounted for changes in crystal face contributions to the surface site density as a function of specific surface area. Further, the model was constrained to consider only two dominant surface complex species for each metal ion adsorption reaction. These two species were selected from 44 possible species through objective curve fitting of single-solute macroscopic adsorption data. While most of the alkaline earth metal surface complexes formed outer-sphere complexes at the goethite surface, an inner-sphere species was utilized for Mg 2+ . With the surface complex species and equilibrium constants obtained from this study, the calibrated model successfully predicted alkaline earth metal ion adsorption over a wide range of solution and surface conditions; the model predictions encompassed a wide range of pH (5–11), solute/solid ratio (1.37 × 10 -5 – 8.33 × 10 -4 mol -solute /g -solid ), ionic strengths (0.01 M – 0.7 M), and background electrolytes (Na + , Cs + , Rb + , Cl - , and NO 3 - ) using the same crystal face contribution methodology for site density, capacitance values, and surface acidity constants adopted for proton and cadmium adsorption in previous work (Han and Katz, 2019). Model simulations for a range of background water chemistries demonstrated the potential for Mg 2+ to reduce Cd 2+ adsorption to goethite in model seawater and oil- and gas-produced waters.

42 ENGINEERING↗
DOE OSTI2023

Surface Complexation/Ion Exchange Hybrid Model for Radionuclide Sorption to Clay Minerals (M4SF-23LL010301062)

This progress report (Level 4 Milestone Number M4SF-23LL010301062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite International Collaborations Activity Number SF-23LL01030106. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. The focus of this milestone is the establishment of international collaborations for surface complexation modeling and the associated impacts of unlocking larger, community-based datasets. More specifically, we are developing a database framework for Spent Fuel and Waste and Science Technology (SFWST) that is aligned with the Helmholtz Zentrum Dresden Rossendorf (HZDR) sorption database development group in support of the database needs of the SFWST program. In our FY22 effort, we described a detailed analysis of U(VI) sorption to quartz through both traditional surface complexation modeling and through a hybrid ML framework. In FY23, effort was placed on publication of these results and expansion of the LLNL surface complexation and ion exchange database (L-SCIE) in order to assess mineral-based radionuclide retardation under a wider variety of geochemical conditions (e.g., ionic strength, varying electrolyte compositions). Efforts were initiated to expand L-SCIE to include radionuclide surface complexation and ion exchange to clays that are relevant to subsurface geochemical processes occurring at nuclear waste repositories. In particular, a large source of sorption data for clays resides at the Paul Scherrer Institute (PSI) (work primarily by Bradbury and Baeyens) and we initiated discussions on how to retrieve those data and apply FAIR principles to those datasets. In addition to L-SCIE development, two hybrid models that incorporate AI/ML were investigated and compared to discern the most promising approaches for accurate and precise estimations of radionuclide retardation. Key considerations for future model development include (1) the ability to reduce computational burden on determining retardation coefficients for PA and (2) the ability to quantify and predict radionuclide-mineral partitioning at a more efficient, rapid pace due to automated workflows. Upon the careful consideration of the most effective modeling approaches, we are identifying ways to implement these approaches into PA. Ultimately, the data science-based workflows will provide a major incentive for other institutions to adopt a FAIR-formatted, interoperable database. LLNL will play a key role in disseminating sorption data and acting as good data stewards by updating the database in a consistent format and assessing the quality of the newly assimilated data in an organized fashion. To this end, all data and workflows are open access and made available on the LLNL Seaborg research website (https://seaborg.llnl.gov/resources/geochemical-databases-modeling-codes).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Adsorption of REEs to Kaolinite via Ion Exchange and Surface Complexation as a Function of Water Chemistry

Rare earth elements (REEs) are critical components of modern technology behind renewable energy, transportation, and electronics but have a limited current supply. A substantial portion of global REE production relies on ion adsorption deposits. A high abundance of kaolinite in REE enrichment zones within these deposits suggests that kaolinite controls the subsurface migration of REEs. This study aimed to improve the current understanding of REE binding to kaolinite under varying water chemistry conditions. We conducted batch experiments with kaolinite (KGa-2) and three REEs (Nd, Dy, and Yb) at varying pH, electrolyte concentration, dissolved inorganic carbon (DIC), low molecular weight organic acids (citric and oxalic acids), and total REE concentration conditions. Increasing electrolyte concentration inhibits REE adsorption at pH < 7, suggesting that ion exchange contributes to adsorption at these pH values. DIC affects adsorption above pH 7–8 by forming strong aqueous complexes with heavy REEs. Citric acid decreases REE adsorption via aqueous complexation of REEs at pH > 5 but does not affect adsorption at pH < 5. The surface complexation model captures the main adsorption trends with two mechanisms: ion exchange on basal planes at pH < ∼6 and inner-sphere surface complexation to edge sites at pH > ∼6. Equilibrium constants for surface complexation increase in the order of Yb > Dy > Nd, indicating a higher strength of adsorption for heavy REEs. This study demonstrates how water chemistry conditions control the adsorption mechanisms that may determine the mobility of REEs in subsurface environments rich in kaolinite.

58 GEOSCIENCES↗
DOE OSTI2025

Adsorption of Rare Earth Elements to Kaolinite [dataset]

Rare earth elements (REEs) are critical components of modern technology behind renewable energy, transportation, and electronics but have a limited current supply. A substantial portion of global REE production relies on ion adsorption deposits. A high abundance of kaolinite in REE enrichment zones within these deposits suggests that kaolinite controls the subsurface migration of REEs. This study aimed to improve the current understanding of REE binding to kaolinite under varying water chemistry conditions. We conducted batch experiments with kaolinite (KGa-2) and three REEs (Nd, Dy, and Yb) at varying pH, electrolyte concentration, dissolved inorganic carbon (DIC), low molecular weight organic acids (citric and oxalic acids), and total REE concentration conditions. Increasing electrolyte concentration inhibits REE adsorption at pH 7, suggesting that ion exchange contributes to adsorption at these pH values. DIC affects adsorption above pH 7−8 by forming strong aqueous complexes with heavy REEs. Citric acid decreases REE adsorption via aqueous complexation of REEs at pH 5 but does not affect adsorption at pH 5. The surface complexation model captures the main adsorption trends with two mechanisms: ion exchange on basal planes at pH ∼6 and inner-sphere surface complexation to edge sites at pH ∼6. Equilibrium constants for surface complexation increase in the order of Yb Dy Nd, indicating a higher strength of adsorption for heavy REEs. This study demonstrates how water chemistry conditions control the adsorption mechanisms that may determine the mobility of REEs in subsurface environments rich in kaolinite.

58 GEOSCIENCES↗