Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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

Hydrogen induced dislocation core reconstruction in bcc tungsten

Dislocation, playing a crucial role in the plastic deformation of metals, can be significantly affected by introducing solute elements. Hydrogen (H) embrittlement is one such example, while the underlying mechanism for H affected dislocation structural stability and mobility remains unclear and the role of H has been controversial. Here, using first-principles calculations, we demonstrate that the effect of H on screw dislocation in bcc tungsten (W) is H concentration-dependent, signified by a H-induced transition of SD core structure. At low concentrations of H segregation, dislocation maintains the intrinsic easy-core structure, and H atoms are attached to the “periphery” of dislocation to enhance dislocation motion. In contrast, at high H concentrations, dislocation transforms into a hard-core, metal hydride-like structure, as H atoms become the “body” of dislocation to significantly reduce the dislocation mobility in W. Further, such local easy-to-hard transition can also be induced by the other solutes, including helium, carbon, nitrogen and oxygen. Our work sheds new light on the H-dislocation interactions in bcc W, having broad implications in the interstitial solute-related phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The atomic-level structure and stability of interfaces of Pt nanoparticles in alumina: An experimental and computational evaluation

The atomic-level structure of interfaces between Pt and a transition form of Al2O3 were studied using a combination of electron microscopy and first principles calculations. A model system of Pt nanoprecipitates in Al2O3 were formed in sapphire wafers via high-energy ion implantation of Pt followed by thermal annealing at 1000 °C in air. The Pt nanoparticles took the form of tetrahedra and truncated tetrahedra primarily bound by {111}Pt facets. The high prevalence of these facets motivated the development of density functional theory (DFT) based models of (111)Pt interfaces with six different chemical terminations of ( 2 ¯ 01 ) θ-alumina. The atomic-level structure of the Pt/Al2O3 interfaces was characterized with aberration-corrected scanning transmission electron microscopy (STEM) and the experimental images were compared to STEM image simulations of the DFT models. The model interface with Pt bonded to oxygen-terminated θ-Al2O3, with the Pt located on top of the O and with an underlying layer of octahedral Al, provided the best match to the experimental images. This interfacial termination is also the most stable for the thermal annealing conditions used based on thermodynamic calculations of the interfacial energy as a function of temperature and oxygen partial pressure. This experimentally verified model provides a basis for improving models of Pt/γ-alumina interfaces.

Clauser, A L↗
DOE OSTI2025

Fundamental microscopic properties as predictors of large-scale quantities of interest: Validation through grain boundary energy trends

Correlations between fundamental microscopic properties computable from first principles, which we term canonical properties, and complex large-scale quantities of interest (QoIs) provide an avenue to predictive materials discovery. Here, we propose that such correlations can be efficiently discovered through simulations utilizing approximate interatomic potentials (IPs), which serve as an ensemble of “synthetic materials”. As a proof of principle we build a regression model relating canonical properties to the symmetric tilt grain boundary (GB) energy curves in face-centered cubic crystals, characterized by the scaling factor in the universal lattice matching model of Runnels et al. (2016), which we take to be our QoI. Our analysis recovers known correlations of GB energy to other properties and discovers new ones. We also demonstrate, using available density functional theory (DFT) GB energy data, that the regression model constructed from IP data is consistent with DFT results, confirming the assumption that the IPs and DFT belong to same statistical pool and thereby validating the approach. Regression models constructed in this fashion can be used to predict large-scale QoIs based on first-principles data and provide a general method for training IPs for QoIs beyond the scope of first-principles calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influence of isoelectronic mass modulation on phonon decay and transport in layered pnictogen-carbon systems

This study systematically investigates the structural, electronic, and thermal transport properties of layered pnictogen–carbon (Pn 2 C 2 ; Pn = P, As, Sb, Bi) monolayers, revealing a profound influence of isoelectronic mass modulation on phonon dynamics and thermal conductivity. Through first-principles calculations and lattice dynamics analysis, we demonstrate that increasing the atomic mass of pnictogens leads to elongated Pn-C bonds, weakened interatomic interactions, and enhanced phonon-phonon scattering, resulting in a dramatic reduction in lattice thermal conductivity (κ). Specifically, κ decreases exponentially from 65.6 W/m-K in P 2 C 2 to an ultralow 0.37 W/m-K in Bi 2 C 2 , driven by suppressed acoustic-optical phonon gaps Δ$^{Γ}_{A-0}$, increased anharmonicity, and reduced phonon lifetimes and group velocities. The transition from semiconducting to metallic behavior in heavier pnictogen-based systems further enhances phonon-electron scattering, contributing to thermal conductivity suppression. Structural analysis highlights the strengthening of out-of-plane C—C bonds and the contraction of C—Pn—C bond angles, which disrupt in-plane phonon transport. Further, these findings establish a design framework for engineering ultralow-κ materials through chemical substitution and structural tuning, offering significant potential for thermoelectric and thermal management applications. This work provides fundamental insights into phonon decay mechanisms and thermal transport in 2D materials, paving the way for advanced phononic and energy-efficient devices.

2D materials↗
DOE OSTI2025

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗
DOE OSTI2026

Chemical complexity suppresses interstitial solute segregation to screw dislocation cores in the bcc NbTaTiHfZr high-entropy alloy

Recent reports suggest that interstitial solute additions to body-centered cubic (bcc) high-entropy alloys (HEAs) may enhance their mechanical properties. However, details of interactions between interstitial atoms and dislocations in these HEAs remain incompletely understood. Using first-principles calculations, we examine the energetics of C, N, and O interstitial solutes in elemental bcc metals and NbTaTiHfZr, focusing on their segregation in screw dislocation cores. We examine two types of core sites and show that the low-energy configuration and associated segregation energy depend on the transition metal group. In NbTaTiHfZr, we find that chemical complexity substantially suppresses interstitial segregation in dislocation cores, as fluctuations in bulk solute energies induce energetically favorable sites that disfavor segregation to core sites. Local chemical order further enhances this effect by promoting solute clustering away from dislocation cores. These findings reveal fundamental differences between the behavior of interstitial atoms in elemental metals and HEAs, with implications for high-temperature plasticity and dynamic strain aging.

Body-centered cubic transition metals↗
DOE OSTI2025

Water effects on NH 3 -SCR over Cu-based small-pore zeolite catalysts: A review

The selective catalytic reduction (SCR) of NO x with NH 3 as the reductant over small-pore Cu-zeolite catalysts is a leading strategy for controlling emissions from diesel engines and mobile sources. While the hydrothermal stability of Cu-zeolites such as Cu-SSZ-13, Cu-SAPO-34, and Cu-SSZ-39 has been extensively studied, the multifaceted role of water vapor under operating conditions has only recently received systematic attention. Water impacts multiple aspects of SCR catalysis, including NO x conversion, N 2 O formation, Cu ion speciation, redox kinetics, and surface reactions. This review comprehensively summarizes recent advances in understanding the effect of water on NH 3 -SCR over Cu-based catalysts. For Cu-SSZ-13, water vapor induces both inhibitory and promotive effects depending on temperature, Cu loading, and Si/Al ratio by altering NH 3 adsorption, Cu mobility, and redox half-cycle kinetics. Studies have revealed that water can suppress ammonium nitrate formation at low temperatures and enhance nitrate pathways at higher temperatures. Over Cu-SAPO-34, water promotes NO conversion by improving Cu(II) reducibility and strengthening Brønsted acidity. First-principles calculations confirm that water-coordinated Cu species lower O 2 activation barriers. In Cu-SSZ-39, water was found to facilitate Cu reduction, migration, and O 2 activation, resulting in superior water tolerance and unexpected catalytic promotion. Both experimental and computational studies highlight that water vapor exerts complex, catalyst-dependent effects on NH 3 -SCR performance, underscoring the need for tailored catalyst designs that leverage beneficial water interactions while minimizing inhibitory impacts.

Ammonia selective catalytic reduction↗
DOE OSTI2020

The potential application of VS 2 as an electrode material for Mg ion battery: A DFT study

Herein, by means of first-principle calculations based on density functional theory (DFT), the electrochemical properties of monolayer VS 2 (M-VS 2 ), double-layer VS 2 (D-VS 2 ) and bulk VS 2 (B-VS 2 ) as electrode materials for Mg-ion batteries (MIBs) were comprehensively explored. The computation results reveal that Mg atom can strongly bind with the three different forms of VS 2 . All of the Mg adsorbed VS 2 systems demonstrate metallic characteristics, which indicates a good electronic conductivity. In addition, crystal orbital hamiltonian population shows that the stability of V-S bond is weakened after adsorption of Mg atom. The low diffusion barriers of Mg give rise to the high rate performance of VS 2 in MIBs. More interestingly, the three types of VS 2 display same storage ability for Mg cations, which can adsorb 0.5 Mg atoms for VS 2 , producing maximum theoretical capacity 233 mA h g -1 for MIBs. Furthermore, the average working voltages results suggest that M-VS 2 can be employed as anode materials, while D-VS 2 , B-VS 2 can be used as cathode materials for MIBs.

25 ENERGY STORAGE↗
DOE OSTI2025

First-Principles study of interface capture of Ni during high-temperature oxidation of NiCr alloys

The formation of precipitated phases of less reactive metals within an oxide scale during alloy oxidation can disrupt its structural integrity, significantly compromising the protective function of the oxide. To better understand this phenomenon, it is crucial to investigate the atomic-level diffusion mechanisms that drive precipitate formation. In this study, we employ first-principles calculations to examine the diffusion and aggregation of Ni atoms across the NiCr/Cr 2 O 3 interface during the high-temperature oxidation of NiCr alloys. Our results reveal that, in the absence of Cr vacancies, Ni atoms predominantly remain within the NiCr substrate, with minimal migration toward the Cr 2 O 3 layer. However, when Cr vacancies are present at both the NiCr/Cr 2 O 3 interface and within the Cr 2 O 3 bulk, Ni atoms exhibit a strong propensity to migrate from the NiCr substrate into the Cr 2 O 3 bulk, where they aggregate into Ni-rich regions with a corundum structure. This study provides valuable atomic-level insights into the diffusion and aggregation of less-reactive metals at the alloy/oxide interface, enhancing our understanding of solute capture and its contribution to the degradation of oxidation resistance in alloys.

Chen, Jianmin [Xiangtan University (China)] (ORCID↗
DOE OSTI2022

Modeling the thermodynamics of the FeTi hydrogenation under para-equilibrium: An ab-initio and experimental study

We report FeTi-based hydrides have recently re-attracted attention as stationary hydrogen storage materials due to favorable reversibility, good sorption kinetics and relatively low costs compared to alternative intermetallic hydrides. Employing the OpenCalphad software, the thermodynamics of the (FeTi) 1-x H x (0 ≤ x ≤ 1) system were assessed as a key basis for modeling hydrogenation of FeTi-based alloys. New thermodynamic data were acquired from our experimental pressure-composition-isotherm (PCI) curves, as well as first-principles calculations utilizing density functional theory (DFT). The thermodynamic phase models were carefully selected based on critical analysis of literature information and ab-initio investigations. Key thermodynamic properties such as dissociation pressure, formation enthalpies and phase diagrams were calculated in good agreement to our performed experiments and literature-reported data. This work provides an initial perspective, which can be extended to account for higher-order thermodynamic assessments and subsequently enables the design of novel FeTi-based hydrides. In addition, the assessed thermodynamic data can serve as key inputs for kinetic models and hydride microstructure simulations.

08 HYDROGEN↗
DOE OSTI2021

A Transient-Response methodology based on experiments and modeling for Cu-Redox Half-Cycle kinetic analysis on a Cu-SSZ-13 SCR catalyst

In this report the Standard SCR reaction catalyzed by Cu-SSZ-13 is a redox process consisting of a reduction half cycle (RHC) and an oxidation half cycle (OHC) that cycle the active Cu sites between the Cu(II) and Cu(I) states. In the current work, a transient-response methodology consisting of experimental transient response Cu-redox (TRCR) measurements and kinetic modeling was developed for detailed study of individual SCR Cu-redox half cycles. The TRCR protocol allows quantification of the reducible Cu density, Cu(II)-Cu(I) partitioning, relative native RHC and OHC rates, and limiting half cycle during SCR. The half-cycle kinetics are studied over a wide (170–450 °C) temperature, and in differential segments along the catalyst length using spatially resolved capillary inlet mass-spectrometry (SpaciMS). The protocol alone provides two independent measures indicating that OHC increases faster than RHC with temperature, and that at all but the lowest temperature SCR is RHC limited. Introducing the transient kinetic models allows half-cycle reaction pathways, orders and activation energies to be studied and determined; e.g., RHC involves surface and gas-phase NH 3 routes, and OHC with O 2 involves dimer formation. Application of the full methodology provides further confirmation and quantification of the protocol conclusions, and specifically the first measurements of RHC and OHC activation energies; specifically for parallel RHC surface and gas-phase routes, and experimentally obtained OHC activation energy vs. first-principle calculations in the literature. Beyond insights related to the degreened catalyst studied here, the methodology provides a broadly available approach for quantifying how aging processes selectively impact RHC and OHC.

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

A perspective on challenges and opportunities in developing high entropy-ultra high temperature ceramics

The emergence of a new class of high entropy ultra-high temperature ceramics (HE-UHTCs) is expected to possess exceptionally superior mechanical, oxidation, and erosion resistance as compared to conventional UHTCs. The improvement in the properties of HE-UHTCs is attributed to the high configurational entropy of the multi-component system stabilized into a single-phase, but with little discussion on entropy contributions. This perspective article presents a viewpoint on the critical scientific question regarding the actual role entropy plays in stabilizing and enhancing the thermo-mechanical properties of a multi-component system. The high throughput first-principle calculations, which can provide the entropy formation ability (EFA) of thermodynamically stable HE-UHTCs, can prove to be valuable toolset to this puzzle. HE-UHTCs show enhancement in the mechanical properties much higher than the predicted rule of the mixture, which cannot traditionally be explained by a mere solid solution formation effect. Additionally, isothermal oxidation studies of HE-UHTCs limited to furnace testing state remarkable oxidation resistance which is attributed to the passivity of complex, layered structure, and composition of the oxide scale. However, whether all the constituents in the HE-UHTC system react simultaneously or preferential oxidation takes place needs to be addressed. To address all these issues, an integrated computational and experimental approach has been postulated towards designing a thermal protection system (TPS) for re-entry applications. Succeeding the approach will provide new opportunities in the compositional space of HE-UHTCs with tailored properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Chemical bonding, phase stability and magnetic property in Sm 2 Fe 17 X 3 (X=H, C, N): A first-principles perspective

As a promising alternative to Nd–Fe–B magnets, the critical rare earth free Sm 2 Fe 17 X 3 (X = C, N) exhibits potential for high-performance magnets. However, their poor phase stability remains a major obstacle to developing bulk magnets. We investigated the phase stability and intrinsic magnetic properties of Sm 2 Fe 17 X 3 (X = H, C, N) using first-principles calculations and chemical bond analysis. The formation energies are negative, while the decomposition energies are −1.53, 0.348, and −0.74 eV per formula unit for X = H, C, and N, respectively, which is responsible for the weak thermal stability. Our chemical bond analysis reveals that the bonding asymmetry between Sm–X and Fe–X interactions creates local structural distortions and degrades the phase stability of Sm 2 Fe 17 X 3 . The project Crystal Orbital Hamilton Population (-pCOHP) analysis indicates that the Sm–X bonding remains positive up to the Fermi level, indicating stable bonding interactions. Here, in contrast, the Fe–X bonding becomes negative near the Fermi level, signifying anti-bonding contributions that reduce structural stability. Interstitial atoms X expand the lattice and enhance Fe magnetic moments, but Fe–X bonding suppresses neighboring Fe moments. Electron transfer from Sm to X modifies the valence state of Sm and the crystal field at the site, contributing to enhanced magnetocrystalline anisotropy in Sm 2 Fe 17 X 3 . Among the interstitial elements, carbon and nitrogen—with their larger atomic radius and higher electronegativity—induce greater lattice expansion and form stronger bonds with neighboring Sm and Fe atoms compared to hydrogen. Consequently, Sm 2 Fe 17 X 3 (X = C and N) exhibits better phase stability and significant improvement in magnetic properties.

Chemical bonding↗
DOE OSTI2020

Migration energy barriers and diffusion anisotropy of point defects on tungsten surfaces

We report tungsten is one of the most promising candidates for plasma-facing materials in future fusion devices, owing to its high performance under extreme irradiation conditions. However, irradiation-induced surface morphology varies significantly, depending on the irradiation type, fluence, and flux. Therefore, it is critical to examine the dynamics of point defects on tungsten (W) surfaces to understand how irradiation affects surface morphology. In this study, we employ the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method to search for potential migration paths of point defects on W (1 0 0), (1 1 0) and (1 1 1) surfaces. The first-principles calculations are then used to accurately determine the migration energy barriers. The obtained paths and barriers are incorporated into a kinetic Monte Carlo (KMC) model to determine trajectories and diffusivities, which are described by diffusion tensors to demonstrate their anisotropic features. Multiple diffusion mechanisms have been identified on different surfaces, with various anisotropy factors. Particularly, point defects on the W (1 1 0) surface have the highest diffusivities, with anisotropy factors independent of temperatures. In comparison, the anisotropy of the W (1 0 0) surface decreases as temperature increases, while the W (1 1 1) surface is isotropic for point defects. This study provides insights into defect transport properties on different surfaces, which are essential for understanding the early stages of irradiation-induced microstructural evolutions and surface morphology of tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2021

First-principles studies of the concentration-dependent tritium diffusion in the zirconium hydrides with and without Sn impurity

We performed first-principles calculations to study the 3 H diffusion in the zirconium hydrides with or without Sn impurity. Our results show that the formation of Zr 3 H x becomes preferable as the 3 H concentration increases. Based on the most stable configurations of Zr 3 H x (x = 0.5, 1.0, 1.5, 2.0), we studied the interstitial 3 H diffusion with or without Sn impurity. The results show that the interstitial 3 H diffusion becomes less probable in the hydrides with a higher 3 H concentration. When introducing a substitutional Sn impurity, the diffusion barriers increase. Therefore, the substitutional Sn in the diffusion pathway hampers the 3 H diffusion in Zr hydrides.

36 MATERIALS SCIENCE↗
DOE OSTI2023

First-principles study of interfaces in Al/SiC metal-matrix composite system

In this work, first-principles calculations were performed on the interfaces between Al and SiC, which is a widely used strengthening agent in aluminum metal-matrix-composites (Al/MMC). C-terminated interfaces have much larger work of adhesion than Si-terminated interfaces, indicating that the former has much stronger interfacial bonding. The electron localization function shows that the chemical bonding between Al and C has a strong covalent character, while the bonding between Al and Si is largely metallic. As a result of the vastly different chemical bonding, the work of adhesion for C-terminated interfaces increases with the number of dangling bonds at the interface, while the opposite trend was observed for Si-terminated interfaces. Additionally, the interface energy for Si-terminated interfaces is comparable to that for C-terminated interfaces, suggesting both types of terminations can coexist in the Al/SiC system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pressure modulated charge transfer and phonon interactions drive phase transitions in uranium–aluminum laves phases

Lanthanide AB 2 intermetallic compounds known as Laves phases have itinerant 3d and localized 4f electrons, which lead to interesting physical properties such as magnetic anisotropy and high Curie temperatures. Actinide Laves phases can display physical properties that are similarly intriguing. However, at reduced A–A spacing the C14 and C15 polytypes may exhibit larger wavefunction overlap for their 5f electron states and distinct characteristics for phases with more delocalized chemical bonding. The C36 polytype, on the other hand, is extraordinarily rare (<5% of known Laves phases). UAl 2 is the only known actinide Laves phase to show a pressure-controllable C15 → C36 transition. Here, we apply first principles calculations to determine the origin of the C15 → C36 phase transition and reveal the differences between the corresponding properties of each phase. Pressure increases lead to bond compression–induced electron transfer from Al to U, which drives dynamic instability in the C15 phonon modes because of the uniform U–U bonding environment. Further, the opposite phenomena is observed in C36: varied U–U bonding environments are vibronically more stable after charge transfer. We find that the interplay between charge transfer, chemical bonding, and phononic stability are central to predicting phase transitions and corresponding changes in physical properties for both C15 and C36 UAl 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of interstitial oxygen on ω transformations and twin formation in bcc NbTaTiHf multi-principal element alloy from first-principles

Transformation- and twinning-induced plasticity (TRIP and TWIP) have been reported to contribute to the low-temperature deformation of some body-centered cubic (bcc) multi-principal element alloys (MPEAs) containing large fractions of group IV transition metals. The influence of interstitial solutes on the mechanisms underlying these forms of plasticity, however, remains unclear. Using first-principles calculations, we study the effects of interstitial O atoms on the relative stability of bcc and ω phases and on unstable and twin boundary stacking fault energy profiles in a representative bcc MPEA with high group-IV elemental fraction: NbTaTiHf. We find that O additions generally promote the relaxation of ω configurations back to their parent bcc structure, therefore inhibiting ω transformation. Calculations of the Rice parameter for bulk bcc and phases, as well as bcc-ω interfaces, further show that ω formation is a potent embrittlement factor, an effect that is enhanced by O additions, suggesting that the formation of bcc-ω interfaces is energetically preferred over the formation of the bulk ω phase. By contrast, the Rice parameter for twin boundaries indicates that these interfaces do not embrittle the material, even with O atoms at twin boundaries, providing a more favorable pathway for plastic deformation compared to ω transformation.

Density functional theory↗