Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orientation mapping”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

In-situ observations of cyclic deformation in an extruded Mg-2Nd-1Y-0.1Zr-0.1Ca alloy

In this study, the evolution of deformation mechanisms during cyclic loading in an extruded, solution-treated Mg–2Nd–1Y–0.1Zr–0.1Ca alloy was investigated using a combination of in-situ loading, scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), and focused ion beam (FIB) nanofabrication. The initial microstructure exhibited a random crystallographic texture with no preferred grain orientation. Flat, rectangular dog-bone specimens were subjected to load-controlled, fully reversed fatigue for 50 cycles, during which the same region was sequentially mapped to track microstructural changes. After 10 cycles of loading deformation twins were observed. During tensile reloading detwinning or narrowing of those twinned regions occurred. After 20 cycles, detwinning ceased and residual twins remained in the material. SEM imaging revealed numerous surface slip traces after cyclic loading. EBSD-assisted slip trace analysis identified the activation of prismatic and pyramidal < c+a> slip systems during low-cycle fatigue. Site-specific scanning transmission electron microscopy (STEM) further revealed that deformation was also accommodated by basal < a> slip and the dissociation of < c+a> dislocations. Center-of-symmetry (COS) analysis confirmed that the dissociation of < c+a> dislocations resulted in the formation of I₁ intrinsic stacking faults after cyclic loading. These findings provide new insights into the complex interplay of dislocation mechanisms governing fatigue deformation in rare-earth-containing Mg alloys.

Cyclic deformation↗
DOE OSTI2025

Kelvin probe force microscopy under ambient conditions

Kelvin probe force microscopy (KPFM) is a technique derived from atomic force microscopy that provides maps of surface potential or work function differences across material systems, with nanometre-scale resolution. KPFM is a useful tool for investigating electrical phenomena such as dipole orientation, interfacial charge transfer, charge accumulation, band bending and doping levels. This Primer aims to provide an overview of typical ambient-condition KPFM measurements, covering their underlying principles, experimental implementations and wide-ranging applications. Key KPFM variants, including amplitude and frequency modulation, heterodyne detection schemes and innovative open loop and pulsed force techniques, are discussed, with practical guidance on optimizing signal acquisition and reducing errors. Specialized approaches, such as time-resolved KPFM and multimodal KPFM, are discussed for their ability to capture dynamic charge processes and chemical information, respectively. Here, we highlight recent advances in KPFM applications, spanning metal alloys, soft matter, ferroelectrics, photovoltaics and 2D materials, showcasing its versatility across research domains. By addressing current limitations and identifying future opportunities, this Primer underscores the transformative potential of KPFM in advancing the understanding of nanoscale electrical phenomena.

Zahmatkeshsaredorahi, Amirhossein [Lehigh Univ., B↗
DOE OSTI2025

Poincaré beams from a free electron laser

Poincaré beams are light beams that have spatially inhomogeneous polarization structure that spans a finite portion of the Poincaré sphere. This feature bestows the beams with intriguing topological properties and has led to a surge in research on their fundamental characteristics, their controlled generation and on emerging applications. Here we present an experimental demonstration of a Poincaré beam generated in the extreme ultraviolet (16.7 nm) at the FERMI free electron laser (FEL). The ‘star’ type Poincaré beam is generated by exploiting the phase and intensity structure intrinsic to FEL radiation without relying on optical elements. Here, we controlled the spatial polarization distribution through a precise overlap and power balance between two FEL pulses, each with different transverse phase distributions and orthogonal circular polarizations. The spatial polarization structure was mapped in detail and shows extensive coverage of the Poincaré sphere, in agreement with analytic predictions. This method of in situ Poincaré beam production in FELs enables straightforward flexibility in the orientation and balance of polarization states, and can readily be extended to other vector beams and to shorter wavelengths enabling novel science applications in modern light sources.

Morgan, Jenny [SLAC National Accelerator Laborator↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Atomic-Scale Scanning of Domain Network in the Ferroelectric HfO 2 Thin Film

Ferroelectric HfO 2 -based thin films have attracted much interest in the utilization of ferroelectricity at the nanoscale for next-generation electronic devices. However, the structural origin and stabilization mechanism of the ferroelectric phase are not understood because the film is typically nanocrystalline with active yet stochastic ferroelectric domains. Here, in this study, electron microscopy is used to map the in-plane domain network structures of epitaxially grown ferroelectric Y:HfO 2 films in atomic resolution. The ferroelectricity is confirmed in free-standing Y:HfO 2 films, allowing for investigating the structural origin for their ferroelectricity by 4D-STEM, high-resolution STEM, and iDPC-STEM. At the grain boundaries of <111>-oriented Pca2 1 orthorhombic grains, a high-symmetry mixed-(R3m, Pnm2 1 ) phase is induced, exhibiting enhanced polarization due to in-plane compressive strain. Nanoscale Pca2 1 orthorhombic grains and their grain boundaries with mixed-(R3m, Pnm2 1 ) phases of higher symmetry cooperatively determine the ferroelectricity of the Y:HfO 2 film. It is also found that such ferroelectric domain networks emerge when the film thickness is beyond a finite value. Furthermore, in-plane mapping of oxygen positions overlaid on ferroelectric domains discloses that polarization is suppressed at vertical domain walls, while it is active when domains are aligned horizontally with subangstrom domain walls. In addition, randomly distributed 180° charged domain walls are confined by spacer layers.

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

radkit base v1.6

The radkit (base) software suite (python) consists of three primary libraries: stark, trajan, and curie. The trajan library provides the tools to analyze and manipulate data from lidar and inertial measurement unit (IMU) devices, cameras, as well as trajectories from algorithms such as simultaneous localization and mapping (SLAM). These components allow reading and writing standard data formats, performing rigid affine transformations, discretizing three-dimensional space, and visualizing data products. The curie library comprises a standard set of object-oriented tools for radiation data and analysis in the following modules: (1) listmode and binmode data classes with methods for manipulation, plotting, slicing and file IO; (2) radiological/nuclear source detection/identification analysis results; (3) source encounters of correlated analyses and (4) energy-dependent angular detector response functions. The stark package provides low-level tools that are leveraged by both curie and trajan. The tools are flexible for offline analysis as well as performant for real-time integrations.

Salathe, Marco [Lawrence Berkeley National Laborat↗
DOE OSTI2024

Automating the detection of hydrological barriers and fragmentation in wetlands using deep learning and InSAR

The loss of hydrological connectivity and fragmentation of natural wetlands is a widespread driver of wetland degradation. Understanding where and how natural connectivity is impaired is essential for managing, protecting and remediating these ecosystems. Wetland Interferometric Synthetic Aperture Radar (Wetland InSAR) can provide information on surface flow orientation in wetlands at a high spatial resolution, which can be used for barrier detection. However, the broad application of this approach is constrained by the labour-intensive manual delineation of barriers based on mapped water levels. This study presents the first deep learning-based methodology for the automated detection of hydrological barriers. We trained a deep convolutional network to segment edge features of hydrological barriers in 25 image pairs captured by ALOS PALSAR-1 L-Band InSAR between 2006 and 2011. The training dataset consists of manually labelled and delineated barriers showing abrupt changes in water surface elevation and wrapped interferograms with high coherence. We tested this method across three wetland sites: the Everglades and southern Louisiana wetlands (United States) and the Cienaga de Zapata (Cuba). Across these sites, the convolutional network detected hydrological barriers with up to 84% accuracy. The model performed particularly well for linear hydrological barriers such as roads, dikes, and channels. Notably, some barriers impede flow only seasonally, appearing during low water levels and disappearing when water levels rise. Our automated approach to detecting and assessing wetland hydrologic connectivity can be applied more broadly to support the effective management of fragmented wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Seasonal Precipitation Classification during Surface Atmosphere Integrated Field Laboratory Campaign

The Surface Atmosphere Integrated Field Laboratory (SAIL) campaign, conducted from September 2021 to June 2023 in Crested Butte, Colorado, aimed to characterize precipitation processes in the Upper Colorado River Basin (UCRB). This increased observations of snowfall accumulation in this hydrologically significant watershed would be useful for quantitative precipitation estimates (QPE). Therefore, the Surface Quantitative Precipitation Estimate (SQUIRE) product was developed using the ARM-supported Colorado State University (CSU) X-band Precipitation Radar. Although SQUIRE will be only released for snowfall, by categorizing precipitation types, users can effectively utilize relevant datasets under diverse meteorological conditions. Moreover, the dataset facilitates validation of the QPE product and the analysis of seasonal variations in precipitation types at the surface. Hydrometeors classes are organized based on their phase and physical characteristics mapping the CSU (both winter Summer) and Py-ART classifications into four groups. 1. Liquid Precipitation: includs drizzle, rain, and large raindrops. 2. Frozen Snow and Ice : Pure Snow, combining ice crystals, aggregates, and vertically oriented ice structures. 3. Dense and Large frozen hydrometeors: including low- and high-density graupel and dry hail. 4.Melting: Wet Snow and Melting Hail, hydrometeors exhibiting both liquid and frozen characteristics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗
DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗
DOE OSTI2025

Software Quality Assurance for the MOOSE-Based Open-Source Multiphysics Code Cardinal - An Expanded CI Testing Suite

Cardinal is a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and the Monte Carlo particle transport code OpenMC within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Cardinal provides high-resolution thermal-hydraulics and/or radiation transport feedback to MOOSE multiphysics simulations. Multiphysics feedback is implemented in a geometry-agnostic manner which eliminates the need for rigid one-to-one mappings. A generic data transfer implementation also allows NekRS and OpenMC to couple to any MOOSE application, enabling a broad set of multiphysics capabilities. Cardinal simulations can also leverage combinations of MPI, OpenMP, and GPU resources. Cardinal continuous development and improvement efforts have led to the software being considered as a high-fidelity design and licensing tool for key areas of nuclear reactor relevant physics, including neutron transport, fluid flow, heat transfer, and mechanical processes. The fast development and expansion of the software from a pure R&D framework towards its application in the nuclear industry and regulation require a focus on developing, enhancing and, maintaining Cardinal’s software quality through strict adherence to a Software Quality Assurance (SQA) framework and SQA program. To facilitate compliance with SQA standards, the Cardinal SQA Program has been initiated during Fiscal Year 2023 (FY23). During the development of the Cardinal SQA Program, multiple gaps have been identified. These gaps are primarily related to model verification and code pedigree as they relate to the use of Cardinal as a safety analysis tool. These gaps have been captured in a report published in 2023. A second report highlighted the progress made during Fiscal Year 2024 (FY24) and described Argonne’s effort to document and integrate software verification within Cardinal’s software development process. This report documents a snapshot of the verification test cases currently available for Cardinal and NekRS in their assimilation into a Continuous Integration (CI) platform. Following the CI practice permits the integrating of source code changes frequently and ensuring that the integrated codebase clears the verification testing for the software. It should be noted that the SQA program itself, including the program plans, procedures, configuration management, and testing strategies, need to be developed in a future step of this task.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Progress Towards NQA-1 for Cardinal in FY25

Cardinal is a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and the Monte Carlo particle transport code OpenMC within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Cardinal provides high-resolution thermal-hydraulics and/or radiation transport feedback to MOOSE multiphysics simulations. Multiphysics feedback is implemented in a geometry-agnostic manner which eliminates the need for rigid one-to-one mappings. A generic data transfer implementation also allows NekRS and OpenMC to couple to any MOOSE application, enabling a broad set of multiphysics capabilities. Cardinal simulations can also leverage combinations of MPI, OpenMP, and GPU resources. Cardinal continuous development and improvement efforts have led to the software being considered as a high-fidelity design and licensing tool for key areas of nuclear reactor relevant physics, including neutron transport, fluid flow, heat transfer, and mechanical processes. The fast development and expansion of the software from a pure R&D framework towards its application in the nuclear industry and regulation require a focus on developing, enhancing,and maintaining Cardinal’s software quality through strict adherence to a Software Quality Assurance (SQA) framework and SQA program. To facilitate compliance with SQA standards, the Cardinal SQA Program was initiated during Fiscal Year 2023 (FY23). During the development of the Cardinal SQA Program, multiple gaps have been identified. These gaps are primarily related to model verification and code pedigree as they relate to the use of Cardinal as an analysis tool. These gaps were captured in a report published in 2023. A second report highlighted the progress made during Fiscal Year 2024 (FY24) and described Argonne’s effort to document and integrate software verification within Cardinal’s software development process. This report documents the progress made towards NQA-1 for Cardinal in the Fiscal Year 2025 (FY25). All cases in the expanded Continuous Integration (CI) suite of NekRS are included in this report which test the solvers and modules available in NekRS exhaustively. The NekRS tests are integrated with the Cardinal CI suite and made available in publicly accessible Github documentation. Following the CI practice permits integrating of source code changes frequently and ensuring that the integrated codebase clears the verification testing for the software. Also in this report is a brief overview of the development of the Cardinal Software Quality Assurance Plan (SQAP) that was done in FY25, though it should be noted that the rest of the documentation for the SQA program needs to be developed in a future step of this task.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Effects of break geometry and orientation on helium-air mixing in simulated reactor cavities of high temperature gas reactors

Here, this study experimentally examined the spatial and temporal variations in air and helium concentrations and temperature fields within simulated reactor cavities of a High Temperature Gas Reactor (HTGR) following helium discharge into an initially air-filled reactor cavity system. Detailed temperature maps were generated using a combination of fiber optics temperature sensor and multiple thermocouple probes within the simulated reactor cavities. The research scenario involved a hypothetical small pipe break in the Reactor Pressure Vessel, resulting in the release of high-temperature helium into the surrounding cavity. A scaled multi-compartment experimental facility, modeled after the General Atomics Modular High Temperature Gas Reactor (GA-MHTGR) design, was constructed for helium and air mixing experiments. Oxygen sensors and thermocouple probes were installed in all five cavities to measure the concentrations of oxygen (or helium) and the temperature distributions of the gas mixture. The experimental findings highlighted the significant impact of the injected helium jet velocity on the gas mixing process and demonstrated how the direction of the helium jet influences the air-helium temperature profiles within the cavities.

Air-ingress↗
DOE OSTI2025

Classical eikonal from Magnus expansion

In a classical scattering problem, the classical eikonal is defined as the generator of the canonical transformation that maps in-states to out-states. It can be regarded as the classical limit of the log of the quantum S-matrix. In a classical analog of the Born approximation in quantum mechanics, the classical eikonal admits an expansion in oriented tree graphs, where oriented edges denote retarded/advanced worldline propagators. The Magnus expansion, which takes the log of a time-ordered exponential integral, offers an efficient method to compute the coefficients of the tree graphs to all orders. We exploit a Hopf algebra structure behind the Magnus expansion to develop a fast algorithm which can compute the tree coefficients up to the 12th order (over half a million trees) in less than an hour. In a relativistic setting, our methods can be applied to the post-Minkowskian (PM) expansion for gravitational binaries in the worldline formalism. We demonstrate the methods by computing the 3PM eikonal and find agreement with previous results based on amplitude methods. Importantly, the Magnus expansion yields a finite eikonal, while the naïve eikonal based on the time-symmetric propagator is infrared-divergent from 3PM on.

Black Holes↗
DOE OSTI2024

Self-Assembly of a Triblock Copolymer in the Presence of a Rigid Conjugated Polyelectrolyte

The properties of conducting polymers are strongly influenced by structural changes induced by long-range order, which can be achieved by using block copolymers that self-assemble into crystalline structures. These blends result in unique mesophases distinct from the pure components with self-assembly behavior modulated by solution conditions and polymer architectures. High-throughput small-angle X-ray scattering data of aqueous Pluronic P123 (PEO 20 –PPO 70 –PEO 20 ) and conjugated polyelectrolyte poly[3-(potassium-4-butanoate) thiophene-2,5-diyl] (PPBT) blends at various concentrations and temperatures were automatically classified into phase maps by autophasemap, an unsupervised statistical analysis algorithm. The outlined phase boundaries revealed that adding PPBT to high P123 concentrations induced a transition from cubicly ordered spherical micelles to hexagonally packed cylindrical micelles. Shear alignment via rheological small-angle neutron scattering of the blends produced monolithic oriented cubic and hexagonal crystal gels. Furthermore, these insights into the self-assembly of conductive polymer blends will aid in the design of soft materials with tunable structural and electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Autonomous Nanoparticle Synthesis Guided by In Situ Multiscale Structural Characterization

Autonomous synthesis platforms promise rapid exploration of vast parameter spaces; yet, integrating in situ structural characterization in closed-loop synthesis optimization remains challenging. We demonstrate a realization of such a closed-loop platform coupled with a droplet-flow microreactor, in situ X-ray scattering methods (SAXS/WAXS), and Gaussian process optimization to synthesize citrate-reduced Au nanoparticles with targeted characteristics. The system efficiently explored ∼19,000 synthesis recipes through 365 experiments, achieving precise control over size (4–60 nm) and polydispersity (σ < 0.11) across large citrate/gold ratios, exceeding traditional synthesis boundaries (1–10). Beyond confirming classical Turkevich–Frens trends, partial-dependence analysis revealed strong nonlinear coupling among precursor, citrate, and pH effects. Combining quantitative SAXS/WAXS analysis with electron microscopy characterization, we uncovered that crystallite size (d c ) and particle size (d) follow d c = 0.18d + β, where synthesis chemistry controls the intercept β while maintaining a universal slope. This parallel-band structure enables independent tuning of crystallite domain size at fixed particle diameter through a combination of chloride, gold precursor, citrate, and pH contributions (cross-validated Spearman ρ = 0.7 ± 0.1). High-resolution electron microscopy shows multiple lattice-fringe orientations within single particles, directly confirming polycrystalline domains and the ability to tune d c at the fixed d. The platform’s validation includes indistinguishable static versus flowing measurements, stable droplet transport at 100 °C, and <5% run-to-run variation, establishing a robust framework for mapping and controlling multiscale nanoparticle structure across expansive chemical spaces. In conclusion, the developed closed-loop platform can be applied to a borad range of nanosyntheis processes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

3D Radio Frequency Mapping and Polarization Observations Show Lightning Flashes Were Ignited by Cosmic‐Ray Showers

Previous 2D radio frequency interferometric observations showed that lightning is commonly started with a positive fast discharge (+FD) before a normal negative leader continues from the origin of the +FD to the ensuing lightning flash. However, the inception and development of the +FD cannot be convincingly explained by existing discharge theories. With our new 3D broadband interferometric mapping and polarization system, we observed that the +FD was sometimes followed by an even faster and more extensive negative discharge (−FD) that propagated backward and overshoot the origin by a few hundred meters, as reported by an earlier 2D study. Surprisingly, the signal polarization, which measures the orientation of the discharge current, systematically slanted from the discharge propagation direction and rotated between the two opposite discharges, showing the +FD and −FD were driven by other storm-independent factors in addition to the storm electric field, or otherwise all would align in the same direction. Assuming a cosmic-ray shower (CRS) piercing through the cloud immediately before the +FD/−FD discharge, we found that their path is consistent with a pre-ionized path by the CRS, and the polarizations for the two opposite discharges are consistent with the respective deflected trajectories of high-energy positrons and electrons in the geomagnetic and an electric field. We further analyzed the more commonly observed +FD and showed that it is consistent with the CRS interpretation, suggesting these flashes in thunderstorms were ignited by cosmic-ray showers.

58 GEOSCIENCES↗
DOE OSTI2026

Learning interpretable surface elasticity properties from bulk properties via neural network equation learners

Surface elasticity is central to understanding the mechanics and stability of surfaces and interfaces. It is characterized by quantities such as surface tension, residual surface stress, and surface stiffness. However their analytical expressions are typically difficult to derive from atomistic data, and depend strongly on modeling choices. This work presents a neural network-based equation learner which combines customized activation functions and connection-based pruning to discover parsimonious, closed-form equations for surface elasticity from atomistic simulations. Applying the method to seven face-centered cubic (FCC) metals, our equation learner uncovers interpretable equations that describe both low-Miller index and high-Miller index surface properties, capturing long-tail property distributions accurately. The discovered expressions are decoupled into two components: a universal, geometry-driven orientation function, and material-specific baseline coefficients. We find that lower-order properties such as surface tension are fundamentally geometry dependent, while higher-order properties such as surface stress and elasticity show more complex geometry and material dependence. We also relate material dependent coefficients to bulk properties, forming a clear map from bulk material properties to surface elasticity. Overall, this approach demonstrates that interpretable neurosymbolic machine learning can bridge the gap between atomistic simulations and physical laws, enabling the discovery of generalizable structure–property relationships for materials science phenomena such as surface elasticity.

Equation learning↗