Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Link Analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Combining Agent Based Modeling and System Dynamics to Investigate the Circularity of Plastics

The United States currently produces about 1 million metric ton of ocean plastic pollution annually. One proposed solution to combat plastic waste is a circular economy (CE), which aims to transition from today's take-make-waste linear pattern of production and consumption to a system where the value of resources is maximized over time. Two key methods in industrial ecology are useful in assessing the viability of CE: (1) System Dynamics (SD) and (2) Agent Based Modeling (ABM). In prior work, the plastic life cycle was modeled with SD and ABM. The two models calculate recycling rates and costs in different ways, making it difficult to pinpoint necessary next steps. We integrate the ABM and SD models - linking the emergent patterns from micro-level human decisions to system level processes - which allows a more comprehensive understanding of feedbacks, costs, and environmental impacts. The integrated model is more accurate, and can be used to visualize recycling rates and human health and environmental impacts over time. The difference between the integrated and original SD model prompts a Sobol sensitivity analysis, which is used to understand which behavioral factors most affect plastic recycling patterns. We find that the habitual component is typically the most influential in promoting positive recycling behavior. Additionally, we utilize the combined model to understand and visualize how various behavioral intervention scenarios, like improved access to recycling programs and cart tagging, influence recycling patterns; these results can guide future policy-making.

agent-based modeling↗
DOE OSTI2021

Search for efficient formulations for Hamiltonian simulation of non-Abelian lattice gauge theories

Hamiltonian formulation of lattice gauge theories (LGTs) is the most natural framework for the purpose of quantum simulation, an area of research that is growing with advances in quantum-computing algorithms and hardware. It, therefore, remains an important task to identify the most accurate, while computationally economic, Hamiltonian formulation(s) in such theories, considering the necessary truncation imposed on the Hilbert space of gauge bosons with any finite computing resources. This paper is a first step toward addressing this question in the case of non-Abelian LGTs, which further require the imposition of non-Abelian Gauss’s laws on the Hilbert space, introducing additional computational complexity. Focusing on the case of SU(2) LGT in 1+1 dimensions coupled to one flavor of fermionic matter, a number of different formulations of the original Kogut-Susskind framework are analyzed with regard to the dependence of the dimension of the physical Hilbert space on boundary conditions, system’s size, and the cutoff on the excitations of gauge bosons. The impact of such dependencies on the accuracy of the spectrum and dynamics obtained from a Hamiltonian computation is examined, and the (classical) computational-resource requirements given these considerations are studied. Besides the well-known angular-momentum formulation of the theory, the cases of purely fermionic and purely bosonic formulations (with open boundary conditions), and the Loop-String-Hadron formulation are analyzed, along with a brief discussion of a Quantum-Link-Model formulation of the same theory. Clear advantages are found in working with the Loop-String-Hadron framework which implements non-Abelian Gauss’s laws a priori using a complete set of gauge-invariant operators. Although small lattices are studied in the numerical analysis of this work, and only the simplest algorithms are considered, a range of conclusions will be applicable to larger systems and potentially to higher dimensions. Future studies will extend this investigation to the analysis of resource requirements for quantum simulations of non-Abelian LGT, with the goal of shedding light on the most efficient Hamiltonian formulation of gauge theories of relevance in nature

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Analysis of Tasks in Autonomous Systems Using the EMRALD Dynamic Risk Assessment Tool

An autonomous system refers to the system that has the power and ability for self-governance in the performance of system functions. Autonomous systems have been actively pursued in a variety of domains such as automotive, aviation, maritime, medicine, and nuclear fields. As an unmanned concept employing the highest automation level, the autonomous system basically performs most of the work in normal operations or emergency situations. However, despite advances in technology, many researchers have noted these systems still require human actions. The nature of human actions on autonomous systems is different than the human actions that are considered in existing systems. Nevertheless, only a few studies have been conducted on 1) characterizing the different types of errors and risks associated with human actions interacting with autonomous systems and 2) how to evaluate human actions in the autonomous operations. As a starting point, this study aims to investigate differences of tasks in autonomous operation compared to those in existing nuclear power plant operation using the Event Modeling Risk Assessment Using Linked Diagram (EMRALD) software. In this paper, insights aspect of human error and time are derived out and discussed based on the output of the EMRALD models.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

A Case for Electron-Astrophysics

Abstract The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the research theme of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. In addition, plasma electrons often play an important role for the spatial transfer of thermal energy due to the high heat flux associated with their velocity distribution. The regulation of this electron heat flux is likewise not understood. By focussing on these and other fundamental electron processes, the research theme of electron-astrophysics links outstanding science questions of great importance to the fields of space physics, astrophysics, and laboratory plasma physics. In this White Paper, submitted to ESA in response to the Voyage 2050 call, we review a selection of these outstanding questions, discuss their importance, and present a roadmap for answering them through novel space-mission concepts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Microscale Mechanical‐Chemical Modeling of Granular Salt: Insights for Creep

Abstract Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical‐chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer‐term creep of salt at larger scales.

58 GEOSCIENCES↗
DOE OSTI2023

Hybrid Simulation and Quasi-linear Theory of Bi-Kappa Proton Instabilities

The quasi-steady states of collisionless plasmas in space (e.g., in the solar wind and planetary environments) are governed by the interactions of charged particles with wave fluctuations. These interactions are responsible not only for the dissipation of plasma waves but also for their excitation. The present analysis focuses on two instabilities, mirror and electromagnetic ion cyclotron instabilities, associated with the same proton temperature anisotropy T ⊥ > T ∥ (where ⊥, ∥ are directions defined with respect to the local magnetic field vector). Theories relying on standard Maxwellian models fail to link these two instabilities (i.e., predicted thresholds) to the proton quasi-stable anisotropies measured in situ in a completely satisfactory manner. Here we revisit these instabilities by modeling protons with the generalized bi-Kappa (bi-κ power-law) distribution, and by a comparative analysis of a 2D hybrid simulation with the velocity-moment-based quasi-linear (QL) theory. It is shown that the two methods feature qualitative and, even to some extent, quantitative agreement. The reduced QL analysis based upon the assumption of a time-dependent bi-Kappa model thus becomes a valuable theoretical approach that can be incorporated into the present studies of solar wind dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Mechanism of Increasing Summer Water Vapor Over the Tibetan Plateau

Rapid warming over the Tibetan Plateau (TP) has been associated with a wetting trend in atmospheric water vapor content, which is critical for recharging the Asian water tower. However, the mechanism associated with the wetting phenomenon remains unclear. Long-term changes in moisture balance and precipitation recycling processes are investigated using the ERA5 reanalysis from 1979 to 2019. The wetting trend over TP is mainly due to the summer water vapor trend with signi?cantly increases over the western TP. Based on the moisture balance analysis, it is found that precipitation, evaporation, and the convergence of moisture are all increasing over the western TP but decreasing over the eastern TP. Based on the dynamical precipitation recycling model, the results suggest that both internal and external cycles contributes to the wetting over TP, with 63.87% and 36.13% contributions respectively. Further analysis found that the atmospheric heating source is also increasing over the western TP, which could shift the moisture transportation from east to west at the southern boundary of TP. The increasing moisture convergence could enhance precipitation, and the enhanced latent heating in the mid-atmosphere can further induced moisture convergence, which forms a positive feedback. However, an opposite situation occurred over the eastern TP. The internal and external cycle of the recycling process can stimulate (suppress) each other through precipitation over the western (eastern) TP. This mechanism linked the changes in the precipitation recycling and atmospheric circulation, and induced the wetting trend over TP in summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
DOE OSTI2022

Discovery of Innovative Polymers for Next-Generation Gas-Separation Membranes using Interpretable Machine Learning

Polymer membranes perform innumerable separations with far-reaching environmental implications. Despite decades of research on membrane technologies, design of new membrane materials remains a largely Edisonian process. To address this shortcoming, we demonstrate a generalizable, accurate machine-learning (ML) implementation for the discovery of innovative polymers with ideal separation performance. Specifically, multitask ML models are trained on available experimental data to link polymer chemistry to gas permeabilities of He, H2, O2, N2, CO2, and CH4. Here, we interpret the ML models and extract chemical heuristics for membrane design, through Shapley Additive exPlanations (SHAP) analysis. We then screen over nine million hypothetical polymers through our models and identify thousands of candidates that lie well above current performance upper bounds. Notably, we discover hundreds of never-before-seen ultrapermeable polymer membranes with O2 and CO2 permeability greater than 104 and 105 Barrer, respectively. These hypothetical polymers are capable of overcoming undesirable trade-off relationship between permeability and selectivity, thus significantly expanding the currently limited library of polymer membranes for highly efficient gas separations. High-fidelity molecular dynamics simulations confirm the ML-predicted gas permeabilities of the promising candidates, which suggests that many can be translated to reality.

Yang, Jason↗
DOE OSTI2024

The Acetylene Bridge in Intramolecular Singlet Fission: A Boon or A Nuisance?

Various analogues of the alkylsilylacetylene group are frequently used as auxiliary groups to enhance the solubility and stability of the acene dimer core, widely used as platforms to investigate intramolecular singlet fission (iSF) mechanisms. However, while in the 2,2’-linked dimers they are primarily auxiliary groups, these are essential fragments of the bridging units in 6,6’/5,5’-linked dimers, the two preferred choices for dimerization. The starkly different iSF dynamics observed in the two variants raise the question of what role the acetylene bridges play. Here, we systematically designed a set of (oligo-)para-phenylene bridged 2,2’-linked pentacene dimers with an additional acetylene fragment in the bridging unit to mimic the structure of 6,6’-linked dimers. Contrasting the results with previously reported analogous 2,2’-linked and 6,6’-linked pentacene dimers reveals that the acetylene bridges contribute to significant conformational freedom. This effect provides a mechanism to promote spin evolution within the triplet pair to achieve free triplets but also offers new parasitic pathways for triplet-pair recombination, revealing that this structural motif can be both a boon and a nuisance. Additionally, our analysis reveals that these bridges directly modify the electronic states, highlighting significant pitfalls of the standard chromophore-bridge-chromophore framework used to design and interpret photophysics of iSF materials.

14 SOLAR ENERGY↗
DOE OSTI2024

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Insights into the dynamics between viruses and their hosts in a hot spring microbial mat

Abstract Our current knowledge of host–virus interactions in biofilms is limited to computational predictions based on laboratory experiments with a small number of cultured bacteria. However, natural biofilms are diverse and chiefly composed of uncultured bacteria and archaea with no viral infection patterns and lifestyle predictions described to date. Herein, we predict the first DNA sequence-based host–virus interactions in a natural biofilm. Using single-cell genomics and metagenomics applied to a hot spring mat of the Cone Pool in Mono County, California, we provide insights into virus–host range, lifestyle and distribution across different mat layers. Thirty-four out of 130 single cells contained at least one viral contig (26%), which, together with the metagenome-assembled genomes, resulted in detection of 59 viruses linked to 34 host species. Analysis of single-cell amplification kinetics revealed a lack of active viral replication on the single-cell level. These findings were further supported by mapping metagenomic reads from different mat layers to the obtained host–virus pairs, which indicated a low copy number of viral genomes compared to their hosts. Lastly, the metagenomic data revealed high layer specificity of viruses, suggesting limited diffusion to other mat layers. Taken together, these observations indicate that in low mobility environments with high microbial abundance, lysogeny is the predominant viral lifestyle, in line with the previously proposed “Piggyback-the-Winner” theory.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Toward Decoding the Relationship between Domain Structure and Functionality in Ferroelectrics via Hidden Latent Variables

Polarization switching mechanisms in ferroelectric materials are fundamentally linked to local domain structure and the presence of the structural defects, which both can act as nucleation and pinning centers and create local electrostatic and mechanical depolarization fields affecting wall dynamics. However, the general correlative mechanisms between domain structure and polarization dynamics are only weakly explored, precluding insight into the associated physical mechanisms. In this work, the correlation between local domain structures and switching behavior in ferroelectric materials is explored using convolutional encoder–decoder networks, enabling image to spectral (im2spec) and spectral to image (spec2im) translations via encoding of latent variables. The latter reflect the assumption that the relationship between domain structure and polarization switching is parsimonious, i.e., is based upon a small number of local mechanisms. The analysis of latent variables distributions and their real-space representations provides insight into the predictability of the local switching behavior and hence associated physical mechanisms. We further pose that the regions where these correlative relationships are violated, i.e., predictability of the polarization dynamics from domain structure is reduced, represent the obvious target for detailed studies, e.g., in the context of automated experiments. This approach provides a workflow to establish the presence of correlation between local spectral responses and local structure and can be universally applied to spectral imaging techniques such as piezoresponse force microscopy (PFM), scanning tunneling microscopy (STM) and spectroscopy, and electron energy loss spectroscopy (EELS) in scanning transmission electron microscopy (STEM).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Regional and Temporal Variability of Atmospheric River Seasonality: Influences of Detection Algorithms and Moisture Transport Dynamics

Abstract Understanding the regional and temporal variability of atmospheric river (AR) seasonality is crucial for preparedness and mitigation of extreme events. While ARs were thought to peak in winter, recent research shows they exhibit region‐specific seasonality and are heavily influenced by the chosen detection algorithm. This study examines the link between the year‐to‐year consistency of peak‐AR activity to the presence of a dominant seasonal pattern, considering both location and algorithm choice. Regions are categorized by their temporal characteristics: consistent patterns (e.g., East Asia), patterns with occasional outliers (e.g., British Columbia coast), and regions lacking a clear dominant peak season (e.g., South Atlantic, parts of Australia). Hence, not all regions display a consistent seasonal cycle of AR activity. This study quantifies the extent to which a region experiences a dominant peak season of AR activity (or lacks one) and offers insights to enhance decision‐making in water management, natural hazard preparedness, and forecasting. Furthermore, given our finding that detection algorithms influence the peak season of AR activity, we also examine two diagnostic variables representative of moisture transport to corroborate our results. Integrated vapor transport, which captures meridional and zonal moisture transport, and Moist Wave Activity, representing moisture intrusions from lower to higher latitudes, are examined. Our analysis indicates that inconsistencies in the seasonal cycle of AR activity are not solely due to discrepancies in detection algorithms but also arise from changes in moisture transport. Plain Language Summary Atmospheric rivers (ARs) are critical weather phenomena that can cause extreme events like heavy rainfall and flooding. Understanding when and where ARs are most likely to occur throughout the year is essential for preparing and responding to these events. Traditionally, ARs were thought to peak in winter, but recent studies show this varies by region. Our study helps address the challenge decision‐makers face in anticipating and preparing for AR events by providing insights into the consistency of peak seasonal patterns across different areas. Some regions, like East Asia, and the British Columbia coast, show a consistent peak season, while others, like the South Atlantic and parts of Australia, have significant year‐to‐year variations, making it hard to identify a dominant season. To better understand these changes over time, the study also examines how moisture moves in the atmosphere, using Integrated Vapor Transport (which looks at moisture movement in various directions) and Moist Wave Activity (which tracks moisture shifts from lower to higher latitudes). The findings suggest that inconsistencies in AR patterns are due not only to detection methods but also due to changes in moisture transport. Key Points The peak season of atmospheric river activity can change depending on the year in some areas Interannual variations in the peak season can make identifying a dominant season challenging for some regions Frequent shifts in peak season across years reflect inconsistencies tied to detection algorithms and to underlying dynamics

Kamnani, Diya↗
DOE OSTI2026

Field-Coupled Water Splitting with Metal-Free Donor–Acceptor Covalent Organic-Framework Junctions

Advancing metal-free electrocatalysts for hydrogen and oxygen evolution reactions (HER/OER) across acidic and alkaline media requires coordinated control of intermediate binding thermodynamics, interfacial charge delivery, and near-electrode transport dynamics. Here, we design amide-linked benzene–triazine covalent organic frameworks (BTA/TzTA-Hz COFs) and integrate them with carbon nanotubes (CNTs) to form COF–CNT junctions that establish a built-in interfacial electric field. Density functional theory (DFT) and electrostatic potential maps indicate complementary active motifs, with benzene-proximal fragments associated with HER and triazine-proximal motifs associated with OER. CNT integration shifts the contact-potential difference by ≈0.20 V, while operando electrochemical impedance spectroscopy suggests partially separable high-frequency junction-charging and lower-frequency Faradaic/transport responses. A 300 mT static magnetic field lowers the HER and OER overpotentials by tens of millivolts. Under anodic bias, the effective interfacial charging capacitance increases, and Mott–Schottky analysis shows an apparent ∼0.15 V flat-band shift with an essentially unchanged slope. Together, these observations are consistent with field-perturbed interfacial charging and altered bias partitioning. Field-dependent impedance and bubble imaging are consistent with magnetohydrodynamic convection that promotes bubble detachment and near-electrode mass transport for both half-reactions, and they reveal an OER-specific high-frequency perturbation under anodic bias. Under field, the heterostructure reaches an OER onset overpotential of ∼261 mV and requires an overpotential of 366 mV at 10 mA cm –2 in alkaline electrolyte. These results illustrate how reticular-framework chemistry, junction engineering, and both built-in and applied fields can program reactivity through interfacial electrostatics and near-electrode transport in organic-framework electrocatalysts.

Garcia-Enriquez, Lissette [Univ. of Texas at El Pa↗
DOE OSTI2021

Grand Challenges of Hydrologic Modeling for Food-Energy-Water Nexus Security in High Mountain Asia

Climate-influenced changes in hydrology affect water-food-energy security that may impact up to two billion people downstream of the High Mountain Asia (HMA) region. Changes in water supply affect energy, industry, transportation, and ecosystems (agriculture, fisheries) and as a result, also affect the region's social, environmental, and economic fabrics. Sustaining the highly interconnected food-energy-water nexus (FEWN) will be a fundamental and increasing challenge under a changing climate regime. High variability in topography and distribution of glaciated and snow-covered areas in the HMA region, and scarcity of high resolution (in-situ) data make it difficult to model and project climate change impacts on individual watersheds. We lack basic understanding of the spatial and temporal variations in climate, surface impurities in snow and ice such as black carbon and dust that alter surface albedo, and glacier mass balance and dynamics. These knowledge gaps create challenges in predicting where and when the impact of changes in river flow will be the most significant economically and ecologically. In response to these challenges, the United States National Aeronautics and Space Administration (NASA) established the High Mountain Asia Team (HiMAT) in 2016 to conduct research to address knowledge gaps. This paper summarizes some of the advances HiMAT made over the past 5 years, highlights the scientific challenges in improving our understanding of the hydrology of the HMA region, and introduces an integrated assessment framework to assess the impacts of climate changes on the FEWN for the HMA region. The framework, developed under a NASA HMA project, links climate models, hydrology, hydropower, fish biology, and economic analysis. The framework could be applied to develop scientific understanding of spatio-temporal variability in water availability and the resultant downstream impacts on the FEWN to support water resource management under a changing climate regime.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Explainable machine learning reveals that local structural motifs encode the thermodynamic state across the CuZr metallic glass-forming range

Metallic glasses derive their properties from the statistics of local atomic motifs rather than from long-range order, yet a quantitative, chemistry-specific link between motif populations and the underlying glassy state has remained elusive. In this work we combine large-scale molecular dynamics, Voronoi tessellation, deep neural networks, and SHapley Additive exPlanations (SHAP) to identify which local structural motifs define the glassy state of Cu—Zr metallic glasses. A dataset of 17,180 atomistic configurations spanning ten compositions (Cu 20 Zr 80 –Cu 80 Zr 20 ) and four quench rates (10 9 –10 12 K/s) is used to train a feed-forward neural network that regresses temperature across the 50–2000 K liquid–supercooled–glass range, achieving a mean absolute error of 19.89 K and R 2 = 0.9974, confirming that the local structural state is faithfully encoded in motif-level structure. SHAP analysis then reveals that a tightly coupled near-icosahedral family of motifs (coordination numbers (CN) 11–13, including the full icosahedron 001200 and its single-atom-perturbation sibling 10930) collectively encodes the thermodynamic state of the system across the full glass-forming range. The CN = 11–13 ordered members carry negative SHAP values at high populations, tracking the most deeply-quenched configurations, while 10930 shows the reversed signature consistent with its role as a soft-spot host whose population shrinks as the icosahedral network deepens. The analysis demonstrates that explainable machine learning can isolate the minimal motif vocabulary defining the glassy state and recovers the near-icosahedral building blocks previously identified by data-driven analyses of Cu—Zr. The approach provides a general, chemistry-specific route for characterizing the structural state of disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Exploring Domain-Wall Pinning in Ferroelectrics via Automated High-Throughput Atomic Force Microscopy

Domain-wall dynamics in ferroelectric materials are strongly position-dependent, since each polar interface is locked into a unique local microstructure. This necessitates spatially resolved studies of wall pinning using scanning-probe microscopy techniques. The pinning centers and pre-existing domain walls are usually sparse within the image plane, precluding the use of dense hyperspectral imaging modes and requiring time-consuming human experimentation. Here, a large-area epitaxial PbTiO 3 film on cubic KTaO 3 was investigated to quantify the electric-field-driven dynamics of the polar–strain domain structures using ML-controlled automated piezoresponse force microscopy. Analysis of 1500 switching events reveals that domain-wall displacement depends not only on field parameters but also on the local ferroelectric–ferroelastic configuration. For example, twin boundaries in polydomains regions, like a 1 – /c+ ∥ a 2 – /c – , stay pinned up to a certain level of bias magnitude and change only marginally as the bias increases from 20 to 30 V, whereas single-variant boundaries, like the a 2 + /c + ∥ a 2 – /c – stack, are already activated at 20 V. These statistics on the possible ferroelectric and ferroelastic wall orientations, together with the automated high-throughput AFM workflow, can be distilled into a predictive map that links domain configurations to pulse parameters. Here, this microstructure-specific rule set forms the foundation for the design of ferroelectric memories.

automated scanning probe microscopy↗