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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Design of super-strong and thermally stable nanotwinned Al alloys via solute synergy

Al alloys have widespread industrial applications. However, their mechanical strength is often much lower than steels. Here, we investigate the influence of solutes on achieving ultrahigh strength and thermal stability of nanotwinned Al alloys. In situ micropillar compression tests show the addition of a small amount of Ti can significantly increase the mechanical strength of Al–Ni alloys to 2 GPa. Deformation induced detwinning, Ni segregation and grain coarsening as discovered in binary Al–Ni alloys are mostly absent in the ternary Al–Ni–Ti alloys. Moreover, the ternary Al–Ni–Ti alloys have outstanding thermal stability. Density function theory calculations reveal the synergetic pinning effect of Ni–Ti solute pairs on incoherent twin boundaries. This study demonstrates that the proper selection of synergistic solute pairs is critical to improve the thermal stability and mechanical properties of nanotwinned Al alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗
DOE OSTI2024

Transition-Metal-Related Quantum Emitters in Wurtzite AlN and GaN

Transition-metal centers exhibit a paramagnetic ground state in wide-bandgap semiconductors and are promising for nanophotonics and quantum information processing. Specifically, there is a growing interest in discovering prominent paramagnetic spin defects that can be manipulated using optical methods. Here, we investigate the electronic structure and magneto-optical properties of Cr and Mn substitutional centers in wurtzite AlN and GaN. We use state-of-the-art hybrid density functional theory calculations to determine level structure, stability, optical signatures, and magnetic properties of these centers. The excitation energies are calculated using the constrained occupation approach and rigorously verified with the complete active space configuration interaction approach. Our simulations of the photoluminescence spectra indicate that $Cr$$^{1+}_{Al}$ in AlN and $Cr$$^{1+}_{Ga}$ in GaN are responsible for the observed narrow quantum emission near 1.2 eV. We compute the zero-field splitting (ZFS) parameters and outline an optical spin polarization protocol for $Cr$$^{1+}_{Al}$ and $Cr$$^{1+}_{Ga}$. Our results demonstrate that these centers are promising candidates for spin qubits.

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

Gas-phase fragmentation of single heteroatom-incorporated Co 5 MS 8 (PEt 3 ) 6 + (M = Mn, Fe, Co, Ni) nanoclusters

Functionalization of metal-chalcogenide clusters by either replacing core atoms or by tuning the ligand is a powerful technique to tailor their properties. Central to this approach is understanding the competition between the strength of the metal-ligand and metal-metal interactions. Here, using collision-induced dissociation of atomically precise metal sulfide nanoclusters, Co 5 MS 8 L 6 + (L = PEt 3 , M = Mn, Fe, Co, Ni) and Co 5-x Fe x S 8 L 6 + (x = 1–3), we study the effect of a heteroatom incorporation on the core-ligand interactions and relative stability towards fragmentation. Sequential ligand loss is the dominant dissociation pathway that competes with ligand sulfide (LS) loss. Because the ligands are attached to metal atoms, LS loss is an unusual dissociation pathway, indicating significant rearrangement of the core prior to fragmentation. Both experiments and theoretical calculations indicate the reduced stability of Co 5 MnS 8 L 6 + and Co 5 FeS 8 L 6 + towards the first ligand loss in comparison with their Co 6 S 8 L 6 + and Co 5 NiS 8 L 6 + counterparts and provide insights into the core-ligand interaction.

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

Finite-Size Scaling on a Digital Quantum Simulator Using Quantum Restricted Boltzmann Machine

The critical point and the critical exponents for a phase transition can be determined using the Finite-Size Scaling (FSS) analysis. This method assumes that the phase transition occurs only in the infinite size limit. However, there has been a lot of interest recently in quantum phase transitions occurring in finite size systems such as a single two-level system interacting with a single bosonic mode e.g., in the Quantum Rabi Model (QRM). Since these phase transitions occur at a finite system size, the traditional FSS method is rendered inapplicable for these cases. For cases like this, we propose an alternative FSS method in which the truncation of the system is done in the Hilbert space instead of the physical space. This approach has previously been used to calculate the critical parameters for stability and symmetry breaking of electronic structure configurations of atomic and molecular systems. We calculate the critical point for the quantum phase transition of the QRM using this approach. We also provide a protocol to implement this method on a digital quantum simulator using the Quantum Restricted Boltzmann Machine algorithm. Our work opens up a new direction in the study of quantum phase transitions on quantum devices.

Khalid, Bilal↗
DOE OSTI2022

Different interfacial reactivity of lithium metal chloride electrolytes with high voltage cathodes determines solid-state battery performance

A deep understanding of the interaction of the surface of cathode materials with solid electrolytes is crucial to design advanced solid-state batteries (SSBs). This is especially true for the new class of lithium metal chloride (Li-M-Cl) solid electrolytes which are receiving rapidly growing attention due to their very high oxidative stability (>4 V) in combination with good ionic conductivity that can enable long cell cycle life. While Li-M-Cl electrolytes typically contain resource-limited metals (M) such as indium or rare earths, work has focused on substituting M with more abundant elements such as zirconium. Via operando resistance measurements using intermittent current interruption we explore the dynamic evolution of the interphase at the surface of Ni-rich NCM85 or NCM111 cathode particles inside a working SSB with three different Li-(M 1 ,M 2 )-Cl catholytes (Li 3 InCl 6 , Li 2 Sc 1/3 In 1/3 Cl 4 and Li 5/2 Y 1/2 Zr 1/2 Cl 6 ) to reveal the impact of the cationic metal substitution on the interfacial chemistry. We show that the metal plays a critical role in determining high voltage stability, contrary to prior assumptions. Using a combination of cyclic voltammetry and ultraviolet photoelectron spectroscopy measurements of the electronic band structure to assess oxidative stability; coupled with DFT calculations and ToF-SIMS to evaluate products formed at the interface at different upper cutoff potentials and degrees of delithiation, we are able to differentiate between electrochemical and chemical degradation. We find that Li 2 Sc 1/3 In 1/3 Cl 4 yields the highest (and Li 3 InCl 6 the lowest) stability against electrochemical oxidation, while Li 5/2 Y 1/2 Zr 1/2 Cl 6 undergoes a detrimental chemical reaction with oxygen released from Ni-rich NCM85 at high potentials, resulting in fast capacity fading. Altogether, our work establishes a platform for the metrics and an approach that can be utilized to efficiently evaluate the stability of new halide SEs in SSB cells.

25 ENERGY STORAGE↗
DOE OSTI2022

An in-line burn-up algorithm: proposal and preliminary analysis

A new burn-up algorithm is proposed for coupling neutronics and isotopic depletion. The algorithm is analogous to those proposed for enforcing xenon equilibrium and converging neutronics and thermal-hydraulics. During the corrector step of standard predictor-corrector methods, the corrector nuclide density is updated as neutronics converges. This is expected to improve stability while reducing computational expense compared to other implicit algorithms. The stability of the algorithm is investigated by Fourier analysis for a simple case whose results have been successfully compared with numerical experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Mechanisms of Ethylene Epoxidation over Silver from Machine Learning-Accelerated First-Principles Modeling and Microkinetic Simulations

In this work we employed machine learning-augmented density functional theory (DFT) thermodynamic calculations to assess the stability of different AgO x structures under catalytic ethylene epoxidation reaction conditions. We found that there are multiple AgO x surface motifs that could co-exist under the relevant conditions. These included Ag surface oxides (e.g., AgO_p(4 × 4) and Ag 1.83 O) and atomic oxygen-covered Ag(111) surfaces. Furthermore, we employed DFT calculations to evaluate the energetics of different reaction mechanisms by which ethylene and oxygen can react on these surfaces. These studies revealed several energetically viable reaction pathways for ethylene epoxidation. Microkinetic modeling analyses, based on the DFT-calculated reaction pathways, showed that ethylene epoxidation can proceed on all surfaces and that multiple pathways, including those involving Langmuir–Hinshelwood and Eley–Rideal mechanisms, could be involved in selective and unselective reactions. The diversity of mechanisms that we discovered in the context of the relatively simple ethylene epoxidation reaction on Ag suggests that the richness and complexity of surface chemistry are most likely a rule rather than an exception in heterogeneous catalytic chemical transformations on metal surfaces and that the concept of a single or even a dominant mechanism and reaction intermediates might need to be revisited for many reactions.

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

Quasi‐Binary Transition Metal Dichalcogenide Alloys: Thermodynamic Stability Prediction, Scalable Synthesis, and Application

Abstract Transition metal dichalcogenide (TMDCs) alloys could have a wide range of physical and chemical properties, ranging from charge density waves to superconductivity and electrochemical activities. While many exciting behaviors of unary TMDCs have been demonstrated, the vast compositional space of TMDC alloys has remained largely unexplored due to the lack of understanding regarding their stability when accommodating different cations or chalcogens in a single‐phase. Here, a theory‐guided synthesis approach is reported to achieve unexplored quasi‐binary TMDC alloys through computationally predicted stability maps. Equilibrium temperature–composition phase diagrams using first‐principles calculations are generated to identify the stability of 25 quasi‐binary TMDC alloys, including some involving non‐isovalent cations and are verified experimentally through the synthesis of a subset of 12 predicted alloys using a scalable chemical vapor transport method. It is demonstrated that the synthesized alloys can be exfoliated into 2D structures, and some of them exhibit: i) outstanding thermal stability tested up to 1230 K, ii) exceptionally high electrochemical activity for the CO 2 reduction reaction in a kinetically limited regime with near zero overpotential for CO formation, iii) excellent energy efficiency in a high rate Li–air battery, and iv) high break‐down current density for interconnect applications. This framework can be extended to accelerate the discovery of other TMDC alloys for various applications.

Hemmat, Zahra↗
DOE OSTI2022

H 2 BZmacropa-NCS: A Bifunctional Chelator for Actinium-225 Targeted Alpha Therapy

Actinium-225 ( 225 Ac) is one of the most promising radionuclides for targeted alpha therapy (TAT). With a half-life of 9.92 days and a decay chain that emits four high-energy α particles, 225 Ac is well-suited for TAT when conjugated to macromolecular targeting vectors that exhibit extended in vivo circulation times. The implementation of 225 Ac in these targeted constructs, however, requires a suitable chelator that can bind and retain this radionuclide in vivo. Previous work has demonstrated the suitability of a diaza-18-crown-6 macrocyclic chelator H 2 macropa for this application. Building upon these prior efforts, in this study, two rigid variants of H 2 macropa, which contain either one (H 2 BZmacropa) or two (H 2 BZ 2 macropa) benzene rings within the macrocyclic core, were synthesized and investigated for their potential use for 225 Ac TAT. The coordination chemistry of these ligands with La 3+ , used as a nonradioactive model for Ac 3+ , was carried out. Both NMR spectroscopic and X-ray crystallographic studies of the La 3+ complexes of these ligands revealed similar structural features to those found for the related complex of H 2 macropa. Thermodynamic stability constants of the La 3+ complexes, however, were found to be 1 and 2 orders of magnitude lower than those of H 2 macropa for H 2 BZmacropa and H 2 BZ 2 macropa, respectively. The decrease in thermodynamic stability was rationalized via the use of density functional theory calculations. 225 Ac radiolabeling and serum stability studies with H 2 BZmacropa showed that this chelator compares favorably with H 2 macropa. Based on these promising results, a bifunctional version of this chelator, H 2 BZmacropa-NCS, was synthesized and conjugated to the antibody codrituzumab (GC33), which targets the liver cancer biomarker glypican-3 (GPC3). The resulting GC33-BZmacropa conjugate and an analogous GC33-macropa conjugate were evaluated for their 225 Ac radiolabeling efficiencies, antigen-binding affinities, and in vivo biodistribution in HepG2 liver cancer tumor-bearing mice. Although both conjugates were comparably effective in their radiolabeling efficiencies, [ 225 Ac]Ac-GC33-BZmacropa showed slightly poorer serum stability and biodistribution than [ 225 Ac]Ac-GC33-macropa. Altogether, these results establish H2BZmacropa-NCS as a new bifunctional chelator for the preparation of 225 Ac radiopharmaceuticals.

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

Solid breeder material crystal structure evolution due to Li burn up–Loss of crystal stability

Solid breeder materials for fusion power plants are subjected to extreme conditions of radiation damage, high temperature and designed component burn-up. Here, the loss of lattice stability of Li 2 TiO 3 due to Li burn-up was investigated using first principles-based approaches to aid evaluation of upper limit of component lifetime. The lattice stability is analyzed using calculated phonon dispersion in the Li 2 TiO 3 supercell, while Li burn-up is modeled by the introduction of Li-vacancies. It has been determined that the studied ceramic can be structurally stable up to ~40% of Li atoms burn-up. At 50% burn-up, the negative frequency branch in the phonon dispersion spectrum appears indicating the loss of lattice dynamical stability. Moreover, the structure obtained by introduction of unstable frozen phonons with minima energy amplitude also results in a phonon dispersion with negative branches in the structure obtained. Therefore, the system completely loses stability when approximately half of Li atoms are consumed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct Mechanochemical Synthesis, Phase Stability, and Electrochemical Performance of α-NaFeO 2

To better understand polymorph control in transition metal oxides, the mechanochemical synthesis of NaFeO 2 was explored. In this paper, we report the direct synthesis of α-NaFeO 2 through a mechanochemical process. By milling Na 2 O 2 and γ-Fe 2 O 3 for 5 h, α-NaFeO 2 was prepared without high-temperature annealing needed in other synthesis methods. While investigating the mechanochemical synthesis, it was observed that changing the starting precursors and mass of precursors affects the resulting NaFeO 2 structure. Density functional theory calculations on the phase stability of NaFeO 2 phases show that the α phase is stabilized over the β phase in oxidizing environments, which is provided by the oxygen-rich reaction between Na 2 O 2 and Fe 2 O 3 . This provides a possible route to understanding polymorph control in NaFeO 2 . Annealing the as-milled α-NaFeO 2 at 700 °C has resulted in increased crystallinity and structural changes that improved electrochemical performance in terms of capacity over the as-milled sample.

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

Equilipy: a python package for calculating phase equilibria

The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Stabilization of N6 and N8 anionic units and 2D polynitrogen layers in high-pressure scandium polynitrides

Abstract Nitrogen catenation under high pressure leads to the formation of polynitrogen compounds with potentially unique properties. The exploration of the entire spectrum of poly- and oligo-nitrogen moieties is still in its earliest stages. Here, we report on four novel scandium nitrides, Sc 2 N 6 , Sc 2 N 8 , ScN 5, and Sc 4 N 3 , synthesized by direct reaction between yttrium and nitrogen at 78-125 GPa and 2500 K in laser-heated diamond anvil cells. High-pressure synchrotron single-crystal X-ray diffraction reveals that in the crystal structures of the nitrogen-rich Sc 2 N 6 , Sc 2 N 8, and ScN 5 phases nitrogen is catenated forming previously unknown N 6 6 − and N 8 6 − units and$${\!\,}_{\infty }{\!\,}^{2}({{{{{\rm{N}}}}}}_{5}^{3-})$$ ∞ 2 ( N 5 3 − ) anionic corrugated 2D-polynitrogen layers consisting of fused N 12 rings. Density functional theory calculations, confirming the dynamical stability of the synthesized compounds, show that Sc 2 N 6 and Sc 2 N 8 possess an anion-driven metallicity, while ScN 5 is an indirect semiconductor. Sc 2 N 6 , Sc 2 N 8 , and ScN 5 solids are promising high-energy-density materials with calculated volumetric energy density, detonation velocity, and detonation pressure higher than those of TNT.

Science & Technology - Other Topics↗
DOE OSTI2024

Cover crops and poultry litter impact on soil structural stability in dryland soybean production in southeastern United States

Abstract This study explored the efficacy of soil aggregate indices in quantifying soil structural development, utilizing 5‐year field experiment data from the Southeastern United States. The experiment utilized a split‐plot design with cover crops (native vegetation as control, cereal rye (Secale cerealeL.), winter wheat (Triticum aestivum), hairy vetch (Vicia villosa), and mustard (Brassica rapa) plus cereal rye as the main factor and fertilizer source (no fertilizer as control, inorganic fertilizer with phosphorus, potassium, and elemental sulfur, and poultry litter) as the secondary factor. Aggregate size fractions were determined using the wet‐sieving method, and aggregate stability index (ASI), mean weight diameter (MWD), geometric mean diameter (GMD), and fractal dimension (FD) were calculated to assess soil structural stability. Main effects results indicated that cereal rye (55.11%) and poultry litter (50.97%) exhibited the highest ASI values. The highest MWD, GMD, and FD were observed under mustard plus cereal rye (1.187 mm), cereal rye (0.462 mm), and hairy vetch (2.573), respectively. Principal component analysis revealed that cover crops significantly improved soil aggregate structure and stability, overcoming limitations of sole fertilization practices. Regression analysis suggested that ASI, MWD, and GWD positively correlated with soil organic carbon, whereas FD negatively correlated with MWD, GMD, and ASI. Principal component analysis exhibited that FD decreased with increasing soil organic carbon, ASI, MWD, and GMD, demonstrating that lower FD values indicate enhanced soil aggregation and structure. Assessed indices, FD included, effectively gauged soil structural stability. These metrics should be prioritized in managerial decisions to support soil productivity and health in agricultural systems.

Agriculture↗
DOE OSTI2024

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗
DOE OSTI2026

Localizing tetrahedral aluminum in nitrate-bearing gibbsite to constrain defect-impurity coupling

The enhanced radiolytic stability of gibbsite (α-Al(OH) 3 ) containing trace nitrate (NO 3 − ) is a phenomenon in nuclear waste management, but its structural origins remain unresolved. Motivated by the detection of minority tetrahedral aluminum (T d ) defects in synthetic gibbsite, we hypothesized that these sites may participate in NO 3 − retention or mediate H 2 suppression. To evaluate this, we combined orthogonal techniques comprised of spatially selective solid-state 27 Al MAS NMR, comparative spectroscopy, and density functional theory (DFT) modeling. Paramagnetic editing and dynamic nuclear polarization (DNP) MAS NMR confirm that T d defects are confined to the particle interior. DFT calculations reveal no energetic stabilization of NO 3 − near T d sites. Comparative NMR analysis shows that T d is also present in chloride-bearing gibbsite, which exhibits high radiolytic hydrogen yields. These three independent disqualifications rule out T d as a structural contributor to nitrate-mediated suppression and narrow the scope of defect-driven explanations. The findings redirect mechanistic attention away from coordination defects and toward redox-active impurity pathways, providing a refined foundation for understanding radiation tolerance in Al(OH) 3 .

Graham, Trent R. [Pacific Northwest National Labor↗