Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials database”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Shuttle environment database

The proceedings of a workshop on the development of a data base for the Space Shuttle environment are summarized. The data base is intended for use by experimenters in designing scientific payloads for Shuttle missions. The status of currently available data on the Shuttle environment is reviewed, with emphasis given to ten specific topics which are believed to be of concern to Shuttle experiment designers. These topics include: flight dynamics, vibration and acoustics, orbiter motion effects, EMI, and thermal and humidity effects. Attention is also given to: the Shuttle particulate environment; surface interactions inside and on the exterior of the spacecraft, and the natural Shuttle environment. Recent experimental results concerning Vehicle Glow on the spacecraft surface material are discussed in detail. A schematic diagram describing the functional organization of computer hardware for the Shuttle environment data base is provided.

Lauriente, M.↗
DOE OSTI2021

Predicting Energetics Materials’ Crystalline Density from Chemical Structure by Machine Learning

To expedite new molecular compound development, a long-sought goal within the chemistry community has been to predict molecules’ bulk properties of interest a priori to synthesis from a chemical structure alone. In this work, we demonstrate that machine learning methods can indeed be used to directly learn the relationship between chemical structures and bulk crystalline properties of molecules, even in the absence of any crystal structure information or quantum mechanical calculations. We focus specifically on a class of organic compounds categorized as energetic materials called high explosives (HE) and predicting their crystalline density. An ongoing challenge within the chemistry machine learning community is deciding how best to featurize molecules as inputs into machine learning models—whether expert handcrafted features or learned molecular representations via graph-based neural network models—yield better results and why. We evaluate both types of representations in combination with a number of machine learning models to predict the crystalline densities of HE-like molecules curated from the Cambridge Structural Database, and we report the performance and pros and cons of our methods. Our message passing neural network (MPNN) based models with learned molecular representations generally perform best, outperforming current state-of-the-art methods at predicting crystalline density and performing well even when testing on a data set not representative of the training data. However, these models are traditionally considered black boxes and less easily interpretable. Here, to address this common challenge, we also provide a comparison analysis between our MPNN-based model and models with fixed feature representations that provides insights as to what features are learned by the MPNN to accurately predict density.

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

User Manual for the Thermal Hydraulic Code SE2RCT

SE2-ANL is a modified version of SUPERENERGY-2 which interfaces the power distribution in the ANL maintained physics codes for the user. SE2RCT is a special version of SE2-ANL developed for specific use on models of the Experimental Breeder Reactor-II (EBR-II). SE2RCT fulfills a specific role in the calculated results as part of the Fuels Irradiation and Physics Database (FIPD) and is built to connect with the Physics and Analysis DataBase (PADB) for metallic fuel research and development in advanced sodium-cooled fast reactors (SFRs) created for EBR-II. SE2RCT is primarily a thermal analysis code that interfaces with the RCT based pin power details instead of the built-in pin power dehomogenization strategy in SE2-ANL. Beyond the thermal hydraulic calculation, the SE2RCT code provides detailed axial distributions of the RCT calculated results. SE2RCT also calculates the the axial distribution of pin wise DPA for the cladding materials which is part of the FIPD calculated results. This report serves as the user manual for SE2RCT and thus eliminates the QA gap of missing documentation on the SE2RCT software. All of the features and limitations of SE2RCT are documented in this report. The input and output specifications of SE2RCT are discussed in detail, and some unique features of SE2RCT, e.g. accurate pin power and temperature results for problems with non-uniform composition within an assembly, are also demonstrated and tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Use of Frit‐Disc Crucible Sets to Make Solution Growth More Quantitative and Versatile

The recent availability of step‐edge, frit‐disc crucible sets (generally sold as Canfield Crucible Sets or CCS) has led to multiple innovations associated with the group's use of solution growth. The use of CCS allows for the clean separation of liquid from solid phases during the growth process. This clean separation enables the reuse of the decanted liquid, either allowing for simple, economic, savings associated with recycling expensive precursor elements or allowing for the fractionation of a growth into multiple, small steps, revealing the progression of multiple solidifications. Clean separation of liquid from solid phases also allows for the determination of the liquidus line (or surface) and the creation, or correction, of composition–temperature phase diagrams. The reuse of clean decanted liquid has also allowed to prepare liquids ideally suited for the growth of large single crystals of specific phases by tuning the composition of the melt to the optimal composition for growth of the desired phase, often with reduced nucleation sites. Finally, it is discussed how solution growth and CCS use can be harnessed to provide a plethora of composition–temperature data points defining liquidus lines or surfaces with differing degrees of precision to either test or anchor artificial intelligence and/or machine‐learning‐based attempts to augment and extend the limited experimentally determined database.

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

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Using FIPD and OPTD to Benchmark Metallic Fuel Performance

This report serves as an introduction, tutorial, and benchmark specification for out-of-pile tests on metallic fuel. It introduces a new user to the EBR-II legacy fuel performance test program and the fast reactor fuel performance databases built to preserve the records. It then details the information stored in each database and how to find it. A benchmark specification is included for a small set of out-of-pile tests on U-10Zr fuel to function as a tutorial demonstrating how the legacy fuel performance data sets stored in the FIPD and OPTD databases can be used together to benchmark fuel performance models for steady-state and transient performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

High fidelity multiphysics tightly coupled model for a lead cooled fast reactor concept and application to statistical calculation of hot channel factors

A tightly coupled multiphysics code system is established using the MOOSE framework for hot channel factor (HCF) evaluation on a Lead Fast Reactor (LFR) concept. The coupled system is driven by the Griffin multiphysics coupling capability under which the MOOSE Heat Transfer module and NekRS computational fluid dynamics solver are coupled for conjugate heat transfer using the Cardinal application. The coupled capability is demonstrated on an LFR assembly model based on materials and geometry of a prototypical lead-cooled fast reactor design by Westinghouse Electric Company, LLC. Moreover, the work integrates the Multiphysics Object Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM) to perform calculations for statistical analysis of HCF. Furthermore, the coupling strategy and workflow demonstrated in this paper is not only useful for predicting accurate hot channel factors for different kinds of advanced reactors but also for other engineering applications such as control rod worth assessment, generation of high-fidelity database for Artificial intelligence (AI)/machine learning (ML) training, design optimization and multi-resolution modeling.

Cardinal↗
NASA NTRS2009

Smoke Point in Co-flow Experiment

The Smoke Point In Co-flow Experiment (SPICE) determines the point at which gas-jet flames (similar to a butane-lighter flame) begin to emit soot (dark carbonaceous particulate formed inside the flame) in microgravity. Studying a soot emitting flame is important in understanding the ability of fires to spread and in control of soot in practical combustion systems space. Previous experiments show that soot dominates the heat emitted from flames in normal gravity and microgravity fires. Control of this heat emission is critical for prevention of the spread of fires on Earth and in space for the design of efficient combustion systems (jet engines and power generation boilers). The onset of soot emission from small gas jet flames (similar to a butane-lighter flame) will be studied to provide a database that can be used to assess the interaction between fuel chemistry and flow conditions on soot formation. These results will be used to support combustion theories and to assess fire behavior in microgravity. The Smoke Point In Co-flow Experiment (SPICE) will lead to a o improved design of practical combustors through improved control of soot formation; o improved understanding of and ability to predict heat release, soot production and emission in microgravity fires; o improved flammability criteria for selection of materials for use in the next generation of spacecraft. The Smoke Point In Co-flow Experiment (SPICE) will continue the study of fundamental phenomena related to understanding the mechanisms controlling the stability and extinction of jet diffusion flames begun with the Laminar Soot Processes (LSP) on STS-94. SPICE will stabilize an enclosed laminar flame in a co-flowing oxidizer, measure the overall flame shape to validate the theoretical and numerical predictions, measure the flame stabilization heights, and measure the temperature field to verify flame structure predictions. SPICE will determine the laminar smoke point properties of non-buoyant jet diffusion flames (i.e., the properties of the largest laminar jet diffusion flames that do not emit soot) for several fuels under different nozzle diameter/co-flow velocity configurations. Luminous flame shape measurements would also be made to verify models of the flame shapes under co-flow conditions. The smoke point is a simple measurement that has been found useful to study the influence of flow and fuel properties on the sooting propensity of flames. This information would help support current understanding of soot processes in laminar flames and by analogy in turbulent flames of practical interest.

Urban, David L.↗
DOE OSTI2022

Calibrating the SPECTACULAR constitutive model using legacy Sandia data for two filled epoxy systems: 828/CTBN/DEA/GMB and 828/DEA/GMB

The SPECTACULAR model is a development extension of the Simplified Potential Energy Clock (SPEC) model. Both models are nonlinear viscoelastic constitutive models used to predict a wide range of time-dependent behaviors in epoxies and other glass-forming materials. This report documents the procedures used to generate SPECTACULAR calibrations for two particulate-filled epoxy systems, 828/CTBN/DEA/GMB and 828/DEA/GMB. No previous SPECTACULAR or SPEC calibration exists for 828/CTBN/DEA/GMB, while a legacy SPEC calibration exists for 828/DEA/GMB. To generate the SPECTACULAR calibrations, a step-by-step procedure was executed to determine parameters in groups with minimal coupling between parameter groups. This procedure has often been deployed to calibrate SPEC, therefore the resulting SPECTACULAR calibration is backwards compatible with SPEC (i.e. none of the extensions specific to SPECTACULAR are used). The calibration procedure used legacy Sandia experimental data stored on the Polymer Properties Database website. The experiments used for calibration included shear master curves, isofrequency temperature sweeps under oscillatory shear, the bulk modulus at room temperature, the thermal strain during a temperature sweep, and compression through yield at multiple temperatures below the glass transition temperature. Overall, the calibrated models fit the experimental data remarkably well. However, the glassy shear modulus varies depending on the experiment used to calibrate it. For instance, the shear master curve, isofrequency temperature sweep under oscillatory shear, and the Young's modulus in glassy compression yield values for the glassy shear modulus at the reference temperature that vary by as much as 15 %. Also, for 828/CTBN/DEA/GMB, the temperature dependence of the glassy shear modulus when fit to the Young's modulus at different temperatures is approximately four times larger than when it is determined from the isofrequency temperature sweep under oscillatory shear. For 828/DEA/GMB, the temperature dependence of the shear modulus determined from the isofrequency temperature sweep under oscillatory shear accurately predicts the Young's modulus at different temperatures. When choosing values for the shear modulus, fitting the glassy compression data was prioritized. The new and legacy calibrations for 828/DEA/GMB are similar and appear to have been calibrated from the same data. However, the new calibration improves the fit to the thermal strain data. In addition to the standard calibrations, development calibrations were produced that take advantage of development features of SPECTACULAR , including an updated equilibrium Helmholtz free energy that eliminates undesirable behavior found in previous work. In addition to the previously mentioned experimental data, the development calibrations require data for the heat capacity during a stress-free temperature sweep to calibrate thermal terms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Comprehensive Approach towards Multi-Objective EMI Filter Design Optimization in High-Frequency SiC-Based Motor Drives

With rapid penetration of high-frequency, highdensity power electronics into the industry, controlling the system generated EMI becomes one of the major design challenges. Designing practical optimized EMI filters require simultaneous consideration of multiple aspects - including electromagnetic couplings, magnetic materials, practical component parasitics and their impacts at high-frequencies, and also the overall manufacturing costs. The multi-dimensional nature of the problem often results in overcompensated designs that hurt the system efficiency and power density. This paper focuses on a comprehensive design platform (developed in MATLAB) towards designing high-density, optimized, highly efficient EMI filters by analyzing the generated noise spectrum and evaluating a vast array of potential solutions before outputting an optimized filter solution. The proposed tool includes database-based component selection, optimal topology selection, multistage filter design, overall optimization for volume, mass, cost and total loss. The proposed filter has been thoroughly evaluated in a PSIM simulation environment emulating the hardware in CRD300DA12EXM3, a 300kW three-phase inverter from Wolfspeed with an RL load representing a motor drive application. Furthermore, a 20kHz switching frequency is chosen and impacts of the motor high-frequency characteristics on the generated EMI noise spectrum of the overall system have been analyzed.

42 ENGINEERING↗
NASA NTRS2004

Trajectory Analysis

The Systems Analysis Branch performs trajectory and systems analysis for next-generation launch vehicle and space transportation technologies. Currently the branch is supporting the Project Prometheus with analysis of nuclear electric and nuclear thermal propulsion missions to a variety of destinations. Within Project Prometheus a proposed mission to Jupiter and three of the planet sized icy moons that orbit it is developing. The Jupiter Icy Moons Orbiter (JIMO) as the project is being called will enable detailed scientific investigation of Jupiter's moons Callisto, Ganymede, and Europa. These moons were choose to orbit for intensive study in particular because they are each believed to have water, energy, and organic material. The JIMO mission will utilize nuclear fission power and electric propulsion in order to allow the spacecraft to orbit the moons at close range for long durations of time. My assignment this summer was to assist in developing a trajectory analysis for the spacecraft system the Jupiter Icy Moons Orbiter by rewriting an inefficient Excel file into a more efficient FORTRAN program. This program has been created for use planning the trajectory of the Jupiter Icy Moons Orbiter Mission. The program uses a database of thousands of data points representing flight time, burn time, thrust, mass of the propellant used, final mass of the spacecraft, ratio of finals mass to initial mass, and change in velocity that a spacecraft experiences during each phase in the Jupiter Icy Moons Orbiter Mission. The trajectory program is capable of taking a specific user entered specific impulse (isp), final mass fraction, and thrust (P(sub J)/m0) and through the use of cubic splines to fit specific data curves, the program can locate the exact flight time linked to the user specified values of specific impulse and final over initial mass fraction. Currently, the database used by the program to calculate flight times for a given thrust is only for isps in the range of 2000 to 9000 seconds. Although the program is specific to the Jupiter Icy Moons Orbiter mission, it can be easily modified to fit any mission with a propulsion system expected isp range between 2000 and 9000 seconds. This program performs the same tasks as the old Excel file but in a more timely fashion. Finally, this FORTRAN program will replace the Excel file as the tool used to calculate trajectory properties in missions and will provide valuable data to mission analysts for years to come.

Burke, Laura M.↗
NASA NTRS2016

Mars 2020 Entry, Descent and Landing Instrumentation 2 (MEDLI2)

The Mars Entry Descent and Landing Instrumentation 2 (MEDLI2) sensor suite will measure aerodynamic, aerothermodynamic, and TPS performance during the atmospheric entry, descent, and landing phases of the Mars 2020 mission. The key objectives are to reduce design margin and prediction uncertainties for the aerothermal environments and aerodynamic database. For MEDLI2, the sensors are installed on both the heatshield and backshell, and include 7 pressure transducers, 17 thermal plugs, and 3 heat flux sensors (including a radiometer). These sensors will expand the set of measurements collected by the highly successful MEDLI suite, collecting supersonic pressure measurements on the forebody, a pressure measurement on the aftbody, direct heat flux measurements on the aftbody, a radiative heating measurement on the aftbody, and multiple near-surface thermal measurements on the thermal protection system (TPS) materials on both the forebody and aftbody. To meet the science objectives, supersonic pressure transducers and heat flux sensors are currently being developed and their qualification and calibration plans are presented. Finally, the reconstruction targets for data accuracy are presented, along with the planned methodologies for achieving the targets.

Aerodynamics↗
DOE OSTI2025

Recent Metallic Fuel Data Recovery in FIPD

The Metallic Fuels Irradiation and Physics Database (FIPD) [1] is an organized collection of metallic fuel test pin data (U-xPu-yZr, = 0 ~ 28; y = 2 ~ 10) and documentation available to industry. FIPD mainly contains three types of data: (1) Fuel pin fabrication data, including fuel slug diameter, fuel slug length, cladding diameter, smear density, etc. (2) Fuel pin operation conditions, including axial distributions for power, temperatures, fluences, burnup, and isotopic densities, etc. and (3) Fuel pin post-irradiation examination (PIE) data, including fission gas release and gas chemistry, profilometry, and neutron radiography, etc. The operating conditions for pins with PIE data available in FIPD span significant ranges across key parameters. The fuel peak burnup extends from less than 5% up to 20 at%. The cladding peak temperature varies from about 490°C to 660°C. Finally, the cladding peak DPA shows a wide range from less than 5 to 120. These broad ranges reflect the diverse testing conditions and operational parameters captured in the available PIE data. More detail about FIPD can be found in ref. [2]. The database development is an ongoing effort covering metallic fuel experiments from the Experimental Breeder Reactor II (EBR-II) and the Fast Flux Test Facility (FFTF). As reported in the ref. [3, 4], most of the PIE data generated during the IFR program [5] has been collected, reviewed, processed, and integrated into FIPD. The most recently added PIE data can be found in ref. [4], which shows the collection of over 95% of the PIE data by the time of this paper. The recent improvements to the database are summarized in this paper.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Imaginary-time correlation function thermometry: A new, high-accuracy and model-free temperature analysis technique for x-ray Thomson scattering data

The accurate interpretation of experiments with matter at extreme densities and pressures is a notoriously difficult challenge. In a recent work [Dornheim et al., Nat. Commun. 13, 7911 (2022)], we have introduced a formally exact methodology that allows extracting the temperature of arbitrary complex materials without any model assumptions or simulations. Here, we provide a more detailed introduction to this approach and analyze the impact of experimental noise on the extracted temperatures. In particular, we extensively apply our method both to synthetic scattering data and to previous experimental measurements over a broad range of temperatures and wave numbers. We expect that our approach will be of high interest to a gamut of applications, including inertial confinement fusion, laboratory astrophysics, and the compilation of highly accurate equation-of-state databases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Gender in Mineral Names

Minerals are the fundamental constituents of Earth, and mineral names appear in scientific literature for disciplines including geology, chemistry, materials science, biology, and medicine, among others. Choosing a name is the full responsibility of the authors of new mineral proposals submitted to the International Mineralogical Association (IMA). Scientific nomenclature and its traditions have evolved over time and, consequently, mineral names track changes in the landscape of mineralogy with respect to language, technology, and culture. To evaluate these changes, the namesake information for all 5896 minerals approved by the IMA or ‘grandfathered’ into use as of December 2022 was recorded and categorized within a workable database. The compiled information yields diverse insights into the intersection of science and culture and could also be used to project future trends. In this study, we used the name database to investigate gender diversity among mineral eponyms. More than half (c. 54%) of all mineral species are named after people, the identities of whom are largely a reflection of the people that have historically been involved, in one way or another, in the geosciences and in the mining industry. Of the 2738 people with minerals named for them, approximately 6.1% are (interpreted to be) women. Nearly all minerals named for women were named during the last sixty years, although the rate of growth in the year-on-year percentage of women among new mineral namesakes has slowed since about 1985. If current and historical trends hold, our model predicts that women will not comprise more than about 10.35% of newly established mineral namesakes in future years. The representation of women among mineral namesakes also differs starkly among countries. For example, Russians comprise 43.11% of women with minerals named for them, but account for only 15.12% of all eponyms. However, there are additional disparities beyond the proportions of namesakes. For scientists who were alive when a mineral was named for them, women were an average of 3.74 years older than men when evaluated over the same timespan (1954–2022). These results demonstrate that gender-based disparities are imprinted into current mineral nomenclature and indicate that gender parity among new mineral namesakes is impossible without unprecedented changes in the upstream demographics that are most likely to affect naming trends.

58 GEOSCIENCES↗
DOE OSTI2023

An Interface–Type Memristive Device for Artificial Synapse and Neuromorphic Computing

Interface-type (IT) metal/oxide Schottky memristive devices have attracted considerable attention over filament-type (FT) devices for neuromorphic computing because of their uniform, filament-free, and analog resistive switching (RS) characteristics. The most recent IT devices are based on oxygen ions and vacancies movement to alter interfacial Schottky barrier parameters and thereby control RS properties. However, the reliability and stability of these devices have been significantly affected by the undesired diffusion of ionic species. Herein, a reliable interface-dominated memristive device is demonstrated using a simple Au/Nb-doped SrTiO 3 (Nb:STO) Schottky structure. The Au/Nb:STO Schottky barrier modulation by charge trapping and detrapping is responsible for the analog resistive switching characteristics. Because of its interface-controlled RS, the proposed device shows low device-to-device, cell-to-cell, and cycle-to-cycle variability while maintaining high repeatability and stability during endurance and retention tests. Furthermore, the Au/Nb:STO IT memristive device exhibits versatile synaptic functions with an excellent uniformity, programmability, and reliability. A simulated artificial neural network with Au/Nb:STO synapses achieves a high recognition accuracy of 94.72% for large digit recognition from MNIST database. These results suggest that IT resistive switching can be potentially used for artificial synapses to build next-generation neuromorphic computing.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Feature engineering descriptors, transforms, and machine learning for grain boundaries and variable-sized atom clusters

Abstract Obtaining microscopic structure-property relationships for grain boundaries is challenging due to their complex atomic structures. Recent efforts use machine learning to derive these relationships, but the way the atomic grain boundary structure is represented can have a significant impact on the predictions. Key steps for property prediction common to grain boundaries and other variable-sized atom clustered structures include: (1) describing the atomic structure as a feature matrix, (2) transforming the variable-sized feature matrix to a fixed length common to all structures, and (3) applying a machine learning algorithm to predict properties from the transformed matrices. We examine how these steps and different combinations of engineered features impact the accuracy of grain boundary energy predictions using a database of over 7000 grain boundaries. Additionally, we assess how different engineered features support interpretability, offering insights into the physics of the structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗