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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Study of Alfvén eigenmode stability in Quasi-Poloidal Stellarator (QPS) plasma using a Landau closure model

The aim of this study is to analyze the linear stability of Alfvén eigenmodes (AE) in the QPS device heated by a tangential neutral beam injector (NBI). The analysis is performed using the gyro-fluid code FAR3d, that solves the reduced MHD equations for the thermal plasma coupled with moments of the kinetic equation for the energetic particles (EP). The AE stability is calculated in several operational regimes of the tangential NBI: EP β between 0.001 and 0.1, EP energy between 12 and 180 keV and different radial locations of the beam. The analysis is performed for vacuum and finite β equilibria as well as QPS configurations with two and three periods. The EP β threshold in the vacuum case is 0.001 and the AE frequency is lower as the energy of the EP population decreases. Toroidal Alfvén eigenmodes with f = 80–120 kHz and elliptical AE between f = 120–350 kHz are triggered between the middle-outer plasma region ($r/a \gt 0.5$). The AE stability improves in the simulations with finite β equilibria and three period configurations with respect to the vacuum case with two periods because the continuum gaps are slender, leading to a higher threshold of the EP β, above 0.03 for the AEs triggered by the helical mode families. Helical effects are not strong enough to destabilize Helical Alfvén eigenmodes, the AEs with the largest growth rates are triggered by the n = 1 and n = 2 toroidal families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Study of the Alfven eigenmodes stability in CFQS plasma using a Landau closure model

The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP β threshold. Also, if the thermal β of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. Overall, the simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Enhanced ionic conductivity and lack of paddle-wheel effect in pseudohalogen-substituted Li argyrodites

Superionic conductors are key to the development of safe and high-energy-density, all-solid-state batteries. Using a combined theoretical and experimental approach, we explore the feasibility of increasing the ionic conductivity through pseudohalogen substitution in the Li argyrodite structure. Under the guidance of calculated thermodynamic stability, BH 4 -substituted Li argyrodite, Li 5.91 PS 4.91 (BH 4 ) 1.09 , was successfully synthesized via a mechanochemical method. As-synthesized BH 4- substituted Li argyrodite displays an ionic conductivity of 4.8 mS/cm at 25°C. Ab initio molecular dynamics simulation trajectory analysis was used to investigate how BH 4 facilitates Li-ion diffusion and indicates only a weak correlation with the B–H bond motion. Here, we find that the enhanced conductivity mainly originates from the weak interaction between Li and BH 4 and find no evidence of a paddle-wheel effect from the polyanion. This work provides insight into how cluster ions enhance Li diffusion and systematically describes how to explore superionic conductors with pseudohalogen substitution.

25 ENERGY STORAGE↗
DOE OSTI2024

THD-C Sheet: A Novel Nonbenzenoid Carbon Allotrope with Tetra-, Hexa-, and Dodeca-Membered Rings

Here, we propose a novel two-dimensional carbon-based structure with tetra-, hexa-, and dodeca-membered rings, which we refer to by the abbreviated name, THD-C. The structure presents a mixture of sp–sp 2 hybridization and can potentially be synthesized by the topological assembly of 4-ethynyldiphenylacetylene molecules. By employing first-principles calculations, the stability and ease of synthesis of the sheet are investigated and compared with various C-allotropes. We predict its metallic behavior and excellent kinetic and dynamic stability. Due to the crystal structure of the sheet, a strong mechanical anisotropy is observed. The effects of functionalization on the electronic properties of the material are also studied, and different semiconducting systems are obtained. The potential of THD-C for energy storage in metal-based batteries, hydrogen storage, and catalysis is also investigated, and we find a superior performance in comparison to graphite and other allotropes. The quantum confinement effect is investigated by constructing nanoribbons and nanotubes of various sizes. For ribbons, we find that tailor-made electronic and magnetic properties can be obtained and explored in potential spintronic devices. Additionally, we observe that nanotubes are conducting irrespective of their chirality and can potentially be used for capture, storage, and separation of industrially relevant small gas molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stiff neural ordinary differential equations

Neural Ordinary Differential Equations (ODEs) are a promising approach to learn dynamical models from time-series data in science and engineering applications. This work aims at learning neural ODEs for stiff systems, which are usually raised from chemical kinetic modeling in chemical and biological systems. We first show the challenges of learning neural ODEs in the classical stiff ODE systems of Robertson’s problem and propose techniques to mitigate the challenges associated with scale separations in stiff systems. We then present successful demonstrations in stiff systems of Robertson’s problem and an air pollution problem. The demonstrations show that the usage of deep networks with rectified activations, proper scaling of the network outputs as well as loss functions, and stabilized gradient calculations are the key techniques enabling the learning of stiff neural ODEs. The success of learning stiff neural ODEs opens up possibilities of using neural ODEs in applications with widely varying time-scales, such as chemical dynamics in energy conversion, environmental engineering, and life sciences.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Electronic and Vibrational Contributions to the Bulk Stabilities of Trivalent 3d Transition Metal Oxyhydroxides from Electronic Structure Calculations

Any understanding of the thermodynamic driving forces active in the global cycles of earth abundant metals requires attendant knowledge of the bulk mineral free energies. This can be particularly challenging for metal (oxyhydr)oxide phases which at low temperature aqueous conditions can be produced at a variety of hydration states and particle sizes. Here we present high quality density functional theory calculations of the relative bulk stabilities of a geologically important family of polymorphic trivalent metal oxyhydroxide crystal structures MOOH (M=Cr, Al, Fe). To improve understanding of their relative stabilities we permute the metal cations into the various possible canonical structure types including goethite (isostructural to diaspore and bracewillite), lepidocrocite (isosctructural to boehmite), guyaniite, and grimaldiite. Electronic, vibrational, and magnetic contributions are considered on a consistent theoretical footing using hybrid functionals and converged plane-wave bases. We find the considered bulk structures to be thermodynamically degenerate within the precision of the available theory, despite the fact that the vibrational and magnetic contributions vary substantially among the cations. A significant conclusion is that Cr3+ supports strong hydrogen bonds, low O-H stretching frequencies and facile H movement relative to the other cations.

Prange, Micah P.↗
DOE OSTI2022

Effects of thin metal contacts on few-layer van der Waals ferrielectric CuInP 2 S 6

Out-of-plane polarized ferroelectric materials in a capacitive structure provide a key component for several technological applications. Furthermore, two-dimensional materials are expected to aid in the quest for both ultrathin and flexible electronics. Of the various two-dimensional ferroelectrics with out-of-plane polarization, CuInP 2 S 6 is special in that the Cu atoms are highly mobile and it has been shown to possess both low- and high-polarization states. Here, using density-functional-theory calculations, we explore the stabilization of the ferroelectric state for several prototypical metal contacts (Gr, Ni, Cu, Au, and Ag). In all cases, we find that the ferroelectric state can be stabilized at fewer layers than in the freestanding case. For all of the considered conventional metal contacts, we also find the existence of a quasi-ferroelectric state that stabilizes a polar phase for thicknesses greater than two layers of CIPS. In the cases of Au and Ag, interfacial alignment and strain can be used to stabilize ferroelectricity at the bilayer limit. Furthermore, we find that the strength of the interaction between the contact and CuInP 2 S 6 also leads to stabilization of the high-polarization state when ferroelectricity is stabilized. Lastly, energy-barrier calculations show that the system is still switchable in the presence of contact doping from the metal contacts.

2D materials↗
DOE OSTI2023

Structural Stability of Graphene-Supported Pt Layers: Diffusion Monte Carlo and Density Functional Theory Calculations

In this study, we used a combination of diffusion Monte Carlo and density functional theory calculations to investigate the stability and interlayer binding of various layered structures of Pt atoms adsorbed on graphene. Our findings show that vertically buckled Pt monolayer and bilayer with (111)-packing order are more energetically favorable than the corresponding buckled (100) or flat (100)-packing structures. This can be attributed to the significant lattice mismatch (>10%) between pristine graphene and a free-standing (100)-packing Pt layer. Additionally, our calculations reveal that among the (100)-packing Pt layers, an incommensurate structure with a Pt/C atomic ratio less than 1/2 may be more stable than the commensurate structures registered at the bridge sites, which aligns with recent experimental findings of incommensurate (100)-packing Pt layers on graphene. The interlayer binding between the Pt layer and graphene is found to be primarily driven by van der Waals interaction, except for the AA-stacked buckled-(100) Pt bilayer where the bottom Pt atoms show chemisorption to the graphene surface. In conclusion, this research offers a comprehensive examination of the stability and interlayer binding of metallic Pt layers, providing valuable insights for potential applications as a next-generation catalyst.

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

Effects of substitutional Ga on the phonons of δ -phase Pu based on density functional theory calculations

The delta phase of Pu is stabilized by Ga doping, but the mechanism of this stabilization remains an open question. Density functional theory calculations focused on how Ga doping affects the phonons sheds some light on the phonons' contribution to the stabilization. Here, the calculated phonon modes of Ga-doped delta phase Pu fall into two distinct types: localized, high frequency Ga-dominated phonon modes, and Pu-dominated modes at lower frequencies. Increasing the Ga concentration has an effect on the Pu-dominated phonon modes opposite to that of compression: higher-frequency modes soften, and lower-frequency modes stiffen. The latter provides an indication that the stabilization mechanism is not due to a thermodynamic contribution from the phonons. Furthermore, the stiffened phonon modes include candidate modes that describe possible pathways into low-temperature phases, suggesting that doping with Ga could impede such pathways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Including the parallel mass flow in calculating the steady-state solutions and stability of the momentum balance equations for a quasisymmetric stellarator

The Helically Symmetric Experiment (HSX) is a quasisymmetric stellarator with minimal parallel viscous damping in a helical direction. The parallel flow (Vǁ) along the magnetic field is similarly weakly damped by viscosity. In this paper, the self-consistent steady-state parallel and poloidal momentum balance equations are used to show that a large Vǁ on the order of the ion thermal velocity can increase the ion resonant radial electric field (Er) beyond the value calculated using the typical approximation that Vǁ is zero. By altering the damping of Vǁ, either by degrading the quasisymmetry or varying the neutral density, the ion resonant Er can shift in a controllable fashion. It is shown explicitly that there exist stable and unstable steady-state solutions in the two-dimensional space of Vǁ and Er. A stability analysis of each solution is performed by calculating the eigenvalues and eigenvectors of the Jacobian. The unstable solution corresponds to a saddle point in which the eigenvalues have opposite signs. The analysis leads to the conclusion that unstable solutions occur when the derivative of the total poloidal damping with respect to Er is positive. A hysteresis in Er and Vǁ is observed when the radial current density is linearly increased to a maximum and then decreased back to zero. Jumps in the radial electric field and the parallel flow are observed as the radial current density drives the evolution from one stable point to the next. This result is similar to experimental data observed on several devices.

Physics↗
DOE OSTI2023

Role of Anions in Stabilizing the [Zn–Al] Layered Double Hydroxides: A Thermodynamic Study

Here, room-temperature acid solution calorimetry, high-temperature oxide melt solution calorimetry, and low-temperature heat capacity measurements were employed to calculate the thermodynamic stabilities of the [Zn–Al–X] layered double hydroxides (LDH) containing different anions (X = Cl – , CO 3 2– , and SO 4 2– ). Cryogenic heat capacity measurements demonstrated a Schottky-type anomaly in the heat capacity of all three LDHs below 11 K. This anomaly is attributed to the tunneling of protons between adjacent oxygen atoms in the LDH interlayer as this creates an energy system similar to a two-level system modeled with a Schottky term. These heat capacity measurements were also used to determine vibrational entropies which, when combined with configurational entropies, provide standard entropies of these LDHs. Enthalpies of formation of LDHs from binary components were determined and combined with the entropies of formation to calculate Gibbs free energies. Based on these values, the order of stability is [Zn–Al–SO 4 ] > [Zn–Al–CO 3 ] > [Zn–Al–Cl]. This trend results from a combination of the interlayer spacing, amount of water in the interlayer, interactions among the interlayer species, and interactions between the metal hydroxide layer and the interlayer.

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

Neoclassical toroidal viscosity torque prediction via deep learning

GPECnet is a densely connected neural network that has been trained on GPEC data, to predict the plasma stability, neoclassical toroidal viscosity (NTV) torque, and optimized 3D coil current distributions for desired NTV torque profiles. Using NTV torque, driven by non-axisymmetric field perturbations in a tokamak, can be vital in optimizing pedestal performance by controlling the rotation profile in both the core, to ensure tearing stability, and the edge, to avoid edge localized modes (ELMs). The generalized perturbed equilibrium code (GPEC) software package can be used to calculate the plasma stability to 3D perturbations and the NTV torque profile generated by applied 3D magnetic fields. These calculations, however, involve complex integrations over space and energy distributions, which takes time to compute. Initially, GPECnet has been trained solely on data representative of the quiescent H-mode (QH) scenario, in which neutral beams are often balanced and toroidal rotation is low across the plasma profile. Lastly, this work provides the foundation for active control of the rotation shear using a combination of beams and 3D fields for robust and high performance QH mode operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Tuning Two-Dimensional Phthalocyanine Dual Site Metal–Organic Framework Catalysts for the Oxygen Reduction Reaction

Metal-organic frameworks (MOFs) offer an interesting opportunity for catalysis, particularly for metal-nitrogen-carbon (M-N-C) motifs by providing an organized porous structural pattern and well-defined active sites for the oxygen reduction reaction (ORR), a key need for hydrogen fuel cells and related sustainable energy technologies. Here, in this work, we leverage electrochemical testing with computational models to study the electronic and structural properties in these systems and their relationship to ORR activity and stability based on dual transitional metal centers. These consists of two M1 metals with amine nodes coordinated to a single M2 metal with a phthalocyanine linker, where M1/M2 = Co, Ni, or Cu. Co-based metal centers, in particular Ni-Co, demonstrate the highest overall activity of all nine tested MOFs. Computationally, we identify the dominance of Co-sites, relative higher importance of the M2 site, and the role of layer M1 interactions on the ORR activity. Selectivity measurements indicate that M1 sites of MOFs, particularly Co, exhibits lowest (< 4%), and Ni demonstrates highest (>46%) two-electron selectivity, in good agreement with computational studies. Direct in-situ stability characterization, measuring dissolved metal ions, and calculations, using an alkaline stability metric, confirm that Co is the most stable metal in the MOF, while Cu exhibits notable instability at the M1. Overall, this study reveals how atomistic coupling of electronic and structural properties affects the ORR performance of dual site MOF catalysts and opens new avenues for tunable design and future development of these systems for practical electrochemical applications.

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

The solution structures and relative stability constants of lanthanide–EDTA complexes predicted from computation

Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La 3+ , Eu 3+ , Gd 3+ , and Lu 3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln–EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ~7 and ~11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln–EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln 3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ~11), as well as for a single Ln 3+ ion (La) with varying ligands (EDTA at pH ~7 and ~11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln–EDTA complex resulting in a reduction of the complex stability. As a result, increased Ln–ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature.

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

Out-of-plane ferroelectricity and robust magnetoelectricity in quasi-two-dimensional materials

Thin-film ferroelectrics have been pursued for capacitive and nonvolatile memory devices. They rely on polarizations that are oriented in an out-of-plane direction to facilitate integration and addressability with complementary metal-oxide semiconductor architectures. The internal depolarization field, however, formed by surface charges can suppress the out-of-plane polarization in ultrathin ferroelectric films that could otherwise exhibit lower coercive fields and operate with lower power. Here, we unveil stabilization of a polar longitudinal optical (LO) mode in the n = 2 Ruddlesden–Popper family that produces out-of-plane ferroelectricity, persists under open-circuit boundary conditions, and is distinct from hyperferroelectricity. Our first-principles calculations show the stabilization of the LO mode is ubiquitous in chalcogenides and halides and relies on anharmonic trilinear mode coupling. We further show that the out-of-plane ferroelectricity can be predicted with a crystallographic tolerance factor, and we use these insights to design a room-temperature multiferroic with strong magnetoelectric coupling suitable for magneto-electric spin-orbit transistors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Photoemission from hybrid states of Cl@C 60 before and after a stabilizing charge transfer

Photoionization calculations of the endofullerene molecule Cl@C 60 with an open-shell chlorine atom are performed in the linear response density functional theory based on a spherical jellium model of C 60 . Cross sections for atom–fullerene hybrid photoemission studied show the effects of the hybridization symmetry, the giant plasmon and the molecular cavity. Comparisons with the results of Ar@C 60 provide insights into the role of a shell-closing electron and its influence on the dynamics. The results for Cl@C 60 are further compared with those of a more stable configuration that results after a C 60 electron transfers to Cl forming Cl–@${{\text{C}}_{60}}^{+}$. This comparison reveals noticeable differences in the ionization properties of the antibonding hybrid state while the bonding hybrid remains nearly unaltered showing a magnification covering the entire giant plasmon energy range. Furthermore, these predictions should raise interest in probing the molecular configuration using photoelectron spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Impact of Anions and Water Content on [Li–Al] Layered Double-Hydroxide Stability

[Li–Al] layered double hydroxides (LDHs) are compounds with potential as sorbents for lithium extraction from brine solutions. Here, in this work, heat capacities were measured from approximately 2.5 to 300 K for six [Li–Al] LDHs with differing anions (Cl – , OH – , and SO 4 2– ) and water content (denoted A for air-dried and O for oven-dried). These measurements were used to calculate the standard entropy at 298.15 K, and the results were combined with previously performed enthalpy measurements to calculate Gibbs energies of formation from the binary compounds. The calculated order of stability based on Gibbs energies of formation was Cl-LDH-O > OH-LDH-O > Cl-LDH-A > SO 4 -LDH-O > SO 4 -LDH-A > OH-LDH-A. Results support previous findings that higher water content generally raises the Gibbs energy of the LDH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗