Search NASA⌕ Search

SEARCH · Search NASA

Results for “STABILITY 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.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Strain Effects in SrHfO 3 Films Grown by Hybrid Molecular Beam Epitaxy

Perovskite oxide heterostructures host a large number of interesting phenomena such as ferroelectricity, which are often driven by octahedral distortions in the crystal that may induce polarization. SrHfO 3 (SHO) is a perovskite oxide with a pseudocubic lattice parameter of 4.08 Å that previous density functional theory (DFT) calculations suggest can be stabilized in a ferroelectric P4mm phase when stabilized with sufficient compressive strain. Additionally, it is insulating and possesses a large band gap and a high dielectric constant, making it an ideal candidate for oxide electronic devices. Here, to test the viability of epitaxial strain as a driver of ferroic phase transitions, SHO films were grown by hybrid molecular beam epitaxy (hMBE) with a tetrakis(ethylmethylamino)hafnium(IV) source on GdScO 3 and TbScO 3 substrates. Strained SHO phases were characterized using X-ray diffraction, X-ray absorption spectroscopy, and scanning transmission electron microscopy to determine the space group of the strained films, with the results compared to those of DFT-optimized models of phase stability versus strain. Contrary to past reports, we find that compressively strained SrHfO 3 undergoes octahedral tilt distortions without associated ferroelectric polarization and most likely takes on the I4/mcm phase with the a 0 a 0 c – tilt pattern.

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

Discovering Ni/Cu Single-Atom Alloy as a Highly Active and Selective Catalyst for Direct Methane Conversion to Ethylene: A First-Principles Kinetic Study

Direct methane conversion to liquid fuels or value-added chemicals is a promising technology to utilize natural resources without resorting to further petroleum extraction. However, discovering efficient catalysts for this reaction is challenging due to either coke formation or unfavorable C-H bond activation. Herein, we design single-atom alloy (SAA) catalysts to simultaneously eliminate the above two bottlenecks based on mechanism-guided strategies: (1) the active single atom enables favorable C-H bond breaking and (2) the less reactive host metal facilitates C-C coupling and thus avoids strong binding of carbonaceous species. Employing electronic structure theory calculations, we screened the stability of multiple SAAs with 3d-5d transition metals atomically dispersed on a copper surface in terms of avoiding dopant aggregation and segregation. We then evaluated reactivities of the stable SAAs as catalysts for direct methane conversion to C2 products, including methane dehydrogenation and C-C coupling mechanisms. Combining selectivity analysis with kinetic modeling, we predicted that nickel dispersed on copper, i.e., Ni/Cu SAA, is a highly active and selective catalyst that can efficiently transform methane to ethylene. This work designs efficient SAA catalysts for direct methane activation and provides chemical insights into engineering compositions of SAAs to tune their catalytic performances.

Kothakonda, Manish↗
DOE OSTI2025

The Role of Cu 3+ in the Oxygen Evolution Activity of Copper Oxides

Cu-based oxides and hydroxides represent an important class of materials from a catalytic and corrosion perspective. In this study, we investigate the formation of bulk and surface Cu 3+ species that are stable under water oxidation catalysis in alkaline media. So far, no direct evidence existed for the presence of hydroxides (CuOOH) or oxides, which were primarily proposed by theory. Here, this work directly places CuOOH in the oxygen evolution reaction (OER) Pourbaix stability region with a calculated free energy of −208.68 kJ/mol, necessitating a revision of known Cu–H 2 O phase diagrams. We also predict that the active sites of CuOOH for the OER are consistent with a bridge O* site between the two Cu 3+ atoms with onset at ≥1.6 V vs the reversible hydrogen electrode (RHE), aligning with experimentally observed Cu 2+/3+ oxidation waves in cyclic voltammetry of Fe-free and Fe-spiked copper in alkaline media. Trace amounts of Fe (2 μg/mL (ppm) to 5 μg/mL) in the solution measurably enhance the catalytic activity of the OER, likely due to the adsorption of Fe species that serve as the active sites . Importantly, modulation excitation X-ray absorption spectroscopy (ME-XAS) of a Cu thin-film electrode shows a distinct Cu 3+ fingerprint under OER conditions at 1.8 V vs RHE. Additionally, in situ Raman spectroscopy of polycrystalline Cu in 0.1 mol/L (M) KOH revealed features consistent with those calculated for CuOOH in addition to CuO. Overall, this work provides direct evidence of bulk electrochemical Cu 3+ species under OER conditions and expands our longstanding understanding of the oxidation mechanism and catalytic activity of copper.

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

Colloidal silver diphosphide (AgP 2 ) nanocrystals as low overpotential catalysts for CO 2 reduction to tunable syngas

Production of syngas with tunable CO/H 2 ratio from renewable resources is an ideal way to provide a carbon-neutral feedstock for liquid fuel production. Ag is a benchmark electrocatalysts for CO 2 -to-CO conversion but high overpotential limits the efficiency. We synthesize AgP 2 nanocrystals (NCs) with a greater than 3-fold reduction in overpotential for electrochemical CO 2 -to-CO reduction compared to Ag and greatly enhanced stability. Density functional theory calculations reveal a significant energy barrier decrease in the formate intermediate formation step. In situ X-ray absorption spectroscopy (XAS) shows that a maximum Faradaic efficiency is achieved at an average silver valence state of +1.08 in AgP 2 NCs. A photocathode consisting of a n + p-Si wafer coated with ultrathin Al 2 O 3 and AgP 2 NCs achieves an onset potential of 0.2 V vs. RHE for CO production and a partial photocurrent density for CO at –0.11 V vs. RHE ( j –0.11, CO ) of –3.2 mA cm –2 .

Heterogeneous catalysis↗
DOE OSTI2022

Double-Bilayer polar nanoregions and Mn antisites in (Ca, Sr)3Mn2O7

Abstract The layered perovskite Ca 3 Mn 2 O 7 (CMO) is a hybrid improper ferroelectric candidate proposed for room temperature multiferroicity, which also displays negative thermal expansion behavior due to a competition between coexisting polar and nonpolar phases. However, little is known about the atomic-scale structure of the polar/nonpolar phase coexistence or the underlying physics of its formation and transition. In this work, we report the direct observation of double bilayer polar nanoregions (db-PNRs) in Ca 2.9 Sr 0.1 Mn 2 O 7 using aberration-corrected scanning transmission electron microscopy (S/TEM). In-situ TEM heating experiments show that the db-PNRs can exist up to 650 °C. Electron energy loss spectroscopy (EELS) studies coupled with first-principles calculations demonstrate that the stabilization mechanism of the db-PNRs is directly related to an Mn oxidation state change (from 4+ to 2+), which is linked to the presence of Mn antisite defects. These findings open the door to manipulating phase coexistence and achieving exotic properties in hybrid improper ferroelectric.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Single-atom Rh/N-doped carbon electrocatalyst for formic acid oxidation

To meet the requirements of potential applications, it is of great importance to explore new catalysts for formic acid oxidation that have both ultra-high mass activity and CO resistance. In this work, we successfully synthesize atomically dispersed Rh on N-doped carbon (SA-Rh/CN) and discover that SA-Rh/CN exhibits promising electrocatalytic properties for formic acid oxidation. The mass activity shows 28- and 67-fold enhancements compared with state-of-the-art Pd/C and Pt/C, respectively, despite the low activity of Rh/C. Interestingly, SA-Rh/CN exhibits greatly enhanced tolerance to CO poisoning, and Rh atoms in SA-Rh/CN resist sintering after long-term testing, resulting in excellent catalytic stability. Density functional theory calculations suggest that the formate route is more favourable on SA-Rh/CN. According to calculations, the high barrier to produce CO, together with the relatively unfavourable binding with CO, contribute to its CO tolerance.

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

Computational discovery of promising new n-type dopable ABX Zintl thermoelectric materials

Computational prediction of good thermoelectric (TE) performance in several n-type doped Zintl phases, combined with successful experimental realization, has sparked interest in discovering new n-type dopable members of this family of materials. However, most known Zintls are typically only p-type dopable; prior successes in finding n-type Zintl phases have been largely serendipitous. Here in this paper, we go beyond previously synthesized Zintl phases and perform chemical substitutions in known n-type dopable ABX Zintl phases to discover new ones. We use first-principles calculations to predict their stability, potential for TE performance as well as their n-type dopability. Using this approach, we find 17 new ABX Zintl phases in the KSnSb structure type that are predicted to be stable. Several of these newly predicted phases (KSnBi, RbSnBi, NaGeP) are found to exhibit promising n-type TE performance and are n-type dopable. We propose these compounds for further experimental studies, especially KSnBi and RbSnBi, which are both predicted to be good TE materials with high electron concentrations due to self-doping by native defects, when grown under alkali-rich conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Higher loadings of Pt single atoms and clusters over reducible metal oxides: application to C–O bond activation

We develop higher loadings of isolated noble metal atoms and clusters on a mildly reducible metal oxide. Here we demonstrate the approach for Pt supported on TiO 2 and confirmed it by XRD, AC-HAADF-STEM, CO-FTIR, XAS, and XPS. Density functional theory calculations rationalize the experimental stability and the IR shifts using mixtures of CH 3 I and CO. The redispersed catalysts are thermally stable in inert gas or H 2 and afford enhanced selectivity and activity in hydrodeoxygenation reactions compared to metal nanoparticles by creating surface oxygen vacancies that promote C–O cleavage without side reactions. Higher metal loadings, e.g., 1%Pt/TiO 2 , on the oxide surface profoundly increase the activity of the bare oxide catalyst tenfold compared to ultra-low loadings typically used for single atom catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modulation of the coordination environment enhances the electrocatalytic efficiency of Mo single atoms toward water splitting

Here, enhancing the catalytic efficiency through engineering active site environments is expected to work pronouncedly for single atom catalysts (SACs) because of intense atomic scale interactions involved between SAs and their coordination environments. Taking Mo SACs for catalyzation of the hydrogen evolution reaction (HER) as an example, three SACs of different coordination environments, namely Mo-O 2 N 2 , Mo-O 2 N 1 C 1 , and Mo-O 2 C 2 , were successfully created for demonstration. The HER performances are in an increasing order of Mo-O 2 N 2 , Mo-O 2 N 1 C 1 , and Mo-O 2 C 2 , exhibiting η 10 of 98, 71, and 61 mV, η 500 of 340, 248, and 200 mV, Tafel slopes of 95.8, 39.6, and 33.8 mV dec -1 , and current density decays of 9, 6, and 6% after a 50 hour operation at an initial current density of 100 mA cm -2 , respectively. Substituting C with N in the coordination environment results in inferior HER catalytic efficiency and stability. Density functional theory calculations reveal that replacing carbon with nitrogen for coordination with the Mo SA on a carbon substrate of a higher N-doping level shifts the d-band center of Mo more negatively from the Fermi level, thereby increasing the hydrogen adsorption energy and thus decelerating the hydrogen desorption kinetics, giving consequent inferior HER activities.

25 ENERGY STORAGE↗
DOE OSTI2024

Phase transitions in HfO2 probed by first-principles computations

Ever since ferroelectricity was discovered in HfO2, the question of its origin remains controversial. Here, we probe this question using a combination of Landau theory of phase transitions and first-principles computations. In such an approach, the energy landscape associated with the phase transition between cubic and different experimentally demonstrated phases of HfO2 (tetragonal, monoclinic, orthorhombic Pbca, orthorhombic Pnma, and orthorhombic Pca21) is explored using density functional theory calculations. Computations revealed that stabilization of all but orthorhombic Pbca phase is driven by a single unstable zone-boundary antipolar mode X2−. When coupled with zone-center modes (Γ1+ and Γ3+), it stabilizes the tetragonal phase. Coupling with four additional modes (Γ5+, X3−, X5−, X5+) results in the monoclinic phase, which is the ground state of the material. If, however, Γ5+ mode is replaced with Γ4− mode, orthorhombic polar phase Pca21 is stabilized. The application of this framework to examine the effect of electric field on the ferroelectric phase of hafnia reveals that the field of 5 MV/cm is capable of stabilizing ferroelectric phase over the monoclinic one at 0 K.

Physics↗
DOE OSTI2020

Quantifying Large Lattice Relaxations in Photovoltaic Devices

Temporal variations of Cu(In,Ga)Se 2 photovoltaic device properties during light exposure at various temperatures and voltage biases for times up to 100 h are analyzed using the kinetic theory of large lattice relaxations. Open-circuit voltage and p-type doping increased with charge injection and decreased with temperature at low injection conditions. Furthermore, lattice relaxation can account for both trends and activation energies extracted from the data are approximately 0.9 and 1.2 eV for devices with lower and higher sodium content, respectively. In these devices, increased sodium content resulted in higher initial p-type doping with greater stability. First-principles calculations providing revised activation energies for the (V Se -V Cu ) complex suggest that this defect does not account for the metastability observed here.

14 SOLAR ENERGY↗
DOE OSTI2022

First-principles study of the T center in silicon

The T center in silicon is a well-known carbon-based color center that has been recently considered for quantum technology applications. Using first-principles computations, we show that the excited state is formed by a defect-bound exciton made of a localized defect state occupied by an electron to which a hole is bound. The localized state is of strong carbon p character and reminiscent of the localization of the unpaired electron in the ethyl radical molecule. The radiative lifetime for the defect-bound exciton is calculated to be on the order of microseconds, much longer than for other quantum defects such as the nitrogen vacancy center in diamond and in agreement with experiments. The longer lifetime is associated with the small transition dipole moment as a result of the very different nature of the localized and delocalized states forming the defect-bound exciton. Finally, we use first-principles calculations to assess the stability of the T center. Here, we find the T center to be stable against decomposition into simpler defects when keeping the stoichiometry fixed. However, we identify that the T center is easily prone to (de)hydrogenation and so requires very precise annealing conditions (temperature and atmosphere) to be efficiently formed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Initial conditions for simulations of beam physics in linear induction accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the DualAxis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Correct Initial Conditions for Simulations of Beam Physics in Linear Induction Accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Machine learning prediction of the mechanical properties of refractory multicomponent alloys based on a dataset of phase and first principles simulation

In this work, a dataset including structural and mechanical properties of refractory multicomponent alloys was developed by fusing computations of phase diagram (CALPHAD) and density functional theory (DFT). The refractory multicomponent alloys, also named refractory complex concentrated alloys (CCAs) which contain 2–5 types of refractory elements were constructed based on Special Quasi-random Structure (SQS). The phase of alloys was predicted using CALPHAD and the mechanical property of alloys with stable and single body-centered cubic (BCC) at high temperature (over 1,500°C) was investigated using DFT-based simulation. As a result, a dataset with 393 refractory alloys and 12 features, including volume, melting temperature, density, energy, elastic constants, mechanical moduli, and hardness, were produced. To test the capability of the dataset on supporting machine learning (ML) study to investigate the property of CCAs, CALPHAD, and DFT calculations were compared with principal components analysis (PCA) technique and rule of mixture (ROM), respectively. It is demonstrated that the CALPHAD and DFT results are more in line with experimental observations for the alloy phase, structural and mechanical properties. Furthermore, the data were utilized to train a verity of ML models to predict the performance of certain CCAs with advanced mechanical properties, highlighting the usefulness of the dataset for ML technique on CCA property prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The theory of kinetic effects on resistive wall mode stability in tokamaks

Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The effect of elastic strains on the adsorption energy of H, O, and OH in transition metals

The influence of elastic strains on the adsorption of H, O, and OH on the (111) surfaces of 8 fcc (Ni, Cu, Pd, Ag, Pt, Au, Rh, Ir) and on the (0001) surfaces of 3 hcp (Co, Zn, Cd) transition metals was analyzed by means of density functional theory calculations. To this end, surface slabs were subjected to different strain states (uniaxial, biaxial, shear, and a combination of them) up to strains dictated by the mechanical stability limits indicated by phonon calculations. It was found that the adsorption energy followed the predictions of the d-band theory but – surprisingly – the variations in the adsorption energy only depended on the area of the adsorption hole and not on the particular elastic strain tensor applied to achieve this area. The analysis of the electronic structure showed that the applied strains did not modify the shape of the Projected Density of States (PDOS) of the d-orbitals of the transition metals but only led to a shift in the energy levels. Moreover, the presence of the adsorbates on the surfaces led to negligible changes in the PDOS. Thus, the adsorption energies were a function of the Fermi energy which in turn was associated with the change of the area of the adsorption through a general linear law that was valid for all metals. Furthermore, the information in this paper allows the immediate and accurate estimation of the effect of any elastic strain on the adsorption energies of H, O, and OH in 11 transition metals with more than half-filled d-orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗