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

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis↗
DOE OSTI2023

High-Precision Calculation of Nanoparticle (Nanocrystal) Potentials of Mean Force and Internal Energies

Thermodynamic stability assessment of nanocrystal systems requires precise free energy calculations. This study highlights the importance of meticulous control over various factors, including the thermostat, time step, potential cutoff, initial configuration, sampling method, and overall simulation duration. Free energy computations in dry (solvent-free) systems are on the order of several hundred $k_BT$ but can be obtained with consistent accuracy. However, calculation of internal energies becomes challenging, as they are typically much larger in magnitude than free energies and exhibit significant noise and reduced reliability. To address this limitation, we propose a new internal energy estimate that drastically reduces the noise. Here, we also present formulas that enable the optimization of the parameters of the harmonic bias potential for optimal convergence. Finally, we discuss the implications of these findings for the computation of free energies in nanocrystal clusters and superlattices.

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

Structure and good piezoelectric performance in the complex system of Pb[(Zn,Ni)Nb]O 3 –Pb[(In,Yb)Nb]O 3 –Pb(Zr,Hf,Ti)O 3

High-performance piezoelectrics are always demanded for the high-end application. In this study, a complex piezoelectric system of 0.49Pb(Zn 1/2 Ni 1/2 ) 1/3 Nb 2/3 O 3 – x Pb(In 1/2 Yb 1/2 ) 1/2 Nb 1/2 O 3 –(0.51 – x )Pb(Zr 1/2 Hf 1/2 ) 0.1 Ti 0.9 O 3 (0.16 ≤ x ≤ 0.23) was fabricated through the solid-state method. The structure, ferroelectric, piezoelectric, and dielectric properties were investigated. The optimum piezoelectric coefficient d 33 of 761 pC/N, high Curie temperature of 169 °C, dielectric permittivity ( ε r ) of 4557, and electromechanical coupling coefficient ( k p ) of 63% were found at the morphotropic phase boundary composition of x = 0.19, which are superior to other complex piezoelectric materials. In particular, a significant large-signal d 33 * of 913 pm/V and low strain hysteresis (6%) was obtained in the temperature range of 20–170 °C. Temperature-dependent x-ray diffraction (XRD) has demonstrated that good temperature stability is put down to the structure stability. The agreement between the calculated lattice strain from in situ high-energy synchrotron XRD data and the macroscopic measurements suggests that the large lattice strain has a dominant contribution to the high piezoelectric response. The high piezoelectric performance and good temperature stability makes it potential for application.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surface stability and H adsorption and diffusion near surfaces of W borides: a first-principles study

Abstract Understanding the behavior of tungsten boride (W x B y ) surfaces in a fusion reactor environment is an important topic since boronization is a common wall conditioning method used in fusion Tokamaks. We report the results of density functional theory calculations that investigate the surface stability of W x B y (tetragonal I4 1 /amd -WB, hexagonal P 6 3 / mmc -WB 2 and tetragonal I4/m- W 2 B) with low-index orientations, as well as hydrogen (H) energetics near W x B y surfaces. For single element terminated W x B y surfaces, B terminated surfaces are more energetically stable than W terminated as a result of significant reconstruction of B. The H surface adsorption energy and activation energy of H diffusion penetration below W x B y surfaces are mainly related to the outer termination. Specifically, the WB(001) surface terminated with two B layers, referred to as WB(001)-T BB , has higher H adsorption affinity and lower H diffusivity on this surface than other terminations, which is controlled by the significant charge transfer from B to H. However, B atoms on the WB 2 (0001)-T BB surface decrease both H adsorption and diffusivity on the surface, but enhance H diffusion below the surface in comparison to W terminated WB 2 (0001) surface. H would be trapped and diffuse within atomic surface gaps on the WB 2 ( 2 1 ˉ 1 ˉ 0 ) surface, while H below the surface layer would jump along the [0001] direction rather than diffuse into bulk. The surface diffusion activation energy of H on the W 2 B(001) surface slightly varies with terminations. Once H crosses the surface layer of W 2 B(001) with either termination, it prefers to diffuse into the bulk, or back towards the surface, rather than move parallel to the surface. Interestingly, WB 2 (0001) and WB 2 ( 2 1 ˉ 1 ˉ 0 ) surfaces will have relatively higher H retention than the other W x B y surfaces evaluated in this work.

Yang, L. (ORCID:0000000322166071)↗
DOE OSTI2025

Stacking polymorphism of PtSe 2 : its implication to layer-dependent metal-insulator transitions

Using diffusion Monte Carlo (DMC) and density functional theory (DFT) calculations, we examine the structural stability and interlayer binding properties of PtSe 2 . Our DMC study reveals that AA and AB-r bilayer stacking modes are nearly degenerate, highlighting the significant role of interlayer hybridization in offsetting the energy cost due to larger interlayer separations in the AB-r mode. Additionally, our DMC-benchmarked DFT calculations with the r 2 SCAN+rVV10 functional uncover pronounced stacking polymorphism in few-layer PtSe 2 , driven by degenerate AA and AB-r interfaces, which leads to substantial band gap variations across different stacking configurations. This polymorphism, along with selenium vacancies, influences a layer-dependent metal-insulator transition observed in few-layer PtSe 2 . Our findings emphasize the importance of both van der Waals interactions and interlayer hybridization in determining the phase stability and electronic properties of PtSe 2 , advancing our understanding of its fundamental properties and refining theoretical models for practical applications in nanoelectronic devices.

Electronic materials↗
DOE OSTI2024

Anomalous Hall crystals in rhombohedral multilayer graphene. II. General mechanism and a minimal model

Here we propose a minimal "three-patch model"for the anomalous Hall crystal (AHC), a topological electronic state that spontaneously breaks both time-reversal symmetry and continuous translation symmetry. The proposal for this state is inspired by the recently observed integer and fractional quantum Hall states in rhombohedral multilayer graphene at zero magnetic field. There, interaction effects appear to amplify the effects of a weak moiré potential, leading to the formation of stable, isolated Chern bands. It has been further shown that Chern bands are stabilized in mean-field calculations even without a moiré potential, enabling a realization of the AHC state. Our model is built on the dissection of the Brillouin zone into patches centered around high-symmetry points. Within this model, the wave functions at high-symmetry points fully determine the topology and energetics of the state. We extract two quantum geometrical phases of the noninteracting wave functions that control the stability of the topologically nontrivial AHC state. The model predicts that the AHC state wins over the topological trivial Wigner crystal in a wide range of parameters, and agrees very well with the results of full self-consistent Hartree-Fock calculations of the rhombohedral multilayer graphene Hamiltonian.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Structural and vibrational properties of methane up to 71 GPa

Single-crystal synchrotron x-ray diffraction, Raman spectroscopy, and first principles calculations have been used to identify the structure of the high-pressure (HP) phase of molecular methane above 20 GPa up to 71 GPa at room temperature. The structure of the HP phase is trigonal R3, which can be represented as a distortion of the cubic phase B, previously documented at 7–15 GPa and confirmed here. The positions of hydrogen atoms in the HP phase have been obtained from first principles calculations, which also demonstrated the stability of this structure above 260 K at 25 GPa. The molecules occupy four different crystallographic sites in phase B and 11 sites in the HP phase, which result in splitting of molecular stretching modes detected in Raman spectroscopy and assigned here based on a good agreement with the Raman spectra calculated from the first principles. Furthermore, our study points out to a single HP phase unlike up to three HP phases proposed previously based on the Raman spectroscopy results only.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

The crystal chemistry of ZrSi

Reported here is a revised crystal structure of β–ZrSi (TlI/CrB structure type), correcting the atomic position and bond distances. The Si–Si bond length has been modified substantially from 2.723(6) Å to 2.4411(8) Å. The β-ZrSi single crystals were grown from an arc-melted button, and were characterized using single crystal X-ray diffraction. A survey of the TlI/CrB structure type shows that changes to the nomenclature would be useful, separating it into four chemically distinct subtypes: TlI, CrB, CaSi, and ThCo. β–ZrSi is an example of the CaSi subtype. α-ZrSi crystallizes in the related FeB structure type, which is also divided here into the subtypes FeB, CeSi, and YNi. The effect of electron count on the relative phase stabilities of the CaSi and CeSi subtypes is rationalized with the aid of the Zintl concept and electronic structure calculations in the Linear Muffin-Tin Orbital (LMTO) basis.

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

Exploring the Coexistence of Spin States in [Fe(tpy-Ph) 2 ] 2+ Complexes on Au(111) Using DFT Calculations

In this work, we systematically study the electronic structure and stability of spin states of the [Fe-(tpy-ph) 2 ] 2+ molecule in both the gas phase and on a Au(111) substrate using density functional theory + U (DFT+ U ) calculations. We find that the stability of the Fe 2+ ion’s spin states predicted by the computations is significantly influenced by the Hubbard U parameter. In the gas phase, the low-spin (LS, S = 0) state is found to be energetically favorable for U (Fe) ≤ 3 eV, whereas the high-spin (HS, S = 2) state is stabilized for U (Fe) > 3 eV. Interaction with the Au(111) substrate is found to elevate the critical U for the spin-state transition to 3.5 eV. Additionally, we perform L-edge X-ray absorption spectroscopy (XAS) calculations for both HS and LS states. The calculated XAS suggests that the HS state more closely aligns with the experimental observations, indicating the potential coexistence of the HS state as the initial state during the X-ray excitation process. These findings enrich our understanding of spin-state dynamics in [Fe(tpy-Ph) 2 ] 2+ .

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

Atomic-scale modeling of C/N kinetic stability descriptors for PGM-free electrocatalysts at finite temperatures

The durability of platinum group metal-free (PGM-free) electrocatalysts is a major barrier to their usage in polymer electrolyte fuel cell cathodes. C and N removal from active sites may play an important role in the catalyst’s ability to maintain high activity. While C degradation mechanisms are kinetically controlled, previous studies have focused on thermodynamic descriptors. In this work, we develop a temperature-dependent kinetic descriptor of C and N stability using an electron beam-damage model. Our approach considers the electron beam energy threshold (EBET) describing the knock-on displacement of C and N atoms as a stability descriptor for atomic structures. The stability of different sites is calculated to be different showing this approach can discriminate between similar sites with varied configurations. Additionally, we provide important insight regarding TEM beam damage of proposed active sites. We calculate 60 keV electrons can damage some proposed active site structures even at room temperature.

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

Interplay between electron localization, magnetic order, and Jahn-Teller distortion dictates LiMnO2 phase stability

The development of manganese (Mn)-rich cathodes for Li-ion batteries promises to alleviate potential supply chain bottlenecks in battery manufacturing. Fundamental challenges in Mn-rich cathodes arise from phenomena such as structural changes due to cooperative Jahn-Teller (JT) distortions of in octahedral environments, Mn migration, and phase transformations to spinel-like order, all of which affect the electrochemical performance. These physically complex phenomena motivate an re-examination of the Li-Mn-O rock-salt space, with a focus on the thermodynamics of the prototypical, polymorphs. It is found that the generalized gradient approximation (GGA-PBEsol) and meta-GGA ( ) density functionals with empirically fitted on-site Hubbard corrections yield spurious stable phases for , such as predicting a phase with -like order ( ) to be the ground state instead of the orthorhombic (Pmmn) phase, which is the experimentally known ground state. Accounting for antiferromagnetic order in each structure is shown to have a substantial effect on the total energies and resulting phase stability. By using hybrid-GGA (HSE06) and GGA with self-consistent Hubbard parameters (on-site and inter-site ) calculated from linear response theory, the experimentally observed phase stability trends are recovered. The calculated on-site between Mn- states in the experimentally observed orthorhombic, layered, and spinel phases are significantly smaller than in and disordered layered structures, by within GGA. The smaller values of are shown to be correlated with a collinear ordering of JT distortions, in which all orbitals are oriented in the same direction. This cooperative JT effect can lead to greater electron delocalization from Mn along the states due to increased Mn-O covalency, which contributes to the greater electronic stability compared to the phases with noncollinear JT arrangements. The structures with collinear ordering of JT distortions also generate greater vibrational entropy, which helps stabilize these phases at high temperature. These phases are shown to be strongly insulating with large calculated band gaps , which are computed using HSE06 and .

Kam, Ronald L↗
DOE OSTI2020

Stabilizing effect of enhanced resistivity on peeling-ballooning instabilities on EAST

Previous stability analysis of NSTX equilibrium with lithium-conditioning demonstrates that enhanced resistivity due to the increased effective charge number Z eff (i.e. increased impurity level) can provide a stabilizing effect on low-n edge localized modes. This paper extends the resistivity stabilizing effect to the intermediate-n peeling-ballooning (PB) instabilities with the linear stability analysis of EAST high-confinement mode equilibria in NIMROD two-fluid calculations. However, in this work, the resistivity stabilizing effect on PB instabilities in the EAST tokamak appears weaker than that found in NSTX. This work may provide better insight into the physical mechanism behind the beneficial effects of impurities on pedestal stability.

EAST tokamak↗
DOE OSTI2024

Stabilization of U 5 f 2 configuration in UTe 2 through U 6 d dimers in the presence of Te2 chains

We investigate the topological superconductor candidate UTe 2 using high-resolution valence-band resonant inelastic x-ray scattering at the U M 4 , 5 edges. We observe atomiclike low-energy excitations that support the correlated nature of this unconventional superconductor. These excitations originate from the U 5 f 2 configuration, which is unexpected since the short Te2-Te2 distances exclude Te2 being 2 − . By utilizing the photoionization cross-section dependence of the photoemission spectra in combination with band structure calculations, we infer that the stabilization of the U 5 f 2 configuration is due to the U 6 d bonding states in the U dimers acting as a charge reservoir. Our results emphasize that the description of the physical properties should commence with a 5 f 2 ansatz. Published by the American Physical Society 2024

Christovam, Denise S. (ORCID:0000000250698107)↗
DOE OSTI2021

Dy adsorption on and intercalation under graphene on 6 H -SiC(0001) surface from first-principles calculations

Previous experimental observations motivate clarification of configuration stabilities and kinetic processes for intercalation of guest atoms into a layered van der Waals material such as a graphene-SiC system. From our first-principles density functional theory (DFT) calculations, we analyze Dy adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that intercalation of a single Dy atom into the gallery between TLG and BLG is more favorable than adsorption on TLG but that intercalation into the gallery underneath BLG is highly unfavorable. We obtain diffusion barriers of ~0.45 and 0.54 eV for a Dy atom diffusing on and under TLG, respectively. We find that the direct penetration of a Dy atom from the graphene top into the gallery under TLG is almost inhibited below a temperature of ~1400 K due to a large global barrier of at least ~3.5 eV. Instead, we find that a single Dy atom on TLG can easily intercalate by crossing a TLG step (e.g., a zigzag step presaturated by a Dy chain or a reconstructed zigzag step zz57). Additionally, we also perform DFT calculations for different Dy coverages to demonstrate how the favorability of Dy intercalation, as well as the corresponding interlayer spacings, depend on the coverage. Consequently, we can provide general insight and guidance for extensively studied systems involving intercalation of foreign atoms into graphene on a SiC substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Energetic Stability of Freestanding and Metal-Supported Borophenes: Quantum Monte Carlo and Density-Functional Theory Calculations

Quantum Monte Carlo (QMC) and density-functional theory (DFT) calculations were carried out to study cohesion energetics of two-dimensional (2D) sheets of boron atoms called borophenes. Our QMC calculations confirmed the polymorphism among free-standing borophenes that was reported in previous DFT studies. Although Perdew-Burke-Ernzerhof (PBE) calculations significantly overestimate the cohesive energies of 2D boron sheets, DFT-PBE relative energetics with respect to each other among various free-standing borophenes are found to be in quantitative agreement with the corresponding QMC results. This suggests that one can make reliable predictions for relative stability of different boron sheets through DFT-PBE calculations. Our analysis of PBE formation energies of borophenes on metal surfaces shows that the polymorphic range is extended for borophenes on the Ag(111) and the Au(111) surfaces beyond that of free-standing borophenes, reflecting recent experimental synthesis of beta(12) and chi(3) boron sheets on the Ag(111) surface. We have also found that a hexagonal borophene can be stabilized through charge transfer from a metal surface and is energetically favored on the AI(111) surface over other borophene structures. Finally, it is found that the bilayer formation could be energetically favored over its monolayer form for the borophene-Au system, especially for borophenes with hexagonal hole densities eta lower than 1/9. This leads to our prediction that in addition to its monolayer form, a bilayer eta = 1/12 borophene can be synthesized on the Au(111) surface, opening a new possibility for borophene-based electronic devices.

ab initio↗
DOE OSTI2023

Coordination Stoichiometry Effects on the Binding Hierarchy of Histamine and Imidazole–M 2+ Complexes

Histidine–M 2+ coordination bonds are a recognized bond motif in biogenic materials with high hardness and extensibility, which has led to growing interest in their use in soft materials for mechanical function. However, the effect of different metal ions on the stability of the coordination complex remains poorly understood, complicating their implementation in metal-coordinated polymer materials. Herein, rheology experiments and density functional theory calculations are used to characterize the stability of coordination complexes and establish the binding hierarchy of histamine and imidazole with Ni 2+ , Cu 2+ , and Zn 2+ . It is found that the binding hierarchy is driven by the specific affinity of the metal ions to different coordination states, which can be macroscopically tuned by changing the metal-to-ligand stoichiometry of the metal-coordinated network. These findings facilitate the rational selection of metal ions for optimizing the mechanical properties of metal-coordinated materials.

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

Structural constraint integration in a generative model for the discovery of quantum materials

Billions of organic molecules have been computationally generated, yet functional inorganic materials remain scarce due to limited data and structural complexity. Here, in this work, we introduce Structural Constraint Integration in a GENerative model (SCIGEN), a framework that enforces geometric constraints, such as honeycomb and kagome lattices, within diffusion-based generative models to discover stable quantum materials candidates. SCIGEN enables conditional sampling from the original distribution, preserving output validity while guiding structural motifs. This approach generates ten million inorganic compounds with Archimedean and Lieb lattices, over 10% of which pass multistage stability screening. High-throughput density functional theory calculations on 26,000 candidates shows over 95% convergence and 53% structural stability. A graph neural network classifier detects magnetic ordering in 41% of relaxed structures. Furthermore, we synthesize and characterize two predicted materials, TiPd 0.22 Bi 0.88 and Ti 0.5 Pd 1.5 Sb, which display paramagnetic and diamagnetic behaviour, respectively. Our results indicate that SCIGEN provides a scalable path for generating quantum materials guided by lattice geometry.

36 MATERIALS SCIENCE↗