Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT”

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

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

At least 1,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evolution of dislocation loops in irradiated α-Uranium: An atomistically-informed cluster dynamics investigation

An atomistically-informed mean field cluster dynamics model has been developed to investigate the nucleation and growth of prismatic loops in irradiated a-Uranium. TEM analysis of neutron irradiated a-Uranium shows the evolution of self-interstitial atom and vacancy loops on (010) and (100) crystallographic planes, respectively, resulting in an anisotropic lattice swelling of its face-centered orthorhombic crystal. To provide model parameters, the crystallography of loops and the binding energy of point defects to these loops were studied using an angular dependent EAM potential and classical molecular dynamics (MD) simulations. Furthermore, using the bond-boost hyperdynamics method, the anisotropic diffusion of interstitials and vacancies in a-Uranium was investigated. Here, the mechanisms of point defect diffusion and the associated migration energies were reported and compared with previous DFT studies. The energetics and kinetic quantities mentioned above were used in the cluster dynamics model to predict the cluster density as a function of dose rate, dose and temperature and the results were compared to the reported neutron irradiation experiments. The model predictions reveal an accumulation of small sized vacancy loops along with a population of large and growing self-interstitial loops, which closely corresponds to the TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Distorted dislocation cores and asymmetric glide resistances in titanium

The determination of Critical Resolved Shear Stress (CRSS) in titanium for basal, prismatic, and pyramidal slip-planes without empirical constants is presented by combining Density Functional Theory (DFT) and anisotropic elasticity. A new mechanism leading to tension-compression (T-C) asymmetry in the CRSS levels has been revealed for the first time. The conditions for this asymmetry are established, involving a complex interplay between the dislocation core-structure and core-advance behavior. The three conditions for T-C asymmetry that need to be simultaneously satisfied can be summarized as: (1) a medium stacking fault width, d, (3 < d/b F < 10, where b F is the magnitude of the Burgers vector of the full dislocation), (2) an asymmetric core-structure of the extended dislocation (ξ 1 ≠ ξ 2 , where ξ 1 and ξ 2 are the core-widths of the first and second partials, respectively), and (3) intermittent motion of the partials (Δd/b F ≠ 0, where Δd is the magnitude of fluctuation in stacking fault width during intermittent motion). Pyramidal-slip in titanium satisfies all three conditions, resulting in significant T-C asymmetry. The CRSS theory considers a Wigner-Seitz (W-S) cell based domain area assigned to each lattice site for the calculations of core-energies accurately capturing the slip-plane lattice. The W-S based approach is essential due to the lower symmetry of the HCP crystal circumventing potential errors associated with the one-dimensional atomic-row approximation. Dislocation core structures are obtained for all the slip-systems in titanium showing significant distortion of the disregistry profile governing the core-advance behavior. The CRSS values predicted from the theory show agreement with the experimental CRSS levels reported in the literature.

Dislocation core↗
DOE OSTI2025

An energetic link between order and strength in metals: A nanocrystalline strength limit in high-entropy alloys and intermetallic compounds

The metallurgy and materials communities have long understood and exploited fundamental links between chemical and structural ordering in metallic solids to tailor their mechanical properties. We extend these ideas to include prediction of the nanocrystalline strength limit in high-entropy alloys and intermetallic compounds, where a breakdown occurs in the classical Hall-Petch strengthening behavior. The highest reported strength achievable through alloying has rapidly climbed and given rise to new classifications of materials with extraordinary properties, with a notable case being nanocrystalline metals. High-entropy alloys (chemically disordered, concentrated solid solutions) and intermetallic compounds are two boundary cases of how tailored order can be used to manipulate mechanical behavior. Here, we show that the complex electronic-structure mechanisms governing the peak strength of alloys and pure metals can be reduced to a few physically meaningful parameters based on their atomic arrangements and used – with no fitting parameters – to predict the maximum strength of these materials. This includes a generalized energy-based accounting for the degree of structural and chemical ordering that allows for rapid and reasonably accurate prediction of peak strength (validated in the nanocrystalline limit) as a function of temperature. Predictions of maximum strength based on the activation energy (with all materials properties derived from DFT calculations or experiments) for a stress-driven transition to an amorphous state is shown to accurately describe the breakdown in Hall-Petch behavior at the smallest crystallite sizes for pure metals, intermetallic compounds, high-entropy alloys, and metallic glasses. Further, this activation energy is also shown to be directly proportional to interstitial electronic charge density, which is a good predictor of ductility, stiffness (moduli), and phase stability in high-entropy alloys and solid metals generally. The proposed framework suggests the possibility of coupling ordering and intrinsic strength to mechanisms like dislocation nucleation, hydrogen embrittlement, and transport properties, such as through correlations between the activation energies for amorphization with stacking-fault and grain boundary energies. It additionally opens the prospect for greatly accelerated structural materials design and development to address materials challenges limiting more sustainable and efficient use of energy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Tailoring microstructures with mild magnetic-field processing: A case study of CuNiFe alloys

Combined experimental and computational investigations of the CuNiFe spinodal system confirm that application of a mild magnetic field during thermal treatment alters elemental redistribution and the resulting microstructure, relative to that obtained from zero-field annealing. Spinodal decomposition of a Cu 40 Ni 42 Fe 18 alloy was initiated during thermal treatment at 773 K, conducted either under zero field or modest (60 mT) magnetic f ield conditions for up to 200 h. Periodic (~10 nm) chemical modulations into Cu-rich and NiFe-rich regions were observed under both conditions, with the amplitude and wavelength of the segregated regions increasing with treatment time. However, magnetic field annealing resulted in a more than twofold increase in the amplitude of elemental modulations relative to zero-field conditions – consistent with enhanced diffusional f luxes during spinodal decomposition – while the modulation wavelength remained largely unaffected. These microstructural differences are reflected in various extrinsic magnetic properties. In parallel, first-principles DFT calculations indicate that long-range ferromagnetic order, as induced by an applied magnetic field, substantially alters the strength and nature of atomic interactions, enhancing the thermodynamic instability of the CuNiFe solid solution. Collectively, these results suggest that incorporating a mild (millitesla-level) magnetic field – distinct from the strong (tesla-level) fields commonly used in prior studies – during thermal processing has the potential to deliver enhanced control of microstructures for targeted engineering outcomes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase diagram of magnetic shape memory alloy Ni 50 Mn $50–x$ In $x$ , 0 < $x$ , 25 from first principles, via spin cluster expansion and phonon vibrational entropies

The metamagnetic shape memory Heusler alloy Ni 50 Mn $50–x$ In $x$ exhibits a rich phase diagram featuring competing magnetic states, coupled magnetic–structural phase transitions, and strong compositional sensitivity. Existing first-principles approaches struggletocapturetheintertwinedchemical, magnetic, andvibrationaleffectsinthesealloys, necessitating a more integrated modeling framework. We develop a spin cluster expansion (spin-CE) framework augmented by a quasi-harmonic phonon model to capture both configurational (chemical and magnetic) and vibrational contributions to the free energy of Ni 50 Mn $50–x$ In $x$ over the full compositional range 0 ≤ x ≤25. The spin-CE includes both chemical clusters and composition-dependent Ising spin interactions, with parameters fit to a first-principles density functional theory (DFT) dataset. Using this approach, we predict the complete magnetostructural phase diagram and transformation temperatures of Ni 50 Mn $50–x$ In $x$ across the composition space. We find that vibrational entropy alone is insufficient to reproduce the martensitic transformation in the magnetic shape memory alloy regime, highlighting the essential role of magnetism. Incorporating both magnetic and vibrational contributions allows us to reproduce all experimentally known phases, including the disappearance of the stable martensite phase at a critical In concentration and the Curie temperature of the austenite phase. The method also captures the transition with increasing In in martensite from antiferromagnetic to ferromagnetic order and predicts re-entrant ferromagnetism, though the latter occurs at higher In content than reported experimentally. We discuss possible sources of this discrepancy and highlight the broader applicability of the method to other magnetostructurally complex materials, where it may offer mechanistic insight and predictive design capabilities.

Cluster expansion↗
DOE OSTI2023

Influence of residual chlorine on Ru/TiO 2 active sites during CO 2 methanation

Titania-supported ruthenium (Ru/TiO 2 ) is an established catalyst for the hydrogenation of carbon dioxide to methane (Sabatier reaction). Chlorine contamination, owed to the RuCl 3 precursor, is demonstrated to have a detrimental impact on methanation activity. After calcination and reduction the catalyst contains residual chlorine, shown by XPS. An aqueous ammonia wash removes Cl without leaching Ru. The washed catalysts exhibit improvements in CH 4 site-time yields. Low Ru loading catalysts encounter the greatest activity enhancements after washing (~4.5-fold). DFT calculations indicate that chlorine and CO 2 directly compete for adsorption on Ru step sites, with Cl impeding the adsorption of CO 2 at under-coordinated sites and at higher Cl coverages. H 2 -chemisorption/TPR show that Cl removal lowers the onset of low temperature H 2 dissociation on Ru. Finally, DRIFTS provide evidence that the removal of Cl facilitates low temperature dissociative binding of CO 2 , indicated by the formation of surface bound linear CO species.

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

Dual-atom Ag 2 /graphene catalyst for efficient electroreduction of CO 2 to CO

Electrochemical reduction of CO 2 into value-added carbon compounds offers a promising strategy to mitigate global warming, but present challenges for chemistry due to the poor selectivity and stability of electrocatalysts. In this work, we report a dual-atom Ag 2 /graphene catalyst featuring well-defined AgN 3 -AgN 3 active site for CO 2 electrochemical reduction. This dual-atom catalyst can drive CO 2 reduction reaction at a potential as high as -0.25 V, and exhibit excellent CO Faradic efficiency up to 93.4 % with a current density of 11.87 mA cm -2 at -0.7 V and long-term stability, far surpassing the single-atom Ag 1 /graphene and the traditional silver nanoparticle catalysts. Finally, DFT calculations reveal that the dual-atom Ag site lowers the barrier for the formation of *COOH by stabilizing the *CO 2 through the concomitant interactions with the C and an O atom of CO 2 , resulting in excellent catalytic performance.

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

Passive NO x adsorption on Pd/H-ZSM-5: Experiments and modeling

A combined experimental and modeling study of the passive NO x adsorber (PNA) is presented that advances the understanding and prediction of the effects of various operating parameters and material properties of H-ZSM-5 and Pd/H-ZSM-5. Experiments reveal the NO x uptake process to be kinetically limited. Two microkinetic schemes are developed to explain the NO x uptake and release and NO 2 generation data. Both schemes are incorporated into a 1+1 D monolith model to predict and validate NO x uptake and temperature programmed desorption (TPD) data. Scheme I involves Z – [PdOH] + , Z – Pd 2+ Z – and Z – Pd + with Z – Pd + as the strongest NO binding site. Scheme II involves reduction of PdO 2 to PdO and nitrate formation. Model tuning utilizes a combination of infrared spectroscopy measurements and density functional theory (DFT) estimates of energy barriers. Lastly, the model is validated at different conditions including feed parameters and Pd loading.

42 ENGINEERING↗
DOE OSTI2020

HCOOH Decomposition on Sub-Nanometer Pd 6 Cluster Catalysts: The Effect of Defective Boron Nitride Supports Through First Principles

The catalytic properties of a hexagonal boron nitride- (h-BN) supported Pd 6 sub-nanometer cluster in the context of formic acid (HCOOH) decomposition were studied by means of periodic Density Functional Theory (DFT) calculations. The effect of support defectivity – boron (h-B v N) and nitrogen (h-BN v ) monovacancies – on the competition between the formate (HCOO)- and carboxyl (COOH)-mediated decomposition pathways was analyzed. Defects are responsible for charge-transfer leading to a positively or negatively charged cluster, and open new reactive channels in which vacancy-mediated dehydrogenation pathways can occur. Pd 6 cluster reconstructions, induced by the adsorption of reaction intermediates and by the presence of monovacancies in the support, greatly stabilize the formation of CO from COOH, which could drastically decrease the selectivity towards hydrogen production. Here, a simplified descriptor-based analysis, based on selected thermochemical quantities calculated on charged cluster models, suggests that Pd 6 sub-nanometer clusters supported on pristine h-BN and h-BN v can be more selective than Pd 6 supported on defective h-B v N towards HCOOH dehydrogenation.

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

NO x adsorption with CO and C 2 H 4 on Pd/SSZ-13: Experiments and modeling

A transient monolith model containing microkinetic schemes for NO uptake and release over Pd/SSZ-13 without and with reductants CO or C2H4 is presented. The scheme involves three cationic Pd species (Z – [PdOH] + , Z – Pd 2+ Z – , Z – Pd + ) as the active sites during uptake of NO, CO and C 2 H 4 , and their desorption and conversion at higher temperature. Kinetic parameters are estimated through a combination of density functional theory (DFT) estimates and a fit of uptake, desorption and conversion data. The tuned model is validated at different uptake temperatures, ramp rates, and flowrates. A “degree of uptake control” parameter is defined that helps to identify the step(s) that are uptake controlling. The model helps to interpret the data features and is used to identify operating conditions to meet application-relevant performance metrics, including NO trapping efficiency and NO release temperature. Furthermore, the model demonstrates higher NO uptake on Pd/SSZ-13 compared to Pd/ZSM-5.

42 ENGINEERING↗
DOE OSTI2022

Efficient Pd on carbon catalyst for ammonium formate dehydrogenation: Effect of surface oxygen functional groups

Formate solution is one of a kind promising liquid organic hydrogen carrier but suffers from kinetic challenges. This study investigates the kinetics of the surfacefunctionalized Pd on carbon catalysts for formate dehydrogenation and the impact of O-functional groups. Here, the fraction of the distinguished O-functional groups was modulated by the different concentrated HNO 3 solutions treatment or by H 2 reduction. This study shows that the O-functional groups play critical roles in dispersing Pd nanoparticles and decreasing the activation energy for dehydrogenation. Density functional theory (DFT) calculation reveals that most O-functional groups enhance formate adsorption on the Pd active site. However, the existence of C=O groups consumes reducing agents and hinders the formation of metallic Pd. The electron transfer from Pd to oxygen functional groups is unfavorable to dehydrogenation. The as-prepared Pd5/re-ACA (reduced activated carbon washed by acid) exhibited significant activity with a higher turnover frequency of 13,511 h -1 than commercial Pd/C.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

In-situ exsolution of bimetallic CoFe nanoparticles on (La,Sr)FeO 3 perovskite: Its effect on electrocatalytic oxidative coupling of methane

This study presents in-situ reduction of lanthanum strontium cobalt ferrite (LSCF) perovskite as an effective method for modifying its surface properties and enhancing its electrocatalytic reactivity for oxidative coupling of methane (OCM). The evolution of hetero-phases during the reduction of LSCF resulted in CoFe nanoparticles being formed at the surface. The in-situ reduced LSCF cell for OCM could be operated in either an ion pump or a fuel cell mode. High selectivity of 63% and 10.2% were reported for C 2+ hydrocarbons and C 3 H 6 , respectively. DFT calculations on LSCF and CoFe revealed that the high selectivity of C 2+ hydrocarbons on the LSCF primarily stems from the presence of CoFe nanoparticles. In-situ DRIFTS conducted under CH 4 proved that complete oxidation of CH 4 can be effectively inhibited by reducing LSCF, and control of oxygen supply is an important parameter for selective conversion of CH 4 to higher order hydrocarbon.

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

Using ethanol and isopropanol as biomass model compounds for understanding bond scission mechanisms over Cu/Mo 2 N catalysts

The conversion of biomass compounds into fuels and chemicals is an important step towards a more sustainable future. This work combines results from model surfaces and powder catalysts to demonstrate Cu-modified mo- lybdenum nitride (Cu/Mo 2 N) as a selective catalyst for dehydrogenation of the biomass model compounds, ethanol and isopropanol. Results from model surfaces showed that while Mo2N led to unselective decomposition via both dehydrogenation and dehydration, the addition of Cu increased the dehydrogenation activity and selectivity. DFT calculations showed how Cu influenced the structures of active sites, adsorbate interactions, and thus the product selectivity. Batch reactor studies on corresponding powder catalysts confirmed the trend that Cu modification increased dehydrogenation activity, and in situ X-ray absorption spectroscopy elucidated the Cu oxidation state under reaction conditions. Further, this work demonstrates a strategy for promoting dehydrogenation over Mo 2 N-based catalysts, as well as the feasibility of using model surfaces to guide the design of industrially relevant catalysts.

09 BIOMASS FUELS↗
DOE OSTI2024

Catalytic production of δ-valerolactone (DVL) from biobased 2-hydroxytetrahydropyran (HTHP) – Combined experimental and modeling study

δ-Valerolactone (DVL) is a five-carbon (C5) cyclic ester that can undergo ring-opening polymerization to yield high-performance, biocompatible polyesters. But current market prices of C5 chemicals like DVL are very high due to poor availability of C5 feedstock in petroleum. Herein, we demonstrate a novel route to DVL synthesis via dehydrogenation of biomass-derived 2-hydroxytetrahydropyran (HTHP) over Cu/SiO 2 without the use of toxic reagents. Since HTHP exists in thermal equilibrium with 3,4-dihydropyran (DHP) via dehydration, and with 2,2’-oxybis(tetrahydropyran) and 5-(tetrahydropyran-2-yloxy)pentanal via acetalization, we have also determined the thermochemistry (ΔH rxn and ΔG rxn ) of each competing reaction using density functional theory (DFT) calculations at the M06-2X/cc-pVTZ level. Finally, by developing a kinetic model of all 8 reactions involved, we have achieved 84% selectivity to DVL at 150°C in a packed bed reactor for over 72 hours of time on stream.

42 ENGINEERING↗
DOE OSTI2025

Electrosynthesis of high purity ethylene using high-index facet Cu 2 O nanocrystals electrocatalyst

Electrochemical CO 2 reduction reaction (eCO2RR) to multi-carbon (C 2+ ) products with copper-based catalysts is often limited by poor selectivity. This challenge arises from the concurrent formation of various intermediates, dictated by the atomic arrangement and electronic properties of surface atoms. In this study, we found that copper (I) oxide (Cu 2 O) nanocrystals with 50 facets (50F-NC), predominantly featuring (211) facets that offers high density of under-coordinated sites, demonstrate superior ethylene (C 2 H 4 ) selectivity of 92% ± 2 with an overall current density of 212 mA/cm 2 at -650 mV vs RHE. Furthermore, after one month of storage in a 1 M KOH electrolyte, this catalyst demonstrated a C 2 H 4 Faradaic efficiency of 87% highlighting its stabile structure under strong alkaline environments. Here, operando electrochemical Raman spectroscopy revealed enhanced CO* intermediate coverage on the 50F-NC catalyst, correlating with improved C-C coupling. SEM, TEM, and XPS analyses, along with DFT calculations, suggested that Cu sites on the (211) facet of 50F-NC and those at the Cu/Cu 2 O interface formed in-situ due to the surface reconstruction during the reaction, are likely active sites for effective C-C coupling and sustained high-rate C 2 H 4 production.

Cu nano particle↗
DOE OSTI2025

Mechanistic insights into N 2 O formation as a side product in NH 3 -SCR over small pore Cu-zeolites

Here, the present contribution provides clarity to N 2 O formation mechanisms and key influencing factors during low temperature NH 3 -SCR, with the goal of enabling the rational design of advanced SCR catalysts with low greenhouse gas impact. By studying more than 50 small pore Cu-exchanged zeolite SCR catalyst samples, including model catalysts synthesized in our laboratories and state-of-the-art industrial catalysts, we explored a wide range of factors affecting N 2 O formation. These factors included Cu loading, support Si/Al ratio, support topology, catalyst aging, reaction temperature and reactant feed composition effects. We probed N 2 O formation under both steady-state SCR, and during NH 4 NO 3 decomposition via temperature programmed desorption (TPD). Finally, we used DFT to probe energetics of possible N 2 O formation pathways. Based on these studies, we confirm that low temperature N 2 O formation occurs via multiple reaction pathways that all involve NH 4 NO 3 and are supported by Cu moieties that facilitate in-situ NO oxidation to NO 2 .

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

Selective oxidation of 3-hydroxypropionic acid to malonic acid over Pd/C: Mechanistic and kinetic study

Malonic acid (MA) is a high-value dicarboxylic acid with strong industrial demand, yet its current production heavily relies on petrochemical feedstocks. Here, we report the first systematic study for sustainable production of MA via oxidation of 3-hydroxypropionic acid (3-HP) with a Pd/Carbon catalyst. The effects of oxidant type (O 2 and H 2 O 2 ), pH, and temperature on the reaction chemistry were comprehensively evaluated. Guided by experimental observation and DFT-calculated thermochemical energetics, reaction networks for 3-HP oxidation with both oxidants are proposed and validated through kinetic modeling. MA was identified as the primary oxidation product, while further oxidation yields acetic acid and oxalic acid. The kinetic model validated the network, displaying excellent agreement (R 2 > 0.95). Kinetic observations also enabled a direct comparison between O 2 and H 2 O 2 and revealed their distinct behaviors. The model was further developed into a temperature-time map, providing insight into conditions that maximize MA production. Malonic acid selectivity of 56.9 % and yield of 50.5 % were achieved at 3 bar oxygen, equimolar NaOH:3-HP ratio, and 50 °C.

3-hydroxypropionic acid↗
DOE OSTI2020

Towards defect engineering in hexagonal MoS 2 nanosheets for tuning hydrogen evolution and nitrogen reduction reactions

Combined computational and experimental approaches were used to evaluate defective 2H-MoS 2 nanosheets for their activity and selectivity for hydrogen evolution reaction (HER) and nitrogen reduction reaction (NRR). Density functional theory calculations were used to understand the relationship between HER and NRR activity on the ideal basal MoS 2 plane, seven grain-boundaries, ten single-/few-atom vacancies and anti-sites, and zigzag and armchair edge sites. The results confirm that 2H-MoS 2 should contain several defects with high activity for HER: armchair and zigzag edges, VS vacancy, MoS 2 anti-site, and S-S and Mo-Mo grain boundaries. Considering Gibbs free energy change for all the steps in the NRR mechanism and kinetic barriers for a key NRR step, we have found that activation of perspective NRR selective sites in 2H-MoS 2 , namely V MoS6 and clusters of S-vacancies, would require large overpotential, conditions at which HER dominates. The DFT conclusions are supported by the electrochemical studies of NRR activity and selectivity under aqueous conditions, which show an increase in NRR activity but a decrease in Faradaic efficiency as applied cell potential becomes more negative. The results of this work therefore highlight the challenges in activating natural 2H-MoS 2 for NRR, which would require additional material engineering or reaction condition optimization as a way to suppress HER, decrease the NRR overpotential or preferentially both.

2H-MoS2↗