Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Conformer-Dependent Chemistry: Experimental Product Branching of the Vinyl Alcohol + OH + O 2 Reaction

The concentration of formic acid in Earth's tropo-sphere is underestimated by detailed chemical models compared to field observations. Phototautomerization of acetaldehyde to its less stable tautomer vinyl alcohol, followed by the OH-initiated oxidation of vinyl alcohol, has been proposed as a missing source of formic acid that improves the agreement between models and field measurements. Theoretical investigations of the OH + vinyl alcohol reaction in excess O 2 conclude that OH addition to the alpha carbon of vinyl alcohol produces formaldehyde + formic acid + OH, whereas OH addition to the β site leads to glycoaldehyde + HO 2 . Furthermore, these studies predict that the conformeric structure of vinyl alcohol controls the reaction pathway, with the anti-conformer of vinyl alcohol promoting α OH addition, whereas the syn-conformer promotes β addition. However, the two theoretical studies reach different conclusions regarding which set of products dominate. We studied this reaction using time-resolved multiplexed photoionization mass spectrometry to quantify the product branching fractions. Our results, supported by a detailed kinetic model, conclude that the glycoaldehyde product channel (arising mostly from syn-vinyl alcohol) dominates over formic acid production with a 3.6:1.0 branching ratio. Additionally, this result supports the conclusion of Lei et al. that conformer-dependent hydrogen bonding at the transition state for OH-addition controls the reaction outcome. As a result, tropospheric oxidation of vinyl alcohol creates less formic acid than recently thought, increasing again the discrepancy between models and field observations of Earth's formic acid budget.

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

Modeling CH 4 and CO 2 cycling using porewater stable isotopes in a thermokarst bog in Interior Alaska: results from three conceptual reaction networks

Quantifying rates of microbial carbon transformation in peatlands is essential for gaining mechanistic understanding of the factors that influence methane emissions from these systems, and for predicting how emissions will respond to climate change and other disturbances. In this study, we used porewater stable isotopes collected from both the edge and center of a thermokarst bog in Interior Alaska to estimate in situ microbial reaction rates. We expected that near the edge of the thaw feature, actively thawing permafrost and greater abundance of sedges would increase carbon, oxygen and nutrient availability, enabling faster microbial rates relative to the center of the thaw feature. We developed three different conceptual reaction networks that explained the temporal change in porewater CO2, CH4, δ13C-CO2 and δ13C-CH4. All three reaction-network models included methane production, methane oxidation and CO2 production, and two of the models included homoacetogenesis — a reaction not previously included in isotope-based porewater models. All three models fit the data equally well, but rates resulting from the models differed. Most notably, inclusion of homoacetogenesis altered the modeled pathways of methane production when the reaction was directly coupled to methanogenesis, and it decreased gross methane production rates by up to a factor of five when it remained decoupled from methanogenesis. The ability of all three conceptual reaction networks to successfully match the measured data indicate that this technique for estimating in-situ reaction rates requires other data and information from the site to confirm the considered set of microbial reactions. Despite these differences, all models indicated that, as expected, rates were greater at the edge than in the center of the thaw bog, that rates at the edge increased more during the growing season than did rates in the center, and that the ratio of acetoclastic to hydrogenotrophic methanogenesis was greater at the edge than in the center. In both locations, modeled rates (excluding methane oxidation) increased with depth. A puzzling outcome from the effort was that none of the models could fit the porewater dataset without generating “fugitive” carbon (i.e., methane or acetate generated by the models but not detected at the field site), indicating that either our conceptualization of the reactions occurring at the site remains incomplete or our site measurements are missing important carbon transformations and/or carbon fluxes. This model–data discrepancy will motivate and inform future research efforts focused on improving our understanding of carbon cycling in permafrost wetlands.

13CH4↗
DOE OSTI2016

Nitrate vulnerability projections from Bayesian inference of multiple groundwater age tracers

Nitrate is a major source of contamination of groundwater in the United States and around the world. We tested the applicability of multiple groundwater age tracers ( 3 H, 3 He, 4 He, 14 C, 13 C, and 85 Kr) in projecting future trends of nitrate concentration in 9 long-screened, public drinking water wells in Turlock, California, where nitrate concentrations are increasing toward the regulatory limit. Very low 85 Kr concentrations and apparent 3 H/ 3 He ages point to a relatively old modern fraction (40–50 years), diluted with pre-modern groundwater, corroborated by the onset and slope of increasing nitrate concentrations. An inverse Gaussian–Dirac model was chosen to represent the age distribution of the sampled groundwater at each well. Model parameters were estimated using a Bayesian inference, resulting in the posterior probability distribution – including the associated uncertainty – of the parameters and projected nitrate concentrations. Three scenarios were considered, including combined historic nitrate and age tracer data, the sole use of nitrate and the sole use of age tracer data. Each scenario was evaluated based on the ability of the model to reproduce the data and the level of reliability of the nitrate projections. The tracer-only scenario closely reproduced tracer concentrations, but not observed trends in the nitrate concentration. Both cases that included nitrate data resulted in good agreement with historical nitrate trends. Use of combined tracers and nitrate data resulted in a narrower range of projections of future nitrate levels. However, use of combined tracer and nitrate resulted in a larger discrepancy between modeled and measured tracers for some of the tracers. In conclusion, despite nitrate trend slopes between 0.56 and 1.73 mg/L/year in 7 of the 9 wells, the probability that concentrations will increase to levels above the MCL by 2040 are over 95% for only two of the wells, and below 15% in the other wells, due to a leveling off of reconstructed historical nitrate loadings to groundwater since about 1990.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Machine learning for detection of 3D features using sparse x-ray tomographic reconstruction

In many inertial confinement fusion (ICF) experiments, the neutron yield and other parameters cannot be completely accounted for with one and two dimensional models. This discrepancy suggests that there are three dimensional effects that may be significant. Sources of these effects include defects in the shells and defects in shell interfaces, the fill tube of the capsule, and the joint feature in double shell targets. Due to their ability to penetrate materials, x rays are used to capture the internal structure of objects. Methods such as computational tomography use x-ray radiographs from hundreds of projections, in order to reconstruct a three dimensional model of the object. In experimental environments, such as the National Ignition Facility and Omega-60, the availability of these views is scarce, and in many cases only consists of a single line of sight. Mathematical reconstruction of a 3D object from sparse views is an ill-posed inverse problem. These types of problems are typically solved by utilizing prior information. Neural networks have been used for the task of 3D reconstruction as they are capable of encoding and leveraging this prior information. We utilize half a dozen, different convolutional neural networks to produce different 3D representations of ICF implosions from the experimental data. Deep supervision is utilized to train a neural network to produce high-resolution reconstructions. These representations are used to track 3D features of the capsules, such as the ablator, inner shell, and the joint between shell hemispheres. Machine learning, supplemented by different priors, is a promising method for 3D reconstructions in ICF and x-ray radiography, in general.

Wolfe, Bradley T. (ORCID:0000000268301614)↗
DOE OSTI2022

Kinetic model development and Bayesian uncertainty quantification for the complete reduction of Fe-based oxygen carriers with CH 4 , CO, and H 2 for chemical looping combustion

In this work, three kinetic models are developed and calibrated for the complete multi-step reduction of an Fe-based oxygen carrier (OC) particle with CH 4 , CO, and H 2 , using data from thermogravimetric analysis. The complete reduction rate profiles exhibit complex dynamics whose trajectory is significantly different depending on the reducing gas. A Bayesian model building and parameter estimation framework is applied for simultaneous parameter and model structure uncertainty quantification. The final models show excellent agreement between model predictions and calibration data, as well as new data not used for calibration (for the reduction of the OC with CH 4 ). Parameter uncertainty is quantified by determining joint posterior distribution, and model structure uncertainty is addressed by incorporating Gaussian process stochastic functions (represented by Bayesian smoothing splines) into the kinetic models. The final kinetic models with discrepancy functions are readily employable in equation-oriented simulation and optimization platforms.

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

Experimental and Theoretical Studies of Actinide and Lanthanide Ion Transport Across Supported Liquid Membranes

In this paper, we propose a new transport mechanism for metal ions relevant for used nuclear fuel separation processes by a supported liquid membrane (SLM). Two SLM extraction systems were investigated where the membrane was impregnated with either di-(2-ethylhexyl)phosphoric acid (HDEHP) or tributyl phosphate (TBP). A HDEHP impregnated membrane was used to extract neodymium (III), representative of a typical trivalent lanthanide. Cerium, which was oxidized by sodium bismuthate from trivalent to tetravalent state, was extracted by TBP. Oxidized cerium was used as a surrogate for oxidized americium to investigate the kinetics and possibility of americium and curium separation by membrane extraction. Both extraction systems were operated at varying nitric acid concentrations, and changes in the kinetics and extraction efficiency of metal ions were investigated. The proposed transport mechanism that was chosen for our studies was modified from the previous works by Danesi et al.[1,2] and Cussler et al.[3] The mechanism was selected due to the ability to accommodate and describe transport phenomena across a SLM when formation of extractant nano-channels in the membrane may exist. We were able to obtain acceptable fit of the models to our overall data trends although chemical and physical conditions must be well established and purity and homogeneity of the membrane are critical. A reverse transport of metal ions was observed when leaving the system for longer times which agrees with our model. The membrane was investigated for degradation and shown to be stable after contact with up to 7 M nitric acid for over 2000 minutes. Finally, we examined the possibility of partitioning americium from curium using a SLM impregnated by TBP. Separation of americium from curium was observed although not to a degree that was expected based on the Ce(IV) transport. Incomplete oxidation of Am(III) to Am(V) and reduction of Am(VI) on the membrane surface are possible causes for this observed discrepancy. Our model was, however, able to accurately predict Cm(III) transport through the membrane.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Assessing Clouds in GFDL's AM4.0 With Different Microphysical Parameterizations Using the Satellite Simulator Package COSP

We evaluate cloud simulations using satellite simulators against multiple observational data sets. These simulators have been run within the Geophysical Fluid Dynamics Laboratory's Atmosphere Model version 4.0 (AM4.0), as well as an alternative configuration where a fully two‐moment Morrison‐Gettelman cloud microphysical parameterization with prognostic precipitation (MG2) is applied, denoted as AM4‐MG2. The modeled cloud spatial distributions, vertical profiles, phase partitioning, cloud‐to‐precipitation transitions, and radiative effects compare reasonably well with satellite observations. Model biases include the under‐prediction of total and low‐level clouds, especially optically thin/intermediate clouds with cloud optical depth of less than 23, but the over‐prediction of thick clouds, indicating “too few, too bright” biases. These biases counteract each other, and give rise to reasonable estimates of cloud radiative effects. The underestimate of low‐level clouds is associated with too early and too frequent drizzle/precipitation formation. The precipitation bias is improved in AM4‐MG2, where the autoconversion scheme initiates the precipitation more realistically. There also exist discrepancies between models and observations for midlevel and high‐level clouds. Additional biases include the underestimate of liquid cloud fraction and the overestimate of ice cloud fraction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Non -degenerate marginal-likelihood calibration with application to quantum characterization

Here, we propose a marginal likelihood strategy within the Kennedy-O’Hagan (KOH) Bayesian framework, where a Gaussian process (GP) models the discrepancy between a physical system and its simulator. Our approach introduces a novel marginalized likelihood by integrating out the degenerate eigenspace of the covariance matrix, rather than approximating the original likelihood. Unlike approximation methods that compromise accuracy for computational efficiency, our method defines an exact likelihood—distinct from the original but preserving all relevant information. This formulation achieves computational efficiency and stability, even for large datasets where the covariance matrix nears degeneracy. Applied to the characterization of a superconducting quantum device at Lawrence Livermore National Laboratory, the approach enhances the predictive accuracy of the Lindblad master equations for modeling Ramsey measurement data by effectively quantifying uncertainties consistent with the quantum data.

general physics↗
DOE OSTI2023

Northern Hemisphere Snow Drought in Earth System Model Simulations and ERA5‐Land Data in 1980–2014

Abstract Low snow levels over the past few decades and predictions of a low‐to‐no snow future have spurred research into snow droughts, which pose a threat to water security and management. Systematic data‐model comparisons of snow drought have been lacking, hindering our understanding of the drivers of snow drought in the past. To address this gap, we analyzed snow drought events using standardized snow water equivalent index derived from monthly results of four numerical experiments using the E3SM Land Model (ELM) and ERA5‐Land data during the period of 1980–2014. Additionally, we compared snow drought duration calculated from models with those from the ERA5‐Land data during selected El Niño‐Southern Oscillation (ENSO) years. The numerical experiments were conducted with ELM driven by two prescribed atmospheric forcings, and with the coupled land‐atmosphere configuration of E3SM with and without plant hydraulics scheme feedback. Analysis reveals that 20%–30% of snow droughts occur due to factors other than above‐normal temperature and low snowfall, such as low soil moisture, warm soil temperature, and low relative humidity, etc., especially in high latitudes (50° North). Furthermore, our study highlights the exacerbating effect of ENSO events on snow drought conditions in various regions, despite some discrepancies between model and ERA5‐Land results. We also identified limitations of the coupled land‐atmosphere models in our current configuration in capturing the spatial patterns of snow droughts. This study underscores the challenge of predicting and mitigating snow drought and the need for a comprehensive understanding of the factors contributing to snow drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Peatland Loss in Southeast Asia Contributing to U.S. Biofuel’s Greenhouse Gas Emissions

Land use change (LUC) induced by biofuel production could lead to greenhouse gas (GHG) emissions, which potentially increase biofuel’s carbon intensity. Among the sources of LUC-related emissions for soy biodiesel, the contribution from peatland loss to agricultural plantations in Southeast Asia remains uncertain. Here, in this study, we analyzed LUC in Malaysia and Indonesia and modeled its impacts on the GHG emissions of soy biodiesel produced in the United States. It shows that oil palm plantations have more than doubled over 2001–2016 and the area of palm-on-peatlands (PoP) has expanded 3.7 times. Over new palm plantations, the share of PoP is about 19% regardless of time and location and the emission factor (EF) for peatland-to-palm conversion is estimated to be 41.5 Mg CO 2 ha –1 yr –1 . With these updates on PoP and EF, the contribution of peatland loss (0.7–5.1 g CO 2 e MJ –1 ) to biodiesel emissions is only 40–65% of previous estimates, which reduces discrepancies among model simulations used by different agencies. Based on emerging evidence on LUC and related carbon changes, our analysis reexamines regional peatland loss and its impacts on LUC emissions modeling and provides new insights into the estimation of LUC impacts on biofuels’ carbon intensity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

CONSTRAINING THE ABUNDANCES OF COMPLEX ORGANICS IN THE INNER REGIONS OF SOLAR-TYPE PROTOSTARS

The high abundances of Complex Organic Molecules (COMs) with respect to methanol, the most abundant COM, detected toward low-mass protostars, tend to be underpredicted by astrochemical models. This discrepancy might come from the large beam of the single-dish telescopes, encompassing several components of the studied protostar, commonly used to detect COMs. To address this issue, we have carried out multi-line observations of methanol and several COMs toward the two low-mass protostars NGC 1333-IRAS 2A and -IRAS 4A with the Plateau de Bure interferometer at an angular resolution of 2″, resulting in the first multi-line detection of the O-bearing species glycolaldehyde and ethanol and of the N-bearing species ethyl cyanide toward low-mass protostars other than IRAS 16293. The high number of detected transitions from COMs (more than 40 methanol transitions for instance) allowed us to accurately derive the source size of their emission and the COM column densities. The COM abundances with respect to methanol derived toward IRAS 2A and IRAS 4A are slightly, but not substantitally, lower than those derived from previous single-dish observations. The COM abundance ratios do not vary significantly with the protostellar luminosity, over five orders of magnitude, implying that low-mass hot corinos are quite chemically rich as high-mass hot cores. Astrochemical models still underpredict the abundances of key COMs, such as methyl formate or di-methyl ether, suggesting that our understanding of their formation remains incomplete.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (N d ) comparisons are highly variable depending on the N d retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and N d values due to a lower cloud effective radius than observed. However, the simulated albedo response to N d is suppressed due to a correlation between the solar zenith angle (SZA) and N d created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln N d distribution, which increases the cloud albedo sensitivity to ln N d to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of N d to CCN is greater than observed for TOA-retrieved N d , and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved N d . The LWP response to N d in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–N d relationship, observed clouds become much more adiabatic as N d increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–N d relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–N d relationship non-causal. N d retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2019

Physics-Informed CoKriging: A Gaussian-Process-Regression-Based Multifidelity Method for Data-Model Convergence

In this work, we propose a new Gaussian process regression (GPR)-based multifidelity method: physics-informed CoKriging (CoPhIK). In CoKriging-based multifidelity methods, the quantities of interest are modeled as linear combinations of multiple parameterized stationary Gaussian processes (GPs), and the hyperparameters of these GPs are estimated from data via optimization. In CoPhIK, we construct a GP representing low-fidelity data using physics-informed Kriging (PhIK), and model the discrepancy between low- and highfidelity data using a parameterized GP with hyperparameters identified via optimization. Our approach reduces the cost of optimization for inferring hyperparameters by incorporating partial physical knowledge. We prove that the physical constraints in the form of deterministic linear operators are satisfied up to an error bound. Furthermore, we combine CoPhIK with a greedy active learning algorithm for guiding the selection of additional observation locations. The eciency and accuracy of CoPhIK are demonstrated for reconstructing the partially observed modified Branin function, reconstructing the sparsely observed state of a steady state heat transport problem, and learning a conservative tracer distribution from sparse tracer concentration measurements.

Yang, Xiu↗
DOE OSTI2026

First study of the nuclear response to fast hadrons via angular correlations between pions and slow protons in electron-nucleus scattering

We report on the first measurement of angular correlations between high-energy pions and slow protons in electron-nucleus (eA) scattering, providing a new probe of how a nucleus responds to a fast-moving quark. The experiment employed the CLAS detector with a 5-GeV electron beam incident on deuterium, carbon, iron, and lead targets. For heavier nuclei, the pion-proton correlation function is more spread-out in azimuth than for lighter ones, and this effect is more pronounced in the πp channel than in earlier ππ studies. The proton-to-pion yield ratio likewise rises with nuclear mass, although the increase appears to saturate for the heaviest targets. These trends are qualitatively reproduced by state-of-the-art eA event generators, including BeAGLE, eHIJING, and GiBUU, indicating that current descriptions of target fragmentation rest on sound theoretical footing. At the same time, the precision of our data exposes model-dependent discrepancies, delineating a clear path for future improvements in the treatment of cold-nuclear matter effects in eA scattering.

Correlations↗
DOE OSTI2024

Synthesis and characterization of isotopically barcoded nickel, molybdenum, and tungsten taggants for intentional nuclear forensics

Intentional nuclear forensics is a concept wherein the deliberate addition of benign and persistent material signatures to nuclear material can be used to reduce the time between the discovery of material outside of regulatory control and determination of its original provenance. One concept within intentional nuclear forensics involves the use of perturbed stable isotopes to generate unique isotope ratio “barcodes” to encode information (e.g., production batch, location, etc.) and track material throughout the nuclear fuel cycle. Synthesis of taggant species of nickel (Ni), molybdenum (Mo), and tungsten (W) was undertaken via a double-spike mechanism, wherein two highly enriched isotopes of interest per elemental taggant were mixed to form an enriched “double-spike” which was subsequently isotopically diluted with bulk material having a natural isotopic composition. Two taggant species perturbing isotopic ratios, alpha (α) and beta (β), for each of Ni, Mo, and W were synthesized. Independent measurements of double spikes and alpha and beta taggant species agreed within uncertainty and are clearly resolvable from natural compositions. High-precision analyses were independently performed by MC-ICP-MS at two U.S. National Laboratories, with consensus values and uncertainties calculated for all samples. Observed isotopic perturbations in the final taggant species measured on the order of hundreds to thousands of permille (‰) with respect to natural for isotope ratios of interest (e.g., 60 Ni/ 58 Ni, 100 Mo/ 98 Mo, 186 W/ 183 W). Discrepancies between modeled and measured isotopic compositions were observed and are largely attributed to imprecise vendor assay values for starting materials. Using measured starting material compositions as inputs for the mixing model improved the level of agreement between predicted and measured α and β taggant isotope ratios. Overall, characterization of all taggant species demonstrates that this “barcode” concept could have viability for use in nuclear forensics. Finally, it is expected that for any two-isotope mixing array dozens of isotopic barcodes could be encoded into a material system and subsequently resolved utilizing modern mass spectrometric methods.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗