Search NASA⌕ Search

SEARCH · Search NASA

Results for “array modeling”

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

21 cm Power Spectrum Analysis of North Celestial Pole Observations with the Tianlai Dish Pathfinder Array

The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization Universe as a means of measuring large-scale cosmic structure. Using nine nights of observations targeting the North Celestial Pole field, totaling approximately 107 hr of integration time, we analyze data in the frequency range 700–800 MHz (corresponding to redshift z ∼ 0.9). We do the data format conversion, radio frequency interference flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition. The spherically averaged power spectrum Δ 2 (k) is obtained. Furthermore, this work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.

cosmology: large-scale structure of universe↗
DOE OSTI2026

Emulating 2D Materials with Magnons

Spin waves (magnons) in two-dimensional (2D) materials have received increasing interest due to their unique states and potential for tunability. However, many interesting features of these systems, including Dirac points and topological states, occur at high frequencies, where experimental probes are limited. Here, we study a crystal formed by patterning a hexagonal array of holes in a perpendicularly magnetized thin film. Through simulation, we find that the magnonic band structure imitates that of graphene, but additionally has some kagomelike character and includes a few flat bands. Surprisingly, its nature can be understood using a nine-band tight-binding Hamiltonian. This clear analogy to 2D materials enables band-gap engineering in 2D, topological magnons along 1D phase boundaries, and spectrally isolated modes at 0D point defects. Interestingly, the 1D phase boundaries allow access to the valley degree of freedom through a magnonic analog of the quantum valley Hall insulator. These approaches can be extended to other magnonic systems, but are potentially more general due to the simplicity of the model, which resembles existing results from electron, phonon, photon, and cold-atom systems. This finding brings the physics of spin waves in 2D materials to more experimentally accessible scales, augments it, and outlines a few principles for controlling magnonic states.

Ferrimagnets↗
DOE OSTI2025

Mapping Spiking Neural Networks to Heterogeneous Crossbar Architectures using Integer Linear Programming

Advances in novel hardware devices and architectures allow Spiking Neural Network (SNN) evaluation using ultra-low power, mixed-signal, memristor crossbar arrays. As individual network sizes quickly scale beyond the dimensional capabilities of single crossbars, networks must be mapped onto multiple crossbars. Crossbar sizes within modern Memristor Crossbar Architectures (MCAs) are determined predominately not by device technology but by network topology; more, smaller crossbars consume less area thanks to the high structural sparsity found in larger, brain-inspired SNNs. Motivated by continuing increases in SNN sparsity due to improvements in training methods, we propose utilizing heterogeneous crossbar sizes to further reduce area consumption. This approach was previously unachievable as prior compiler studies only explored solutions targeting homogeneous MCAs. Our work improves on the state-of-the-art by providing Integer Linear Programming (ILP) formulations supporting arbitrarily heterogeneous architectures. By modeling axonal interactions between neurons, our methods produce better mappings while removing inhibitive a priori knowledge requirements. We first show a 16.7-27.6% reduction in area consumption for square-crossbar homogeneous architectures. Then, we demonstrate 66.9-72.7% further reduction when using a reasonable configuration of heterogeneous crossbar dimensions. Next, we present a new optimization formulation capable of minimizing the number of inter-crossbar routes. When applied to solutions already near-optimal in area, an 11.9-26.4% routing reduction is observed without impacting area consumption. Finally, we present a profile-guided optimization capable of minimizing the number of runtime spikes between crossbars. Compared to the best-area-then-route optimized solutions, we observe a further 0.5-14.8% inter-crossbar spike reduction while requiring 1–3 orders of magnitude less solver time.

Pohl, Devin [ORNL] (ORCID:0009000040149027)↗
DOE OSTI2026

Locating the missing large-scale emission in the jet of M87* with short EHT baselines

In very-long baseline interferometric arrays, nearly co-located stations probe the largest scales and typically cannot resolve the observed source. In the absence of a large-scale structure, closure phases constructed with these stations are zero and, since they are independent of station-based errors, they can be used to probe data issues. Here, we show how these trivial closure phases become nonzero with a brightness distribution on smaller scales than their short baseline would suggest. When applied to sources that are made up of a bright compact and large-scale diffuse component, the trivial closure phases directly measure the centroid relative to the compact source and higher-order image moments. We present a technique to measure these image moments with minimal model assumptions and validate it on synthetic Event Horizon Telescope (EHT) data. We then apply this technique to 2017 and 2018 EHT observations of M87* and find a weak preference for extended emission in the direction of the large-scale jet. We also apply it to 2021 EHT data and measure the source centroid about 1 mas northwest of the compact ring, which is consistent with the jet observed at lower frequencies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Gravitational wave spectrum of chain inflation

Chain inflation is an alternative to slow-roll inflation in which the inflaton tunnels along a large number of consecutive minima in its potential. In this work we perform the first comprehensive calculation of the gravitational wave (GW) spectrum of chain inflation. In contrast to slow-roll inflation the latter does not stem from quantum fluctuations of the gravitational field during inflation, but rather from the bubble collisions during the first-order phase transitions associated with vacuum tunneling. Our calculation is performed within an effective theory of chain inflation which builds on an expansion of the tunneling rate capturing most of the available model space. The effective theory can be seen as chain inflation’s analog of the slow-roll expansion in rolling models of inflation. The near scale-invariance of the scalar power spectrum translates to a quasiperiodic shape of the inflaton potential in chain inflation, with the tunneling rate changing very slowly during the e-folds leading to cosmic microwave background observables. We show that chain inflation produces a very characteristic double-peak GW spectrum: a faint high-frequency peak associated with the gravitational radiation emitted during inflation, and a strong low-frequency peak associated with the graceful exit from chain inflation (marking the transition to the radiation-dominated epoch). There exist very exciting prospects to test the gravitational wave signal from chain inflation at the aLIGO-aVIRGO-KAGRA network, at LISA and /or at pulsar timing array experiments. A particularly intriguing possibility we point out is that chain inflation could be the source of the stochastic gravitational wave background recently detected by NANOGrav, PPTA, EPTA, and CPTA. We also show that the gravitational wave signal of chain inflation is often accompanied by running/ higher running of the scalar spectral index to be tested at future cosmic microwave background experiments. Published by the American Physical Society 2024

Freese, Katherine↗
DOE OSTI2026

Exploring the keV-scale physics potential of CUORE

We present the analysis techniques developed to explore the keV-scale energy region of the Cryogenic Underground Observatory for Rare Events (CUORE) experiment, based on more than 2 metric ton yr of data collected over five years. By prioritizing a stricter selection over a larger exposure, we are able to optimize data selection for thresholds at 10 keV and 3 keV with 691 kg yr and 11 kg yr of data, respectively. We study how the performance varies among the 988-detector array with different detector characteristics and data-taking conditions. We achieve an average baseline resolution of 2.54 ±0.14 keV FWHM and 1.18 ±0.02 keV FWHM for the data selection at 10 keV and 3 keV, respectively. The analysis methods employed reduce the overall background by about an order of magnitude, reaching 2.06±0.05 counts/(keV kg days) and 16±2 counts/(keV kg days) at the thresholds of 10 keV and 3 keV. We evaluate for the first time the near-threshold reconstruction efficiencies of the CUORE experiment, and find these to be 50 ±2% and 26 ±4% at 10 keV and 3 keV, respectively. This analysis provides crucial insights into rare decay studies, new physics searches, and keV-scale background modeling with CUORE. We demonstrate that ton-scale cryogenic calorimeters can operate across a wide energy range, from keV to MeV, establishing their scalability as versatile detectors for rare event and dark matter physics. These findings also inform the optimization of future large mass cryogenic calorimeters to enhance the sensitivity to low-energy phenomena.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Quantum annealing-aided design of an ultrathin-metamaterial optical diode

Abstract Thin-film optical diodes are important elements for miniaturizing photonic systems. However, the design of optical diodes relies on empirical and heuristic approaches. This poses a significant challenge for identifying optimal structural models of optical diodes at given wavelengths. Here, we leverage a quantum annealing-enhanced active learning scheme to automatically identify optimal designs of 130 nm-thick optical diodes. An optical diode is a stratified volume diffractive film discretized into rectangular pixels, where each pixel is assigned to either a metal or dielectric. The proposed scheme identifies the optimal material states of each pixel, maximizing the quality of optical isolation at given wavelengths. Consequently, we successfully identify optimal structures at three specific wavelengths (600, 800, and 1000 nm). In the best-case scenario, when the forward transmissivity is 85%, the backward transmissivity is 0.1%. Electromagnetic field profiles reveal that the designed diode strongly supports surface plasmons coupled across counterintuitive metal–dielectric pixel arrays. Thereby, it yields the transmission of first-order diffracted light with a high amplitude. In contrast, backward transmission has decoupled surface plasmons that redirect Poynting vectors back to the incident medium, resulting in near attenuation of its transmission. In addition, we experimentally verify the optical isolation function of the optical diode.

36 MATERIALS SCIENCE↗
DOE OSTI2024

FiberFlex: Real-time FPGA-based Intelligent and Distributed Fiber Sensor System for Pedestrian Recognition

In recent years, security monitoring of public places and critical infrastructure has heavily relied on the widespread use of cameras, raising concerns about personal privacy violations. To balance the need for effective security monitoring with the protection of personal privacy, we explore the potential of optical fiber sensors for this application. This article proposes FiberFlex, an intelligent and distributed fiber sensor system. Ultizing Field Programmable Gate Arrays (FPGA) high-level synthesis (HLS) acceleration, FiberFlex offers real-time pedestrian detection by co-designing the entire pipeline of optical signal acquisition, processing, and recognition networks based on the principles of optical fiber sensing. As a promising alternative to traditional camera-based monitoring systems, FiberFlex achieves pedestrian detection by analyzing the vibration patterns caused by pedestrian footsteps, enabling security monitoring while preserving individual privacy. FiberFlex comprises three modules: First , fiber-optic sensing system: A fiber-optic distributed acoustic sensing (DAS) system is built and used to measure the ground vibration waves generated by people walking. Second , algorithms: We first collect the training data by measuring the ground vibration waves, label the data, and use the data to train the neural network models to perform pedestrian recognition. Third , hardware accelerators: We use HLS tools to design hardware modules on FPGA for data collection and pre-processing and integrate them with the downstream neural network accelerators to perform in-line real-time pedestrian detection. The final detection results are sent back from FPGA to the host CPU. We implement our system FiberFlex with the in-house built DAS system and AMD/Xilinx Kintex7 FPGA KC705 board and verify the whole system using the real-world collected data. We conduct recognition tests on five test subjects of varying ages, heights, and weights in a fixed sensing area. Each subject experienced 20 real-time recognition tests using their daily walking habits, and the subjects were given adequate rest between tests. After 100 tests on five test subjects, the overall real-time recognition accuracy exceeded \(88.0\%\) . The whole system uses 55 W of power, 33 W in the optical DAS system and 22 W in the FPGA. Relying on its end-to-end interdisciplinary design, FiberFlex seamlessly combines fiber-optic sensors with FPGA accelerators to enable low-power real-time security monitoring without compromising privacy, making it a valuable addition to the existing security monitoring network. According to FiberFlex, more valuable research can be conducted in the future, such as fall monitoring for the elderly, migration of identification networks between different application scenarios, and improvement of anti-interference performance in more complex environments. In future perception networks, where the “eyes” are not feasible, let’s use fiber optic touch instead.

Distributed↗
DOE OSTI2025

A Parameter-masked Mock Data Challenge for Beyond-two-point Galaxy Clustering Statistics

The past few years have seen the emergence of a wide array of novel techniques for analyzing high-precision data from upcoming galaxy surveys, which aim to extend the statistical analysis of galaxy clustering data beyond the linear regime and the canonical two-point (2pt) statistics. We test and benchmark some of these new techniques in a community data challenge named “Beyond-2pt,” initiated during the Aspen 2022 Summer Program “Large-Scale Structure Cosmology beyond 2-Point Statistics,” whose first round of results we present here. The challenge data set consists of high-precision mock galaxy catalogs for clustering in real space, in redshift space, and on a light cone. Participants in the challenge have developed end-to-end pipelines to analyze mock catalogs and extract unknown (“masked”) cosmological parameters of the underlying ΛCDM models with their methods. The methods represented are density-split clustering, nearest neighbor statistics, BACCO power spectrum emulator, void statistics, LEFTfield field-level inference using effective field theory (EFT), and joint power spectrum and bispectrum analyses using both EFT and simulation-based inference. In this work, we review the results of the challenge, focusing on problems solved, lessons learned, and future research needed to perfect the emerging beyond-2pt approaches. The unbiased parameter recovery demonstrated in this challenge by multiple statistics and the associated modeling and inference frameworks supports the credibility of cosmology constraints from these methods. The challenge data set is publicly available, and we welcome future submissions from methods that are not yet represented.

Krause, Elisabeth [Univ. of Arizona, Tucson, AZ (U↗
DOE OSTI2026

Solvent selection for a biomass-to-bioproduct pipeline through integrated reductive catalytic fractionation and microbial funneling

The growing significance of lignin-first biorefineries, which focus on upgrading the aromatics resulting from lignin depolymerization, presents opportunities for bioproduct synthesis using microbial strains capable of funneling a diverse array of phenolics into a single commodity chemical. In this study, we evaluated a biomass-to-bioproduct pipeline involving the reductive catalytic fractionation (RCF) of poplar biomass followed by biological funneling with a Novosphingobium aromaticivorans strain that produces 2-pyrone-4,6-dicarboxylic acid (PDC), a potential bioplastic precursor. Considering the impact of solvent on RCF reactor operating pressure, and the potential inhibitory effects of solvent on downstream microbial funneling, we performed an analysis of six pure solvents, namely methanol, ethanol, isopropanol, isobutanol, 1,4-dioxane and ethylene glycol, and different variations of their aqueous mixtures comprising 5 to 50 vol% water. For each pure solvent and solvent/water system, we measured phenolic monomer yields in the RCF process and PDC yields from the phenolic monomers. We then developed correlation models that relate phenolic monomer yields from RCF-derived samples to Hansen solubility parameters to determine solvent descriptors that contribute to high yields. Furthermore, we developed an integrated biorefinery system to estimate the minimum selling price (MSP) of PDC and the associated carbon footprint to identify solvent systems with better costs and sustainability metrics. These analyses resulted in the 50 vol% methanol/water system being identified as optimal because it reduces RCF reactor pressure and is compatible with microbial funneling with N. aromaticivorans. This solvent system produced 63 g PDC per kg biomass (264 g PDC per kg lignin) from 85 g phenolic monomers per kg biomass at a reduced reactor pressure of 48 bar (reduced by 26% compared to our previous poplar-to-PDC pipeline). The MSP for this system is $\$$13.98 per kg of purified PDC (carbon footprint of 1.47 kg CO 2 e per kg), which is about 24% lower than a previously described poplar-to-PDC pipeline and 46% lower than a lignin-to-PDC pipeline that used pure methanol as the solvent. The results from this study illustrate improvements that can be made in lignocellulosic biorefineries that are compatible with the hybrid chemical and biological processes needed to gain value from lignin.

Sripada, Sarada [Great Lakes Bioenergy Research Ce↗
DOE OSTI2024

Probing the Impact of Radio-mode Feedback on the Properties of the Cool Circumgalactic Medium

We explore the influence of radio-mode feedback on the properties of the cool circumgalactic medium (CGM). To this end, we assemble a statistical sample of approximately 30,000 radio galaxies with background quasars by combining optical spectroscopic measurements of luminous red galaxies and quasars from the year 1 data set of the Dark Energy Spectroscopic Instrument and radio sources from the LOw-Frequency ARray Two-metre Sky Survey (LoTSS) DR2 catalog and the Very Large Array Sky Survey (VLASS) quick-look catalog. Galaxies with similar optical properties but with no radio counterparts in LoTSS and VLASS are selected as the control group. We measure the cool CGM properties of radio galaxies and their control samples traced by Mg ii absorption lines, including covering fraction, rest equivalent width, and gas kinematics. Our results show no significant difference in the properties of gas around radio galaxies and their control sample, indicating that the operating radio-mode feedback of massive galaxies does not produce detectable effects on the properties of the cool CGM. Finally, we show that the CGM of radio galaxies contains a nonnegligible amount of cool gas with approximately 10 10 M ⊙ . This abundance can place a stringent constraint on the radio-mode feedback models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Even Higher-Level Synthesis: An Exploration of AI Hardware Accelerators using HLS4ML

With the rise of artificial intelligence, the popularization of deep learning, and a constantly evolving industry, the demand for flexible and efficient tools has never been greater. As algorithms grow more complex, their runtime and energy consumption increase exponentially. Customized hardware accelerators, long used for specific mathematical operations, remain essential for managing modern applications' computational and power demands. Hardware accelerators can speed up complex computations by orders of magnitude, but their manual design and verification processes are often challenging and time-consuming. High-Level Synthesis (HLS) provides a solution by transforming high-level algorithm descriptions, typically written in C++ or SystemC, into synthesizable RTL suitable for hardware implementation. This approach reduces development time for RTL engineers while offering flexibility beyond what traditional handwritten RTL can provide. We extended this capability to the machine-learning domain with the open-source framework hls4ml, which allows neural networks trained in Python frameworks like Tensorflow or PyTorch to be synthesized into efficient hardware representations for the traditional FPGA and ASIC flows. This breakthrough addresses the growing need for reduced design turnaround and easy verification of ML hardware accelerators with low latency and power efficiency constraints. During this tutorial, we will demonstrate how Python complements HLS by simplifying the ML design process, bridging the gap between software and hardware development. Attendees will explore how we translate neural networks modeled in Python into fixed-point C++ models suitable for HLS workflows. We will dive into strategies like Value-Range Analysis and Quantization-Aware Training, which optimize these designs for deployment and evaluate their accuracy, power consumption, and energy efficiency. To exemplify these concepts, experts from Fermilab will share their experiences applying this technology to high-energy physics experiments, where real-time, low-latency processing is critical. Over the years, Fermilab engineers have demonstrated how deep neural networks, optimized for hardware using hls4ml, can meet the stringent requirements of trigger systems at the CERN Large Hadron Collider. These systems rely on rapid decision-making to process immense data volumes while retaining only the most relevant events for further analysis. The application of hls4ml has also been extended to innovative technologies like smart pixel arrays. These smart pixels integrate ML inference capabilities directly into sensor devices, enabling localized data processing at the pixel level. This approach drastically reduces the need to transmit raw data to external processing units, significantly decreasing power consumption and latency. By embedding neural networks within the pixel architecture, the smart pixels can identify and prioritize relevant data in real time, providing a highly efficient solution for edge computing in scenarios such as particle detectors and imaging systems. Fermilab's work highlights the potential of hardware-accelerated ML in scenarios where both speed and power efficiency are mission-critical. Through this tutorial, attendees will gain valuable insights into the challenges and solutions of deploying ML in hardware. Understanding how HLS and hls4ml streamline the development of neural network-based hardware accelerators is fundamental for the industry's future. Participants will learn how these technologies are shaping the future of AI and scientific computing.

Di Guglielmo, Giuseppe [Fermilab]↗
DOE OSTI2024

Reference Site Conditions for Floating Wind Arrays in the United States

Floating offshore wind farm design is highly site-specific, requiring detailed information about the specific conditions of a project area for realistic design studies. Unfortunately, publicly available site condition data for potential floating offshore wind project sites in the United States is scarce. To support U.S. offshore wind research, we developed reference site condition datasets, including metocean and seabed information, for four potential floating wind project areas in the U.S.: Humboldt Bay, Morro Bay, the Gulf of Maine, and the Gulf of Mexico. These datasets were compiled using publicly available data. Our metocean analysis, covering wind, waves, and surface currents, utilized measurement data from 2000 to 2020. Sources included the National Renewable Energy Laboratory’s National Offshore Wind Dataset for wind data, National Data Buoy Center buoys for wave data, and the High Frequency Radar Network for surface currents. These data were integrated into hourly time series used to compute extreme return periods up to 500 years, monthly statistics, and joint probability clusters for fatigue analysis. Soil conditions were evaluated using the usSEABED database and bathymetry grids were interpolated from the NCEI Digital Elevation Model Global Mosaic. In addition to providing curated reference site condition datasets for four U.S. areas, our assessment highlights the need for more publicly available metocean and soil condition data.

17 WIND ENERGY↗
DOE OSTI2024

Machine learning-driven descriptions of protein dynamics at solid-liquid interfaces

This chapter has described how ML has enabled quantitative analysis of HS-AFM data to discover the physical phenomena governing protein dynamics and ordering at solid-liquid interfaces. The research detailed in this chapter modeled the rotation models of protein nanorods, the discovery of which would otherwise not be possible. By tracking the trajectories of individual protein rods from frame to frame, it was possible to model Brownian type motion and behaviors and Levy-flight dynamics that had not previously been shown. We also described the application of the Python package AtomAI, which has been developed specifically to analyze and extract physical phenomena, providing exemplar code for training an ensemble of deep neural networks to produce the semantic segmentation of AFM data and functions for encoding and decoding local environments. We last described a combinatorial approach to analyze very noisy data with a densely covered substrate where the emergence of order for the protein liquid crystals could be elucidated. By combining the methods from Case 1 and 2, it was possible to obtain the center of mass and angle for each rod in the images and track the assembly of the rods over time into a 2D liquid crystal array on the surface of mica.

protein dynamics, solid-liquid interfaces, atomic ↗
DOE OSTI2026

Towards Spatial Control of Reaction Selectivity on Photocatalysts using Area Selective Atomic Layer Deposition on Model Dual Site Electrocatalyst Platform

Photocatalytic water splitting is a promising route to low-cost H2; however, this approach is currently limited by solar-to-hydrogen (STH) conversion efficiencies in the sub-10% range. Z-scheme water splitting, in which H2 and O2 evolving particles are operated in separate compartments and electronically coupled by a soluble redox mediator, offers improvements in STH efficiency but introduces high rates of undesired side and back reactions. Nanoscopic oxide coatings (e.g., CrOx, and SiOx) have previously been used to selectively block undesired reactants from reaching active sites; however, a coating encapsulating the entire photocatalyst particle limits activity as the particle can no longer facilitate both half reactions. In contrast to photodeposition, which may produce non-uniform overlayers due to the difficulty of controlling local electrochemical reactivity on the surface, area-selective atomic layer deposition (AS-ALD) can be used to deposit conformal ultrathin oxide coatings while preserving access to non-growth areas of a substrate surface. To develop this technique for Z-scheme photocatalysts, we performed AS-ALD of TiO2 on a dual site planar electrocatalyst based on interdigitated arrays of Pt and Au. Self-assembled monolayers of 1-octadecanethiol (ODT) were used to block growth on the Au array, resulting in a patterned surface in which only the Pt sites were encapsulated with TiO2. These electrocatalysts demonstrated localized reaction selectivity for hydrogen evolution reaction (HER) over the TiO2/Pt sites, while the uncoated Au sites retained activity towards Fe(II)/Fe(III) redox (FeRR). Under independent potential control, these microelectrodes showed that selectively deposited TiO2 coatings can suppress the rate of back reactions on neighboring active sites by an order of magnitude compared to uncoated control samples.

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

Alignments in the orientation of radio jets from AGN and their host galaxies

It is well established that active galactic nuclei (AGNs) play an important role in the evolution of galaxies. These AGNs can be linked to the accretion processes onto massive black holes and past merger events in their host galaxies, which may lead to different alignments of the jets with respect to the host galaxies. This paper presents a study of the position angle (PA) differences between radio and optical images of radio-AGNs based on the second data release (DR2) of the Low Frequency Array (LOFAR) Two-Meter Sky Survey (LoTSS), theKarl G. JanskyVery Large Array Faint Images of the Radio Sky at Twenty-Centimeters Survey (FIRST), the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys, and the Sloan Digital Sky Survey (SDSS). We assessed PA measurement biases in the data and classified the radio-AGNs based on the radio luminosity and infrared colour from the Wide-field Infrared Survey Explorer (WISE). This resulted in the largest sample of radio AGNs with reliable radio and optical PA measurements published to date, with a total of 3682 AGNs. The PA difference (dPA) distributions for the radio-AGN sample show a prominent minor-axis alignment tendency. Based on some simple assumptions, we simulated the projection effect to estimate the intrinsic jet-galaxy alignment. The observed dPA distribution can be well described by a two-component jet-alignment model in which one component is more aligned with the minor axis of the host galaxy than the other. The fitting results indicate that the jet alignment is dependent on radio luminosity and the shape of the host galaxies, with the jets being more likely to be aligned with the minor axis of the galaxy for lower radio luminosity and for optically more elongated radio-AGNs. The minor-axis alignment of the entire sample may suggest a coherent accretion model present in most AGN host galaxies, while a considerable number of luminous radio-AGNs with massive host galaxies might have undergone accretion – according to the chaotic model – or past merger events.

Astronomy & Astrophysics↗
DOE OSTI2025

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗
DOE OSTI2026

Mechanics of pore array collapse and interaction in shock-compressed polymethyl methacrylate (PMMA)

Recent studies on dynamic pore collapse have revealed significant development of shear localization, which can lead to material failure in porous structures and hot spot generation in energetic materials. These findings have dramatically improved the understanding of failure mechanisms during pore collapse but also prompt further investigation of realistic porous materials. In particular, porous media consist of many pores and porous networks. Even in low-porosity materials, pores can form in close proximity during the manufacturing process, leading to the critical question of pore–pore interaction during collapse under dynamic loading conditions. This study investigates, via plate impact experiments coupled with high-speed internal digital image correlation and shadowgraphy techniques, the collapse of two pores in shock-compressed PMMA at stresses between 0.4 and 1 GPa. The results of these experiments provide new insights into shear localization in pore collapse, in addition to distinct interactions between pores. Shadowgraphy measurements reveal novel, direct visualization of shear band development and crack evolution from pore surfaces. Spacing between adiabatic shear bands is measured over a range of impact stresses and is predicted accurately by the Grady–Kipp model. Pore interactions are found to effect a transition in the impact stress threshold at which different failure mechanisms initiate and are also found to possibly influence preferential sites for shear cracking. Throughout the study, numerical and theoretical models are leveraged to understand shear localization behavior. The role of baroclinicity and wave interactions between the pores is used to elucidate interaction mechanisms between pores.

Lawlor, Barry P. [California Institute of Technolo↗