Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Graph properties”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

RivGraph: Automatic extraction and analysis of river and delta channel network topology

River networks sustain life and landscapes by carrying and distributing water, sediment, and nutrients throughout ecosystems and communities. At the largest scale, river networks drain continents through tree-like tributary networks. At typically smaller scales, river deltas and braided rivers form loopy, complex distributary river networks via avulsions and bifurcations.In order to model flows through these networks or analyze network structure, the topology, or connectivity, of the network must be resolved. Additionally, morphologic properties of each river channel as well as the direction of flow through the channel inform how fluxes travel through the network’s channels. Riv Graphis a Python package that automates the extraction and characterization of river channel networks from a user-provided binary image, or mask, of a channel network (Fig. 1). Masks may be derived from (typically remotely-sensed) imagery, simulations, or even hand-drawn. RivGraph will create explicit representations of the channel network by resolving river centerlines as links, and junctions as nodes. Flow directions are solved for each link of the network without using auxiliary data, e.g., a digital elevation model (DEM). Morphologic properties are computed as well, including link lengths, widths, sinuosities, branching angles,and braiding indices. If provided,RivGraph will preserve georeferencing information of the mask and will export results as ESRI shapefiles, GeoJSONs, and GeoTIFFs for easy import into GIS software.RivGraph can also return extracted networks as networkx objects for convenient interfacing with the full-featured networkx package (Hagberg et al., 2008). Finally, RivGraph offers a suite of topologic metrics that were specifically designed for river channel network analysis (Tejedor et al., 2015b).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

LatticeAnalytics: Strut-Level Visualization and Inspection of Additively Manufactured Lattice Structures

Additive manufacturing (AM) is revolutionizing the production of custom components with complex internal geometries, essential for high-performance applications in diverse fields such as medicine and defense. These AM parts optimize strength while minimizing weight by utilizing internal lattice structures consisting of large quantities of small interconnected struts. However, the complexity of these structures, combined with the challenges of using X-ray Computed Tomography (XCT) data, makes validation of part reliability difficult. This ultimately inhibits the development of novel parts for our collaborating material scientists. Here, we introduce LatticeAnalytics, a novel framework specifically designed for visual inspection of defects in these lattice structures. Our framework offers an end-to-end solution that includes the data management of XCT scans, enables remote access for geographically dispersed teams through a web-based dashboard, and incorporates novel visualizations. Our analysis is facilitated by a coarse alignment between the lattice’s nominal model, a spatial graph, and the XCT data. We employ a simple VR-based approach for fast and rough alignment, followed by an offline registration and identification of the struts. With the nodes and struts aligned and identified in the volume, our framework allows querying of subvolumes containing a single strut at multiple resolutions. This avoids computation over the entire lattice and also allow for easy parallelization of down-stream computations, such as strut-specific metrics. To depict a fast overview of the strut quality, we introduce two innovative visual encodings, crucial for our collaborators’ research in creating novel AM parts: the Contour View and the Roughness Map, which depict critical geometrical and surface features of individual struts in standardized two 2D views. We evaluated the integrated system through expert interviews. The feedback confirms the framework’s practicality and its effectiveness in enhancing current inspection workflows. It solves major bottlenecks for our collaborators, ultimately helping them create novel parts with advanced properties.

Miao, Haichao [Lawrence Livermore National Laborat↗
DOE OSTI2026

Machine Learning Vacancy Formation Energy in Nickel-Based Superalloys

Thermal vacancies play a critical role in high-temperature Ni-based superalloys and influence various properties such as creep resistance, oxidation, etc. This study systematically investigates the impact of commonly used transition metals (Cr, Co, Fe), refractory metals (Nb, Ta, Mo, W) and other elements (Al, Cu, Ti, Mn) on the thermodynamic stability of 36 binary, 20 ternary, 11 quaternary, 9 quinary, and 3 senary FCC Ni-based alloys covering various elemental combinations. Density functional theory-based studies on Ni-X binary alloys show that higher concentrations of Cr, Nb, Ta, Al, and Ti introduce significant lattice distortions and broaden the distribution of vacancy formation energies (standard deviation up to 0.15 eV). These elements partially donate electrons, reducing their self-consistent chemical potentials relative to single-element reference values and lowering vacancy formation energies, while Co, Fe, Mo, and W show lower charge localization. These trends extend from 3-6 element alloys, where Cr, Nb, and Ta-rich compositions have low-energy states (~0.5 eV) that increase vacancy concentrations. Finally, graph neural network models are developed to screen over 5000 virtual alloys. Eleven leading compositions are identified with mean vacancy formation energy higher than 1.75 eV and vacancy concentration ~2 orders of magnitude lower than pure Ni at 1000 K. These results provide valuable guidelines to achieve controlled defect engineering in structural alloys.

DFT↗
DOE OSTI2017

Occupation Time of a Randomly Accelerated Particle on the Positive Half Axis: Results for the First Five Moments

In the random acceleration process a point particle is accelerated by Gaussian white noise with zero mean. Although several fundamental statistical properties of the motion have been analyzed in detail, the statistics of occupation times is still not well understood. We consider the occupation or residence time T{sub +} on the positive x axis of a particle which is randomly accelerated on the unbounded x axis for a time t. The first two moments of T{sub +} were recently derived by Ouandji Boutcheng et al. (J Stat Mech 053213:1–10, 2016). With an alternate approach utilizing basis functions which have proved useful in other studies of randomly accelerated motion, results for the first five moments are obtained in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Utah FORGE: Stress Logging Data

This spreadsheet consist of data and graphs from deep well 58-32 stress testing from 6900 - 7500 ft depth. Measured stress data were used to correct logging predictions of in situ stress. Stress plots shows pore pressure (measured during the injection testing), the total vertical in situ stress (determined from the density logging) and the total maximum and minimum horizontal stresses. The horizontal stresses were determined from the DSI (Dipole Sonic Imager) and corrected to match the direct measurements.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Determining Compression Characteristics of Honeycomb Material - 19674

The objective is to determine compression test characteristics of stainless steel honeycomb material to be able to represent honeycomb structures accurately in analytical models used to simulate hypothetical accident scenarios of shipping packages. Honeycomb is a material used primary in the aerospace industry due to its high strength-to-weight ratio. Because it doesn't have a shelf life and can withstand high heats, it is an excellent candidate for a structural material in package designs. Honeycomb has orthotropic material properties. The honeycomb currently being evaluated is made of metal ribbons spot welded together to form hexagon pattern between metal plates. The hexagon ribbons are brazed to the metal plates. The complexity of the honeycomb significantly lengthens the simulation time, which is further complicated when the brazing, welding, and imperfections of the material is considered. This means to affectively represent the honeycomb in simulation programs, such as Abaqus, the material needs to be approximated as a uniform orthotropic material. This requires material properties in each of the three directions, T, W, and L. The T direction is defined as the 'strong' direction, perpendicular to the honeycomb sheet. The L direction is parallel the ribbon and the W direction is perpendicular to the ribbon. Honeycomb has 3 different compression stages. Stage 1 is the initial compression where the honeycomb maintains its structural integrity and does not permanently deform from loads from normal operations. Stage 2 is when the initial buckling causes irreversible damage to the honeycomb. Stage 2 is the one we are most interested in because it absorbs the most energy from a hypothetical accident scenario. During Stage 2, the honeycomb fails layer-by-layer, indicating that the more layers, the longer material crushing is sustained, and thus the more energy that is absorbed. Stage 3 is final compression, similar to compressing solid metal. During stage 3 the stress increases with a diminishing rate of elongation and the honeycomb is completely failed where the plates between the honeycomb core sandwich the crushed honeycomb ribbon. The stress-strain graph bellow shows compression test of 3 different samples in the T direction. The red lines divide the three stages. The far left is Stage 1, the middle stage 2, and the right stage 3. The data collected is displayed in the table below. The expressions in the table represent the regression formula that represents each stage of the compression on a stress-strain graph. Stage 1 is assumed to intersect with the origin; data for stage 3 for the crush in W and L directions where unable to be gathered due to the nature of failure for those directions. Stage 1 and 2 are linear regressions while stage 3 is a degree 2, to best match the curve. Linear regression is used on stage 2 because ultimately that represents the energy absorbed. With more data, a degree n x 2 regression would be a more appropriate regression for stage 2, where n is the number of layers. The table below is the approximation of each scenario and stage. Each stage starts and ends at the intersection of the next stage. There where two main objectives for these tests: firstly to understand how honeycomb performs under extreme compression, and secondly to be able to numerically represent the honeycomb structure. Both of these where completed to varying degrees. It is important understood how the honeycomb fails. If it is crushed in the T direction it retains its integrity even after being crushed, however when crushed in the L and W crush direction, if it fails, it disintegrates, and loses all integrity. Also, the more layers the more time the material spends in stage 2. Failure is started by buckling, thus if there is any imperfection, the stress will not spike but transition straight into stage 2. The data collected and aggregated can be used for initial simulation of honeycomb used in packages. The initial testing has set the ground work for more data to be collected in order to verify results and to allow more confidence in the simulation results. This is only the initial data collected. The next steps is to continue to collect more data to verify results and to test more variants of honeycomb, with different brazing, layers, and shape.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Rare Higgs Processes at CMS and Precision Timing Detector Studies for HL-LHC CMS Upgrade

This thesis describes the search for two rare Higgs processes. The first analysis describes the CMS Run 2 search for $H$ $\rightarrow$ $\mu$$\mu$ decays, with 137.3 fb$^{-1}$ of data at $\sqrt{s}$ = 13 TeV. The analysis targeted four different Higgs production modes: the gluon fusion (ggH), the vector boson fusion (VBF), the Higgs-strahlung process (VH), and the production in association with a pair of top quarks (ttH). Each category used a dedicated machine learning based classifier to separate the signal from the background processes. A combined fit from all these categories saw a slight excess in the data corresponding to 3.0 standard deviations at $M$$_{H}$ = 125.38 GeV, and gave the first evidence for the Higgs boson decay to second-generation fermions. The best-fit signal strength and the corresponding 68% CL interval was found to be +0.17?????? = 1.19 $_{-0.39}^{+0.41}$ (stat)$_{-0.16}^{+0.17}$(syst) at $M$$_{H}$ = 125.38 GeV. The second analysis describes the CMS Run 2 search for 𝐻𝐻 → 𝑏𝑏𝑏𝑏 with highly boosted Higgs bosons. This analysis used a dedicated jet identification algorithm based on graph neural networks (ParticleNet) to identify boosted H→ bb jets. This search targeted the gluon fusion and the vector boson fusion HH production modes, and put constraints on the allowed values of the various Higgs couplings as: 𝜅𝜆 ∈ [−9.9, 16.9] when 𝜅𝑉 = 1, 𝜅2𝑉 = 1; 𝜅𝑉 ∈ [−1.17, −0.79] ∪ [0.81, 1.18] when 𝜅𝜆 = 1, 𝜅2𝑉 = 1; 𝜅2𝑉 ∈ [0.62, 1.41] when 𝜅𝜆 = 1, 𝜅𝑉 = 1. A scenario with 𝜅2𝑉 = 0 was excluded with a significance of 6.3 standard deviations for the first time, when other H couplings are fixed to their SM values. The combined observed (expected) 95% upper limit on the HH production cross section was found to be 9.9 (5.1) × SM. Finally, this thesis also discusses the planned MIP Timing Detector (MTD) upgrade for CMS at the HL-LHC. The MTD will be a time-of-flight (TOF) detector, designed to provide a precision timing information for charged particles using SiPMs + LYSO scintillating crystals, with a time resolution of ∼30 ps. This thesis describes several R&D tests that have been performed for characterizing the sensor properties (time resolution, light yield, etc.) and optimizing the sensor design geometry. This thesis also contains a description of mock test setups for cooling the sensors, since it is known to be an effective way of mitigating the increased dark current rates in the sensors due to radiation damage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Provenance infrastructure for multimodal x-ray experiments and reproducible analysis

Exploring relationships between material structures and desired properties for material discovery and synthesis requires mining machine-readable databases and linking experimental data with sufficiently rich annotations. There are vast gaps both in these data sources and in the connections between them. We present a provenance-based data management and analysis framework that enables capturing, persisting, and reanalyzing experimental data from the BNL National Synchrotron Light Source II (NSLS-II). Our system leverages NSLS-II Bluesky for experimental data acquisition, captures analysis parameters, and enriches the search space with results from the scientific literature. We describe in detail the challenges and benefits of data-driven approaches when comparing computational and experimental data, such as the multiple possible scenarios of interest in often ill-posed inverse problems and missing data. Our infrastructure utilizes the NSLS-II event model and includes discovery mechanisms across multiple beamlines. A text and data mining module using natural language processing classification techniques extracts relevant information from the literature that is added to data sources. We describe a provenance acquisition and analysis rerun workflow based on a directed acyclic graph that keeps track of process and data provenance and allows to replay analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

GraphCH: A Deep Framework for Assessing Cyber-Human Aspects in Insider Threat Detection

Insider threat is one of the most damaging cyber attacks that could cause the loss of intellectual property and enterprise data security breaches. Action sequence data such as host logs are used to investigate such threats and develop anomaly-based AI detectors. However, insider threat actions are similar to legitimate user activities, causing AI detectors to fail and suffer from high false alarm rates. Therefore, user cyber activity logs are inadequate to fully unfold insider threats. In this study, we adopt human psychological principles of risk-taking and impulsiveness along with host data to assess the influence and usefulness of human behavioral aspects in insider threat detection. Here, we hypothesize that individuals' impulsive and risk-taking behavior correlates with cyberspace activities. To validate our hypothesis, we conducted an IRB-approved study recruiting 35 participants who work in a large U.S. university and collected their cyber and psychological data for 90 days. Host and human-behavioral data analysis and mapping indicate that impulsive and risk-taking users trigger more system errors causing (un)intentional insider threats and are susceptible to attackers' social engineering and cognitive hacking. Utilizing cyber-human aspects, we introduce a Cyber-Human Graph Neural Network (GNN) based framework GraphCH to identify abnormal user behaviors and detect insider threats.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Self-stabilizing Connected Components

For the problem of computing the connected components of a graph, this paper considers the design of algorithms that are resilient to transient hardware faults, like bit Hips. More specifically, it applies the technique of self-stabilization. A system is self-stabilizing if, when starting from a valid or invalid state, it is guaranteed to reach a valid state after a finite number of steps. Therefore on a machine subject to a transient fault, a self-stabilizing algorithm could recover if that fault caused the system to enter an invalid state. We give a comprehensive analysis of the valid and invalid states during label propagation and derive algorithms to verify and correct the invalid state. The self-stabilizing label-propagation algorithm performs O (V log V) additional computation and requires O (V) additional storage over its conventional counterpart (and, as such, does not increase asymptotic complexity over conventional label propagation). When run against a battery of simulated fault injection tests, the self-stabilizing label propagation algorithm exhibits more resilient behavior than a triple modular redundancy (TMR) based fault-tolerant algorithm in 80% of cases. From a performance perspective, it also outperforms TMR as it requires fewer iterations in total. Beyond the fault tolerance properties of self-stabilizing label-propagation, we believe, they are useful from the theoretical perspective; and may have other use-cases.

Sao, Piyush↗
DOE OSTI2023

Interaction energy and isosteric heat of adsorption between hydrogen and magnesium diboride

Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).

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

Identifying the Unique Properties of α-Bi 2 Mo 3 O 12 for the Activation of Propene

In order to understand the remarkable activity of α-Bi 2 Mo 3 O 12 for selective oxidation and ammoxidation of propene, the propene activation ability of four molybdenum-based mixed metal oxides - Bi 2 Mo 3 O 12 , PbMoO 4 , Bi 2 Pb 5 Mo 8 O 32 , and MoO 3 - was investigated using density functional theory. Propene activation is considered to occur via abstraction of a hydrogen atom from the methyl group of physisorbed propene by lattice oxygen. For each material, the apparent activation energy was estimated by summing the heat of adsorption of propene, the C-H bond dissociation energy, and the hydrogen attachment energy (HAE) for hydrogen addition to lattice oxygen; this sum provides a lower bound for the apparent activation energy. It was found that two structural features of oxide surfaces are essential to achieve low activation barriers: under-coordinated surface cation sites enable strong propene adsorption, and suitable 5- or 6-coordinate geometries at molybdenum result in favorable HAEs. The impact of molybdenum coordination on HAE was elucidated by carrying out a molecular orbital analysis using a cluster model of the molybdate unit. This effort revealed that, in 5- and 6-coordinate molybdates, oxygen donor atoms trans to molybdenyl oxo atoms destabilize the molybdate prior to H addition but stabilize the molybdate after H addition, thereby providing an HAE ~15 kcal/mol more favorable than that on 4-coordinate molybdate oxo atoms. Bi 3+ cations in Bi 2 Mo 3 O 12 thus promote catalytic activity by providing both strong adsorption sites for propene and forcing molybdate into 5-coordinate geometries that lead to particularly favorable values of the HAE. (Graph Presented).

36 MATERIALS SCIENCE↗
NASA NTRS1999

Modeling of THz Lasers Based on Intersubband Transitions in Semiconductor Quantum Wells

In semiconductor quantum well structures, the intersubband energy separation can be adjusted to the terahertz (THz) frequency range by changing the well width and material combinations. The electronic and optical properties of these nanostructures can also be controlled by an applied dc electric field. These unique features lead to a large frequency tunability of the quantum well devices. In the on-going project of modeling of the THz lasers, we investigate the possibility of using optical pumping to generate THz radiation based on intersubband transitions in semiconductor quantum wells. We choose the optical pumping because in the electric current injection it is difficult to realize population inversion in the THz frequency range due to the small intersubband separation (4-40 meV). We considered both small conduction band offset (GaAs/AlGaAs) and large band offset (InGaAs/AlAsSb) quantum well structures. For GaAs/AlGaAs quantum wells, mid-infrared C02 lasers are used as pumping sources. For InGaAs/AlAsSb quantum wells, the resonant intersubband transitions can be excited by the near-infrared diode lasers. For three- and four-subband quantum wells, we solve the pumpfield-induced nonequilibrium distribution function for each subband of the quantum well system from a set of rate equations that include both intrasubband and intersubband relaxation processes. Taking into account the coherent interactions between pump and THz (signal) waves, we calculate the optical gain for the THz field. The gain arising from population inversion and stimulated Raman processes is calculated in a unified manner. A graph shows the calculated THz gain spectra for three-subband GaAs/AlGaAs quantum wells. We see that the coherent pump and signal wave interactions contribute significantly to the gain. The pump intensity dependence of the THz gain is also studied. The calculated results are shown. Because of the optical Stark effect and pump-induced population redistribution, the maximum THz gain saturates at larger pump intensities.

Liu, Ansheng↗
NASA NTRS2004

Particulate Titanium Matrix Composites Tested--Show Promise for Space Propulsion Applications

Uniformly distributed particle-strengthened titanium matrix composites (TMCs) can be manufactured at lower cost than many types of continuous-fiber composites. The innovative manufacturing technology combines cold and hot isostatic pressing procedures to produce near-final-shape components. Material stiffness is increased up to 26-percent greater than that of components made with conventional titanium materials at no significant increase in the weight. The improved mechanical performance and low-cost manufacturing capability motivated an independent review to assess the improved properties of ceramic titanium carbide (TiC) particulate-reinforced titanium at elevated temperature. Researchers at the NASA Glenn Research Center creatively designed and executed deformation and durability tests to reveal operating regimes where these materials could lower the cost and weight of space propulsion systems. The program compares the elevated-temperature performance of titanium alloy Ti-6Al-4V matrix material to an alloy containing 10 wt% of TiC particles. Initial experiments showed that at these relatively low particle concentrations the material stiffness of the TMC was improved 20 percent over that of the plain Ti-6Al-4V alloy when tested at 427 C. The proportional limit and ultimate strength of the composite in tension are 21- and 14-percent greater than those of the plain alloy. Compression tests showed that the proportional limit is about 30 percent greater for TMC than for the plain alloy. The enhanced deformation resistance of the TMC was also evident in a series of tensile and compressive stress relaxation tests that were made. Specimens were subjected to tensile or compressive strain amplitudes of 0.75 percent for 24 hr followed by a return to zero strain imposed for 24 hr. The stress relaxation data were normalized with respect to the maximum stress for each case and plotted as a function of time in the following graph. Tensile stresses relaxed 19 percent for the TMC and 25 percent for the plain Ti-6Al-4V alloy. Compressive stresses relaxed 25 percent for the TMC and 39 percent for the plain Ti-6Al-4V alloy. The superior deformation resistance of the TMC extends to a creep rate that is 28-percent slower for the TMC when it is loaded to stress levels that are 26-percent higher than for the plain Ti-6Al-4V alloy.

Lerch, Bradley A.↗
DOE OSTI2019

Materials design of perovskite solid solutions for thermochemical applications

Perovskites are excellent candidate materials as oxygen carriers in thermochemical processes. Due to their versatile composition, it is possible to fine-tune the perovskite's properties. Here, we present a method for the rational design of AMO 3– δ perovskite solid solutions with two different species on the Msite in order to tune their redox behavior. To account for the different ionic radii of different transition metal species M, two distinct ions are used in a solid solution on the A site, allowing tolerance-factor adjusted materials design. Using this methodology, we can create stable perovskites over a large range of different compositions. Leveraging the infrastructure of Materials Project, we calculate redox enthalpies for the reduction of over 240 of these $(\text {A}^{\prime}_x\text {A}^{\prime\prime}_{1-x})(\text {M}^{\prime}_y\text {M}^{\prime\prime}_{1-y})$O 3 perovskites to their brownmillerite phases based on density functional theory (DFT). We compare this data to experimentally measured data on thermodynamics of 24 of these materials to verify our theoretical framework. An empirical model is formulated for predicting the enthalpy and entropy changes as a function of the perovskites non-stoichiometry δ , which can be used to simulate the equilibrium composition as a function of temperature and oxygen partial pressure and to create a perovskite search engine based on an energetic analysis of the redox cycles. The data has been added as a contribution to MPContribs, which now includes publicly available user-controlled interactive graphs based on our theoretical and experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Can a deep-learning model make fast predictions of vacancy formation in diverse materials?

The presence of point defects, such as vacancies, plays an important role in materials design. Here, we explore the extrapolative power of a graph neural network (GNN) to predict vacancy formation energies. We show that a model trained only on perfect materials can also be used to predict vacancy formation energies (E vac ) of defect structures without the need for additional training data. Such GNN-based predictions are considerably faster than density functional theory (DFT) calculations and show potential as a quick pre-screening tool for defect systems. To test this strategy, we developed a DFT dataset of 530 E vac consisting of 3D elemental solids, alloys, oxides, semiconductors, and 2D monolayer materials. We analyzed and discussed the applicability of such direct and fast predictions. We applied the model to predict 192 494 E vac for 55 723 materials in the JARVIS-DFT database. Our work demonstrates how a GNN-model performs on unseen data.

2D materials↗
DOE OSTI2018

Representing local atomic environment using descriptors based on local correlations

Statistical learning of material properties is an emerging topic of research and has been tremendously successful in areas such as representing complex energy landscapes as well as in technologically relevant areas, like identification of better catalysts and electronic materials. However, analysis of large data sets to efficiently learn characteristic features of a complex energy landscape, for example, depends on the ability of descriptors to effectively screen different local atomic environments. Thus, discovering appropriate descriptors of bulk or defect properties and the functional dependence of such properties on these descriptors remains a difficult and tedious process. To this end, we develop here a framework to generate descriptors based on many-body correlations that can effectively capture intrinsic geometric features of the local environment of an atom. These descriptors are based on the spectrum of two-body, three-body, four-body, and higher order correlations between an atom and its neighbors and are evaluated by calculating the corresponding two-body, three-body, and four-body overlap integrals. They are invariant to global translation, global rotation, reflection, and permutations of atomic indices. By systematically testing the ability to capture the local atomic environment, it is shown that the local correlation descriptors are able to successfully reconstruct structures containing 10-25 atoms which was previously not possible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

Antiproton Trapping for Advanced Space Propulsion Applications

The Summary of Research parallels the Statement of Work (Appendix I) submitted with the proposal, and funded effective Feb. 1, 1997 for one year. A proposal was submitted to CERN in October, 1996 to carry out an experiment on the synthesis and study of fundamental properties of atomic antihydrogen. Since confined atomic antihydrogen is potentially the most powerful and elegant source of propulsion energy known, its confinement and properties are of great interest to the space propulsion community. Appendix II includes an article published in the technical magazine Compressed Air, June 1997, which describes CERN antiproton facilities, and ATHENA. During the period of this grant, Prof. Michael Holzscheiter served as spokesman for ATHENA and, in collaboration with Prof. Gerald Smith, worked on the development of the antiproton confinement trap, which is an important part of the ATHENA experiment. Appendix III includes a progress report submitted to CERN on March 12, 1997 concerning development of the ATHENA detector. Section 4.1 reviews technical responsibilities within the ATHENA collaboration, including the Antiproton System, headed by Prof. Holzscheiter. The collaboration was advised (see Appendix IV) on June 13, 1997 that the CERN Research Board had approved ATHENA for operation at the new Antiproton Decelerator (AD), presently under construction. First antiproton beams are expected to be delivered to experiments in about one year. Progress toward assembly of the ATHENA detector and initial testing expected in 1999 has been excellent. Appendix V includes a copy of the minutes of the most recently documented collaboration meeting held at CERN of October 24, 1997, which provides more information on development of systems, including the antiproton trapping apparatus. On February 10, 1998 Prof. Smith gave a 3 hour lecture on the Physics of Antimatter, as part of the Physics for the Third Millennium Lecture Series held at MSFC. Included in Appendix VI are notes and graphs presented on the ATHENA experiment. Portable antiproton trap has been under development. The goal is to store and transport antiprotons from a production site, such as Fermilab near Chicago, to a distant site, such as Huntsville, AL, thus demonstrating the portability of antiprotons.

Smith, Gerald A.↗