Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRENGTH CALCULATIONS”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Insensitivity of Magnetic Coupling to Ligand Substitution in a Series of Tetraoxolene Radical-Bridged Fe 2 Complexes

The elucidation of magnetostructural correlations between bridging ligand substitution and strength of magnetic coupling is essential to the development of high-temperature molecule-based magnetic materials. Toward this end, we report the series of tetraoxolene-bridged Fe II 2 complexes [(Me 3 TPyA) 2 Fe 2 ( R L)] n + (Me 3 TPyA = tris(6-methyl-2-pyridylmethyl)amine; n = 2: OMe LH 2 = 3,6-dimethoxy-2,5-dihydroxo-1,4-benzoquinone, Cl LH 2 = 3,6-dichloro-2,5-dihydroxo-1,4-benzoquinone, Na 2 [ NO 2 L] = sodium 3,6-dinitro-2,5-dihydroxo-1,4-benzoquinone; n = 4: SMe 2 L = 3,6-bis(dimethylsulfonium)-2,5-dihydroxo-1,4-benzoquinone diylide) and their one-electron-reduced analogues. Furthermore, variable-temperature dc magnetic susceptibility data reveal the presence of weak ferromagnetic superexchange between Fe II centers in the oxidized species, with exchange constants of J = +1.2(2) (R = OMe, Cl) and +0.3(1) (R = NO 2 , SMe 2 ) cm –1 . In contrast, X-ray diffraction, cyclic voltammetry, and Mössbauer spectroscopy establish a ligand-centered radical in the reduced complexes. Magnetic measurements for the radical-bridged species reveal the presence of strong antiferromagnetic metal–radical coupling, with J = -57(10), -60(7), -58(6), and -65(8) cm –1 for R = OMe, Cl, NO 2 , and SMe 2 , respectively. The minimal effects of substituents in the 3- and 6-positions of R L x –• on the magnetic coupling strength is understood through electronic structure calculations, which show negligible spin density on the substituents and associated C atoms of the ring. Finally, the radical-bridged complexes are single-molecule magnets, with relaxation barriers of U eff = 50(1), 41(1), 38(1), and 33(1) cm –1 for R = OMe, Cl, NO 2 , and SMe 2 , respectively. Taken together, these results provide the first examination of how bridging ligand substitution influences magnetic coupling in semiquinoid-bridged compounds, and they establish design criteria for the synthesis of semiquinoid-based molecules and materials.

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

Molecular Insight into the β-Sheet Twist and Related Morphology of Self-Assembled Peptide Amphiphile Ribbons

Self-assembly of high-aspect-ratio filaments containing β-sheets has attracted much attention due to potential use in bioengineering and biomedicine. However, precisely predicting the assembled morphologies remains a grand challenge because of insufficient understanding of the self-assembly process. We employed an atomistic model to study the self-assembly of peptide amphiphiles (PAs) containing valine–glutamic acid (VE) dimeric repeats. By changing of the sequence length, the assembly morphology changes from flat ribbon to left-handed twisted ribbon, implying a relationship between β-sheet twist and strength of interstrand hydrogen bonds. The calculations are used to quantify this relationship including both magnitude and sign of the ribbon twist angle. Interestingly, a change in chirality is observed when we introduce the RGD epitope into the C-terminal of VE repeats, suggesting arginine and glycine’s role in suppressing right-handed β-sheet formation. Furthermore, this study provides insight into the relationship between β-sheet twist and self-assembled nanostructures including a possible design rule for PA self-assembly.

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

Prediction of O and OH Adsorption on Transition Metal Oxide Surfaces from Bulk Descriptors

In the search for stable and active catalysts, density functional theory and machine learning (ML) based models can accelerate the screening of materials. While stability is conveniently addressed on the bulk level of computation, the modelling of catalytic activity requires expensive surface simulations. Here, in this work, we develop models for the surface adsorption energy of O and OH intermediates across a consistent and extensive dataset of pure transition metal oxide surfaces. We show that adsorption energies across metal oxidation states of +2 to +6 are well captured from the metal-oxygen bond strength extracted from the bulk level calculation. Specifically, we calculate the integrated crystal orbital Hamiltonian population (ICOHP) of the metal-oxygen bond in the bulk oxide and employ a simple normalization scheme to obtain a strong correlation with adsorption energetics. By combining our ICOHP descriptor with non DFT features in a Gaussian Process regression (GPR) model, we achieve high model accuracy with mean absolute errors of 0.166 and 0.219 eV for OH and O adsorption, respectively. By targeting the O-OH adsorption energy difference with our GPR model, we predict the the oxygen evolution reaction (OER) activity from bulk descriptors only. Furthermore, we utilize the strong correlation between the COHP and metal oxygen bond lengths to rapidly predict adsorption energetics and catalytic activity from the optimized bulk geometry. Our approach can enable an efficient search for active catalysts by eliminating the need for surface calculations in the initial screening phase.

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

Valence shell electronically excited states of norbornadiene and quadricyclane

The absolute photoabsorption cross sections of norbornadiene (NBD) and quadricyclane (QC), two isomers with chemical formula C 7 H 8 that are attracting much interest for solar energy storage applications, have been measured from threshold up to 10.8 eV using the Fourier transform spectrometer at the SOLEIL synchrotron radiation facility. The absorption spectrum of NBD exhibits some sharp structure associated with transitions into Rydberg states, superimposed on several broad bands attributable to valence excitations. Sharp structure, although less pronounced, also appears in the absorption spectrum of QC. Assignments have been proposed for some of the absorption bands using calculated vertical transition energies and oscillator strengths for the electronically excited states of NBD and QC. Natural transition orbitals indicate that some of the electronically excited states in NBD have a mixed Rydberg/valence character, whereas the first ten excited singlet states in QC are all predominantly Rydberg in the vertical region. In NBD, a comparison between the vibrational structure observed in the experimental 1 1 B 1 –1 1 A 1 (3sa 1 ← 5b 1 ) band and that predicted by Franck–Condon and Herzberg–Teller modeling has necessitated a revision of the band origin and of the vibrational assignments proposed previously. Similar comparisons have encouraged a revision of the adiabatic first ionization energy of NBD. Simulations of the vibrational structure due to excitation from the 5b 2 orbital in QC into 3p and 3d Rydberg states have allowed tentative assignments to be proposed for the complex structure observed in the absorption bands between ~5.4 and 7.0 eV.

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

Crust Composition and the Shallow Heat Source in KS 1731–260

The presence of a strong shallow heat source of unknown origin in accreting neutron star crusts has been inferred by analyzing x-ray observations of their cooling in quiescence. Here, we model the cooling of KS 1731–260 using realistic crust compositions and nuclear heating and cooling sources from detailed nuclear reaction network calculations. We find that the required strength of the shallow heat source in KS 1731–260 is reduced by more than a factor of 3 compared to previous analysis, a 5⁢𝜎 difference that alleviates the need for exotic solutions. Our analysis also suggests the existence of an impure nuclear pasta layer in the inner crust of KS 1731–260, though future observations will provide more stringent constraints. In addition, we obtain constraints on the dominant surface burning modes of KS 1731–260 over its history.

neutron star cooling↗
DOE OSTI2023

Cross-section measurement of the 82 Kr ($p,γ$) 83 Rb reaction in inverse kinematics

The total cross section of the 82 Kr(p,γ) 83 Rb reaction was measured for the first time at effective center-of-mass energies between 2.4 and 3.0 MeV, within the relevant Gamow window for the astrophysical γ process. The experiment took place at the National Superconducting Cyclotron Laboratory at Michigan State University using the ReA facility. A 82 Kr beam was directed onto a hydrogen gas cell located at the center of the Summing NaI(Tl) (SuN) detector. Here, the obtained spectra were analyzed using the γ-summing technique and the extracted cross section was compared to standard statistical model calculations using the non-smoker and talys codes. The comparison indicates that standard statistical model calculations tend to overproduce the cross section of the 82 Kr(p,γ) 83 Rb reaction relative to the experimentally measured values. Furthermore, the experimental data were used to provide additional constraints on the nuclear level density and the γ-ray strength function used in the statistical model calculations.

59 ≤ A ≤ 89↗
DOE OSTI2025

Oakridge PPU Magnets: Results and Measurements

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) is being upgraded from 1.0 GeV to 1.3 GeV (or 1.4 to 2.8 MW). Several water-cooled magnets have been upgraded to transport 30% higher beam energy. Fermilab contributed the magnet design for the new chicane magnets and injection/extraction septum. Designing the magnets was a challenging task because the new magnets required good combined integrated field quality and needed to occupy the old magnets space but with about 20% greater integrated magnetic field. Additional strong requirements applied to the magnets fringe field so as not to disturb the circulating beam. After fabrication of the magnets, an extensive measurement campaign was developed and performed at Fermilab’s Magnet Test Facility. The measurements needed to assess magnet performance and provide comparison to design calculations. These included verification of field strength and harmonics along an 8 m length and 200 mm good field diameter for the chicane dipoles, end-field Hall probe mapping of these magnets, and measurements along two differently curved trajectories within the ∼3 m septum gradient magnet. Details of the measurements and systems are presented along with results and comparison to field models.

DiMarco, J. [Fermilab] (ORCID:0000000228400036)↗
DOE OSTI2025

Solenoidal Fields from Permanent Magnets for a Polarized Electron Preinjector

This technical note describes the design of a compact focusing solenoid constructed from permanent bar magnets, intended for the low-energy section of a spin-polarized electron injector. Typically, quadrupole triplets provide appropriate focusing for the low-energy spin rotator section; however, their large footprint precludes the insertion of sufficient beam diagnostics. The primary motivation of this study is to provide beam focusing while preserving the vertical spin orientation of the electron beam. We propose a solenoidal field composed of multiple permanent magnets arranged to naturally yield a zero longitudinal magnetic field integral ($\int B_zdz = 0$). This design serves as a compact, passive, and non-coupling alternative to traditional quadrupole triplets. We present the magnet arrangement, field calculation methodology, spin tracking simulations using space charge codes, and calculations demonstrating that the required field strength is achievable.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

VENTSAR V3.0: A Python GUI for Estimating Contaminant Concentrations on or Near Buildings Due to Building Effects and Plume Rise and For Calculating Inhalation and Plume Shine Doses

VENTSAR: Originated as VENTAX (Smith and Weber 1983) on the IBM Mainframe at the Savannah River Site (SRS) as a Fortran program. Updated to VENTSAR XL V1.0 (Simpkins 1997) as a spreadsheet version using macros created in Microsoft Excel. Later updated to VENTSAR XL V2.0 (Dixon 2018) as a spreadsheet executing a modified version of the original VENTAX Fortran code. Estimates contaminant concentrations on or near a building from a release at a nearby location. Calculates concentrations for a given meteorological exceedance probability or for a given stability and wind speed combination. Can model a single building with or without a penthouse on top or a ground location from either a stack or ground release. Plume rise can be considered. Contaminant releases can be chemical or radioactive with downwind concentrations determined at user-specified distances. Wind passing over and around buildings creates a complicated dispersion pattern. Air-intake vents may be located on building roofs or near the ground downwind of a release source. An estimation of pollutant concentrations on or near a structure is important in determining expected pollutant levels. Meteorological data are selected based on a specific area of the site. Fortran-based VENTAX able to make fast, complex mathematical calculations. Not user-friendly VENTSAR XL V1.0 no longer supported due to obsolete Excel macros. VENTSAR XL V2.0 an interim solution using an Excel spreadsheet to interface with the modified VENTAX Fortran code. Requires user to have Microsoft Excel. Not an intuitive interface for someone unfamiliar with VENTSAR. Does not perform dose calculations. VENTSAR V3.0 goal to combine the strengths of previous versions of VENTSAR into a single, powerful yet user-friendly program with a Graphical User Interface (GUI). Not dependent upon a specific program or plug-ins. Self-contained package to be used on any computer running Microsoft Windows. User-friendly, intuitive GUI created in Python for parameter input. Executes same modified VENTAX Fortran code as VENTSAR XL V2.0. Outputs formatted text file of completed calculations for easy review and dissemination. Performs inhalation and plume shine dose calculations for up to 11 user-selected nuclides from a nuclide dose factor library of almost 500 nuclides. 12 test cases were created for verification of V1.0 and V2.0. Same test cases run in V3.0 to verify correct performance. V3.0 produced similar results to V1.0. Confirmed GUI did not alter calculations in any way. Comparisons of the test case results confirm the V3.0 GUI does not influence the VENTSAR calculations. Simply passes same input parameters to VENTAX Fortran code for execution. Provides end user an easy-to-use, intuitive tool to quickly make building effect and plume rise calculations, as well as, inhalation and plume shine dose calculations. No experience with or working knowledge of Fortran, Excel, command line, or macros required. Self-contained package allows VENTSAR be deployed to any Windows computer without the need for specific programs or plug-ins.

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

Probing the 𝑗 dependence of angular distributions and 𝑁=20 shell rigidity via the 36 S (𝑝,𝑑) 35 S reaction

An investigation of the N=20 36 S nucleus has been performed through a detailed study of the 3 6S (p,d)3 5S neutron-removal reaction, employing a 66-MeV proton beam at iThemba Laboratory for Accelerator Based Sciences and an innovative target design. Using the high-resolution K = 600 magnetic spectrometer, 98 states in 3 5S were identified up to 16-MeV excitation energy, including 47 previously unobserved states. Angular distributions and spectroscopic factors, including isobaric analog-state contributions, were extracted for 81 levels. A pronounced j dependence in the angular distributions of ℓ=2 states provides refined insights into the spin-orbit splitting. Finite-range adiabatic distorted-wave approximation calculations qualitatively reproduce the observed j dependence. Comparisons of the measured 1d 5/2 spectroscopic strength distribution with large-scale shell-model and ab initio calculations show good agreement overall, and the robustness of the N=20 shell closure in 36 S is confirmed when comparing the relatively low fp orbital occupancies in 40 Ca and 3 6S across the Fermi surface. This study underscores the utility of neutron-removal reactions in probing nuclear structure and the Fermi surface of sd nuclei and beyond. The findings advance our understanding of shell evolution and offer constraining data for theoretical models.

20 ≤ A ≤ 38↗
DOE OSTI2020

Ab initio polaritonic potential-energy surfaces for excited-state nanophotonics and polaritonic chemistry

Advances in nanophotonics, quantum optics, and low-dimensional materials have enabled precise control of light–matter interactions down to the nanoscale. Combining concepts from each of these fields, there is now an opportunity to create and manipulate photonic matter via strong coupling of molecules to the electromagnetic field. Toward this goal, here we demonstrate a first principles framework to calculate polaritonic excited-state potential-energy surfaces, transition dipole moments, and transition densities for strongly coupled light–matter systems. In particular, we demonstrate the applicability of our methodology by calculating the polaritonic excited-state manifold of a formaldehyde molecule strongly coupled to an optical cavity. This proof-of-concept calculation shows how strong coupling can be exploited to alter photochemical reaction pathways by influencing avoided crossings with tuning of the cavity frequency and coupling strength. Therefore, by introducing an ab initio method to calculate excited-state potential-energy surfaces, our work opens a new avenue for the field of polaritonic chemistry.

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

Magnetic states of the quasi-one-dimensional iron chalcogenide Ba 2 FeS 3

Quasi-one-dimensional iron-based ladders and chains, with the 3d iron electronic density n=6, are attracting considerable attention. Recently, a new iron chain system Ba 2 FeS 3 , also with n=6, was prepared under high-pressure and high-temperature conditions. Here the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered ↑-↓-↑-↓ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard U coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large U/W (W=bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba 2 FeS 3 provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material.

1-dimensional systems↗
DOE OSTI2023

High-throughput exploration of the WMoVTaNbAl refractory multi-principal-element alloys under multiple-property constraints

Development of next-generation gas turbines requires the design and fabrication of novel high-temperature structural materials capable of operating beyond 1300°C. Here, we propose a high-throughput alloy design framework under multiple-property constraints to discover new refractory multi-principal element alloys (MPEAs) for high-temperature applications. The framework treats the development of MPEAs as a composition-agnostic constraint satisfaction problem, i.e., no prescriptions are made concerning the design space before performing investigatory calculations. We target alloys in the WMoVTaNbAl chemistry space that are predicted to meet constraints on the following properties simultaneously: single-phase stability, density, solidus temperature, yield strength at 1300°C, and ductile-to-brittle-transition temperature. These properties are relevant to both applications in gas turbines and manufacturability. A set of 214 MoNbV-rich alloys meet these relevant constraints. These feasible alloys are investigated with density functional theory (DFT) to provide a fundamental electronic basis for their superior properties. Three compositionally representative alloys from the feasible design space (Mo 45 Nb 35 Ta 5 V 15 , Mo 25 Nb 50 V 20 W 5 , and Mo 30 Nb 35 Ta 5 V 25 W 5 ) are selected with a k-medoids-based design scheme for detailed DFT analysis and experimental characterization. The DFT analysis predicted a single-phase BCC at high temperatures with a high yield strength for all three MPEAs, in agreement with CALPHAD (CALculation of PHAse Diagrams) and experiments, respectively. These three alloys are benchmarked against a public database of 1546 MPEAs. Concerning the aforementioned constraints, the Mo 30 Nb 35 Ta 5 V 25 W 5 alloy outperforms these 1546 MPEAs. The present work demonstrates the ability of the proposed design methodology to identify candidate alloys for a given application under multiple property constraints in a combinatorically vast design space.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Theoretical and experimental studies of compression and shear deformation behavior of Osmium to 280 GPa

Abstract The compression behavior of osmium metal was investigated up to 280 GPa (volume compression V/V o = 0.725) under nonhydrostatic conditions at ambient temperature using angle dispersive axial x-ray diffraction (A-XRD) with a diamond anvil cell (DAC). In addition, shear strength of osmium was measured to 170 GPa using radial x-ray diffraction (R-XRD) technique in DAC. Both diffraction techniques in DAC employed platinum as an internal pressure standard. Density functional theory (DFT) calculations were also performed, and the computed lattice parameters and volumes under compression are in good agreement with the experiments. DFT predicts a monotonous increase in axial ratio ( c/a ) with pressure and the structural anomalies of less than 1% in ( c/a ) ratio reported below 150 GPa were not reproduced in theoretical calculations and hydrostatic measurements. The measured value of shear strength of osmium ( τ ) approaches a limiting value of 6 GPa above a pressure of 50 GPa in contrast to theoretical predictions of 24 GPa and is likely due to imperfections in polycrystalline samples. DFT calculations also enable the studies of shear and tensile deformations. The theoretical ideal shear stress is found along the (001)[1–10] shear direction with the maximal shear stress ∼24 GPa at critical strain ∼0.13.

Lin, Chia-Min↗
DOE OSTI2020

Nonrandom two‐liquid activity coefficient model with association theory

Abstract Considered one of the most versatile and widely used classical thermodynamic models to correlate phase equilibria behavior of nonideal systems, the nonrandom two‐liquid (NRTL) theory does not explicitly account for specific chemical associations such as hydrogen bonding. This deficiency has been recognized as the cause for unsatisfactory representation of association systems like methanol–alkanes binaries. This work presents a practical approach to integrate Wertheim's perturbation theory for association contribution with the classical NRTL model. Specifically, the association contribution is calculated from pre‐determined molecule‐specific association strengths while the physical interaction contribution is captured with NRTL binary interaction parameters. The resulting association NRTL model correlates fluid phase equilibria of association systems with few adjustable parameters and offers improved predictive capability for higher order systems.

Hao, Yifan↗
DOE OSTI2022

Superhydrophilicity of $α$-alumina surfaces results from tight binding of interfacial waters to specific aluminols

Understanding the microscopic driving force of water wetting is challenging and important for design of materials. The relations between structure, dynamics and hydrogen bonds of interfacial water can be investigated using molecular dynamics simulations. Here, contact angles at the alumina (0001) and ($11\bar{2}0$) surfaces are studied using both classical molecular dynamics simulations and experiments. To test the superhydrophilicity, the free energy cost of removing waters near the interfaces are calculated using the density fluctuations method. The strength of hydrogen bonds is determined by their lifetime and geometry. Both surfaces are superhydrophilic and the (0001) surface is more hydrophilic. Interactions between surfaces and interfacial waters promote a templating effect whereby the latter are aligned in a pattern that follows the underlying lattice of the surfaces. Translational and rotational dynamics of interfacial water molecules are slower than in bulk water. Hydrogen bonds between water and both surfaces are asymmetric, water-to-aluminol ones are stronger than aluminol-to-water ones. Molecular dynamics simulations eliminate the impacts of surface contamination when measuring contact angles and the results reveal the microscopic origin of the macroscopic superhydrophilicity of alumina surfaces: strong water-to-aluminol hydrogen bonds.

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

Recovering high-purity uranyl nitrate from simulated used nuclear fuel dissolver solutions by crystallization: rejecting technetium

The separation of U from Tc and other problematic fission product elements like Mo and Ru, along with Sr, Zr, Cs, and Nd, has been achieved via the crystallization of uranyl nitrate hexahydrate (UNH). Rejection of technetium as pertechnetate anion ( 99 TcO 4 – ) is an especially important feature of this system, as it otherwise tends to follow U (VI) in extractive separations. It also raises the salient question regarding why this oxoanion cannot replace nitrate within the crystalline lattice of UNH. Results showed high-yield (>90 %), high-purity (>99 %) recovery of U as UNH from solutions containing 99 TcO 4 – by simple reduction of temperature from 60°C to 20°C. There was no observable interaction of 99 TcO 4 – with UO 2 2+ . The addition of other cations like, Sr 2+ , Zr 4+ , Cs + , and Nd 3+ , also did not form secondary, contaminant solid phases, leaving the > 99 % of the fission product elements in the mother liquor, while the U was recovered at > 90 %. Similarly, Mo and Ru, when added to the mixture, were shown to behave as the other fission-product elements, remaining in the mother liquor during crystallization. As a result, DFT calculations showed that, despite the higher binding strength of TcO 4 – , HMoO 4 – , and BiO 3 – with the UO 2 2+ cation compared to NO 3 – , the hydrogen-bonding network of the two coordinated ions and four waters of hydration in the UNH crystal structure is the driving force for the high specificity of this separation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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↗