Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “response surface methods”

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

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗
DOE OSTI2022

Experimental-based mechanistic study and optimization of hydrothermal liquefaction of anaerobic digestates

Valorization of agricultural and food waste digestates is crucial for sustainable waste management to reduce environmental impacts and improve the economics of commercial farms. Hydrothermal liquefaction (HTL) of anaerobic digestates was evaluated to recover resources by converting them into carbon-dense biocrude oil and a nutrient-rich HTL aqueous phase (HTL-AP) coproduct. The effects of HTL temperature (280–360 °C), reaction time (10–50 min), feedstock pH (2.5–8.5), digestate salt content (1–5 wt%), and digestate cellulose-to-lignin ratio (0.2–1.8) on energy and nutrient recovery were systematically investigated in a set of well-designed experiments following a half-fractional central composite protocol. Response surface analysis combined with HTL product characterization and comparative literature study produced a comprehensive reaction pathway for HTL of anaerobic digestates. Moreover, this analysis revealed the importance of acidic feedstocks (pH 3.00–5.53), high reaction temperatures (337–360 °C), and reaction times <45 or 45–50 min for digestates with Cel/Lig >1 or <1, for maximizing the energy recovered in biocrude (high carbon yield and low heteroatom content) and the amounts of P, NH 3 –N, and Mg distributed in the HTL-AP. Acidic conditions catalyzed biocrude production, inhibited the Maillard reaction (lowering the nitrogen content in biocrude), and partitioned nutrients into the HTL-AP. Higher reaction temperatures coupled with longer reaction times activated hydro-denitrogenation and deoxygenation reactions to improve biocrude quality. Finally, this work provides not only validated methods to achieve targeted resource recovery for specific feedstock compositions using HTL, but also a comprehensive mechanistic understanding of the HTL of biomass waste for controlling target product characteristics.

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

Polymer nanoparticle photocatalysts realized in non-aqueous solvents

Colloidal organic nanoparticles (oNPs) have emerged as a promising category of photocatalyst, thanks to their long-lived surface-bound charges, electronic tunability, and strong absorption in the visible spectrum. Our previous research has established a direct correlation between charge generation in oNPs and their photocatalytic activity, highlighting their effectiveness as a framework for stable, long-lived free carriers. However, oNPs have been restricted to use only in aqueous environments as a result of being synthesized via either nano-emulsion or nano-precipitation procedures. Herein, we present a method for transferring oNP photocatalysts from water into polar non-aqueous solvents while retaining their long-term colloidal stability. We observed that the polymer chains in the solvent-transferred oNPs rearrange from a predominantly H-aggregate structure in water to a combination of H- and J-aggregate characteristics in N,N-dimethylformamide, suggesting a dynamic rearrangement in response to the new solvent environment. Importantly, transient absorption and time-resolved microwave conductivity measurements confirm that the solvent-transferred oNPs maintain their ability to generate free charges at an internal heterojunction. This development opens unique opportunities for eventually leveraging light-generated, long-lived electrons and holes in synthetic redox chemistry across diverse solvent environments, a direction that will be explored in future studies.

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

A machining digital twin for hybrid manufacturing

Hybrid manufacturing consisting of metal additively manufactured preforms and computer numerical control (CNC) machining has been established to be an effective method for high material use rates. However, hybrid manufacturing introduces unique challenges. Near-net shape designs are typically selected, which result in a smaller margin for part placement within the stock and stringent requirements for work coordinate system identification. Additionally, less stock material reduces the preform stiffness, which limits the material removal rates during machining. This paper demonstrates a digital twin for CNC machining of a wire arc additively manufactured preform that implements: 1) structured light scanning for stock model identification and tool path generation; 2) a fused filament fabrication apparatus to attach temporary fiducials and scan targets to the preform that enable coordinate system definition for both the CAM and CNC machine; 3) preform and tool tip frequency response function measurements to enable stable milling parameter selection; and 4) post-manufacturing measurements of geometry, surface finish, and structural dynamics to confirm designer intent. These efforts define key components of the machining digital twin for hybrid manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm

Retreat and advance of ice sheets perturb the gravitational field, solid surface and rotation of the Earth, leading to spatially variable sea-level changes over a range of timescales $\textit{O}$ (10 0–6 years), which in turn feed back onto ice-sheet dynamics. Coupled ice-sheet–sea-level models have been developed to capture the interactive processes between ice sheets, sea level and the solid Earth, but it is computationally challenging to capture short-term interactions $\textit{O}$ (10 0–2 years) precisely within longer $\textit{O}$ (10 3–6 years) simulations. The standard forward sea-level modelling algorithm assigns a uniform temporal resolution in the sea-level model, causing a quadratic increase in total CPU time with the total number of input ice history steps, which increases with either the length or temporal resolution of the simulation. In this study, we introduce a new “time window” algorithm for 1D pseudo-spectral sea-level models based on the normal mode method that enables users to define the temporal resolution at which the ice loading history is captured during different time intervals before the current simulation time. Utilizing the time window, we assign a fine temporal resolution $\textit{O}$ (10 0–2 years) for the period of ongoing and recent history of surface ice and ocean loading changes and a coarser temporal resolution $\textit{O}$ (10 3–6 years) for earlier periods in the simulation. This reduces the total CPU time and memory required per model time step while maintaining the precision of the model results. We explore the sensitivity of sea-level model results to the model temporal resolution and show how this sensitivity feeds back onto ice-sheet dynamics in coupled modelling. We apply the new algorithm to simulate sea-level changes in response to global ice-sheet evolution over two glacial cycles and the rapid collapse of marine sectors of the West Antarctic Ice Sheet in the coming centuries and provide appropriate time window profiles for each application. The time window algorithm reduces the total CPU time by ~ 50 % in each of these examples and changes the trend of the total CPU time increase from quadratic to linear. This improvement would increase with longer simulations than those considered here. Our algorithm also allows for coupling time intervals of annual temporal scale for coupled ice-sheet–sea-level modelling of regions such as West Antarctica that are characterized by rapid solid Earth response to ice changes due to the thin lithosphere and low mantle viscosities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Carbon decorated Li 3 V 2 (PO 4 ) 3 for high-rate lithium-ion batteries: Electrochemical performance and charge compensation mechanism

Fast charging and high-power delivering batteries are highly demanded in mobile electronics, electric vehicles and grid energy storage, but there are full of challenges. The star-material Li 3 V 2 (PO 4 ) 3 is demonstrated as a promising high-rate cathode material meeting the above requirements. Herein, we report the carbon decorated Li 3 V 2 (PO 4 ) 3 (LVP/C) cathode prepared via a facile method, which displays a remarkable high-rate capability and long-term cycling performance. Briefly, the prepared LVP/C delivers a high discharge capacity of 122 mAh g -1 (~93% of the theoretical capacity) at a high rate up to 20 C and a superior capacity retention of 87.1% after 1000 cycles. Importantly, by applying a combination of X-ray absorption spectroscopy and full-range mapping of resonant inelastic X-ray scattering, we clearly elucidate the structural and chemical evolutions of LVP upon various potentials and cycle numbers. Furthermore, we show unambiguous spectroscopic evidences that the evolution of the hybridization strength between V and O in LVP/C as a consequence of lithiation/delithiation is highly reversible both in the bulk and on the surface during the discharge-charge processes even over extended cycles, which should be responsible for the remarkable electrochemical performance of LVP/C. Our present study provides not only an effective synthesis strategy but also deeper insights into the surface and bulk electrochemical reaction mechanism of LVP, which should be beneficial for the further design of high-performance LVP electrode materials.

25 ENERGY STORAGE↗
DOE OSTI2020

Comparison of Thin Film Heat Flux Gauge Technologies Emphasizing Continuous-Duration Operation

Abstract Thin-film heat flux gauges (HFGs) have been used for decades to measure surface temperatures and heat flux in test turbines with the majority being used in facilities that are short-duration. These gauges are typically composed of two resistive temperature devices deposited on opposing sides of a dielectric. However, because these sensors have been traditionally applied for measurements in transient-type facilities, the challenges facing adaptation of this technology for a steady facility warrant investigation. These challenges are highlighted, and the solutions are presented throughout the paper. This paper describes the nanofabrication process for heat flux gauges and a new calibration method to address the potential deterioration of gauges over long runtimes in continuous-duration facilities. Because the primary uncertainty of these sensors arises from the ambiguity of the thermal properties, the emphasis is placed on the property determination. Also, this paper presents a discussion on the use of impulse response theory to process the data showing the feasibility of the method for steady-duration facilities after an initial settling time. The latter portion of the paper focuses on comparing well-established heat flux gauges developed for short-duration turbine test facilities to recently developed gauges fabricated using modern nanofabrication techniques for a continuous turbine test facility. The gauges were compared using the test case of an impinging jet over a range of Reynolds numbers. The comparison between the PSU gauge and the reference device indicated agreement within 14%, and similar results were achieved through comparison with established sensors from partner institutions.

Engineering↗
DOE OSTI2024

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗
DOE OSTI2026

Time-domain terahertz emission spectroscopy on van der Waals materials

Time-domain terahertz (THz) emission spectroscopy provides a direct method to probe transient photo-currents by recording the emitted terahertz electric field. Although the basic principles of THz surface emission have been understood for more than 30 years, the constant progress in ultrafast laser science to ever shorter pulses, the development of new materials and enhanced sensitivity promote THz emission spectroscopy as a reliable method to gain insights into charge carrier dynamics with unprecedented precision. It provides a versatile tool to study ultrafast processes, such as plasmon-driven hot carriers, dynamics of Dirac fermions, interfacial charge transfer, coherent phonon emission and quantum beating, to name only a few. However, despite the rapidly growing body of research on van der Waals materials, especially in their low-dimensional limit, THz emission spectroscopy has only been applied to a limited extent in these material systems. In this prospective, we review time-domain THz emission spectroscopy as a complementary approach to probe ultrafast charge carrier dynamics and the material’s nonlinear response. After a description of the experimental method, we report on THz emission spectroscopy of bulk and 2D van der Waals materials with special focus on graphene and transition metal dichalcogenide layers.

2D materials↗
DOE OSTI2021

Detection, monitoring, and determination of location of changes in metallic structures using multimode acoustic signals

Acoustic transducers generate and receive acoustic signals at multiple locations along a surface of rigid structure, wherein longitudinal spacing between transducer locations define measurement zones. Acoustic signals with chosen amplitude-time-frequency characteristics excite multiple vibration modes in the structure within each zone. Small mechanical changes in inspection zones lead to scattering and attenuation of broadband acoustic signals, which are detectable as changes in received signal characteristics as part of a through-transmission technique. Additional use of short, narrowband pulse acoustic signals as part of a pulse-echo technique allows determination of the relative location of the mechanical change within each zone based on the differential delay profiles. For accurate acoustic modeling and simulation, the mesh size, time step, time delay, and time-window size are optimized. Frequency normalization of the Short-Time Fourier Transform of acoustic response output improves experiment-simulation cross-validation. Applications of the method to structures with arbitrarily complex geometries are also demonstrated.

Findikoglu, Alp Tugrul↗
DOE OSTI2021

Application of an automated machine learning-genetic algorithm (AutoML-GA) coupled with computational fluid dynamics simulations for rapid engine design optimization

In recent years, the use of machine learning-based surrogate models for computational fluid dynamics (CFD) simulations has emerged as a promising technique for reducing the computational cost associated with engine design optimization. However, such methods still suffer from drawbacks. One main disadvantage is that the default machine learning (ML) hyperparameters are often severely suboptimal for a given problem. This has often been addressed by manually trying out different hyperparameter settings, but this solution is ineffective in case of a high-dimensional hyperparameter space. Besides this problem, the amount of data needed for training is also not known a priori. In response to these issues that need to be addressed, the present work describes and validates an automated active learning approach, AutoML-GA, for surrogate-based optimization of internal combustion engines. In this approach, a Bayesian optimization technique is used to find the best machine learning hyperparameters based on an initial dataset obtained from a small number of CFD simulations. Subsequently, a genetic algorithm is employed to locate the design optimum on the ML surrogate surface. In the vicinity of the design optimum, the solution is refined by repeatedly running CFD simulations at the projected optima and adding the newly obtained data to the training dataset. It is demonstrated that AutoML-GA leads to a better optimum with a lower number of CFD simulations, compared to the use of default hyperparameters. The proposed framework offers the advantage of being a more hands-off approach that can be readily utilized by researchers and engineers in industry who do not have extensive machine learning expertise.

Owoyele, Opeoluwa↗
DOE OSTI2024

Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies

Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.

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

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigating Soil Organic Matter Complexation using Spectral Induced Polarization

Spectral induced polarization (SIP) laboratory experiments were conducted to determine the sensitivity of this method to the formation of soil organic matter (SOM) complexes, with a long-term goal of field-scale monitoring. There are few SIP experiments that have explored this topic, yet understanding the dynamic behavior and interactions of SOM at the field scale could provide insight into soil fertility and health which influences crop yields, microorganisms that degrade organic pollutants, and carbon stabilization. We present the results of three experiments where the iron oxide, ferrihydrite (Fhy), was used to coat different media, and then the OM compound pentaglycine (PG) was pulse injected to form SOM complexes. SIP data was collected during these injections to capture any surface complexation changes. These experiments were performed in 1) a fluidic cell containing a micromodel, 2) a column containing Fhy coated ceramic beads and 3) a column containing Fhy coated Accusand®. Our results show a higher frequency response (defined here as > 1 Hz) in all three experiments, with the largest amplitude response after the first PG injection (Figure S.1). The repeatability of this response is encouraging and supporting data collected on the Accusand® experiment provides preliminary insight into the mechanisms controlling the SIP signatures. Sampling of fluid conductivity $σ_w$ and pH may indicate deprotonation of SOM occurring or rapid adsorption and release of protons from the Fhy sites. However additional experiments are needed to identify and confirm the primary and secondary reactions impacting the SIP response. We are looking towards other opportunities to continue this work, particularly to repeat experiments while collecting supporting datasets.

58 GEOSCIENCES↗
DOE OSTI2020

Ambient Mass Spectrometry Rapid Analysis Methods for Trace Inorganics

The ability to rapidly detect and characterize the chemistry of trace quantities of inorganic species has multiple potential applications, including nuclear safeguards, nuclear forensics, and environmental monitoring. Some mass spectrometric methods currently exist, but they are costly and time-intensive, often destroying chemical information in the process of ionization. Ambient mass spectrometry is a relatively new analytical method in which samples are introduced to the mass spectrometer from atmospheric pressure, rather than from high vacuum. Ambient MS can operate with gentler ionization methods such as electrospray ionization (ESI) and paper-spray ionization (PSI), reducing fragmentation and preserving the chemical information of the analyte. This may be useful for determining the intended purpose of material in a sample. Strontium, part of nuclear fallout, is analyzed here to test these new analytical methods for viability. Electrospray ionization mass spectrometry functions by spraying the sample as a solution through a needle charged to a few kV. A nebulizing gas flows out along with the sample to create the spray before droplets enter a drying chamber with heated drying gas to remove the solvent molecules. From there bare analyte ions enter the mass spectrometer. This method has the advantage of relatively fast sample preparation and stable ionization. More stable ionization produces higher-quality data. Paper-spray ionization mass spectrometry operates by applying sample to a paper triangle. The paper triangle is then electrified in a metal clamp to a few kV before spray solvent is applied. The spray solvent transports analyte in the sample to the tip, where it is ionized and enters the mass spectrometer. One advantage of this method is reduced sample preparation - an analyst can even swipe a surface and cut a triangle directly from the swipe. The data analysis serves to rapidly process tens of samples all at once. The mass spectrometer outputs a file containing a few hundred spectra for each sample. The pipeline processes each one to identify peaks of interest. A full palette of data analysis metrics are generated, including plots of the raw data and plots of calculated values such as isotopic ratios. The pipeline filters the raw data to remove noise and background peaks, as well as peaks representing detector saturation. A heat map is generated to show the overall spectrum in mass-time- intensity space. Post-processing, we obtained the desired calibration curves, demonstrating that both ESI-MS and PSI-MS are viable methods for quantifying trace amounts of inorganics while retaining their chemical speciation information. ESI-MS demonstrated a linear response region of 100 ppb to 20 ppm, with a limit-of-detection of around 100 ppb. At top left is an example spectrum with the {sup 88}SrNO{sub 3} peak identified to show the retention of chemical speciation. PSI-MS showed a linear response region of around 1 ppb to 1 ppm, with a limit-of-detection of near 1 ppb. PSI-MS produced more erratic results than ESI-MS, likely due to the variability in ionization resulting from manual positioning of the paper triangle near the MS inlet. Future efforts will focus on improving the repeatability of PSI-MS as well as analyzing other species of interest such as cerium and uranium. Improvements will also be made to create more reliable data extraction and intelligent data filtering.

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

Substantial and overlooked greenhouse gas emissions from deep Arctic lake sediment - supporting data and code

This data package contains data, descriptions, and code-based analyses that were used to support conclusions drawn in “Substantial and overlooked greenhouse gas emissions from deep Arctic lake sediment”, by Freitas et al. (2025) (https://doi.org/10.1038/s41561-024-01614-y). The study evaluated greenhouse gas production along a deep sediment core (20 m) taken in 2018 from below Goldstream Lake, a field site approximately 15 km north of Fairbanks, Alaska.The file “ESSDive_NFreitas_2024_flmd.csv” includes an overview of all other csv files in this data package, namely: sediment descriptions (depth and type of sediment), sediment characterizations (bulk density, gravimetric water content, total carbon, etc.), calculated respiration and temperature sensitivity values associated with year-long incubations of the sediment core, and the R code used to process the dataset. Additional details regarding the content of these files and how the calculations were performed are described in the Methods section of this archive. The “data_dictionary_ESSDive_NFreitas_2024_dd.csv” is a data dictionary for all files included in the data package. Each row in the data dictionary represents a column name in a given file. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).Data package updates- August 2024: Title was updated, R code now includes a process for calculating potential production across the whole sediment column, the "sediment_column_production_potentials_ESSDive_NFreitas_2024.csv" is the associated data output file, and the file level metadata and data dictionaries were updated to reflect the contents of the additional csv file.- October 2024: Title was updated, and a statement was added to Step 2 of the Methods that describes that samples were collected (and exported) in a responsible manner and in accordance with relevant permits and local laws.- January 2025: Updated the associated manuscript details in the Abstract (manuscript title, publication date, DOI link) and in the Related References section (full citation).

54 ENVIRONMENTAL SCIENCES↗