Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability characterization”

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

The Effect of Material Properties on Oxygen Evolution Activity and Assessing Half-Cell Screening as a Predictive Tool in Electrolysis

Iridium-based oxygen evolution catalysts are screened in this study for activity and stability with rotating disk electrode (RDE) half-cells. This study focuses on the electrochemical and materials approaches needed to characterize oxygen evolution catalysts, and include testing for activity, stability, composition, oxide content, and structure. Findings also discuss recommendations for data interpretation and detail the difficulties of comparing catalysts across materials sets with different elemental and oxide compositions, and linking RDE activity to device-level performance. The materials evaluated are a mixture of oxides and metals, and several methods are used to quantify metal content, qualitatively assess oxide content, and determine total surface area. Oxygen evolution activities and stabilities are compared, where a wide range of results are reported. In general, higher RDE performances are found for catalysts that contained larger amounts of ruthenium and metals. Higher durability, however, is found for catalysts that only contained iridium and a higher proportion of oxides. Additionally, catalysts are evaluated for performance in membrane electrode assemblies to assess RDE as a predictive tool in electrolysis. While activity trends within individual material sets generally held between ex- and in-situ testing, RDE tends to overestimate the activity of more metallic catalysts when compared to device-level performance. These results stress the need for multiple metal/oxide baselines for mixed catalysts, to better project in-situ kinetics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multifunctional Tunable Polymethacrylates for Enhanced Shear Stability and Wear Prevention

This paper describes the design, synthesis, and characterization of linear, polar, and non-polar polymethacrylates designed for use as lubricant viscosity index improvers with enhanced shear and anti-wear resistance. The polymers were prepared via Reversible Addition-Fragmentation Chain Transfer (RAFT) to obtain relatively low, narrow, molecular weights, with random polar moieties. The resulting polymers were evaluated as oil solutions to determine their viscosity index, shear stability performance, and antiwear characteristics. The polymers had molecular weights (Mw) ranging from 120-170 kDa, and the resulting viscosity indexes (VIs) in general tracked the Mw. Although the VIs were modest (140-189), all polymers displayed higher values than a commercial product designed for similar applications, both at 2% (w/w) and when kinematic viscosities (KVs) were normalized to the lowest KV100 of ~5cSt. One of the imidazole-containing polymers displayed an anomalously high VI for a low apparent Mw. As hypothesized, many of the polar polymers demonstrated substantial wear reduction compared to the non-polar homopolymers or a commercial benchmark. The low-molecular-weight imidazole-containing polymer produced less than 4% of the wear displayed by the commercial standard. The higher molecular weight polymers containing imidazole, hydroxy group and amino group also produced as little as 10% of the wear shown by the benchmark. The strategy demonstrated the benefit of low, narrow, molecular weight polymers to achieve good shear stability in viscosity modifiers while minimizing wear, rendering them suitable for fluid-power applications.

friction modifiers↗
DOE OSTI2023

Regulating Anionic Redox via Mg Substitution in Mn-Rich Layered Oxide Cathodes Enabling High Electrochemical Stability for Sodium-Ion Batteries

Here, with the limited resources and high cost of lithium-ion batteries (LIBs) and the ever-increasing market demands, sodium-ion batteries (SIBs) gain much interest due to their economical sustainability, and similar chemistry and manufacturing processes to LIBs. As cathodes play a vital role in determining the energy density of SIBs, Mn-based layered oxides are promising cathodes due to their low cost, environmental friendliness, and high theoretical capacity. However, the main challenge is structural instability upon cycling at high voltage. Herein, Mg is introduced into the P2-type Na 0.62 Ni 0.25 Mn 0.75 O 2 cathode to enhance electrochemical stability. By combining electrochemical testing and material characterizations, it is found that substituting 10 mol% Mg can effectively alleviate the P2–O2 phase transition, Jahn-Teller distortion, and irreversible oxygen redox. Moreover, structural integrity is greatly improved. These lead to enhanced electrochemical performances. With the optimized sample, a remarkable capacity retention of 92% in the half cell after 100 cycles and 95% in the full cell after 170 cycles can be achieved. Altogether, this work provides an alternative way to stabilize P2-type Mn-based layer oxide cathodes, which in turn, put forward the development of this material for the next-generation SIBs.

25 ENERGY STORAGE↗
DOE OSTI2025

Baseband Digital Network Analyzer Upgrade for LLRF Controllers

Digital Network Analyzers (DNA) have been implemented in many Low-Level Radio Frequency (LLRF) systems, notably NSLS-II and CERN, to help tune feedback loops. DNA characterizes feedback loops by measuring the frequency-dependent magnitude and phase transfer functions. It enables the measurement of open loop gains, gain/phase margins, and loop delays to help fine-tune feedback loops. An FPGA-based DNA has been developed and integrated into the current Relativistic Heavy Ion Collider (RHIC) LLRF infrastructure. Its performance has been tested with an implementation of one-turn delay feedback (OTFB) on the bench to maximize gain and stability. The DNA has been used to characterize a RHIC 28 MHz cavity in a RHIC Accelerator Physics Experiment (APEX) to test transient beam loading compensation strategies.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Spectroscopic Identification of Carbamate Formation and Synergistic Binding in Amide–CO$^{–}_{2}$ Complexes

Carbamate formation is an elementary step that governs nitrogencentered nucleophilic CO 2 capture across diverse environments, yet a direct, structure-specific understanding has been lacking. Here, we report the first gas-phase characterization of closed-shell carbamate formation by deprotonated amides, using cryogenic ion trap vibrational spectroscopy combined with quantum chemical calculations. Reactions of deprotonated benzamide and isophthalamide anions with CO 2 form carbamate species, with the amide nitrogen serving as the nucleophilic site. Diagnostic, strongly red-shifted antisymmetric CO 2 stretching vibrational bands, supported by a bonding analysis, establish chemisorption with substantial charge transfer. In the multiamide system, an intramolecular N−H···O hydrogen bond provides synergistic stabilization, correlating with larger red shifts and more exothermic binding. These structure-assigned benchmarks provide molecular-level insights into the binding motif, charge redistribution, and hydrogen bond-mediated stabilization of amide carbamates, which aid the characterization of carbamate formation in condensed phase.

Amides↗
DOE OSTI2023

Operation of Layered LiCoO 2 to Higher Voltages with a Localized Saturated Electrolyte

The commercialization of lithium-ion batteries started with the layered LiCoO 2 (LCO) cathode for portable electronics, but only 50% of its theoretical capacity can be used in practical cells due to detrimental surface and bulk degradations when charged to high voltages. Here we demonstrate here that the stability of electrolyte plays a critical role in the performance of LCO at high voltages by employing a localized saturated electrolyte (LSE). With a cutoff voltage of 4.5 V, LCO achieves an initial 1C discharge capacity of 176 mA h g -1 and a capacity retention of 80% over 230 cycles. Even with a cutoff voltage of 4.6 V, LCO with 2% aluminum doping displays an initial discharge capacity of 189 mA h g -1 at 1C rate with 80% capacity retention over 137 cycles. With extensive analytical characterization, we show that the improved cycling stability stems from a suppression of the O3 to H1-3 phase transition as well as a robust inorganic-rich cathode-electrolyte interphase (CEI) facilitated by the LSE. This work highlights the importance of protecting the surface as well as the bulk with appropriate electrolytes for the high-voltage, higher-energy-density operation of LCO.

25 ENERGY STORAGE↗
DOE OSTI2021

Understanding Degradation Mechanisms in SrIrO 3 Oxygen Evolution Electrocatalysts: Chemical and Structural Microscopy at the Nanoscale

Designing acid-stable oxygen evolution reaction electrocatalysts is key to developing sustainable energy technologies such as polymer electrolyte membrane electrolyzers but has proven challenging due to the high applied anodic potentials and corrosive electrolyte. This work showcases advanced nanoscale microscopy techniques supported by complementary structural and chemical characterization to develop a fundamental understanding of stability in promising SrIrO 3 thin film electrocatalyst materials. Cross-sectional high-resolution transmission electron microscopy illustrates atomic-scale bulk and surface structure, while secondary ion mass spectrometry imaging using a helium ion microscope provides the nanoscale lateral elemental distribution at the surface. After accelerated degradation tests under anodic potential, the SrIrO 3 film thins and roughens, but the lateral distribution of Sr and Ir remains homogeneous. A layer-wise dissolution mechanism is hypothesized, wherein anodic potential causes the IrO x -rich surface to dissolve and be regenerated by Sr leaching. Finally, the characterization approaches utilized herein and mechanistic insights into SrIrO 3 are translatable to a wide range of catalyst systems.

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

Implementational Aspects for the Characterization of Low-Inertia in Power Systems

The replacement of conventional generating units with non-conventional renewable energy resources impacts the stability of the power system. While the general goal of supplying demand could be achieved from both generation resources, inverter-based energy resources are different from traditional synchronous machines in the sense that they do not inherently provide inertia to the system. The latter obviously worsen the system ability to maintain the balance between supply and demand after certain contingencies such as generator outages. With high integration of renewable-based generation resources, there is a reduction in total system inertia, which could lead to systems with less time to react for frequency restoration in case of such unexpected unit failures. While conceptually it is true, there is a missing gap in the field of power system stability community in terms of test-bed power systems for characterizing low-inertia conditions. In this paper, a framework for low-inertia power systems is presented. The purpose of the test-bed is to establish an environment for testing and conducting research on power systems experiencing low-inertia issues. The frequency nadir, the rate of change of frequency, and the under-frequency load shedding scheme are used as a measure of the healthiness of power grids.

Baquedano-Aguilar, Mario D.↗
DOE OSTI2025

Initial Characterization of the NREL Large-Amplitude Motion Platform

The Large Amplitude Motion Platform (LAMP) at NREL represents a significant advancement in the controlled testing of Wave Energy Converters (WECs) under laboratory conditions. Originally designed by E2M as a six-degree-of-freedom (DOF) Stewart platform for flight simulation, LAMP has been adapted by NREL to facilitate the mounting and evaluation of WECs. This adaptation enables dry testing of WECs using motion profiles similar to the ocean, facilitating the iterative design, testing, and validation of WEC performance prior to ocean deployments. This report presents the initial work completed to characterize LAMP, with particular emphasis on its stability and operational capabilities across various single and multi-degree-of-freedom (DOF) motion profiles. The report includes planned comparisons at three distinct mass payloads, aimed at assessing the platform's positional accuracy, frequency response, and endurance over extended runtime periods. These experimental tests are critical for establishing the platform's limitations and ensuring that the data generated during WEC validation is both accurate and reproducible. The outcomes of this study not only contribute to a deeper understanding of LAMP's capabilities but also lay the groundwork for future advancements in WEC testing methodologies. By providing robust and reliable performance data within a controlled laboratory setting, the findings are expected to significantly enhance the development and commercialization of marine energy technologies. This report presents the initial findings of the LAMP Characterization work and proposed steps to further understand and characterize LAMP. Data collected during this work can be found on MHKDR at: https://mhkdr.openei.org/submissions/602 Note that this report shares the measured/found instantaneous maximum operating range of LAMP. For most applications the maximum operating range cannot be used for system health and longevity. The operating range and capabilities of LAMP will be evaluated on a case-by-case basis, single-DOF position, velocity, and acceleration values presented in Table 4a-c, and Table 5 should be taken as instantaneous absolute maximum values. Future use of LAMP will likely be limited to smaller values.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

13 C ENDOR Characterization of the Central Carbon within the Nitrogenase Catalytic Cofactor Indicates That the CFe 6 Core Is a Stabilizing “Heart of Steel”

Substrates and inhibitors of Mo-dependent nitrogenase bind and react at Fe ions of the active-site FeMo-cofactor [7Fe–9S–C–Mo–homocitrate] contained within the MoFe protein α-subunit. The cofactor contains a CFe 6 core, a carbon centered within a trigonal prism of six Fe, whose role in catalysis is unknown. Targeted 13 C labeling of the carbon enables electron-nuclear double resonance (ENDOR) spectroscopy to sensitively monitor the electronic properties of the Fe–C bonds and the spin-coupling scheme adopted by the FeMo-cofactor metal ions. Here this report compares 13 CFe 6 ENDOR measurements for (i) the wild-type protein resting state (E 0 ; α-Val 70 ) to those of (ii) α-Ile 70 , (iii) α-Ala 70 -substituted proteins; (iv) crystallographically characterized CO-inhibited “hi-CO” state; (v) E 4 (4H) Janus intermediate, activated for N 2 binding/reduction by accumulation of 4[ e– /H + ]; (vi) E 4 (2H)* state containing a doubly reduced FeMo-cofactor without Fe-bound substrates; and (vii) propargyl alcohol reduction intermediate having allyl alcohol bound as a ferracycle to FeMo-cofactor Fe6. All states examined, both S = 1/2 and 3/2 exhibited near-zero 13 C isotropic hyperfine coupling constants, $^Ca$ = [-1.3 ↔ +2.7] MHz. Density functional theory computations and natural bond orbital analysis of the Fe-C bonds show that this occurs because a (3 spin-up/3 spin-down) spin-exchange configuration of CFe 6 Fe-ion spins produces cancellation of large spin-transfers to carbon in each Fe–C bond. Previous X-ray diffraction and DFT both indicate that trigonal-prismatic geometry around carbon is maintained with high precision in all these states. The persistent structure and Fe–C bonding of the CFe 6 core indicate that it does not provide a functionally dynamic (hemilabile) “beating heart”-instead it acts as “a heart of steel”, stabilizing the structure of the FeMo-cofactor-active site during nitrogenase catalysis.

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

Reactive Extrusion of Waste Plastics with Compatibilizer and Lightly Pyrolyzed Crumb Rubber for Asphalt Modification

Modifying asphalt is a potentially high-value application for reusing waste plastics because of the high-volume usage of asphalt in highway construction. However, simply blending hot plastics and asphalt encounters difficulties related to the poor solubility of polymers, which limits the formation of a swollen network with asphalt molecules. The polymer phases also tend to coalesce and separate from asphalt during high-temperature storage in static conditions. The present study developed an innovative process to stabilize waste plastics in asphalt and improve binder storage stability by using lightly pyrolyzed crumb rubber together with a chemical compatibilizer. Both polymers were extruded to produce a thermoplastic elastomer (TPE) for asphalt modification. The mechanical performance and chemical reactions of TPEs were characterized via tension test and Fourier transform infrared spectroscopy. The storage stability and rheological properties of modified binder blends were evaluated through laboratory asphalt stability test and dynamic shear rheometer test. Polymer phases and network structures were characterized through optical microscopy. It was found that the pyrolyzed and reactive extrusion process improved the rubber solubility and polymer interaction, and therefore the storage stability of modified binder blends. The co-existence of rigid plastic and soft rubbery regimes in an entangled network provided a promising pathway to improve the mechanical performance of asphalt binders in both high- and low-temperature domains.

Ma, Yuetan↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

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

Commissioning the H i Observing Mode of the Beam Former for the Cryogenically Cooled Focal L-band Array for the GBT (FLAG)

We present the results of commissioning observations for a new digital beam-forming back end for the Focal plane L-band Array for the Robert C. Byrd Green Bank Telescope (FLAG), a cryogenically cooled Phased Array Feed (PAF) with the lowest measured T {sub sys}/η of any PAF outfitted on a radio telescope to date. We describe the custom software used to apply beam-forming weights to the raw element covariances to create research-quality spectral-line images for the new fine-channel mode, study the stability of the beam weights over time, characterize FLAG’s sensitivity over a frequency range of 150 MHz, and compare the measured noise properties and observed distribution of neutral hydrogen emission from several extragalactic and Galactic sources with data obtained with the current single-pixel L-band receiver. These commissioning runs establish FLAG as the preeminent PAF receiver currently available for spectral-line observations on the world’s major radio telescopes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Cu(II) Stability and UV-Induced Electron Transfer in a Metal–Organic Hybrid: An EPR, DFT, and Crystallographic Characterization of Copper-Doped Zinc Creatininium Sulfate

Single-crystal X-ray diffraction and electron paramagnetic resonance (EPR) spectroscopic experiments, complemented by quantum chemical DFT calculations, were carried out on the copper-doped metal–organic hybrid and Tutton salt analogue zinc creatininium sulfate to determine its crystal structure, to characterize the electronic structure of the doped Cu(II) binding site, and to propose a pathway for an excited-state, proton-coupled electron transfer (PCET) process in UV-exposed crystals. The crystal structure is isomorphous to that of cadmium creatininium sulfate, which has the transition ion, not in direct coordination with the creatinine, but forming a hexahydrate complex, which is bridged to a creatininium through an intervening sulfate ion. The EPR g (2.446, 2.112, 2.082) and copper hyperfine (A Cu : -327, -59.6, 10.8 MHz) tensor parameters are consistent with doped copper replacing host zinc in the metal–hexahydrate complex. These parameters are similar to those observed for copper hexahydrate in doped Tutton salt systems at low temperature, where the unpaired electron occupies mainly the copper 3d x 2 –y 2 orbital. At room temperature in the Tutton systems, vibration couplings stemming from a dynamic Jahn–Teller effect cause tensor averaging which results in a reduction in their maximum g-tensor and hyperfine tensor values. However, like for the doped isomorphous Cd creatinine crystal, the Cu(II) EPR exhibits little, or no room temperature averaging compared to its low temperature pattern. Samples exposed to 254 nm UV light generate a carbon-centered free radical species, characterized by an isotropic g-tensor (g = 2.0029) and an alpha-proton hyperfine coupling (-24 -14 +4 G). These parameters identify it as a creatinine radical cation formed by the oxidative release of one of its C2 methylene hydrogens. DFT calculations confirm the unpaired electronic structures of both the Cu(II) site and free radical. The growth in radical concentration with an increase in the UV exposure time coincides with a decrease in the copper EPR signal, indicating a coupled light-induced oxidation reduction process. A comparison of the crystal structure with the EPR parameters and DFT results provides evidence for a UV-induced PCET.

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

Atomically Dispersed Dual‐Metal Site Catalysts for Enhanced CO 2 Reduction: Mechanistic Insight into Active Site Structures

Abstract Carbon‐supported nitrogen‐coordinated single‐metal site catalysts (i.e., M−N−C, M: Fe, Co, or Ni) are active for the electrochemical CO 2 reduction reaction (CO 2 RR) to CO. Further improving their intrinsic activity and selectivity by tuning their N−M bond structures and coordination is limited. Herein, we expand the coordination environments of M−N−C catalysts by designing dual‐metal active sites. The Ni‐Fe catalyst exhibited the most efficient CO2RR activity and promising stability compared to other combinations. Advanced structural characterization and theoretical prediction suggest that the most active N‐coordinated dual‐metal site configurations are 2N‐bridged (Fe‐Ni)N 6 , in which FeN 4 and NiN 4 moieties are shared with two N atoms. Two metals (i.e., Fe and Ni) in the dual‐metal site likely generate a synergy to enable more optimal *COOH adsorption and *CO desorption than single‐metal sites (FeN 4 or NiN 4 ) with improved intrinsic catalytic activity and selectivity.

Li, Yi↗
DOE OSTI2024

SEC ‐ SAXS / MC Ensemble Structural Studies of the Microtubule Binding Protein Cdt1 Show Monomeric, Folded‐Over Conformations

ABSTRACT Cdt1 is a mixed folded protein critical for DNA replication licensing and it also has a “moonlighting” role at the kinetochore via direct binding to microtubules and the Ndc80 complex. However, it is unknown how the structure and conformations of Cdt1 could allow it to participate in these multiple, unique sets of protein complexes. While robust methods exist to study entirely folded or unfolded proteins, structure–function studies of combined, mixed folded/disordered proteins remain challenging. In this work, we employ orthogonal biophysical and computational techniques to provide structural characterization of mitosis‐competent human Cdt1. Thermal stability analyses shows that both folded winged helix domains1 are unstable. CD and NMR show that the N‐terminal and linker regions are intrinsically disordered. DLS shows that Cdt1 is monomeric and polydisperse, while SEC‐MALS confirms that it is monomeric at high concentrations, but without any apparent inter‐molecular self‐association. SEC‐SAXS enabled computational modeling of the protein structures. Using the program SASSIE, we performed rigid body Monte Carlo simulations to generate a conformational ensemble of structures. We observe that neither fully extended nor extremely compact Cdt1 conformations are consistent with SAXS. The best‐fit models have the N‐terminal and linker disordered regions extended into the solution and the two folded domains close to each other in apparent “folded over” conformations. We hypothesize the best‐fit Cdt1 conformations could be consistent with a function as a scaffold protein that may be sterically blocked without binding partners. Our study also provides a template for combining experimental and computational techniques to study mixed‐folded proteins.

Cell Biology↗
DOE OSTI2025

Cross-scale modeling and experimental integration for advancing cathode electrolyte interphase studies in high energy density lithium-ion batteries

Electrochemical interfaces are critical to the performance and durability of lithium-ion batteries (LIBs). The solid electrode-electrolyte interphase (SEI and CEI) structures that form during cycling can passivate reactive surfaces, ensuring safe operation, but also may contribute to performance degradation. Understanding the microscopic factors influencing interphase formation, growth, and evolution is essential for balanced battery design. While significant research has focused on the anode-electrolyte interphase (SEI), the cathode-electrolyte interphase (CEI) remains less explored, despite its importance in high-voltage and advanced battery technologies. Challenges in conducting in-situ or operando experiments arise from the occluded nature of these interfaces and the long timescales involved, often leading to biased interpretations. A validated multi-scale, multi-physics modeling approach, integrated with advanced characterization techniques, can effectively elucidate the intrinsic stability of electrolyte and cathode surfaces, the impact of chemical heterogeneity, and the role of microstructural features on CEI performance. In conclusion, this article reviews current modeling and simulation strategies for studying CEI in advanced LIBs and highlights opportunities for future methodological advancements and experimental integration.

Cathode-electrolyte interphase↗
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↗