Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electron localization function”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Atomically Revealing Bulk Point Defect Dynamics in Hydrogen‐Driven γ‐Fe 2 O 3 → Fe 3 O 4 → FeO Transformation

Understanding how point defects in the bulk govern redox transformations is essential for advancing hydrogen-based metal production and designing high-performance oxide materials. This study reveals the atomic-scale mechanisms driving hydrogen-induced reduction of γ-Fe 2 O 3 to Fe 3 O 4 , focusing on how bulk vacancy dynamics dictate structural evolution and reaction kinetics. A key finding is the pronounced contrast in defect behavior between the two oxides: in γ-Fe 2 O 3 , intrinsic Fe vacancies promote oxygen vacancy clustering, destabilizing the local lattice and driving nanopore formation. In contrast, Fe 3 O 4 exhibits a higher oxygen vacancy formation energy and lacks intrinsic Fe vacancies, suppressing vacancy aggregation and maintaining a dense, pore-free structure. This divergence governs distinct reduction pathways—γ-Fe 2 O 3 undergoes an interface-reaction-limited transformation confined to the γ-Fe 2 O 3 /Fe 3 O 4 boundary, while Fe 3 O 4 supports a uniform increase in oxygen vacancy concentration, enabling bulk-phase reduction to lower-oxide FeO. Integrated in situ electron microscopy and density functional theory modeling uncover a vacancy-mediated mechanism, where synergistic cation-anion vacancy dynamics steer microstructure evolution and phase progression. These insights highlight the critical role of vacancy dynamics in controlling oxide reactivity and offer a pathway toward vacancy engineering to enhance reduction kinetics in hydrogen metallurgy and to tailor porosity, reactivity, and structural resilience in oxide-based catalysts and energy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Confinement of quasi-atomic structures in Ti 2 N and Ti 3 N 2 MXene electrides

Metal carbides, nitrides, or carbonitrides of early transition metals, better known as MXenes, possess notable structural, electrical, and magnetic properties. Analyzing electronic structures by calculating structural stability, band structure, density of states, Bader charge transfer, and work functions utilizing first principle calculations, we revealed that titanium nitride MXenes, namely Ti 2 N and Ti 3 N 2 , have excess anionic electrons in their lattice voids, making them MXene electrides. Bulk Ti 3 N 2 has competing antiferromagnetic (AFM) and ferromagnetic (FM) configurations with slightly more stable AFM configuration, while the Ti 2 N MXene is nonmagnetic. Although Ti 3 N 2 favors AFM configuration with hexagonal crystal systems having 6/ mmm point group symmetry, Ti 3 N 2 does not support altermagnetism. The monolayer of the Ti 3 N 2 MXene is a ferromagnetic electride. These unique properties of having non-nuclear interstitial anionic electrons in the electronic structure of titanium nitride MXene have not yet been reported in the literature. Density functional theory calculations show TiN is neither an electride, MXene, or magnetic.

Anionic electrons↗
DOE OSTI2025

Holistic Microstructure Control Strategies in Photopolymerization‐Induced Phase Separation of Acrylate Systems

Open porous materials, known for their large surface area and interconnected structures, are essential in various applications, including batteries, ion exchange, catalysis, filtration, and electronic waste recycling. A critical aspect of the functionality of porous membranes is the precise control of pore size and morphology. Photopolymerization-induced phase separation (photo-PIPS) offers a convenient and versatile methods for creating porous structures. However, controlling the porous morphology remains challenging due to the complex interplay between thermodynamics, polymerization kinetics, and monomer structures, which makes it difficult to establish the relationship between processing conditions and resulting morphology in photo-PIPS. Herein, a physics-based phase-field model capable of generating and characterizing the microstructures of porous materials based on both average and localized features is developed. Using the phase-field simulations as test bed, the effects of polarity, light intensity, and curing temperature, as well as the previously unexplored roles of chain transfer agents and substrates, on the morphology of the resulting porous microstructure are investigated. Experiments are performed to verify the results predicted by the simulations. This work lays out a comprehensive guide for designing PIPS-derived porous microstructures and offers practical engineering strategies for tailoring microstructure-level topology and size of pores for application-specific needs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Straightforward Model for Quantifying Local pH Gradients Governing the Oxygen Evolution Reaction

The production and consumption of protons by an electrocatalyst will, under certain conditions, generate localized microenvironments with properties distinct from those of the bulk solution. These local properties are particularly impactful for reactions involving proton-coupled electron transfer, where the generation of locally basic or acidic environments may significantly influence the energy efficiency and reaction selectivity of the electrocatalyst. Whereas local pH environments have been observed and characterized in reductive half-reactions, including the CO 2 reduction and hydrogen evolution reactions, the incompatibility of conventional techniques and materials has limited studies in oxidative half-reactions, including the oxygen evolution reaction (OER), which provides the reducing equivalents for solar-to-fuels electrolysis. With the straightforward parameters bulk pH, buffer composition and pK a , and mass transport, we develop a model for describing local pH as a function of current density regardless of the microscopic details of the mechanism. Using an acid-stable PbO x OER catalyst, we observe the formation and dissipation of pH gradients during the OER and validate the model with voltammetric and potentiometric studies. Here, the model predicts how local acidic environments can develop over a narrow OER current density window, thus providing further motivation for the development of OER catalysts that are stable to acid, even when operating in basic aqueous conditions. More generally, the model is not restricted to the OER and is useful for determining the onset of local pH gradients for other electrocatalytic reactions that involve the consumption or generation of protons in energy conversion reactions.

Anions↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2026

Multireference Methods for Chemistry and Materials Science: Automated Active Spaces, Efficient Dynamic Correlation, and Extended Systems

While multiconfigurational approaches have long been relegated to expert practitioners working on a case-by-case basis, recent developments have increasingly made these methods more routine and applicable to broader sets of systems. This article outlines the state-of-the-art in multiconfigurational approaches, with an emphasis on moving from delicate hand-selected pathways through configuration space toward more robust and efficient approaches to treating a host of challenging chemical systems accurately. First, we overview recent work in automated active-space selection, which has enabled increasingly large-scale applications of multireference methods to modeling vertical excitations and reactivity. Second, we highlight the increasingly efficient methods for recovering correlation energy beyond the active space, as headlined by extensions of pair-density functional theory and its role in accurate and efficient treatment of excited-state dynamics and its utilization to train machine-learned potentials. Finally, we highlight recent efforts to treat extended systems that until recently have lied beyond the traditional limits of active-space methods, giving center stage to product-form wave functions of the localized active space family of methods that allow for the computation of multiconfigurational band structures. These recent advancements point to a broader use of multireference approaches for high-impact chemical and materials science applications.

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

Tracking Local pH Dynamics during Water Electrolysis via In-Line Continuous Flow Raman Spectroscopy

The performance of electrochemical devices, which play a critical role in decarbonization efforts, is often governed by proton-coupled electron transfer reactions at the electrode–electrolyte interface. These reactions are highly sensitive to the complex and dynamic microenvironment present at the electrode surface. However, characterizing this environment─particularly monitoring interfacial pH and its evolution under reaction conditions─remains challenging, necessitating the development of advanced analytical tools. Here, in this study, we introduce in-line continuous flow Raman spectroscopy (CFRS) as a spectroelectrochemical platform for quantifying interfacial pH swings generated during water-splitting. By monitoring phosphate ion speciation and controlling the hydrodynamics with a flow cell, we measure pH swings as a function of current density, flow rate, and distance from the electrode. Comparison with theoretical models reveals the impact of bulk pH, boundary layer thickness, and bubble dynamics at high current densities. Collectively, these findings establish CFRS as a platform for quantitatively investigating pH dynamics, offering critical insights for advancing electrochemical energy conversion technologies.

Marquez, Raul A. [Univ. of Texas, Austin, TX (Unit↗
DOE OSTI2025

2D Vacancy Confinement in Anatase TiO 2 for Enhanced Photocatalytic Activities

Abstract Light‐driven energy conversion devices call for the atomic‐level manipulation of defects associated with electronic states in solids. However, previous approaches to produce oxygen vacancy (V O ) as a source of sub‐bandgap energy levels have hampered the precise control of the distribution and concentration ofV O . Here, a new strategy to spatially confineV O at the homo‐interfaces is demonstrated by exploiting the sequential growth of anatase TiO 2 under dissimilar thermodynamic conditions. Remarkably, metallic behavior with high carrier density and electron mobility is observed after sequential growth of the TiO 2 films under low pressure and temperature (L‐TiO 2 ) on top of high‐quality anatase TiO 2 epitaxial films (H‐TiO 2 ), despite the insulating properties of L‐TiO 2 and H‐TiO 2 single layers. Multiple characterizations elucidate that theV O layer is geometrically confined within 4 unit cells at the interface, along with low‐temperature crystallization of upper L‐TiO 2 films; this 2DV O layer is responsible for the formation of in‐gap states, promoting photocarrier lifetime (≈300%) and light absorption. These results suggest a synthetic strategy to locally confine functional defects and emphasize how sub‐bandgap energy levels in the confined imperfections influence the kinetics of light‐driven catalytic reactions.

Chemistry↗
DOE OSTI2025

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

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

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spectroscopic Signatures of Phonon Character in Molecular Electron Spin Relaxation

Spin–lattice relaxation constitutes a key challenge for the development of quantum technologies, as it destroys superpositions in molecular quantum bits (qubits) and magnetic memory in single molecule magnets (SMMs). Gaining mechanistic insight into the spin relaxation process has proven challenging owing to a lack of spectroscopic observables and contradictions among theoretical models. Here, we use pulse electron paramagnetic resonance (EPR) to profile changes in spin relaxation rates (T 1 ) as a function of both temperature and magnetic field orientation, forming a two-dimensional data matrix. For randomly oriented powder samples, spin relaxation anisotropy changes dramatically with temperature, delineating multiple regimes of relaxation processes for each Cu(II) molecule studied. We show that traditional T 1 fitting approaches cannot reliably extract this information. Single-crystal T 1 anisotropy experiments reveal a surprising change in spin relaxation symmetry between these two regimes. We interpret this switch through the concept of a spin relaxation tensor, enabling discrimination between delocalized lattice phonons and localized molecular vibrations in the two relaxation regimes. Variable-temperature T 1 anisotropy thus provides a unique spectroscopic method to interrogate the character of nuclear motions causing spin relaxation and the loss of quantum information.

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

Plasmonic Hot-Carrier Generation and Catalysis in Ti 3 C 2 O 2 from Real-Time TDDFT Simulations

Photoinduced hot electrons are central to plasmon-driven catalysis. Atomically thin Ti 3 C 2 O 2 , with high carrier density and broad optical absorption, offers a promising platform for plasmon-driven reactions. However, comprehensive investigations of its plasmon resonance, hot-carrier generation, and plasmonic catalytic performance remain limited. In this work, real-time time-dependent density functional theory (rt-TDDFT) was employed to study Ti 3 C 2 O 2 ’s plasmon excitation and hot-carrier generation from nonradiative plasmon damping. The temporal evolution of the dipole moment reveals plasmon resonance in Ti 3 C 2 O 2 , followed by strong plasmon damping that redistributes the stored energy to generate hot carriers. Ti 3 C 2 O 2 with low oxygen vacancy concentration (O v -Ti 3 C 2 O 2 ) exhibits plasmonic behavior resembling the pristine surface, and the plasmon-generated hot electrons can markedly reduce the dissociation barrier of CO 2 at the oxygen vacancy. These findings provide fundamental insights into the plasmonic properties of Ti 3 C 2 O 2 and how they drive its catalytic performance in surface reactions, which is valuable for advancing plasmon-driven catalysis.

CO2 reduction↗
DOE OSTI2025

Low‐Loss Far‐Infrared Surface Phonon Polaritons in Suspended SrTiO 3 Nanomembranes

Phonon polaritons (PhPs), excitations arising from the coupling of light with lattice vibrations, enable light confinement and local field enhancement, which is essential for various photonic and thermal applications. To date, PhPs with high confinement and low loss are mainly observed in the mid-infrared regime and mostly in manually exfoliated flakes of van der Waals (vdW) materials. In this work, the existence of low-loss, thickness-tunable phonon polaritons in the far-infrared regime within transferable freestanding SrTiO 3 membranes synthesized through a scalable approach, achieving high figures of merit is demonstrated, which are comparable to the previous record values from the vdW materials. Leveraging atomic precision in thickness control, large dimensions, and compatibility with mature oxide electronics, functional oxide membranes present a promising large-scale 2D platform alternative to vdW materials for on-chip polaritonic technologies in the infrared regime.

functional oxide membranes↗
DOE OSTI2024

Defect-Limited Mobility and Defect Thermochemistry in Mixed A-Cation Tin Perovskites: (CH3NH3)1-xCsxSnBr3

Hybrid organic-inorganic semiconductors crystallizing in the perovskite structure present a significant opportunity for realizing defect-tolerant semiconductors. In this work, we examine the solid solution, (CH3NH3)1-xCsxSnBr3, and identify the thermochemistry dictating the intrinsic carrier concentrations and how local structural distortions influences this electronic behavior. This family of compounds exhibits the expected systematic trend in decreasing optical gap with the cesium to methylammonium A-site mixing ratio, x, in the visible region. However, the carrier mobility, as determined from time-resolved microwave conductivity measurements, trends opposite to that expected from first-principles calculations combined with Boltzmann scattering theory calculations of the carrier mobility. We propose that this is a result of increasing carrier scattering with x in (CH3NH3)1-xCsxSnBr3 due to a significant increase in the carrier density with x. By examining the dependence of the carrier density as a function of x, we infer the compositional-dependence of the average enthalpy and nonconfigurational entropy per defect. While diffraction reveals a cubic aristotypic perovskite structure as a function of x at room temperature, the pair distribution functions obtained from synchrotron X-ray total scattering from these materials are better described by symmetry-adapted displacement modes of the Pm3m crystal structure, which we attribute to a large degree of anharmonic dynamics of the atom positions. This analysis shows that CH3NH3+-rich compositions retain mostly linear Sn-Br-Sn bonding environments through the displacements, while Cs+-rich compositions lead to more significantly bent Sn-Br-Sn environments. We propose that these bent bonding arrangements in Cs-rich compositions yield a higher propensity for defect formation. This also provides a rationale for carrier trapping that gives rise to anomalous microwave transients. Together, these results provide insight into the structure-dynamics-properties relationships in this highly anharmonic system with high amplitude atomic motions and low defect formation energies.

defect thermochemistry↗
DOE OSTI2025

Impact of structural distortions on the correlated electronic structure of orbital-selective Mott insulating Na 3 ⁢Co 2 ⁢SbO 6 under strain

Na 3 Co 2 SbO 6 is a promising candidate to realize the Kitaev spin liquid phase since the large Kitaev spin exchange interaction is tunable via the change in electronic structure, such as the trigonal crystal field splitting (Δ TCF ). Here, in this study, we show that the uncorrelated electronic structure of Na 3 Co 2 SbO 6 is rather insensitive to the strain effect due to the low crystal symmetry accompanied by oxygen displacements and the presence of Sb s orbitals. Using density functional theory plus dynamical mean field theory, we find that the correlated electronic structure of Na 3 Co 2 SbO 6 is an orbital-selective Mott insulating state where the trigonal a 1g orbital is insulating due to the nearly full occupation, while other d orbitals behave as typical Mott insulators, resulting in the effective tunability of Δ TCF under the strain effect. The sign change of Δ TCF can occur as the in-plane tensile strain is applied, and the Kitaev spin liquid phase could possibly be realized due to the strongly suppressed Δ TCF under tensile strain. Our results show that the local Co-site symmetry and dynamical correlation effects will play an important role in engineering the novel magnetic phase in this and related materials.

Nguyen, Nam [Univ. of Illinois, Chicago, IL (Unite↗
DOE OSTI2025

Solution-Processed Temperature-Adaptive Radiative Paint as a Thermal Imaging Sensitizer

Thermography detects mid-infrared radiation from surfaces based on the Stefan-Boltzmann law, mapping surface temperatures and potentially revealing subsurface thermal activity. Recent developments offer an alternative strategy to traditional camera-based improvements: a thermal imaging sensitizer (TIS) coating whose emissivity increases sharply with local temperature, amplifying small thermal variations into high-contrast signals. Existing TIS structures are nanofabricated solid membranes and face significant challenges in fabrication complexity. Here, in this study, we present a solution-processed, liquid form of TIS, termed temperature-adaptive radiative paint (TARP), to address these limitations. TARP offers drastically reduced fabrication costs, scalability to large areas, applicability to curved surfaces, and an extended operating temperature range, while maintaining the function of temperature amplification. Application of TARP enhances small temperature contrast by more than 3 times, substantially improving ambient thermography and enabling broader applications such as detection of structural defects and hot spots in electronic components.

paint↗
DOE OSTI2026

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗
DOE OSTI2025

Morphology Effects on Free Energies of Proton-Coupled Electron Transfer in Polyoxotungstates

Polyoxotungstates have previously been established to facilitate the hydrogenation of small molecule substrates via hydrogen atom transfer from reactive hydroxyl groups formed at the assembly surface. Understanding structure−function relationships that dictate the thermochemistry and kinetics of protoncoupled electron transfer is key to controlling this chemistry. In this work, we combine comprehensive electrochemical experiments and density functional theory calculations to address how different polyoxotungstate morphologies, specifically W 6 O 19 −2 , W 10 O 32 −4 , SiW 12 O 40 −4 , and P 2 W 18 O 62 −6 , affect the bond dissociation free energies of surface hydroxides (BDFE(O−H)) formed upon reduction of the assembly in acidic media. Our results reveal increasing hydroxide bond strengths with increasing cluster size, and that anisotropic cluster geometries result in substantial thermodynamic differentiation of H-binding sites. We demonstrate an excellent agreement between theory and experiments on the reported BDFE(O−H) values and, importantly, we elucidate how cluster size and shape affect electronic properties (local charges and frontier molecular orbitals), giving rise to sites with increased preference for hydrogen binding, demonstrated in higher BDFE(O−H). Overall, this work aids the understanding and design of polyoxometalates exhibiting surface sites with tailored interaction strengths.

anions↗