Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear Stability Theory”

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.

44 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Canted antiferromagnetism and spin reorientation in corner-shared single chain quasi-one-dimensional Ba 2 ⁢FeSe 3

Here, we report the canted antiferromagnetic (AFM) structure together with a spin reorientation in a single chain quasi-one-dimensional (Q-1D) iron chalcogenide Ba 2⁢ FeSe 3 . Ba 2 ⁢FeSe 3 crystallizes in Pnma (No. 62) orthorhombic structure with linear single iron chains consisting of corner-shared distorted FeSe 4 tetrahedra along the 𝑏 axis. Ba 2 ⁢FeSe 3 is a narrow-gap semiconductor and orders AFM below 60 K. Modeling of neutron powder diffraction data reveals a canted AFM ground state of magnetic space group 𝑃⁢𝑎⁢21/𝑐 (BNS No. 14.80) with commensurate propagation vector 𝐤 =(0, $\frac{1}{2}$, 0), where the Fe ion spins are AFM aligned with up-down-up-down (↑−↓−↑−↓) sequence along the Q-1D chain direction of the 𝑏 axis. In the magnetically ordered state, the canting of magnetic moments reorients from the 𝑎⁢𝑐 plane to the 𝑎⁢𝑏 plane below 30 K, with a 10° tilting angle toward the 𝑎 axis, and the magnetic moment does not induce a net moment in either orientation. The density functional theory results indicate that an ↑−↓−↑−↓ AFM state is stabilized along the chain direction. In this work, we elucidate the unique canted AFM of the iron chalcogenide and pave the way for searching exotic physics in Q-1D Ba 2⁢ FeSe 3 .

Gao, Fei [Univ. of Texas at Dallas, Richardson, TX↗
DOE OSTI2024

Parametric Instability of Alfvén Waves and Wave Packets in Periodic and Open Systems

The parametric decay instability of Alfvén waves has been widely studied, but few investigations have examined wave packets of finite size and the effect of different boundary conditions on the growth rate. In this paper, we perform a linear analysis of circular and arc-polarized wave trains and wave packets in periodic and open boundary systems in a low- β plasma. We find that both types of wave are 3–5 times more stable in open boundary conditions compared to periodic. Additionally, once the wave packet width ℓ becomes smaller than the system size L , the growth rate decreases nearly with a power law γ ∝ ℓ / L . This study demonstrates that the stability of a pump wave cannot be separated from the laboratory settings, and that the growth rate of daughter waves depends on the conditions downstream and upstream of the pump wave and on the fraction of volume it fills. Our results can explain simulations and experiments of localized Alfvén waves. They also suggest that Alfvénic fluctuations in the solar wind, including sharp impulses known as switchbacks, can be more stable than traditional theory suggests depending on wind conditions.

Alfven waves↗
DOE OSTI2025

Structural Evolution of Mixed-Addenda Keggin Polyoxometalate Anions with Atom-by-Atom Substitution

Polyoxometalates (POMs) are molecular metal oxides with distinctive electronic properties that make them promising materials for applications in energy, sensors, and memory devices. One of the most promising methods of tuning the stability, photochromic, redox, and electron-spin properties of POMs is through the substitution of the metal “addenda” atoms that, along with oxygen, constitute their cage-like structures. Because traditional synthesis methods typically produce a distribution of POMs, the isolation and characterization of multimetallic POMs with predetermined stoichiometry remains challenging. The presence of multiple energetically accessible isomers further complicates the experimental characterization and theoretical modeling of multimetallic POMs. Herein, we leverage the distinguishing mass-selection capabilities of ion soft landing to prepare stoichiometrically selected Keggin PMo x W 12-x O 40 3- (x = 0 – 6, 8, 10, and 12) POMs on self-assembled monolayer surfaces free of the solvent molecules and counterions that often confound characterization of complex species at interfaces. The structures of the supported POMs are characterized with atom-by-atom precision using in situ infrared (IR) reflection absorption spectroscopy complemented by detailed density functional theory calculations. Our joint experimental and theoretical results reveal an almost linear shift in the positions of the IR bands towards lower wavenumbers with an increase in the number of lighter molybdenum atoms compared to heavier W atoms in PMoxW 12-x O 40 3- . The theoretical calculations also indicate that numerous isomeric structures may be populated at the experimental conditions and, consequently, contribute to the overall IR spectra. In conclusion, our findings indicate that in addition to the number of substituted addenda atoms and the presence of multiple isomeric structures, interactions with the surface play an important role in determining the IR spectra and structure of supported bimetallic POMs.

Prabhakaran, Venkateshkumar [Pacific Northwest Nat↗
DOE OSTI2024

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deep Koopman Neural Network for Analyzing High-Energy-Density Simulations of Electrical Wire Explosions

Megaampere-scale electrical wire experiments (EWEs) provide a platform for studying magnetohydrodynamic (MHD) instability growth in magneto-inertial fusion (MIF) devices. Even when nonlinear simulations of these experiments can digitally reproduce much of the experimentally observed instability growth, interpreting the results and understanding mode growth and evolution can be non-trivial. As a first step toward providing better interpretation of these simulation features, this work investigates the use of a deep neural network that uses Koopman operator theory to analyze the dynamics of pulsed-power-driven explosions of EWEs. This deep neural network is trained on 1-D resistive MHD simulations of EWEs. This neural network learns to transform the nonlinear data into a lower-dimensional representation where the time dynamics are linear. Layers of this neural network are shown to learn features of the simulations, including the locations of shock waves and different physical regimes of the simulation. Using the learned features, the network can compress a time state of the simulation consisting of 5120 data point into a 36-parameter lower-dimensional latent space embedding. Furthermore, these embeddings are shown to be clustered in the latent space by initial radius and time state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Non–linear bonding trends in maleonitrile-1,2–dithiolate complexes of the transuranium actinides

The trivalent actinides are produced in the nuclear fuel cycle during power production and provide the largest long-term radiation dose in used nuclear fuel. It is ideal for these elements to be removed from used nuclear fuel for disposal and a necessity for fuel recycling. A key challenge to this is the similarity of chemical behavior of the trivalent actinides to the lanthanides that are also present as fission products in used fuel. Thus far, some of the most effective separations of actinides from lanthanides utilise chelating agents containing sulfur moieties such as dithiophosphinates that selectively bind to actinide ions because of a greater bond covalency relative to lanthanide ions. Typically, greater differences between actinide and lanthanide ions are observable the more ligands and chelators bonds have a covalent character. Here, a series of complexes of the trivalent actinides Np(III) through Cf(III) (excluding Bk(III)) with maleonitrile-1,2-dithiolate (mnt 2– ) are synthesized along with their lanthanide counterparts (La(III) – Nd(III), Sm(III) – Gd(III), Dy(III)), in order to characterize the nature of chemical bonds with these metal ions and a polarizable, non-innocent, sulfur-donor ligand. The metal-sulfur bonds in these complexes trend shorter than measured for lanthanides with equivalent ionic radii. However, particularly large deviations are observed in the neptunium and plutonium complexes in both structure and bonding, resulting in a nonlinear bond length trendline for the actinide series. Density Functional Theory (DFT) calculations with Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Order (NBO) analyses indicate that for the neptunium and plutonium complexes, the presence of increased 5f-orbital participation, energy degeneracy of the metal and ligand orbitals, and the structure packing result in shortened M–S bonds. The stabilization of the energy of the 5f-orbitals and the decrease in f-contribution to bonding orbitals in the later actinides results in structural properties more similar to the lanthanide complexes.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Neutrino flavor instabilities in neutron star mergers with moment transport: Slow, fast, and collisional modes

Determining where, when, and how neutrino flavor oscillations must be included in large-scale simulations of hot and dense astrophysical environments is an enduring challenge that must be tackled to obtain accurate predictions. Here, using an angular moment-based linear stability analysis framework, we examine the different kinds of flavor instabilities that can take place in the context of the postprocessing of a neutron star merger simulation, with a particular focus on the collisional flavor instability and a careful assessment of several commonly used approximations. First, neglecting anisotropies of the neutrino field, we investigate the extent to which commonly used monoenergetic growth rates reproduce the results obtained from a full multienergy treatment. Contrary to the large discrepancies found in core-collapse supernova environments, we propose a simple combination of energy-averaged estimates that reproduces the multienergy growth rates in our representative simulation snapshot. We then quantify the impact of additional physical effects, including nuclear many-body corrections, scattering opacities, and the inclusion of the vacuum term in the neutrino Hamiltonian. Finally, we include the neutrino distribution anisotropies, which allows us to explore, for the first time in a multienergy setting, the interplay between collisional, fast, and slow modes in a moment-based neutron star merger simulation. We find that, despite a dominance of the fast instability in most of the simulation volume, certain regions exhibit only a collisional instability, while others, especially at large distances, exhibit a slow instability that is largely underestimated if anisotropic effects are neglected.

neutrino oscillations↗