Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dependence of predicted bulk properties of hexagonal hydroxyapatite on exchange-correlation functional

Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Vibrations and Phase Stability in Mixed Valence Antimony Oxide

α-Sb 2 O 4 (cervantite) and β-Sb 2 O 4 (clinocervantite) are mixed valence compounds with equal proportions of Sb III and Sb V as represented in the formula Sb III Sb V O 4 . Their structure and properties can be difficult to calculate owing to the Sb III lone-pair electrons. Here, we present a study of the lattice dynamics and vibrational properties using a combination of inelastic neutron scattering, Mössbauer spectroscopy, nuclear inelastic scattering, and density functional theory (DFT) calculations. DFT calculations that account for lone-pair electrons match the experimental densities of phonon states. Mössbauer spectroscopy reveals the β phase to be significantly harder than the α phase. Calculations with O vacancies reveal the possibility for nonstoichiometric proportions of Sb III and Sb V in both phases. An open question is what drives the stability of the α phase over the β phase, as the latter shows pronounced kinetic stability and lower symmetry despite being in the high-temperature phase. Since the vibrational entropy difference is small, it is unlikely to stabilize the α phase. In conclusion, our results suggest that the α phase is more stable only because the material is not fully stoichiometric.

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

Stability, electronic quantum states, and magnetic interactions of Er 3+ ions in Ga 2 ⁢O 3

Here, we report an ab initio study of phase stability, defect formation, electronic structure, and multiple magnetic, Dzyaloshinskii-Moriya, optical, hyperfine, and crystal field interactions in erbium (Er)-doped wide band gap 𝛼- and 𝛽-gallium oxides (Ga 2 ⁢O 3 ), critically important to make a foundation for both optoelectronic and quantum information applications. The chemical, structural, mechanical, and dynamical stabilities of the pristine phases are confirmed from respective negative formation energies, negative cohesive energies, favorable elastic constants, and positive phonon frequencies. The phonon dispersions indicate that the Ga-O bonds are uniform in the 𝛼-phase, while they vary in the 𝛽-phase due to the anisotropic polyhedral movement. The defect formation energy analysis confirms that both Er-doped 𝛼- and 𝛽−Ga 2 ⁢O 3 prefer Er 3+ (neutral) state. The underestimated band gaps of the pristine phases from standard density functional theory (DFT) calculations as compared to experimental values are corrected by employing the hybrid functional calculations, resulting in the indirect band gaps of 5.21 eV in 𝛼−Ga 2 ⁢O 3 and 4.94 eV in 𝛽−Ga 2 ⁢O 3 . The site preference energy analysis indicates partial occupation of Er in the octahedral site of Ga. The anisotropic nature of hyperfine tensor coefficients of Er is similar in both phases, which may be due to the occupation of Er in the same octahedral Ga site. On the other hand, the calculated magnetic exchange interaction between two Er dopants is negative for 𝛼 and positive for 𝛽, indicating an antiferromagnetic (AFM) ground state in the former and a ferromagnetic (FM) ground state in the latter. Large values of Dzyaloshinskii-Moriya interactions (DMIs) are obtained along the 𝑥 direction in the 𝛼 and along the 𝑦 direction in the 𝛽. The large DMI may support exotic magnetic textures, a promising direction for spintronic applications. The analysis of dielectric constants and refractive indices of both pristine and Er-doped phases shows a good agreement with available experimental values. The calculated optical anisotropy is slightly higher in 𝛽 than those in 𝛼, which is due to the involvement of lower symmetry in 𝛽. The crystal field coefficients (CFCs) calculated from DFT are used to analyze 4⁢𝑓 multiplets and 4⁢𝑓 −4⁢𝑓 transitions. Thus calculated lowest energy level of the first excited state to the lowest energy level of the ground state is about 1.53 µ⁢m, which is in a good agreement with available experiments, and it falls within the quantum telecommunication wavelength range.

3-dimensional systems↗
DOE OSTI2022

Multiphase equation of state for Ta 2 O 5

A new equation of state for Ta 2 O 5 is presented. The EOS is constructed using the OpenSesame software and is referred to as SESAME 3530. The EOS uses a combination of density functional theory (DFT) calculations and experimental data. DFT calculations include cold curves and phonons of the solid phases, as well as DFT-based molecular dynamics simulations of the liquid phase. Experimental data includes isobaric, diamond anvil cell, and porous shock Hugoniot data. To fit the data, we create a multiphase EOS consisting of two solid phases and the liquid. Overall agreement with experimental data is shown, and we provide some suggestions for future experiments that could improve our knowledge of the phase diagram.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

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

Heteroleptic Square Planar Cobalt(I/II) Complexes

Reduction of the cobalt(II) chloride complex, Ph 2 B( t BuIm) 2 Co(THF)Cl (1) in the presence of t BuN≡C affords the diamagnetic, square planar cobalt(I) complex Ph 2 B( t BuIm) 2 Co(C≡N t Bu) 2 (2). This is a rare example of a 16-electron cobalt(I) complex that is structurally related to square planar noble metal complexes. Accordingly, the electronic structure of 2, as calculated by DFT, reveals that the HOMO is largely d z 2 in character. Complex 2 is readily oxidized to its cobalt(II) congener [Ph 2 B( t BuIm) 2 Co(C=N t Bu) 2 ]BPh 4 (3-BPh 4 ), whose EPR spectral parameters are characteristic of low-spin d 7 with an unpaired electron in an orbital of d z 2 parentage. Further, this is also consistent with the results of DFT calculations. Despite its 16-electron configuration and the d z 2 parentage of the HOMO, the only tractable reactions of 2 involve one electron oxidation to afford 3.

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

Development and application of a uranium mononitride (UN) potential: Thermomechanical properties and Xe diffusion

We report that atomic-scale modeling of thermophysical and defect properties of uranium mononitride (UN) plays an important role in establishing a better understanding and improved models of UN fuel performance. Having an accurate interatomic potential is crucial for generating reliable data at finite temperatures using molecular dynamic simulations. We report a new interatomic potential for UN, based on a combination of many-body and pairwise interactions, a simple form that we later show could be easily adapted to include Xe-U and Xe-N interactions, i.e., generating a U-N-Xe interatomic potential. The potential was fitted to experimental thermal expansion and single crystal elastic constants, as well as Frenkel, Schottky, anti-Schottky, and antisite pair reaction energies from density functional theory (DFT) calculations. Using the potential, we successfully reproduced experimental lattice parameters, thermal expansion, single crystal elastic constants, and temperature dependent heat capacity. The potential also performs reasonably well in reproducing the energy of the aforementioned stoichiometric defect reactions and defect migration barriers calculated using DFT. However, the potential underestimates the energy difference between the tetrahedral and dumbbell uranium interstitials, and a more complex potential form might be needed to overcome this issue. The potential was also used to predict UN single crystal elastic constants and elastic properties at different temperatures, showing that UN becomes softer and more compressible with increasing temperature. We also compare our potential against literature data from previous empirical potentials, demonstrating similar or better behavior depending on the property of interest. To enable the simulation of Xe in UN a Buckingham potential has been fitted to DFT-derived Xe incorporation energies. The potential was then used to determine the activation energy for Xe diffusion due to various Xe-containing defects, with {XeU : VU} exhibiting the lowest activation energy.

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

Engineering defect energy landscape of CoCrFeNi high-entropy alloys by the introduction of additional dopants

The concept of high-entropy alloys (HEAs) focusing on tuning the overall chemical complexity represents a novel alloy design strategy. In contrast, alloying of a metallic matrix with minor doping elements with limited and localized tunability has been a common practice to improve material performance. Combining the idea of globally engineering defect energy landscape in HEAs and the localized doping strategy in dilute alloys, in this work, we explore doping effects of minor elements in a HEA matrix to further enhance the overall and localized chemical tunability, aiming to improve its irradiation resistance. Specifically, we study the influence of minor Al, Cu, Ti, and Pd substitutional doping elements on defect energetics in a CoCrFeNi model HEA based on density-functional theory (DFT) calculations. The DFT results indicate that the formation and migration energies of vacancies can be strongly influenced when a dopant is introduced at the first nearest neighbor shells around a vacancy. On the other hand, interstitial energetics are only slightly affected. Among the four elements, Ti and Pd generally decrease vacancy formation energies and increase vacancy migration energies more significantly than Al and Cu. The doping effects become more pronounced when the concentration of the substitutional dopants increases. Based on the energy distributions obtained from DFT, we build a kinetic Monte Carlo (kMC) model to assess the impact of dopants on vacancy-mediated diffusivity in the doped HEAs. Our results suggest that Ti and Pd can lower the tracer diffusivity in the considered HEAs and act as trapping sites, whereas Cu may enhance the atomic transport. This work indicates that substitutional doping in HEAs is an effective strategy in metallurgy to further tune the defect and transport properties of complex alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

A Model for Atomic Precision p-Type Doping with Diborane on Si(100)-2×1

Diborane (B 2 H 6 ) is a promising molecular precursor for atomic precision p-type doping of silicon that has recently been experimentally demonstrated. In this study, we use density functional theory (DFT) calculations to determine the reaction pathway for diborane dissociating into a species that will incorporate as electrically active substitutional boron after adsorbing onto the Si(100)-2×1 surface. Our calculations indicate that diborane must overcome an energy barrier to adsorb, explaining the experimentally observed low sticking coefficient (<1 × 10 –4 at room temperature) and suggesting that heating can be used to increase the adsorption rate. Upon sticking, diborane has an ≈50% chance of splitting into two BH 3 fragments versus merely losing hydrogen to form a dimer such as B 2 H 4 . As boron dimers are likely electrically inactive, whether this latter reaction occurs is shown to be predictive of the incorporation rate. The dissociation process proceeds with significant energy barriers, necessitating the use of high temperatures for incorporation. Using the barriers calculated from DFT, we parameterize a Kinetic Monte Carlo model that predicts the incorporation statistics of boron as a function of the initial depassivation geometry, dose, and anneal temperature. Our results suggest that the dimer nature of diborane inherently limits its doping density as an acceptor precursor and furthermore that heating the boron dimers to split before exposure to silicon can lead to poor selectivity on hydrogen and halogen resists. This suggests that, while diborane works as an atomic precision acceptor precursor, other non-dimerized acceptor precursors may lead to higher incorporation rates at lower temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2021

A Uranium(II) Arene Complex That Acts as a Uranium(I) Synthon

We report two-electron reduction of the amidate-supported U(III) mono(arene) complex U(TDA) 3 (2) with KC 8 yields the anionic bis(arene) complex [K[2.2.2]cryptand][U(TDA) 2 ] (3) (TDA = N-(2,6-di-isopropylphenyl)pivalamido). EPR spectroscopy, magnetic susceptibility measurements, and calculations using DFT as well as multireference CASSCF methods all provide strong evidence that the electronic structure of 3 is best represented as a 5f 4 U(II) metal center bound to a monoreduced arene ligand. Reactivity studies show 3 reacts as a U(I) synthon by behaving as a two-electron reductant toward I 2 to form the dinuclear U(III)–U(III) triiodide species [K[2.2.2]cryptand][(UI(TDA) 2 ) 2 (μ-I)] (6) and as a three-electron reductant toward cycloheptatriene (CHT) to form the U(IV) complex [K[2.2.2]cryptand][U(η 7 -C 7 H 7 )(TDA) 2 (THF)] (7). The reaction of 3 with cyclooctatetraene (COT) generates a mixture of the U(III) anion [K[2.2.2]cryptand][U(TDA) 4 ] (1-crypt) and U(COT) 2 , while the addition of COT to complex 2 instead yields the dinuclear U(IV)–U(IV) inverse sandwich complex [U(TDA) 3 ] 2 (μ-η 8 :η 3 -C 8 H 8 ) (8). Two-electron reduction of the homoleptic Th(IV) amidate complex Th(TDA) 4 (4) with KC 8 gives the mono(arene) complex [K[2.2.2]cryptand][Th(TDA) 3 (THF)] (5). The C–C bond lengths and torsion angles in the bound arene of 5 suggest a direduced arene bound to a Th(IV) metal center; this conclusion is supported by DFT calculations.

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

Spectroscopic signature of negative electronic compressibility from the Ti core-level of titanium carbonitride MXene

Two-dimensional transition metal carbides, carbonitrides, and nitrides called MXenes exhibit high metallic conductivity, ion intercalation capability and reversible redox activity, prompting their applications in energy storage and conversion, electromagnetic interference (EMI) shielding, and electronics, among many other fields. It has been shown that replacement of about 50% of carbon atoms in the most popular MXene family member, titanium carbide (Ti 3 C 2 T x ), by nitrogen atoms, forming titanium carbonitride (Ti 3 CNT x ), leads to drastically different properties, such as very high negative charge in solution and extreme EMI shielding effectiveness, exceeding all known materials, even metals at comparable thicknesses. Here, by using ultraviolet photoemission spectroscopy (UPS), the electronic structures of Ti 3 CNT x and Ti 3 C 2 T x are systematically investigated and compared as a function of charge carrier density. We observe that, in contrast to Ti 3 C 2 T x , the Ti 3p core-level of Ti 3 CNT x exhibits a counterintuitive shift to a lower binding energy of up to approximately 250 meV upon increasing the electron density, which is a spectroscopic signature of negative electronic compressibility (NEC). These experimentally measured chemical potential shifts are well-captured by the density functional theory (DFT) calculation. The DFT results also further suggest that the hybridization of titanium-nitrogen bonding in Ti 3 CNT x helps promoting the available states of Ti atoms for receiving more electron above the Fermi level and leads to the observed NEC. Furthermore, our findings explain the differences in electronic properties between the two very important and widely studied MXenes and also suggest a new strategy to apply the NEC effect of Ti 3 CNT x in energy and charge storage applications.

2D materials↗
DOE OSTI2020

Magnetotransport and electronic structure of the antiferromagnetic semimetal YbAs

A number of rare-earth monopnictides have topologically nontrivial band structures together with magnetism and strong electronic correlations. In order to examine whether the antiferromagnetic (AFM) semimetal YbAs (TN=0.5 K) exhibits such a scenario, we have grown high-quality single crystals using a flux method, and characterized the magnetic properties and electronic structure using specific heat, magnetotransport, and angle-resolved photoemission spectroscopy (ARPES) measurements, together with density functional theory (DFT) calculations. Both ARPES and DFT calculations find no evidence for band inversions in YbAs, indicating a topologically trivial electronic structure. From low-temperature magnetotransport measurements, we map the field-temperature phase diagram, where we find the presence of a field stabilized phase distinct from the AFM phase at low temperatures. Additionally, an extremely large magnetoresistance (XMR) for both YbAs and the nonmagnetic counterpart LuAs is also observed, which can consistently be accounted for by the presence of electron-hole compensation. Moreover, an angle-dependent study of the Shubnikov-de Haas effect oscillations reveals very similar Fermi surfaces between YbAs and LuAs, with light effective masses down to at least 0.5 K, indicating that the Yb-4f electrons are well localized, and do not contribute to the Fermi surface. However, the influence of the localized Yb-4f electrons on the magnetotransport of YbAs can be discerned from the distinct temperature dependence of the XMR compared to that of LuAs, which we attribute to the influence of short-ranged spin correlations that appear well above TN.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Low-energy interband transition in the infrared response of the correlated metal SrVO 3 in the ultraclean limit

We studied the low-energy electronic response of the prototypical correlated metal SrVO 3 in the ultraclean and disordered limit using infrared spectroscopy and density functional theory plus dynamical mean field theory calculations (DFT+DMFT). A strong optical excitation at 70 meV is observed in the optical response of the ultraclean samples but is hidden by the low-energy Drude-like response from intraband excitations in the more disordered samples. DFT+DMFT calculations reveal that this optical excitation originates from interband transitions between the bands split by orbital off-diagonal hopping, which has often been ignored in cubic systems, such as SrVO 3 . A memory function analysis of the optical data shows that this interband transition can lead to deviations of optical self-energy from the expected Fermi-liquid behavior. Our findings demonstrate that analysis schemes employed to extract many-body effects from optical spectra may be oversimplified to study the true electronic ground state and that improvements in material quality can guide efforts to refine theoretical approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Revealing the electronic structure of van der Waals antiferromagnetic NiPS 3 through synchrotron-based 𝜇-ARPES and alkali metal dosing

Antiferromagnetic NiPS 3 has recently emerged as a quantum material of considerable interest, thanks to the discovery of multiple new couplings involving electrons, spins, orbitals, phonons, and magnons. However, controversies and open questions persist concerning the fundamental origins of these couplings. A critical piece of information required to advance the understanding is the precise electronic band structure of NiPS 3 . Angle-resolved photoemission spectroscopy (ARPES), combined with alkali metal dosing (AMD), can enable us to directly observe the subtle electronic states that appear around the Fermi surface, offering valuable insights into the intriguing quantum properties and interplays of the examined material. Here, in this study, we present a comprehensive characterization and analysis of the band structure of van der Waals layered antiferromagnet NiPS 3 , leveraging state-of-the-art μ-ARPES measurements supported by density functional theory (DFT) calculations. Theoretical DFT results identify the orbital contributions to the observed bands, providing a precise understanding of the experimental ARPES data. Crucially, AMD enables the observation of conduction band and defect-related states above the valence band maximum in NiPS 3 . Furthermore, temperature dependent ARPES results across the Néel transition temperature of NiPS 3 reveal that the paramagnetic and antiferromagnetic phases have nearly identical band structures, underlining the highly localized character of Ni d states. These findings substantially deepen our understanding of the electronic properties of NiPS 3 and lay a vital foundation for exploring the intriguing quantum phenomena it exhibits.

Cao, Yifeng [Boston Univ., MA (United States); Law↗
DOE OSTI2026

Machine learning guided prediction of solute segregation at coherent and semi-coherent metal/oxide interfaces

Investigation of semi-coherent metal/oxide interfaces with misfit dislocations using density functional theory (DFT) is computationally intensive to the point of being prohibitive, as it involves several hundreds to many thousands of atoms. In this study, we examined the solute segregation behavior at the Fe/Y 2 O 3 interface—a model interface for cladding applications in nuclear fission reactors—using a combination of DFT calculations and machine learning (ML) approaches. Both coherent and semi-coherent interfaces were considered. ML models were trained on DFT-calculated segregation energies to identify the key chemical, geometric and strain energy related features that govern solute segregation behavior at coherent Fe/Y 2 O 3 interfaces. Furthermore, it was found that ML models when trained on DFT calculated segregation energy of elements at a coherent interface, comprising of about a hundred-atom supercell, can predict the segregation energy of elements at a semi-coherent Fe/Y 2 O 3 interface (with multiple hundreds of atoms) at a fraction of computational cost (1/35th), with an accuracy comparable to DFT calculations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Prediction of Solute Segregation at Metal/Oxide Interfaces Using Machine Learning Approaches

The atomic structure and chemistry at metal/oxide interfaces play a crucial role in determining their properties. However, studying semi-coherent metal/oxide interfaces that include misfit dislocations through density functional theory (DFT) is often computationally expensive due to the large number of atoms involved, ranging from hundreds to thousands. In this study, we explore solute segregation behavior at the Fe/Y 2 O 3 interface—an important model interface for cladding applications in nuclear fission reactors—by combining DFT calculations with a machine learning (ML) approach. ML models are trained using DFT-calculated segregation energies (𝐸 𝑆𝑒𝑔 ) to identify the key chemical and geometric factors influencing solute segregation at metal/oxide interfaces, revealing the competition between these features in determining 𝐸 𝑆𝑒𝑔 . Moreover, the segregation behavior at a specific Fe/Y 2 O 3 interface is predicted with high accuracy using ML models trained on data from this interface. Furthermore, it is found that the ML models could also predict solute segregation at a different Fe/Y 2 O 3 interface with a new orientation relationship (OR), at a computational cost of less than 1/45 of that required for similar DFT calculations.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Selective butene formation in direct ethanol-to-C3+-olefin valorization over Zn-Y/Beta and single-atom alloy composite catalysts using in situ generated hydrogen

The selective production of C3+ olefins from renewable feedstocks, especially via C1 and C2 platform chemicals, is a critical challenge for obtaining economically viable low-carbon middle distillate transportation fuels (i.e., jet and diesel). Here, we report a multifunctional catalyst system composed of Zn-Y/Beta and “single-atom” alloy (SAA) Pt-Cu/Al2O3 which selectively catalyzes ethanol-to-olefins (C3+, ETO) valorization in the absence of cofed hydrogen, forming butenes as the primary olefin products. Beta zeolites containing predominately isolated Zn and Y metal sites catalyze ethanol upgrading steps (588 K, 3.1 kPa ethanol, ambient pressure) regardless of cofed hydrogen partial pressure (0-98.3 kPa H2), forming butadiene as the primary product (60% selectivity at 87% conversion). The Zn-Y/Beta catalyst possesses site-isolated Zn and Y Lewis acid sites (at ~7 wt% Y) and Brønsted acidic Y sites, the latter of which has been previously uncharacterized. A secondary bed of SAA Pt-Cu/Al2O3 selectively hydrogenates butadiene to butene isomers at a consistent reaction temperature using hydrogen generated in situ from ethanol-to-butadiene (ETB) conversion. This unique hydrogenation reactivity at near-stoichiometric hydrogen and butadiene partial pressures is not observed over monometallic Pt or Cu catalysts, highlighting these operating conditions as a critical SAA catalyst application area for conjugated diene selective hydrogenation at high reaction temperatures (>573 K) and low H2/diene ratios (e.g., 1:1). Single-bed steady state selective hydrogenation rates, associated apparent hydrogen and butadiene reaction orders, and DFT calculations of the Horiuti-Polanyi reaction mechanisms indicate that the unique butadiene selective hydrogenation reactivity over SAA Pt-Cu/Al2O3 reflects lower hydrogen scission barriers relative to monometallic Cu surfaces and limited butene binding energies relative to monometallic Pt surfaces. DFT calculations further indicate the preferential desorption of butene isomers over SAA Pt-Cu(111) and Cu(111) surfaces while Pt(111) surface favors subsequent butene hydrogenation reactions to form butane over butene desorption events. Under operating conditions without of hydrogen cofeeding, this combination of Zn-Y/Beta and SAA Pt-Cu catalysts can selectively form butenes (65% butenes, 78% C3+ selectivity at 94% conversion) and avoid butane formation using only in situ generated hydrogen, avoiding costly hydrogen cofeeding requirements that hinder many renewable energy processes.

Cordon, Michael↗