Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional Data Analysis”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

BayFlux: A Bayesian method to quantify metabolic Fluxes and their uncertainty at the genome scale

Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13 C Metabolic Flux Analysis ( 13 C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13 C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13 C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13 C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in “non-gaussian” situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13 C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13 C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P- 13 C MOMA and P- 13 C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Phenomenological Model for Cavitation

A phenomenological model of cavitation is presented, based on observations that both large relative negative pressures and large negative time derivatives of pressure are required for cavitation onset. We simulated two cavitation experiments to generate cavitation scaling parameters for relative pressure drop and rate of pressure drop. Our results show the model, while simple, is effective at reproducing results from laboratory experiments of cavitation. The parameters were then used in conjunction with a human surrogate computational model to predict, at any position within the head, the probability of intracranial cavitation caused by exposure to a blast event. The results suggest that the magnitude of blast overpressure observed in field data is sufficient to cause intracranial cavitation. Our analysis indicates that the helmeted head, when compared to the unhelmeted head configuration, results in a decrease but not elimination of cavitation exposure. When density functions of cavitation probability versus cumulative brain volume are combined with an injury severity model, the results show helmet efficacy at low and moderate risk levels. However, the convergence of unhelmeted and helmeted probability density functions at high-to-excessive risk thresholds indicates the helmet offers diminishing protection at elevated exposure levels, relative to the unhelmeted baseline. Future investigation and collaboration with neuroscience subject matter experts are needed to contextualize the current computational results. While the present work contributes specific and quantified predictions of intracranial cavitation location and severity, more research is required to apply our results to clinical settings with population-based brain injury subjects and controls. The relationship between our intracranial cavitation predictions with their anticipated clinical sequelae remains a topic in need of exploration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Analysis of Air-Purifying Respirator (APR) and Powered Air-Purifying Respirator (PAPR) Cartridge Performance Testing on a Hanford AP Tank Farm Exhauster Slipstream: Volume 2 Raw Analytical Data

As the Tank Operations Contractor for U.S. Department of Energy operations at the Hanford site in Washington State, Washington River Protection Solutions (WRPS) is responsible for managing highly radioactive wastes stored in tanks at Hanford. WRPS tests air-purifying respirator (APR) and powered air-purifying respirator (PAPR) chemical cartridges commonly used at Hanford Tank Farms to determine the period of time that the cartridges would provide adequate performance for APRs and PAPRs used to protect workers when exposed to a mixture of Chemicals of Potential Concern (COPC) from any vapors exiting headspaces in the storage tanks. The Occupational Safety and Health Administration (OSHA) Standard promulgated in Title 29 of the Code of the Federal Regulations (CFR) 1910.134(d)(3)(iii)(b)(2) specifies that for protection against gases and vapors, employers shall implement a schedule for cartridges to ensure that change-outs occur before the end of service life. The change schedule can be based on objective information or data that ensures cartridge change-outs occur before the end of their service life.[2-5] The primary function of the WRPS Cartridge Test Program is to obtain objective data to determine service lives for the APR and PAPR cartridges used at Hanford Tank Farms. WRPS contracted with Pacific Northwest National Laboratory to analyze the test data and offer an independent analysis and any recommendations. Volume 1 of this report summarizes data analysis of APR and PAPR cartridge testing on vapors from the AP tank farm exhauster. Previous testing of APR cartridges was conducted on the AP exhauster in June 2016. However, an AP exhauster upgrade was completed in September 2016. Two different APR cartridges from SCOTT Safety (Monroe, North Carolina) were assessed for the new AP exhauster source, along with two different PAPR cartridges—one from MSA Safety Inc. (Pittsburgh, Pennsylvania) and another from 3M (Maplewood, Minnesota). Volume 2 provides an introduction to the raw data, including analytical laboratory analysis results that supported the analysis and conclusions documented in Volume 1.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Analysis of Air-Purifying Respirator (APR) and Powered Air-Purifying Respirator (PAPR) Cartridge Performance Testing on Hanford Tanks BY-108 and BY-110 Volume 2: Raw Analytical Data

As the Tank Operations Contractor for U.S. Department of Energy operations at the Hanford site in Washington State, Washington River Protection Solutions (WRPS) is responsible for managing highly radioactive wastes stored in tanks at Hanford. WRPS tests air-purifying respirator (APR) and powered air-purifying respirator (PAPR) chemical cartridges commonly used by workers at Hanford Tank Farms to determine the period of time that the cartridges would provide adequate performance for APRs and PAPRs used to protect workers when exposed to a mixture of Chemicals of Potential Concern (COPC) from any vapors exiting headspaces in the storage tanks. The Occupational Safety and Health Administration (OSHA) Standard promulgated in Title 29 of the Code of the Federal Regulations (CFR) 1910.134(d)(3)(iii)(b)(2) specifies that for protection against gases and vapors, employers shall implement a schedule for cartridges to ensure that change-outs occur before the end of service life. The change schedule can be based on objective information or data that ensures cartridge change-outs occur before the end of their service life. The primary function of the WRPS Cartridge Test Program is to obtain objective data to determine service lives for the APR and PAPR cartridges used at Hanford Tank Farms. WRPS contracted with Pacific Northwest National Laboratory to analyze the test data and offer an independent analysis and any recommendations. This report summarizes data analysis of APR and PAPR cartridge testing on BY-110 and BY-108 headspace vapors, respectively. Two different APR cartridges from SCOTT Safety (Monroe, North Carolina) were assessed for the BY-110 headspace vapors, and two different PAPR cartridges—one from MSA Safety Inc. (Pittsburgh, Pennsylvania) and another from 3M (Maplewood, Minnesota)—were assessed for the BY-108 headspace vapors. Volume 1 of this report documents the testing, data analysis, results, conclusions, and recommendations resulting from the cartridge testing of vapors from the BY-110 and BY-108 headspaces. Volume 2 provides an introduction to the raw data, including analytical laboratory analysis results that supported the analysis and conclusions documented in Volume 1.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Analysis of Air-Purifying Respirator (APR) and Powered Air-Purifying Respirator (PAPR) Cartridge Performance Testing on a Hanford AX Tank Farm Exhauster Slipstream Volume 2: Raw Analytical Data

As the Tank Operations Contractor for U.S. Department of Energy operations at the Hanford site in Washington State, Washington River Protection Solutions (WRPS) is responsible for managing highly radioactive wastes stored in tanks at Hanford. To protect workers at Hanford Tank Farms, WRPS tests air-purifying respirator (APR) and powered air-purifying respirator (PAPR) chemical cartridges commonly used at the tank farms. The tests were conducted to determine the period of time the cartridges would provide adequate performance for APRs and PAPRs when workers are exposed to a mixture of Chemicals of Potential Concern (COPC) from any vapors exiting headspaces in the storage tanks. Occupational Safety and Health Administration (OSHA) Standard 29 Code of the Federal Regulations (CFR) 1910.134(d)(3)(iii)(b)(2) specifies that for protection against gases and vapors, employers shall implement a schedule for cartridges to ensure that change-outs occur before the end of service life. The change schedule can be based on objective information or data that ensures cartridge change-outs occur before the end of their service life. The primary function of the WRPS Cartridge Test Program is to obtain objective data to determine service lives for the APR and PAPR cartridges used at Hanford Tank Farms. WRPS contracted with Pacific Northwest National Laboratory to analyze the test data and offer an independent analysis and any recommendations. This report summarizes data analysis of APR and PAPR cartridge testing on a vapor slipstream from the Hanford AX tank farm exhauster. Volume 1 of this report documents the testing, data analysis, results, conclusions, and recommendations resulting from the cartridge testing on AX exhauster slipstream vapors. Volume 2 provides an introduction to the raw data, including analytical laboratory analysis results that supported the analysis and conclusions documented in Volume 1.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Assessment of Mechanical Draft Cooling Tower Thermal Emissions from Visual Images of Plumes

Using a one-dimensional code, we computed the power (enthalpy discharge rate) of a twelve-cell mechanical draft cooling tower (MDCT) using over two hundred visible condensed water vapor plume volume measurements derived from images, weather data, and tower operating conditions. The plume images were simultaneously captured by multiple stationary digital cameras surrounding the cooling tower. An analysis technique combining structure from motion (SfM), a neural-network-based image segmentation algorithm, and space carving was used to quantify the volumes. Afterwards, the power output was computed using novel techniques in the one-dimensional code that included cooling tower exhaust plume adjacency effects implemented with a modified version of the entrainment function, weather data averaged from eleven stations, and fan operations at the times when plume volumes were measured. The model was then compared with the averaged observed power output, and it validated well with an average error ranging from 6 to 12%, depending on the meteorological data used in the simulations. This methodology can possibly determine power plant fuel consumption rates by applying visible imagery.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

COMPASO : A new halo finder for competitive assignment to spherical overdensities

We describe a new method (compaso) for identifying groups of particles in cosmological N-body simulations. compaso builds upon existing spherical overdensity (SO) algorithms by taking into consideration the tidal radius around a smaller halo before competitively assigning halo membership to the particles. In this way, the compaso finder allows for more effective deblending of haloes in close proximity as well as the formation of new haloes on the outskirts of larger ones. This halo-finding algorithm is used in the abacussummit suite of N-body simulations, designed to meet the cosmological simulation requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. compaso is developed as a highly efficient on-the-fly group finder, which is crucial for enabling good load-balancing between the GPU and CPU and the creation of high-resolution merger trees. In this paper, we describe the halo-finding procedure and its particular implementation in abacus, accompanying it with a qualitative analysis of the finder. We test the robustness of the compaso catalogues before and after applying the cleaning method described in an accompanying paper and demonstrate its effectiveness by comparing it with other validation techniques. We then visualize the haloes and their density profiles, finding that they are well fit by the NFW formalism. Finally, we compare other properties such as radius–mass relationships and two-point correlation functions with that of another widely used halo finder, rockstar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Constraining dark matter microphysics with the annihilation signal from subhalos

In the cold dark matter scenario, galactic dark matter halos are populated with a large number of smaller subhalos. Previous work has shown that dark matter annihilations in subhalos can generate a distinctive, non-Poisson signal in the gamma-ray photon counts probability distribution function (PDF). Here we show that the gamma-ray PDF also carries information about the velocity dependence of the dark matter annihilation cross section. After calculating the PDF assuming s-wave and Sommerfeld-enhanced annihilation, we perform a mock data analysis to illustrate how current and future observations can constrain the microphysics of the dark matter annihilation. We find that, with current Fermi data, and assuming a dark matter annihilation cross section roughly at the limit of current bounds from annihilation in dwarf spheroidal galaxies, one can potentially distinguish the non-Poissonian fluctuations expected from dark matter annihilation in subhalos from Poisson sources, as well as from dark matter models with an incorrect velocity-dependence. Here we explore how robust these results are to assumptions about the modeling of the galactic gamma-ray background, but further work is needed to determine the impact of realistic astrophysical source populations on our results. We also point out a four-parameter degeneracy between the velocity dependence of the dark matter annihilation, the minimum subhalo mass, the power law index of the subhalo mass function, and the normalization of the dark matter signal. This degeneracy can be broken with priors from N-body simulations or from observational constraints on the subhalo mass function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Characterization of nuclear breakup as a function of hard-scattering kinematics using dijets measured by ATLAS in p+Pb collisions

This Letter analyzes the sensitivity of event geometry estimators to the initial-state kinematics of hard scattering in proton–lead collisions. This analysis uses dijets as a proxy for the parton–parton scattering configuration, correlating it with event geometry estimators, namely the energy deposited in the Zero-Degree Calorimeter and the transverse energy recorded in the Forward Calorimeter in the Pb-going direction. The analysis uses data recorded by the ATLAS detector at the Large Hadron Collider with a nucleon–nucleon center-of-mass energy of 8.16 TeV, corresponding to an integrated luminosity of 56 nb − 1 . The jets are measured within the pseudorapidity interval −2.8 < η < 4.5, where positive η values correspond to the direction of the proton beam. Results are presented as a function of the Bjorken-x of the parton originating from the proton, xp . Both event geometry estimators are found to be dependent on xp , with the energy deposited in the Zero-Degree Calorimeter about six times less sensitive to xp compared with the transverse energy deposited in the Forward Calorimeter.

Aad, G↗
DOE OSTI2023

First test of the consistency relation for the large-scale structure using the anisotropic three-point correlation function of BOSS DR12 galaxies

ABSTRACT We present, for the first time, an observational test of the consistency relation for the large-scale structure (LSS) of the Universe through a joint analysis of the anisotropic two- and three-point correlation functions (2PCF and 3PCF) of galaxies. We parameterize the breakdown of the LSS consistency relation in the squeezed limit by Es, which represents the ratio of the coefficients of the shift terms in the second-order density and velocity fluctuations. Es ≠ 1 is a sufficient condition under which the LSS consistency relation is violated. A novel aspect of this work is that we constrain Es by obtaining information about the non-linear velocity field from the quadrupole component of the 3PCF without taking the squeezed limit. Using the galaxy catalogues in the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12, we obtain $E_{\rm s} = -0.92_{-3.26}^{+3.13}$, indicating that there is no violation of the LSS consistency relation in our analysis within the statistical errors. Our parameterization is general enough that our constraint can be applied to a wide range of theories, such as multicomponent fluids, modified gravity theories, and their associated galaxy bias effects. Our analysis opens a new observational window to test the fundamental physics using the anisotropic higher-order correlation functions of galaxy clustering.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2023

Structure-thermal property relationships of aikinite PbCuBiS 3

Understanding the structural properties and chemical bonding in complex multinary chalcogenides is highly valuable, as it provides vital information for controlling their physical properties to suit various technological applications. Here, we report on the structural, thermal, and electronic properties of PbCuBiS 3 by incorporating an experimental and theoretical investigation of this material. Analysis of the experimental temperature-dependent heat capacity and thermal conductivity data show that weak bonding and strong lattice anharmonicity give rise to ultralow thermal conductivity. Chemical bonding information obtained from the density functional theory electronic calculations reveal significant antibonding interactions that underscore the dominant role of the Cu-S tetrahedra on occupied antibonding states, as well as a highly distorted Pb/Bi-S environment inducing weak bonding and strong lattice anharmonicity. This study broadens our understanding of the structure-property relationships for this material, and can be useful for developing this and similar chalcogenide materials for potential applications of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multioutput Convolutional Neural Network for Improved Parameter Extraction in Time-Resolved Electrostatic Force Microscopy Data

Time-resolved scanning probe microscopy methods, like time-resolved electrostatic force microscopy (trEFM), enable imaging of dynamic processes ranging from ion motion in batteries to electronic dynamics in microstructured thin film semiconductors for solar cells. Reconstructing the underlying physical dynamics from these techniques can be challenging due to the interplay of cantilever physics with the actual transient kinetics of interest in the resulting signal. Previously, quantitative trEFM used empirical calibration of the cantilever or feed-forward neural networks trained on simulated data to extract the physical dynamics of interest. Both these approaches are limited by interpreting the underlying signal as a single exponential function, which serves as an approximation but does not adequately reflect many realistic systems. Here, we present a multi-branched, multi-output convolutional neural network (CNN) that uses the trEFM signal in addition to the physical cantilever parameters as input. The trained CNN accurately extracts parameters describing both single-exponential and bi-exponential underlying functions, and more accurately reconstructs real experimental data in the presence of noise. This article demonstrates an application of physics-informed machine learning to complex signal processing tasks, enabling more efficient and accurate analysis of trEFM.

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

Fluoropolymer Composites from Partially Perfluoroalkylated Waste Polyethylene

Chemically modified plastics have emerged as practical solutions to plastic waste increases. Here, the inherent novelty of decorating polymer chains with chemical functionality results in distinct properties that expand the available application space. Nevertheless, developing designer materials for specific applications beyond compatibilization or mild property enhancement is difficult due to the synergistic effects of both the polar functionality imparted and the parent materials' intrinsic properties. By incorporating perfluoro-alkyl side-chains onto the backbone of dehydrogenated waste HDPE, unique surface properties intermediate between polytetrafluoroethylene (PTFE, the model fluoropolymer) and HDPE become apparent, while the overall material mechanical and thermal properties result in more LLDPE-like materials. This is demonstrated through moderate decreases in the surface free energy of the perfluoroalkylated polyolefin surface (increase in H 2 O contact angle of ~ 6°) and increased ordering under shear when blended with PTFE nanoparticles where the crossover point occurred at higher strains. Critically, perfluoroalkylated HDPE possesses improved rheological modification properties at elevated temperatures with PTFE nanoparticles, resulting in more thermally robust and stable composite materials.

fluoropolymer↗
DOE OSTI2023

Autonomous x-ray scattering

Autonomous experimentation (AE) is an emerging paradigm that seeks to automate the entire workflow of an experiment, including—crucially—the decision-making step. Beyond mere automation and efficiency, AE aims to liberate scientists to tackle more challenging and complex problems. We describe our recent progress in the application of this concept at synchrotron x-ray scattering beamlines. We automate the measurement instrument, data analysis, and decision-making, and couple them into an autonomous loop. We exploit Gaussian process modeling to compute a surrogate model and associated uncertainty for the experimental problem, and define an objective function exploiting these. We provide example applications of AE to x-ray scattering, including imaging of samples, exploration of physical spaces through combinatorial methods, and coupling to in situ processing platforms These uses demonstrate how autonomous x-ray scattering can enhance efficiency, and discover new materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Measurement of the Ar ( e , e ' p ) and Ti ( e , e ' p ) cross sections in Jefferson Lab Hall A

The E12-14-012 experiment, performed in Jefferson Lab Hall A, has collected exclusive electron-scattering data ( e , e ' p ) in parallel kinematics using natural argon and natural titanium targets. In this study we report the first results of the analysis of the data set corresponding to beam energy 2222 GeV, electron scattering angle 21 . 5 ° , and proton emission angle – 50 ° . The differential cross sections, measured with ≈ 4 % uncertainty, have been studied as a function of missing energy and missing momentum, and compared to the results of Monte Carlo simulations, obtained from a model based on the distorted-wave impulse approximation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Fire, insect and disease‐caused tree mortalities increased in forests of greater structural diversity during drought

Abstract Structural diversity is an emerging dimension of biodiversity that accounts for size variations in organs among individuals in a community. Previous studies show significant effects of structural diversity on forest growth, but its effects on forest mortality are not known, particularly at a large scale. To address this knowledge gap, we quantified structural diversity using stem structural diversity (SSD) based on both tree diameter and height. We obtained U.S. Forest Service Forest Inventory and Analysis (FIA) data from over 2400 plots across southcentral U.S. forests that have suffered a recent drought. Using data from multiple sampling times, we calculated SSD and compared the relative importance of SSD, species diversity, functional diversity and other stand attributes in determining tree mortalities caused by fire, insects and diseases. We also used FIRETEC, a physics‐based fire model, to test the effect of SSD on canopy consumption by fire. Our results showed that (1) SSD was positively associated with tree mortalities caused by all three disturbances; (2) species richness was negatively associated with insect‐ and disease‐caused mortalities; (3) functional diversity was negatively associated with fire‐ and disease‐caused mortalities and (4) more phylogenetically related species had more similar mortality rates by insect and disease but not fire. Moreover, the FIRETEC model showed increasing canopy consumption by fire in stands with greater SSD. Together, the different tree mortalities during drought associated with SSD more consistently than the other biodiversity metrics were evaluated. Synthesis . Our results suggest that SSD could be considered in modelling forest dynamics and planning management to sustain forest health under disturbances.

54 ENVIRONMENTAL SCIENCES↗