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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The microscopic mechanisms of high temperature oxidation of Haynes 282

Nickel-based superalloy finds widespread applications in aerospace and extreme environments. They are known for their high temperature oxidation resistance under extended periods. However, the oxide formation and evolution which sets the stages for the later parabolic oxidation kinetics is not fully understood. This paper aims to provide new insights into the transient stage oxidation mechanism and kinetics of a typical Ni-based superalloy (Haynes 282). While tracking the chemical and microstructural evolution under micrometer scale at 800°C, we show that the oxide scale and its grain boundary species change significantly during the initial stages of oxidation and can have a profound impact on the oxidation kinetics. Cr 2 O 3 forms initially at the grain boundary along with minor amount of Al 2 O 3 . Then, the grain boundary region is enriched with copious amounts TiO 2 while Cr migrates away from the grain boundary. A noticeable change in the isothermal oxidation kinetics observed likely results from a mechanistic change from uniform surface oxidation to preferential outward diffusion of Ti 4+ ions through the grain boundary. Through first-principles calculations combined with energy dispersive spectroscopy, we confirmed the preferential outward diffusion of Ti 4+ ions as the diffusion barrier is lower for Ti than Cr along the grain boundaries. In conclusion, these findings highlight the critical role of early-stage grain boundary oxidation dynamics in dictating the long-term oxidation resistance of Ni-based superalloys and provide a foundation for future strategies to enhance their performance in extreme environments.

DFT↗
DOE OSTI2021

$\mathrm{Perturbo}$: A software package for ab initio electron–phonon interactions, charge transport and ultrafast dynamics

We report Perturbo is a software package for first-principles calculations of charge transport and ultrafast carrier dynamics in materials. The current version focuses on electron–phonon interactions and can compute phonon-limited transport properties such as the conductivity, carrier mobility and Seebeck coefficient. It can also simulate the ultrafast nonequilibrium electron dynamics in the presence of electron–phonon scattering. Perturbo uses results from density functional theory and density functional perturbation theory calculations as input, and employs Wannier interpolation to reduce the computational cost. It supports norm-conserving and ultrasoft pseudopotentials, spin–orbit coupling, and polar electron–phonon interactions for bulk and 2D materials. Hybrid MPI plus OpenMP parallelization is implemented to enable efficient calculations on large systems (up to at least 50 atoms) using high-performance computing. Taken together, Perturbo provides efficient and broadly applicable ab initio tools to investigate electron–phonon interactions and carrier dynamics quantitatively in metals, semiconductors, insulators, and 2D materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

VTAnDeM: A python toolkit for simultaneously visualizing phase stability, defect energetics, and carrier concentrations of materials

Phase stability, defect formation energies, and carrier concentrations are closely interrelated features of semiconductors. Due to their joint dependence on the multidimensional chemical potential space, it is challenging to quantitatively establish patterns between these quantities in a given semiconductor, especially when the semiconductor is comprised of multiple elements. To enable synchronous visualization and analysis of these complementary material properties and their interdependence, we developed the Visualization Toolkit for Analyzing Defects in Materials (VTAnDeM). This python-based toolkit allows users to interactively explore how defect formation energies and carrier concentrations vary across the composition and chemical potential spaces of multicomponent semiconductors. Here, we illustrate the computational workflow that employs VTAnDeM as a post-processing tool for first-principles calculations and describe the data organization and theory underlying the visualization scheme. Furthermore, we believe that this software will serve as a useful tool for simultaneously visualizing the often complex and non-intuitive chemical potential – defect – carrier concentration phase space of semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanistic understanding of pH effects on the oxygen evolution reaction

The oxygen-evolution reaction (OER) is pivotal in many energy-conversion technologies as it is an important counter reaction to others that convert stable chemicals to higher-value products using electrochemistry. The local microenvironment and pH for the anode OER can vary from acidic to neutral to alkaline depending on the system being explored, making definitive mechanistic insights difficult. In this paper, we couple experiments, first-principles calculations based on density functional theory, microkinetics, and transport modeling to explore the entire pH range of the OER. At low current densities, neutral pH values unexpectedly perform better than the acidic and alkaline conditions, and this trend is reversed at higher current densities (> 20 mA cm -2 ). Using multiscale modeling, this switch is rationalized by a change from a dual-reaction mechanism to a single rate-determining step. The model also shows how the alkaline reaction rates dominate in the middle to high pH range. Furthermore, we explore that the local pH for near-neutral conditions is much different (e.g., 2.4 at the reaction surface vs. 9 in the bulk) than the pH extremes, demonstrating the criticality that transport phenomena plays in kinetic activity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Continuum-Level Modeling of Li-Ion Battery SEI by Upscaling Atomistically Informed Reaction Mechanisms

Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.

continuum-level model↗
DOE OSTI2025

Thermal decomposition pathways of bulk electrolytes on vanadium oxide nanocrystals

The thermal stability of electrolytes at an elevated temperature induced by battery charge-discharge cycling is critical for the long cycling performance of a rechargeable battery. For many multivalent systems, such as rechargeable magnesium batteries, which offer great potential for high energy density and utilize earth-abundant resources, electrolyte instability and electrode surface passivation, arising from electrolyte decomposition, remain as major roadblocks. Understanding the electrolyte decomposition pathways at the electrode-electrolyte interface is essential to provide guidance in overcoming this challenge. In this work, in situ 13 C magic angle spinning nuclear magnetic resonance (MAS NMR) and first-principles calculations were used to investigate the thermal decomposition of the electrolyte in a system consisting of MgV 2 O 4 , a novel cathode for magnesium batteries, mixed with a bulk electrolyte consisting of magnesium bis(trifluoromethanesulfonyl)imide (Mg(TFSI) 2 ) in diglyme (G2). We show that significant electrolyte decomposition is observed in bulk 1.0 M Mg(TFSI) 2 in G2 mixed with nanometer sized MgV 2 O 4 powder at elevated temperatures. This observation is to mimic the possible thermal decomposition that might happen during battery cycling. We demonstrate that the MgV 2 O 4 surface is covered by a layer of decomposed G2 products. We conclude that the dominant reaction pathway for electrolyte decomposition is the thermal decomposition of the pure electrolytes at elevated temperatures, followed by adsorption of G2 decomposition products to the MgV 2 O 4 surface. The activation energy for the major decomposition pathway is obtained. In conclusion, this work highlights the importance of studying thermal decomposition of electrolytes for overall system stability and explores electrolyte stability at significantly elevated temperatures.

25 ENERGY STORAGE↗
DOE OSTI2024

Computing the solubility of argon and xenon in molten sodium chloride and potassium chloride salts

Molten salt reactors (MSRs) offer significant advancements in nuclear reactor safety and efficiency by operating at higher temperatures and lower pressures compared to traditional reactors. A critical aspect of MSR operation involves understanding the solubility of fission byproducts, particularly noble gases, in the molten salts used. This study employs molecular dynamics (MD) simulations to compute Henry’s law constants and enthalpies of solvation for argon and xenon in molten sodium chloride (NaCl) and potassium chloride (KCl). We developed a new pairwise potential for the noble gas and salt interactions based on first principles calculations. We then used this potential to calculate Henry’s law constants of the two gases in the molten salts, which were modeled using both a rigid ion model (RIM) and a polarizable ion model (PIM). The solubility calculations, performed using the Widom insertion method, show qualitative agreement with limited experimental data, highlighting the temperature dependence and greater solubility of both gases in KCl compared to NaCl. Additionally, free volume analysis elucidated the role of available space within the molten salts in governing solubility trends. Our findings suggest that PIM trajectories provide more reliable predictions for noble gas solubility than RIM due to their accurate density representation. Furthermore, these results enhance understanding of gas solubility in MSR environments, and the methods can be readily extended to other systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Anomalous elastic behavior of tantalum at high pressures: Experimental and theoretical studies

Elastic wave velocities of polycrystalline tantalum have been measured up to 13.6 GPa at room temperature using ultrasonic interferometry technique in a multi-anvil apparatus. The bulk and shear moduli, as well as their pressure derivatives are obtained using Eulerian 3rd order finite strain equations, yielding K S0 = 193.9(29) GPa, G 0 = 69.3(10) GPa, K S0' = 3.18(5) and G 0' = 0.88(1). Here the Steinberg-Guinan yield stress model is examined using the current shear modulus and its pressure derivative, where a slightly lower A 0 = 0.0127 GPa –1 compared to previous studies is obtained. First-principles calculations on the elasticity of tantalum have also been performed up to 200 GPa. Anomalous softening of the elastic shear constant C 44 and shear wave velocity V S is observed between 80– 180 GPa, whereas the anisotropy shows a softening-stiffening behavior with increasing pressure and reaches a minimum at around 150 GPa. Such anomalies might commonly originate from the nesting of Fermi surface for Group VB metals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Primary solidification of ternary compounds in Al-rich Al–Ce–Mn alloys

Primary solidification of ternary compounds Al 20 Mn 2 Ce and Al 10 Mn 2 Ce were analyzed through the coupling of the thermodynamic modeling and classic nucleation theory. Thermodynamic models of Al 20 Mn 2 Ce and Al 10 Mn 2 Ce were developed using the CALPHAD approach based on first-principles calculated enthalpy of formation and experimental data obtained from this work and the literature. The analysis suggested that despite the larger thermodynamic driving force for nucleation of Al 10 Mn 2 Ce, nucleation is preferred for the Al 20 Mn 2 Ce phase in the highly undercooled liquid due to its smaller interfacial energy. Therefore, manufacturing methods with rapid cooling rates will favor primary solidification of Al 20 Mn 2 Ce for Al-rich Al–Ce–Mn alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Exploring the elastic and electronic properties of chromium molybdenum diboride alloys

We perform first-principles calculations to study the structural, mechanical, thermal, electronic, and magnetic properties of Cr 1-x Mo x B 2 for x = 0.25, 0.33, 0.50, 0.67 and 0.75. Based on structural search methods, we determine the ground-state structure for each concentration. The ternaries are either monoclinic (x = 0.25, 0.75) or trigonal (x = 0.33, 0.50, 0.67). The calculated mechanical properties reveal that the strength of Cr 1-x Mo x B 2 is maximized for x = 0.50. Cr 0.5 Mo 0.5 B 2 exhibits excellent mechanical properties (B = 298 GPa, Y = 558 GPa, G = 235 Gpa, ν = 0.19, H v =27 GPa), surpassing those of β-MoB 2 at a lower cost. All of these ternaries are hard alloys with Vickers hardness greater than 24 GPa. Chemical bonding analysis demonstrates that the strength of the new compounds is related to the alternating planar and buckled B-B layers, as well as the strong TM-B bonds. Finally, the enhanced strength of Cr 0.5 Mo 0.5 B 2 is a consequence of the high density of strong interlayer Cr-Mo metallic bonds around the Fermi level.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strain induced electronic structure, and magnetic and structural properties in quaternary Heusler alloys ZrRhTiZ (Z = Al, In)

In this work, electronic structure, and magnetic and structural properties of quaternary Heusler alloys ZrRhTiZ (Z = Al, In) have been analyzed from first principles calculations. The ferromagnetic ground state and thermodynamical stability of these compounds are supported by relative total energies, and by derived formation and cohesive energies, respectively. The investigated Curie temperatures of the compounds exceed room temperature indicating that these compounds are promising candidates for beyond room temperature spintronics and magneto-electronics applications. Both compounds follow the Slater-Pauling 18 electron rule and exhibit half-metallic ferromagnetic behavior with magnetic moments of 2 μB with 100% spin polarization in one spin channel at the Fermi-level. The effect of uniform strain on half-metallicity of these compounds inverts the band gap from one spin channel to another spin channel with a mixed spin (magnetic moment) regime that can be mapped to an exchange bias like effect observed in some of the magnetic Heusler compounds. The half metallic ferromagnetic behavior is retained within –2 to 2% strain contrasting with the anomalous spin flip in –1% and 1% strain in ZrRhTiZ (Z = Al, In).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Additively manufactured Ni-20Cr to V functionally graded material: Computational predictions and experimental verification of phase formations

Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Angle-resolved polarized Raman study of layered b-As x P x-1 alloys: Identification of As-P vibrational modes

In this study, the polarization-resolved Raman spectra of b-As x P x-1 flakes with varying arsenic concentration (x = 0, 0.4, and 0.8) are systematically investigated. A clear polarization dependence is observed in the Raman intensity of all the P-P, As-P, and As-As modes and confirmed by the first-principles calculations. The observed angle-dependence can be related to the structure of the Raman tensors of each peak. In particular, the Raman intensities corresponding to As-P vibration modes are studied in detail. The out of plane As-P mode ($A$$^{1}_{g}$) is clearly identified, but the butterfly-shaped behavior observed in the polar plot around 350 cm -1 is interpreted as the combination of the two in-plane vibrations, one along the armchair and the other along the zigzag direction. This work allows unambiguous assignment of the $A$$^{1}_{g}$, B 2g , and $A$$^{2}_{g}$ modes of the As-P vibrational modes and the importance of the consideration of the superposition of the in-plane B 2g and $A$$^{2}_{g}$ modes. Such study can be utilized to further characterize the anisotropic structural properties and analyze the local distribution of the As and P in b-As x P x-1 alloys which will play key roles in electronic and optical properties. Furthermore, this methodology can be further extended to alloyed two-dimensional material systems with in-plane anisotropy in general.

2D materials↗
DOE OSTI2020

Investigation of the reaction mechanism of the hydrodeoxygenation of propionic acid over a Rh(1 1 1) surface: A first principles study

Microkinetic models based on first principles calculations have been used to study the vapor and liquid phase hydrodeoxygenation of propionic acid on a Rh(111) surface. Calculations suggest that both decarboxylation and decarbonylation do not occur at an appreciable rate under all reaction environments. Propanol and propionaldehyde are the main products on this surface and they are produced at similar rates in both vapor and liquid phase environments. While a condensed phase can shift the reaction rate, the dominant pathways and selectivity towards the various products are hardly affected. At 473 K, the turnover frequency is increased by about a factor of 1.5 in liquid water relative to the gas phase. In liquid 1,4-dioxane, the turnover frequency is also slightly increased relative to the gas phase. Here, given the uncertainty in Rh cavity radius in the liquid phase calculations, computations with different cavity radii have been performed. With larger Rh cavity radius, the promotional effect of the solvents on the turnover frequency becomes more significant, while practically no changes are observed for a smaller cavity radius.

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

Selective dehydra-decyclization of cyclic ethers to conjugated dienes over zirconia

ZrO 2 provides high selectivity (>90%) to conjugated pentadienes through the dehydra-decyclization of C-5 cyclic ethers, even at high conversions. 1,3-Pentadiene was the major product in both reaction of 2-methyltetrahydrofuran and tetrahydropyran over ZrO 2 . The reaction of 3-methyltetrahydrofuran produced nearly stoichiometric amounts of isoprene. Other catalysts, including TiO 2 , γ-Al 2 O 3 , and H-ZSM-5, were generally much less selective and produced a mixture of diene isomers. A combination of TPD and steady-state measurements revealed that both piperylenes are exclusively produced through primary catalytic pathways from 2-methyltetrahydrofuran, avoiding any isomerization once formed. First-principle calculations on ZrO 2 imply the presence of an energetically favored, surface isomerization of ring-opened intermediates to conjugated alkenolates that selectively dehydrate to conjugated dienes, providing high selectivity to the desired products. Finally, the stabilization of the conjugated alkenolate is key for understanding the ability of ZrO 2 to selectively produce conjugated dienes from cyclic ethers, without the need for the thermo-limited diene isomerization.

1,3 -pentadiene↗
DOE OSTI2025

Tunable broadband luminescence in lead-free hybrid copper halides

Metal halides are an important class of optoelectronic materials combining exceptional optical and electronic properties. An inherent advantage of metal halides is their solution synthesis and processability, which render them as low-cost and environmentally friendly materials for a range of applications from photovoltaics and photodetection to solid-state lighting (SSL). Here, in this study, we synthesized three previously unreported lead-free organic–inorganic hybrid copper halides: (OA) 4 CuX 5 (X = Br, I; OA + = C 8 H 17 NH 3+ , n-octylammonium cation) and (HA) 2 CuI 3 (HA + = C 6 H 13 NH 3 + , n-hexylammonium cation), all of which exhibit broadband emissions arising from self-trapped excitons (STEs). Among these compounds, (OA) 4 CuI 5 demonstrates tunable dual-band white-light emission with a high color rendering index value of 91 at room temperature. Temperature-dependent photoluminescence measurements and first-principles calculations reveal distinct behaviors between the two emission states in (OA) 4 CuI 5 . These findings highlight the potential of copper halide compounds for optoelectronic applications, particularly in the development of environmentally friendly solid-state lighting technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reinvestigation of the intrinsic magnetic properties of (Fe 1–$x$ Co $x$ ) 2 B alloys and crystallization behavior of ribbons

In this paper, new determination of the magnetic anisotropy from single crystals of ( Fe 1 - x Co x ) 2 B alloys are presented. The anomalous temperature dependence of the anisotropy constant is discussed using the standard Callen-Callen theory, which is shown to be insufficient to explain the experimental results. A more material specific study using first-principles calculations with disordered moments approach gives a much more consistent interpretation of the experimental data. Since the intrinsic properties of the alloys with x = 0.3 - 0.35 are promising for permanent magnets applications, initial investigation of the extrinsic properties are described, in particular the crystallization of melt spun ribbons with Cu, Al, and Ti additions. Previous attempts at developing a significant hysteresis have been unsuccessful in this system. Our melt-spinning experiment indicates that this system shows rapid crystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronic, structural and magnetic properties of Mn (1+x) Pt (1-x) Sb

Electronic and magnetic properties of half-metallic Heusler alloys can be modified by tuning their chemical compositions. We have carried out a combined theoretical and experimental investigation of Mn (1+x) Pt (1-x) Sb (0 ≤ x ≤ 0.5) alloys. Our first-principles calculations indicate that the stoichiometric MnPtSb exhibits nearly half-metallic band structure, but a robust half-metallicity can be achieved in Mn-rich compositions Mn (1+x) Pt (1-x) Sb (0 ≤ x ≤ 0.5) with x = 0.25 and higher in their cubic structures. In addition, while MnPtSb exhibits ferromagnetic alignment, Mn (1+x) Pt (1-x) Sb are ferrimagnetic for all non-zero values of x. We have also synthesized cubic MnPtSb and Mn 1.25 Pt 0.75 Sb alloys using arc melting and annealing. The magnetic properties of these alloys are consistent with our theoretical predictions. Furthermore, these results indicate that the Mn-rich Mn (1+x) Pt (1-x) Sb alloys have potential for spin-transport-based devices.

36 MATERIALS SCIENCE↗