Search NASA⌕ Search

SEARCH · Search NASA

Results for “barrier 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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

SEED: Semantic Energy Exploration and Discovery

The Bioenergy Knowledge Discovery Framework (KDF) hosts a vast repository of specialized data, yet traditional keyword-based search methods often struggle to provide direct answers, requiring significant domain expertise and manual effort to filter through raw documents. To overcome these barriers, this software introduces a semantic search engine that enables both specialists and non-specialists to query the KDF using natural language. By shifting from rigid keyword matching to intent-based retrieval, the tool automatically identifies and ranks the most relevant sources within the database. The system functions by processing natural language queries to extract the most pertinent information, delivering an AI-generated plain-language summary alongside exact supporting quotes from retrieved documents. This integrated approach provides users with immediate, evidence-based answers while eliminating the need for exhaustive manual review. By surfacing direct insights and contextual evidence, the software enhances the usability of existing KDF resources and democratizes access to complex bioenergy data. Ultimately, this semantic search solution accelerates the discovery process and supports faster, more informed decision-making across the bioenergy sector.

Pan, Meiyu (Melrose) [Oak Ridge National Laborator↗
DOE OSTI2025

Vacancy-Dependent Diffusion Mechanism in Oxygen-Defective SrFeO 3 Perovskite Materials: First-Principles Density Functional Theory and Experimental Approach

Understanding oxygen diffusion at the atomic scale in SrFeO 3−δ perovskites is crucial for developing oxygen storage materials with optimal performance. Such materials are required to have high stability, corrosion resistance, and acceptable oxygen storage capacity at moderate operating temperatures and pressures. Here, in this study, we used first-principles density functional theory and thermogravimetric analysis to study the vacancy-dependent oxygen diffusion in oxygen-deficient SrFeO 3−δ (δ = 0, 0.065, 0.125, 0.25, 0.5) perovskites. The electronic structures, including the partial- and spin-resolved density of states, for different SrFeO 3−δ phases were calculated and compared with available experimental and theoretical results. By mapping the migration pathways, we investigated diffusion mechanisms and calculated the energy barriers for oxygen diffusion in cubic, orthorhombic, and brownmillerite phases of SrFeO 3−δ perovskites. Using the calculated energy barriers, we deduced the diffusion time scales and diffusion coefficients within SrFeO 3−δ . A diffusion coefficient on the order of 10 –8 m 2 /s was obtained for SrFeO 2.875 . We experimentally investigated the roles of temperature and oxygen partial pressures on the redox kinetics and deduced the kinetics rate and diffusion density, which agreed well with the calculated values for the density of diffusing oxygen vacancy in the lattice. Our results showed that the energy barrier tends to reduce at higher oxygen concentrations. Our results serve as an important guideline for designing oxygen storage materials with optimal redox kinetics.

chemical looping with oxygen uncoupling (CLOU)↗
DOE OSTI2026

A framework and tool for designing cost-effective, resilient, and circular net-zero supply chains under uncertainty with an application to multilayer plastic films

While 55% of Fortune 500 companies have committed to achieving net-zero emissions and/or zero-waste operations by 2035, only 2% are currently on track, revealing a critical gap between ambition and action. Designing supply chains that reduce both emissions and waste is a complex non-intuitive, multi-objective challenge, compounded by the high costs of new technologies and the need for resilient, profitable solutions. This paper aims to address this challenge by presenting a generic framework and multi-objective optimization formulation for designing cost-effective, circular, and resilient supply chains under uncertainty, implemented through a user-friendly decision-support tool with intuitive data visualization capabilities, enabling communication of results to both technical and non-technical stakeholders. We demonstrate the application of this framework in the context of multilayer plastic films (barrier films), which are widely used in food packaging and composite materials. The model quantifies trade-offs across three objectives: minimizing global warming potential, maximizing circularity, and minimizing cost. A key contribution of this work is the explicit modeling of technological resilience, the ability of supply chains to maintain function under disruption. In the cost-minimization case, the resilience constraint makes the design approximately three times more expensive in the short-term metric, but shifts the system from relying on a single recovery pathway to a portfolio of four recovery pathways, improving the robustness of the optimization solution under uncertainty. Lastly, we introduce TranZero, a decision-support tool that integrates material flow analysis, hotspot identification, and optimization-based scenario planning to support net-zero and circularity decisions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Metal–Ligand Cooperativity via Exchange Coupling Promotes Iron- Catalyzed Electrochemical CO 2 Reduction at Low Overpotentials

Biological and heterogeneous catalysts for the electrochemical CO 2 reduction reaction (CO 2 RR) often exhibit a high degree of electronic delocalization that serves to minimize overpotential and maximize selectivity over the hydrogen evolution reaction (HER). Here, we report a molecular iron(II) system that captures this design concept in a homogeneous setting through the use of a redox non-innocent terpyridine-based pentapyridine ligand (tpyPY2Me). As a result of strong metal-ligand exchange coupling between the Fe(II) center and ligand, [Fe(tpyPY2Me)] 2+ exhibits redox behavior at potentials 640 mV more positive than the isostructural [Zn(tpyPY2Me)] 2+ analog containing the redox-inactive Zn(II) ion. This shift in redox potential is attributed to the requirement for both an open-shell metal ion and a redox non-innocent ligand. The metal-ligand cooperativity in [Fe(tpyPY2Me)] 2+ drives the electrochemical reduction of CO 2 to CO at low overpotentials with high selectivity for CO 2 RR (>90%) and turnover frequencies of 100,000 s -1 with no degradation over 20 h. The decrease in the thermodynamic barrier engendered by this coupling also enables homogeneous CO 2 reduction catalysis in water without compromising selectivity or rates. Synthesis of the two-electron reduction product, [Fe(tpyPY2Me)] o , and characterization by X-ray crystallography, Mössbauer spectroscopy, X-ray absorption spectroscopy (XAS), variable temperature NMR, and density functional theory (DFT) calculations, support assignment of an open-shell singlet electronic structure that maintains a formal Fe(II) oxidation state with a doubly reduced ligand system. Furthermore, this work provides a starting point for the design of systems that exploit metal-ligand cooperativity for electrocatalysis where the electrochemical potential of redox non-innocent ligands can be tuned through secondary metal-dependent interactions.

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

How the Bioinspired Fe 2 Mo 6 S 8 Chevrel Breaks Electrocatalytic Nitrogen Reduction Scaling Relations

The nitrogen reduction reaction (NRR) is a renewable alternative to the energy- and CO 2 -intensive Haber–Bosch NH 3 synthesis process but is severely limited by the low activity and selectivity of studied electrocatalysts. The Chevrel phase Fe 2 Mo 6 S 8 has a surface Fe–S–Mo coordination environment that mimics the nitrogenase FeMo-cofactor and was recently shown to provide state-of-the-art activity and selectivity for NRR. Here, we elucidate the previously unknown NRR mechanism on Fe 2 Mo 6 S 8 via grand-canonical density functional theory (GC-DFT) that realistically models solvated and biased surfaces. Fe sites of Fe 2 Mo 6 S 8 selectively stabilize the key *NNH intermediate via a narrow band of free-atom-like surface d-states that selectively hybridize with p-states of *NNH, which results in Fe sites breaking NRR scaling relationships. These sharp d-states arise from an Fe–S bond dissociation during N 2 adsorption that mimics the mechanism of the nitrogenase FeMo-cofactor. Furthermore, we developed a new GC-DFT-based approach for calculating transition states as a function of bias (GC-NEB) and applied it to produce a microkinetic model for NRR at Fe 2 Mo 6 S 8 that predicts high activity and selectivity, in close agreement with experiments. Furthermore, our results suggest new design principles that may identify effective NRR electrocatalysts that minimize the barriers for *N 2 protonation and *NH 3 desorption and that may be broadly applied to the rational discovery of stable, multinary electrocatalysts for other reactions where narrow bands of surface d-states can be tuned to selectively stabilize key reaction intermediates and guide selectivity toward a target product. Furthermore, our results highlight the importance of using GC-DFT and GC-NEB to accurately model electrocatalytic reactions.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Ion diffusion coefficients in poly(3-alkylthiophenes) for energy conversion and biosensing: role of side-chain length and microstructure

Conductive polymers are promising materials as active elements for energy storage and conversion devices due to mixed ion–electron conduction. The ion diffusion coefficient is a relative measure of the efficacy of ion transport, allowing for comparison between materials and electrochemical conditions. In this work, diffusion coefficients of hexafluorophosphate (PF6–) counterions in poly(3-alkylthiophene) (P3AT) materials are measured as a function of both side-chain length and microstructure using electrochemical impedance spectroscopy (EIS). For semi-crystalline films, the diffusion coefficient is found to be anomalous and nearly independent of applied electrochemical potential. Here, the anomalous behavior of diffusion indicates that spin casting yields compact films with an enthalpic barrier to ion transport, attributed to ionic trapping. Diffusion coefficient values ~10 –11 cm 2 s –1 were measured for all films, indicating interchain spacing, in the absence of strong intermolecular interactions with the electrolyte, is not a viable design strategy to control ion transport. For the prototypical system of poly(3-hexylthiophene), we observe almost no potential dependence in ion transport for regioregular and regiorandom films of comparable molecular weight, with both exhibiting anomalous diffusion. Alternatively, changing the microstructure of poly(3-hexylthiophene) to a mostly amorphous, ion-imprinted structure yields ~500× increase in the diffusion coefficient to ~2 × 10 –8 cm 2 s –1 at 0.8 V vs. Ag/Ag + with behavior closer to ordinary diffusion. Collectively, these results indicate new insight into ion transport in conductive polymers, where ionic trapping effects can be mitigated through electrodeposition protocols over post-synthesis processing (i.e. spin coating).

14 SOLAR ENERGY↗
DOE OSTI2026

Leveraging chiral induced spin selectivity to improve the efficiency and selectivity of electrocatalytic reduction

Energy production is a fundamental pillar of modern society, but it remains heavily reliant on fossil fuel combustion, which leads to environmental degradation and resource depletion. Electrocatalytic reduction reactions (ERRs) offer a cleaner alternative by functioning under mild conditions, eliminating the need for high temperatures or pressures. However, their widespread adoption is hindered by challenges such as low efficiency, inadequate selectivity, and sluggish kinetics. Here, the introduction of chiral-induced spin selectivity (CISS) offers a promising strategy to address these limitations by reducing energy barriers in reactions involving spin-sensitive intermediates or transition states. By leveraging CISS, ERR pathways can be manipulated to favor desired reactions, enhance product yield, and suppress unwanted byproducts, improving selectivity. Beyond reducing the impact of competing side reactions, CISS also enhances the kinetics of ERRs by accelerating reaction rates and lowering the energy barriers associated with multi-electron transfer steps. Thanks to its straightforward and versatile application, CISS shows significant potential in optimizing well-studied ERRs, such as the hydrogen evolution reaction (HER) and carbon dioxide reduction reaction (CO2RR), facilitating the future transition of these technologies from research laboratories to commercial energy plants worldwide.

14 SOLAR ENERGY↗
DOE OSTI2025

16 O + 16 O fusion at sub-coulomb energies

A simple optical model (OPM) method using non-monotonic (NM) potentials characterized by a repulsive core from the microscopic theory of the Pauli-led energy-density functional (EDF) has been developed to investigate the 16 O + 16 O fusion at sub-Coulomb energies relevant to the oxygen burning. Here, the study involves the analysis of the experimental fusion cross-section (FCS) data in energy range 6.92 ≤ E cm ≤ 13.83 MeV, which includes the Coulomb barrier region E cm = 10.0–11.92 MeV. Apart from an excellent description of the existing FCS data in the studied energy range, the associated hindrance, characterized by the S-factor reaching a maximum and then gradually falling off at lower energies, so far suggested empirically for the system, is reproduced down to 4 MeV for the first time in the simple OPM. An S-factor maximum of S 0 = 3.15 × 10 25 MeV.b at E 0 = 7.14 MeV is observed with T ≃ 2.6 GK, which conforms to the values reported for quiescent and explosive burning. Our reaction rate, deduced from the NM potential, compares well with the Caughlan and Fowler data. Dominant partial waves implicit in the observed maximum S-factor in the studied Gamow range are also explored. Our present findings, with the success of NM potentials, suggest that the nucleus–nucleus potential is non-monotonic.

16O+16O↗
DOE OSTI2022

The impacts of charge transfer, localization, and metallicity on hydrogen retention and transport capacity

Solid state hydrides such as early transition metal hydrides are of inestimable importance for the future of hydrogen energy and are actively being investigated for energy conversion and storage applications such as fuel cells, solid-state batteries and neutron moderators. The retention and transport behavior of hydrogen in these hydrides has a huge role on the extended performance of components. While early transition-metal-based compounds exhibit many peculiar properties due to their unique correlated electronic signatures arising from -orbital electrons, the fundamental chemistry and transport behavior of hydrogen in such hydrides is not well understood. In the present work, using density functional theory, a highly intricate bonding feature is revealed through the theoretical investigation of the electronic structure of early transition metal hydrides YH 2 and ZrH 2 . In particular, a pronounced charge transfer from the transition element to H, results in localized electron densities at deep energy levels. The interplay between intrinsic charge transfer, charge localization, and metallicity in YH 2 and ZrH 2 leads to strong chemical bonding between metal and hydrogen atoms and large energy barriers for the migration of hydrogen vacancies. Specifically, hydrogen vacancies are found to be stable in the neutral state due to electron screening effects, accompanied by substantially high migration barriers between 0.8–1.2 eV along different crystallographic directions. In contrast, recent literature shows the migration barrier for charged H vacancies in insulating s-block metal hydrides lie between 0.1–0.4 eV, which is suitable for fast conduction applications. This pivotal electron structure difference exploited between early transition metal hydrides and alkali/alkaline earth metal hydrides determines extended hydrogen retention in these early transition metal hydrides. Finally, this work explains fundamental differences between the electronic structure of s-block and d-block metal hydrides, and its impact on the mobility of hydrogen vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Revealing Robust Room Temperature Ferromagnetism in Gd‐Doped Few‐Layered MoS 2 Thin Films

2D MoS 2 holds great promise for spintronics, yet is limited by intrinsic diamagnetism. This study demonstrates inducing ferromagnetic behavior in MoS 2 films doped with 0.47% Gd, achieving an ultrahigh saturation magnetization of 454 emu/cm 3 in a few‐layered film over 11‐times higher than bulk films (40 nm). Raman spectroscopy, X‐ray photoelectron spectroscopy, X‐ray magnetic circular dichroism, and density functional theory (DFT) calculations reveal an interplay between Gd dopants and Mo, S vacancies (V 1Mo+2S ), leading to the formation of bound magnetic polarons (BMPs) that drive ferromagnetic ordering. H 2 S annealing and DFT calculations reveal that defect healing reduces the saturation magnetization by 83%. High sulfur migration barrier in few‐layered films helps preserve BMPs, thereby sustaining ferromagnetism, whereas lower migration barriers in bulk films lead to suppression. These findings highlight the synergy between Gd doping and defect engineering in achieving ultrahigh room‐temperature ferromagnetism, offering a scalable strategy for developing high‐performance 2D magnetic materials for spintronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Mechanistic Insights into Dinitrogen Reduction to Ammonia in Light-Controlled Nanocrystal:Nitrogenase Complexes

Developing systems that can efficiently capture photon energy and convert this energy into fuels and chemicals requires understanding how to assemble molecular components with diverse functions into complete systems possessing selectivity and efficiency in directing charge carriers to catalytic reactions. There are many challenges to achieving this goal. One promising approach is the development of hybrid systems that combine semiconductor nanocrystals (NCs) for light capture and enzymes as efficient catalysts. Such biohybrid systems capitalize on the tunable electronic and optical properties of NCs while leveraging the unmatched specificity and efficiency of enzymes in catalyzing chemical reactions, thereby offering opportunities to surpass the limitations of each component alone. Here, we focus on recent progress in developing a biohybrid system that combines CdS NCs for photon capture with the enzyme nitrogenase to accomplish light-driven dinitrogen (N 2 ) reduction to ammonia (NH 3 ). Integrating light-harvesting materials with biological catalysts requires a deep understanding of NC properties, protein stability, and electron transfer (ET), making it an inherently multidisciplinary problem. The reduction of N 2 to NH 3 is a challenging reaction, with a high demand in both agriculture and industrial chemical production. This reaction is intrinsically energy intensive, due to the need to activate the N≡N triple bond. The current standard industrial approach to N 2 reduction, the Haber−Bosch reaction, obtains the necessary energy input from fossil fuels, whereas biological systems capable of N 2 reduction utilize the hydrolysis of ATP as their energy source. Replacing these costly, energy-intensive inputs with renewable light energy represents a critical step toward sustainable NH 3 production. Recent progress has demonstrated that semiconductor CdS NCs can be coupled to the catalytic component of nitrogenase, the MoFe protein, to form a biohybrid CdS NC:MoFe protein complex, enabling light-driven N 2 reduction rather than energy input from fossil fuels or ATP. This illustrates how inorganic NCs can functionally replace the natural Fe protein partner, yielding a biohybrid catalyst that enables controlled electron delivery and provides not only light-driven NH 3 production but also new approaches for probing enzyme catalytic function. The CdS NC:MoFe protein biohybrid system enables light-initiated electron delivery at ambient temperature, as well as temperatures below freezing, allowing for stabilization and spectroscopic characterization of key reaction intermediates. These findings highlight how photochemical biohybrids can serve as both functional catalysts and mechanistic probes. Beyond studies of the nitrogenase mechanism, studies of the CdS NC:MoFe system reveal how variables such as NC size, electrostatic binding interactions, and sacrificial electron donors (SEDs) govern complex stability, charge transfer efficiency, and catalytic performance. In addition, studies of nitrogenase and the high activation barrier for N 2 reduction are enabling investigations of new and interesting questions regarding the properties and limitations of NC biocatalysis. In this Account, we describe the key features of CdS NC:MoFe protein biohybrids and the parameters for optimal light-driven N 2 reduction, and how controlling ET with light illuminates the path to new insights into the nitrogenase mechanism.

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

Computational Design of Metal-Exchanged MFI Zeolites for Coupling CO 2 –Ethylene to Form Acrylic Acid

The production of acrylic acid by coupling CO 2 and ethylene is highly desirable, as it uses a C 1 molecule to introduce the carboxylic group into a large commodity chemical, forming an unsaturated carboxylic acid. Molecular complexes have been used as catalysts in the literature for this reaction; however, these homogeneous studies were all limited by low activity and stability, and efficient heterogeneous catalysts remain to be explored. Here, we report density functional theory calculations, through which we investigated metal-exchanged MFI zeolites with metals from 3d, 4d, and 5d series for this reaction. Screening of reaction thermodynamics has shown that Sc, Y, Cr, Mo, W, Mn, Tc, and Re dispersed into MFI zeolites are potential catalysts. We further calculated the activation barriers for the β–H transfer, a key elementary step that normally determines the overall rates, and found the lowest barrier over Sc- and Y-MFI; such a feasible kinetics is attributed to the low electronegativity of Sc and Y that strongly affects the nature of the transition states and the charge redistribution during the β–H transfer. In conclusion, this work thus provides guidelines for designing a zeolite-based catalyst for this valuable reaction and identifies promising metal centers to be used.

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

Inferring relative dose-dependent color center populations in proton irradiated thoria single crystals using optical spectroscopy

In this work, we have utilized photoluminescence spectroscopy and optical ellipsometry to characterize the dose-dependence of the photoluminescence emission intensity and changes in optical absorption of thoria single crystals subject to irradiation with energetic protons at room- and elevated-temperatures. The photoluminescence peaks and the optical absorption bands are attributed to creation of new electronic states emerging from defects resulting from displacement damage. These bands are most likely associated with electrons trapped at the oxygen vacancy sites similar to color centers formed in other irradiated oxides and halides. Our experimental observations are supported by a standard density functional theory calculation of the electronic structure in pristine and oxygen vacancy-bearing thoria crystals. The dose-dependence of the intensity of the photoluminescence peaks is used to parameterize a rate theory model that estimates the concentration of color centers in the irradiated crystals. This parameterization provides optimized migration barrier parameters for oxygen interstitials and vacancies that simultaneously capture the optical response of the crystals irradiated at room- and elevated-temperature. These optical spectroscopy techniques offer a promising pathway to characterize the population of color centers formed at the sites of oxygen anion vacancies, particularly in irradiated nuclear fuels, where atomic-level defects cannot be readily imaged using electron microscopy. When combined with other direct and indirect characterization tools, our approach can provide new insight into defect formation and accumulation in energy materials over single atomic to extended length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, deep-ultraviolet optoelectronics, and other applications due to its ultra-wide bandgap and low projected fabrication cost of large-size and high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles due to deep character of acceptor levels and polaron transport of resulting holes. However, successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, and free carrier concentration, providing opportunities for predictive process simulation. We use first-principles defect theory and defect equilibrium calculations to simulate a 3-step growth-annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen-oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. This process is similar to the Mg acceptor activation by H removal in GaN. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. The resulting free hole concentration in quenched samples is very low (109 – 1012 cm-3) due to deep energy level of these Mg acceptors, but this type converted Ga2O3 material can create a significant built-in field in a p-n junction with an adjoining n-type material. Additionally, the electron concentration is greatly suppressed in such acceptor-doped Ga2O3, which could enable the use as a current blocking layer to fabricate normally-off (enhancement-mode) vertical Ga2O3 based metal-oxide-semiconductor field effect transistors (MOSFETs), and as a guard ring for edge termination in MOSFETs and Schottky barrier diodes (SBDs) with increased breakdown voltage.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗
DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

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

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

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

Electronegativity-Guided Site Differentiation in High-Entropy Alloy for pH-universal Hydrogen Evolution Reactions

Enhancing the intrinsic activity of transition metal catalysts for the hydrogen evolution reaction (HER) remains a critical challenge in sustainable energy conversion. Herein, we report an electronegativity-guided site differentiation strategy in a single-phase CoNiCuMoW high-entropy alloy (HEA) via electrodeposition by incorporating high-electronegativity 4d/5d orbital transition metals (Mo, W) into the face-centered cubic (fcc) matrix (CoNiCu). The as-synthesized HEA demonstrates exceptional HER performance in all pH conditions, delivering an outstanding overpotential of 65 mV (alkaline), 28 mV (acidic), and 155 mV (neutral) at a current density of 100 mA cm−2, showing performance comparable to commercial Pt/C and has excellent long-term stability at high current density (1 A cm−2, 1000 h). X-Ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations reveal that the incorporation of Mo/W simultaneously alters the local coordination environment and induces element-dependent charge redistribution, accompanied by a system-level d-band center downshift, thereby optimizing the hydrogen binding strength across multimetallic sites. Meanwhile, oxophilic Mo/W sites lower the water dissociation energy barrier. These synergistic effects collectively enable efficient and durable pH-universal HER performance.

Wu, Yutong↗