Search NASA⌕ Search

SEARCH · Search NASA

Results for “NEUTRAL ELEMENT”

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 397 records · Page 22

NASA Tech Briefs, September 2003

Topics include: Oxygen-Partial-Pressure Sensor for Aircraft Oxygen Mask; Three-Dimensional Venturi Sensor for Measuring Extreme Winds; Swarms of Micron-Sized Sensors; Monitoring Volcanoes by Use of Air-Dropped Sensor Packages; Capacitive Sensors for Measuring Masses of Cryogenic Fluids; UHF Microstrip Antenna Array for Synthetic- Aperture Radar; Multimode Broad-Band Patch Antennas; 164-GHz MMIC HEMT Frequency Doubler; GPS Position and Heading Circuitry for Ships; Software for Managing Parametric Studies; Software Aids Visualization of Computed Unsteady Flow; Software for Testing Electroactive Structural Components; Advanced Software for Analysis of High-Speed Rolling-Element Bearings; Web Program for Development of GUIs for Cluster Computers; XML-Based Generator of C++ Code for Integration With GUIs; Oxide Protective Coats for Ir/Re Rocket Combustion Chambers; Simplified Waterproofing of Aerogels; Improved Thermal-Insulation Systems for Low Temperatures; Device for Automated Cutting and Transfer of Plant Shoots; Extension of Liouville Formalism to Postinstability Dynamics; Advances in Thrust-Based Emergency Control of an Airplane; Ultrasonic/Sonic Mechanisms for Drilling and Coring; Exercise Device Would Exert Selectable Constant Resistance; Improved Apparatus for Measuring Distance Between Axles; Six Classes of Diffraction-Based Optoelectronic Instruments; Modernizing Fortran 77 Legacy Codes; Active State Model for Autonomous Systems; Shields for Enhanced Protection Against High-Speed Debris; Scaling of Two-Phase Flows to Partial-Earth Gravity; Neutral-Axis Springs for Thin-Wall Integral Boom Hinges.

Source record↗

Acid Sulfate Alteration on Mars

A variety of mineralogical and geochemical indicators for aqueous alteration on Mars have been identified by a combination of surface and orbital robotic missions, telescopic observations, characterization of Martian meteorites, and laboratory and terrestrial analog studies. Acid sulfate alteration has been identified at all three landing sites visited by NASA rover missions (Spirit, Opportunity, and Curiosity). Spirit landed in Gusev crater in 2004 and discovered Fe-sulfates and materials that have been extensively leached by acid sulfate solutions. Opportunity landing on the plains of Meridiani Planum also in 2004 where the rover encountered large abundances of jarosite and hematite in sedimentary rocks. Curiosity landed in Gale crater in 2012 and has characterized fluvial, deltaic, and lacustrine sediments. Jarosite and hematite were discovered in some of the lacustrine sediments. The high elemental abundance of sulfur in surface materials is obvious evidence that sulfate has played a major role in aqueous processes at all landing sites on Mars. The sulfate-rich outcrop at Meridiani Planum has an SO3 content of up to 25 wt.%. The interiors of rocks and outcrops on the Columbia Hills within Gusev crater have up to 8 wt.% SO3. Soils at both sites generally have between 5 to 14 wt.% SO3, and several soils in Gusev crater contain around 30 wt.% SO3. After normalization of major element compositions to a SO3-free basis, the bulk compositions of these materials are basaltic, with a few exceptions in Gusev crater and in lacustrine mudstones in Gale crater. These observations suggest that materials encountered by the rovers were derived from basaltic precursors by acid sulfate alteration under nearly isochemical conditions (i.e., minimal leaching). There are several cases, however, where acid sulfate alteration minerals (jarosite and hematite) formed in open hydrologic systems, e.g., in Gale crater lacustrine mudstones. Several hypotheses have been suggested for the aqueous formation of sulfate-bearing phases under acidic conditions on the surface of Mars including (1) sulfuric acid weathering of basaltic materials; (2) oxidative weathering of ultramafic igneous rocks containing sulfides; (3) acid fog weathering of basaltic materials, and (4) near-neutral pH subsurface solutions rich in Fe2(+) that were rapidly oxidized to Fe3(+), which produced excess acidity as iron was oxidized on exposure to O2 or photo-oxidized by ultraviolet radiation at the martian surface. Next, we briefly describe evidence for these hypothesis.

Ming, D. W.↗

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗

Stability analysis of a reinforced carbon carbon shell

This paper presents the development of a stability analysis for the nose cap of the NASA Space Shuttle Orbiter. Stability is evaluated by the differential stiffness analysis of the NASTRAN finite-element computer code, addressing those nonstandard characteristics in the nose cap such as nonuniform curvature, asymmetrical and nonuniform loads, support fixity, and various combinations of membrane and bending stresses. A full-sized nose cap, thinner than production, was statically tested and stability analyzed. The failing load level correlated to within 30%. The region and mode of buckling that occurred during test was accurately predicted by analysis. The criterion for predicting instability is based on the behavior of the nonlinear deflections. The deflections are nonlinear elastic in that the stresses are well within the elastic range of the material, but the geometry-load relationship produces nonlinear deflections. The load-deflection relationship is well defined by differential stiffness analysis up to the zero-slope portion of the curve, the point of neutral stability or where the shell 'snaps through' just prior to general instability.

Agan, W. E.↗

Potentially Habitable Ancient Environments in Gusev Crater, Mars

Habitable environments must sustain liquid water at least intermittently and also provide both chemical building blocks and useful sources of energy for life. Observations by Spirit rover indicate that conditions have probably been too dry to sustain life, at least since the emplacement of the extensive basalts that underlie the plains around the Columbia Memorial Station landing site. Local evidence of relatively minor aqueous alteration [1] probably occurred under conditions where the activity of water was too low to sustain biological processes as we know them. In contrast, multiple bedrock units in West Spur and Husband Hill in the Columbia Hills have been extensively altered. Patterns of elemental abundances are consistent with aqueous processes involving migrating fluids [2]. Fe in several of these units has been extensively oxidized [3]. Conceivably any microbiota present during the aqueous alteration of these rocks might have obtained energy from Fe oxidation. Spirit discovered olivine-rich ultramafic rocks during her descent from Husband Hill southward into Inner Basin. Alteration of similar ultramafic rocks on Earth can yield H2 that can provide both energy and reducing power for microorganisms. Spirit's discovery of deposits rich in ferric sulfate is consistent with the aqueous dissolution and/or alteration of olivine under acidic conditions [2] such as those associated with hydrothermal activity. The oxidation of iron and sulfur that can accompany such activity can be an energy source for life. Hydrothermal systems on Earth that sustain either acidic [4] or neutral to alkaline fluids [5] have been shown to provide this energy. Collectively the observations by Spirit rover are consistent with the possibility that habitable environments existed in Gusev crater at least intermittently in the distant geologic past.

Fe oxidation↗

Search for heavy neutral leptons decaying into muon-pion pairs in the MicroBooNE detector

We present upper limits on the production of heavy neutral leptons (HNLs) decaying to $\mu \pi$ pairs using data collected with the MicroBooNE liquid-argon time projection chamber (TPC) operating at Fermilab. This search is the first of its kind performed in a liquid-argon TPC. We use data collected in 2017 and 2018 corresponding to an exposure of $2.0 \times 10^{20}$ protons on target from the Fermilab Booster Neutrino Beam, which produces mainly muon neutrinos with an average energy of $\approx 800$ MeV. HNLs with higher mass are expected to have a longer time-of-flight to the liquid-argon TPC than Standard Model neutrinos. The data are therefore recorded with a dedicated trigger configured to detect HNL decays that occur after the neutrino spill reaches the detector. We set upper limits at the $90\%$ confidence level on the element $\lvert U_{\mu4}\rvert^2$ of the extended PMNS mixing matrix in the range $\lvert U_{\mu4}\rvert^2<(6.6$-$0.9)\times 10^{-7}$ for Dirac HNLs and $\lvert U_{\mu4}\rvert^2<(4.7$-$0.7)\times 10^{-7}$ for Majorana HNLs, assuming HNL masses between $260$ and $385$ MeV and $\lvert U_{e 4}\rvert^2 = \lvert U_{\tau 4}\rvert^2 = 0$.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Delamination growth analysis in quasi-isotropic laminates under loads simulating low-velocity impact

A geometrically nonlinear finite-element analysis has been developed to calculate the strain energy released by delaminating plates during impact loading. Only the first mode of deformation, which is equivalent to static deflection, was treated. Both the impact loading and delamination in the plate were assumed to be axisymmetric. The strain energy release rate in peeling, GI, and shear sliding, GII, modes were calculated using the fracture mechanics crack closure technique. Energy release rates for various delamination sizes and locations and for various plate configurations and materials were compared. The analysis indicated that shear sliding was the primary mode of delamination growth. The analysis also indicated that the midplane (maximum transverse shear stress plane) delamination was more critical and would grow first before any other delamination of the same size near the midplane region. The delamination growth rate was higher (neutrally stable) for a low toughness (brittle) matrix and slower (stable) for high toughness matrix. The energy release rate in the peeling mode, GI, for a near-surface delamination can be as high as 0.5GII, and can contribute significantly to the delamination growth.

Shivakumar, K. N.↗

Ratio of Sodium to Potassium in the Mercurian Exosphere

Sodium (Na) and Potassium (K) atoms can be seen in the exosphere of Mercury and the Moon because they are extremely efficient at scattering sunlight. These species must be derived from surface materials, so that we might expect the ratio of sodium to potassium to reflect the ratio of these elements in the surface crust. This expectation is approximately born out for the Moon, where the ratio of sodium to potassium in the lunar exosphere averages to be about 6, not too far from the ratio in lunar rocks of 2 to 7. However, the ratio in the Mercury exosphere was found to be in the range 80 to 190, and at least once, as high as 400. The sodium and potassium atoms seen in the Mercury exosphere represent a balance between production from the surface and loss to space. Only if the production efficiencies and loss rates for Na and K were equal, would the ratio of Na to K in the exosphere reflect the ratio in the surface rocks. Since a value of 100 or more for the ratio of sodium to potassium in the surface rocks seems very unlikely, the high values of the observed ratios suggests that either production efficiencies or loss processes for the two elements are not equivalent. It does not seem likely that source processes should be different on the Moon and Mercury by an order of magnitude. This suggests that loss processes rather than source processes are the cause of the difference between the two. The major loss processes for sodium and potassium on Mercury are radiation pressure and trapping of photoions by the solar wind. Radiation pressure can reach 50-70% of surface gravity, and can sweep sodium and potassium atoms off the planet, provided they are sufficiently hot. Photoionization followed by trapping of the ions in the solar wind is the other major loss process. Photoions are accelerated to keV energies in the magnetosphere, and may either intercept the magnetopause, and be lost from the planet, or impact the planetary surface. Ions that impact the surface are neutralized, and are then available for resupply to the exosphere. The loss efficiency depends on characteristics of the magnetosphere that determine the fraction of the ions that are recycled by neutralization on the surface. Over the preceding decade, we have collected sodium and potassium data for Mercury at irregular intervals. We analyzed these data to extract values for the Na/K ratio at a variety of conditions on Mercury. Additional information is contained in the original extended abstract.

Potter, A. E.↗

Lick optical spectra of quasar HS 1946+7658 at 10 kilometers per second resolution Lyman-alpha forest and metal absorption systems

We present optical spectra of the most luminous known quasi stellar object (QSO) HS 1946+7658 (z(sub em) = 3.051). Our spectra have both full wavelength coverage, 3240-10570 A, and in selected regions, either high signal-to-noise ratio, SNR approximately equals 40-100, or unusually high approximately 10 km/sec resolution, and in parts of the Ly alpha forest and to the red of Ly alpha emission they are among the best published. We find 113 Ly alpha systems and six metal-line systems, three of which are new. The metal systems at z(sub abs) = 2.844 and 3.050 have complex velocity structure with four and three prominent components, respectively. We find that the system at z(sub abs) = 2.844 is a damped Ly alpha absorption (DLA) system, with a neutral hydrogen column density of log N(H I) = 20.2 +/- 0.4, and it is the cause of the Lyman limit break at lambda approximately equals 3520 A. We believe that most of the H I column density in this system is in z(sub abs) = 2.8443 component which shows the strongest low-ionization absorption lines. The metal abundance in the gas phase of the system is (M/H) approximately equals -2.6 +/- 0.3, with a best estimate of (M/H) = -2.8, with ionizaion parameter log gamma = -2.75, from a photoionization model. The ratios of the logarithmic abundances of C, O, Al, and Si are all within a factor of 2 of solar, which is important for two reasons. First, we believe that the gas abundances which we measure are close to the total abundances, because the ratio of aluminum to other elements is near cosmic, and Al is a refractory element which depletes very readily like chromium, in the interstellar medium. Second, we do not see the enhancement of O with respect to C of (O/C) approximately equals 0.5-0.9 reported in three partial Lyman limit systems by Reimers et al. (1992) and Vogel & Reimers (1993); we measure (O/C) = -0.06 for observed ions and (O/C) approximately equals 0.2 after ionization corrections, which is consistent with solar abundances. We see C II*(lambda 1335) offset by 15 km/sec with respect to C II(lambda 1334), presumably because the gas density varies from 2 to 8 cm(exp -3) with changing velocity in the DLA system. These densities imply that the damped component is 6-25 pc thick, which is reasonable for a single cloud in a cold spiral disk. They also imply that the cloud is relatively highly ionized with more C III than C II, more O III than O I, and log N(H I) = 20.72, which is 3 times the H I column. The system at z(abs) = 1.7382 is also believed to be damped with N(H I) approximately equals 10(exp 21) cm(exp -2), because we see Cr II, but its Ly alpha line will never be seen because it is below the Lyman limit of the other DLA system. We see a 2.6 sigma lack of Ly alpha forest lines well away from the QSO redshift, which may be a chance fluctuation. We also see a correlation between column density N(H I) and Doppler parameter b for 96 unsaturated Ly Alpha forest absorption lines, and although this correlation persists in the 36 Ly alpha lines which lie in regions where the SNR approximately equals 8-16, we agree with Rauch et al. (1993) that it is probably a bogus effect of low supernova remnant (SNR). The same applies to lines with very low b values: in regions where SNR less than or equal to 8 we see many Ly alpha lines which appear to have 10 less than or equal to b less than or equal to 20, but when 8 less than or equal to SNR less than or equal to 16 we see only one line with b less than or equal to 15 km/sec, and two others which we believe have b less than or equal to 20, with values of 20 and 16 km/sec. Traditional Ly alpha line samples which include all lines which have W/sigma(W) greater than or equal to 4 are not adequate to explore the distribution of the properties of individual clouds, because we need much higher (W/sigma(W)) and SNR to avoid the strong biases.

Fan, Xiao-Ming↗

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" 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 "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗

Delamination growth analysis in quasi-isotropic laminates under loads simulating low-velocity impact

A geometrically nonlinear finite-element analysis was developed to calculate the strain energy released by delamination plates during impact loading. Only the first mode of deformation, which is equivalent to static deflection, was treated. Both the impact loading and delamination in the plate were assumed to be axisymmetric. The strain energy release rate in peeling, G sub I, and shear sliding, G sub II, modes were calculated using the fracture mechanics crack closure technique. Energy release rates for various delamination sizes and locations and for various plate configurations and materials were compared. The analysis indicated that shear sliding (G sub II) was the primary mode of delamination growth. The analysis also indicated that the midplane (maximum transverse shear stress plane) delamination was more critical and would grow before any other delamination of the same size near the midplane region. The delamination growth rate was higher (neutrally stable) for a low toughness (brittle) matrix and slower (stable) for high toughness matrix. The energy release rate in the peeling mode, G sub I, for a near-surface delamination can be as high as 0.5G sub II and can contribute significantly to the delamination growth.

Shivakumar, K. N.↗

Reionization of the Universe and the Photoevaporation of Cosmological Minihalos

The first sources of ionizing radiation to condense out of the dark and neutral Intergalactic Medium (IGM) sent ionization fronts sweeping outward through their surroundings, overtaking other condensed objects and photoevaporating them. This feedback effect of universal reionization on cosmic structure formation is demonstrated here for the case of a cosmological minihalo of dark matter and baryons exposed to an external source of ionizing radiation with a quasar-like spectrum, just after the passage of the global ionization front created by the source. We model the pre-ionization minihalo as a truncated, nonsingular isothermal sphere in hydrostatic equilibrium following its collapse out of the expanding background universe and virialization. Results are presented of the first, gas dynamical simulations of this process, including radiative transfer. A sample of observational diagnostics is also presented, including the spatially-varying ionization levels of C, N, and O in the flow if a trace of heavy elements is present and the integrated column densities of H I, He I and He II, and C IV through the photoevaporating gas at different velocities, which would be measured in absorption against a background source like that responsible for the ionization.

Shapiro, Paul R.↗

Impact assessment of geomagnetic induced current neutral blocking devices with power electronics sources on distance relay for 230 kV transmission lines

Electrical grids increased the integration of power electronics sources onto transmission lines. In addition, recent reports show a significant rise in geomagnetic storms occurring in 2024. These geomagnetic storms can affect transmission lines by inducing electrical currents within them, potentially causing power outages due to overloaded power transformers. Geomagnetic induced current neutral blocking devices (GIC-NBDs) are capacitors on the ground of wye power transformers, to avoid damage caused by geomagnetic storms. The integration of distance relays and transmission lines with power electronics sources and GIC-NBDs needs to be studied to observe if GIC-NBDs and power electronics sources could adversely affect operation of the electrical grid. In this analysis, the effect of 2 MW power electronics sources and GIC-NBDs on distance relays is assessed for different 230 kV transmission line lengths in radial and nonradial power systems. The simulations measured the apparent impedances, with different electrical faults. Results for a typical 2,650 μF GIC-NBD application were presented on impedance plots, and the distance relay model operations were assessed. In conclusion, the tests for radial and nonradial power systems showed the behavior of the distance elements and source impedance ratios, assessing the effect of power electronics sources and GIC-NBDs on 230 kV transmission lines.

Distance relays↗

Nonthermal Ionization of Kilonova Ejecta: Observable Impacts

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass (≲0.05 M ⊙ ) and high ejecta velocity (≳0.1c). We show that, even at very early times (≲2 days), regions of ejecta fall below critical density and temperature thresholds at which nonlocal thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized r-process elements. We find that incorporating ionization from high-energy electrons produces an “inverted” and “blended” ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta, and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light-curve properties of AT 2017gfo, such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity (∼0.3 c) ejecta components, our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

High energy astrophysics↗

On the Language of Reliability: A System Engineer Perspective

In its classical definition, risk is defined by three elements: what can go wrong, what are its consequences, and how likely is it to occur? While this definition makes sense in a regulatory-based framework where for the current fleet of operating light water reactors (LWRs), the risks associated with nuclear power plants typically are characterized in terms of core damage and large early release frequency (LERF), this approach does not provide a useful snapshot of the health of the plant from a broader perspective. This is due to the very narrow context in which the term “risk” typically is defined as nuclear safety aspects that have the potential to impact public health. In this paper, we take the viewpoint of nuclear safety that is reflective of the current fleet of operating LWRs for which core damage frequency and LERF are appropriate metrics. For other advanced reactor designs, other more applicable technology neutral metrics of reactor safety metrics would be specified. A possible alternate path would start by redefining the word risk with a broader meaning that better reflects the needs of a system health and asset management decision-making process. Rather than asking how likely an event could occur (in probabilistic terms), we can ask how far this event is from occurring. Our approach starts by defining and quantifying component and system health in terms of a “distance” between its actual and limiting conditions, i.e., determination of the margin that exists between the current state/condition and the state where the component/system is no longer capable of achieving its intended function. A margin is a measure that is more reflective of the current state or performance of a component, and therefore more closely tied to decisions that are made on an ongoing basis. Further, we will show how, given the data available from plant equipment reliability and monitoring (e.g., pump vibration data) and prognostic (e.g., component remaining useful life estimation) data, a margin can be described and determined for all types of maintenance approaches (e.g., corrective or predictive maintenance). We show how classical reliability models (e.g., fault trees) can be used to quantify the system margin provided component margin values. In the approach described in this paper, the propagation of margin values through classical reliability models are not performed using classical probabilistic calculations applied to sets (as performed in a typical plant probabilistic risk assessment). Instead, we show how it is possible to propagate margin values through Boolean logic gates (i.e., AND and OR operators) through distance-based operations.

97 MATHEMATICS AND COMPUTING↗

Radiative recombination of complex ions

The rates of radiative capture by atomic systems is computed for the first four ionization stages of the abundant elements C, N, O, Ne, Mg, Si, S, and Ar. A simple prescription is given for calculating the rates for systems in higher ionization stages. Results for capture to atomic helium are also given. In the general case, recent calculations of photoionization cross sections have been used to compute capture rates to levels of the ground-state configuration of the recombined species. Captures to other levels of the valence shell and to the higher shells are evaluated using hydrogenic formulae but with a correction factor which takes into account the nonhydrogenic nature of these states. This correction factor is due to the incomplete shielding by the inner electrons and represents essentially an effective charge for the recombining ion; it is determined semiempirically from an appropriate weighting of values derived from the observed level structure for each species. In some cases, for the capture to neutral atomic systems, its application increases the resulting computed recombination rate by a significant amount (about 50%). Computed recombination rates are tabulated for the four ionization stages at an electron temperature of 10,000 K and for the first stage at 100 K. A convenient procedure is outlined for evaluating the rates at other temperatures in the general neighborhood of these values.

Gould, R. J.↗

Unraveling the solvation phenomena of Ln 3+ cations in room-temperature ionic liquids: A computational study

Rare-earth elements (REEs), classified as critical materials, are difficult to separate due to their similar chemical properties and slight differences in ionic radius. Room-temperature ionic liquids (RTILs) have gained significant attention for REE separation because of their unique physicochemical properties and environmental advantages. In this study, we have investigated the solvation mechanisms of two RE cations (Ln 3+ ), Nd 3+ (light REE) and Yb 3+ (heavy REE), in two 1-butyl-3-methylimidazolium ([BMIM] + )-based RTILs using bis(trifluoromethylsulfonyl)imide ([NTf 2 ] − ) and acetate ([OAc] − ) as the respective anions, employing classical molecular dynamics (MD) simulations and density functional theory (DFT) techniques. From MD simulations, it is evident that Ln 3+ cations are primarily solvated by RTIL anions in the first solvation shell, while [BMIM] + cations form the second solvation shell with their numbers depending on the size and composition of the first solvation shell. The neutralizing [NO 3 − ] counterions are mainly solvated by [BMIM] + cations in their first solvation shell. Relative free energy change of solvation calculations using thermodynamic integration (TI) method indicate that Yb 3+ is more strongly solvated than Nd 3+ in both RTILs, a trend further supported by DFT calculations. While both methods predict consistent qualitative behavior, differences in models and energy evaluations lead to variations in absolute values. Overall, this study provides a comprehensive understanding of the solvation behavior of Ln 3+ cations in RTILs, demonstrating a stronger solvation preference for the heavy Ln 3+ cation (Yb 3+ ). In conclusion, these findings have implications for the design of RTIL-based separation processes for REEs.

Ash, Tamalika [Ames Lab., and Iowa State Univ., Am↗

Flowability of Crumbler Rotary Shear Size-Reduced Granular Biomass: An Experiment-Informed Modeling Study on the Angle of Repose

Biomass has potential as a carbon-neutral alternative to petroleum for chemical and energy products. However, complete replacement of fossil fuel is contingent upon efficient processes to eliminate undesirable characteristics of biomass, e.g., low bulk density, variability, and storage-induced quality problems. Mechanical size reduction via comminution is a processing operation to engineer favorable biomass flowability in handling. Crumbler rotary shear mill has been empirically demonstrated to produce more uniformly shaped particles with higher flowability than hammermilled biomass. This study combines modeling and experimentation to unveil fundamental understandings of the relation between granular particle characteristics and biomass flow behavior, which elucidate underlying mechanisms and guide selection of critical processing parameters. For this purpose, the impact of critical material attributes, including particle size (2–6 mm), particle shape (briquette, chip, clumped-sphere, cube, etc.), and surface roughness, on the angle of repose (AOR) of milled pine chips were investigated using discrete element method (DEM) simulations. Forest Concepts Crumbler rotary shear system is used to produce milled pine particles within the same size range considered in DEM simulations. AOR of different sets of these particles were measured experimentally to benchmark DEM results against experimental data. Specific energy consumption for the comminution of biomass with different particle size and moisture content are measured for technoeconomic analysis. Our results show that the smaller size (2 mm) of pine particle achieves better followability (i.e., smaller AOR) while the energy cost of comminution is significantly higher and bulk density is almost the same as the 6-mm pine particles. For the 2-mm particle size, Crumbles from veneer have better flow properties than Crumbles from chips. Contrarily, no significant difference was observed between the AOR of the two materials for the 6-mm particle size. Furthermore, from DEM simulations, mechanical interlocking between particles was found as a dominant factor in determining AOR of complex-shaped particles such as milled pine, which cannot be accurately captured by using simple particle shapes (e.g., mono-sphere) with a rolling resistance model. Conversely, clumped-sphere model alleviates this limitation without increasing computational cost significantly and can be used for accurate representation of biomass granular particles when simulating free-flow behavior.

09 BIOMASS FUELS↗