Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “importance sampling”

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 OSTI2023

MetaboDirect: an analytical pipeline for the processing of FT-ICR MS-based metabolomic data

Background: Microbiomes are now recognized as the main drivers of ecosystem function ranging from the oceans and soils to humans and bioreactors. However, a grand challenge in microbiome science is to characterize and quantify the chemical currencies of organic matter (i.e., metabolites) that microbes respond to and alter. Critical to this has been the development of Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS), which has drastically increased molecular characterization of complex organic matter samples, but challenges users with hundreds of millions of data points where readily available, user-friendly, and customizable software tools are lacking. Results: Here, we build on years of analytical experience with diverse sample types to develop MetaboDirect, an open-source, command-line-based pipeline for the analysis (e.g., chemodiversity analysis, multivariate statistics), visualization (e.g., Van Krevelen diagrams, elemental and molecular class composition plots), and presentation of direct injection high-resolution FT-ICR MS data sets after molecular formula assignment has been performed. When compared to other available FT-ICR MS software, MetaboDirect is superior in that it requires a single line of code to launch a fully automated framework for the generation and visualization of a wide range of plots, with minimal coding experience required. Among the tools evaluated, MetaboDirect is also uniquely able to automatically generate biochemical transformation networks (ab initio) based on mass differences (mass difference network-based approach) that provide an experimental assessment of metabolite connections within a given sample or a complex metabolic system, thereby providing important information about the nature of the samples and the set of microbial reactions or pathways that gave rise to them. Finally, for more experienced users, MetaboDirect allows users to customize plots, outputs, and analyses. Conclusion: Application of MetaboDirect to FT-ICR MS-based metabolomic data sets from a marine phage-bacterial infection experiment and a Sphagnum leachate microbiome incubation experiment showcase the exploration capabilities of the pipeline that will enable the research community to evaluate and interpret their data in greater depth and in less time. It will further advance our knowledge of how microbial communities influence and are influenced by the chemical makeup of the surrounding system. The source code and User’s guide of MetaboDirect are freely available through (https://github.com/Coayala/MetaboDirect) and (https://metabodirect.readthedocs.io/en/latest/), respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bayesian Inference with Latent Hamiltonian Neural Networks (L-HNNs)

When sampling for Bayesian inference, one popular approach is to use Hamiltonian Monte Carlo (HMC) and the No-U-Turn Sampler (NUTS). However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice. We propose Hamiltonian neural networks (HNNs) with HMC and NUTS for solving Bayesian inference problems [1, 2]. Once trained, HNNs do not require gradients of the target density while sampling. Moreover, they satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well suited for use within HMC and NUTS because stationarity can be shown. We also propose an HNN extension called latent HNNs (L-HNNs), which predict latent variable outputs. Compared to HNNs, L-HNNs offer improved expressivity and a reduction in integration errors. Finally, we propose employing L-HNNs in NUTS with an online error monitoring scheme to prevent degeneracy of the sampling in regions of low probability density. We demonstrate L-HNNs in NUTS with online error monitoring by using several example cases involving complex, heavy-tailed, and high local curvature probability densities. Overall, L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. Compared to traditional NUTS, L-HNNs in NUTS with online error monitoring improved the effective sample size (ESS) per gradient by an order of magnitude.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Local cation order and ferrimagnetism in compositionally complex spinel ferrites

We present an exploration of a family of compositionally complex cubic spinel ferrites featuring combinations of Mg, Fe, Co, Ni, Cu, Mn, and Zn cations, systematically investigating the average and local atomic structures, chemical short-range order, magnetic spin configurations, and magnetic properties. All compositions result in ferrimagnetic average structures with extremely similar local bonding environments; however, the samples display varying degrees of cation inversion and, therefore, differing apparent bulk magnetization. Additionally, first-order reversal curve analysis of the magnetic reversal behavior indicates varying degrees of magnetic ordering and interactions, including potentially local frustration. Finally, reverse Monte Carlo modeling of the spin orientation demonstrates a relationship between the degree of cation inversion and the spin collinearity. Collectively, these observations correlate with differences in synthesis procedures. This work provides a framework for understanding magnetic behavior reported for “high-entropy spinels,” revealing many are likely compositionally complex oxides with differing degrees of chemical short-range order—not meeting the community established criteria for high or medium entropy compounds. Moreover, this work highlights the importance of reporting complete sample processing histories and investigating local to long-range atomic arrangements when evaluating potential entropic mixing effects and assumed property correlations in high entropy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Small Angle X-Ray Scattering and Polymers

Polymers are some of the most versatile and useful materials on this planet. Properties of polymers are a result of their processing history, chemical makeup, and physical structure. Therefore, characterization is necessary to develop new polymer formulations and processes. Many established techniques for polymer characterization are destructive, and alternative methods that don’t destroy samples are in high demand. Thermal analysis techniques such as differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) yield important characteristic information for polymer samples while destroying the sample. Non-destructive strategies such as X ray diffraction (XRD) or small angle X-ray scattering (SAXS) can also find polymer characteristics while preserving a sample’s attributes. Sample preservation is desirable when investigating novel formulation and processing methods. SAXS analysis will be done on two different polymers: polyethylene, and a triblock copolymer. Using these two polymers allows for confidence between tests as polyethylene is very common, and well understood. Polyethylene analysis will be done first, and when methods have been refined so that the SAXS data matches accepted thermal analysis values for HDPE, there will be confidence to use the same analysis on block copolymers. Block copolymers remain relatively unexplored, so SAXS data will yield valuable information on their characteristics. Understanding the mechanical properties of these polymers is paramount to their application in upcycling, a process by which waste plastic can be converted into a higher value commodity. For this project, the upcycling process will be done via 3-dimensional (3d) printing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nitrogen fixation rates in the Guinea Dome and the equatorial upwelling regions in the Atlantic Ocean

Abstract Biological nitrogen fixation is a key process balancing the loss of combined nitrogen in the marine nitrogen cycle. Its relevance in upwelling or high nutrient regions is still unclear, with the few available studies in these regions of the ocean reporting rates that vary widely from below detection limit to > 100 nmol N L −1 d −1 . In the eastern tropical Atlantic Ocean, two open ocean upwelling systems are active in boreal summer. One is the seasonal equatorial upwelling, where the residual phosphorus associated with aged upwelled waters is suggested to enhance nitrogen fixation in this season. The other is the Guinea Dome, a thermal upwelling dome. We conducted two surveys along 23° W across the Guinea Dome and the Equator from 15° N to 5° S in September 2015 and August–September 2016 with high latitudinal resolution (20–60 nm between stations). The abundance of Trichodesmium colonies was characterized by an Underwater Vision Profiler 5 and the total biological nitrogen fixation in the euphotic layer was measured using the 15 N 2 technique. The highest abundances of Trichodesmium colonies were found in the area of the Guinea Dome (9°–15° N) with a maximum of 3 colonies L −1 near the surface. By contrast, colonies were almost absent in the Equatorial band between 2° N and 5° S. The highest nitrogen fixation rate was measured at the northern edge of the Guinea Dome in 2016 (ca. 31 nmol N L −1 d −1 ). In this region, where diazotrophs thrived on a sufficient supply of both phosphorus and iron, a patchy distribution was unveiled by our increased spatial resolution scheme. In the Equatorial band, rates were considerably lower, ranging from below detection limit to ca. 4 nmol N L −1 d −1 , with a clear difference in magnitude between 2015 (rates close to zero) and 2016 (average rates around 2 nmol N L −1 d −1 ). This difference seemed triggered by a contrasting supply of phosphorus between years. Our study stresses the importance of surveys with sampling at fine-scale spatial resolution, and shows unexpected high variability in the rates of nitrogen fixation in the Guinea Dome, a region where diazotrophy is a significant process supplying new nitrogen into the euphotic layer.

Fernández-Carrera, Ana (ORCID:0000000183470996)↗
DOE OSTI2024

Multi-level Monte Carlo methods in chemical applications with Lennard-Jones potentials and other landscapes with isolated singularities

We describe and compare outcomes of various Multi-Level Monte Carlo (MLMC) method variants, motivated by the potential of improved computational efficiency over rejection based Monte Carlo, which scales poorly with problem dimension. With an eye toward its application to computational chemical physics, we test MLMC's ability to sample trajectories on two problems — a familiar double-well potential, with known stationary distributions, and a Lennard-Jones solid potential (a Galton Board). By sampling Brownian motion trajectories, we are able to compute expectations of observable averages. These multi-basin potential energy problems capture the essence of the challenges with using MLMC, namely, maintaining correspondence of sample paths as time-resolution is varied. Addressing this challenge properly can lead to MLMC significantly outperforming standard Monte Carlo path sampling. We describe the essence of this problem and suggest strategies that circumvent diverging multilevel sample paths for an important class of problems. In the tests we also compare the computational cost of several, “adaptive,” variants of MLMC. Our results demonstrate that MLMC overcomes the collision, time scale limitation of the more familiar Brownian path MC samplers, and our implementation provides tunable error thresholds, making MLMC a promising candidate for application to larger and more complex molecular systems.

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

Extended X-ray absorption fine structure of dynamically-compressed copper up to 1 terapascal

Abstract Large laser facilities have recently enabled material characterization at the pressures of Earth and Super-Earth cores. However, the temperature of the compressed materials has been largely unknown, or solely relied on models and simulations, due to lack of diagnostics under these challenging conditions. Here, we report on temperature, density, pressure, and local structure of copper determined from extended x-ray absorption fine structure and velocimetry up to 1 Terapascal. These results nearly double the highest pressure at which extended x-ray absorption fine structure has been reported in any material. In this work, the copper temperature is unexpectedly found to be much higher than predicted when adjacent to diamond layer(s), demonstrating the important influence of the sample environment on the thermal state of materials; this effect may introduce additional temperature uncertainties in some previous experiments using diamond and provides new guidance for future experimental design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The X-ray Focusing System at the Time-Resolved AMO Instrument

The emergence of X-ray Free Electron Lasers (XFELs) in recent years has opened up new fields of research, with the merging of short wavelengths and ultrafast timescales. As XFEL facilities including the Linac Coherent Light Source (LCLS) have matured, there is renewed interest in pushing the performance boundaries of experimental endstations, as well as improving ease of operation for experiments. One of the major areas of improvement is that of X-ray micro- and nano-focusing. The focus quality is directly linked to the peak intensity at the interaction point with a sample, and is especially important for intensity-dependent (e.g. x-ray pump) studies and imaging efforts such as single particle imaging.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Synthetic spectra for Lyman- α forest analysis in the Dark Energy Spectroscopic Instrument

Synthetic data sets are used in cosmology to test analysis procedures, to verify that systematic errors are well understood and to demonstrate that measurements are unbiased. In this work we describe the methods used to generate synthetic datasets of Lyman-α quasar spectra aimed for studies with the Dark Energy Spectroscopic Instrument (DESI). In particular, we focus on demonstrating that our simulations reproduces important features of real samples, making them suitable to test the analysis methods to be used in DESI and to place limits on systematic effects on measurements of Baryon Acoustic Oscillations (BAO). We present a set of mocks that reproduce the statistical properties of the DESI early data set with good agreement. Additionally, we use a synthetic dataset to forecast the BAO scale constraining power of the completed DESI survey through the Lyman-α forest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nontrivial critical behavior at magnetic transitions: A case study of Sm 7 ⁢Pd 3

We present a comprehensive analysis of the critical behavior of Sm 7 ⁢Pd 3 in the vicinity of its second-order magnetoelastic transition at 𝑇 c =173 K. The critical exponents (CEs) 𝛽 and 𝛾, determined using both the standard convergence procedure and the average normalized slope (ANS) method, diverge at 𝑇 c –-a characteristic typically associated with first-order transitions. Notably, none of the established universality classes satisfactorily describe the critical behavior of Sm 7 ⁢Pd 3 , and we discuss the possible origins of this deviation in the context of the strong spin-lattice coupling intrinsic to the sample. We emphasize the importance of accurately selecting the critical temperature and magnetic field ranges to ensure robust critical behavior analysis, and propose a quantitative approach to assess the reliability of the extracted CEs. Additionally, we demonstrate that in the ANS method, the critical exponents 𝛽 and 𝛾 should be calculated separately using data for 𝑇 ⩽ 𝑇 c and 𝑇 ⩾ 𝑇 c , respectively. In conclusion, our findings underscore the need for a revised theoretical framework to accurately describe second-order magnetoelastic transitions.

Ferrimagnets↗
DOE OSTI2020

Thermal dissipation sensors enter a new age: Navigating frontiers in transpiration and hydrologic function

Thermal dissipation (TD) sensors have been used extensively for a wide array of sap flow applications, and yet recently several limitations of this approach have been identified, including a tendency to underestimate absolute flows, especially in small stems with very high flows. However, when used properly, the TD method remains the most cost effective and reliable approach for determining relative flows in comparison studies. Major advantages of this method include the ability to replicate maximally within trees, across trees of different sizes, and across multiple species to capture the full range of natural variation in the study system. From early studies, such variation has been deemed a major challenge for representing stand-level dynamics. In a series of recent comparison studies, we determined mean sap flux density dynamics within and among groups of trees. These studies were conducted within temperate coniferous forest (n=46), temperate oak savanna (n=36), temperate pine forest (n=49), temperate bottomland hardwood forest (n=23), tropical premontane rainforest (n=43), semiarid subtropical shrubland (n=25), and tropical dry forest (n=15), with a total of 237 sensors. Two study sites were also equipped with eddy covariance systems for determination of stand-level evapotranspiration. Sites subject to high background thermal gradients used the transient thermal dissipation method. Comparisons varied between studies but tended to focus on relative flow differences between species, between wetter and drier microsites, between understory and overstory components, or between wetter and drier periods of the growing season. Except for the shrubland sites, study trees tended to be at least 15-cm diameter and frequently exceeded 50-cm diameter. Results of this cross-site comparison highlight the inherent variation in natural stands and the importance of using large sample sizes. The TD approach remains a valuable and preferred method for the future determination of water use in trees. However, sensor replication can fundamentally impact a study's outcomes and should be more carefully considered in study design. Customized error mitigation strategies are best to address the sources of variation most problematic for a particular study. In conclusion, the development of “smart” calibration approaches that correct for fundamental effects of radial variation inherent in thermal dissipation studies is discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Phase Picking Beyond Local Distances: Where Waveform Filtering Still Matters for Deep Learning Models

Waveform filtering is a standard step in traditional seismic phase picking but often receives little attention in deep learning workflows, where models are typically trained on raw or minimally processed waveforms. Although this strategy performs well for local events, we show that performance can degrade substantially at regional distances. To address this limitation, we introduce two ways to incorporate multiband-filtered waveforms into deep learning phase pickers. The stacking approach concatenates filtered inputs along the channel dimension, while the branching approach processes each frequency band through a dedicated network branch before feature fusion. Both approaches can substantially improve performance across epicentral distances of 0° to 20°, but their effectiveness depends strongly on the selected frequency bands. Tests with multiple filter banks show that filter-bank design should be treated as part of model optimization rather than as a fixed preprocessing choice. Grad-CAM analysis of the branching model indicates that band importance varies among waveform samples and across training realizations, with only a weak overall preference for the 0.25 to 0.5 Hz band. These results show that no single filter band is consistently optimal and demonstrate that explicit feature engineering remains valuable for robust deep learning-based seismic phase picking.

58 GEOSCIENCES↗
DOE OSTI2020

Use of Remote Sensing and In-Situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the ACME Model

The overachieving goal of the whole CMDV-MCS project is to improve understanding of warm season continental convection and to develop treatments of convection and microphysics capable of representing mesoscale convective systems (MCSs) features in large-scale models. Our tasks for this project contributing to the overachieving goal include: (1) Improve the ice nucleation formulation for MG2 and P3 cloud microphysics schemes; (2) Improve the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM; and (3) Test the performance of improved ice microphysics in E3SM with observation data. In this project, we have (1) Improved the ice nucleation parameterization for MG2 and P3 in E3SM by implementing two advanced empirical parameterizations with connection to aerosols. The two deterministic heterogeneous ice nucleation parameterizations (i.e., DeMott et al., 2015; Niemand et al., 2012) were merged with the MG2 and P3 cloud microphysics schemes in E3SM. Long-term simulations were conducted to examine the impacts of these new parameterizations on simulated cloud properties; (2) Improved the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM. We evaluated the double Gaussian PDF of vertical velocity simulated by the Cloud Layers Unified By Binormals (CLUBB) and the sub-column vertical velocity sampled from the Subgrid Importance Latin Hypercube Sampler (SILHS) in E3SM. We introduced the vertical velocity variance induced by topographic gravity waves for ice nucleation and droplet activation; and (3) Tested the performance of improved ice microphysics in E3SM with observation data. We tested the new treatments of ice nucleation in the single column model (SCM) mode for the stratiform mixed-phase clouds observed during 9-10 October 2004 in the DOE ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) and for the convective clouds observed on 20 May 2011 in the Midlatitude Continental Convective Clouds Experiment (MC3E). Modeled ice nucleating particles (INPs) concentrations were compared against observations collected around the globe.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Imaging 192/193m Ir sources using digital autoradiography for nuclear forensic applications [Slides]

This work aims to develop a new methodology for assaying 192/193m Ir-containing materials. SDDs cannot determine activity spatially, thus it is important that the prepared samples are uniformly distributed. Autoradiography is used to image radioactive samples using imaging media by direct exposure. Using this technique, radio-iridium samples will be imaged to determine uniformity and self-attenuation as a function of mass.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

MethodOpt: a Shiny-based graphical user interface for multivariate optimization of sampling and analytical instrumentation

Method optimization is an important step in producing useful data in various experimental settings involving the use of sampling and analytical instrumentation, such as gas-chromatography mass-spectrometry or other analytical techniques. However, traditional optimization techniques often lack the sophistication of more modern optimization techniques developed in areas of applied mathematics. A graphical user interface has been developed that implements a multivariate, multi-objective optimization technique for spectra-generating sampling and analytical instrumentation, which saves substantial time and resources compared to the more traditional approaches to method development.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Giant kelp genetic monitoring before and after disturbance reveals stable genetic diversity in Southern California

Given the impacts of climate change and other anthropogenic stressors on marine systems, there is a need to accurately predict how species respond to changing environments and disturbance regimes. The use of genetic tools to monitor temporal trends in populations gives ecologists the ability to estimate changes in genetic diversity and effective population size that may be undetectable by traditional census methods. Although multiple studies have used temporal genetic analysis, they usually involve commercially important species, and rarely sample before and after disturbance. In this study, we run a temporal analysis of giant kelp, Macrocystis pyrifera , genetic diversity over the scope of 10 years (2008-2018) using the same microsatellite marker panel to assess the genetic consequences of disturbance in several populations of giant kelp ( Macrocystis pyrifera ) in the Southern California Bight. The study is a rare pre- and post-disturbance microsatellite analysis that included declines to giant kelp caused by the 2015/16 El Nino Southern Oscillation event. We used canopy biomass estimated by remote sensing (Landsat) to quantify the extent of disturbance to kelp beds, and sea surface temperature data to understand how kelp was pushed towards its temperature limits during this period. Despite prolonged periods with decreased canopy at several sites, no changes in genetic structure and allelic richness were observed. We argue that giant kelp in the region is best described as a “patchy population” system where true extinctions are rare. We discuss how deep refugia of subsurface sporophytes and cryptic microscopic life stages could have kept genetic diversity through disturbance. Given the increasing effects of climate change and uncertainty in modeling impacts of species with cryptic life history stages, we suggest further investigation to reveal the role such stages play in species resilience. Genetic monitoring studies of sites selected by remote census demographic and climate surveys should be continued in the future given the predicted impacts of climate change.

Klingbeil, III, William H.↗
DOE OSTI2023

CRITICAL MINERAL PARTITIONING IN COAL-HOSTED CLAYS OF THE POWDER RIVER BASIN, WY DETERMINED BY SEQUENTIAL EXTRACTION

Increasing demand for a more robust domestic supply of rare earth elements (REE) and critical minerals (CM) has led to significant investigation into unconventional sources. Coal and coal byproducts from the Powder River Basin (PRB) of Wyoming are potential sources of REE and CM. The association between strategic metals and their mineral hosts has important implications for extractability. Samples of overburden, underclays, and clay-rich partings were selected for analysis from core and bucket samples from three PRB coal mines. Quartz, illite, and kaolinite-group minerals are abundant in clay-rich portions of the PRB coal strata. We performed a five-step sequential extraction to constrain critical metal partitioning in samples of overburden and underclays (n=10). Unlike the occurrences of REE in other clay deposits, insignificant amounts of total REE + Y (TREY) were recovered in the ion-exchangeable fraction (< 0.5 ppm). A significant amount of TREY was leached from the acid-soluble fraction by 0.1M HCl, suggesting REE may be tied up within crystal lattices rather than sorbed to clay surfaces. The acid-soluble fraction contained between 3 and 250 ppm TREY, with an average TREY concentration of 60 ppm. Relative standard deviations for each extraction step were generally < ±10%. Whole-rock samples (n=72) contained greater TREY enrichment in overburden than in underclays and partings, however, the opposite is true for other CM, such as Ti and V. Other locations in the PRB show different enrichment trends. Our findings supplement a broader effort to couple REE and CM extraction with existing domestic coal production and contribute to the fundamental understanding of how these metals concentrate in low-temperature basin environments.

Stuart, Sophia↗
DOE OSTI2021

A new long-term sampling approach to viruses on surfaces

The importance of virus disease outbreaks and its prevention is of growing public concern but our understanding of virus transmission routes is limited by adequate sampling strategies. While conventional swabbing methods provide merely a microbial snapshot, an ideal sampling strategy would allow reliable collection of viral genomic data over longer time periods. This study has evaluated a new, paper-based sticker approach for collection of reliable viral genomic data over longer time periods up to 14 days and after implementation of different hygiene measures. In contrast to swabbing methods, which sample viral load present on a surface at a given time, the paper-based stickers are attached to the surface area of interest and collect viruses that would have otherwise been transferred onto that surface. The major advantage of one-side adhesive stickers is that they are permanently attachable to a variety of surfaces. Initial results demonstrate that stickers permit stable recovery characteristics, even at low virus titers. Stickers also allow reliable virus detection after implementation of routine hygiene measures and over longer periods up to 14 days. Overall, results for this new sticker approach for virus genomic data collection are encouraging, but further studies are required to confirm anticipated benefits over a range of virus types.

59 BASIC BIOLOGICAL SCIENCES↗