Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Common data models”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Extending CF Conventions to Enhance Data FAIRness for Atmospheric Composition Observations

The Hierarchical Data Format (HDF) and Network Common Data Form (NetCDF) are data file formats created to aid users in the creation or use of scientific data. These file formats are useful for handling large data volumes and hosting extensive metadata as global, group, or variable attributes and are popular with the modeling community. HDF and NetCDF files are widely used with atmospheric remote sensing data and have been used to support measurements from numerous field campaigns, from satellite to aircraft or ground and mobile based measurements. The files from airborne field studies, however, vary greatly in terms of the file structure and the amount and content of their metadata. Information relevant to the file that can be useful to the user such as the data producer, location where data was taken, variable descriptions, or information about the instrument might not be included in the file. Recently, the Measurements of Aerosols, Clouds, and their Interactions for Earth System Models (MACIE) group started a grassroots effort to develop a CF-based template for the HDF and NetCDF files for field studies, with the aim of making the data products more interoperable and usable. This template seeks to make the files more compliant to Climate and Forecast (CF) metadata conventions and to standardize the file structure and the global and variable attributes. The template would help to ensure that HDF and NetCDF files contain adequate metadata to better support their use for research, e.g., the modeling community, and to enhance the usability and interoperability of data for research communities at large. The draft template has been applied to recent field studies for various instruments and their merge files in support of the Atmosphere Observing System (AOS) project. The details of the revised template are to be presented, as well as examples of the implementation of these requirements for merge files and lidar observation data files and issues revealed during the implementation process.

Sean Leavor↗
NASA NTRS2006

Exploring the Production of NOx by Lightning and Its Impact on Tropospheric Ozone

Our quantitative understanding of free tropospheric (FT) chemistry is quite poor. State-of-the-art regional air quality models (e.g., US EPA's CMAQ) perform very poorly in simulating FT chemistry, with Uniform ozone around 70 ppb throughout the FT in summer, while ozonesonde data show much higher levels of ozone and much spatial-temporal structure. Such models completely neglect lightning-NOx (LNOx) emissions (the most significant source of NOx in the FT), and also contain large uncertainties in the specifications of intercontinental transport, stratosphere-troposphere exchange (STE) and PBLFT exchange (PFTE). Global air chemistry models include LNOx, but in very crude fashion, with the frequency and distribution of lightning being based on modeled cloud parameters (hence large uncertainty), lightning energetics being assumed to be constant for all flashes (literature value, while in reality there is at least a two-orders of magnitude variability from flash-to-flash), and the production of NOx in the surrounding heated air, per Joule of heating, being assumed to be constant also (literature value, while in fact it is a non-linear function of the dissipated heat and local air density, p). This situation is commonly blamed on paucity of pertinent observational data, but for the USA, there is now a wealth of surface- and satellite-based data of lightning available to permit much improved observation-based estimation of LNOx emissions. In the FT, such NOx has a long residence time, and also the ozone production efficiency from NOx there is considerably higher than in the PBL. It is, therefore, of critical importance in FT chemistry. This paper will describe the approach and data products of an ongoing NSSTC project aimed at a much-improved quantification of not only LNOx production on the scale of continental USA based on local and regional lightning observations, but also of intercontinental transport, STE and PFTE, all in upgraded simulations of tropospheric transport and chemistry. In our approach for LNOx, (a) we utilize continuous observed lightning information from the NLDN ground network and from satellite imagers (OTD and LIS) to quantify lightning frequency and distribution at the spatial-temporal scales of models such as CMAQ; (b) we develop new methodologies to quantify flash-specific lightning energy dissipation as heat (epsilon) using data from the research-grade lightning measurement facility at NASA-KSC, and to parameterize epsilon based on regional lightning monitoring data (ground- and satellite-based); and, (c) we develop a new parameterization of NOx production as a function of epsilon and rho. Based on such observation-based information, we are working to develop a gridded, episodic LNOx emissions inventory for the USA for use in models like CMAQ. We are also developing approaches for global(MOZART)- regional(CMAQ) chemistry coupling to improve intercontinental transport and STE. Finally, we are developing new methodologies for assimilation of satellite-observed (GOES) clouds into meteorological modeling (MM5), to improve PFTE and to optimize co-location of cloud convection and observed lightning. We will incorporate these improvements in CMAQ simulations over the USA to better understand FT processes and chemistry, and its impact on ground-level ozone.

Gillani, Noor↗
DOE OSTI2026

Microscopic theory of angular momentum distributions across the full range of fission fragments

Modern nuclear theory provides qualitative insights into the fundamental mechanisms of nuclear fission and is increasingly capable of making reliable quantitative predictions. Most quantities of interest pertain to the primary fission fragments, whose subsequent decay is typically modeled using statistical reaction models. Consequently, a key objective of fission theory is to inform these models by predicting the initial conditions of the primary fragments. Here, in this work, we employ a framework that combines joint angular momentum and particle number projection with time-dependent configuration mixing to calculate the angular momentum distributions of primary fragments. Focusing on the benchmark cases of neutron-induced fission of 235 U and 239 Pu , we predict—for the first time—microscopic angular momentum distributions for all fragments observed in experiments. Our results reveal a pronounced sawtooth pattern in the average angular momentum as a function of fragment mass, consistent with recent measurements. Additionally, we observe substantial variations in angular momentum distributions along isobaric chains, indicating that commonly used empirical formulas lack sufficient accuracy. We also quantify a strong correlation between the angular momentum and the deformation of the fragments at scission, and a weak correlation in the magnitude of the angular momentum between fragment partners. The generated data will enable estimation of the impact of microscopic distributions on fission spectra, paving the way toward fission modeling based on microscopic inputs.

Physics - Nuclear physics and radiation physics↗
NASA NTRS2011

Global Intercomparison of 12 Land Surface Heat Flux Estimates

A global intercomparison of 12 monthly mean land surface heat flux products for the period 1993-1995 is presented. The intercomparison includes some of the first emerging global satellite-based products (developed at Paris Observatory, Max Planck Institute for Biogeochemistry, University of California Berkeley, University of Maryland, and Princeton University) and examples of fluxes produced by reanalyses (ERA-Interim, MERRA, NCEP-DOE) and off-line land surface models (GSWP-2, GLDAS CLM/ Mosaic/Noah). An intercomparison of the global latent heat flux (Q(sub le)) annual means shows a spread of approx 20 W/sq m (all-product global average of approx 45 W/sq m). A similar spread is observed for the sensible (Q(sub h)) and net radiative (R(sub n)) fluxes. In general, the products correlate well with each other, helped by the large seasonal variability and common forcing data for some of the products. Expected spatial distributions related to the major climatic regimes and geographical features are reproduced by all products. Nevertheless, large Q(sub le)and Q(sub h) absolute differences are also observed. The fluxes were spatially averaged for 10 vegetation classes. The larger Q(sub le) differences were observed for the rain forest but, when normalized by mean fluxes, the differences were comparable to other classes. In general, the correlations between Q(sub le) and R(sub n) were higher for the satellite-based products compared with the reanalyses and off-line models. The fluxes were also averaged for 10 selected basins. The seasonality was generally well captured by all products, but large differences in the flux partitioning were observed for some products and basins.

Jimenez, C.↗
DOE OSTI2020

Tracking Seasonal and Interannual Variability in Photosynthetic Downregulation in Response to Water Stress at a Temperate Deciduous Forest

Abstract The understanding and modeling of photosynthetic dynamics affected by climate variability can be highly uncertain. In this paper, we examined a well‐characterized eddy covariance site in a drought‐prone temperate deciduous broadleaf forest combining tower measurements and satellite observations. We find that an increase in spring temperature usually leads to enhanced spring gross primary production (GPP), but a GPP reduction in late growing season due to water limitation. We evaluated how well a coupled fluorescence‐photosynthesis model (SCOPE) and satellite data sets track the interannual and seasonal variations of tower GPP from 2007 to 2016. In SCOPE, a simple stress factor scaling of Vcmax as a linear function of observed predawn leaf water potential ( ψ pd ) shows a good agreement between modeled and measured interannual variations in both GPP and solar‐induced chlorophyll fluorescence (SIF) from the Global Ozone Monitoring Experiment‐2 (GOME‐2). The modeled and satellite‐observed changes in SIF yield are ~30% smaller than corresponding changes in light use efficiency (LUE) under severe stress, for which a common linear SIF to GPP scaling would underestimate the stress reduction in GPP. Overall, GOME‐2 SIF tracks interannual tower GPP variations better than satellite vegetations indices (VIs) representing canopy “greenness.” However, it is still challenging to attribute observed SIF variations unequivocally to greenness or physiological changes due to large GOME‐2 footprint. Higher‐resolution SIF data sets (e.g., TROPOMI) already show the potential to well capture the downregulation of late‐season GPP and could pave the way to better disentangle canopy structural and physiological changes in the future.

He, Liyin↗
DOE OSTI2022

What you get is not always what you see—pitfalls in solar array assessment using overhead imagery

Effective integration planning for small, distributed solar photovoltaic (PV) arrays into electric power grids requires access to high quality data: the location and power capacity of individual solar PV arrays. Unfortunately, national databases of small-scale solar PV do not exist; those that do are limited in their spatial resolution, typically aggregated up to state or national levels. While several promising approaches for solar PV detection have been published, strategies for evaluating the performance of these models are often highly heterogeneous from study to study. The resulting comparison of these methods for practical applications for energy assessments becomes challenging and may imply that the reported performance evaluations overly optimistic. The heterogeneity comes in many forms, each of which we explore in this work: the degree of diversity of the locations and sensors (e.g. different satellites, aerial photography) from which the training and validation data originate, the validation of ground truth (manual annotation of imagery vs known solar PV locations), the level of spatial aggregation (e.g. array-level vs regional estimates), and inconsistencies in the training and validation datasets (e.g. different datasets are used for each study and those data are not always made accessible). For each, we discuss emerging practices from the literature to address them or suggest directions of future research. As part of our investigation, we evaluate solar PV identification performance in two large regions: the entire state of Connecticut and the city of San Diego, CA. In Connecticut, we also use 33,114 known parcel-level solar PV installations from Berkeley Lab’s Tracking the Sun dataset to evaluate parcel-level performance and evaluate capacity estimates using 169 municipalities. We also make our code (which we call SolarMapper), pre-trained models, training data, and predictions publicly available and provide a web portal for interactively inspecting each prediction that was made. Here our findings suggest that traditional performance evaluation of the automated identification of solar PV from satellite imagery may be optimistic due to common limitations in the validation process. The takeaways from this work are intended to inform and catalyze the large-scale practical application of automated solar PV assessment techniques by energy researchers and professionals.

14 SOLAR ENERGY↗
DOE OSTI2023

Kinetic model-based group contribution method for derived cetane number prediction of oxygenated fuel components and blends

A four-step autoignition model was used to derive an expression for the ignition delay (ID) measured in ignition quality testers (IQT) with the derived cetane number (DCN) determined from this ID using the ASTM D6890 standard correlation. The model predicts DCN for individual compounds and blends as a function of each compound's global initiation and net chain branching rate constants. Expressions for these values were determined assuming they could be related to the functional groups present in each compound. Measurement data for 125 compounds and 94 binary and ternary blends (including oxygenates: alcohols, aldehydes, esters, ethers, and ketones), from literature and from measurements performed in our lab, were used to obtain the dependence of the measured ignition delay on each functional group. The new kinetic model-based group contribution method was able to predict the ignition delay of both pure compounds and blends, with an average DCN error of 4.4 (19%) and 2.8 (11%), respectively. Here, the blend model was also used to develop an ID mixing rule by incorporating existing IQT ignition delay data for each compound. Use of the mixing rule gave an average DCN error of 3.6 (16%) for blends. Both the blend model and mixing rule were found to be superior compared to standard linear by volume fraction or mole fraction mixing rules commonly used to estimate the DCN of mixtures.

42 ENGINEERING↗
DOE OSTI2023

Early Alarm: Robust Event Analysis for Power Systems using 1-D Fully Convolutional Network

This work presents a novel deep learning model for early, accurate, and robust detection, recognition, and temporal localization of multi-type events in large-scale power systems. The proposed method develops a unified 1-D fully convolutional network (FCN) model that takes time series of raw frequency signals measured from a power system as input, extracts distinguishing features, and predicts at every temporal point in the time series if an event is happening and what the type of the event is. Compared to existing methods, the proposed model eliminates the necessity for hand-crafted feature extraction or complicated data pre-processing, can flexibly handle input signals of arbitrary length, and precisely infer the event occurrence time. Most importantly, the model is intentionally trained with incomplete patterns, such that it is more robust to partial features of an event which is common in real-world online recognition, resulting in early alarm for power system failures. Extensive experimental results demonstrate that the proposed method achieves superior performance to the state-of-the-art, and also shows strong robustness to noise and system oscillations.

Li, Chengcheng↗
DOE OSTI2021

Correlating Time-Resolved Pressure Measurements With Rim Sealing Effectiveness for Real-Time Turbine Health Monitoring

Purge flow is bled from the upstream compressor and supplied to the under-platform region to prevent hot main gas path ingress that damages vulnerable under-platform hardware components. A majority of turbine rim seal research has sought to identify methods of improving sealing technologies and understanding the physical mechanisms that drive ingress. While these studies directly support the design and analysis of advanced rim seal geometries and purge flow systems, the studies are limited in their applicability to real-time monitoring required for condition-based operation and maintenance. As operational hours increase for in-service engines, this lack of rim seal performance feedback results in progressive degradation of sealing effectiveness, thereby leading to reduced hardware life. To address this need for rim seal performance monitoring, the present study utilizes measurements from a one-stage turbine research facility operating with true-scale engine hardware at engine-relevant conditions. Time-resolved pressure measurements collected from the rim seal region are regressed with sealing effectiveness through the use of common machine learning techniques to provide real-time feedback of sealing effectiveness. Two modelling approaches are presented that use a single sensor to predict sealing effectiveness accurately over a range of two turbine operating conditions. Results show that an initial purely data-driven model can be further improved using domain knowledge of relevant turbine operations, which yields sealing effectiveness predictions within three percent of measured values.

Compressors↗
DOE OSTI2026

SetGo: Metadata Readiness for Scientific AI Datasets

Scientific datasets intended for AI use require both computational readiness for model training and metadata readiness for discovery, sharing, and reuse. The Readiness Engine for Data Integration (REDI) addresses computational readiness, but no corresponding tool evaluates whether a dataset’s metadata are sufficiently complete, governed, and standards-compliant for publication and agent-based consumption. Existing FAIR assessors operate only on published repository records, and no single system covers FAIR compliance, licensing, provenance, governance, reproducibility, and catalog readiness together. We present SetGo, an open-source Python toolkit that assesses and repairs metadata readiness across these six dimensions before a dataset is published or archived. Applied to four scientific corpora, SetGo surfaces deficiencies that general-purpose tools do not detect: ERA5 climate metadata scores 4% on ACDD 1.3 compliance; materials datasets fail OPTIMADE species-definition requirements; and PDB-derived proteomics data carries licensing terms incompatible with standard SPDX identifiers. Guided enrichment raises overall FAIR scores from 52–57% to 81–91%, and a single setgo publish command pushes to Hugging Face Hub, CKAN, or OpenMetadata with ML Commons Croissant 1.0 metadata sidecars. To support interactive and automated workflows, SetGo integrates with coding agents powered by large language models (LLMs) through a /setgo skill that enables natural-language execution of the full assess–enrich–publish loop, with user involvement limited to supplying missing metadata values.

Wilkinson, Sean [ORNL] (ORCID:0000000214437479)↗
DOE OSTI2023

Automatic Point Cloud Building Envelope Segmentation (AutoCuBES)

The Auto-CuBES algorithm is based on unsupervised machine learning that automatically labels 3D point cloud data and reduces the time spent in manual segmentation. The algorithm can process high-resolution point clouds and generate a wire-frame building envelope model with a small set of calibration parameters. The algorithm inputs a 3D point cloud generated by commonly available surveying equipment and outputs a wire-frame model of the building envelope. Unsupervised machine learning methods were used to identify facades, windows, and doors while minimizing the number of calibration parameters.

Puente, BryanMaldonado↗
DOE OSTI2020

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

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

FAIR Surrogate Benchmarks Supporting AI and Simulation Research (Final Report)

Computational Science is being revolutionized by integrating AI and simulation and, in particular, by deep learning surrogate models that can replace all or part of traditional large‐scale HPC computations. Such surrogates can achieve remarkable performance improvements, as much as several orders of magnitude, and save both compute time and energy. The Surrogate Benchmark Initiative (SBI) project creates a community repository and FAIR (Findable, Accessible, Interoperable, and Reusable) data ecosystem for HPC application surrogate benchmarks. The SBI team comes from Argonne National Laboratory (ANL), Indiana University (IU), Rutgers University, the University of Tennessee, Knoxville (UTK), and the University of Virginia (UVA). SBI repositories include data, code, and all relevant collateral artifacts that the science and engineering community need to use and reuse these data sets and surrogates. SBI repositories generate active research from both the participants in SBI and the broad community of AI and domain scientists. This project develops surrogates that use several different neural nets to learn and quickly infer the results of simulations and data systems and captures them as surrogate benchmarks with a rich set of metadata covering: Data; Model; Metrics specification; Machine specification; and Science, Speed, and Power Results. We research FAIR metadata for these benchmarks. We develop application surrogate examples as benchmarks across many fields (ANL, UTK, IU, UVA). We also study non-Surrogate benchmarks that have many common features and similar issues as regards FAIRness. We work with MLCommons (UVA, UTK), which is a major machine learning benchmarking activity where we get metadata ontologies, software, and benchmarks, Benchmarks have datasets, models, and metadata and they need a technical framework developed by UTK and Rutgers and deployed by UVA. We study features of Surrogates including performance, training set size, and uncertainty quantification (Rutgers, UVA and IU).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Confined turbulent fluid-particle flow modeling using multiple-realization particle trajectory schemes

A multiple-realization particle trajectory scheme has been developed and applied to the numerical prediction of confined turbulent fluid-particle flows. The example flows investigated include the vertical pipe upflow experimental data of Tsuji et al. and the experimental data of Leavitt for a coaxial jet flow, comprising a particle-laden central jet and a clean annular jet, into a large recirculation chamber. The results obtained from the numerical scheme agree well with the experimental data, lending confidence to the modeling approach. The multiple-realization particle trajectory turbulent flow modeling scheme is believed to be a more elegant and accurate approach to the extension of single-particle hydrodynamics to dilute multi-particle systems than the more commonly employed two-fluid modeling approach. It is also better able to incorporate additional force items such as lift, virtual mass and Bassett history terms directly into the particle equation of motion as appropriate. This makes it a suitable candidate for particle migration studies and an extension to situations involving liquid particulate phases with possible propulsion applications, such as in spray combustion, follows naturally.

Adeniji-Fashola, A. A.↗
NASA NTRS1994

General purpose architecture for intelligent computer-aided training

An intelligent computer-aided training system having a general modular architecture is provided for use in a wide variety of training tasks and environments. It is comprised of a user interface which permits the trainee to access the same information available in the task environment and serves as a means for the trainee to assert actions to the system; a domain expert which is sufficiently intelligent to use the same information available to the trainee and carry out the task assigned to the trainee; a training session manager for examining the assertions made by the domain expert and by the trainee for evaluating such trainee assertions and providing guidance to the trainee which are appropriate to his acquired skill level; a trainee model which contains a history of the trainee interactions with the system together with summary evaluative data; an intelligent training scenario generator for designing increasingly complex training exercises based on the current skill level contained in the trainee model and on any weaknesses or deficiencies that the trainee has exhibited in previous interactions; and a blackboard that provides a common fact base for communication between the other components of the system. Preferably, the domain expert contains a list of 'mal-rules' which typifies errors that are usually made by novice trainees. Also preferably, the training session manager comprises an intelligent error detection means and an intelligent error handling means. The present invention utilizes a rule-based language having a control structure whereby a specific message passing protocol is utilized with respect to tasks which are procedural or step-by-step in structure. The rules can be activated by the trainee in any order to reach the solution by any valid or correct path.

Loftin, R. Bowen↗
NASA NTRSDate not supplied

Icing Formation Strains on Airfoil Skin

Aircraft icing has long presented a safety hazard to the aerospace industry. One problem inherent to aircraft icing is predicting the adhesive properties of the ice-aircraft interface; values commonly reported span multiple orders of magnitude, making shedding prediction and surface evaluation difficult. An effort has been undertaken at the NASA Glenn Research Center to understand the cause of the scatter in the data. Geometry dependencies, variation in testing procedures, and residual stresses in ice samples have previously been identified as significant contributors to the problem. Residual stresses were singled out for this study, which investigated the strains at the leading edge of a Common Research Model (CRM) wind tunnel model in the NASA Glenn Icing Research Tunnel. Sprays were conducted at a particle size median volumetric diameter (MVD) of 20 µm, with liquid water contents (LWC) spanning from 0.4 to 3.0 g/m3. Total temperatures were varied between 0 and -25 °C, and test section velocity between 51 and 118 m/s. The leading edge was instrumented with custom strain gages, resistance temperature detectors (RTDs), optical fibers, and a pressure mat. Results indicate Δ strains, which were the difference between recorded strains and those predicted from isotropic thermal expansion of aluminum, below -400 µε with consistent patterns of strain across a wide variety of tunnel conditions. Isotropic, linear elastic finite element analysis (FEA) was conducted by varying the temperature of simplified airfoil and ice geometries, producing similar strains to those observed in testing. This allowed induced stresses to be predicted. Predicted stresses were above the expected tensile strength for ice, and suggest a source for previously observed cracking in tunnel ice.

Ice Adhesion↗
NASA NTRS2022

Rock Abundance on the Lunar Mare on Surfaces of Different Age: Implications for Regolith Evolution and Thickness

The growth of lunar regolith over time affects surface rock abundance, because larger, less frequent impacts are needed to penetrate thicker regolith developed on older surfaces and excavate rocks. On younger surfaces with thinner regolith, smaller, more frequent impacts are sufficient to excavate rocks. We quantify the correlation between observed rock abundances and age on the lunar surface by comparing Diviner rock abundance data to the surface ages of inter-crater parts on the maria. Our observations show the expected negative correlation between age and rock abundance. The commonality of non-zero rock abundance values on ancient surfaces, combined with a simple Monte Carlo model of the rock excavation process, suggest that rocks re-excavated from the regolith volume contribute to the presently observed rock population on the lunar surface. The half-life of meter-scale surface rocks most consistent with our observations is 80±20 Myr.

Sashank Vanga↗
DOE OSTI2026

Constraining Hamiltonians from chiral effective field theory with neutron-star data

Multi-messenger observations of neutron stars (NSs) and their mergers have placed strong constraints on the dense-matter equation of state (EOS). The EOS, in turn, depends on microscopic nuclear interactions that are described by nuclear Hamiltonians. These Hamiltonians are commonly derived within chiral effective field theory (EFT). Ideally, multi-messenger observations of NSs could be used to directly inform our understanding of EFT interactions, but such a direct inference necessitates millions of model evaluations. This is computationally prohibitive because each evaluation requires us to calculate the EOS from a Hamiltonian by solving the quantum many-body problem with methods such as auxiliary-field diffusion Monte Carlo (AFDMC), which provides very accurate and precise solutions but at a significant computational cost. Additionally, we need to solve the stellar structure equations for each EOS which further slows down each model evaluation by a few seconds. In this work, we combine emulators for AFDMC calculations of neutron matter, built using parametric matrix models, and for the stellar structure equations, built using multilayer perceptron neural networks, with the PyCBC data-analysis framework to enable a direct inference of coupling constants in an EFT Hamiltonian using multi-messenger observations of NSs. We find that astrophysical data can provide informative constraints on two-nucleon couplings despite the high densities probed in NS interiors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗