Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Computational Thermochemistry of Jet Fuels and Rocket Propellants

The design of new high-energy density molecules as candidates for jet and rocket fuels is an important goal of modern chemical thermodynamics. The NASA Glenn Research Center is home to a database of thermodynamic data for over 2000 compounds related to this goal, in the form of least-squares fits of heat capacities, enthalpies, and entropies as functions of temperature over the range of 300 - 6000 K. The chemical equilibrium with applications (CEA) program written and maintained by researchers at NASA Glenn over the last fifty years, makes use of this database for modeling the performance of potential rocket propellants. During its long history, the NASA Glenn database has been developed based on experimental results and data published in the scientific literature such as the standard JANAF tables. The recent development of efficient computational techniques based on quantum chemical methods provides an alternative source of information for expansion of such databases. For example, it is now possible to model dissociation or combustion reactions of small molecules to high accuracy using techniques such as coupled cluster theory or density functional theory. Unfortunately, the current applicability of reliable computational models is limited to relatively small molecules containing only around a dozen (non-hydrogen) atoms. We propose to extend the applicability of coupled cluster theory- often referred to as the 'gold standard' of quantum chemical methods- to molecules containing 30-50 non-hydrogen atoms. The centerpiece of this work is the concept of local correlation, in which the description of the electron interactions- known as electron correlation effects- are reduced to only their most important localized components. Such an advance has the potential to greatly expand the current reach of computational thermochemistry and thus to have a significant impact on the theoretical study of jet and rocket propellants.

Crawford, T. Daniel↗
DOE OSTI2023

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low-rank approximation of the dielectric matrix

We present an efficient preconditioning technique for accelerating the fixed-point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low-rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed-point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous self-consistent field iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from ~100 to 1100 atoms (~500–20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized d orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20–30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed-point iteration in some cases. Lastly, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4-fold reduction in computational cost for the total ground-state calculation compared to other preconditioners.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A Rate Theory Model of Radiation-Induced Swelling in an Austenitic Stainless Steel

Many rate theory models of cavity (void) swelling have been published over the past 50 years, all having the same, or similar, structures. A rigorous validation of the models has not been possible because of the dearth of information concerning the microstructures that correspond with the swelling data. Whereas the lack of microstructure information is still an issue for historical swelling data, in the past 10–20 years data have been published on the evolution of the microstructure (point defect yields from collision cascades, cavity number densities, and dislocation densities/yield strengths) allowing certain gaps in information to be filled when considering historic swelling data. With reasonable estimates of key microstructure parameters, a standard rate theory model can be applied, and the model parameter space explored, in connection with historical swelling data. By using published data on: (i) yield strength as a function of dose and temperature (to establish an empirical expression for dislocation density evolution); (ii) cavity number densities as a function of temperature; and (iii) freely migrating defect (FMD) production as a function of primary knock-on atom (PKA) spectrum, the necessary parameter and microstructure inputs that were previously unknown can be used in model development. This paper describes a rate-theory model for void swelling of 316 stainless steel irradiated in the EBR-2 reactor as a function of irradiation temperature and neutron dose.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-Performance, Long-Life, Rechargeable Li–CO 2 Batteries based on a 3D Holey Graphene Cathode Implanted with Single Iron Atoms

A highly efficient cathode catalyst for rechargeable Li-CO 2 batteries is successfully synthesized by implanting single iron atoms into 3D porous carbon architectures, consisting of interconnected N,S-codoped holey graphene (HG) sheets. The unique porous 3D hierarchical architecture of the catalyst with a large surface area and sufficient space within the interconnected HG framework can not only facilitate electron transport and CO 2 /Li+ diffusion, but also allow for a high uptake of Li 2 CO 3 to ensure a high capacity. Consequently, the resultant rechargeable Li-CO 2 batteries exhibit a low potential gap of approximate to 1.17 V at 100 mA g -1 and can be repeatedly charged and discharged for over 200 cycles with a cut-off capacity of 1000 mAh g -1 at a high current density of 1 A g -1 . Density functional theory calculations are performed and the observed appealing catalytic performance is correlated with the hierarchical structure of the carbon catalyst. This work provides an effective approach to the development of highly efficient cathode catalysts for metal-CO 2 batteries and beyond.

(2) batteries↗
DOE OSTI2025

Achieving Multimodal and Multicolor Luminescence in LaAlO 3 :Pr 3+ , Gd 3+ via Trap Engineering and Energy Transfer

Achieving multimodal luminescence within a single phosphor is vital for multifunctional applications but remains challenging due to complex color tuning and trap engineering. In this study, we report Pr 3+ and Gd 3+ co‐doped LaAlO 3 (LAO:PG) phosphors, designed through careful modulation of multilevel traps and Pr 3+ → Gd 3+ energy transfer dynamics. These materials exhibit diverse luminescence modes, including down‐conversion luminescence (DCL), up‐conversion luminescence (UCL), persistent luminescence (PersL), optically stimulated luminescence (OSL), and thermally stimulated luminescence (TSL) across a wide spectral range. Unlike previously studied Pr 3+ ‐doped LAO, the co‐doped LAO:PG shows DCL in both UV‐visible and NIR regions and displays ultraviolet‐C UCL under visible excitation. Notably, we observe, for the first time, PersL lasting several minutes in these phosphors—an improvement over the non‐PersL behavior of Pr 3+ ‐only doped LAO. Additionally, the LAO:PG phosphors exhibit strong OSL response. TSL analysis reveals five distinct trap levels linked to these properties. Density functional theory calculations further correlate intrinsic defects to these traps, supporting a proposed mechanism for the observed multimodal luminescence. These findings highlight LAO:PG as a promising platform for developing advanced phosphors with integrated luminescence modes, paving the way for future applications in data storage, phototherapy, and anti‐counterfeiting technologies.

Chemistry↗
DOE OSTI2022

A Low-Voltage Layered Na 2 TiGeO 5 Anode for Lithium-Ion Battery

Titanium-based anode materials have achieved much progress with the wide studies in lithium-ion batteries. However, these known materials usually possess high discharge voltage platforms and limited energy densities. In this work, a titanium-based oxide of Na 2 TiGeO 5 with layered structure, two-dimensional lamellar frame and exposed highly active (001) facet, exhibiting good electrochemical performance in terms of high capacity (410 mAh g –1 with a current density of 50 mA g –1 ), excellent rate capability and cycling stability with no obvious capacity attenuation after 4000 cycles, is reported. The appropriate discharge voltage plateau at around 0.2 V endows the Na 2 TiGeO 5 anode material high security compared with graphite and high energy density compared with spinel Li 4 Ti 5 O 12 . Combining the electrochemical tests and the density functional theory calculations, the Li + storage mechanism of Na 2 TiGeO 5 is elucidated and the conversion reaction process is revealed. More importantly, this study provides a way to develop low-voltage and high-capacity titanium-based anode materials for efficient energy storage.

25 ENERGY STORAGE↗
DOE OSTI2023

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗
DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

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

Distinct behavior of localized and delocalized carriers in anatase TiO 2 (001) during reaction with O 2

Two-dimensional (2D) metallic states induced by oxygen vacancies $(V_{O}\text{s})$ at oxide surfaces and interfaces provide opportunities for the development of advanced applications, but the ability to control the behavior of these states is still limited. As such, we used angle resolved photoelectron spectroscopy combined with density-functional theory (DFT) to study the reactivity of $V_O$-induced states at the (001) surface of anatase TiO 2 , where both 2D metallic and deeper lying in-gap states (IGs) are observed. The 2D and IG states exhibit remarkably different evolutions when the surface is exposed to molecular O 2 : while IGs are almost completely quenched, the metallic states are only weakly affected. DFT calculations indeed show that the IGs originate from surface $V_{O}\text{s}$ and remain localized at the surface, where they can promptly react with O 2 . In contrast, the metallic states originate from subsurface vacancies whose migration to the surface for recombination with O 2 is kinetically hindered on anatase TiO 2 (001), thus making them much less sensitive to oxygen dosing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

FY20 Proxy App Suite Release: Report for ECP Proxy App Project Milestone ADCD-504-10

Version 4.0 of the ECP Proxy App Suite is practically unchanged from the previous release. The current set of proxies has proven useful for many aspects of benchmarking and co-design and we see little reason to alter the suite. Although there have been few changes to the ECP suite, the team has been hard at work in other areas. In the area of Machine Learning (ML) we have now created a separate proxy suite dedicated to this scientific applications of ML. The suite includes: miniGAN (Generative Adversarial Networks), miniRL (Reinforcement Learning), CRADL (inline inference), Cosmoflow-Benchmark (Convolutional Neural Network), and MLPerf-DeepCam (Climate Segmentation Benchmark). Section 3 contains more information about these proxies as well as the principles that are guiding the development of the suite. We have surveyed available proxies for several application domains including Computational Fluid Dynamics, Quantum Chemistry, Quantum Computing Simulation, Molecular Dynamics, Monte Carlo Transport, and Density Functional Theory to identify gaps in proxy coverage. Several new proxy apps are either already available or will be released soon to fill these gaps. Section 4 provides full details. Finally, section 5 reports on our continued collaboration with the ECP Continuous Integration (CI) effort to use proxy apps to help identify problems and roadblocks to cross-lab CI. We have also assisted the El Capitan Center of Excellence (COE) to stand up a CI system that will be used to test releases of HPE and AMD software stacks. We hope that the COE effort can serve as a model for ECP by showing how the Proxy App Team can work with AD and ST teams identify critical features, kernels, patterns, etc. and incorporate them into a CI system that will help ensure that Frontier and Aurora will provide those needed capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Density functional theory (DFT) study of UF 6 hydrolysis: reaction pathways, spectroscopy, and chemical kinetics

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the gas-phase reaction mechanism and associated chemical kinetics are poorly understood. During the performance period we undertook development of a state-of-the-art ab initio gas-phase chemical kinetics simulation workflow to model the hydrolysis of uranium hexafluroride (UF 6 ). In doing so, we addressed several outstanding issues in the theoretical treatment of uranium-containing systems. At the outset it was unclear how to generate accurate estimates of kinetic and thermodynamic data for U-containing chemical reactions. Generation of such data has been made routine. Prior to our work, the literature associated with UF 6 hydrolysis were disparate and inaccurate. This body of work provides a modern and comprehensive theoretical assessment of the reaction mechanism, molecular clustering towards deposition, and chemical kinetics. New methodological implementations and software integrations resulting from this work are also highlighted. As much as possible, our predictions were validated against experimental data including particle morphologies, vibrational spectroscopy, atomization enthalpies, and kinetic rate constants. Nevertheless, we were unable to reconcile kinetic measurements with high-accuracy simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Grain‐Boundary‐Rich Noble Metal Nanoparticle Assemblies: Synthesis, Characterization, and Reactivity

Abstract Here, a comprehensive study on the synthesis, characterization, and reactivity of grain‐boundary (GB)‐rich noble metal nanoparticle (NP) assemblies is presented. A facile and scalable synthesis of Pt, Pd, Au, Ag, and Rh NP assemblies is developed, in which NPs are predominantly connected via Σ3 (111) twin GBs, forming a network. Driven by water electrolysis, the random collisions and oriented attachment of colloidal NPs in solution lead to the formation of Σ3 (111) twin boundaries and some highly mismatched GBs. This synthetic method also provides convenient control over the GB density without altering the crystallite size or GB type by varying the NP collision frequency. The structural characterization reveals the presence of localized tensile strain at the GB sites. The ultrahigh activity of GB‐rich Pt NP assembly toward catalytic hydrogen oxidation in air is demonstrated, enabling room‐temperature catalytic hydrogen sensing for the first time. Finally, density functional theory calculations reveal that the strained Σ3(111) twin boundary facilitates oxygen dissociation, drastically enhancing the hydrogen oxidation rate via the dissociative pathway. This reported large‐scale synthesis of the Σ3 (111) twin GB‐rich structures enables the development of a broad range of high‐performance GB‐rich catalysts.

Geng, Xin↗
DOE OSTI2025

Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning

Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.

Adsorption↗
DOE OSTI2022

Fluorination Enables Simultaneous Improvements of a Dialkoxybenzene-Based Redoxmer for Nonaqueous Redox Flow Batteries

Redoxmers or redox-active organic materials, arc one critical component for nonaqueous redox flow batteries (RFBs), which hold high promise in enabling the time domain of the grid. While tuning redox potentials of redoxmers is a very effective way to enhance energy densities of NRFBs, those improvements often accompany accelerated kinetics of the charged species, undermining stability and cycling performance. In this paper, a strategy for designing redoxmers with simultaneous improvements in redox potential and stability is proposed. Specifically, the redoxmer 1,4-di-tert-butyl-2,5-bis(2,2,2-trifluoroethoxy)benzene (ANL-C46) is developed by incorporating fluorinated substitutions into the dialkoxybenzene-based platform. Compared to the non-fluorinated analogue, ANL-C46 demonstrates not only an increased (similar to 0.41 V) redox potential but also much enhanced stability (1.6 times) and cyclability (4 times) evidenced by electron paramagnetic resonance kinetic study, H-cell and flow cell cycling. In fact, the cycling performance of ANL-C46 is among the best of high potential (>1.0 V vs Ag/Ag + ) redoxmers ever reported. Density functional theory calculations suggest that while the introduced fluorine substitutions elevate the redox potentials, they also help to depress the decomposition reactions of the charged redoxmers, affording excellent stability. The findings represent an interesting strategy for simultaneously improving energy density and stability, which could further prompt the development of high-performance redoxmers.

25 ENERGY STORAGE↗
DOE OSTI2024

17 O NMR Spectroscopy Reveals CO 2 Speciation and Dynamics in Hydroxide-Based Carbon Capture Materials

Carbon dioxide capture technologies are set to play a vital role in mitigating the current climate crisis. Solid-state 17 O NMR spectroscopy can provide key mechanistic insights that are crucial to effective sorbent development. In this work, we present the fundamental aspects and complexities for the study of hydroxide-based CO 2 capture systems by 17 O NMR. We perform static density functional theory (DFT) NMR calculations to assign peaks for general hydroxide CO 2 capture products, finding that 17 O NMR can readily distinguish bicarbonate, carbonate and water species. However, in application to CO 2 binding in two test case hydroxide-functionalised metal-organic frameworks (MOFs) – MFU-4l and KHCO 3 -cyclodextrin-MOF, we find that a dynamic treatment is necessary to obtain agreement between computational and experimental spectra. We therefore introduce a workflow that leverages machine-learning force fields to capture dynamics across multiple chemical exchange regimes, providing a significant improvement on static DFT predictions. In MFU-4l, we parameterise a two-component dynamic motion of the bicarbonate motif involving a rapid carbonyl seesaw motion and intermediate hydroxyl proton hopping. For KHCO 3 -CD-MOF, we combined experimental and modelling approaches to propose a new mixed carbonate-bicarbonate binding mechanism and thus, we open new avenues for the study and modelling of hydroxide-based CO 2 capture materials by 17 O NMR.

NMR spectroscopy↗
DOE OSTI2021

Electronic tuning of confined sub-nanometer cobalt oxide clusters boosting oxygen catalysis and rechargeable Zn–air batteries.

Reasonable design of robust bifunctional oxygen catalysts from an electronic structure perspective is intriguing and challenging for the development of high active rechargeable zinc-air batteries (ZABs). In this study, the favorable regulation of the electronic structure of the cobalt oxide nanoclusters was firstly predicted by density functional theory (DFT) simulation, and then experimentally verified by confining sub-nanometer CoOx clusters (0.86 nm) into the small pore of ZIF-8 derived N-doped nanomaterials (PNC) using a microporous MOFs confinement strategy. The confined effect of the MOF micropores not only enhanced the stability of the subnanometer cobalt oxide clusters, but also make it coupled with Co-Nx to further regulate the electronic structure of the former, synergistic resulting in enhanced ORR/OER actives. As a result, the optimized 0.05CoOx@PNC catalyst demonstrates outstanding bifunctional oxygen performance with a smaller potential gap of 0.67 V. Moreover, the rechargeable Zn-air batteries integrated 0.05CoOx@PNC air cathode displays encouraging performance with a peak power density of 157.1 mW cm(-2), a specific capacity of 887 mAh g(Zn)(-1)at 10 mA cm(-2) and long-term cyclability for over 200 h, significantly outperforming the benchmark electrode couple consisted of Pt/C/RuO2. DFT calculation further revealed that reducing particle size and coupling with Co-N could effectively regulate the charge distribution of CoOx nanoclusters and downshift the D-band center of Co adsorption sites in CoOx nanoclusters, which reduced the reaction barrier of intermediate O-2* and OH* and ORR/ OER over potential, thus accelerating the overall ORR/OER kinetic process. This work offers a novel reference for the construction of a robust sub-nanometer cluster catalysts in the field of ZABs.

Bifunctional oxygen electrocatalysts↗
DOE OSTI2023

Molecular Simulation of Lithium Carbonate Reactive Vapor–Liquid Equilibria Using a Deep Potential Model

We developed a first-principles machine learning model for the reactive vapor–liquid phase behavior of molten Li 2 CO 3 . The model was trained on ab initio electronic density functional theory data using the Deep Potential (DP) methodology, and its accuracy was evaluated by comparing model predictions of density and viscosity to experimental measurements. Direct coexistence simulations with the DP model over time scales of tens of nanoseconds were used to observe equilibrium dissociation of Li 2 CO 3 into CO 2 residing primarily in the vapor phase and Li 2 O which remains dissolved in the liquid. The simulations covered a range of temperatures, overall system sizes, and vapor-to-liquid volume ratios. Results were analyzed in terms of the observed chemical composition of the liquid and vapor phases, product structure, and CO 2 partial pressures. In addition, we calculated equilibrium constants for the dissociation reaction by assuming ideal-solution behavior for the liquid. As expected on the basis of thermodynamic arguments and prior experiments for this system, the observed partial pressure of CO 2 in the gas phase depends on both the temperature and the ratio of vapor to liquid volumes, while the calculated equilibrium constants only depend on temperature. DP model predictions for the equilibrium constant of the reaction are generally consistent with the available experimental measurements. Furthermore, the present study establishes the validity of the DP methodology for the description of reactive, multiphase equilibria from first principles, with possible applications to many other systems of scientific and technological interest even in the absence of relevant experimental measurements.

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

Quantifying the Impact of Electric Fields on the Local Structure and Migration of Potassium Ions at the Orthoclase (001) Surface

Developing an understanding of the response of mineral/water interfaces to applied electric fields is central to detecting and interpreting signatures of interfacial processes in the subsurface. Here we focus density functional theory calculations on understanding K + cation binding and diffusion across the (001) surface of orthoclase feldspar under various applied electric fields with and without surface hydration. The calculations reveal how water ligands labilize surface K + cations for diffusion while also increasing their sensitivity to electric field effects on the binding energy at different surface sites. The calculations also show how the direction and strength of the electric field systematically affect surface cation mobility, sorption, and hydration behavior. Specifically, electric fields directed toward the surface increase the potassium binding energies, facilitate K + diffusion across multiple binding sites, and increase the propensity of surface K + ions to lose their hydration shells. Furthermore, the findings help fill a major knowledge gap into the impact of electric fields on mineral/water interface structure and dynamics more generally, featuring in this case a commonly found type of feldspar involved in a multitude of atmospheric and geochemical processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗