Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RESPONSE FUNCTIONS”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗
DOE OSTI2025

time-resolved spectroscopy fit (trspecfit) v0.01

Analyze 2D time- and energy-resolved data, such as from a pump-probe spectroscopy experiment. User can select and input different peak shapes/ functions and background types to first fit a ground state/ unperturbed spectrum. This would be similar to how standard spectroscopy data is fit. Subsequently, to describe the time domain, users can choose functions that describe the temporal dynamics of one or more spectral features, such as a peak amplitude, peak position, etc. These time dynamics functions can be added or convoluted (e.g. describing an instrument response function) with each other. Functionality to integrate implicit variables leading to distributions of certain parameters/ functions is in development. Alternatively, 2D data can be analyzed one time step at a time to get an idea of the time dynamics of the system before deploying the global 2D fit described above. Typically people write custom software for this purpose. During my PhD I've seen five internal LBL and external researchers write one-off code in different languages to analyze time- and energy-resolved spectra. While this was specifically was for a laser pump - X-ray probe spectroscopy experiment, I'm trying to write a general package for the time-resolved spectroscopy community.

Mahl, Johannes [Lawrence Berkeley National Laborat↗
DOE OSTI2023

Demonstrating Safeguards Applications of DRiFT Gas Detector Capabilities [Slides]

DRiFT, a Detector Response Function Toolkit, post-processes MCNP output to model detector response. Three primary components: scintillator (focus of FY 21 TED work), gas detectors (focus of this FY), and semiconductor detectors. Capability to model scintillator detectors has previously been demonstrated.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Forward Modeling of Gamma Reaction History Signatures From Anticipated Deuterium-Tritium Filled MagLIF Implosions on Sandia’s Z-Machine

Nuclear reaction history measurements provide a bang time and burn width of Inertial Confinement Fusion (ICF) implosions and are essential for understanding implosion performance to constrain ICF capsule design. When fusion fuel contains Deuterium (D) and Tritium (T) gas, reaction history is informed by measuring the 16.75 MeV gamma rays generated from the D(T,γ) 5 He reaction. Such DT based reaction history measurements have not been made on the Magnetized Laser Inertial Fusion (MagLIF) platform on Sandia’s Z-Machine due to the lack of Tritium being used. The recent development of ICF implosions with tritiated fuel will open the possibility of measuring the gamma reaction history on the Z-Machine. A forward model of the Gamma Reaction History diagnostic on Z (GRH-Z) has been developed using the MCNP6.3 (Monte-Carlo N-Particle) radiation transport code. The model included the Z-Machine geometry of interest to characterize the impact of neutron induced gamma rays on the DT reaction history measurements. In addition, the impulse response functions of the GRH-Z diagnostic to understand the temporal response of the detector, and the minimum yields required to make a reaction history measurement were calculated. This approach also predicted that with T 2 gas doping of MagLIF implosions a reaction history may be made for high performance shots >8e12-2.4e13 depending on the chosen threshold for the detector, with a maximum signal to background ratio of 25%. It was found that for long duration ICF implosions that additional collimation will be needed to prevent the neutron induced gamma rays from modifying the shape of the measured DT reaction history curve.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Evaluation of the GaAs displacement damage metric using updated nuclear data

The emerging use of the physics-based athermal recombination-corrected displacement per atom (arc-dpa) model for the displacement damage efficiency has motivated a re-evaluation of the historical empirically-derived GaAs damage response function with the purpose of highlighting needs for future analytical and experimental work. The 1-MeV neutron damage equivalence methodology used in the ASTM E-722 standard for GaAs has been re-evaluated using updated nuclear data. This yielded a higher fidelity representation of the GaAs displacement kerma and, through the use of the refined PKA recoil energy-dependent damage efficiency model, an updated 1-MeV(GaAs) displacement damage function. This re-evaluation included use of the Norgett-Robinson-Torrens (NRT) model for an updated threshold treatment, rather than the sharp-threshold Kinchin-Pease model used in the current ASTM standard. The underlying nuclear data evaluations have been updated to use the ENDF/VIII.0 {sup 75}As and TENDL-2019 {sup 71}Ga/{sup 69}Ga evaluations. The displacement kerma and 1-MeV-equivalent damage responses were calculated using a modified NJOY-2016 code which allowed for refinements in some of the damage models. This paper shows that an updated displacement damage function, based upon the latest nuclear data, is consistent with the experimental data used to develop the current ASTM E-722 GaAs standard. Using a double ratio approach to compare the available experimental data with the calculated response, the average legacy double ratio was found to be 0.97 ± 0.05 and the average updated double ratio was found to be 0.94 ± 0.05. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2024

Neural Networks to Find the Optimal Forcing for Offsetting the Anthropogenic Climate Change Effects

Abstract Of great relevance to climate engineering is the systematic relationship between the radiative forcing to the climate system and the response of the system, a relationship often represented by the linear response function (LRF) of the system. However, estimating the LRF often becomes an ill-posed inverse problem due to high-dimensionality and nonunique relationships between the forcing and response. Recent advances in machine learning make it possible to address the ill-posed inverse problem through regularization and sparse system fitting. Here, we develop a convolutional neural network (CNN) for regularized inversion. The CNN is trained using the surface temperature responses from a set of Green’s function perturbation experiments as imagery input data together with data sample densification. The resulting CNN model can infer the forcing pattern responsible for the temperature response from out-of-sample forcing scenarios. This promising proof of concept suggests a possible strategy for estimating the optimal forcing to negate certain undesirable effects of climate change. The limited success of this effort underscores the challenges of solving an inverse problem for a climate system with inherent nonlinearity. Significance Statement Predicting the climate response for a given climate forcing is a direct problem, while inferring the forcing for a given desired climate response is often an inverse, ill-posed, problem, posing a new challenge to the climate community. This study makes the first attempt to infer the radiative forcing for a given target pattern of global surface temperature response using a deep learning approach. The resulting deeply trained convolutional neural network inversion model shows promise in capturing the forcing pattern corresponding to a given surface temperature response, with a significant implication on the design of an optimal solar radiation management strategy for curbing global warming. This study also highlights the technical challenges that future research should prioritize in seeking feasible solutions to the inverse climate problem.

Ren, Huiying↗
DOE OSTI2022

Physics Discovery in Nanoplasmonic Systems via Autonomous Experiments in Scanning Transmission Electron Microscopy

Abstract Physics‐driven discovery in an autonomous experiment has emerged as a dream application of machine learning in physical sciences. Here, this work develops and experimentally implements a deep kernel learning (DKL) workflow combining the correlative prediction of the target functional response and its uncertainty from the structure, and physics‐based selection of acquisition function, which autonomously guides the navigation of the image space. Compared to classical Bayesian optimization (BO) methods, this approach allows to capture the complex spatial features present in the images of realistic materials, and dynamically learn structure–property relationships. In combination with the flexible scalarizer function that allows to ascribe the degree of physical interest to predicted spectra, this enables physical discovery in automated experiment. Here, this approach is illustrated for nanoplasmonic studies of nanoparticles and experimentally implemented in a truly autonomous fashion for bulk‐ and edge plasmon discovery in MnPS 3 , a lesser‐known beam‐sensitive layered 2D material. This approach is universal, can be directly used as‐is with any specimen, and is expected to be applicable to any probe‐based microscopic techniques including other STEM modalities, scanning probe microscopies, chemical, and optical imaging.

42 ENGINEERING↗
DOE OSTI2021

Quantifying the impacts of compound extremes on agriculture

Agricultural production and food prices are affected by hydroclimatic extremes. There has been a growing amount of literature measuring the impacts of individual extreme events (heat stress or water stress) on agricultural and human systems. Yet, we lack a comprehensive understanding of the significance and the magnitude of the impacts of compound extremes. This study combines a fine-scale weather product with outputs of a hydrological model to construct functional metrics of individual and compound hydroclimatic extremes for agriculture. Then, a yield response function is estimated with individual and compound metrics, focusing on corn in the United States during the 1981–2015 period. Supported by statistical evidence, the findings suggest that metrics of compound hydroclimatic extremes are better predictors of corn yield variations than metrics of individual extremes. The results also confirm that wet heat is more damaging than dry heat for corn. This study shows the average yield damage from heat stress has been up to four times more severe when combined with water stress.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2024

Mortality attributable to PM 2.5 from wildland fires in California from 2008 to 2018

In California, wildfire risk and severity have grown substantially in the last several decades. Research has characterized extensive adverse health impacts from exposure to wildfire-attributable fine particulate matter (PM 2.5 ), but few studies have quantified long-term outcomes, and none have used a wildfire-specific chronic dose-response mortality coefficient. Here, we quantified the mortality burden for PM 2.5 exposure from California fires from 2008 to 2018 using Community Multiscale Air Quality modeling system wildland fire PM 2.5 estimates. We used a concentration-response function for PM 2.5 , applying ZIP code–level mortality data and an estimated wildfire-specific dose-response coefficient accounting for the likely toxicity of wildfire smoke. We estimate a total of 52,480 to 55,710 premature deaths are attributable to wildland fire PM 2.5 over the 11-year period with respect to two exposure scenarios, equating to an economic impact of $\$432$ to $\$456$ billion. These findings extend evidence on climate-related health impacts, suggesting that wildfires account for a greater mortality and economic burden than indicated by earlier studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Activity-based protein profiling identifies alternating activation of enzymes involved in the bifidobacterium shunt pathway or mucin degradation in the gut microbiome response to soluble dietary fiber

While deprivation of dietary fiber has been associated with adverse health outcomes, investigations concerning the effect of dietary fiber on the gut microbiome have been largely limited to compositional sequence-based analyses or utilize a defined microbiota not native to the host. To extend understanding of the microbiome’s functional response to dietary fiber deprivation beyond correlative evidence from sequence-based analyses, approaches capable of measuring functional enzymatic activity are needed. In this study, we use an activity-based protein profiling (ABPP) approach to identify sugar metabolizing and transport proteins in native mouse gut microbiomes that respond with differential activity to the deprivation or supplementation of the soluble dietary fibers inulin and pectin. We found that the microbiome of mice subjected to a high fiber diet high in soluble fiber had increased functional activity of multiple proteins, including glycoside hydrolases, polysaccharide lyases, and sugar transport proteins from diverse taxa. The results point to an increase in activity of the Bifidobacterium shunt metabolic pathway in the microbiome of mice fed high fiber diets. In those subjected to a low fiber diet, we identified a shift from the degradation of dietary fibers to that of gut mucins, in particular by the recently isolated taxon “Musculibacterium intestinale”, which experienced dramatic growth in response to fiber deprivation. When combined with metabolomics and shotgun metagenomics analyses, our findings provide a functional investigation of dietary fiber metabolism in the gut microbiome and demonstrates the power of a combined ABPP-multiomics approach for characterizing the response of the gut microbiome to perturbations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Analytical continuation of matrix-valued functions: Carathéodory formalism

Finite-temperature quantum field theories are formulated in terms of Green's functions and self-energies on the Matsubara axis. In multiorbital systems, these quantities are related to positive semidefinite matrix-valued functions of the Carathéodory and Schur class. Analysis, interpretation, and evaluation of derived quantities such as real-frequency response functions requires analytic continuation of the off-diagonal elements to the real axis. Here we derive the criteria under which such functions exist for given Matsubara data and present an interpolation algorithm that intrinsically respects their mathematical properties. For small systems with precise Matsubara data, we find that the continuation exactly recovers all off-diagonal and diagonal elements. In real-materials systems, we show that the precision of the continuation is sufficient for the analytic continuation to commute with the Dyson equation, and we show that the truncation of the off-diagonal self-energy elements leads to considerable approximation artifacts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Scale dependence in functional equivalence and difference in the soil microbiome

Climatic history can shape the functioning of soil microbial communities and thus rates of ecosystem processes such as organic matter decomposition. For example, broad spatial scale differences in climatic history, such as contrasting precipitation regimes, have been shown to generate unique microbial functional responses to contemporary moisture conditions. Yet it is an open question as to whether local differences in soil microclimate similarly influence the functional potential of decomposer communities. Here, we use a multi-scale approach within and among two temperate forest field sites to investigate this question. Soils from fifty-four microsites, that vary in their soil moisture climate-regimes, were used as inocula for a common leaf litter (Quercus rubra L.) in a controlled, laboratory microcosm study. Microcosms were placed under dry, mesic and wet lab-moisture conditions and the rate of carbon (C) mineralization of the litter was measured over 202 days. Our results reveal differences in decomposition rates under controlled conditions that highlight broad-scale functional differences between the soil communities at each site. Specifically, we found that C mineralization differed by as much as two-fold for soil communities when compared between the sites. Our results also show that functional differences of soil communities are observable within one site but not the other. In the site where local-scale functional legacies were apparent, the historical soil moisture microclimate-regimes generated as much as an 89% change in C mineralization rates of the leaf litter under the same contemporary, lab-imposed moisture conditions. Here, a similar pattern was not observable in the other site; instead, laboratory moisture conditions explained almost all variation in C mineralization. Notably, for the site with pronounced local-scale functional legacies, there was much greater within-site variation in field-soil microsite moisture than at the site which did not exhibit functional legacies, suggesting that the extent of local-scale variation in microclimate may act as control on whether local-scale functional legacies are observed. Regardless of whether this mechanism does explain our findings, our observations do confirm those from prior studies where regional-scale moisture-regime differences shape microbial function, and extend this prior work by providing evidence that pronounced local-scale differences in soil moisture microclimate-regimes are associated with microbial functional legacies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of moisture on microbial decomposition phenotypes and enzyme dynamics

Soil organic matter decomposition is a complex process reflecting microbial composition and environmental conditions. Moisture can modulate the connectivity and interactions of microbes. Due to heterogeneity, a deeper understanding of the influence of soil moisture on the dynamics of organic matter decomposition and resultant phenotypes remains a challenge. Soils from a long-term field experiment exposed to high and low moisture treatments were incubated in the laboratory to investigate organic matter decomposition using chitin as a model substrate. By combining enzymatic assays, biomass measurements, and microbial enrichment via activity-based probes, we determined the microbial functional response to chitin amendments and field moisture treatments at both the community and cell scales. Chitinolytic activities showed significant responses to the amendment of chitin, independent of differences in field moisture treatments. However, for other measurements of carbon metabolism and cellular functions, soils from high moisture field treatments had greater potential enzyme activity than soils from low moisture field treatments. A cell tagging approach was used to enrich and quantify bacterial taxa that are actively producing chitin-degrading enzymes. By integrating organism, community, and soil core measurements we show that (i) a small subset of taxa compose the majority (>50%) of chitinase production despite broad functional redundancy, (ii) the identity of key chitin degraders varies with moisture level, and (iii) extracellular enzymes that are not cell-associated account for most potential chitinase activity measured in field soil.

activity-based probes↗
DOE OSTI2024

Uranium measurements in the field using high-resolution cadmium zinc telluride detectors

A new generation of cadmium zinc telluride (CZT) detectors has become available and is being evaluated by the International Atomic Energy Agency (IAEA) for safeguards verifications in the field. The new CZT detector, model M400, is a room temperature spectrometer manufactured by H3D, Inc. The M400 demonstrates superior energy resolution, effective isotope identification capabilities, and convenient usability features when tested in a controlled laboratory environment. These characteristics define the M400 as a potential platform for IAEA field detection applications, which could become suitable for nuclear material characterization (e.g., enrichment verification) and nuclear safeguards missions. The capabilities of gamma spectrometry codes including Fixed energy, Response function Analysis with Multiple efficiencies (FRAM) from Los Alamos National Laboratory, CZT for Uranium (CZTU) from Lawrence Livermore National Laboratory, and Gamma Detector Response and Analysis Software (GADRAS) from Sandia National Laboratories were adapted for M400 spectra, and the performance of the codes has been validated. This was reported in a prior work. To further validate the performance of the high-energy resolution CZT detector and the isotopic analysis codes, a field measurement campaign consisting of uranium hexafluoride (UF 6 ) cylinder measurements was conducted at a fuel fabrication facility. A total of 34 Type 30B cylinders containing UF6 were measured using three different M400 CZT detectors. Each detector was outfitted with a custom rectangular collimator and shield made out T-Flex®, a tungsten-impregnated polymer. Measurements were performed at three different locations of the cylinder, ensuring that the measurement geometry satisfied the infinite thickness criterion. The spectra from the M400 CZT were analyzed using the code General Enrichment Meter (GEM). For analyzing the gamma-ray spectra from UF 6 cylinder, the GEM code is the appropriate tool since it relies only on the gamma-ray emissions from 235 U and not from other isotopes. Results from the spectral analysis were compared with the known abundance of 235 U in the cylinders, as well as with the International Target Values 2020 (ITV2020). The suitability of the different underlying techniques used by the various codes for UF 6 analysis is discussed. The challenges of measuring UF 6 contained in cylinders and mitigation strategies are highlighted.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

DRiFT - Release 1.1.1: Organic Scintillators

DRiFT (a Detector Response Function Toolkit) is LANL-developed software that postprocesses output from the extensively validated radiation transport code, MCNP, and generates realistic nuclear instrumentation response. DRiFT is designed to be flexible, enabling users to specify detector type and many experimental settings, as well as accommodating the addition of their own desired features. Although DRiFT development has included scintillator, gas, and semiconductor features, the focus of this release is on organic scintillator and associated capabilities. Organic scintillators are widely used in the areas of nuclear safeguards and nuclear non-proliferation efforts. DRiFT has several diagnostic and detector physics features relevant to detailed scintillator simulations including: tracking source particle information, scintillation light production, the effects of PMT quantum efficiency and gain, and digitizer settings. Users can select responses from many scintillator and PMT types supported natively by DRiFT, or add their own by following the instructions in this document. We acknowledge that DRiFT is under active development, bug reports and general questions and comments should be directed to Madison Andrews, madison@lanl.gov. This manual is divided into four parts: I) An overview of DRiFT, including how to obtain and install the executable, II) A description of the detector physics related to scintillators available, III) a description of more general DRiFT features the user may find useful, and IV) a description of the test suite and examples made available with the code release.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗