Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory development”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Development of efficient solar cells using combination of QSPR and DFT approaches

The awarded research project aimed to design and determine the effective organic dye-sensitizers from the ensemble of chromophores for the dye-sensitized solar cell (DSSC) and fullerene-derivatives (FDs) as acceptor for polymer solar cell (PSC) followed by explore the electron transfer mechanism using the combined quantitative structure-property relationship (QSPR) analysis in conjunction with density functional theory (DFT) and time-dependent DFT (TDDFT)-based calculations. Based on our proposal we have employed different in silico approaches to design photo-efficient organic dye-sensitizers for DSSCs and FDs as acceptor for PSCs with exhaustive screening methodology.

14 SOLAR ENERGY↗
DOE OSTI2020

Anion Association Strength as a Unifying Descriptor for the Reversibility of Divalent Metal Deposition in Nonaqueous Electrolytes

Developing next-generation battery chemistries that move beyond traditional Li-ion systems is critical to enabling transformative advances in electrified transportation and grid-level energy storage. Here, we provide the first evidence for common descriptors for improved reversibility of metal plating/stripping in nonaqueous electrolytes for multivalent ion batteries. Focusing first on the specific role of chloride (Cl – ) in promoting electrochemical reversibility in multivalent systems, rotating disk (RDE) and ring-disk electrode (RRDE) investigations were performed utilizing a variety of divalent cations (Mg 2+ , Zn 2+ , and Cu 2+ ) and the bis-(trifluoromethane sulfonyl) imide (TFSI – ) anion. By introducing varying concentrations of Cl – , a cooperative effect is observed between TFSI – and Cl – that yields the more reversible behavior of mixed electrolytes relative to electrolytes containing only TFSI – . This effect is shown to be general for Mg, Zn, and Cu electrodeposition, and mechanistic understanding of the role of Cl – in improving reversibility of TFSI-based electrolytes is obtained through the combination of R(R)DE experimental results and density functional theory (DFT) evaluation of the redox activity and thermodynamic stability of various TFSI- and Cl-based solution complexes of metal ions. The cooperative anion effect is further generalized to other mixed-anion systems, where systematic variations in anion association strength predicted from DFT (i.e., Cl – > OTf – ≈ TFSI – > BF 4 – > PF 6 – ) yield corresponding trends in redox potentials and improvements of ≥200 mV in the reversibility of metal deposition/dissolution. These results identify anion association strength as a common descriptor for the reversibility of divalent metal anodes and suggest a set of general design principles for developing new electrolytes with improved activity and stability.

25 ENERGY STORAGE↗
DOE OSTI2020

Tunable Magnetic Exchange between Rare-Earth Metal 5d and Iron 3d States: A Case Study of the Multiple Magnetic Transitions in Gd 6 FeBi 2 and the Solid Solutions Dy 6– x Gd x FeBi 2 (1 ≤ x ≤ 5) with Curie Temperatures in the Range 120–350 K

In this paper, the room-temperature magnet Gd 6 FeBi 2 was comprehensively characterized by means of temperature-dependent single-crystal X-ray diffraction, magnetization measurements, and experimental electron density and electronic structure computations. This work explores the electron-spin effects of Fe on this structure and suggests that Gd6FeBi2 shows structural features and exchange interactions, which are clearly distinguishable from other analogues, including the solid solutions Dy 6–x Gd x FeBi 2 (0 ≤ x ≤ 5). The unique traits of Gd 6 FeBi 2 and its derivatives encompass abnormal variations of lattice parameters with the temperature, as well as the presence of multiple magnetic transitions of complex nature. Based on the comprehensive analyses, it can be suggested that the unique magnetic response in this material is the result from tunable Fe spin states, which are coupled with strong Gd 5d–Fe 3d interactions, which differ from other members of this large, structural family. In the ground state, according to density functional theory calculations, the Fe atom carries two net electrons, whose spins are coupled in an antiparallel fashion to the spins of the Gd electrons. Near the Curie temperature, the Fe net moment is reduced intermittently, coupled with a multiple of magnetic transitions. The magnetic correlations are manifested in unexpected variations of the lattice parameters as a function of temperature, suggesting spin–lattice interactions. Furthermore, this study emphasizes the important role of Gd 5d electrons in tuning Fe 3d-based magnetic contribution, which enables better understanding on the spin coupling in other related compounds and sheds light on the development of new magnetic and spintronic compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Scalable neutral H 2 O 2 electrosynthesis by platinum diphosphide nanocrystals by regulating oxygen reduction reaction pathways

Despite progress in small scale electrocatalytic production of hydrogen peroxide (H 2 O 2 ) using a rotating ring-disk electrode, further work is needed to develop a non-toxic, selective, and stable O 2 -to-H 2 O 2 electrocatalyst for realizing continuous on-site production of neutral hydrogen peroxide. We report ultrasmall and monodisperse colloidal PtP 2 nanocrystals that achieve H 2 O 2 production at near zero-overpotential with near unity H 2 O 2 selectivity at 0.27 V vs. RHE. Density functional theory calculations indicate that P promotes hydrogenation of OOH* to H 2 O 2 by weakening the Pt-OOH* bond and suppressing the dissociative OOH* to O* pathway. Atomic layer deposition of Al 2 O 3 prevents NC aggregation and enables application in a polymer electrolyte membrane fuel cell (PEMFC) with a maximum r(H 2 O 2 ) of 2.26 mmol h -1 cm -2 and a current efficiency of 78.8% even at a high current density of 150 mA cm -2 . Catalyst stability enables an accumulated neutral H 2 O 2 concentration in 600 mL of 3.0 wt% (pH = 6.6).

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

Tuning electrochemically driven surface transformation in atomically flat LaNiO 3 thin films for enhanced water electrolysis

Structure–activity relationships built on descriptors of bulk and bulk-terminated surfaces are the basis for the rational design of electrocatalysts. However, electrochemically driven surface transformations complicate the identification of such descriptors. In this work, we demonstrate how the as-prepared surface composition of (001)-terminated LaNiO 3 epitaxial thin films dictates the surface transformation and the electrocatalytic activity for the oxygen evolution reaction. Specifically, the Ni termination (in the as-prepared state) is considerably more active than the La termination, with overpotential differences of up to 150 mV. A combined electrochemical, spectroscopic and density-functional theory investigation suggests that this activity trend originates from a thermodynamically stable, disordered NiO 2 surface layer that forms during the operation of Ni-terminated surfaces, which is kinetically inaccessible when starting with a La termination. Our work thus demonstrates the tunability of surface transformation pathways by modifying a single atomic layer at the surface and that active surface phases only develop for select as-synthesized surface terminations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Noble Gas Transport in the MSRE

This study explores the relationship between the physicochemical parameter known as Henry’s gas constant and gas transport across a two-layer film interface. The investigation utilized the Gibbs free energy, incorporating surface and volume terms to elucidate trends in enthalpy and entropy. Notably, our findings align with experimental data and offer predictive insights into the Henry’s gas constants for helium and krypton, which hold significance for future experiments and theoretical developments. Furthermore, this study enhances the Gibbs free energy theory pertaining to the liquid–gas interface. It underscores the substantial contribution of noble gases in this region to volumetric energy as temperature increases. Additionally, we employed Monte Carlo simulations to analyze the effective thermal neutron multiplication factor, denoted as k eff . Our analysis reveals a linear correlation between graphite density and uniform density as a function of temperature. For the 1D Molten Salt Reactor Experiment (MSRE) system, we employed the Mole code to conduct heat and mass transfer calculations. These computations enable us to ascertain the distribution of fuel temperature based on coefficients and thermal properties. We also studied delayed neutron precursors during fuel cycling, taking into account the drift of cycling fuel through Mole–Griffin coupling. Our model represents k eff and β eff across various volume flow rates and salinity compositions. Finally, this study leveraged xenon-135 for continuous on-line monitoring of fuel salts to investigate the impact of steady-state xenon-135 on the MSRE and to better understand its distribution. These efforts build upon previous research related to removal processes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Triggering C-N Coupling on Metal Oxide Nanocomposite for the Electrochemical Reduction of CO 2 and NO x ⁻ to Formamide

The co-electroreduction of CO 2 and NO x ⁻ (NO 3 ⁻ /NO 2 ⁻ ) to generate formamide (HCONH 2 ) offers an opportunity for downstream chemical and polymer manufacturing decarbonization; however, significant challenges lie in the C-N coupling and the associated low product selectivity. Herein, p-block metal oxides are incorporated in copper oxides to provide more accessible active sites for reactant adsorption and activation, tuning the reaction selectivity toward the formamide production. Through in situ Raman and synchrotron-based infrared spectroscopy measurements, C-N bond formation is demonstrated in real-time with the CuO x /BiO x catalyst, where the C-N bond is detected via a *CHO and *NH 2 intermediates formation, in agreement with the density functional theory calculations. When tested in a flow electrolyzer, a formamide yield rate of 134 ± 11 mmol h -1 g cat -1 is reported, the first report of co-electroreduction of CO 2 and NO x ⁻ to formamide beyond conventional H-cell measurements. These new insights on the C-N coupling mechanisms and scale-up capability provide directions for further development of electrocatalysts for the formamide production.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Common workflows for computing material properties using different quantum engines

The prediction of material properties based on density-functional theory has become routinely common, thanks, in part, to the steady increase in the number and robustness of available simulation packages. This plurality of codes and methods is both a boon and a burden. While providing great opportunities for cross-verification, these packages adopt different methods, algorithms, and paradigms, making it challenging to choose, master, and efficiently use them. We demonstrate how developing common interfaces for workflows that automatically compute material properties greatly simplifies interoperability and cross-verification. We introduce design rules for reusable, code-agnostic, workflow interfaces to compute well-defined material properties, which we implement for eleven quantum engines and use to compute various material properties. Each implementation encodes carefully selected simulation parameters and workflow logic, making the implementer’s expertise of the quantum engine directly available to non-experts. All workflows are made available as open-source and full reproducibility of the workflows is guaranteed through the use of the AiiDA infrastructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine learning potential assisted exploration of complex defect potential energy surfaces

Abstract Atomic-scale defects generated in materials under both equilibrium and irradiation conditions can significantly impact their physical and mechanical properties. Unraveling the energetically most favorable ground-state configurations of these defects is an important step towards the fundamental understanding of their influence on the performance of materials ranging from photovoltaics to advanced nuclear fuels. Here, using fluorite-structured thorium dioxide (ThO 2 ) as an exemplar, we demonstrate how density functional theory and machine learning interatomic potential can be synergistically combined into a powerful tool that enables exhaustive exploration of the large configuration spaces of small point defect clusters. Our study leads to several unexpected discoveries, including defect polymorphism and ground-state structures that defy our physical intuitions. Possible physical origins of these unexpected findings are elucidated using a local cluster expansion model developed in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unconventional surface state pairs in a high-symmetry lattice with anti-ferromagnetic band-folding

Many complex magnetic structures in a high-symmetry lattice can arise from a superposition of well-defined magnetic wave vectors. These “multi-q” structures have garnered much attention because of interesting real-space spin textures such as skyrmions. However, the role multi-q structures play in the topology of electronic bands in momentum space has remained rather elusive. Here we show that the type-I anti-ferromagnetic 1q, 2q and 3q structures in an face-centered cubic sublattice with band inversion, such as NdBi, can induce unconventional surface state pairs inside the band-folding hybridization bulk gap. Our density functional theory calculations match well with the recent experimental observation of unconventional surface states with hole Fermi arc-like features and electron pockets below the Neel temperature. We further show that these multi-q structures have Dirac and Weyl nodes. Our work reveals the special role that band-folding from anti-ferromagnetism and multi-q structures can play in developing new types of surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Controlling Structural, Electronic, and Energy Flow Dynamics of Catalytic Processes through Tailored Nanostructures

MoS 2 (molybdenum disulfide) is a highly-versatile catalyst material for support of numerous reactions from hydrodesulfurization and denitrogenation to the focus of this renewal proposal: hydrogenation of CO/CO 2 towards (higher) alcohols. At the same time, MoS 2 is a non-toxic, environmentally-benign and rather inert material- which under ambient conditions for some time has served as a lubricant and, more recently, as next-generation electronic material. The apparent contrast between inertness and stability in ambient, and catalytic activity under reactive conditions is puzzling and calls for a synergistic theoretical and experimental investigation with the long-term objective of enabling the rational design of MoS 2 -based catalysts for alcohol-formation reactions by providing a microscopic understanding of the environmental factors that determines site activity and selectivity. Our research project seeks answers to the questions (a) what conformation does MoS 2 adopt under reaction conditions (as opposed to that under ultrahigh vacuum and low temperatures)?; (b) what reaction pathways exist on such a material?; (c) how can the local environment of the active sites be manipulated so as to make MoS 2 an efficient catalyst for production of higher alcohol from syngas? In particular, research strategies will explore how the basal plane composed of sulfur atoms can be activated so as to exhibit a reactivity of its own, by addressing three research targets and building on extensive preliminary and enabling work: (1) vacancies and vacancy aggregates on the basal plane; (2) non-local catalyst transformation through alkali doping, hydrogenation and phase transition; (3) fabrication of a metal-nanoparticle-activated MoS 2 system, in which particle anchoring, reactive sites and pathways as well as selectivity are controlled by design, as an example of predictive development of a catalyst material,. All strategies are directed to improve the efficacy of the key reactive sites and selectivity of chemical pathways by design, to replace the inefficient methodology of trial and error in catalyst development. This research project represents a synergistic combination of computational guidance, foundational surface-science-based experiments and validation under reactive conditions that aims at transformative new insights into the working of MoS 2 -based hydrogenation-catalysts, a topic squarely at the center of the interest of DOE BES. Guided and led by Talat Rahman, a computational physicist, this project will apply density functional theory to understand structure, reaction pathways and chemical potentials associated with MoS 2 -based CO/CO 2 hydrogenation, augmented by kinetic Monte Carlo methods for reaction rates and prefactors as well as ab-initio molecular dynamics for evaluation of thermal stability. Complementary experimental input and validation will originate from co-PI Ludwig Bartels, a physical chemist and materials scientist, whose group focuses on local imaging and preparation of MoS 2 materials, from co-PI Peter Dowben, an experimental physicist, whose group is expert in the spectroscopy of occupied and unoccupied electronic states, and from senior collaborator Michael White of Brookhaven National Laboratory, whose group generates high-resolution electronic and activity information on size-selected well-defined metal chalcogenide clusters. This collaborative effort will enable a comprehensive understanding of the correlation of structural integrity and catalytic activity of MoS 2 in forms ranging from extended films to individual particles with known geometries and binding sites. Alcohol formation from syngas is a rapidly emerging application that has great potential through facile, economic, and decentralized biomass gasification. CO 2 activation is one of the most pressing concerns of our time: increasing CO 2 levels in the atmosphere change the climate and expose the globe to environmental transformations with the potential for enormous economic and societal impact. We will investigate CO/CO 2 hydrogenation via an interdisciplinary research collaboration with established synergy – one that, in accordance with the mission of the DOE, involves accredited Hispanic-Serving Institutions and that, through student exchange with international collaborators and National Labs, will directly benefit a broad spectrum of communities and generate human resources in sciences essential to their future.

2D materials↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a 10 multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory (DFT). The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress-strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston-Tonks-Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability 20 of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid-liquid but never fully liquid state, never fully lost strength, 25 stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗
DOE OSTI2021

Scaling of Transition State Vibrational Frequencies and Application of d -Band Theory to the Brønsted–Evans–Polanyi Relationship on Surfaces

Semiempirical energy relations provide a means of estimating thermodynamic properties. Specifically, linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) relationships correlate adsorption energies between adsorbates across surfaces and reaction energies with activation barriers, respectively. Although vibrational scaling relations (VSRs) exist between adsorbates at identical sites, scaling between vibrational frequencies of adsorbed local minima and transition states is lacking. Here, we present density functional theory calculations for AH X (A = C, N, O) diffusions on transition metal surfaces and reveal linear scaling between frequencies of local minima and the transition state between those minima. Using d-band theory and linear muffin tin orbital theory (LMTO), we derive the slopes of these transition state vibrational scaling relations (TSVSRs) and, in so doing, provide a rigorous theory extending the original BEP relations developed for solution chemistry to surface chemistry. Furthermore, with a single reference DFT calculation, we predict the slopes and quantify uncertainty in the predictions.

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

pyGWBSE: a high throughput workflow package for GW-BSE calculations

Abstract We develop an open-source python workflow package, py GWBSE to perform automated first-principles calculations within the GW-BSE (Bethe-Salpeter) framework. GW-BSE is a many body perturbation theory based approach to explore the quasiparticle (QP) and excitonic properties of materials. GW approximation accurately predicts bandgaps of materials by overcoming the bandgap underestimation issue of the more widely used density functional theory (DFT). BSE formalism produces absorption spectra directly comparable with experimental observations. py GWBSE package achieves complete automation of the entire multi-step GW-BSE computation, including the convergence tests of several parameters that are crucial for the accuracy of these calculations. py GWBSE is integrated with Wannier90 , to generate QP bandstructures, interpolated using the maximally-localized wannier functions. py GWBSE also enables the automated creation of databases of metadata and data, including QP and excitonic properties, which can be extremely useful for future material discovery studies in the field of ultra-wide bandgap semiconductors, electronics, photovoltaics, and photocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Experimental trends and theoretical descriptors for electrochemical reduction of carbon dioxide to formate over Sn-based bimetallic catalysts

The electrochemical carbon dioxide reduction reaction (CO 2 RR) using renewable energy sources is a promising solution for mitigating CO 2 emissions. In particular, CO 2 RR to formate represents a commercially profitable target. However, a comprehensive understanding of the catalytic mechanisms of Sn-based catalysts under reaction conditions, including the real-time structural evolution of catalysts and the role of all key reaction intermediates in influencing the CO 2 RR selectivity, is still lacking. The current study reports a framework to study the selectivity preference of Sn-based bimetallic catalysts using a combination of electrochemical measurements, in situ characterization, and density functional theory (DFT) calculations. The addition of a second metal (Co, Ni, Ag, Zn, Ga, Bi) was found to play a vital role in affecting the CO 2 RR performance. In situ X-ray absorption near edge structure (XANES) measurements revealed a dynamic evolution in the Sn valence state induced by different secondary metals. A multidimensional descriptor involving all the key reaction intermediates was developed to assess formate selectivity using a 2-dimensional volcano plot. Finally, this research offers an effective framework for understanding CO 2 RR catalytic selectivity by considering both the real-time structural evolution of catalysts and all the key intermediates involved.

58 GEOSCIENCES↗
DOE OSTI2022

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations

The thermal neutron scattering law (TSL) for crystalline beryllium hydride (BeH 2 ) is developed from first-principles ab initio lattice dynamics calculations and the impact of neutron thermalization in this material on critical mass is estimated. BeH 2 has a body-centered orthorhombic crystal structure with 12 molecules per unit cell and a theoretical density of 0.755 g/cm 3 . The vibrational (phonon) densities of states for H and Be bound in BeH 2 are determined using VASP density functional theory and PHONON lattice dynamics calculations. The TSLs for H bound in BeH 2 , H(BeH 2 ), and Be bound in BeH 2 , Be(BeH 2 ), are then evaluated in the incoherent approximation from the calculated H and Be partial phonon density of states using FLASSH. Finally, critical mass as a function of 235 U loading density for bare and reflected BeH 2 moderated spheres is predicted from MC21 Monte Carlo neutron transport calculations using ENDF/B-VIII.0 cross sections and the H(BeH 2 ) and Be(BeH 2 ) TSL evaluations. Comparisons are made to water (H 2 O), polyethylene (CH 2 ), and beryllium oxide (BeO) as moderators. These critical mass predictions are a refinement upon the prior work by Rao and Srinivasan that neglected thermal neutron scattering effects. The minimum critical mass of a BeH 2 moderated assembly is estimated to be 0.207 kg 235 U for a 0.20 m thick BeO reflected sphere and 0.178 kg 235 U for a 0.40 m thick BeO reflected sphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗