Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (DFT)”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Orbital-Free Quantum Simulation Methods for Application to Warm Dense Matter (Final Technical Report)

Predictive simulations for prediction of condensed system behavior in state conditions far from ambient is increasingly crucial to DOE priorities. Warm dense matter (WDM) is the paradigm: temperature T > 1-15 eV, pressures P to 1 Mbar or greater. Experiments under such state conditions are difficult and costly. We summarize work driven by the need and opportunity to make free-energy density functional theory (DFT) as powerful a tool for ab initio simulation of matter under such extreme conditions as ground state DFT is for ordinary matter Advancing orbital-free DFT (OF-DFT) to eliminate the Kohn-Sham (KS) scaling bottleneck in such simulations is the other priority. The concurrent challenge for both goals is the intrinsic complexity of WDM. We summarize 15 years of successes and major progress on (1) free energy exchange-correlation functionals; (2) non-interacting free energy functionals (counterpart to T=0 Kohn-Sham kinetic energy density functionals); (3) rigorous results and constraints for free-energy DFT; (4) software for free energy DFT calculations in both conventional Kohn-Sham and OF-DFT form; (5) de-orbitalization of advanced orbital-dependent ground state functionals for use in OF-DFT; (6) demonstration calculations; (7) ancillary achievements (e.g. major review articles, secondary explorations motivated by primary goals).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hydrogen irradiation-driven computational surface chemistry of lithium oxide and hydroxide

We have investigated, using molecular dynamics, the surface chemistry of hydrogen incident on the amorphous and crystalline lithium oxide and lithium hydroxide surfaces upon being slowed down by a collision cascade and retained in the amorphous surface of either Li 2 O or LiOH. We looked for the bonding of H to the resident structures in the surface to understand a possible chain of chemical reactions that can lead to surface transformation upon H atom impact. Our findings, using Density-Functional Theory (DFT) trained ReaxFF force field/electronegativity equalization method potentials, stress the importance of inclusion of polarization in the dynamics of a Li–O–H system, which is also illustrated by DFT energy minimization and quantum–classical molecular dynamics using tight binding DFT. The resulting polar-covalent chemistry of the studied systems is complex and very sensitive to the instantaneous positions of all atoms as well as the ratio of concentrations of various resident atoms in the surface.

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

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tetracene Functionalized Si(111) Achieves Enhanced Solar-to-Chemical Energy Conversion via Molecular Acceptor States

The properties of semiconductor|liquid interfaces play a critical role in determining the efficiency of solar-to-hydrogen (STH) conversion. Here, we investigate how molecular functionalization of Si(111) and Si(111)|TiO 2 surfaces impacts photoelectrochemical (PEC) hydrogen production efficiency. We find that functionalization of ∼3% of the atop sites of Si(111) with either 9-anthracene (Anth) or 5-tetracene (Tet), with the remaining sites passivated by methyl groups, provides substrates with high electronic quality and low surface oxide densities, as determined by X-ray photoelectron spectroscopy (XPS) measurements. Surface photovoltage (SPV) spectroscopy shows that surfaces modified with Anth or Tet exhibit an increased photovoltage, with Tet-functionalized surfaces yielding an additional 192 meV relative to methyl-terminated Si(111), indicating improved charge separation for Si-Tet. Further improvement in onset potential was achieved by replacing a nitrogen-containing TiO 2 atomic layer deposition (ALD) precursor (TDMAT) with a precursor lacking nitrogen (TTIP), which eliminates the parasitic defect band in the TiO 2 overlayer (p-Si(111)-Tet|TTIP-TiO 2 |Pt: V OC = +0.283 ± 0.041 V vs RHE). Density functional theory (DFT) analysis demonstrates that compared with Anth-modified Si(111), the Tet-modified surface exhibits more hybridized Si(111)-Tet states closer to the silicon band edges. Mercury contact current–voltage (I–V, dark) measurements quantified the relative interfacial density of states of Si-Tet, Si-Anth and Si-Me surfaces─revealing that the interfacial state density was highest for Si-Tet. This suggests that such hybridized interfaces serve to capture better photoexcited charge, which enables facile electron transfer to molecular acceptors in solution. Altogether, the data indicate that beneficial hybrid molecular LUMO surface states interacting with the Si conduction band edge results in improved hydrogen evolution (HER) performance for p-Si devices.

Group theory↗
DOE OSTI2023

Multiple redox mechanisms for water-gas shift reaction on Fe 3 O 4 (1 1 1) surface: A density functional theory and mean-field microkinetic modeling study

In this work, we investigated redox mechanisms for the high temperature (HT) water–gas shift reaction (WGSR) on the Fe 3 O 4 (1 1 1) surface using density functional theory (DFT) calculations and mean-field microkinetic modeling (MKM). The redox pathways branched into three Langmuir-Hinshelwood processes (LH1, LH2 and LH3) and one Mars-van-Krevelen (MVK) process (in the presence of oxygen vacancy) from multidentate binding of CO and CO 2 over four available reactive sites (Fe oct2 , Fe tet1 , Fe bridge, and O1). We found that the LH1 and LH2 processes have CO adsorption at a single iron site (Fe oct2 or Fe tet1 ), while the LH3 and MVK processes have stronger chemisorption of CO or CO 2 by both the Fe oct2 site and the O 1 site. From the mechanistic study of these reaction path, we recognized that availability of O 1 sites was key to proceed to either CO oxidation by the single site (LH1, LH2) or by the dual site (LH3, MVK). We observed that the reaction energetics were significantly different in the CO oxidation steps, where the single or dual site results in markedly different apparent activation energy (73–281 kJ/mol) and reaction rates (10 -3 –10 -10 mol∙m -2 ∙s -1 ) among the four reaction mechanisms. The utilization of mean field MKM with DFT reaction energetics helps explain the experimental debates for the catalytic reaction details as well as providing a possible direction to engineer the catalyst with higher activity.

42 ENGINEERING↗
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Impact of solvation on the GW quasiparticle spectra of molecules

First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

CO and H 2 adsorption on Au-Ni bimetallic surfaces: a combined experimental and DFT theoretical study

Au-Ni bimetallic thin films were grown on refractory metal substrates. CO and H 2 adsorption on Au-Ni bimetallic surfaces have been studied by a combination of in situ polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS), temperature-programmed desorption (TPD), and density functional theory (DFT) calculations. It is found CO desorption peak shifts from 413 K on pure Ni surfaces to 293 K on the isolated Ni atoms formed by alloying with Au atoms. The sharp decrease of CO desorption temperature on Au-Ni surfaces with increasing Au coverage is caused by the change of the favored CO adsorption sites from bridge/hollow sites on pure Ni surfaces to Ni top sites on Au-Ni bimetallic surfaces. In situ PM-IRRAS shows two CO adsorption bands on Au sites at 2119 cm −1 and 2103 cm −1 on Au-Ni surfaces at 80 K, which are due to CO bound on under-coordinated Au atoms and electron negatively charged Au sites modified with nearby Ni atoms, respectively. Even with the Au-Ni surface temperature at as low as 100 K, CO adsorption induced Ni surface segregation has been observed by in situ PM-IRRAS. Furthermore, DFT calculation results discover the adsorption energy of CO on Ni top sites continues to decrease with increasing Au coverage due to the geometric ensemble effect and the lowered d-band center after Ni alloying with Au. H 2 desorption temperature decreases from 363 K on pure Ni thin films to 302 K with increasing Au coverage to 0.6 ML. A new H 2 peak appears at around 170 K on the Au-Ni surfaces with Au coverages between 0.6 ML and 0.9 ML. This new H 2 TPD peak is assigned to H 2 desorption from the totally isolated Ni sites. With Au coverage above 1.5 ML, there is no any H 2 desorption detected. Finally, the combined surface science studies and DFT calculations provide new insights into the surface structure-activity correlation of Ni-base bimetallic surface alloys.

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

Novel ceramic capacitors with ultrahigh energy density and efficiency (Final Technical Report)

Antiferroelectric ceramics are a special class of material that have shown great potential as the dielectric in electrical capacitors due to their high energy- and power-density. During each charge-discharge cycle, the ceramic undergoes transformation to a ferroelectric phase and resumes its antiferroelectric phase. The hysteresis associated with the transitions leads to a mediocre energy efficiency and service lifetime of antiferroelectric capacitors and, hence, their almost absence in commercial products. Under the support of this research project, we first formulated a universal lattice-compatibility theory that included electrostatic polarization energy along with elastic energy and thermal energy to understand the origin of the hysteresis in antiferroelectric oxides. Guided by this compatibility theory, we conducted high-throughput density functional theory (DFT) calculations to assess chemical modifiers and their effect on crystal structures of 400+ PbZrO 3 -based compositions. Down-selected compositions were experimentally validated for their suppressed hysteresis and higher energy efficiency. The verified low-hysteresis compositions were then expanded to an antiferroelectric ceramic library with nearly 500 new compositions (more than 1,500 samples) using high-throughput experiments involving ceramic synthesis and property screening. The large quantity of data generated (theory and experimental) in these tasks were processed by machine-learning techniques and identified trends were fed to the next iteration. In the end, we successfully discovered four compositions with near-zero hysteresis, yielding a world-record energy efficiency of 98.2% at an energy density of 3.0 J/cm 3 . Furthermore, our antiferroelectric ceramic capacitor reaches 79.5 million charge-discharge cycles lifetime, a factor of 80 enhancement over previous antiferroelectric ceramics with large hysteresis. These research accomplishments have not only met the milestones set in the SOPO, but also led to two patent filings, three journal publications (one of them was in Advanced Materials, impact factor 29.4), and nine oral presentations at various venues. Through the course of the project, three postdocs, four Ph.D. students, and one M.S. student were trained. In short, our project established a new methodology in searching next-generation functional ceramics on the fundamental side and discovered several high-efficiency antiferroelectric compositions for capacitors on the applied side. Once fabricated into the multilayer form for commercial applications, these ceramic capacitors can potentially enable the high temperature high power density DC-link capacitors that are critical for the next generation inverters in electric vehicles. The project also significantly contributed to the nation’s workforce development in the STEM fields.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Design and fabrication of TiO 2 -based dye sensitized solar cells using plant-derived organic dyes

The utilization of solar energy, an abundant and renewable resource, shows great promise. Various technologies have been employed to capture this vast energy potential. In this study, the effectiveness of natural dye-sensitized solar cells (DSSCs) was explored through experimental and computational methods. Organic dyes derived from specific plant species were examined, with a solvent mixture of ethanol, methanol, and tetrahydrofuran used for extraction. A UV–vis spectrophotometer was utilized to measure the sensitizers’ optical characteristics. The study employed density functional theory (DFT), implemented in Gaussian 09 W software, to perform computational calculations. The B3LYP method and a basis set of 6-31G++(d,p) were chosen for optimizing the geometries and energetics of dye molecules. The DFT results indicated that the dye molecules had a bandgap in the range of 2.16–2.38 eV. The photochemical performance of all four fabricated DSSC-based devices was evaluated to be between 0.33% and 1.04%, and the Euryops pectinatus sensitizer demonstrated the highest efficiency of 1.04% among the studied sensitizers.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning prediction of the mechanical properties of refractory multicomponent alloys based on a dataset of phase and first principles simulation

In this work, a dataset including structural and mechanical properties of refractory multicomponent alloys was developed by fusing computations of phase diagram (CALPHAD) and density functional theory (DFT). The refractory multicomponent alloys, also named refractory complex concentrated alloys (CCAs) which contain 2–5 types of refractory elements were constructed based on Special Quasi-random Structure (SQS). The phase of alloys was predicted using CALPHAD and the mechanical property of alloys with stable and single body-centered cubic (BCC) at high temperature (over 1,500°C) was investigated using DFT-based simulation. As a result, a dataset with 393 refractory alloys and 12 features, including volume, melting temperature, density, energy, elastic constants, mechanical moduli, and hardness, were produced. To test the capability of the dataset on supporting machine learning (ML) study to investigate the property of CCAs, CALPHAD, and DFT calculations were compared with principal components analysis (PCA) technique and rule of mixture (ROM), respectively. It is demonstrated that the CALPHAD and DFT results are more in line with experimental observations for the alloy phase, structural and mechanical properties. Furthermore, the data were utilized to train a verity of ML models to predict the performance of certain CCAs with advanced mechanical properties, highlighting the usefulness of the dataset for ML technique on CCA property prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2026

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗
DOE OSTI2025

Revealing symmetry-broken superconducting configurations by density functional theory

A coherent theory for the superconductivity of both conventional and unconventional superconductors is currently lacking. Here we show that superconductivity arises from the formation of a symmetry-broken superconducting configuration (SCC) due to atomic perturbation of the normal conducting configuration (NCC). This electron–phonon interaction creates straight one-dimensional tunnels (SODTs) for charge density of electrons and/or holes as revealed by the calculations based on density functional theory (DFT). The SODTs act as resistance-free superhighways and are correlated to the Cooper pairs in the Bardeen–Cooper–Schrieffer (BCS) theory. The formation of SODTs implies that the electron–phonon interaction in the BCS theory can be represented by the difference in charge densities between SCC and NCC predicted by DFT. The present work highlights that in conventional superconductors, SODTs are embedded within the bulk materials and are easily destroyed by phonon vibrations, resulting in a low critical superconducting temperature (T C ). Conversely, in unconventional superconductors such as YBa 2 Cu 3 O 7 (YBCO 7 ), SODTs are protected by a layered pontoon structure with very weak bonding to the bulk materials, maintaining SODTs’ stability at higher temperatures and leading to a much higher T C . The present approach is validated for 14 conventional superconductors of 18 pure elements and MgB 2 examined in this work, including the presently predicted superconductivity in Cu, Ag, Au, Sb, Bi, and MgB 2 at 0 K and 0 GPa, and one unconventional superconductor of YBCO 7 . Our discovery indicates that DFT can be a practical tool for predicting superconductors, enabling a systematic search for new superconducting materials in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding the influence of solvents on the Pt-catalyzed hydrodeoxygenation of guaiacol

Bio-oil derived from fast pyrolysis of lignocellulosic biomass needs to be deoxygenated to become a substitute for petroleum fuels. Here, we study the hydrodeoxygenation mechanism of guaiacol, a bio-oil model compound derived from the lignin fraction of biomass, on Pt(1 1 1) terrace sites in the presence of water, diethyl ether, 1-butanol, and n-hexane as solvent. Using first-principles periodic density functional theory (DFT) calculations and mean-field microkinetic reactor modeling, a detailed reaction mechanism is investigated targeting various products such as catechol, phenol, anisole, benzene, cyclohexanone, and cyclohexanol. Solvent phase DFT outcomes are mostly similar to that of the vapor phase; however, microkinetic modeling results suggest that rate controlling species and transition states differ somewhat in the various reaction environments. Catechol was found to be the major aromatic product across all reaction environments. Over Pt(1 1 1), unsaturated monooxygenate production from catechol is unlikely.

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

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Machine-learning enabled thermodynamic model for the design of new rare-earth compounds

We employ a descriptor based machine-learning approach to assess the effect of chemical alloying on formation-enthalpy of rare-earth intermetallics. Application of machine-learning approaches in rare-earth intermetallic design have been sparse due to limited availability of reliable datasets. In this work, we developed an ‘in-house’ rare-earth database with more than 600 + compounds, each entry was populated with formation enthalpy and related atomic features using high-throughput density-functional theory (DFT). The SISSO (sure independence screening and sparsifying operator) based machine-learning method with meaningful atomic features was used for training and testing the formation enthalpies of rare earth compounds. The complex lattice function coupled with the machine-learning model was used to explore the effect of transition metal alloying on the energy stability of Ce based cubic Laves phases (MgCu 2 type). The SISSO predictions show good agreement with high-fidelity DFT calculations and X-ray powder diffraction measurements. Our study provides quantitative guidance for compositional considerations within a machine-learning model and discovering new metastable materials. The electronic-structure of Ce-Fe-Cu based compound was also analyzed to get an in-depth understanding of the electronic origin of phase stability. The interpretable analytical models in combination with density-functional theory and experiments provide a fast and reliable design guide for discovering technologically useful materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High temperature elastic properties of sub-stoichiometric yttrium dihydrides

Yttrium hydrides are considered as candidate materials for neutron moderation applied in microreactors (akin transportable tiny nuclear reactors) owing to their superior thermal stability and hydrogen retention. The evolution of elastic properties of these materials at elevated temperatures, needed for predicting the thermomechanical response and performance of the moderator during in-service reactor conditions, however, is lacking. Here, we report the Young’s and shear elastic moduli of three stoichiometries of bulk yttrium hydride (YH x , x = 1.61, 1.82, and 1.84) from room temperature to 1000°C. In situ temperature-dependent measurements of the longitudinal and shear wave velocities were performed using a laser ultrasonic technique while heating the sample in a vacuum-pumped heating stage. The elastic moduli increased linearly with increasing hydrogen content and decreased by ~10% during heating from room temperature to 1000°C in the three YH x compositions. The linear relationship between the elastic moduli and the hydrogen content in yttrium hydride was verified by atomistic calculations based on density functional theory (DFT). The absence of abrupt changes in the temperature-dependent measurements of elastic modulus of the YH x samples suggested negligible loss of hydrogen at elevated temperatures. Excellent agreement was found between the measured and calculated dependence of the elastic moduli on the stoichiometry, thereby providing a new approach for investigating the effects of fabrication-induced parameters (such as porosity) on the elastic moduli. Furthermore, this study demonstrates the utility of the combined approach involving DFT-based atomistic calculations and measurements of the elastic moduli for the informative development of metal hydrides and can be used as a metric for novel moderator materials investigations for emerging microreactors and beyond.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗