Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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

Effect of synthesis temperature on the phase formation of NiTiAlFeCr compositionally complex alloy thin films

Here, the synthesis temperature dependent phase formation of Ni 10 Ti 10 Al 25 Fe 35 Cr 20 thin films is compared to a bulk processed sample of identical composition. The as-cast alloy exhibits a dual-phase microstructure which is composed of a disordered BCC phase and AlNiTi-based B2- and/or L2 1 -ordered phase(s). Formation of the BCC phase as well as an ordered AlNi-based B2 phase is observed for a thin film synthesised at 500 °C (ratio of synthesis temperature of thin film to melting temperature of bulk alloy: T/T m = 0.49), which is attributed to both surface and bulk diffusion mediated growth. Post deposition annealing at 900 °C (T/T m = 0.75) of a thin film deposited without intentional heating results in the formation of NiAlTi-based B2 and/or L2 1 -phase(s) similar to the bulk sample, which is attributed to bulk diffusion. Depositions conducted at room temperature without intentional substrate heating (T/T m = 0.20) resulted in the formation of an X-ray amorphous phase, while a substrate temperature increase to 175 °C (T/T m = 0.28) causes the formation of a BCC phase. Atom probe tomography of the thin films deposited without intentional substrate heating and at 175 °C indicates the formation of ~5 nm and ~10 nm FeAl-rich domains, respectively. This can be rationalized based on the activation energy for surface diffusion, as Ti and Ni exhibt 2.5 to 4 times larger activation energy barriers than Al, Fe and Cr. It is evident from the homologous temperature that the phase formation observed at 500 °C (T/T m = 0.49) is a result of both surface and bulk diffusion. As the temperature is reduced, the formation of FeAl-rich domains can be understood based on the differences in activation energy for surface diffusion and is consistent with kinetically limited thin film growth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

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

Charge-density based convolutional neural networks for stacking fault energy prediction in concentrated alloys

A descriptor-less machine learning (ML) model based only on charge density images extracted from density functional theory (DFT) is developed to predict stacking fault energies (SFE) in concentrated alloys. The model is based on convolutional neural networks (CNNs) as one of the promising ML techniques for dealing with complex images and data. Identification of correct descriptors is a key bottleneck to develop ML models for predicting materials properties. Often, in most ML models, textbook physical descriptors such as atomic radius, valence charge and electronegativity are used as descriptors which have limitations because these properties change in concentrated alloys when multiple elements are mixed to form a solid solution. Here, we illustrate that, within the scope of DFT, the search for descriptors can be circumvented by electronic charge density, which is the backbone of the Kohn-Sham DFT and describes the system completely. The performance of our model is demonstrated by predicting SFE of concentrated alloys with an RMSE and R 2 of 6.18 mJ/m 2 and 0.87, respectively, validating the accuracy of the proposed approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electron Transfer Theory Elucidates the Hidden Role Played by Triethylamine and Triethanolamine during Photocatalysis

Triethylamine (TEA) and triethanolamine (TEOA) are renowned, in part, for their ability to reductively quench excited states by outer-sphere electron transfer with vast and still growing applications as sacrificial electron donors for photocatalytic systems. Upon amine oxidation, the resulting TEA •+ and TEOA •+ radical cations undergo proton transfer (or hydrogen atom transfer), resulting in the formation of a chemical reductant that has an α-carbon centered radical adjacent to the nitrogen center (TEA • and TEOA • ). In this contribution, we have electrochemically and spectroscopically characterized a set of electron acceptors which, upon accepting an electron, are a series of photocatalysts, [ReCl(R 1 R 2 -bpy)(CO) 3 ] •– , where R 1 and R 2 are electron-donating and electron-withdrawing groups in the 4,4′- and 5,5′-positions on the bipyridyl ligand. Here, we substantiated the formation of the electron donors, TE(O)A • , by spin trapping electron paramagnetic resonance spectroscopy, where TE(O)A • reacts with 2,4,6-tri-tert-butylnitrosobenzene to generate N-centered and O-centered radical adducts. Having established the chemical behaviors of the electron acceptors and donors individually, the electron transfer rate constants were determined across a 1.43 V range in driving force. The redox potential of TEA • was benchmarked to within ±80 mV on an absolute scale in V vs Fc + /Fc in CH 3 CN by using an empirical rate vs free-energy correlation, electron transfer theory, and density functional theory calculations. The equilibrium potentials for TEA • and TEOA • were determined to be −1.98 V and −1.76 V, respectively. Based on the kinetic and thermochemical analysis presented for TEA • and TEOA • , these transient radicals can be broadly considered strong homogeneous chemical reductants within the wider context of photoredox potentials. Thus, this work clarifies a frequently unnoticed secondary function for these sacrificial electron donors during photocatalysis and rationalizes the possibility of a one-photon/two-electron conversion process that is dependent on the free-energy exchange between TE(O)A • and photocatalysts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Efficient method for estimation of fission fragment yields of $\textit{r}$-process nuclei

Background: More than half of all the elements heavier than iron are made by the rapid neutron capture process (or $\textit{r}$ process). For very-neutron-rich astrophysical conditions, such at those found in the tidal ejecta of neutron stars, nuclear fission determines the $\textit{r}$-process endpoint, and the fission-fragment yields shape the final abundances of 110 ≤ $\textit{A}$ ≤ 170 nuclei. The knowledge of fission-fragment yields of hundreds of nuclei inhabiting very-neutron-rich regions of the nuclear landscape is thus crucial for the modeling of heavy-element nucleosynthesis. Purpose: In this study, we propose a model for the fast calculation of fission-fragment yields based on the concept of shell-stabilized prefragments defined with help of the nucleonic localization functions. Methods: To generate realistic potential-energy surfaces and nucleonic localizations, we apply Skyrme density-functional theory. In this work, the distribution of the neck nucleons among the two prefragments is obtained by means of a statistical model. Results: We benchmark the method by studying the fission yields of 178 Pt, 240 Pu, 254 Cf, and 254,256,258 Fm and show that it satisfactorily explains the experimental data. We then make predictions for 254 Pu and 290 Fm as two representative cases of fissioning nuclei that are expected to significantly contribute during the $\textit{r}$-process nucleosynthesis occurring in neutron-star mergers. Conclusions: The proposed framework provides an efficient alternative to microscopic approaches based on the evolution of the system in a space of collective coordinates all the way to scission. It can be used to carry out global calculations of fission-fragment distributions across the $\textit{r}$-process region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Accelerating the development of new solar absorbers by photoemission characterization coupled with density functional theory

Abstract The expectation to progress towards Terawatts production by solar technologies requires continuous development of new materials to improve efficiency and lower the cost of devices beyond what is currently available at industrial level. At the same time, the turnaround time to make the investment worthwhile is progressively shrinking. Whereas traditional absorbers have developed in a timeframe spanning decades, there is an expectation that emerging materials will be converted into industrially relevant reality in a much shorter timeframe. Thus, it becomes necessary to develop new approaches and techniques that could accelerate decision-making steps on whether further research on a material is worth pursuing or not. In this review, we will provide an overview of the photoemission characterization methods and theoretical approaches that have been developed in the past decades to accelerate the transfer of emerging solar absorbers into efficient devices.

14 SOLAR ENERGY↗
DOE OSTI2021

Role of quadrupole deformation and continuum effects in the “island of inversion” nuclei 28,29,31 F

The peculiar properties of nuclei in the so-called “island of inversion” around Z = 10 and N = 20 are the focus of current nuclear physics research. Recent studies showed that 28 F has a negative-parity ground state and thus lies within the southern shore of the island of inversion, and 29 F presents a halo structure in its ground state, but it is unclear which effects, such as deformation, shell evolution due to tensor forces, or couplings to the continuum, lead to this situation. Here, we investigate the role of quadrupole deformation and continuum effects on the single-particle structure of 28,29,31 F from a relativistic mean-field approach and show how both phenomena can lead to a negative-parity ground state in 28 F and halo structures in 29,31 F. We solve the Dirac equation in the complex-momentum (Berggren) representation for a potential with quadrupole deformation at the first order obtained from relativistic mean-field calculations using the NL3 interaction and calculate the continuum level densities using the Green’s function method. We extract single-particle energies and widths from the continuum level densities to construct the Nilsson diagrams of 28,29,31 F in the continuum and analyze the evolution of both the widths and occupation probabilities of relevant Nilsson orbitals in 28 F and find that some amount of prolate deformation must be present. In addition, we calculate the density distributions for bound Nilsson orbitals near the Fermi surface in 29,31 F and reveal that, for a quadrupole deformation 0.3 ≤ β 2 ≤ 0.45 (prolate), characteristic halo tails appear at large distances. Using the relativistic mean-field approach in the complex-momentum representation with the Green’s function method, we demonstrate that in neutron-rich fluorine isotopes, while in the spherical case the pf shells are already inverted and close to the neutron emission threshold, a small amount of quadrupole deformation can dramatically reduce the gap between positive- and negative-parity states and increase the role of continuum states, ultimately leading to the negative parity in the ground state of 28 F and the halo structures in 29,31 F.

20 ≤ A ≤ 38↗
DOE OSTI2025

Atomic Energy Accuracy of Neural Network Potentials: Harnessing Pretraining and Transfer Learning

Machine learning-based interatomic potentials (MLIPs) have transformed the prediction of potential energy surfaces (PESs), achieving accuracy comparable to ab initio calculations. However, atomic energy predictions, often assumed to lack physical meaning, remain underexplored. In this study, we demonstrate that inaccuracies in atomic energy predictions reduce the robustness and transferability of Neural Network Potentials (NNPs) and atomic energy error can be masked in total energy predictions due to error cancellation. Here, we validate this finding using challenging configurations involving deformation and failure under tensile loading. By pretraining atomic energy predictions using empirical potentials and applying transfer learning with density functional theory (DFT) data, we achieve notable improvements in the accuracy of total energy, forces, and stress predictions. Furthermore, this approach enhances the robustness and transferability of NNPs, emphasizing the importance of atomic energy predictions in developing high-quality and reliable MLIPs.

Active Learning↗
DOE OSTI2024

A First-Principles Approach to Modeling Surface Site Stabilities on Multimetallic Catalysts

The study of multimetallic alloys and the multitude of possible surface compositions have sparked a tremendous interest in engineering low-cost materials with high activity and selectivity in heterogeneous catalysis. Multimetallic systems provide complementary functionalities and an unprecedented tunability when designing catalyst formulations. However, due to their immense structural and compositional complexity, the investigation and identification of an optimal catalyst is a tedious and time-consuming process, both experimentally and theoretically. Therefore, theoretical design principles are highly desirable to accelerate the screening of catalyst structures across the vast compositional space. In this paper, we introduce a simple and general model for predicting the site stability of multimetallic surfaces and nanoparticles, which is based on physical principles. The model requires only a small set of density functional theory (DFT) calculations of metal atom binding energies on monometallic and dilute alloy surface slabs to optimize the parameters in the simple model. The resulting model allows for the quantification of the stability of any particular atom site in any conceivable chemical environment across a wide range of morphologies, sizes, and arrangements by interpolating the derived parameters from a monometallic system to a completely diluted alloyed system. Herein, we demonstrate the robustness of the model across an extensive data set of transition metal alloy surfaces and 147-atoms cuboctahedral nanoparticles (NPs) composed of IrRhRu and PtPdRu. In conclusion, our approach yields mean absolute errors of ≈0.15 (IrRhRu), 0.20 (PtPdRu), 0.19 (IrRhRu NP), and 0.26 (PtPdRu NP) eV relative to site binding energies calculated using DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atoms and interatomic bonding synergism inherent in molecular electronic wave functions

The chemical model of matter consists of atoms held together by bonds. The success of this model implies that the physical interactions of the electrons and nuclei in molecules combine into compound interactions that create the bonding. In the quantum mechanical description, the modified atoms in molecules and the bonding synergism are contained in the molecular electronic wave function. So far, only part of this information has been recovered from the wave function. Notably, the atoms have remained unidentified in the wave function. One reason is that conventional energy decomposition analyses formulate separate model wave functions, independent of the actual wave function, to represent “prepared atoms” and preconceived interactions and, then, intuitively catenate the parts. In the present work, the embedded modified atoms and the inherent physical synergisms between them are determined by a unified derivation entirely from the actual molecular valence space wave function. By means of a series of intrinsic orbital and configurational transformations of the wave function, the energy of formation of a molecule is additively resolved in terms of intra-atomic energy changes, interference energies, and quasi-classical, non-classical, and charge-transfer Coulombic interactions. Here, the analysis furnishes an algorithm for the quantitative resolution of the energy of formation, which enables analyses elucidating reaction energies.

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

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso–meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4- phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi’s golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

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

Atomistic Insight into Cation Effects on Binding Energies in Cu-Catalyzed Carbon Dioxide Reduction

Electrolyte properties in general and electrolyte ions in particular have been shown to have significant effects on the kinetics of electrochemical CO 2 and CO reduction at metal electrodes, but these effects have not yet been fully understood. Here, we investigate the effects of cations at the aqueous electrolyte–Cu(211) interface on adsorbate binding energies and the electrolyte structure using density functional theory (DFT). Charging the interface via explicit Na + has systematic effects on adsorbate–electrolyte interactions and conformations. We describe specific local adsorbate–ion interactions, including direct alkali ion–adsorbate coordination and hydrogen bonding via ion-coordinated water molecules. The relative importance of these specific interactions and purely electrostatic field–adsorbate interactions is investigated by comparing the DFT-calculated ion effects to those predicted by purely electrostatic models of the interface and the adsorbates. We find that the trend in ion effects among different adsorbates at constant surface charge density is well explained by a purely electrostatic interaction model. The binding energy of OCCHO is found to depend strongly on the surface charge density as well as the spatial distribution of ions at constant surface charge density. These effects are also explained by a purely electrostatic local electrostatic field–adsorbate model. This indicates that alkali ion effects can be mainly attributed to purely electrostatic field interactions and that the local field at ion-stabilized active sites can depend significantly on both the overall charge density and on the spatial distribution of ions at constant charge density. This work provides new insight on alkali ion effects on a variety of adsorbates relevant to the CO 2 and CO reduction reactions by describing specific local ion–adsorbate interactions, systematic changes in adsorbate–electrolyte interactions induced by explicit surface charge, and comparisons of the calculated effects on binding energies to simple electrostatic field–adsorbate models to explain trends among various adsorbates and at varying surface charge densities and supercell configurations.

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

Impact of symmetry breaking and spin-orbit coupling on the band gap of halide perovskites

Halide perovskite (HP) materials have recently emerged as a class of semiconductors with immense promise for various optoelectronic applications, ranging from solar cells to light-emitting diodes. One of the unique attributes of HPs is their tunable band gaps with different factors governing their value. Furthermore, the first factor is related to relativistic corrections [“mass-Darwin,” connected to the ns 2 lone pairs, and spin-orbit coupling (SOC)] that induce an orbital shift or degeneracy splitting, resulting in a band-gap reduction. The second factor involves the structural configuration: in HPs the local symmetry of each Wyckoff position tends to be broken, inducing an opening of the band gap. Based on high-throughput density functional theory calculations, this paper systematically studies a possible self-cancelation on the band-gap correction for HPs when the polymorphous configuration—structural effects—and the SOC—electronic effects—are included. Our results indicate that the nature of interplay between SOC and symmetry breaking (SB) is that they are independent decoupling effects to describe the band-gap magnitude in halide perovskites. As a result of that, we observe a transitivity of the band-gap description; i.e., if we know the band gap of halide perovskites without SB and SOC, we can independently add the effects of band-gap reduction due to SOC and band-gap opening due to SB, regardless of the order in which these effects are considered.

14 SOLAR ENERGY↗
DOE OSTI2023

Exchange correlation potentials from full configuration interaction in a Slater orbital basis

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. As a result, examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

CHARGED POINT DEFECTS AND DEFECT COMPLEXES IN PHOSPHORENE AND STRUCTURAL PREDICTION OF TAU-11, A DENSITY FUNCTIONAL THEORY APPROACH

With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase Stability, Band Gap Tuning, and Rashba Splitting in Selenium-Alloyed Bournonite: CuPbSb(S 1– x Se x ) 3

Recently, bournonite (CuPbSbS 3 ) has been identified as a potential ferroelectric photovoltaic (PV) material with a 1.3 eV band gap, which falls in an appropriate range for single junction PV devices. Progress in applying bournonite has yielded several studies reporting successful thin-film processing, the most recent of which demonstrated a 2.65% power conversion efficiency for a bournonite-based PV device. In an effort to explore bournonite band gap engineering for PV and other prospective applications (e.g., thermoelectric, spintronic), we here consider the solid-state synthesis of selenium-alloyed bournonite, CuPbSb(S 1-x Se x ) 3 , across the full range of x (0.0 ≤ x ≤ 1.0), and report phase purity for 0.0 ≤ x ≤ 0.5. We characterize the crystal structure and band gap of the samples using X-ray diffraction and diffuse reflectance spectroscopy, determining a band gap decrease for single-phase samples from 1.25 eV at x = 0.0 to 1.06 eV at x = 0.5. Formation energies determined using dispersion corrected hybrid density functional theory (DFT) support experimental findings and are consistent with the stability (instability) of the x = 0.5 (x = 1.0) structure relative to starting materials. The computed band structures show a decreasing trend in the band gap with increasing x, again consistent with experiment, with the states at the conduction (valence) band minima (maxima) being principally derived from Pb (S/Se) states. Spin texture analysis and detailed comparison of spin splitting parameters show that Se-alloying modifies the band gap while maintaining Rashba spin splitting character within the electronic structure. Rashba splitting is most pronounced for the conduction band minimum along the Γ-X path in the Brillouin zone. Energy separation between spin states is maintained at ΔE~0.2 eV for the various Se contents, with σz as the spin direction. Furthermore, a smaller spin splitting occurs along Γ-Z (decreasing value with increasing x), with σ x as the spin direction.

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

Critical Assessment of Electronic Structure Descriptors for Predicting Perovskite Catalytic Properties

The discovery and design of materials which can efficiently catalyze the oxygen reduction and evolution reactions at reduced temperatures is important for facilitating the widespread adoption of fuel cell and electrolyzer technologies. Numerous studies have produced correlations between catalytic properties, such as oxygen surface exchange or electrode area specific resistance (ASR), and properties of the catalyst material. However, correlations have historically been limited in scope (e.g., using only a few materials or at a single temperature) and it has been difficult to provide detailed assessments of their robustness. Here, in this study, we assess the ability of the O p-band center electronic structure descriptor, obtained from density functional theory (DFT) calculations, to correlate with oxygen surface exchange rates, diffusivities, and area specific resistances for a large database of perovskite oxide catalytic properties. By data mining the literature, we obtain 747 catalytic property value data points spanning 299 unique perovskite compositions from 313 studies. We assess linear correlations of each property with the O p-band center and find generally modest correlations that are qualitatively useful (prediction mean absolute errors of about 0.5 log units are typical), where the correlations are improved at higher temperatures (e.g., 800 °C vs. 500 °C) and significantly improve when considering fits to the subset of materials which have multiple independent measurements. These findings suggest that the spread of property data is significantly influenced by experimental uncertainty, and subsequent measurements of additional materials will likely improve the O p-band center correlations.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Accurate Electron-Phonon Interactions from Advanced Density Functional Theory

Electron-phonon coupling (EPC) is key to understanding charge transport, band renormalization, and superconductivity in energy materials, including correlated transition-metal oxides, ferroelectric perovskites, optoelectronic semiconductors, and phonon-mediated superconductors. Although first-principles density functional theory (DFT)-based EPC calculations are used widely, their predictive power is limited by the accuracy, transferability, and efficiency of the underlying exchange-correlation functionals. These limitations become exacerbated in complex 𝑑- and 𝑓-electron materials, where beyond-DFT approaches and additional corrections, such as the Hubbard 𝑈, are commonly invoked. Here, using the examples of Co⁢O and Ni⁢O, we show how the 𝑟 2 SCAN density functional correctly captures strong EPC effects in transition-metal oxides without requiring the introduction of Hubbard 𝑈 parameters. We also find that 𝑟 2 SCAN successfully describes the subtle interplay between ionic and covalent bonding, and strong EPC effects in the low-temperature rhombohedral phase of the prototypical ferroelectric Ba⁢Ti⁢O 3 , without requiring Hubbard 𝑈 and intersite 𝑉 corrections. We further demonstrate the ability of 𝑟 2 SCAN to accurately model the EPC of the main-group semiconductor Ga⁢As and the phonon-mediated superconducting properties of Mg⁢B 2 , with reliable electronic bands and phonons. Our study establishes 𝑟 2 SCAN-based EPC as a transferable, parameter-free, and computationally efficient framework for predictive material-specific modeling of EPC in energy materials, and opens a practical route toward high-throughput screening of superconductors, thermoelectrics, optoelectronic semiconductors, and oxide electrodes and catalysts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗