Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT calculations”

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

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

250 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗
DOE OSTI2026

Reducing experimental uncertainty in the calculation of cavity swelling in a transmission electron microscope through crystallographically aided void volume tracking

Accurate measurement of cavity swelling in a transmission electron microscope is essential to define material performance under irradiation, and the conventionally applied spherical assumption for the calculation of cavity volumes in irradiated materials can result in errors ranging between a 25% underestimation and 72% overestimation of volume purely based on the assumption of shape. This assumption is undeniably expedient for calculation but does not sufficiently account for the 3D nature of polyhedral cavities and their shape projection in the transmission electron microscope, and therefore presents too large of an associated uncertainty in swelling determination for faceted cavities. This uncertainty has been defined for many common cavity shapes in FCC and BCC crystal systems, and has been tabulated across the cubic fundamental region. A revised methodology for crystallographically aided void volume tracking, or CAVV-T, is presented and demonstrated on a specimen of neutron irradiated Ni. In-depth discussion on the application of this technique is provided along with resources to allow for conversion between the spherical assumption and this revised method. This work seeks to increase experimental confidence in the characterization and quantification of critical aspects of irradiation damage in materials by applying a crystallographically-resolved approach for cavity swelling calculation.

Haag IV, James V. [Pacific Northwest National Labo↗
NASA NTRS1977

Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting

Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.

O I Ornatskaya↗
DOE OSTI2024

Thermodynamically informed priors for uncertainty propagation in first-principles statistical mechanics

Here, this work demonstrates how first-principles statistical mechanics approaches within a Bayesian framework can quantify and propagate uncertainties to downstream thermodynamic calculations. To address the issue of Bayesian prior selection, knowledge of 0 K ground states in the material system of interest is incorporated into the prior. The effectiveness of this framework is shown by creating a phase diagram for the fcc zirconium nitride system, including confidence intervals on order-disorder transition temperatures.

Bayesian methods↗
DOE OSTI2025

Distinguishing Thermal Fluctuations from Polaron Formation in Halide Perovskites

Recent angle-resolved photoelectron spectroscopy (ARPES) measurements of the hole effective mass in CsPbBr$_3$ revealed an enhancement of $\sim$50 % compared to the bare mass computed from first principles for CsPbBr$_3$ at $T = 0 K$. This large enhancement was interpreted as evidence of polaron formation. Employing accurate finite-temperature first-principles calculations, we show that the calculated hole effective mass of CsPbBr$_3$ at $T = 300 K$ can explain experimental results without invoking polarons. Thermal fluctuations are particularly strong in halide perovskites compared to conventional semiconductors such as Si and GaAs, and cannot be ignored when comparing with experiment. We not only resolve the debate on polaron formation in halide perovskites, but also demonstrate the general importance of including thermal fluctuations in first-principles calculations for strongly anharmonic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗
NASA NTRSDate not supplied

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation↗
DOE OSTI2025

Quantum mechanical dataset of 836k neutral closed-shell molecules with up to 5 heavy atoms from C, N, O, F, Si, P, S, Cl, Br

Abstract We introduce the Vector-QM24 (VQM24) dataset comprehensively covering all possible neutral closed-shell small organic and inorganic molecules with up to five heavy (p-block) atoms: C, N, O, F, Si, P, S, Cl, Br. All valid stoichiometries, Lewis-rule-consistent graphs, and stable conformers (identified via GFN2-xTB) were enumerated combinatorially, yielding 577k conformational isomers spanning 258k constitutional isomers and 5,599 unique stoichiometries. DFT (ωB97X-D3/cc-pVDZ) optimizations were performed for all, and diffusion quantum Monte Carlo (DMC@PBE0(ccECP/cc-pVQZ)) energies are provided for 10,793 lowest-energy conformers with up to 4 heavy atoms. VQM24 includes structures, vibrational modes, rotational constants, thermodynamic properties (Gibbs free energies, enthalpies, ZPVEs, entropies, heat capacities), and electronic properties such as atomization, electron interaction, exchange-correlation, dispersion energies, multipole moments (dipole to hexadecapole), alchemical potentials, Mulliken charges, and wavefunctions. Machine learning models of atomization energies on this dataset reveal significantly higher complexity than QM9, with none achieving chemical accuracy. VQM24 offers a rigorous, high-fidelity benchmark for evaluating quantum machine learning models.

Science & Technology - Other Topics↗
DOE OSTI2025

Fast lithium ion diffusion in brownmillerite Li x Sr 2 Co 2 O 5

Transition metal oxides not only exhibits novel magnetic properties but also provides outstanding ionic transports. Ionic conductors have great potential for interesting tunable physical properties via ionic liquid gating and novel energy storage applications such as all-solid-state lithium batteries. In particular, low migration barriers and high hopping attempt frequency are the keys to achieve fast ion diffusion in solids. Taking advantage of the oxygen-vacancy channel in Li x Sr 2 Co 2 O 5 , we show that migration barriers of lithium ion are as small as 0.28–0.17 eV depending on the lithium concentration rates. Our first-principles calculation also investigated hopping attempt frequency and concluded the room temperature ionic diffusivity and ion conductivity are high as 10 −7 –10 −6 cm 2 s −1 and 10 −3 –10 −2 Scm −1 , respectively, which outperform most of perovskite-type, garnet-type, and sulfide Li-ion solid-state electrolytes. This work proves Li x Sr 2 Co 2 O 5 as a promising super-ionic conductor.

Crystallographic defects↗
NASA NTRS1995

A Case Study of View-Factor Rectification Procedures for Diffuse-Gray Radiation Enclosure Computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Robert P Taylor↗
DOE OSTI2024

Phase stability in the Hf-N and Zr-N systems

Hf and Zr nitrides are promising compounds for many technologically important areas, including high-temperature structural applications, quantum computing, and solar and optical applications. Here, this article reports on a comprehensive first-principles statistical mechanics study of phase stability in the Hf-N and Zr-N binary systems. A high solubility of nitrogen in the hcp forms of Hf and Zr is predicted. The rocksalt forms of HfN and ZrN can also tolerate a high degree of off-stoichiometry through the introduction of nitrogen and metal vacancies. The Hf-N binary favors a family of stacking faulted parent crystal structures at intermediate nitrogen concentrations that host a unique form of short-range order among nitrogen interstitials and vacancies. These phases can accommodate some degree of configurational entropy and remain ordered to temperatures as high as 1200 K.

Monte Carlo methods↗
NASA NTRS1975

Interpretation of Various Radiation Backgrounds Observed in the Gamma-Ray Spectrometer Experiments Carried on the Apollo Missions and Implications for Diffuse Gamma-Ray Measurements

Since the report of a preliminary analysis of cosmic gamma-ray measurements made during the Apollo 15 mission, an improved calculation of the spallation activation contribution has been made including the effects of short-lived spallation fragments, which can extend the correction to 15 MeV. In addition, a difference between Apollo 15 and 16 data enables an electron bremsstrahlung contribution to be calculated. A high level of activation observed in a crystal returned on Apollo 17 indicates a background contribution from secondary neutrons. These calculations and observations enable an improved extraction of spurious components and suggest important improvements for future detectors.

C S Dyer↗
DOE OSTI2025

Energy filtering–induced ultrahigh thermoelectric power factors in Ni 3 Ge

Traditional thermoelectric materials rely on low thermal conductivity to enhance their efficiency but suffer from inherently limited power factors. Innovative pathways to optimize electronic transport are thus crucial. Here, we achieve ultrahigh power factors in Ni 3 Ge-based systems through an unconventional thermoelectric materials design principle. When overlapping flat and dispersive bands are engineered to the Fermi level, charge carriers can undergo intense interband scattering, yielding an energy filtering effect similar to what has long been predicted in certain nanostructured materials. Via a multistep DFT-based screening method developed here, we find a family of L1 2 -ordered binary compounds with ultrahigh power factors up to 11 mW m −1 K −2 near room temperature, which are driven by an intrinsic phonon-mediated energy filtering mechanism. Our comprehensive experimental and theoretical study of these intriguing materials paves the way for understanding and designing high-performance scattering-tuned metallic thermoelectrics.

Science & Technology - Other Topics↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

First-Principles Investigation of Phase Transitions in RuNb, RuTa, and Ru-Nb-Ta Ultrahigh Temperature Shape Memory Alloys

Ultrahigh temperature shape memory alloys (UHT-SMAs) have transition temperatures above 600 C, and they have found applications for sensing and actuating devices in aerospace industry. Among very few such UHT-SMAs currently known are Ru-based alloys such as RuNb and RuTa, whose martensite structures and phase transitions are totally different from those of NiTi-based SMAs and were poorly understood. In this work, we carried out a systematical study of RuNb using first-principles total energy calculations and molecular dynamics (MD) simulations. The transition paths and mechanisms in cubic → tetragonal →monoclinic transitions are revealed. and the transition sequence and martensitic transition temperatures are determined (MTTs) by evaluating the Gibbs free energies using thermodynamic integration. The calculated MTTs are in very good agreement with the experimental data. We found that the monoclinic phase at the second transition has the 𝑃21/m symmetry instead of experimentally identified 𝑃2/m. Our calculations demonstrate that RuTa has very similar phase transitions to those of RuTa. Furthermore, we studied the Ru0.5Nb0.25Ta0.25 ternary. Our results can explain the measured significant drop of MTT (~ 200 C) in the second transition for Ru0.5Nb0.25Ta0.25 compared with those of RuNb and RuTa, while in the first transition its MTT is between MTTs of RuNb and RuTa. The insights gained by this study and the verified ab initio methods for accurate MTT calculations can be applied to fast screen and quantitatively design novel UHT-SMAs having similar properties with desirable MTTs and much reduced cost.

Zhigang Wu↗
DOE OSTI2026

Mu2e Production Target Thermal Test Project

The Mu2e experiment requires a production target that is capable of operating under extreme thermal conditions caused by an 8 GeV proton beam. This project’s objective supports the development of the Mu2e Pro- duction Target by testing the Stickman model’s thermal behavior. Angel Flores Luviano has assisted Jonathan Williams in progressing this project by contributing to the development of a Radiative Cooling Test Fixture (RCTF) that will be used to evaluate the thermal behavior of the new production target model. Engineering calculations were performed to an- alyze thermal performance and pressure drop within the cold well and wa- ter cooling circuit. These calculations also determined the optimal sizing for key components of the water system. An engineering note was made to document these calculations. CAD models of the water circuit piping, thermocouple mounting bars, and radiator chimney were developed, and prototype test rig components were fabricated using 3D printing. Future work will focus on continued RCTF development which includes control system integration, heater hardware design, and interfaces that can be scaled up to increase thermal capacity.

Flores Luviano, Angel Augusto [Unlisted, US]↗
NASA NTRS1977

The Chemical Composition of the Cores of the Terrestrial Planets and the Moon

Using models of the quasi-chemical theory of solutions, the activity coefficients of silicon are calculated in the melts Fe-Si, Ni-Si, and Fe-Ni-Si. The calculated free energies of solution of liquid nickel and silicon in liquid iron in the interval 0 to 1400 kbar and 1500 to 4000 K, shows that Fe-Ni-Si alloy is stable under the conditions of the outer core of the Earth and the cores of the terrestrial planets. The oxidation-reduction conditions are studied, and the fugacity of oxygen in the mantles of the planets and at the core-mantle boundary are calculated. The mechanism of reduction of silicon is analyzed over a broad interval of P and T. The interaction between the matter of the core and mantle is studied, resulting in the extraction of silicon from the mantle and its solution in the material of the core. It is concluded that silicon can enter into the composition of the outer core of the Earth and Venus, but probably does not enter into the composition of the cores of Mercury, Mars, and the Moon, if in fact the latter possesses one.

O L Kuskov↗